Porosity

Porosity may be fixed at a constant value, or it may be a function of the effective porepressure, the volumetric strain, the temperature and/or chemical precipitates

Available porosity formulations include:

Constant: PorousFlowPorosityConst

The simplest case where porosity does not change during the simulation. A single value of porosity can be used, or a spatially varying porosity AuxVariable can be used to define a heterogeneous porosity distribution.

Linear: PorousFlowPorosityLinear

In this case

(1)

Here the "ref" values are reference values, is the effective fluid pressure, is the temperature and is the total volumetric strain. These may be spatially-varying AuxVariables to define a heterogeneous porosity distribution. , and are real-valued coefficients. A lower-bound on porosity, , may also be defined.

Exponential: PorousFlowPorosity

Using PorousFlowPorosity with appropriately set flags, porosity can depend on:

  • total strain, with mechanical = true

  • effective porepressure, with fluid = true

  • temperature, with thermal = true

  • precipitated minerals, with chemical = true

A combination of these may be used, to simulate, for instance, THM or HM coupling.

Theoretical evolution of porosity leading to the exponential form

The evolution of the porosity is governed by (Detournay and Cheng, 1993; Chen and Zhang, 2009) (2) Here is the bulk modulus of the drained porous skeleton: . Also is the volumetric coefficient of thermal expansion. denotes the total precipitated mineral concentration: where is the user-defined weight for each mineral and is the concentration (m/m) of precipitated mineral.

Special form for fracture flow

It is sometimes convenient to extend Eq. (2) to include a new Biot coefficient, : (3) This form is useful when modelling flow through fractures, where the porosity increases with porepressure, which can be simulated with .

Solution

The equation for porosity has solution (4) where is the porosity at reference porepressure, zero elastic strain, reference temperature and reference mineral concentration. Note this porosity can become negative, and an option for ensuring positivity is detailed below.

Time-lag in mineralisation problems

With mineralisation, now depends on total mineral concentration, . However, the evolution of is governed by , where is a reaction rate which is independent of porosity (but dependent on the primary chemical species, temperature, etc) and is the aqueous phase saturation. Therefore a circular dependency exists: depends on , and depends on . This could be broken by promoting porosity to a MOOSE Variable, and solving for it. Instead, PorousFlow replaces in Eq. (4) by the approximate form Note that the old value of porosity is used on the right-hand-side, which breaks the cyclic dependency problem.

Thermal-mechanical simulations

Without porepressure and mineralisation effects, the correct expression for porosity as a function of volumetric strain and temperature is (5)

These expressions may be modified to include the effects of plasticity.

Preventing negative porosity

The exponential expressions Eq. (4) and Eq. (5) can yield negative porosity values, which are unphysical. To ensure positivity of , PorousFlow offers the following option. First write both equations in the form In deriving specific forms for , and , above, it has been assumed that is small. Since physically , for the porosity will be physically meaningful, but for , can become negative. For example, for positive volumetric strain () the porosity is always physically meaningful. However, when the porous material is squashed with negative and large volumetric strain () the porosity can become negative.

Define Then, for , the above expression for is replaced by At first glance this expression may appear rather obscure. It has been constructed to satisfy the following requirements:

  • As , , which is necessary for continuity at .

  • As , , which is physically desirable.

  • The expression is monotonically decreasing, which is both physically desirable and computationally desirable (otherwise there may be non-unique solutions in a PorousFlow simulation).

  • Finally, at , its derivative is , as desired from continuity of the derivative at .

Evolution of porosity for an isothermal, mineral-free, situation

The evolution of porosity is fundamental to the coupling between flow and solid mechanics. Consider the isothermal situation with no plasticity. The following presentation is mostly drawn from (Detournay and Cheng, 1993).

Denote the change of a quantity, q, by . Recall that the porosity is defined by , where is an arbitrary volume of the porous material, and is the porevolume within that volume. Also, by definition of the effective stress, Taking the trace of this equation, and using yields In most instances it is appropriate to write this equation as (6) (which is Detournay and Cheng Eqn(20a)) where the total mechanical pressure is and is the so-called drained bulk modulus . To find the evolution equation for porosity, a similar equation for must be derived.

Assuming linearity (7) (This is Detournay and Cheng Eqn(20b).) The Betti-Maxwell reciprocal theorem yields and , as is now shown.

The work increment is So during some deformation that takes from to , and from to , the total work is

Now consider two experiments:

  • First take from to with fixed at . Then, leaving fixed, take from to . The first takes work , while the second takes work .

  • First take from to with fixed at . Then, leaving fixed, take from to . The first takes work , and the second takes work .

The two experiments must give the same work done (this is called the Betti-Maxwell reciprocal theorem), which yields (8) (This is Detournay and Cheng Eqn(22).)

Now to identify . Consider a so-called ideal porous material, which is characterised by a fully-connected pore space and a homogeneous and isotropic matrix material. In this case, applying a uniform porepressure, , and an equal mechanical pressure, , the solid material will experience a uniform pressure throughout its skeleton. This means it will deform uniformly without any shape change, and Substituting this equation, this specific pressure condition, and Eq. (8) into Eq. (6) and Eq. (7), yields (This is Detournay and Cheng Eqns(22), (24b), (37) and the definition of in Table 1.)

Now that and have been identified, they may be substituted into Eq. (7). Rearranging yields Using the expression for yields Now , so using the definition of yields as written in Eq. (2).

Evolution of porosity due to mineral precipitation

Consider now the case without any fluid porepressure or temperature, but with mineral precipitation.

The concentration of a mineral species, is defined to be , where the numerator is the volume of the mineral species and the denominator is the reference volume. Suppose that and the definition of porosity is .

Now consider volume changes via (where is the volumetric strain) and changes in the volume of mineral, . Under these changes, the porosity evolves as This yields, after a little algebra, as written in Eq. (2).

References

  1. Z. Chen and Y. Zhang. Well flow models for various numerical methods. International Journal of Numerical Analysis and Modelling, 6:375–388, 2009.[BibTeX]
  2. E. Detournay and A. Cheng. “Fundamentals of poroelasticity,” Chapter 5 in Comprehensive Rock Engineering: Principles, Practice and Projects, Vol II, Analysis and Design Method. Also available online. Pergamon Press, New York, 1993.[BibTeX]