Theory for Moisture Transport in Porous Media
In a porous material, moisture may be present under vapor and liquid state, and moisture transport occurs through convection and diffusion, liquid water being fixed on the pores walls mainly in two ways:
The predominance of adsorption or capillarity depends both on the amount of liquid water present in the porous medium, and on the size and shape of the pores, and should be defined in function of the material and process under study. Two moisture regions can be distinguished:
The relative humidity and capillary pressure are called state variables: they describe the thermodynamic state of water in the pore space under two phases in equilibrium. At equilibrium, they are related through Kelvin’s law, which is a relation of thermodynamical origin. As long as equilibrium is satisfied, the state variables are continuous between two media (at the interface between two porous media with different porosities, or between a porous and a free medium), and Kelvin’s law is valid on both sides of the interface. In opposition, the liquid saturation or the moisture content depend on the porous matrix.
Note that capillary transport of liquid water should be considered in both regions.
The Moisture Transport in Porous Media Interface is dedicated to the modeling of porous media in the hygroscopic moisture region. The total moisture content is thus defined by a sorption curve as
where the liquid saturation, sl, describes the amount of liquid water within the pores, and εp (dimensionless) is the porosity.
By summing the mass conservation equations for liquid water and vapor, and equilibrium between the liquid and gas phases, a single equation for the transport of the total moisture can be written, in which the evaporation and condensation source terms of each equation cancel out:
(4-137)
with the following material properties, fields, and source:
ww) (SI unit: kg/m3) is the total moisture content.
  ϕw (dimensionless) is the relative humidity.
ρg (SI unit: kg/m3) is the moist air density, defined from the dry air and vapor densities, in function of the amount of vapor.
ug (SI unit: m/s) is the moist air velocity field, that should be interpreted as the Darcy velocity, that is, the volume flow rate per unit cross sectional area.
ωv (dimensionless) is the vapor mass fraction in moist air, defined as
Μv (SI unit: kg/mol) is the molar mass of water vapor.
csat (SI unit: mol/m3) is the vapor saturation concentration.
T (SI unit: K) is the temperature.
gw is the moisture diffusive flux is defined as:
Deff (SI unit: m2/s) is the effective vapor diffusion coefficient in the porous medium, computed from the diffusion coefficient in a free medium, and accounting for the porosity and tortuosity of the porous medium.
ul (SI unit: m/s) is the liquid water velocity field, defined from the absolute pressure gradient by the Darcy’s law as:
κrl (dimensionless) is the relative liquid water permeability, that may be defined as a function of the liquid saturation.
κ (SI unit: m2) is the porous medium permeability.
pA (SI unit: Pa) is the absolute pressure.
ρl (SI unit: kg/m3) is the liquid water density, defined as a function of the temperature.
glc is the liquid water capillary flux, defined either from the capillary pressure gradient by a Darcy’s law:
or by a diffusion equation using the relative humidity:
Dw (SI unit: m2/s) is the moisture diffusivity.
pc is the capillary pressure, that can be related to the relative humidity using Kelvin’s law:
G (SI unit: kg/(m3s)) is a moisture source (or sink). See the Moisture Source node.
The moisture source due to evaporation may be obtained from the transport equation for vapor:
This moisture source is added as a mass source in the fluid flow equations solved for the computation of the moist air velocity field, ug.
The corresponding latent heat source is defined as:
where Lv is the latent heat of evaporation.
For a steady-state problem, the relative humidity does not change with time and the first term in the moisture transport equation disappears.