[논문 리뷰] Palindromic discontinuous Galerkin method for kinetic equations with stiff relaxation

David Coulette, Emmanuel Franck외 3명|arXiv (Cornell University)|2016. 12. 30.

Lattice Boltzmann Simulation Studies참고 문헌 5인용 수 3

한 줄 요약

이 논문은 경직성 있는 붕괴를 가진 운동 방정식을 위한 고차수, 점근 보존성 피아노형 비연속 갈레르킨 라티스 볼츠만 방법을 소개한다. 음성적 상향 흐름 DG 공간 이산화와 피아노형 시간 통합을 조합하여 고차수 정확도(최대 6차), 높은 CFL 수치에서의 안정성, 그리고 작은 붕괴 시간에 대한 강건성을 달성한다. 이는 도전적인 시험 케이스, 특히 $\tau > 0$인 리만 문제에서 최적의 수렴 속도와 안정성을 입증한다. 방법은 점근 보존성을 유지하며 삼각형 선형 시스템을 통해 효율적이고 저메모리 계산을 가능하게 한다.

ABSTRACT

We present a high order scheme for approximating kinetic equations with stiff relaxation. The objective is to provide efficient methods for solving the underlying system of conservation laws. The construction is based on several ingredients: (i) a high order implicit upwind Discontinuous Galerkin approximation of the kinetic equations with easy-to-solve triangular linear systems; (ii) a second order asymptotic-preserving time integration based on symmetry arguments; (iii) a palindromic composition of the second order method for achieving higher orders in time. The method is then tested at orders 2, 4 and 6. It is asymptotic-preserving with respect to the stiff relaxation and accepts high CFL numbers.

연구 동기 및 목표

• 경직성 있는 붕괴를 가진 운동 방정식을 위한 고차수 수치적 방법을 개발하여, 높은 CFL 수치에서도 안정성과 정확도를 유지한다.
• 명시적 방법과 표준 LBM의 한계를 극복하기 위해, 삼각형 선형 시스템을 가진 음성적 상향 흐름 비연속 갈레르킨(DG) 공간 이산화를 도입한다.
• 대칭적인 두 번째 차수 시간 통합기법을 피아노형 조합하여 고차수 시간 정확도를 확보함으로써, $\tau \to 0$일 때 점근 보존성을 보장한다.
• 경직성 붕괴 영역, 예를 들어 점성 및 접촉 불연속 문제를 포함하여, 방법의 강건성과 수렴 특성을 시험한다.
• 병렬 및 하이브리드 아키텍처에 적합한 효율적이고 저메모리 계산을 가능하게 하며, 향후 고차원으로의 확장 가능성을 고려한다.

제안 방법

• 방법은 음성적 상향 흐름 DG 공간 이산화를 사용하며, 유도된 선형 시스템은 삼각형이 되어 명시적으로 해석 가능하여 높은 CFL 수치와 안정성을 보장한다.
• 대칭성 논리를 기반으로 두 번째 차수, 대칭적인 시간 통합기법을 구성하여 $\tau \to 0$일 때 점근적 행동을 유지한다.
• 두 번째 차수 방법의 피아노형 조합을 통해 시간에 대한 고차수 정확도를 확보하며, 수식은 수지와 카한-리 방법을 사용해 4차 및 6차로 시험하였다.
• 운동 모델은 그라일의 벡터형 LBM에 기반하며, $2m$개의 이산 속도 $v_i = \pm\lambda$와 보존 변수 $w$ 및 유량 $q(w)$를 복원하는 균형 분포를 정의한다.
• 해법은 점근 보존성이다: $\tau \to 0$일 때 해는 원래의 쌍곡 보존 법칙 시스템으로 수렴한다.
• 1D 리만 문제 및 점성 시험 케이스에서 수치 시험을 수행하였으며, 2차, 4차, 6차에서의 수렴 분석과 $\tau > 0$에 대한 안정성 분석을 실시하였다.

실험 결과

연구 질문

• RQ1 음성적 DG 방법과 상향 흐름을 사용한 해법이 경직성 붕괴를 가진 운동 방정식에서 고차수 정확도와 높은 CFL 수치에서의 안정성을 달성할 수 있는가?
• RQ2 대칭적인 두 번째 차수 시간 통합기법의 피아노형 조합이 작은 붕괴 시간 $\tau$에서 점근 보존성을 유지하는가?
• RQ3 $\tau \approx \Delta t$ 또는 $\tau < \Delta t$일 때, 복소수 시간 단계를 사용한 해법과 실수 단계 해법의 안정성과 정확도는 어떻게 비교되는가?
• RQ4 고 CFL 조건과 경직성 붕괴 하에서 4차 및 6차 해법의 수렴 특성은 어떠한가?
• RQ5 점성 시험 케이스에서 $\tau > 0$일 때, 특히 기존 방법이 실패하는 임계 값 근처에서 정확도와 안정성을 유지할 수 있는가?

주요 결과

• CFL 수치 100에서 저마하 수치 시험 케이스에서 4차 및 6차 수렴 속도가 최적으로 달성되었으며, 4차 수지 및 6차 카한-리 방법이 기대한 수렴 행동을 보였다.
• 높은 CFL 수치($\beta = 100$)에서 복소수 시간 단계를 사용한 해법은 불안정해지지만, 피아노형 해법은 안정적이고 정확하게 유지된다.
• $\tau = 0.000519$일 때, 6차 카한-리 방법은 특이 조건($2\tau + \gamma_2\Delta t/2 = 0$)으로 인해 실패하지만, 복소수 단계 해법은 안정적이며 $N_x=100$ 및 $N_x=1000$ 메esh에서 기준 해와 구분 불가능한 결과를 생성한다.
• 해법은 점근 보존성을 유지한다: $\tau \to 0$일 때 해는 원래의 쌍곡 보존 법칙으로 수렴하며, $\tau > 0$일 땐 채프먼-엔스코그 이론과 일치하는 점성 효과를 갖는다.
• 음성적 DG 해법은 삼각형 선형 시스템을 생성하며, 전진 대체를 통해 효율적으로 해석 가능하여 저메모리 구현과 높은 효율성을 가능하게 한다.
• 해법은 $\tau > 0$에서도 강건하며, 복소수 단계 해법은 중간 수준의 CFL 수치에서 $2\tau + \Delta t$ 안정성 장벽을 넘을 수 있어, 경직성 영역에서 더 넓은 적용 가능성을 시사한다.