Back to Input Deck Cards

Back to SUBSURFACE_FLOW

Back to MODE

WIPP_FLOW

Defines options for the WIPP_FLOW subsurface flow mode.

SIMULATION Options

TIMESTEPPER Options

NEWTON_SOLVER Options

Examples

SIMULATION Options

(under SUBSURFACE_FLOW in SIMULATION PROCESS_MODELS block)

Basic Settings

2D_FLARED_DIRICHLET_BCS <string>

Enables Dirichlet boundary conditions to be specified at cell centers through a list of cells ids and the initial condition. The entries are read from the file named by string, with each entry providing the cell ID and flags for whether the pressure and saturation are held constant.

ALLOW_NEGATIVE_GAS_PRESSURE

Allows the flow solution to converge with a negative gas pressure. The only reason for this is to match what is coded in BRAGFLO. (default = True) It is NOT recommended.

ALPHA_DATASET <string>

Specifies the dataset containing ALPHA values for all grid cells. Only needed for the 2D flared grid scenario.

AUTO_PRESSURE_C_B <float>

Parameter for calculating initial pressures internally.

AUTO_PRESSURE_MATERIAL_IDS <int1> … <intN>

Parameter for calculating initial pressures internally.

AUTO_PRESSURE_PB_0 <float>

Parameter for calculating initial pressures internally.

AUTO_PRESSURE_PB_REF <float>

Parameter for calculating initial pressures internally.

AUTO_PRESSURE_RHO_B0 <float>

Parameter for calculating initial pressures internally.

AUTO_PRESS_SHALLOW_ORIGIN <float> <float> <float>

Parameter for calculating initial pressures internally.

BRAGFLO_RESIDUAL_UNITS

Calculates the residual in units of [kg/m³ ] instead of [kmol/sec].

DEFAULT_ALPHA

Allows the users to run a simulation without providing ALPHA through a dataset.

DIP_ROTATION_ANGLE <float>

Dip rotation angle for calculating ELEVATION internally.

DIP_ROTATION_BASEMENT <float>

Dip rotation basement for calculating ELEVATION internally.

DIP_ROTATION_CEILING <float>

Dip rotation ceiling for calculating ELEVATION internally.

DIP_ROTATION_ORIGIN <float>

Dip rotation origin for calculating ELEVATION internally.

DIP_ROTATION_REGIONS <string1> … <stringN>

Specifies regions for which ELEVATION will be calculated internally.

ELEVATION_DATASET <string>

Specifies a dataset containing elevation value for all grid cells.

GAS_COMPONENT_FORMULA_WEIGHT <float>

Molecular weight of gas [g/mol]. (default = 2.01588 for H_2(g))

HARMONIC_PERMEABILITY_ONLY

Performs a harmonic mean weighting on permeability alone, not the lumped density * permeability / viscosity.

LOGGING_VERBOSITY <integer>

Prints additional convergence logging to the screen when the integer is greater than 0 (for process models that have logging verbosity implemented).

MATCH_BRAGFLO_OUTPUT

Reports convergence information in units that match BRAGFLO for comparison purposes.

NO_CREEP_CLOSURE

Turns off creep closure without removing it from the input file.

NO_FRACTURE

Turns off fracture without removing it from the input file.

NO_GAS_GENERATION

Turns off gas generation without removing it from the input file.

REPLACE_INIT_PARAMS_ON_RESTART

Forces the overwrite of base (0) porosity and permeabilities off which solution dependent parameters may be calculated.

REVERT_PARAMETERS_ON_RESTART

Reverts permeability and porosity back to the values specified in the input file instead of the checkpointed values.

SKIP_RESTART

On a restarted simulation, employs the initial condition for process model instead of the restarted solution.

WITH_RADIOLYSIS

Toggles on radiolysis.

---

Expert Settings

COUNT_UPWIND_DIRECTION_FLIP

Turns on and increments a counter whenever the upwind direction changes. This is only useful for non-parallel execution.

FIX_UPWIND_DIRECTION

Limits the update of the upwind direction (based on pressure gradients) to UPWIND_DIR_UPDATE_FREQUENCY Newton iterations.

UNFIX_UPWIND_DIRECTION

Turns off fixing of upwind direction.

UPWIND_DIR_UPDATE_FREQUENCY

Frequency at which the upwind direction is updated if FIX_UPWIND_DIRECTION is included.

