Sorption

Sorption:

This page discusses:

Sorption

A porous medium becomes partially saturated when the total pore liquid pressure, uw, becomes negative (see Effective stress principle for porous media). Negative values of uw represent capillary effects in the medium. For uw<0 it is known that the saturation lies within certain limits that depend on the value of the capillary pressure, -uw (see Continuity statement for the wetting liquid phase in a porous medium). Typical forms of these limits are shown in Figure 1. We write these limits as sasse, where sa(uw) is the limit at which absorption will occur (so that s˙>0), and se(uw) is the limit at which exsorption will occur (so that s˙<0). The transition between absorption and exsorption and vice versa takes place along “scanning” curves (discussed below). These curves are approximated by the single straight line shown in Figure 1.

Typical absorption and exsorption behaviors.

When partial saturation is included in the analysis of flow through a porous medium, the absorption behavior, the exsorption behavior, and the scanning behavior (between absorption and exsorption) should each be defined. Each of these behaviors is discussed below. If sorption is not defined at all, Abaqus/Standard assumes fully saturated flow (s=1.0) for all values of uw.

Strongly unsymmetric partially saturated flow coupled equations result from the definition of sorption. Therefore, Abaqus/Standard automatically uses its unsymmetric matrix storage and solution scheme (see Defining an Analysis) if you request partially saturated analysis (i.e., if sorption is defined).

Defining Absorption and Exsorption

Absorption and exsorption behaviors are defined by specifying the pore liquid pressure, uw (negative “capillary tension”), as a function of saturation. In most physical cases the wetting liquid cannot be driven to zero saturation; to achieve zero saturation, the data would have to define uw- as s0.0. Absorption and exsorption data can be defined in either a tabular form or an analytical form.

Tabular Form

By default, you define the absorption and exsorption behaviors by specifying uw as a tabular function of s, where 0+s1.0.

Analytical Form

The absorption and exsorption behaviors can be defined by the following analytical form:

uw=1Bln[(s-s0)(1-s0)+A(1-s)]    for    s1s<1,uw=uw|s1-duwds|s1(s1-s)            for    s0s<s1,

where A,B are positive material constants and s0,s1 are parameters used to define the lower bound of the saturation values of interest (see Figure 2).

Logarithmic form of absorption and exsorption behaviors.

Defining the Behavior between Absorption and Exsorption

The behavior between absorption and exsorption is defined by a scanning line of user-specified constant slope, (duw/ds)|s. This slope should be larger than the slope of any segment of the absorption or exsorption behaviors.

If absorption and exsorption behaviors are defined with no scanning line, the slope of the scanning line is taken as 1.05 times the largest value of duw/ds given in the absorption and exsorption behavior definitions.