MOMOP

MOMOP describes additional options as shown in Table 17.

Free format for 26 parameters. Or Format(26I1)

MOP2(I), I = 1,26

In TOUGH4, input with keyword "MOMOP" becomes obsolete. Its Input can be replaced by input of IE array through keyword "SELEC". if any element in IE(101)-IE(126) is not 0, it will replace corresponding element in MOP2(1)-MOP2(26). Therefore, using MOMOP is not encouraged in TOUGH4.

Table 17. More Options

MOP2

Value

Description

1

0, 1

2

Minimum number of Newton-Raphson iterations

Allow convergence in a single Newton-Raphson iteration

Perform at least two iterations; primary variables are always updated.

2

Void

3

Void

4

0

>0

Vapor pressure reduction

No vapor pressure reduction at low liquid saturation

Reduces vapor pressure for Sl < 0.02 to prevent liquid disappearance by evaporation (only certain EOS modules)

5

0

>0

Active Fracture Model (Appendix C)

Active Fracture Model applied to liquid phase only

Active Fracture Model applied to all phases

6

0

>0

Leverett scaling of capillary pressure

No Leverett scaling

Rescale capillary pressure: if element-specific permeabilities are specified

7

Void

8

Void

9

0

>0

Time stepping after time-step reduction to honor printout time

Continue with time step imposed by forced time-step reduction

Continue with time step used before forced time-step reduction (may conflict with time-step size for time-dependent honoring source/sink or first-type boundary.)

10

0

1

2

Writing SAVE file

Write SAVE file only at the end of a forward run

Write SAVE file after each printout time

Write separate SAVE files after each printout time

11

0

1 2

Water properties

International Formulation Committee (1967)

IAPWS-IF97

IAPWS-95

12

0

>0

Enthalpy of liquid water

Potential energy not included in enthalpy of liquid water

Potential energy included in enthalpy of all phases (Stauffer et al., 2014)

13

0

>0

Adjustment of Newton-Raphson increment weighting

No adjustment

Reduce WNR by MOP2(13) percent if Newton-Raphson iterations oscillate and time step is reduced because ITER=NOITE

14

Void

15

0

>0

Porosity used for calculation of rock energy content

Use porosity of block ROCKS; this assumes that the porosities provided in block INCON were the result of a pore compressibility/expansivity calculation; the “original” porosity from block ROCKS is used to compensate for equivalent rock-grain density changes.

Use porosity from block INCON; this assumes that these porosities were not the result from a pore compressibility/expansivity calculation; changes in rock-grain density due to pore compressibility/expansivity are not compensated.

16

Void

17

0

1

Variables printed on FOFT files (time series)

Print variables according to input data in block FOFT

Print variables according to input data in block OUTPU

18

0

1

Adjustment of time step size to coincide with transient boundary times.

No adjustment

Do the adjustment

19

0

1

Use anisotropic permeability values kx, ky and kz of rocks to construct

the diagonal of the permeability tensor.

No

Yes

20

Void

21

0

>0

Honoring generation times in blocks TIMBC.

Time stepping ignores times specified in block TIMBC

Time steps adjusted to match times in block TIMBC

22

Void

23

Void

24

Void

25

0

1

Check mesh for duplicate element names

Do not check the mesh

Check the mesh

26

0

1

2

3

Printout format

Both traditional output and CSV formats

CSV/TECPLOT format (types separated; times combined)

CSV/TECPLOT format (types separated; times separated)

CSV/TECPLOT format (types combined; times combined)

Used in: All EOS modules

Last updated