Solution Control Card Options (EOR-BO)
For STOMP-EOR-BO, four classes of options are specified via the Solution Control Card: 1) Execution Mode Options, 2) Operational Mode Options, 3) Fluid Diffusion Options, and 4) Interfacial Averaging Options. Execution Mode Options are used to specify the state of the simulation. Operational Mode Options are used to specify the solved equations and active processes. Fluid Diffusion Options are used to specify the models for component diffusion through fluids. Interfacial Averaging Options are used to specify the model for computing state variables at the centroids of grid-cell surfaces.
Execution Mode Options
Two Execution Modes are recognized: Normal, and Restart. In the Normal mode, initial state conditions are declared through the Initial Conditions Card. In the Restart mode, initial state conditions are assigned via a restart file from a previous execution or declared through the Initial Conditions Card, using the special overwrite option for selected parameters. Unless specified through the Output Control Card, restart files (i.e., restart.n) are generated at each plot.n write event, and have name extensions that correspond to the generating time step (e.g., the file restart.0028 would have been generated at the conclusion of time step 28). Restart files are text files that contain simulation time and control information, and a collection of field variables needed to redefine the simulation state for the operational mode.
Normal Mode
In the Normal mode, STOMP-W executes from a declared start time, using an initial state declared through the Initial Conditions Card, until the declared stop time, the declared number of time steps, an execution error, or a sequence of convergence failures. The following keywords may also be used after specifying the "Normal" Execution Mode.
No Flow
This option results in the coupled flow and transport equations only being computed once, eliminating the flow calculations each time step for a reactive transport problem with a steady flow field.
No Flow Option
The No Flow option, when used in conjunction with the Normal or Restart execution modes, results in the coupled flow and transport equations only being computed once. This option can be used to eliminate the flow calculations each time step for a solute transport problem with a steady flow field.
Scaling Factor
Inverse (UCode)
Restart Mode
In the Restart mode, STOMP-W executes from either a declared start time or the start time specified in the restart file, using an initial state defined by a previous execution, until the declared stop time, the declared number of time steps, an execution error, or a sequence of convergence failures. The following keywords may be after specifying the "Restart" Execution Mode.
File
This option triggers STOMP-W to read an additional character string, which is the name of the restart file.
No Flow
This option results in the coupled flow and transport equations only being computed once, eliminating the flow calculations each time step for a reactive transport problem with a steady flow field.
No Flow Option
The No Flow option, when used in conjunction with the Normal or Restart execution modes, results in the coupled flow and transport equations only being computed once. This option can be used to eliminate the flow calculations each time step for a solute transport problem with a steady flow field.
Scaling Factor
Inverse (UCode)
Overwrite
When the keyword "overwrite" is included in the initial or boundary conditions cards with any of the above options during a restart simulation, the specified values will overwrite those from the restart file.
Operational Mode Modifiers
The Operational Mode is EOR. This identifier is used to make certain that the operational mode of the STOMP executable matches the operational mode declared in the input file. The activation of solute transport, and reactive species transport is specified via keyword modifiers to the operational mode. Models for transporting passive solutes and reactive species is additionally controlled via keyword modifiers to the operational mode. For example, solute transport is solved using the Patankar method, unless the keywords TVD or Roe Superbee also appear.
w/Transport
This modifier activates solute transport.
w/ECKEChem
This modifier activates reactive species transport.
ECKEChem
Reactive transport is specified by including the keyword ECKEChem in the Operational Mode field. The reactive transport algorithms use the same transport schemes as the solute transport model, and therefore are controlled through the keyword options TVD and Roe Superbee.
ECKEChem
Using this keyword requires the ECKEChem Module to be implemented in the simulator.
Leonard-TVD Solute Transport
If either TVD or Leonard keywords appear with the Eckechem keyword, then the Leonard-TVD transport scheme is implemented. This third-order scheme using a TVD technique (Datta Gupta et al. 1991) is most appropriate for advection-dominated flow (high Peclet numbers). Conventional techniques, like the one discussed by Patankar (1980), suffer from artificial diffusion that smears otherwise sharp fronts. The smearing is a result of the first-order approximation of the advective term in the transport equation. Datta Gupta et al. (1991) proposed and successfully tested a third-order differencing scheme with an appropriate flux limiting function which significantly minimizes numerical diffusion while, at the same time, avoiding oscillations that commonly affect classical higher-order schemes.
Roe Superbee
First-Order Upwind
Flow Velocity or Flux Velocity for Hydraulic Dispersion
flow-hyd or flux-hyd hydraulic dispersion coefficient
Guess
Guess species concentration on first call.
Porosity Alteration with Precipitation
Effective Reaction Area
keywords - Area scales with water saturation
area + constant surface = w/ Mineral surface area held constant at initial value'
Fluid Diffusion Options
Models for diffusion coefficients for components in the three mobile fluids: 1) aqueous, 2) nonaqueous-liquid, and 3) gas are specified via the Fluid Diffusion Options
Interfacial Averaging Options
State variables at the centroids of grid-cell surfaces are required to compute fluxes between grid-cell centroids. Models for computing the state variables on grid-cell surfaces are referred to as interfacial averaging schemes and use the state variables at adjacent grid-cells to compute the surface variables. Default interfacial averaging schemes or various variable types have been selected for STOMP. Schemes other than the defaults can be specified.
References
[1] Peng, D-Y, and DB Robinson. 1976. “A new two-constant equation of state.” Ind. Eng. Chem. Fundam., 15:59-64.
[2] Redlich, O and JNS. Kowg. 1949. “On the thermodynamics of solutions.” Chem. Rev., 44:233-244.
[3] Soave, G. 1972. “Equilibrium constants from a modified Redlich-Kwong equation of state.” Chem. Eng. Sci., 27:1197-1203.
[4] Peneloux, A, E Rauzy, and R Fréze. 1982. “A consistent correction of Redlich-Kwong-Soave volumes.” Fluid Phase Equilibria, 8:7-23.
[5] Michelsen, ML and JM Mollerup. 2007. Thermodynamic Models: Fundamentals & Computational Aspects, Second Edition, Tie-Line Publications, Holte, Denmark, ISBN 87-989961-3-4.
[6] Wilson, GM 1969. “A modified Redlich-Kwong equation of state. Application to general physical data calculation.” Paper No. 15C presented at the 1969 AIChE 65th National Meeting. Cleveland, OH, March 4-7, 1969.
STOMP User Guide Home
- Simulation Title Card
- Solution Control Card
- Grid Card
- Inactive Nodes Card
- Rock/Soil Zonation Card
- Mechanical Properties Card
- Hydraulic Properties Card
- Saturation Function Card
- 3 Phase Relative Permeability Card
- Salt Transport Card
- Black Oil Properties Card
- Solute/Fluid Interaction Card
- Solute/Porous Media Interaction Card
- Initial Conditions Card
- Boundary ConditionsCard
- SourceCard
- Coupled Well Card
- Output Control Card
- Surface Flux Card
ECKEChem Input Cards