[논문 리뷰] A topology-motivated mixed finite element method for dynamic response of porous media

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

Numerical methods in engineering참고 문헌 52인용 수 3

한 줄 요약

이 논문은 포화 다孔성 매체에서의 동적 다孔성 탄성 거동을 모델링하기 위해 위상 구조에 기반한 혼합 유한요소 방법을 제안한다. 유체의 유량과 압력(w-p 필드)에 대해 Raviart-Thomas 요소를 사용하고, 고체의 변위(u 필드)에 대해 표준 Galerkin 유한요소를 사용한다. 이 방법은 세 필드 약한 공식화를 통해 안정성과 정확한 응력-공극압력 결합을 보장하며, 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의 동적 다孔성 탄성 방정식에서 유도된다.
• 유량과 압력 필드(w-p)에 대해 Raviart-Thomas(RT) 요소를 사용하여 국소 질량 보존과 정확한 유량 근사가 보장되며, 특히 Darcy 법칙을 고려할 때 중요하다.
• 고체 뼈대의 연속성과 최적 수렴을 확보하기 위해 변위 필드 u에 대해 표준 노드 기반 Galerkin 유한요소를 사용한다.
• 효과적 응력 개념 σᵗ = σ - αpI 를 적용하며, 토양의 경우 α = 1 으로 설정하여 총 응력과 뼈대 응력, 공극압력 간의 연결을 보장한다.
• 고체(우선 조건: Dirichlet, 자연 조건: Neumann) 및 유체(우선 조건: w·n, p) 필드에 대해 경계 조건을 분리하여 설정함으로써 문제의 잘 정의됨을 보장한다.
• 최종적으로 유도된 시스템은 공간 이산화된 미분-대수방정식(DAE)이며, 시간 적분 기법을 통해 해를 구한다.

실험 결과

연구 질문

• RQ1 고체와 유체의 응답이 결합된 동적 다孔성 탄성에 대해 안정적인 혼합 유한요소 방법을 어떻게 구성할 수 있는가?
• RQ2 변위, 유량, 공극압력을 근사할 때 물리적 일관성을 유지하면서 최적의 유한요소 공간은 무엇인가?
• RQ3 w와 p에 대해 Raviart-Thomas 요소를 사용할 경우, 표준 Galerkin 요소에 비해 정확도와 안정성은 어떻게 향상되는가?
• RQ4 표준 기준 조건 하에서 혼합 형식의 수렴 거동은 어떠한가?
• RQ5 유한요소 공간의 위상은 해의 안정성과 정확성에 어떤 영향을 미치는가?

주요 결과

• 혼합 유한요소 형식은 혼합 방법 이론에 따라 기대되는 바와 같이 변위 및 유체 변수에 대해 최적 수렴률을 달성한다.
• w-p 필드에 Raviart-Thomas 요소를 사용함으로써 국소 질량 보존과 정확한 유량 근사가 보장되며, 다공성 매체 내의 유체 흐름을 모델링하는 데 핵심적이다.
• u, w, p를 독립 변수로 하는 세 필드 형식은 락킹과 진동을 방지하는 안정적이고 일관된 프레임워크를 제공하여 동적 다孔성 탄성에 적합하다.
• 효과적 응력 모델 σᵗ = σ - αpI 는 성공적으로 구현되었으며, α = 1 인 경우 기계적 응답과 수리적 응답 간의 정확한 결합이 보장된다.
• 경계 조건은 도메인 경계상의 서로소인 Dirichlet 및 Neumann 부분집합을 통해 적절히 강제되며, 약한 공식화의 잘 정의됨을 보장한다.
• 이 방법은 수치적 불안정성이 Biot 방정식의 표준 Galerkin 형식에서 흔한 문제를 피하기 때문에 동적 문제에 대해 강력하다.