A Stabilized Mortar Method for Discontinuities in Geological Media with Non-Conforming Grids

1. 연구 배경 및 문제 제기 (Problem Statement)

• 배경: 지열 에너지 생산, 지하 가스 저장, 탄소 포집 및 저장 (CCS) 등 지하 시스템의 안전성 및 성능 평가에는 지질 매체 내 단층 (fault) 과 균열 (fracture) 의 정확한 수치 시뮬레이션이 필수적입니다.
• 주요 문제:
  1. 기하학적 비정합성 (Non-conformity): 산업 표준 격자 형식인 '코너 포인트 그리드 (corner-point grid)'를 사용할 경우, 단층 양쪽의 격자가 정렬되지 않는 (비정합) 상황이 빈번하게 발생합니다.
  2. 접촉 조건 enforcement: 단층의 마찰 접촉 거동 (미끄러짐, 고정, 개방) 을 정확히 모델링하는 것은 여전히 큰 도전 과제입니다.
  3. 수치적 불안정성: Lagrange 승수법 (LMM) 을 사용하여 접촉 조건을 엄격하게 부과할 때, 단층 양쪽 격자의 세분화 비율 (mesh refinement ratio) 이 다르면 inf-sup 안정성 조건이 위반되어 해가 진동하거나 발산할 수 있습니다. 특히, 비 모르터 (non-mortar) 측의 격자가 더 미세할 때 이러한 진동이 심화됩니다.
  4. 기존 방법의 한계: 페널티 방법은 제약 조건이 근사적이고 조건수가 나빠지며, 기존 안정화된 LMM (예: 듀얼 노드, 버블 함수) 은 복잡한 지질 기하학 (교차점, 경계 조건) 에서 구현이 어렵거나 여전히 진동 문제를 겪을 수 있습니다.

2. 제안된 방법론 (Methodology)

이 논문은 비정합 격자 (Non-conforming grids) 환경에서 단층 접촉 문제를 해결하기 위해 안정화된 모르터 방법 (Stabilized Mortar Method) 을 제안합니다.

• 기본 프레임워크:

  • 이산 불연속 모델 (DFM): 단층을 차원이 낮은 격자로 명시적으로 이산화합니다.
  • 모르터 방법 (Mortar Method): 단층 한쪽 (모르터 측) 과 다른 쪽 (비 모르터 측) 의 격자가 정렬되지 않아도 되도록 약한 형식 (weak formulation) 으로 접촉 조건을 부과합니다.
  • Lagrange 승수: 인터페이스의 힘 (traction) 역할을 하는 Lagrange 승수를 도입하여 접촉 제약을 정확히 처리합니다.

• 핵심 혁신: 안정화 기법 (Stabilization Strategy)

  • 공간 선택: 변위장은 조각별 선형 (piecewise linear, $Q_1$), Lagrange 승수장은 조각별 상수 (piecewise constant, $P_0$) 로 선택합니다. $P_0$ 승수는 구현의 용이성 (경계 조건 적용, 교차점 처리) 이 뛰어나지만, 전통적으로 inf-sup 조건을 만족하지 못합니다.
  • 트랙션 점프 안정화 (Traction-Jump Stabilization): inf-sup 조건을 회복하기 위해 제약 방정식에 트랙션 점프 (traction-jump) 항을 추가합니다. 이는 인접 요소 간의 힘 불연속을 패널티화하여 안정성을 확보합니다.
  • 자동 스케일링 알고리즘:
    • 기존 방법과 달리 사용자가 추가적인 안정화 파라미터를 선택할 필요가 없습니다.
    • 매크로 요소 (Macro-element) 분석을 기반으로 로컬 Schur 여분 (Schur complement) 을 근사하여 안정화 항의 적절한 스케일링 계수 ($S_{EKL}$) 를 자동으로 계산합니다.
    • 이는 진동 모드를 제거하면서도 주 변수의 정확도를 과도하게 떨어뜨리지 않는 균형을 이룹니다.

• 해법 전략:

  • Coulomb 마찰 법칙을 적용하여 접촉 상태를 '고정 (stick)', '미끄러짐 (slip)', '개방 (open)' 영역으로 분류합니다.
  • Active-set 전략을 사용하여 부등식 제약 조건을 반복적으로 풀며, Newton-Raphson 방법을 통해 비선형 시스템을 해결합니다.

3. 주요 기여 (Key Contributions)