---

Debugging Options

RESIDUAL_TEST

Toggle on printing of residual information at a specific cell. RESIDUAL_TEST_CELL must be defined.

RESIDUAL_TEST_CELL

Cell at which residual information will be printed when RESIDUAL_TEST is present.

JACOBIAN_TEST

Toggles on testing of numerical Jacobian usign full residual evaluation.

JACOBIAN_TEST_RDOF

Residual equation that will be printed for JACOBIAN_TEST (X in dR/dX).

JACOBIAN_TEST_XDOF

Unknown that will be printed for JACOBIAN_TEST (R in dR/dX).

NO_ACCUMULATION

Skip accumulation term calculation.

NO_FLUX

Skip internal flux calculation.

NO_BCFLUX

Skip boundary flux calculation.

NO_FRACTURE

Skip fracture.

NO_CREEP_CLOSURE

Skip creep closure.

NO_GAS_GENERATION

Skip gas generation.

PRINT_RESIDUAL

Print the residual to a file pf_residual.txt at each Newton iteration.

PRINT_SOLUTION

Print the solution to a file pf_solution.txt at each Newton iteration.

PRINT_UPDATE

Print the update to a file pf_update.txt at each Newton iteration.

DEBUG

Toggles on increasing verbose output for debugging.

DEBUG_GAS_GENERATION

Increasingly verbose information for gas generation from pm_wipp_srcsink.

DEBUG_FIRST_ITERATION

Stops the simulation after the first Newton iteration when debugging is toggled on.

DEBUG_OSCILLATORY_BEHAVIOR

Turns on increasingly verbose information for a cell where the residual is oscilating.

DEBUG_TS_UPDATE

Prints dtime(1:2) ramping factors used in updating the time step size.

USE_BRAGFLO_CC

Toggles the use of characteristic curves exactly as coded in BRAGFLO. The code was lifted from BRAGFLO and wrapped for use in PFLOTRAN for debugging purposes.

TIMESTEPPER Options

Basic Settings

GAS_SAT_CHANGE_TS_GOVERNOR <float>

Time step size will be reduced on the subsequent time step when the maximum change in gas saturation is above the governor value. Default = 0.3

GAS_SAT_GOV_SWITCH_ABS_TO_REL <float>

When the gas saturation is greater than this value, a relative change in gas saturation will be used to calculate the new time step size based on the GAS_SAT_CHANGE_TS_GOVERNOR, while absolute change is used below the value. Default = 0.01

LIQ_PRES_CHANGE_TS_GOVERNOR <float>

Time step size will be reduced on the subsequent time step when the maximum change in liquid pressure is above the governor value. Default = 5.e5

MINIMUM_TIMESTEP_SIZE <float>

Value below which the time step size is truncated to MINIMUM_TIMESTEP_SIZE. Default = 8.64e-4 sec.

NEWTON_SOLVER Options

Basic Settings

CONVERGENCE_TEST <string>

BOTH - Both the residual infinity tolerance and absolute/relative change criteria must be met to declare convergence. EITHER - Either the residual infinity tolerance or absolute/relative change criter ia must be met to declare convergence. Default = BOTH

DO_NOT_SCALE_JACOBIAN

Turns off Jacobian scaling.

GAS_RESIDUAL_INFINITY_TOL <float>

Infinity norm tolerance for gas residual equation. Units depend on the units of the residual: [kg/m³] with BRAGFLO_RESIDUAL_UNITS and [kmol/sec] without. Default = 1.e-2

GAS_SAT_THRESH_FORCE_EXTRA_NI <float>

Gas saturation threshold forcing an extra Newton iteration. If the gas saturation is outside the bounds of -GAS_SAT_THRESH_FORCE_EXTRA_NI and 1+GAS_SAT_THRESH_FORCE_EXTRA_NI, an extra Newton iteration will be required. Default = 1.e-3

GAS_SAT_THRESH_FORCE_TS_CUT <float>

Gas saturation threshold forcing a time step cut. If the gas saturation is outside the bounds of -GAS_SAT_THRESH_FORCE_TS_CUT and 1+GAS_SAT_THRESH_FORCE_TS_CUT, the time step will be cut immediately. Default = 0.2

JACOBIAN_PRESSURE_DERIV_SCALE

Scalign factor for the liquid pressure derivative. Default = 1.e7

LIQUID_RESIDUAL_INFINITY_TOL <float>

