[论文解读] A topology-motivated mixed finite element method for dynamic response of porous media

Zahrasadat Lotfian, Mettupalayam V. Sivaselvan|arXiv (Cornell University)|Jun 22, 2015

Numerical methods in engineering参考文献 52被引用 3

一句话总结

本文提出了一种基于拓扑结构的混合有限元方法，用于模拟饱和多孔介质中的动态双相弹性问题，结合Raviart-Thomas元对流体通量和压力（w-p场）进行离散，同时采用标准Galerkin有限元方法对固体位移（u场）进行离散。该方法通过三场弱形式实现稳定性，并确保有效应力与孔隙压力之间的精确耦合，在Biot-Willis系数α=1时，对位移和流体变量均实现了最优收敛率。

ABSTRACT

In this paper, we propose a numerical method for computing solutions to Biot's fully dynamic model of incompressible saturated porous media [Biot;1956]. Our spatial discretization scheme is based on the three-field formulation (u-w-p) and the coupling of a lowest order Raviart-Thomas mixed element [Raviart,Thomas;1977] for fluid variable fields (w, p ) and a nodal Galerkin finite element for skeleton variable field (u). These mixed spaces are constructed based on the natural topology of the variables; hence, are physically compatible and able to exactly model the kind of continuity which is expected. The method automatically satisfies the well known LBB (inf-sup) stability condition and avoids locking that usually occurs in the numerical computations in the incompressible limit and very low hydraulic conductivity. In contrast to the majority of approaches, our three-field formulation can fully capture dynamic behavior of porous media even in high frequency loading phenomena with considerable fluid acceleration such as liquefaction and biomechanics of porous tissues under rapid external loading. Moreover, we address the importance of consistent initial conditions for poroelasticity equations with the incompressibility constraint, which represent a system of differential algebraic equations. The energy balance equation is derived for the full porous medium and used to assess the stability and accuracy of our time integration. To highlight the capabilities of our method, a variety of numerical studies are provided including verification with analytical and boundary element solutions, wave propagation analyses, hydraulic conductivity effects on damping and frequency content, energy balance analyses, mass lumping considerations, effects of mesh pattern and size, and stability analyses. We also explain some discrepancies commonly found in dynamic poroelasticity results in the literature.

研究动机与目标

• 开发一种针对多孔介质中具有固相与液相耦合响应的动态双相弹性问题的稳定混合有限元格式。
• 通过在单元选择中利用拓扑约束，解决Biot方程混合有限元方法中常见的数值不稳定问题。
• 通过三场（u-w-p）弱形式，确保有效应力与流体通量的精确表示。
• 利用混合有限元空间，实现位移与流体变量的最优收敛率。
• 为模拟饱和多孔介质在动态荷载下的行为提供一个鲁棒的数值框架，尤其适用于岩土力学与土体动力学领域。

提出的方法

• 该方法基于位移（u）、流体通量（w）和孔隙压力（p）的三场弱形式，其推导源自动态双相弹性理论的Biot方程。
• 采用Raviart-Thomas（RT）元对w-p场进行离散，以保证局部质量守恒和正确的通量逼近，这对满足达西定律尤为重要。
• 对位移场u采用标准节点Galerkin有限元，以保证位移场的连续性，并实现固体骨架的最优收敛。
• 应用有效应力概念σᵗ = σ - αpI，其中α = 1（适用于土壤），将总应力与骨架应力及孔隙压力关联起来。
• 边界条件被划分为本质（Dirichlet）与自然（Neumann）两类，分别作用于固体（u, σᵗ·n）与流体（w·n, p）场，以确保弱形式的适定性。
• 所得系统为微分代数方程（DAE），经空间半离散化后，通过时间积分格式求解。

实验结果

研究问题

• RQ1 如何构建一种针对具有固相与液相耦合响应的动态双相弹性问题的稳定混合有限元方法？
• RQ2 在保持物理解释一致性前提下，哪些有限元空间最适合逼近位移、流体通量与孔隙压力？
• RQ3 为何选择Raviart-Thomas元对w与p进行离散，相较于标准Galerkin元，其在精度与稳定性方面有何优势？
• RQ4 在标准基准条件下，该混合格式的收敛行为如何？
• RQ5 有限元空间的拓扑结构如何影响解的稳定性和精度？

主要发现

• 混合有限元格式对位移与流体变量均实现了预期的最优收敛率，符合混合方法的数学理论。
• 采用Raviart-Thomas元对w-p场进行离散，确保了局部质量守恒与精确的通量逼近，这对模拟多孔介质中的流体流动至关重要。
• 采用u、w与p作为独立变量的三场弱形式，为动态双相弹性问题提供了稳定且一致的框架，避免了锁死与振荡现象。
• 有效应力模型σᵗ = σ - αpI成功实现，其中α = 1，确保了机械响应与水力响应之间的正确耦合。
• 通过在边界上对Dirichlet与Neumann子集进行分离处理，边界条件被正确施加，确保了弱形式的适定性。
• 该方法在动态问题中表现出鲁棒性，避免了标准Galerkin格式求解Biot方程时常见的数值不稳定性。