FORCH

FORCH Invokes flow calculation according to the Forchheimer equation; introduces information on non-Darcy flow coefficient model and weighting scheme.

Record FORCH.1

Free format for 7 parameters, or FORMAT (2I5, 6E10.4)

IFORCH, MFORCH, (FORCH(I), I=1,6)

IFORCH Integer parameter to choose type of non-Darcy flow coefficient model (see Table x, in Non-Darcy Flow section). If a negative number is given, the Forchheimer equation is applied to gas flow only; liquid flow is governed by Darcy’s law.

MFORCH Determines interface weighting scheme for non-Darcy flow coefficient β :

MFORCH = 0: Upstream weighting (recommended)

MFORCH = 1: Upstream weighting with WUP (see record PARAM.3)

MFORCH = 2: Harmonic weighting

MFORCH = 3: Geometric average

MFORCH = 4: Arithmetic average

FORCH(I) I = 1, …, 6; user-specified parameters for non-Darcy flow coefficient model (see Table 16; only needed for models IFORCH = 1 and 14).

Record FORCH.2

Free format or FORMAT (8A10)

ND_Rocks1, ND_Rocks2, ......

ND_Rocks Input a list of rocks in which the flow will be treated as non-Darcy flow. These rocks must be defined in the ROCKS data section. If no rock is specified in this record, the whole flow system will be treated as non-Darcy flow.

Used in: All EOS modules

Example:

FORCH

2, 0 //use the second model for the coefficient calculation; use upstream weighting.

SAND1, SAND2 //flow in rock SAND1 and SAND2 is treated as non-Darcy flow.