For eSTOMP-W, three classes of options are specified via the Solution Control Card:
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.28 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.
In the Normal mode, eSTOMP executes from a declared start time, using an initial state declared through the Initial Conditions Card, until the declared stop time or the declared number of time steps is reached, or an execution error or a sequence of convergence failures occur. 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 being computed only once, eliminating the flow calculations at each time step for a solute or reactive transport problem with a steady flow field. | |
with PETSc |
This option allows the user to specify PETSc absolute and relative tolerances. | |
with inverse |
This option triggers the generation of UCODE output files and allows for the user to specify the Observed Data Card to compare eSTOMP results with observed data. |
The “W” in eSTOMP-W signifies the water operational mode of the simulator, which solves a mass balance equation for water. This operational mode is therefore specified using the keyword water. Solution of transport equations for passive solutes or reactive species can also be specified by adding the keywords solute, e.g., “water with solute,” or reactions, e.g., ”water with reactions,”. Modeling of water uptake by plant roots can enabled by adding the keyword plants, e.g., “water with plants”.
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 eSTOMP. Schemes other than the defaults can be specified.
Datta-Gupta, A, LW Lake, GA Pope, and K Sepehnoori. 1991. High-Resolution Monotonic Schemes for Reservoir Fluid Flow Simulation. In Situ, 15(3):289-317.
Patankar, SV. 1980. Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing Corporation, Washington, D. C.
Poeter, EP and MC Hill. 1998. Documentation of Ucode: A Computer Code for Universal Inverse Modeling, 98-4080, USGS Water-Resources Investigations Report.