HydroGeoSphere/Porous Medium

Default porous media transport values
By default, all porous media zones (and elements) in the domain will be assigned default porous media transport properties which are listed in Table 5.21. Included here are parameters for modifying the porous medium so that it acts as a double-porosity medium for simulating transport, as described in Section 2.6.1.3 and also for isotopic fractionation, as described in Section 2.6.1.4.

The following instructions can be applied to porous media, as discussed in Section 5.8.1, to modify the default transport parameters. For each instruction we will indicate its scope (i.e. .grok, .mprops). Recall that if an instruction is used in the prefix.grok file, it will affect the current set of chosen zones, while in a properties (e.g. .mprops) file, it will only affect the named material of which it is a part.

Longitudinal dispersivity
Scope: .grok .mprops
 * 1) val Longitudinal dispersivity [L], $${\alpha}_l$$ in Equation 2.97.
 * &bull; &bull; &bull;

Transverse dispersivity
Scope: .grok .mprops
 * 1) val Horizontal component of the transverse dispersivity [L], $${\alpha}_t$$ in Equation 2.97.
 * &bull; &bull; &bull;

Vertical transverse dispersivity
Scope: .grok .mprops
 * 1) val Vertical component of the transverse dispersivity [L], $${\alpha}_t$$ in Equation 2.97.
 * &bull; &bull; &bull;

Tortuosity
Scope: .grok .mprops
 * 1) val Tortuosity, $${\tau}$$ in Equation 2.97.
 * &bull; &bull; &bull;

Anisotropic tortuosity ratio
<tt>Scope: .grok .mprops</tt>
 * 1) y_tortratio Tortuosity ratio in the y-direction. Default value is 1.
 * 2) z_tortratio Tortuosity ratio in the z-direction. Default value is 1.

By default, tortuosity is isotropic, since the ratio values are set to 1 in both the y- and z-directions. You may make tortuosity anisotropic by entering a value greater than 0 and less than 1. These values will be used to multiple the tortuosity, $${\tau}$$ in the y- and z-directions respectively, to obtain the directional values.
 * &bull; &bull; &bull;

Bulk density
<tt>Scope: .grok .mprops</tt>
 * 1) val Bulk density [M L−3], $${\rho}_b$$ in Equation 2.96. If this instruction is used, the value of the density of solids previously saved for this material is overwritten by the density of solids computed from the bulk density, the density of water and the porosity $$[{\rho}_s=({\rho}_b-{\theta}_s{\rho})/(1-{\theta}_s)]$$.
 * &bull; &bull; &bull;

By default, the porous medium acts a single-porosity medium (i.e. the immobile zone is inactive) because the porosity and mass transfer coefficient are set to zero. In order to activate the double-porosity feature, you can enter non-zero values for these parameters using the following two instructions:

Immobile zone porosity
<tt>Scope: .grok .mprops</tt>
 * 1) val Immobile zone porosity, $${\theta}_{Imm}$$ in Equations 2.101, 2.135 and 2.136.
 * &bull; &bull; &bull;

Immobile zone mass transfer coefficient
<tt>Scope: .grok .mprops</tt>
 * 1) val Immobile zone mass transfer coefficient [T−1], $${\alpha}_{Imm}$$ in Equation 2.134.
 * &bull; &bull; &bull;

Isotope fractionation data...End
<tt>Scope: .mprops</tt>

Causes grok to begin reading a group of isotope fractionation instructions until it encounters an End instruction.

If no further instructions are issued, the default isotopic fractionation parameter values listed in Table 5.21 will be used.
 * &bull; &bull; &bull;

The following three instructions can be used to change these values:

Reverse rate
<tt>Scope: .mprops</tt>
 * 1) val Reverse fractionation rate [L−1], $$k_r$$ in Equation 2.143.
 * &bull; &bull; &bull;

Fractionation factor
<tt>Scope: .mprops</tt>
 * 1) val Fractionation factor, $${\alpha}_r$$ in Equation 2.143.
 * &bull; &bull; &bull;

Rock-water mass ratio
<tt>Scope: .mprops</tt>
 * 1) val Isotopic rock-water mass ratio, $$x_r$$ in Equation 2.102.
 * &bull; &bull; &bull;

The next four instructions can be used to change the thermal properties of the porous medium:

Thermal conductivity of solid
<tt>Scope: .grok .mprops</tt>
 * 1) val Temperature invariant thermal conductivity of the solids [W L−1 K−1]. The bulk thermal conductivity is computed internally from the volume fractions of the solid and liquid phases.
 * &bull; &bull; &bull;

Temperature-dependent thermal conductivity of solid
<tt>Scope: .grok .mprops</tt>
 * 1) k_s1 Thermal conductivity [W L−1 K−1] at temperature $$t_{s1}$$.
 * 2) t_s1 Temperature [&deg;C] at which the thermal conductivity is equal to $$k_{s1}$$.

If that instruction is specified, the thermal conductivity of the solid phase is temperature-dependent. The bulk thermal conductivity is also temperature-dependent and is computed internally from the volume fractions of the solid and liquid phases. It is assumed here that the thermal conductivity of the solids decreases at a constant rate of 1% per 10&deg;C increase in temperature and the relationship between thermal conductivity and temperature is defined with the pair of values $$(k_{s1},t_{s1})$$.
 * &bull; &bull; &bull;

Specific heat capacity of solid
<tt>Scope: .grok .mprops</tt>
 * 1) val Specific heat capacity of the solid phase [J kg−1 K−1]. The default value is 730.0 J kg−1 K−1.
 * &bull; &bull; &bull;

Note that the density of the solid phase of the porous medium is now computed automatically from the bulk density and porosity.