A scalable preconditioning framework for stabilized contact mechanics with hydraulically active fractures

기술 요약: 수력학적으로 활성화된 균열을 포함한 안정화된 접촉 역학을 위한 확장 가능한 프리컨디셔닝 프레임워크

문제 정의
본 논문은 다공성 매질 내에 삽입된 균열 네트워크 내에서의 마찰 접촉 역학 및 유체 흐름의 결합된 물리학을 시뮬레이션하는 데 따르는 계산적 과제를 다룬다. 이 결합된 문제는 지열 에너지 생산, 수압 파쇄 및 CO2 격리와 같은 지하 응용 분야에서 매우 중요하다. 물리적 과정은 밀접하게 결합되어 있다: 균열의 개구부(aperture)와 미끄러짐(slippage)은 유체 흐름을 구동하며, 유체 압력의 변화는 응력 상태를 교란하고 접촉 역학에 영향을 미친다.

저자들은 균열이 명시적으로 표현되는 이산 균열 모델(Discrete Fracture Model, DFM)을 채택한다. 이산화에는 변위장(displacement field)을 위한 저차 연속 유한 요소와 유체 압력을 위한 셀 중심 유한 체적법이 사용된다. 접촉 제약 조건은 압력장과 동일한 격자를 공유하는 셀 중심 라그랑주 승수(traction forces를 나타냄)를 사용하여 부과되며, 이는 필드 간의 보간 필요성을 제거한다. 그러나 이 특정 정식화는 균일하게 inf-sup 안정적(L-B-B 불안정)이지 않으며 전역 대수적 안정화(global algebraic stabilization)를 필요로 한다.

결과적으로 발생하는 선형 시스템은 비선형 문제의 뉴턴-크릴로프(Newton-Krylov) 선형화로부터 유도된 크고 희소한 비대칭 $3 \times 3$ 블록 자코비안(Jacobian) 행렬이다. 행렬 구조는 동적이며, 균열 네트워크가 스틱(stick), 슬립(slip), 개방(open) 영역을 거치며 진화함에 따라 그 성질(대칭 세들 포인트에서 비대칭 부정-정부호 행렬로)이 변한다. 표준 직접 솔버(direct solvers)는 대규모 문제에 대해 계산적으로 불가능하므로, 견고하고 확장 가능한 반복 솔버가 필요하다.

방법론
핵 핵심 기여는 물리 기반의 블록 분할과 최신 대수적 다중격자(Algebraic Multigrid, AMG) 기술에 기초한 확장 가능한 프리컨디셔닝 프레임워크를 설계한 것이다. 이 프레임워크는 전체 결합 다물리 시스템을 단일 물리 서브 문제로 제한하고, 이를 근사적으로 해결한 다음, 솔루션을 다시 전체 시스템으로 확장(prolonging)하는 방식으로 작동한다.

두 가지 별개의 축소 시퀀스(접근법)가 제안되고 분석되었다:

1. $t-p-u$ 접근법 (Traction-Pressure-Displacement):

   • 전략: 변위($u$)에 대한 단일 평형 방정식을 줄이기 위해 트랙션($t$)과 압력($p$) 변수를 동시에 제거한다.
   • 슈르 보완(Schur Complement): 이는 $S = A + C_1 \tilde{H}^{-1}C_2 - Q_1 T^{-1}Q_2$ 형태의 슈르 보완을 결과로 낳는다.
   • 근사: 트랙션 블록 $H$는 블록 대각 행렬 $\tilde{H}$로 근사되고, 압력 블록 $T$는 그 대각 성분인 $\tilde{T}$로 근사된다. 결과적인 슈르 보완 $S$의 역행렬은 AMG 연산자(구체적으로 PETSc의 GAMG)를 사용하여 부정확하게(inexactly) 적용된다.
   • 특징: 이 접근법은 반복당 계산 비용이 저렴하지만(한 번의 AMG 호출), $S$의 양의 정부호성(positive definiteness)에 의존한다. 논문은 개방 영역이 존재할 때(특히 균열 상태에 따라) $S$가 부정-정부호(indefinite)가 될 수 있으며, 이는 AMG 솔버의 효과를 저해할 수 있다고 언급한다.

