[논문 리뷰] Asymptotic preserving Implicit-Explicit Runge-Kutta methods for non linear kinetic equations

Giacomo Dimarco, Lorenzo Pareschi|arXiv (Cornell University)|2012. 05. 04.

Gas Dynamics and Kinetic Theory참고 문헌 30인용 수 25

한 줄 요약

이 논문은 강한 비선형 운동학 방정식, 특히 볼츠만 방정식에 대해 펜리티드 음성-음성(IMEX) 룬게-쿠타 스킴을 제안하며, 충돌 연산자의 비용이 큰 암시적 역행렬 계산 없이도 점근적 보존성과 정확성을 보장한다. 충돌 항을 보상 기법을 통해 평형 및 비평형 부분으로 분해함으로써, 다양한 수축 시간에 걸쳐 안정성과 효율성을 달성하였으며, 3차 스킴과 수치 실험을 통해 검증되었다.

ABSTRACT

We discuss Implicit-Explicit (IMEX) Runge Kutta methods which are particularly adapted to stiff kinetic equations of Boltzmann type. We consider both the case of easy invertible collision operators and the challenging case of Boltzmann collision operators. We give sufficient conditions in order that such methods are asymptotic preserving and asymptotically accurate. Their monotonicity properties are also studied. In the case of the Boltzmann operator, the methods are based on the introduction of a penalization technique for the collision integral. This reformulation of the collision operator permits to construct penalized IMEX schemes which work uniformly for a wide range of relaxation times avoiding the expensive implicit resolution of the collision operator. Finally we show some numerical results which confirm the theoretical analysis.

연구 동기 및 목표

• 스모올 쿠드센 수가 작은 유체역학적 극한에서의 강한 비선형 운동학 방정식에 대해 효율적인 수치 스킴을 개발하기 위해.
• 강한 영역에서 비선형 볼츠만 충돌 연산자의 직접 역행렬 계산을 피하여 암시적 해법의 계산 부담을 줄이기 위해.
• 보상 전략을 사용하여 IMEX 룬게-쿠타 방법을 전체 볼츠만 방정식에 확장하여 점근적 보존성 및 점근적 정확성 특성을 유지하기 위해.
• 특히 동질적 경우에서의 단조성과 안정성 분석 및 보장하기 위해.
• 강한 ODE 시스템에 적용 가능한 일반적인 프레임워크 제공을 위해, 특히 수축 연산자와 보존량을 포함한 대규모 시스템에 대해.

제안 방법

• 충돌 연산자를 $ R(Y) = N(Y) + L(Y) $ 로 분해하며, 여기서 $ L(Y) $ 는 평형 상태에서의 자코비안을 근사하는 선형화된 수축 항이다.
• 보상 기법을 도입하여 볼츠만 충돌 적분의 양성 항을 재구성함으로써, 비선형 부분은 명시적으로 처리하고 선형화된 수축 항은 암시적으로 처리할 수 있도록 한다.
• IMEX 스킴을 구성하여 암시적 단계에서 보상된 연산자를 포함하는 선형 시스템을 풀게 하여 전체 비선형 충돌 항의 역행렬 계산이 필요 없도록 한다.
• 기본 ODE 시스템의 점근적 분석에서 유도된 충분 조건을 만족시킴으로써 점근적 보존성을 확보한다.
• 고차 정확도와 강력한 안정성 특성을 갖춘 스킴을 설계하며, GSA(일반화된 단계 순서) 및 단조성 조건을 포함한다.
• 부처 테이블을 사용하여 특정 스킴을 구성하며, 일차, 이차, 삼차 방법을 포함하며, 매개변수 조정을 통해 단조성 간격을 명시적으로 제어한다.

실험 결과

연구 질문

• RQ1 전체 비선형 볼츠만 방정식에 적용할 때, IMEX 룬게-쿠타 스킴이 점근적 보존성과 점근적 정확성을 유지할 수 있는가?
• RQ2 강한 영역에서 암시적 충돌 해법의 계산 비용을 줄일 수 있는 방법은 무엇인가, 안정성이나 정확성을 희생하지 않고서도?
• RQ3 어떤 보상 전략이 충격 항의 효율적인 암시-명시적 통합을 가능하게 하며, 동시에 점근적 행동을 유지하는가?
• RQ4 어떤 조건에서 보상된 IMEX 스킴이 단조성을 유지하는가, 특히 동질적 경우에서?
• RQ5 3차까지의 고차 스킴을 구성할 수 있는가, 이는 점근적 보존성과 무조건적 안정성을 동시에 확보하는가?

주요 결과

• 보상된 IMEX 스킴은 부처 테이블과 보상 매개변수에 대한 충분 조건을 만족할 경우 전체 볼츠만 방정식에서 점근적 보존성과 점근적 정확성을 달성한다.
• 보상 기법을 통해 전체 비선형 충돌 연산자의 고비용 암시적 해법을 피할 수 있어 다양한 수축 시간 범위에서 효율적인 계산이 가능하다.
• 일차 스킴인 DP-A $(1,2,1)$ 과 DP-ARS $(1,2,1)$ 는 AA(점근적 정확성) 조건을 만족하며 $ ho \neq 1 $ 일 때 단조성을 확보할 수 있다.
• 이차 스킴인 DP2-A $(2,4,2)$ 는 DIRK(대각선 암시적 룬게-쿠타) 부분에서 3차 정확도를 달성하며, $ ho \to \rho_{\text{min}} $ 일 때 $ \rho \neq 1 $ 이면 단조성을 유지한다.
• 삼차 스킴인 DP1-A $(2,4,2)$ 는 네 단계로 구성되며 GSA 및 AA 조건을 모두 만족하여 강한 영역에서 고차 정확도를 달성한다.
• 수치 결과는 이론적 분석을 확인하며, $ \rho = 10^{-6} $ 를 포함한 다양한 $ \rho $ 값에서 수렴성과 안정성을 보이며, 유체역학적 극한을 정확하게 해석한다.