Infinity norm tolerance for liquid residual equation. Units depend on the units of the residual: [kg/m³ ] with BRAGFLO_RESIDUAL_UNITS and [kmol/sec] without. Default = 1.e-2

MAX_ALLOW_GAS_SAT_CHANGE_TS <float>

Maximum absolute change in gas saturation during a time step. Default = 1

MAX_ALLOW_LIQ_PRES_CHANGE_TS <float>

Maximum absolute change in liquid pressure during a time step. Default = 1e7

MAX_ALLOW_REL_GAS_SAT_CHANGE_NI <float>

Maximum allowable relative change in gas saturation between Newton iterations. Default = 1.e-3

MAX_ALLOW_REL_LIQ_PRES_CHANG_NI <float>

Maximum allowable relative change in liquid pressure between Newton iterations. Default = 1.e-2

MIN_GAS_SATURATION_PERTURBATION <float>

Minimum gas saturation perturbation for derivatives.

MIN_LIQ_PRESSURE_PERTURBATION <float>

Minimum liquid pressure perturbation for derivatives.

MIN_LIQ_PRES_FORCE_TS_CUT <float>

The minimum liquid pressure below which Newton iteration will cease and the time step will be cut. Default = -1.e8

REL_GAS_SATURATION_PERTURBATION <float>

Relative gas saturation perturbation for derivatives.

REL_LIQ_PRESSURE_PERTURBATION <float>

Relative liquid pressure perturbation for derivatives.

SCALE_JACOBIAN

Toggles on scaling of the linear system. The Jacobian pressure derivatives are first right-hand scaled by JACOBIAN_PRESSURE_DERIV_SCALE, and then the Jacobian is left-hand scaled by the reciprocal of the maximum absolute value in each row. Default = TRUE

Examples

SIMULATION
  ...
  PROCESS_MODELS
    SUBSURFACE_FLOW flow
      MODE WIPP_FLOW
      OPTIONS
        GAS_COMPONENT_FORMULA_WEIGHT 2.01588d0 #hardwired
        2D_FLARED_DIRICHLET_BCS
          EXTERNAL_FILE ../dirichlet_bcs.txt
        /
        ALLOW_NEGATIVE_GAS_PRESSURE
        ALPHA_DATASET alpha
        BRAGFLO_RESIDUAL_UNITS
        DIP_ROTATION_ANGLE 1.d0
        DIP_ROTATION_ORIGIN 23495.7d0 0.d0 378.685d0
        DIP_ROTATION_CEILING 779.69d0
        DIP_ROTATION_BASEMENT 178.07d0
        DIP_ROTATION_REGIONS rSHFTU
      /
    /
  /
END
...
SUBSURFACE
  NUMERICAL_METHODS FLOW
    TIMESTEPPER
      LIQ_PRES_CHANGE_TS_GOVERNOR 5.d5       ! PRESNORM
      GAS_SAT_CHANGE_TS_GOVERNOR 0.3d0       ! SATNORM
      GAS_SAT_GOV_SWITCH_ABS_TO_REL 1.d0     ! TSWITCH
    /
    NEWTON_SOLVER
      CONVERGENCE_TEST BOTH                  ! ICONVTEST 1
      SCALE_JACOBIAN                         ! LSCALE
      GAS_RESIDUAL_INFINITY_TOL 1.d-2        ! FTOL_PRES
      GAS_SAT_THRESH_FORCE_EXTRA_NI 1.d-3    ! SATLIMIT
      GAS_SAT_THRESH_FORCE_TS_CUT 0.2d0      ! DSATLIM
      LIQUID_RESIDUAL_INFINITY_TOL 1.d-2     ! FTOL_SAT
      MAX_ALLOW_GAS_SAT_CHANGE_TS 1.d0       ! DSAT_MAX
      MAX_ALLOW_LIQ_PRES_CHANGE_TS 1.d7      ! DPRES_MAX
      MAX_ALLOW_REL_GAS_SAT_CHANGE_NI 1.d-3  ! EPS_SAT
      MAX_ALLOW_REL_LIQ_PRES_CHANG_NI 1.d-2  ! EPS_PRES
      MIN_LIQ_PRES_FORCE_TS_CUT -1.d8        ! DPRELIM
      MINIMUM_TIMESTEP_SIZE 8.64d-4          ! DELTMIN
    /
  /
  ...