[논문 리뷰] A finite element method with strong mass conservation for Biot's linear consolidation model

Béatrice Rivière, Guido Kanschat|arXiv (Cornell University)|2017. 12. 20.

Advanced Numerical Methods in Computational Mathematics참고 문헌 10인용 수 46

한 줄 요약

이 논문은 Biot의 선형 응고 모델을 위한 H(div)-형식 유한요소법을 제안하며, 혼합 유한요소 형식과 내부 페널티 불연속 Galerkin 항을 조합하여 강력한 질량 보존을 보장한다. 이 방법은 유체 속도와 변위에 대해 최적의 L² 수렴률을 달성하며, 유체 및 고체의 압축성에 관계없이 안정적인 성능을 보이며, 메esh 크기나 시간 간격과 무관하게 정확한 질량 균형을 유지한다.

ABSTRACT

An H(div) conforming finite element method for solving the linear Biot equations is analyzed. Formulations for the standard mixed method are combined with formulation of interior penalty discontinuous Galerkin method to obtain a consistent scheme. Optimal convergence rates are obtained.

연구 동기 및 목표

• 다공성 탄성 해석에서 강력한 질량 보존을 보장하는 유한요소법을 개발한다.
• 기존의 Taylor-Hood 혼합 방법에서 발생하는 유체 속도와 변위 간의 수렴률 격차를 제거한다.
• 일반적인 물성 매개변수, 특히 압축성 없는 유체(cₛ = 0) 조건에서도 유체 속도 및 변위 필드에 대해 최적의 L² 수렴률을 달성한다.
• Biot-Willis 상수 α와 저장계수 cₛ에 대해 강건성을 확보하며, 특히 압축성 없는 유체의 극한 조건에서 성능을 유지를 한다.
• 투영 기반 연결을 통한 유체 및 고체 속도 간의 일관된 결합을 제공하여, 다양한 매개변수 영역에서 안정성과 정확성을 확보한다.

제안 방법

• 유체 속도와 고체 변위 모두에 대해 H(div)-형식 Raviart-Thomas 유한요소를 사용하여 국소 질량 보존을 보장한다.
• 유체 압력과 속도에 대한 표준 혼합 형식과 내부 페널티 불연속 Galerkin 항을 조합하여 변위 필드의 H¹ 연속성을 강제한다.
• 속도 및 변위에 동일한 벡터 유한요소 공간을 사용하여 L²에서 동일한 순서의 근사가 가능하도록 한다.
• 최적의 오차 추정을 유도하기 위해 변위 필드의 발산에 대한 초수렴 가정(superapproximation assumption)을 사용한다.
• 시간 이산화에 θ-스키마를 적용하며, θ = 0.501을 사용하여 강력한 A-안정성을 확보하고 시간 이산화 오차를 감소시킨다.
• 메쉬가 거친 경우에도 기계 정밀도 수준(마이크로스코픽 정밀도) 내에서 이산 질량 균형 잔여항 Δm(t)을 통해 질량 보존을 검증한다.

실험 결과

연구 질문

• RQ1 투영 기반 결합에 의존하지 않고 Biot의 응고 모델에서 강력한 질량 보존을 보장하는 유한요소법을 구축할 수 있는가?
• RQ2 H(div)-형식 혼합 방법과 내부 페널티 DG 항을 조합하면 유체 속도 및 변위 필드에 대해 최적의 수렴률을 달성할 수 있는가?
• RQ3 특히 압축성 없는 유체 극한 조건(cₛ = 0)에서 Biot-Willis 상수 α와 저장계수 cₛ에 대해 강건한가?
• RQ4 안정성과 일관성을 유지하면서 L²에서 속도 및 변위에 대해 동일한 순서 근사가 가능하고 최적의 수렴률을 달성할 수 있는가?
• RQ5 메쉬 크기나 시간 간격과 무관하게 이산 수준에서 질량 균형을 정확히 유지하는가?

주요 결과

• RT₁/Q₁ 요소를 사용할 경우, 유체 속도 및 변위에 대해 최적의 L² 수렴률 h²를 달성하며, 수치 실험으로 확인되었다.
• RT₂/Q₂ 요소를 사용할 경우, L² 수렴률이 h³에 도달하여 고차 근사에서 최적 수렴을 보였다.
• 변위 필드의 발산 오차는 순서 h¹로 수렴하며, 이는 이론적 가정과 수치적 검증에 부합한다.
• 이산 질량 균형 잔여항 Δm(t)는 cₛ = 0 및 α ≠ 1를 포함한 모든 테스트 조건에서 기계 정밀도 수준(약 10⁻¹⁷) 내에 유지되었다.
• 이전 방법이 성능이 저하되거나 최적 수렴을 달성하지 못하는 cₛ = 0(압축성 없는 유체) 조건에서도 이 방법은 강건한 성능을 유지한다.
• 수치 실험 결과, 다양한 물성 매개변수 하에서 방법이 안정적이고 정확하며, 메쉬 세분화에 따라 오차가 예상된 비율로 감소하는 것으로 확인되었다.