Solution Control Card Syntax
~Solution Control Card
Choose One Execution Mode OptionNormal
,Restart (expects filename to be "restart")
,Restart with specified restart filename
,Restart Filename
Choose One Operational Mode OptionIsothermal Mode
+ Modifiers(Optional), One or More Operational Mode modifiers may be specifed after the Operational Mode keyword,,If this option is specified, the following must be provided on the line following the Gas Diffusion Option:
Aqueous Density
,units(kg/m^3),Aqueous Viscosity,units(Pa s),Gas Density,units(kg/m^3),Gas Viscosity,units(Pa s),
,
If this option is specified, the following must be provided on the line following the Gas Diffusion Option, and after the density and viscosity values if the Invariant Fluid Density modifier is specified
Fractional CO2 Solubility Factor
,Additional Reactive Transport Modifiers (Optional)
,,,,,,,,,,,,,,
ECKEChem ModifiersNote that Patankar is the default option for reactive transport, and does not need to be specified, but can be
Non-Isothermal Mode
+ Modifiers(Optional), One or More Operational Mode modifiers may be specifed after the Operational Mode keyword,,,If this option is specified, the following must be provided on the line following the Gas Diffusion Option:
Aqueous Density,units(kg/m^3),Aqueous Viscosity,units(Pa s),Gas Density,units(kg/m^3),Gas Viscosity,units(Pa s),
,
If this option is specified, the following must be provided on the line following the Gas Diffusion Option, and after the density and viscosity values if the Invariant Fluid Density modifier is specified
Fractional CO2 Solubility Factor
,Additional Reactive Transport Modifiers (Optional)
,,,,,,,,,,,,,,ECKEChem ModifiersNote that Patankar is the default option for reactive transport, and does not need to be specified, but can be
Number of Execution Time Periods
,
Repeat for Each Execution Time Period Choose One Time-step Reduction OptionFour Successive Time Step Reductions before Execution Stop
,Unit(s),Final Time ,Unit(s) ,Initial Time Step ,Units ,Maximum Time Step ,Units ,Time Step Acceleration Factor ,Maximum number of Newton-Raphson iterations per time step ,Convergence Criterion Time Step Reductions to Minimum Time Step before Execution Stop
,Unit(s),Final Time ,Unit(s) ,Initial Time Step ,Units ,Maximum Time Step ,Units ,Time Step Acceleration Factor ,Maximum number of Newton-Raphson iterations per time step ,Convergence Criterion Minimum Time Step ,Units ,Time Step Reduction Factor ,Maximum Number of Time Steps,
Choose One Aqueous Diffusion Option,,CO2 Diffusion Coefficient,Unit(s),Salt Diffusion Coefficient,Unit(s),, Choose One Gas Diffusion Option,,Gas Diffusion Coefficient,
,Clay Mass Fraction,
,Number of Interfacial Averaging Variables,
Interfacial Averaging SchemesWhile STOMP allows the user to specify alternative interfacial averaging schemes, it is strongly suggested that novice users in particular specify "0" for the Number of Interfacial Averaging Variables and use the default schemes, as not all schemes are appropriate for all field variables.
For details see
Repeat for Number of Interfacial Averaging VariablesSurface Variable Option,Interfacial Averaging Scheme,
Execution Time PeriodsExecution periods refer to a period of simulation time. STOMP allows the user to specify a single or multiple execution periods. For each execution period, the user can control the initial time step, maximum time step, time step acceleration factor, maximum number of Newton-Raphson iterations, and convergence criterion independently.
If no initial time record is read for a restart simulation, the initial time record is obtained from the restart file.
Two additional solution controls can be specified for each Execution Period input line: 1) the minimum time step, and 2) the time-step cut factor. These optional solution controls must be specified together. The minimum-time-step parameter sets the minimum time step for an execution period. Without this control, STOMP quits after four time-step reductions. Cyclic injection well schedules can yield simulations that execute with large time steps during steady injection periods, but require small time steps during the injection startup. The minimum-time-step parameter allows the code to cut the time step to the minimum value before quitting. The time-step cut factor is set to 0.2 by default, which means that time steps are cut to 20% of their value when a convergence failure occurs. Execution performance can be improved for some simulations by altering this factor.
Recommended InputsRecommended values for the Time Step Acceleration Factor, Maximum Number of Newton-Raphson Iterations, and Convergence Criterion are 1.25, 8, and 1.e-06, respectively. Except under special circumstances, it is not recommended to change the value for the Convergence Criterion from its recommended value. This value has proven through numerous applications to achieve a good balance between accuracy and execution speed.
Zero Time Step SimulationsBy specifying 0 as the maximum number of time steps, initial conditions can be established for a set of conditions. This can be useful for calculating properties of a fluid that can be then used to determine the pressure gradient to be specified in the initial conditions card.