Monolithic Multi-level Overlapping Schwarz Solvers for Fluid Problems

기술 요약: 유체 문제를 위한 모놀리식 다중 레벨 오버랩핑 슈바르츠 솔버

문제 제기
본 논문은 비압축성 유체 역학의 수치 시뮬레이션, 특히 대규모 안장점 (saddle point) 시스템을 초래하는 이산화된 스토크스 (Stokes) 방정식의 해법을 다룹니다. 기하학적 복잡성이 증가하고 경계층을 해결하기 위해 메시 해상도가 세밀해짐에 따라, 직접 솔버는 prohibitive 한 계산 비용과 메모리 소비로 인해 실행 불가능해집니다. 저자들은 표준 2 레벨 도메인 분해 방법, 특히 일반화된 Dryja–Smith–Widlund (GDSW) 접근법의 확장성 한계에 초점을 맞춥니다. 2 레벨 방식에서 부분 도메인의 수가 증가함에 따라 조상 문제 (coarse problem, $A_0$) 의 크기가 커져 결국 그 직접 인수분해가 계산적으로 처리 불가능하게 됩니다.

방법론
저자들은 2 레벨 GDSW 접근법의 확장성 한계를 극복하기 위해 **모놀리식 다중 레벨 오버랩핑 슈바르츠 전구조건부 (preconditioner)**를 제안하고 구현합니다.

• 수학적 공식화: 이 방법은 2 레벨 GDSW 전구조건부에 재귀적 적용을 통해 확장됩니다. 조상 문제 $A_0$를 직접 역행렬화하는 대신, 저자들은 $A_0^{-1}$을 더 조밀한 레벨 (3 레벨 또는 다중 레벨 계층 구조) 에서 유도된 근사치로 대체합니다. 전구조건부는 다음과 같이 정의됩니다:
  $$B_{GDSW, 3-level} = \Phi^{(2)} \left( \Phi^{(3)} A_0^{(3)-1} \Phi^{(3)T} + \sum_{i=1}^{N^{(2)}} R_i^{(2)T} A_i^{(2)-1} R_i^{(2)} \right) \Phi^{(2)T} + \sum_{i=1}^{N^{(1)}} R_i^{(1)T} A_i^{(1)-1} R_i^{(1)}$$
  이 구조를 통해 조상 문제를 직접 솔버가 아닌 도메인 분해의 추가 레벨을 통해 해결할 수 있습니다.

• 구현: 솔버는 Trilinos 생태계의 ShyLU 패키지의 일부인 FROSch(Fast and Robust Overlapping Schwarz) 라이브러리 내에서 구현되었습니다.

  • 연동: 확장 가능한 인터페이스가 FROSch와 FEATFLOW 유한 요소 라이브러리를 연결합니다. FEATFLOW 는 속도에 대해 $Q_2$ 요소와 압력에 대해 불연속 $P_1$ 요소를 사용하여 시스템 행렬과 우변 벡터의 조립을 처리합니다.
  • 대수적 성질: 구현은 순수하게 대수적이며, 자유도 (dofs) 의 병렬 분포에 대한 정보만 필요합니다. 이는 로컬 부분 도메인 행렬, 제한/보간 연산자, 그리고 조상 문제를 자동으로 구성합니다.
  • 조상 공간 구성: 부분 도메인 간에 불연속인 압력 변수의 경우, 저자들은 압력 인터페이스 $\Gamma_p$의 필요성을 피하기 위해 전체 부분 도메인에 걸친 상수 함수를 조상 기저 함수로 활용합니다. 속도의 경우, 영공간 (null space) 제한의 에너지 최소화 확장을 사용합니다.

주요 기여

