An efficient preconditioner for mixed-dimensional contact poromechanics based on the fixed stress splitting scheme

기술 요약: 고정 응력 분리 기법(Fixed Stress Splitting Scheme)에 기반한 혼합 차원 접촉 포로메카닉스를 위한 효율적인 프리컨디셔너

문제 정의
균열 접촉 포로메카닉스(fracture contact poromechanics)의 수치 시뮬레이션은 CO2 격리, 지열 에너지 생산, 지하 가스 저장과 같은 응용 분야에서 매우 중요합니다. 이러한 문제는 유체 흐름, 매트릭스 변형(Biot poromechanics), 그리고 균열 내 마찰 접촉 역학이 강력하게 결합된 물리 현상을 포함합니다. 균열의 복잡한 기하학적 구조(종종 고차원 도메인에 삽입된 저차원 객체로 표현됨)와 접촉 조건(쿨롱 마찰 및 비침투 제약 조건에 의해 제어되는 stick, slide, open 상태)의 비매끄러운 특성으로 인해 정확한 모델링은 매우 까다롭습니다.

완전 암시적 시간 단계(fully implicit time-stepping)와 뉴턴 방법(Newton's method)을 사용하여 해결할 때, 이 문제들은 거대한 완전 결합 선형 시스템을 생성합니다. 솔버의 효율성과 견고함은 프리컨디셔너(preconditioner)에 크게 의존합니다. 기존 방법들은 다음과 같은 구체적인 난관에 봉착해 있습니다:

1. 사들 포인트 구조(Saddle-Point Structure): 접촉 역학의 증강 라그랑주(Augmented Lagrange) 정식화는 특이 행렬(singular submatrix)을 가진 사들 포인트 구조를 도입하며, 이는 표준 프리컨디셔닝 기법의 적용을 어렵게 만듭니다.
2. 결합 복잡성(Coupling Complexity): 포로메카닉스의 결합을 분리하기 위해 잘 알려진 방법인 고정 응력 분리 기법은 균열 매질로 확장하기가 어렵습니다. 표준 사들 포인트 솔버는 연산자가 결합된 포로메카닉 문제를 포함하지 않는 특정 행렬 구조(예: 라플라시안 유사 연산자)를 가정하는 경우가 많습니다.
3. 확장성(Scalability): 프리컨디셔너는 실제 저류층 시뮬레이션에서 볼 수 있는 대규모 균열 네트워크를 포함하여 문제의 크기가 커지더라도 효율성을 유지해야 합니다.

방법론
저자들은 중첩된 슈르 보완(nested Schur complement) 근사와 선형 변환을 결합한 견고하고 확장 가능한 프리컨디셔너를 제안합니다. 이 접근 방식은 고정 응력 기법을 매트릭스와 균열 서브도메인 모두로 확장하면서 접촉 역학의 특이성을 해결합니다.

방법론은 다음 단계로 진행됩니다:

1. 선형 변환 (전처리): 접촉 역학 서브행렬의 특이성(특히 sticking 상태의 균열)을 처리하기 위해, 자코비안(Jacobian)에 우측 선형 변환 행렬($Q_r$)을 적용합니다. 이 변환은 다른 사들 포인트 문제에서 사용되는 기법에서 영감을 얻어, 접촉 역학 서브행렬을 인터페이스 힘 균형 방정식과 결합함으로써 이를 가역적으로 만듭니다.
2. 중첩된 슈르 보완 분해: 프리컨디셔너는 변환된 자코비안의 블록 LDU 분해의 상삼각 인수를 세 단계의 슈르 보완 제거를 통해 근사합니다:
   • 레벨 1 (접촉 역학): 접촉 역학 자유도를 제거합니다. 변환된 접촉 서브행렬은 인터페이스 힘으로부터 접촉 조건을 분리하기 위해 정확하게 역행렬을 구하거나 블록 대각 근사를 사용합니다.
   • 레벨 2 (운동량 균형): 운동량 균형 방정식을 분리하기 위해 고정 응력 기반 근사를 적용합니다. 여기에는 주변 매트릭스($l_{mat}$)와 균열 서브도메인($l_{frac}$)에 대한 특정 계수를 가진 안정화 행렬($D_{55}$)을 구성하는 과정이 포함됩니다. 결합된 운동량 및 인터페이스 힘 균형 방정식은 대수적 다중 격자(Algebraic Multigrid, AMG)를 사용하여 단일론적(monolithically)으로 해결됩니다.
   • 레벨 3 (인터페이스 유동): 인터페이스 유동을 제거하기 위해 인터페이스 유동 서브행렬을 그 대각 성분으로 근사합니다(인터페이스 유동 시스템에 ILU(0)가 사용됨). 결과적인 슈르 보완은 유체 질량 균형을 나타내며, AMG와 함께 단일론적으로 해결됩니다.
3. 이론적 분석: 저자들은 이 맥락에서 고정 응력 분리 기법에 대한 이론적 정당성을 제공합니다. 단순화된 선형화 모델을 분석함으로써, 재료 특성, 마찰 계수 및 접촉 상태의 안정성에 관한 특정 가정 하에서 순차적 반복 스킴이 수축(contraction)임을 증명합니다. 이 분석은 프리컨디셔너에 사용된 안정화 계수를 검증합니다.

주요 기여

• 새로운 프리컨디셔너 설계: 본 논문은 혼합 차원 포로메카닉스에서 접촉 역학을 고정 응력 분리 기법과 성공적으로 통합하는 프리컨디셔너를 소개합니다. 선형 변환을 통해 접촉 서브행렬의 특이성을 극복하고, 중첩된 슈르 보완을 통해 복잡한 결합을 처리합니다.
• 이론적 수렴 증명: 저자들은 균열 접촉 포로메카닉스에 적용되는 고정 응력 분리 기법에 대한 수렴 보증을 도출했습니다. 이 스킴이 물리적으로 타당한 가정 하에서 수축임을 입증하였으며, 수렴 속도를 마찰 계수 및 접촉 탄성 매개변수와 연결했습니다.
• 확장성 및 견고성: 이 방법은 탄성(elasticity) 및 유체 질량 균형과 같은 주요 하위 문제에 AMG를 사용하고, 인터페이스 유동에 ILU를 사용하여 확장 가능하도록 설계되었습니다. 이는 복잡한 균열 네트워크를 가진 대규모 문제에 적합합니다.

결과
프리컨디셔너는 PorePy 프레임워크를 사용하여 일련의 수치 실험을 통해 검증되었습니다:

• 매개변수 민감도: 단일 균열 모델에 대한 테스트 결과, 반복 횟수는 주로 마찰 계수($F$)에 따라 달라지며, 마찰이 높을수록 더 많은 반복이 필요함을 보여주었습니다. 이는 이론적 분석과 일치합니다. 균열 법선 강성($K_n$)은 수렴에 거의 영향을 미치지 않았습니다.
• 격자 세밀화 (2D): 단일 균열이 있는 2D 모델에서, 직접 서브솔버(direct subsolvers)를 사용할 때 GMRES 반복 횟수가 안정적이었습니다. 직접 솔버를 AMG로 교체했을 때, 자유도가 3단계(three orders of magnitude) 증가함에 따라 반복 횟수가 약 2배 정도 약간 증가하였으며, 이는 좋은 확장성을 입증합니다.
• 기하학적 복잡성 (2D): 교차하는 여러 균열이 있는 실험에서도 프리컨디셔너는 효율성을 유지하였으며, 격자가 세밀해짐에 따라 AMG 서브솔버의 GMRES 반복 횟수가 소폭 증가하는 데 그쳤습니다.
• 3D 확장성: 여러 개의 균열과 거의 $10^6$개의 자유도를 가진 3D 시뮬레이션에서, 격자 세밀화에 따라 GMRES 반복 횟수가 거의 일정하게 유지되어 대규모 문제를 위한 방법의 확장성을 확인했습니다. 3D 케이스에서 선형 솔버의 붕괴(breakdown)는 관찰되지 않았습니다.

의의
본 논문은 제안된 프리컨디셔너가 균열 매질 내 마찰 접촉 포로메카닉스라는 도전적인 문제에 대해 견고하고 확장 가능한 솔루션을 제공한다고 주장합니다. 선형 변환과 중첩된 슈르 보완을 통해 접촉 매카닉스의 사들 포인트 구조를 처리하도록 고정 응력 기법을 확장함으로써, 이 방법은 복잡한 지하 프로세스의 효율적인 시뮬레이션을 가능하게 합니다. 저자들은 이 접근 방식이 라그랑주 승수 기반 접근법을 사용하는 모든 접촉(포로)역학 이산화에 적용 가능하며, 유도 지진, 유체 이동 및 균열된 저류층의 변형을 정확하고 계산 가능한 방식으로 모델링할 수 있는 경로를 제공한다는 점을 강조합니다.