2. $t-u-p$ 접근법 (Traction-Displacement-Pressure):

   • 전략: 먼저 트랙션($t$)을 제거하여 첫 번째 레벨의 슈르 보완 $S_1 = A + C_1 \tilde{H}^{-1}C_2$(가상 강성을 가진 탄성 행렬)를 형성한다. 그 다음, 변위($u$)를 제거하여 두 번째 레벨의 슈르 보완 $S_2 = T - Q_2 S_1^{-1} Q_1$을 형성한다.
   • 근사: $S_1$의 역행렬은 AMG를 사용하여 적용된다. 두 번째 레벨의 슈르 보완 $S_2$는 대수적으로 계산되거나 국부 부피 추정을 통해 계산되는 대각 고정 응력(fixed-stress) 행렬 $S_K$로 근사된다.
   • 특징: 이 접근법은 프리컨디셔너 적용당 두 번의 AMG 호출을 필요로 하므로 비용이 더 높다(반복당 높은 비용). 그러나 첫 번째 레벨의 슈르 보완 $S_1$이 균열 상태와 관계없이 이론적으로 양의 정부호임이 보장되므로 더 높은 견고성을 제공한다. 잠재적인 부정-정부호성은 대칭이며 중첩된 직접 솔버(nested direct solvers)에 의해 효과적으로 처리될 수 있는 두 번째 레벨의 슈르 보완 $S_2$에 국한된다.

주요 기여

• 프레임워크 설계: 저자들은 수력학적으로 활성화된 균열을 포함한 안정화된 접촉 역학에서 발생하는 비대칭, 부정-정부호 및 시변적(time-varying) $3 \times 3$ 블록 자코비안을 처리할 수 있는 통합된 프리컨디셔닝 프레임워크를 제시한다.
• 비교 분석: 본 논문은 $t-p-u$와 $t-u-p$ 축소 시퀀스를 엄밀한 이론적 및 수치적 비교를 제공한다. $t-p-u$가 이론적으로 효율적이지만, 복잡한 실제 시뮬레이션에서 균열 상태를 예측할 수 없는 경우 $t-u-p$ 접근법이 견고성 측면에서 우월함을 입증한다.
• 확장성 입증: 제안된 프리컨디셔너가 문제의 크기가 증가하거나 격자 내 균열 기반 미지수와 벌크 미지수의 비율이 변하더라도 크릴로프 반복 횟수를 거의 일정하게 유지하는 알고리즘적 확장성을 달성함을 보여준다.

수치 결과
프레임워크는 세 가지 실험 세트에서 테스트되었다:

1. 견고성 (테스트 1): 단일 균열 문제를 통해, $t-p-u$ 접근법은 직접 솔버와 함께 사용할 때는 잘 작동하지만, 슈르 보완이 부정-정부호가 될 경우 AMG를 사용하면 수렴에 실패할 수 있음을 밝혔다. $t-u-p$ 접근법은 테스트된 모든 구성에서 수렴을 유지했다.
2. 약한 확장성 (테스트 2): 균일하게 세분화된 격리(최대 약 760,000개의 미지수) 및 다양한 2D-to-3D 변수 비율에 대한 테스트에서, 선형 솔버의 반복 횟수는 격자 세분화나 균열 네트워크의 상대적 크기와 관계없이 안정적으로 유지되었다(평균 17–24회의 GMRES 반복).
3. 실제 응용 (테스트 3): 9개의 균열과 교차하는 기울어진 우물을 포함한 수압 파쇄 시나리오(100만 개 이상의 미지수) 시뮬레이션을 통해 계산 효율성을 확인했다. 복잡한 스틱/슬립/개방 전이로 인해 비선형 솔버가 많은 액티브 세트(active-set) 단계와 뉴턴 반복을 필요로 했음에도 불구하고, 선형 솔버는 안정적이고 효율적이었다(단계당 평균 24.2회의 GMRES 반복).

의의 및 주장
본 논문은 제안된 프리컨디셔닝 프레임워크가 균열 매질 내의 결합된 마찰 접촉 및 유체 흐름 시뮬레이션의 계산적 병목 현상을 성공적으로 해결한다고 주장한다. 물리 기반의 블록 분할과 다중격자 기술을 활용함으로써, 이 방법은 대규모 병렬 시뮬레이션을 위한 필요한 확장성을 달 achieve 한다.

저자들은 $t-p-u$ 접근법이 이론적 효율성과 낮은 적용 비용 측면에서 약간의 이점이 있지만, 복잡한 균열 진화 중에 발생하는 슈르 보완 행렬의 부정-정부호성에 대한 우수한 견고성 때문에 실무 적용에는 $t-u-p$ 접근법이 더 적합하다고 겸허히 결론짓는다. 이 연구는 점차 커지고 복잡해지는 지질 역학 문제를 해결하기 위한 토대를 마련하였으며, 향후 고성능 컴퓨팅 구현 및 다상 유동(multi-phase flow)으로의 확장을 목표로 한다.