On A Saturated Poromechanical Framework and Its Relation to Abaqus Soil Mechanics and Biot Poroelasticity Frameworks

Lei Jin

2023-04-05

Abstract:We introduce a conservational and constitutive framework for a closed and isothermal two-phase material system consisting of a deformable porous solid matrix and a fully saturating, single-phase, and compressible pore fluid without inter-phase mass exchange. We re-derive a generalized fluid mass balance law using fundamental transport rules. We also summarize from the literature a fundamental force balance law for the fluid-solid mixture that does not require any effective stress law a priori. We show that the two conservation laws are coupled naturally to second-order without any constitutive prerequisites. This differs from Biot poroelasticity, which first postulates first-order fluid-solid coupling as two linearized constitutive relationships and then enforces them into simple Eulerian form of conservation laws. Next, we examine a limiting-case unsaturated soil mechanics framework implemented in Abaqus, by assuming isothermal conditions, full saturation, and no adsorption, and then relate it to our framework. We prove that (1) the two mass balance laws are always equivalent regardless of fluid constitutive behaviors, and (2) the two force balance laws are equivalent in their specific forms with a linearly elastic solid skeleton. Finally, taking advantage of a fundamental pore constitutive law, we show how our framework, and by extension the limiting-case Abaqus framework, naturally gives rise to the distinction between drained and undrained settings, and reduces to Biot poroelasticity under simplifying conditions. Notably, our framework indicates the presence of an additional solid-to-fluid coupling term when the solid particle velocity is non-orthogonal to the Darcy velocity.

Fluid Dynamics,Geophysics

What problem does this paper attempt to address?

The problems that this paper attempts to solve are: in a multiphase material system, how to establish a general mechanical framework for saturated porous media without a priori assumptions, and explore its relationship with Abaqus soil mechanics and Biot poroelasticity theory. Specifically:

1. **Natural Emergence of Fluid - Solid Coupling**:
- The author re - derives the law of fluid mass conservation and summarizes the basic force - balance laws for fluid - solid mixtures in the literature. These laws do not require any a priori effective - stress laws.
- The author proves that these two conservation laws are naturally coupled at the second order without any constitutive pre - conditions. This is different from Biot poroelasticity theory, in which first - order fluid - solid coupling is assumed as two linearized constitutive relations and then they are substituted into the conservation laws in simple Euler form.

2. **Relationship with the Abaqus Soil Mechanics Framework**:
- The author examines the limiting case (i.e., isothermal conditions, fully saturated and without adsorption) of an unsaturated soil mechanics framework implemented in Abaqus and compares it with the proposed framework.
- The author proves that under these two frameworks, the law of mass conservation is always equivalent, regardless of the constitutive behavior of the fluid; and in a specific form, the force - balance laws are also equivalent, especially in the case of a linearly elastic solid skeleton.

3. **Differences between Drained and Undrained Conditions**:
- Using the basic pore constitutive laws, the author shows how their framework (and extended to the Abaqus framework) naturally leads to the distinction between drained and undrained conditions and reduces to Biot poroelasticity theory under simplified conditions.

4. **Existence of the Solid - Liquid Coupling Term**:
- The author points out that when the solid - particle velocity is non - orthogonal to the Darcy velocity, there is an additional solid - liquid coupling term.

By solving these problems, the author aims to help researchers from different fields understand the common problems behind these frameworks, promote clear understanding, and facilitate the appropriate use of commercial tools.

### Summary of Mathematical Formulas

1. **Law of Fluid Mass Conservation**:
\[
\frac{\partial (\rho_f \phi_f)}{\partial t}+\nabla\cdot(\rho_f \phi_f \mathbf{v}_f) = 0
\]

2. **Force - Balance Law**:
\[
\nabla\cdot(\sigma' + p_f \mathbf{I})+\rho_f \phi_f \mathbf{g}+\rho_s(1 - \phi_f)\mathbf{g} = 0
\]
where \(\sigma'\) is the effective Cauchy stress tensor, \(p_f\) is the pore - fluid pressure, \(\mathbf{I}\) is the unit tensor, and \(\mathbf{g}\) is the gravitational acceleration.

3. **Effective - Stress Law**:
\[
\sigma'=\sigma-\alpha p_f \mathbf{I}
\]
where \(\alpha = 1-\frac{K}{K_s}\) is the Biot - Willis coefficient, and \(K\) and \(K_s\) are the bulk moduli of the porous solid skeleton and solid particles, respectively.

4. **Darcy's Law**:
\[
\mathbf{v}_f = -\frac{k}{\eta}(\nabla p_f-\rho_f \mathbf{g})
\]
where \(\mathbf{v}_f\) is the Darcy velocity, \(k\) is the permeability tensor, and \(\eta\) is the fluid viscosity.

5. **Hooke's Law (Linearly Elastic Solid Skeleton)**:
\[
\sigma'=\mathbb{C}_e:\varepsilon
\]
where \(\mathbb{C}_e\) is the fourth - order elastic stiffness tensor.