1. 다중 레벨로의 확장: 이 작업은 안장점 유체 문제를 위해 이전에 [7] 에서 소개된 모놀리식 2 레벨 GDSW 전구조건부를 3 레벨 (그리고 잠재적으로 다중 레벨) 프레임워크로 성공적으로 확장했습니다.
2. 확장 가능한 연동: 기하학적 조립과 대수적 도메인 분해 솔버의 효율적인 연동을 가능하게 하는 FEATFLOW 라이브러리와 Trilinos/FROSch 프레임워크 간의 인터페이스 개발.
3. 대규모 병렬 결과: 저자들은 고성능 컴퓨팅 (HPC) 아키텍처에서 약한 확장성 (weak scaling) 결과를 입증하여, 단위 큐브에서의 포아죄유 (Poiseuille) 유동의 경우 32,768 MPI 랭크까지, 복잡한 압출 다이 (extrusion die) 기하학의 경우 4,096 MPI 랭크까지 확장성을 달성했습니다.

수치 결과
저자들은 단위 큐브에서의 포아죄유 유동과 복잡한 압출 다이 기하학에서의 포아죄유 유동이라는 두 가지 테스트 사례를 사용하여 솔버를 평가했습니다.

• 단위 큐브 (구조화된 분해):

  • 2 레벨 솔버: 32,768 개의 부분 도메인에서 2 레벨 솔버는 조상 문제 (407,837 개의 자유도) 를 인수분해할 수 없었기 때문에 실패했습니다. 조상 문제 인수분해를 위한 설정 시간이 병목 현상이 되었습니다.
  • 3 레벨 솔버: 32,768 개의 부분 도메인에서 문제를 성공적으로 해결했습니다. 부분 도메인이 64 배 증가함에 따라 크릴로프 (Krylov) 반복 횟수가 41 회에서 66 회로 증가했지만, 솔버는 견고성을 유지했습니다. 3 레벨 접근법은 대규모 부분 도메인 수에 대해 2 레벨 접근법에 비해 설정 시간을 크게 단축하여 거대한 조상 문제의 실행 불가능한 인수분해를 피했습니다.

• 압출 다이 (비구조화된 분해):

  • 2 레벨 솔버: 4,096 개의 부분 도메인에서 설정 시간은 조상 문제 인수분해 (84.3 초) 에 의해 지배되었습니다.
  • 3 레벨 솔버: 성능 향상을 보여주었습니다. 수치 설정 시간은 2 레벨 솔버에 비해 약 6 배 감소했고, 해결 시간은 거의 2 배 감소했습니다. 3 레벨 접근법은 2 레벨 방식이 확장성 한계에 직면했던 비구조화된 분해를 성공적으로 처리했습니다.

의의 및 주장
본 논문은 이 작업의 주요 의의가 유체 문제에 내재된 2 레벨 GDSW 방법의 확장성 병목 현상을 극복할 수 있는 다중 레벨 모놀리식 오버랩핑 슈바르츠 전구조건부의 성공적인 시연에 있다고 주장합니다. 구체적으로:

• 3 레벨 접근법은 2 레벨 조상 문제가 인수분해가 실행 불가능할 정도로 커지는 경우에도 (최대 32,768 개까지) 훨씬 더 많은 수의 부분 도메인을 가진 문제의 해결을 가능하게 합니다.
• 이 방법은 특히 조상 공간이 커질 때 대규모 문제의 설정 및 해결 시간 모두에서 성능 향상을 제공합니다.
• 이 작업은 이질적인 CPU-GPU 시스템을 갖춘 미래 엑사스케일 아키텍처에서의 효율적인 전산 유체 역학 (CFD) 시뮬레이션을 향한 한 걸음으로, 독일의 "StroemungsRaum" 프로젝트의 일부입니다.

저자들은 매우 높은 부분 도메인 수에서 크릴로프 반복 횟수의 증가에 관해 겸손하게 언급하며, 전구조건부가 아직 점근적 범위에 도달하지 않았을 수 있음을 지적하지만, 반복 횟수의 성장이 막대한 병렬성 증가에 비추어 수용 가능하다고 강조합니다.