[논문 리뷰] Multiphysics Finite Element Methods for a Poroelasticity Model

Xiaobing Feng, Zhihao Ge외 1명|arXiv (Cornell University)|2014. 11. 27.

Advanced Numerical Methods in Computational Mathematics참고 문헌 22인용 수 6

한 줄 요약

이 논문은 두 개의 의사압력을 도입하여, 계면체성 모델의 비정적 다상태 유한요소 방법을 제안한다. 이 방법은 시스템을 일반화된 스토크스 문제와 확산 방정식으로 분리함으로써, 이산 에너지 안정성, 최적 수렴성, 그리고 급격한 초기 압력 변화나 제약된 특정 저장계수의 영향을 받지 않는 수치적 시뮬레이션에서의 락킹 현상의 제거를 보장한다.

ABSTRACT

This paper concerns with finite element approximations of a quasi-static poroelasticity model in displacement-pressure formulation which describes the dynamics of poro-elastic materials under an applied mechanical force on the boundary. To better describe the multiphysics process of deformation and diffusion for poro-elastic materials, we first present a reformulation of the original model by introducing two pseudo-pressures, one of them is shown to satisfy a diffusion equation, we then propose a time-stepping algorithm which decouples (or couples) the reformulated PDE problem at each time step into two sub-problems, one of which is a generalized Stokes problem for the displacement vector field (of the solid network of the poro-elastic material) along with one pseudo-pressure field and the other is a diffusion problem for the other pseudo-pressure field (of the solvent of the material). In the paper, the Taylor-Hood mixed finite element method combined with the $P_1$-conforming finite element method is used as an example to demonstrate the viability of the proposed multiphysics approach. It is proved that the solutions of the fully discrete finite element methods fulfill a discrete energy law which mimics the differential energy law satisfied by the PDE solution and converges optimally in the energy norm. Moreover, it is showed that the proposed formulation also has a built-in mechanism to overcome so-called "locking phenomenon" associated with the numerical approximations of the poroelasticity model. Numerical experiments are presented to show the performance of the proposed approach and methods and to demonstrate the absence of "locking phenomenon" in our numerical experiments.

연구 동기 및 목표

• 급격한 초기 압력 변화 상황에서 발생하는 유한요소 근사에서 알려진 수치적 불안정성인 '락킹 현상'을 해결하기 위해.
• 시간 스텝 알고리즘을 개발하여 다상태성 모델을 두 개의 해결 가능한 부분 문제로 분리하기 위해: 이는 이동도와 하나의 의사압력에 대한 일반화된 스토크스 문제, 나머지 하나의 의사압력에 대한 확산 문제이다.
• 보존량에서 유도된 적절한 경계 조건을 통해 일반화된 스토크스 문제의 유일한 해를 보장하기 위해.
• 완전 이산 유한요소 스킴이 이산 에너지 법칙을 만족하고 에너지 노름에서 최적 수렴성을 보임을 입증하기 위해.
• 제약된 특정 저장계수의 극한 과정에서의 거동을 보여주며, 비트의 응고 모델로 복원됨을 보여주기 위해.

제안 방법

• 원래의 다상태성 모델을 두 개의 의사압력 도입을 통해 재구성하여, 한 쪽은 확산 방정식을 만족함으로써 시스템의 분리를 가능하게 한다.
• 시간 이산화된 시스템을 두 부분 문제로 분리한다: 이동도와 하나의 의사압력에 대한 일반화된 스토크스 문제, 나머지 의사압력에 대한 확산 문제.
• PDE 해의 보존량을 이용하여 일반화된 스토크스 문제의 유일한 해를 보장하는 경계 조건을 구성한다.
• 일반화된 스토크스 문제의 해에서 유도된 경계 조건을 확산 문제에 적용하여 일관성과 유일성을 확보한다.
• 속도와 압력에 대해 타일러-후드 혼합 유한요소 방법을 사용하고, 의사압력에는 P1-연속 유한요소 방법을 사용하여 각 부분 문제에 대한 기존의 해법기를 활용할 수 있도록 한다.
• 완전 이산 스킴이 연속 에너지 법칙을 모방하는 이산 에너지 법칙을 만족하고 에너지 노름에서 최적 수렴성을 달성함을 증명한다.

실험 결과

연구 질문

• RQ1 다상태 유한요소 접근법을 설계하여 다상태성 시스템을 두 개의 해결 가능한 부분 문제로 분리하면서도 안정성과 정확성을 유지할 수 있는가?
• RQ2 시간 스텝 알고리즘에서 일반화된 스토크스 문제와 확산 문제의 유일한 해를 보장하는 방법은 무엇인가?
• RQ3 제안된 방법은 이전의 혼합 유한요소 방법에서 관찰된 급격한 초기 압력 변화 상황에서의 락킹 현상을 제거하는가?
• RQ4 제약된 특정 저장계수가 0으로 수렴하는 극한에서 이 방법은 어떻게 행동하며, 비트의 응고 모델을 회복하는가?
• RQ5 기존의 스토크스 문제와 확산 문제 해법기의 코드를 수정하지 않고도 이 방법을 구현할 수 있는가?

주요 결과

• 보존량에서 유도된 경계 조건과 부분 문제 간의 결합을 통해 제안된 방법은 두 부분 문제의 유일한 해를 보장한다.
• 완전 이산 유한요소 스킴은 연속적인 에너지 법칙을 반영하는 이산 에너지 법칙을 만족한다.
• строго한 오차 추정을 통해 에너지 노름에서 최적 수렴 속도가 확보됨을 입증하였다.
• 수치 실험을 통해 압력 진동의 부재와 락킹 현상의 제거를 확인하였으며, 급격한 초기 압력 변화가 있는 매우 작은 시간 단계에서도 동일한 결과를 얻었다.
• 제약된 특정 저장계수가 0으로 수렴하는 극한에서도 방법이 안정적으로 유지되며, 정확히 비트의 응고 모델로 복원됨을 보였다.
• 압력의 정규성에 민감하지 않다. 압력은 더 이상 주요 변수가 아니라 의사압력에서 유도된 양이기 때문이다.