Retention and Permeability Relationships
Richards’ equation appears deceptively similar to the saturated flow equation set out in the Darcy’s Law interface, but it is notoriously nonlinear (Ref. 3). Nonlinearities arise because the material and hydraulic properties θ, Se, Cm, and kr vary for unsaturated conditions (for example, negative pressure) and reach a constant value at saturation (for example, pressure of zero or above). The volume of liquid per porous medium volume, θ, ranges from a small residual value θr to the total porosity θs. Its value is given in a constitutive relation in the model. The effective saturation, Se, amounts to θ normalized to a maximum value of 1. The specific moisture capacity, Cm, describes the change in θ as the solution progresses, the slope on a plot of θ versus pressure (or pressure head). The relative permeability, kr, increases with moisture content and varies from a nominal value to its maximum value at saturation, which reveals that the fluid moves more readily when the porous medium is fully wet.
The Richards’ Equation interface in COMSOL Multiphysics includes three retention models to define θ, Se, Cm, and kr:
The analytic formulas of van Genuchten (Ref. 4) and Brooks and Corey (Ref. 5) are so frequently used that they are synonymous with this variably saturated flow modeling. Posed in terms of pressure head Hp = p/(ρg), the analytic expressions require data for the saturated θs and residual θr liquid volume fractions as well as constants α, n, m, and l, which specify a particular medium type.
The van Genuchten equations define saturation when the fluid pressure is atmospheric (that is, Hp = 0). These equations are
(4-15)
here, the constitutive parameter m is equal to 11 / n. With the Brooks and Corey approach, an air-entry pressure distinguishes saturated (Hp > −1) and unsaturated (Hp < −1) flow so that
(4-16)
Here, the constitutive parameter m is equal to 1 − 1 / n. COMSOL Multiphysics also provides user-defined options for those who want to incorporate experimental data or arbitrary expressions to define these relationships.