STOMP

Boundary Conditions Card Options (W)

This card is used to specify boundary conditions for the aqueous phase system, passive solutes or reactive species.  For the aqueous system, boundary conditions are one of two types: 1) Dirichlet or 2) Neumann.  Dirichlet-type boundary conditions specify the aqueous pressure at the boundary surface centroid. Neumann-type boundary conditions specify an aqueous flux (i.e., Darcy velocity) across a boundary surface. The mechanism for specifying these two types of boundary types, however, vary from simple to complex. Solutes are nonreactive passive tracers; where, passive indicates that aqueous-phase properties are independent of solute concentrations. A boundary condition type must be specified for every solute.  Species are reactive passive tracers; where, passive indicates that aqueous-phase properties are independent of species concentrations. A single boundary condition type is specified for all species.

Aqueous Boundary Condition Options

Dirichlet Type
Dirichlet

Dirichlet-type boundary where the aqueous pressure is specified at the boundary surface centroid.

Dirichlet Outflow

Dirichlet-type boundary where the aqueous pressure is specified at the boundary surface centroid, but flow is only permitted in the direction of the surface normal (i.e., out of the computational domain).

Hydraulic Gradient

Dirichlet-type boundary where the aqueous pressure is specified at the boundary surface centroid of the lowest I, J, K indexed node.  The aqueous pressure at all other boundary surface centroids within the range of I, J, K indexed nodes is set to be in hydrostatic equilibrium with the aqueous pressure at the boundary surface centroid of the lowest I, J, K indexed node. This boundary type should be applied only to a column or plane of vertical surfaces.

 
Initial Condition

Dirichlet-type boundary where the aqueous pressure is set and to and held at the initial value specified at the node adjacent to the boundary surface. This boundary type is invariant with time. If used with a restart simulation, the aqueous pressure for each specified boundary node is set at the value that is present in the restart file.

 
Saturated

Dirichlet-type boundary where the aqueous pressure is set to the gas pressure at the boundary surface centroid.  The gas pressure is invariant in the simulation, but can be specified via the Initial Conditions Card. This boundary type is available only for two-phase conditions and imposes total-liquid saturation conditions (e.g., water table) at the boundary surface. 

 
Seepage Face

Dirichlet-type boundary similar to a Hydraulic Gradient boundary, but limited to pressure boundaries of the local gas pressure. The aqueous pressure is specified at the boundary surface centroid of the lowest I, J, K indexed node. This boundary type is designed to model an exposed vertical face that seeps liquids. Liquid can enter a seepage face only for phase pressures that exceed the local gas pressure.  

Unit Gradient

Dirichlet-type boundary where hydrostatic conditions are imposed across the boundary surface for the specified phase.

X-Y-Z Hydraulic Gradient

Dirichlet-type boundary which is equivalent to the Hydraulic Gradient boundary condition, but with the option to specify additional gradients in the x-direction (global coordinates), y-direction (global coordinates), and z-direction (global coordinates).

X-Y-Z Seepage Face

Dirichlet-type boundary which is equivalent to the Seepage Face, but with the option to specify additional gradients in the x-direction (global coordinates), y-direction (global coordinates), and z-direction (global coordinates).

Neumann Type
Neumann

Neumann-type boundary where the aqueous volumetric flux is specified across the boundary surface.  Positive flux is in the direction of the surface normal.

Zero Flux

Neumann-type boundary where the aqueous volumetric flux is zero.

Solute Boundary Condition Options

Dirichlet Type
Aqueous Concentration

Dirichlet-type boundary where the aqueous concentration, in terms of solute per mass of aqueous phase, of the solute is specified at the boundary surface centroid. Solutes migrate across the boundary surface via diffusion through the aqueous phase or advection with the aqueous phase.

Inflow Aqueous Concentration

Dirichlet-type boundary where the aqueous concentration, in terms of solute per mass of aqueous phase, of the solute is specified at the boundary surface centroid. Solutes migrate across the boundary surface only via aqueous phase advection in the direction opposite of the boundary-surface normal (i.e., into the domain).

Inflow-Outflow Aqueous Concentration

