Material for unsaturated flow in lattice models - LatticeTransMat

Lattice transport elements have to be used with this material model. A positive sign is assumed for liquid tension, unlike the convention of soil mechanics which assumes compression positive.

These transport elements are idealised as one-dimensional conductive pipes. The gradient of hydraulic head, which governs flow rate along each transport element, is determined from the capillary pressures at the two nodes.

The mass balance equation describes the change in moisture inside a porous element as a consequence of liquid flow and solid-liquid retention. It leads to the following partial differential equation

div (265)

where is the capillary pressure, $ c$ is the mass capacity function( ), $ k$ is the Darcy hydraulic conductivity (), is the fluid mass density, is the acceleration of gravity, is the capillary height and t is the time.

The hydraulic conductivity $ k$ consists of


where is the hydraulic conductivity of the intact material and is the additional conductivity due to cracking.

Darcy hydraulic conductivity is defined as


where is the dynamic viscosity (), $ \kappa$ is the permeability also called intrinsic conductivity(), and is the relative conductivity. is a function of the effective degree of saturation.

The cracking part is


where $ \xi$ is a tortuosity factor taking into account the roughness of the crack surface, is the equivalent crack opening of the dual mechanical lattice and $ h$ is the length of the dual mechanical element.

Borek Patzak