1. 새로운 안정화된 모르터 공식화: 조각별 상수 Lagrange 승수 ($P_0$) 와 트랙션 점프 안정화 항을 결합하여, 비정합 격자 환경에서도 inf-sup 안정성을 보장하는 새로운 형식을 제시했습니다.
2. 파라미터 없는 자동 스케일링: 사용자 정의가 필요 없는 알고리즘을 개발하여, 로컬 Schur 여분 분석을 통해 안정화 항의 크기를 자동으로 결정했습니다.
3. 과도한 진동 제거: 비 모르터 측 격자가 더 미세한 경우 ($r_h < 1$) 에 발생하는 전형적인 힘 (traction) 진동을 성공적으로 억제하고 매끄러운 힘 분포를 복원했습니다.
4. 실제 적용성 입증: 복잡한 기하학 (교차점, 경계 조건) 과 산업 표준 격자 (코너 포인트 그리드) 를 사용하는 실제 지질 문제에 적용 가능성을 검증했습니다.

4. 수치 실험 및 결과 (Numerical Results)

논문은 5 가지 테스트 케이스를 통해 방법론을 검증했습니다.

• Test 1 (패치 테스트): 비정합 격자 교차점에서 일정한 힘 전달 능력을 검증. 제안된 방법은 표준 노드 승수법에서 발생하는 진동 없이 정확한 상수 힘 분포를 재현했습니다.
• Test 2 (일정 미끄러짐): 경사 단층에서의 미끄러짐 접촉. 해석해와 완벽하게 일치하며, 단일 액티브 세트 업데이트로 수렴했습니다.
• Test 3 (압축 하중 하의 단일 균열): 균열 끝단에서의 특이점 (singularity) 을 포함한 해석해와 비교. 격자 비정합성과 안정화 기법의 결합이 진동을 완화하여 해석해와 잘 일치함을 보였습니다.
• Test 4 (압축 및 전단 하중 하의 두 블록):
  • 핵심 결과: 비 모르터 측 격자가 더 미세한 경우 ($r_h = 1/2$), 기존 inf-sup 안정 이론을 가진 방법들 (듀얼 승수, 버블 안정화) 은 여전히 힘의 진동을 보였습니다.
  • 반면, 제안된 트랙션 점프 안정화 방법은 매끄러운 힘 분포를 제공하여 기준 해 (fine conforming mesh) 와 거의 동일한 결과를 얻었습니다. 이는 마찰 법칙 하에서 진동이 접촉 상태 판별 (active-set) 의 안정성을 해치는 것을 방지합니다.
• Test 5 (지질 매체 내 대수층 물 추출):
  • 실제 적용: 코너 포인트 그리드를 사용하여 10 년 간의 유체 추출에 따른 단층 활성화 시뮬레이션.
  • 결과: 유체 압력 강하로 인한 지반 수축이 단층의 미끄러짐과 개방을 유발하는 과정을 매끄럽게 시뮬레이션했습니다. 복잡한 지질 구조와 비정합 격자에서도 방법론이 견고하게 작동함을 입증했습니다.

5. 의의 및 결론 (Significance and Conclusion)

• 기술적 의의: 이 연구는 지질 역학 시뮬레이션에서 가장 까다로운 문제 중 하나인 비정합 격자와 마찰 접촉의 결합을 해결했습니다. 특히, 산업 표준 격자 형식을 사용할 때 발생하는 기하학적 불일치로 인한 수치적 불안정성을 안정화 기법 없이 해결할 수 있는 자동화된 알고리즘을 제시했습니다.
• 실용적 가치: 제안된 방법은 복잡한 지하 구조물 (단층, 교차 균열) 을 가진 실제 지하 자원 개발 및 지질 재해 평가 문제에 직접 적용 가능하여, 기존 방법론의 한계를 극복하고 더 정확하고 신뢰할 수 있는 예측을 가능하게 합니다.
• 미래 전망: 더 복잡한 기하학적 시나리오와 고급 접촉 알고리즘을 결합하여 방법론의 견고성을 더욱 향상시킬 계획입니다.

요약하자면, 이 논문은 $P_0$ Lagrange 승수와 자동 스케일링 트랙션 점프 안정화를 결합하여, 비정합 격자 환경에서도 진동 없이 안정적인 접촉 해석을 가능하게 하는 획기적인 수치 기법을 제시했습니다.