[논문 리뷰] Monolithic and splitting based solution schemes for fully coupled quasi-static thermo-poroelasticity with nonlinear convective transport

Mats K. Brun, Elyes Ahmed외 3명|arXiv (Cornell University)|2019. 02. 15.

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

한 줄 요약

이 논문은 대류 열수송을 수반하는 비선형이고 완전히 결합된 정적 열다孔탄성 문제를 해결하기 위한 6개의 단일체적이고 분할 기반 반복적 방법을 제안한다. 3필드 모델을 5필드 체계로 확장하고 고정 응력 분할을 적용한 L-스키마 선형화를 통해, 다양한 결합 강도와 메esh 정밀도에서 안정적인 수렴을 달성하였으며, 안정화 기법이 비안정화 대비 수렴성과 안정성을 크게 향상시켰다.

ABSTRACT

This paper concerns splitting-based iterative procedures for the coupled nonlinear thermo-poroelasticity model problem. The thermo-poroelastic model problem we consider is formulated as a three-field system of PDE's, consisting of an energy balance equation, a mass balance equation and a momentum balance equation, where the primary variables are temperature, fluid pressure, and elastic displacement. Due to the presence of a nonlinear convective transport term in the energy balance equation, it is convenient to have access to both the pressure and temperature gradients. Hence, we introduce these as two additional variables and extend the original three-field model to a five-field model. For the numerical solution of this five-field formulation, we compare three approaches that differ by how we treat the coupling/decoupling between the flow and/from heat and/from mechanics; these approaches have in common a simultaneous application of the fixed-stress splitting scheme on both the non-linearity and the coupling structure of the problem. More precisely, the derived procedures transform a nonlinear and fully coupled problem into a set of simpler subproblems to be solved sequentially in an iterative fashion. We provide a convergence proof for the derived algorithms, and validate our results through several numerical examples.

연구 동기 및 목표

• 대류 열수송을 포함하는 열, 유량, 기계적 현상이 완전히 결합된 비선형 열다孔탄성 문제를 해결하는 데 도전한다.
• 강한 결합이나 비선형성으로 인해 기존 방법이 실패하는 문제를 해결하기 위해 안정화된 반복적 분할 프레임워크를 도입한다.
• 열, 유량, 기계적 현상의 세 가지 물리적 과정에 대한 모든 가능한 결합/분리 조합을 다루는 6개의 독립된 반복 알고리즘을 개발한다.
• 강한 결합과 비선형 대류 항이 존재하는 상황에서도 L-스키마 선형화와 고정 응력 분할을 통해 수렴성과 안정성을 확보한다.
• 다양한 수치 벤치마크를 통해 제안된 방법의 안정성과 효율성을 검증하며, 다양한 매개변수 영역과 메쉬 정밀도에서의 성능을 입증한다.

제안 방법

• 온도, 압력, 변형률의 원래 3필드 열다孔탄성 모델을 압력 및 온도 기울기 변수를 추가함으로써 5필드 체계로 확장한다.
• 비선형 대류 항과 결합된 구조 양쪽에 동시에 고정 응력 분할 기법을 적용하여 하위 문제를 순차적으로 해결할 수 있도록 한다.
• 에너지 방정식 내 비선형 대류 항의 선형화 및 안정화를 위해 L-스키마를 사용하여 수렴성과 안정성을 보장한다.
• 모든 6개 알고리즘을 열, 유량, 기계적 현상의 서로 다른 결합 순서로 하위 문제를 번갈아가며 해결하는 반복 절차로 기술한다.
• 이론적 분석에서 유도된 시간단계 제약 조건 하에 수렴 증명을 유도하였으며, 실질적인 성능은 이론적 한계를 초월해 안정적으로 유지됨을 관찰하였다.
• 공간 이산화를 위해 혼합 유한요소 방법을 구현하여 수치 해에서 국소 질량 및 운동량 보존을 보장한다.

실험 결과

연구 질문

• RQ1 대류 열수송에 의해 유도되는 강한 결합과 비선형성 문제를 효과적으로 다룰 수 있는 반복 분할 기법은 어떻게 설계할 수 있는가?
• RQ2 L-스키마를 통한 안정화 기법이 다양한 결합 강도와 매개변수 영역에서 분할 알고리즘의 수렴성과 안정성에 미치는 영향은 어떠한가?
• RQ3 열, 유량, 기계적 현상의 모든 6가지 가능한 결합/분리 구성이 통합된 반복 프레임워크를 통해 효과적으로 구현되고 검증될 수 있는가?
• RQ4 안정화된 알고리즘과 비안정화된 알고리즘 간의 반복 횟수 및 매개변수 변화에 대한 민감도는 어떻게 비교될 수 있는가?
• RQ5 이론적 수렴 한계로 도출된 시간단계 크기 제약 조건이 복잡한 비선형 다물리적 시뮬레이션에서 실질적인 성능과 얼마나 일치하는가?

주요 결과

• 제안된 6개 알고리즘은 메쉬 정밀도(최대 40×40 요소)와 다양한 매개변수 영역에서 강한 결합 상태를 포함해 안정적인 수렴성과 안정성을 모두 확보하였다.
• 다양한 메쉬 크기와 결합 강도에서 반복 횟수가 일관되게 낮게 유지됨(12–18회)하여 뛰어난 확장성과 안정성을 나타냈다.
• L-스키마를 통한 안정화 기법은 매개변수 변화에 대한 민감도를 감소시켰다: 비안정화 알고리즘은 강한 결합 영역에서 실패하거나 과도한 반복 횟수를 요구하는 반면, 안정화된 버전은 여전히 효과적으로 작동하였다.
• 일정한 열원(z = 2×10⁻⁴ W m⁻³ K⁻¹) 조건에서 온도 필드가 압력과 변형률에 상당한 영향을 미치며, 이는 시스템 내 강한 결합 효과를 확인한다.
• 이론적 시간단계 제약 조건은 지나치게 보수적인 것으로 밝혀졌으며, 실질적인 성능은 유량이 유한한 범위 내에 있을지라도 합리적인 시간단계로 안정적으로 유지되었다.
• 기울기 변수를 추가한 5필드 체계는 비선형 대류 항을 정확하게 처리할 수 있게 하여 수치 안정성과 해의 정확도를 향상시켰다.