[논문 리뷰] Purely algebraic domain decomposition methods for the incompressible Navier-Stokes equations

Pawan Kumar|arXiv (Cornell University)|2011. 04. 17.

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

한 줄 요약

이 논문은 희소 행렬의 그래프 구조를 기반으로 한 대칭 다중수준 기법을 사용하여 슈어 여백을 근사화함으로써, 비압축성 라우아이-스토크스 방정식을 해결하기 위한 순수 대수적 도메인 분할 방법을 제안한다. 이 방법은 고레이놀즈수(최대 3000)에서도 강력한 수렴성을 확보하며, PCD 및 LSC 조절자보다 뛰어난 성능을 보이며, 설정 및 적용 단계 모두에서 높은 병렬성을 가능하게 한다.

ABSTRACT

In the context of non overlapping domain decomposition methods, several algebraic approximations of the Dirichlet-to-Neumann (DtN) map are proposed in [F. X. Roux, et. al. Algebraic approximation of Dirichlet- to-Neumann maps for the equations of linear elasticity, Comput. Methods Appl. Mech. Engrg., 195, 2006, 3742-3759]. For the case of non overlapping domains, approximation to the DtN are analogous to the approximation of the Schur complements in the incomplete multilevel block factorization. In this work, several original and purely algebraic (based on graph of the matrix) domain decomposition techniques are investigated for steady state incompressible Navier-Stokes equation defined on uniform and stretched grid for low viscosity. Moreover, the methods proposed are highly parallel during both setup and application phase. Spectral and numerical analysis of the methods are also presented.

연구 동기 및 목표

• 비압축성 라우아이-스토크스 방정식의 이산화에서 발생하는 크고 희소하며, 부정부호이자 비대칭인 선형 연립방정식의 해결 과제를 다루기.
• 레이놀즈 수가 100을 초과할 경우 성능이 저하되는 기존의 조절자인 PCD 및 LSC의 한계를 극복하기.
• 기하학적 또는 메쉬에 의존하는 가정 없이, 오직 행렬의 그래프 구조에만 기반한 순수 대수적이고 겹치지 않는 도메인 분할 방법을 개발하기.
• 현대의 다핵 및 분산 컴퓨팅 아키텍처에 적합하게 설정 및 적용 단계 모두에서 높은 병렬성을 확보하기.
• 제안된 조절자들의 강건성과 효율성을 검증하기 위해 스펙트럼 분석 및 수치 분석을 제공하기.

제안 방법

• 희소 행렬의 그래프 구조를 이용하여 도메인-에서-노이만(DtN) 지도를 근사화함으로써 대수적 도메인 분할 조절자를 구성하기.
• 그래프 분할을 통해 겹치지 않는 부분영역을 정의하여 기하학적 정보 없이도 하위구조화를 가능하게 하기.
• 대수적 다중수준 방법을 사용하여 슈어 여백을 근사화하며, 행렬을 그래프로 간주하여 군집 수준의 연산자를 정의하기.
• MSCN, LUM, MSCE, OMSCN, OMSCNR, MPCD 및 LSC 변형 등 새로운 조절자를 도입하며, 이들은 모두 슈어 여백의 대수적 근사에서 유도된다.
• 수렴성을 반복 횟수 및 해법 시간으로 모니터링하기 위해, 결과로 생긴 선형 연립방정식을 해결하기 위해 조절된 GMRES를 적용하기.
• 이산 라우아이-스토크스 시스템의 블록 구조(2×2의 안장점 형태)를 활용하여 블록 기반의 조절 전략을 정의하기.

실험 결과

연구 질문

• RQ1 희소 행렬의 그래프 구조에 기반한 순수 대수적 도메인 분할 방법이 고레이놀즈수에서 비압축성 라우아이-스토크스 방정식에 대해 강건한 조절을 제공할 수 있는가?
• RQ2 슈어 여백의 대수적 근사 방법이 전통적인 PCD 및 LSC 조절자와 비교해 수렴성 및 확장성 측면에서 어떻게 성능을 발휘하는가?
• RQ3 제안된 방법이 설정 및 적용 단계 모두에서 얼마나 높은 병렬성을 유지할 수 있는가?
• RQ4 스트레칭된 격자에서도 대수적 접근 방식이 효과를 유지하는가? (이러한 격자는 경계층 시뮬레이션에서 흔하다.)
• RQ5 제안된 방법이 메쉬 독립적인 수렴성을 달성하고, 레이놀즈 수가 증가함에 따라 여전히 강건한가?

주요 결과

• 제안된 대수적 조절자들(예: OMSCN, OMSCNR)은 레이놀즈 수 최대 3000까지 메쉬 독립적인 수렴성을 확보하며, PCD 및 LSC 방법보다 뛰어난 성능을 보인다.
• 스트레칭된 격자에서 Re = 3000일 경우, OMSCNR 방법은 128개 프로세서에서 119회의 GMRES 반복과 28.4초의 해법 시간을 기록했고, LSC는 1127회의 반복과 362.3초의 해법 시간이 소요되었다.
• Leaky lid-driven cavity 문제에서 Re = 1000일 경우, OMSCNR는 128개 프로세서에서 114회의 반복과 51.0초의 해법 시간을 기록했고, MPCD는 133회의 반복과 164초의 해법 시간이 소요되었다.
• MSCN 및 OMSCN 방법은 모든 시험 케이스에서 뛰어난 성능을 보였으며, 반복 횟수와 해법 시간이 모두 PCD 및 LSC보다 일관되게 낮았다.
• 이 방법들은 강력한 병렬 확장성을 보이며, 설정 및 적용 단계 모두가 매우 높은 병렬성을 확보하여 현대 고성능 컴퓨팅 환경에 적합하다.
• 스펙트럼 분석 결과, 조절된 시스템의 고유값들이 1에 가까이 뭉쳐져 있고 0에서 멀리 떨어져 있음을 확인하여, 빠른 킬로프 하위공간 수렴이 보장된다.