[논문 리뷰] Preconditioning of weighted H(div)-norm and applications to numerical simulation of highly heterogeneous media

Johannes Kraus, Raytcho Lazarov외 3명|arXiv (Cornell University)|2014. 06. 17.

Advanced Numerical Methods in Computational Mathematics참고 문헌 21인용 수 4

한 줄 요약

이 논문은 고도로 이방성 매질에서 2계 타원 문제의 혼합 유한요소 해석에 대해, 추가적인 슈어 여집합 근사와 함께 보조 공간 multigrid(ASMG) 방법을 사용하는 강건한 대수적 다수준 조절자(algebraic multilevel preconditioner)를 제안한다. 이 방법은 메쉬 크기와 물성 대비에 관계없이 최적 수렴성을 확보하며, 최대 787,456 자유도를 가진 SPE10 기준 문제에서 수치적으로 검증되었다.

ABSTRACT

In this paper we propose and analyze a preconditioner for a system arising from a finite element approximation of second order elliptic problems describing processes in highly het- erogeneous media. Our approach uses the technique of multilevel methods and the recently proposed preconditioner based on additive Schur complement approximation by J. Kraus (see [8]). The main results are the design and a theoretical and numerical justification of an iterative method for such problems that is robust with respect to the contrast of the media, defined as the ratio between the maximum and minimum values of the coefficient (related to the permeability/conductivity).

연구 동기 및 목표

• 높은 이방성 다孔성 매질에서 타원 문제의 혼합 유한요소 이산화로 인해 발생하는 대규모 비정칙 손대점 시스템을 해결하는 데 도전한다.
• 물성 계수의 극도로 큰 변동(고대비)이 존재하더라도 효과적이고 효율적인 조절자를 개발한다.
• 메쉬 세분화(h)와 계수 대비(κ)에 대해 강건성을 확보한다. 이는 석유공학 및 환경 모델링에서의 현실적인 시뮬레이션에 매우 중요하다.
• 채널과 내포물 등을 포함한 다스케일 특징을 가진 실제 복잡한 기준 문제, 예를 들어 SPE10 문제에 대해 방법을 검증한다.
• 다양한 메쉬 수준에서의 확장성과 균일한 수렴성을 보장하며, 총 계산 작업이 자유도에 대해 선형적으로 증가함을 확인한다.

제안 방법

• 속도에 대해 Raviart-Thomas RT0 요소, 압력에 대해 조각별 상수 요소를 사용하여 혼합 유한요소로 문제를 수립한다.
• 블록-구조 선형 시스템을 손대점 형태로 구성: [M_α -B_div^T; -B_div 0][u; p] = [0; f], 여기서 M_α는 가중 H(div) 내적을 나타낸다.
• 추가 슈어 여집합 근사 기반의 대수적 다수준 반복(AMLI)-사이클 보조 공간 multigrid(ASMG) 조절자를 적용한다.
• ASMG 솔버에서 W-사이클 스무스러이저를 사용하며, 전처리 및 후처리 스무스링 단계를 각각 한 번 적용하고, 전역 시스템에 대해 MinRes를 Krylov 솔버로 사용한다.
• 계수 k(x)를 통해 가중 H(div) 노름을 도입하고, 이 가중치의 변동에 대해 조절자가 강건하도록 설계한다.
• 이산 Poincaré 부등식과 다수준 전이 연산자를 활용하여 다양한 척도에서 균일한 수렴성을 보장한다.

실험 결과

연구 질문

• RQ1 메쉬 세분화와 극도로 높은 물성 대비에 대해 수렴성이 강건하게 유지되는 조절자를 설계할 수 있는가?
• RQ2 제안된 ASMG 조절자는 채널과 내포물 등의 다스케일 특징을 가진 실제 복잡한 기준 문제, 예를 들어 SPE10 문제에서 어떻게 성능을 발휘하는가?
• RQ3 추가 슈어 여집합 근사와 AMLI-사이클 multigrid의 조합이 고대비 H(div) 문제에 대해 균일한 수렴성을 제공하는가?
• RQ4 비선형 스무스링과 V-사이클 대비 W-사이클 구성이 잔차 감소와 반복 횟수에 어떤 영향을 미치는가?
• RQ5 우리의 이론적 가정을 초월하여 우측항이 0이 아닐 경우 조절자가 얼마나 효과적으로 유지되는가?

주요 결과

• 제안된 ASMG 조절자는 MinRes 솔버와 함께 강건한 수렴성을 확보하며, MinRes와 ASMG 반복 횟수의 곱이 메쉬 수준에 관계없이 거의 일정하게 유지된다.
• SPE10 기준 문제(case [c])에서 자유도가 3,136에서 787,456으로 증가함에 따라 MinRes 반복 횟수는 13에서 18로 뿐만 아니라 증가하여 확장성이 확인된다.
• 한 개의 스무스링 단계를 가진 ASMG W-사이클은 모든 테스트 케이스, 특히 고대비 케이스 [a]와 [b]에서도 균일한 수렴성을 달성하였으며, 반복 횟수가 수준 간에 안정되어 있었다.
• 스무스링 단계 없이(m=0)도 비선형 ASMG V-사이클은 SPE10 케이스에서 강건한 수렴성을 보였으며, 이는 방법의 내성 강건성을 확인한다.
• 조절된 시스템의 조건수는 메쉬 크기와 대비에 관계없이 유계임이 수치적 근사값 ‖π_D̃‖²_Ã를 통해 확인되었다.
• 총 계산 작업은 미세 메쉬 행렬-벡터 곱의 수로 측정되었으며, 자유도에 대해 선형적으로 증가함을 확인하여 최적 복잡도를 입증하였다.