[논문 리뷰] Efficient $p$-multigrid method based on an exponential time discretization for compressible steady flows

Shujie Li|arXiv (Cornell University)|2018. 06. 30.

Advanced Numerical Methods in Computational Mathematics참고 문헌 14인용 수 3

한 줄 요약

이 논문은 고차분리간점유분해법을 사용한 압축성 정 steady 유동을 해결하기 위한 효율적인 p-multigrid 프레임워크를 제안한다. 지수 시간 적분과 s단계의 전처리된 룬게쿠타 스무딩기를 결합하여, 완전히 암시적 방법에 비해 메모리 사용을 크게 줄이고 p-독립적 수렴성을 달성한다. 이는 3차원 고차 시뮬레이션에 매우 적합하다.

ABSTRACT

An efficient multigrid framework is developed for the time marching of steady-state compressible flows with a spatially high-order ($p$-order polynomial) modal discontinuous Galerkin method. The core algorithm that based on a global coupling, exponential time integration scheme provides strong damping effects to accelerate the convergence towards the steady state, while high-frequency, high-order spatial error modes are smoothed out with a $s$-stage preconditioned Runge-Kutta method. Numerical studies show that the exponential time integration substantially improves the damping and propagative efficiency of Runge-Kutta time-stepping for use with the $p$-multigrid method, yielding rapid and $p$-independent convergences to steady flows in both two and three dimensions.

연구 동기 및 목표

• 고차분리간점유분해법을 사용한 정 steady 압축성 유동 시뮬레이션에 대해 강력하고 효율적인 p-multigrid 프레임워크를 개발하는 것.
• 고차 공간 이산화에 대해 수렴 효율성과 메모리 소비 측면에서 전통적 시간적분 방법의 한계를 극복하는 것.
• 지수 시간 적분과 전처리된 룬게쿠타 스무딩기를 통합하여 고속의 p-독립적 수렴을 달성하는 것.
• 특히 고차 다항식 순서를 가진 3차원 문제에서 완전히 암시적 방법에 비해 메모리 사용량을 줄이는 것.
• 다양한 p-순서에서 2차원 및 3차원 기준 문제에 대해 이 방법의 효율성과 확장성을 입증하는 것.

제안 방법

• 강한 수치적 감쇠와 정 steady 상태로의 빠른 수렴을 위해 1차 지수 시간 적분(EXP1)을 사용한다.
• 고주파수, 고차원 공간 오차 모드를 효과적으로 스무딩하기 위해 s단계의 전처리된 룬게쿠타(PRK) 방법을 사용한다.
• 지수 시간 적분과 PRK를 V-사이클 p-multigrid 프레임워크에 통합하여 다항식 순서 전반에서 수렴 속도를 가속화한다.
• 지수 시간 적분 기법에서 행렬 지수를 효율적으로 계산하기 위해 Krylov 부분공간 방법을 적용한다.
• 강성 있는 자코비안 부분은 해석적으로 다루고 비선형 잔여항은 수치적으로 근사하는 전역 결합 지수 방법을 사용한다.
• 안정성과 효율성을 유지하기 위해 국소 고유값 추정 기반의 시간단계 전략을 구현한다.

실험 결과

연구 질문

• RQ1 지수 시간 적분이 고차분리간점유분해법 이산화에 대한 p-multigrid 방법의 수렴 속도를 크게 향상시킬 수 있는가?
• RQ2 제안된 eMG 프레임워크가 2차원 및 3차원 압축성 유동 문제 모두에서 p-독립적 수렴을 달성하는가?
• RQ3 고차 다항식 순서에서 eMG 방법의 메모리 소비는 완전히 암시적 방법과 비교해 어떻게 되는가?
• RQ4 CPU 시간과 반복 횟수 측면에서 eMG 방법은 ILU-GMRES 기반 암시적 솔버에 비해 얼마나 효율적인가?
• RQ5 eMG 프레임워크는 3차원 고차 시뮬레이션에서 메모리 사용량을 줄이면서도 높은 정확도를 유지할 수 있는가?

주요 결과

• eMG 방법은 2차원 및 3차원 비점성 유동에서 다항식 순서 p = 1에서 p = 3까지 p-독립적 수렴 속도를 달성한다.
• 10^5개 요소와 p = 3을 가진 3차원 문제에서, eMG 방법은 메모리 사용량이 0.03 GB에 그치지만, 완전히 암시적 방법은 45 GB를 차지한다.
• 반복 횟수가 더 많지만, p = 3 경우에서 eMG 방법은 ILU-GMRES 암시적 솔버보다 더 빠른 CPU 시간을 기록한다.
• 지수 시간 적분은 강력한 감쇠 효과를 제공하여 행렬 분해 없이도 정 steady 상태로의 수렴을 가속화한다.
• 이 방법은 3차원에서도 높은 효율성과 낮은 메모리 사용량을 유지하여 대규모 고차 시뮬레이션에 실용적이다.
• 비슷하거나 더 좋은 수렴 성능를 달성하면서도 계산 비용과 메모리 포트폴리오 측면에서 완전히 암시적 방법보다 뛰어난 성능을 보인다.