Experiment
Below, the different panels of the Experiment tab are described in detail.
Select one or two Transport Mechanisms by checking the Advection and/or Diffusion boxes in the upper left part of the Experiment tab. Checking Simulate Diffusion activates the computation of transport by diffusion. If the button is unchecked, the diffusivity is set to in the advection-diffusion equation. Checking Simulate Advection activates the computation of advective transport by the fluid. If unchecked, the velocity of the fluid is zero, i.e. in the advection-diffusion equation, and therefore, the concentration field is transported by diffusion only. In this case, the Flow Boundary Conditions, the Flow Solver tab, and the Inflow Concentration are inactive. |
In the Simulated Time panel, the physical End Time of the simulation can be specified. The initial time is always zero. By clicking the Approximate End Time button, the End Time is automatically determined such that, for the specified average fluid velocity (defined in the Flow Boundary Conditions panel by Mean Velocity or by Flow Rate), most of the solute passes through the entire structure and a complete breakthrough curve is computed. When it is unclear how to choose the End Time, it is highly recommended that the Approximate End Time used before starting the simulation. This works best for high Péclet numbers, i.e., when advection dominates the transport. The simulation time can be divided into intervals (so-called "time batches"), after which the concentration field and other information is written to disk. The total number of batches is defined via Number of Batches. Depending on the time period to be resolved more finely, one of the following options can be chosen:
|
Estimated Dimensionless Numbers
In this panel, three Dimensionless Numbers are automatically estimated based on the defined parameters, provided a structure is loaded. They give the user an idea of how strongly advection and diffusion contribute to transport before the simulation is started. After the simulation, more precise values are listed in the report. The Courant number quantifies the extent to which particles, molecules, or concentration fields are transported through the structure by advection during the simulation. For instance, if , then the solute (particles, molecules, or the concentration field) is advectively transported by the flowing solvent across half of the domain. The Courant number is defined as where is the characteristic velocity estimated from the flow boundary conditions, is the simulation time (which equals the End Time, given that the initial time is zero), and is the length of the structure in the through direction. In contrast, the Fourier number quantifies the extent to which the solute (particles, molecules, or the concentration field) is transported through the structure by diffusion during the simulation. The Fourier number is defined as where is the characteristic diffusivity, taken as the maximum effective diffusivity across all existing materials. Finally, the Péclet number is the ratio of advective transport to diffusive transport. A Péclet number of infinity implies purely advective transport, whereas a value equal to zero corresponds to purely diffusive transport. The Péclet number is defined as |
Transport Initial and Boundary Conditions
In the Transport Initial and Boundary Conditions, the initial concentration distribution in the structure and the Inflow Concentration can be defined. Inflow Concentration is only available if Simulate Advection is checked as otherwise the inflow boundary is an empty set. When the check box next to Initial Concentration Field .hht File is checked, a concentration field can be loaded by clicking Browse. Through this option, an initial distribution of the solute in the structure can be defined. For each time batch, AddiDict automatically generates a *.hht file, so the output of a previous transport simulation can easily serve as the input for subsequent runs. Alternatively, *.hht files can be manipulated with Python scripting. In the Inflow Concentration field, a concentration value that represents the amount of solute carried into the structure by the inflowing solvent fluid can be set. Currently, the distribution of one solute species can be simulated in a Transport Concentration Field simulation. |
In the Flow Boundary Conditions panel, the Pressure Drop, the mean flow Velocity, or the Flow Rate can be entered. In AddiDict, the flow is always computed in the Z-direction. If another direction is wanted, the 3D model can be rotated with ProcessGeo prior to using AddiDict. Note that this part of the dialog is accessible only if Simulate Advection is checked. From the Flow PDE menu, the user specifies the equations to be solved to compute the fluid flow field required to simulate the transport of the solute concentration field. (Navier–)Stokes equations describe fluid flow in structures containing only solid materials and pores. (Navier–)Stokes–Brinkman equations describe the flow through structures also containing porous materials. Use the Stokes equation when fluid velocities are slow or the viscosities are high, i.e. at low Reynolds numbers. The Navier–Stokes equation should be used for high fluid velocities when inertial effects cannot be neglected. The used flow boundary conditions in X and Y direction are reported and the boundary condition in Z direction can be defined. The boundary conditions in X and Y direction are fixed to Symmetric. The boundary condition in Z direction can be set to Periodic, Symmetric (default), or VinPout. Please refer to the FlowDict user guide for more details. Additionally, the Slip Length can be specified in the Flow Boundary Conditions panel. The Slip Length allows including sliding effects in the simulation. The default Slip Length of zero corresponds to a flow velocity of zero along the structure. A non-zero Slip Length simulates the sliding of the fluid along the structure’s solid parts, increasing the fluid mean flux and thus, the permeability. This option might be used when it is realistic for a given physical material. Currently, the same slip length value must be set for all materials in the structure. The slip length is defined as an extrapolated distance inside the wall where the tangential velocity component vanishes. Find more details on the Slip Length here. For the No-Slip boundary condition two additional options are available:
Find more details on these settings in the FlowDict user guide. |
©2025 created by Math2Market GmbH / Imprint / Privacy Policy