Dirichlet-type boundary where the aqueous concentration, in terms of solute per mass of aqueous phase, of the solute is specified at the boundary surface centroid.  Solutes migrate across the boundary surface only via aqueous phase advection.  When aqueous phase flux is in the direction of the boundary-surface normal (i.e., out of the domain), the solute concentration is that at the grid-cell centroid, and when aqueous phase flux is in the direction opposite of the boundary-surface normal (i.e., into the domain), the solute concentration is that specified via the input.

Volumetric Concentration

Dirichlet-type boundary where the aqueous concentration, in terms of solute per volume of aqueous phase, of the solute is specified at the boundary surface centroid.  Solutes migrate across the boundary surface via diffusion through the aqueous phase or advection with the aqueous phase..

Inflow Volumetric Concentration

Dirichlet-type boundary where the solute concentration, in terms of solute per grid-cell volume, of the solute is specified at the boundary surface centroid.  Solutes migrate across the boundary surface only via aqueous phase advection in the direction opposite of the boundary-surface normal (i.e., into the domain).

Inflow-Outflow Volumetric Concentration

Dirichlet-type boundary where the solute concentration, in terms of solute per grid-cell volume, of the solute is specified at the boundary surface centroid.  Solutes migrate across the boundary surface only via aqueous phase advection. When aqueous phase flux is in the direction of the boundary-surface normal (i.e., out of the domain), the solute concentration is that at the grid-cell centroid, and when aqueous phase flux is in the direction opposite of the boundary-surface normal (i.e., into the domain), the solute concentration is that specified via the input.

Initial Condition

Dirichlet-type boundary where the aqueous concentration is set to be equal to the aqueous concentration at the grid-cell centroid at the start of the simulation. Solutes migrate across the boundary surface via diffusion through the aqueous phase or advection with the aqueous phase.

Outflow

Solutes migrate across the boundary surface only via aqueous phase advection in the direction of the boundary-surface normal (i.e., out of the domain). As diffusive transport across the boundary surface is not considered the solute concentration entry is ignored.

Neumann Type
Zero Flux

Neumann-type boundary where the solute flux is zero.  Solutes are prevented from crossing the boundary surface regardless of their aqueous-phase concentration.

Species Boundary Condition Options 

Dirichlet Type
Aqueous Concentration

Dirichlet-type boundary where the aqueous concentration, in terms of species per mass of aqueous phase, of the species is specified at the boundary surface centroid.  Species migrate across the boundary surface via diffusion through the aqueous phase or advection with the aqueous phase.

Inflow Aqueous Concentration

Dirichlet-type boundary where the aqueous concentration, in terms of species per mass of aqueous phase, of the species is specified at the boundary surface centroid.  Species migrate across the boundary surface only via aqueous phase advection in the direction opposite of the boundary-surface normal (i.e., into the domain).

Inflow-Outflow Aqueous Concentration

Dirichlet-type boundary where the aqueous concentration, in terms of species per mass of aqueous phase, of the species is specified at the boundary surface centroid.  Species migrate across the boundary surface only via aqueous phase advection.  When aqueous phase flux is in the direction of the boundary-surface normal (i.e., out of the domain), the species concentration is that at the grid-cell centroid, and when aqueous phase flux is in the direction opposite of the boundary-surface normal (i.e., into the domain), the species concentration is that specified via the input.

 
Outflow

Species migrate across the boundary surface only via aqueous phase advection in the direction of the boundary-surface normal (i.e., out of the domain). As diffusive transport across the boundary surface is not considered the species concentration entry is ignored.

Initial Condition

Dirichlet-type boundary where the aqueous concentration is set to be equal to the aqueous concentration at the grid-cell centroid at the start of the simulation. Species migrate across the boundary surface via diffusion through the aqueous phase or advection with the aqueous phase.


Warning

ECKEChem 

Species are reactive and require the ECKEChem Module to be implemented in the simulator.
Neumann Type
Zero Flux

Neumann-type boundary where the species flux is zero.  Species are prevented from crossing the boundary surface regardless of their aqueous-phase concentration.

 
Warning

ECKEChem 

Species are reactive and require the ECKEChem Module to be implemented in the simulator.

STOMP User Guide Home

Modules