[논문 리뷰] Parameter-robust convergence analysis of fixed-stress split iterative method for multiple-permeability poroelasticity systems

Qingguo Hong, Johannes Kraus외 2명|arXiv (Cornell University)|2018. 12. 31.

Advanced Numerical Methods in Computational Mathematics참고 문헌 23인용 수 5

한 줄 요약

이 논문은 다중투수성 다중투수성 포름탄성계(MPET) 시스템을 위한 고정 스트레스 분할 반복 방법을 제안하며, 모든 물리적 매개변수에 영향을 받지 않는 수렴 속도를 갖는 선형 수렴성을 증명한다. 이 방법은 모든 유체 압력의 합을 통한 안정화를 통해 매개변수에 강건한 수렴성을 확보하며, 조건부 최적화된 MinRes 해법과의 수치적 비교를 통해 광범위한 매개변수 변화에서 뛰어난 효율성과 강건성을 입증한다.

ABSTRACT

We consider flux-based multiple-porosity/multiple-permeability poroelasticity systems describing multiple-network flow and deformation in a poro-elastic medium, sometimes also referred to as MPET models. The focus of the paper is on the convergence analysis of the fixed-stress split iteration, a commonly used coupling technique for the flow and mechanics equations in poromechanics. We formulate the fixed-stress split method in the present context and prove its linear convergence. The contraction rate of this fixed-point iteration does not depend on any of the physical parameters appearing in the model. The theory is confirmed by numerical results which further demonstrate the advantage of the fixed-stress split scheme over a fully implicit method relying on norm-equivalent preconditioning.

연구 동기 및 목표

• 넓은 물리적 매개변수 변화 범위에서도 효과적인 다중투수성 다중투수성 포름탄성계(MPET) 시스템을 위한 강건한 반복적 결합 방법을 개발하는 것.
• 모델 매개변수에 의존하지 않는 고정 스트레스 분할 방법의 선형 수렴성을 MPET 맥락에서 증명하는 것.
• 고정 스트레스 방법이 범위 등가 조건부 최적화를 갖는 완전히 암시적 방법보다 우월함을 이론적 및 계산적 증거로 확립하는 것.
• 비토의 모델에서부터 더 복잡한 다중유체 네트워크를 가진 MPET 시스템으로 고정 스트레스 분할 프레임워크를 확장하는 것.
• 다양한 메쉬 해상도와 매개변수 영역에서 이론적 수렴 속도를 수치 실험을 통해 검증하는 것.

제안 방법

• 고정 스트레스 분할 방법은 MPET 시스템에 대한 고정점 반복으로 공식화되며, 유동 및 기계적 방정식을 분리한다.
• 모든 유체 압력의 합을 통한 안정화를 도입하여 강건성 확보 및 매개변수에 강건한 수렴 분석이 가능하도록 한다.
• 수렴 분석은 유량 기반 MPET 시스템에 대한 이전 연구의 안정성 결과에 기반하며, 에너지 노름을 사용하여 수축 추정치를 유도한다.
• 변분 공식을 사용하여 MPET 시스템의 약한 형태를 유도함으로써 유한요소 이산화 및 반복적 해법이 가능해진다.
• 이 방법은 변위 및 압력에 대해 표준 유한요소를 사용하며, 유량은 다르시의 법칙에 따라 계산된다.
• 수치 실험에서는 고정 스트레스 방법을 조건부 최적화된 MinRes 해법과 비교하며, 잔차를 10^8 배 감소시키는 데 소요된 반복 횟수를 측정한다.

실험 결과

연구 질문

• RQ1 고정 스트레스 분할 방법은 다중유체 네트워크를 가진 MPET 시스템으로 확장 가능하며, 물리적 매개변수에 영향을 받지 않고 수렴함을 증명할 수 있는가?
• RQ2 고정 스트레스 반복에서 MPET 시스템의 매개변수에 강건한 수렴을 보장하기 위한 안정화 전략은 무엇인가?
• RQ3 고정 스트레스 방법의 수렴 속도는 범위 등가 조건부 최적화를 갖는 완전히 암시적 방법과 비교해 볼 때 어떻게 되는가?
• RQ4 고정 스트레스 방법은 투수성, 강성, 저장 계수의 극단적인 변화가 있을 때에도 효율성과 강건성을 유지하는가?
• RQ5 이론적 수축 속도 추정치는 다양한 메쉬 해상도와 매개변수 영역에서 수치적으로 검증 가능한가?

주요 결과

• 고정 스트레스 분할 방법은 MPET 모델의 모든 물리적 매개변수, 즉 투수성, 강성, 저장 계수에 영향을 받지 않는 선형 수렴을 보이며, 이는 수축 속도가 매개변수에 독립적임을 의미한다.
• 네 네트워크를 가진 MPET 모델에서 고정 스트레스 방법은 잔차를 10^8 배 감소시키는 데 항상 10~14회의 반복만을 요구하였으며, 투수성, 강성, 저장 계수의 극단적인 값(예: λ·10^8 또는 K₃·10^10)이 적용된 경우에도 동일하게 유지되었다.
• 조건부 최적화된 MinRes 해법은 더 많은 반복 횟수를 필요로 하였으며, 특히 투수성 또는 강성이 넓은 범위로 변화할 경우 최대 45회까지 도달하여 낮은 강건성을 보였다.
• 모든 테스트된 매개변수 범위, 특히 λ(기본값의 10^8 배), K₃(기본값의 10^10 배), K(기본값의 10^2 배)의 극단적 값에서도 방법은 안정적이고 효율적이었다.
• 수치 결과는 이론적 수축 속도 추정치를 확인하였으며, 모든 매개변수 변화와 메쉬 해상도에서 반복 횟수가 낮은 값(10~14)으로 안정화됨을 보였다.
• 모든 테스트 시나리오에서 고정 스트레스 방법은 범위 등가 조건부 최적화를 갖는 MinRes 방법을 능가하였으며, 더 뛰어난 강건성과 효율성을 입증하였다.