[논문 리뷰] On the uniform accuracy of implicit-explicit backward differentiation formulas (IMEX-BDF) for stiff hyperbolic relaxation systems and kinetic equations

Jingwei Hu, Ruiwen Shu|arXiv (Cornell University)|2019. 12. 02.

Gas Dynamics and Kinetic Theory참고 문헌 33인용 수 9

한 줄 요약

이 논문은 강한 비선형 초월 시스템과 운동 방정식에 대해 유도된 은폐-명시적 역행성 차분 공식(IMEX-BDF)의 일관된 정확도를 이론적으로 증명한다. 이는 이완 파arameter ε의 모든 영역에서 성립한다. 새로운 다중변수 기법을 사용한 독창적인 에너지 추정을 통해 저자들은 ε에 종속되지 않는 최적의 사전 오차 경계를 확립하여, 기존 방법이 자주 순서 감소 문제를 겪는 중간 영역에서도 고차 정확도를 보장한다.

ABSTRACT

Many hyperbolic and kinetic equations contain a non-stiff convection/transport part and a stiff relaxation/collision part (characterized by the relaxation or mean free time $\varepsilon$). To solve this type of problems, implicit-explicit (IMEX) multistep methods have been widely used and their performance is understood well in the non-stiff regime ($\varepsilon=O(1)$) and limiting regime ($\varepsilon ightarrow 0$). However, in the intermediate regime (say, $\varepsilon=O(Δt)$), uniform accuracy has been reported numerically without a complete theoretical justification (except some asymptotic or stability analysis). In this work, we prove the uniform accuracy -- an optimal {\it a priori} error bound -- of a class of IMEX multistep methods, IMEX backward differentiation formulas (IMEX-BDF), for linear hyperbolic systems with stiff relaxation. The proof is based on the energy estimate with a new multiplier technique. For nonlinear hyperbolic and kinetic equations, we numerically verify the same property using a series of examples.

연구 동기 및 목표

• 수치적 증거는 있었지만 완전한 이론적 근거가 없었던 중간 영역(ε = O(Δt))에서 IMEX-BDF 스킴의 일관된 정확도에 대한 엄밀한 이론적 기반을 마련하는 것.
• ε가 작지도 않고 O(1)이지도 않은 경우, 강한 초월 및 운동 방정식에 대한 IMEX 스킴에서의 순서 감소 현상에 대한 이론적 이해의 격차를 메우는 것.
• 초월 시스템의 대류-이완 구조에 특화된 새로운 다중변수 기법을 개발하고 적용하여 안정성과 수렴성을 증명하는 것.
• 선형 경우에서 관찰된 것과 동일한 일관된 정확도 행동을 보장하기 위해 비선형 운동 방정식에 분석을 확장하고 수치적 검증을 수행하는 것.

제안 방법

• 저자들은 강한 초월 이완 시스템의 구조에 특화된 새로운 다중변수 기법을 적용한 에너지 추정을 사용하여, 이전 다중변수 기법의 한계를 극복한다.
• 이 방법은 강한 이완을 포함한 선형 초월 시스템, 예를 들어 식 (1.1)의 모델 시스템에 적용되며, 계수 변화가 있는 경우로도 확장된다.
• 운동 방정식의 경우, BGK 모델을 사용하여 수치적으로 검증되며, IMEX-BDF 스킴은 비선형 시스템의 모멘트 기반 해법을 통해 구현된다.
• 고차 정확도를 달성하기 위해 스펙트럼 공간 이산화 기법이 스킴과 결합되어 엄밀한 오차 분석이 가능해진다.
• 증명 과정에서 ε에 종속되지 않는 최적의 사전 오차 경계 ‖u(T,·)−U(T,·)‖L² + ‖v(T,·)−V(T,·)‖L² ≤ CΔt^q 가 확립된다.
• 기호 계산을 통해 q=4인 IMEX-BDF 스킴의 고차 계수를 도출하였으며, 안정성과 일致성 조건을 매우 정밀하게 만족시키도록 한다.

실험 결과

연구 질문

• RQ1 강한 초월 이완 시스템에 대해 IMEX-BDF 스킴의 일관된 정확도가 ε의 모든 값에서 엄밀히 증명될 수 있는가?
• RQ2 기존 다중변수 기법이 실패하는 가운데, 중간 영역에서 IMEX-BDF 스킴의 안정성을 확보하기 위해 어떤 새로운 수학적 기법이 필요한가?
• RQ3 IMEX-BDF 스킴의 일관된 정확도 성질은 BGK 모델과 같은 비선형 운동 방정식으로까지 확장되는가?
• RQ4 ε가 시간 단계 Δt와 같은 주기일 때, IMEX-BDF 스킴의 정확도 순서는 어떻게 행동하며, 그 이론적 근거는 무엇인가?

주요 결과

• 논문은 선형 초월 이완 시스템에 적용된 IMEX-BDF 스킴에 대해, 상수 C가 ε에 종속되지 않는 최적의 사전 오차 경계 ‖u(T,·)−U(T,·)‖L² + ‖v(T,·)−V(T,·)‖L² ≤ CΔt^q 를 증명한다.
• 이중 정확도는 4차까지의 IMEX-BDF 스킴에 대해 확립되며, ε가 비강성, 강성, 또는 중간 영역에 있든 간에 고차 수렴성을 보장한다.
• 대류-이완 구조를 정확히 반영할 수 있는 새로운 다중변수 기법이 개발되어, 일관된 안정성과 수렴성의 증명이 가능해졌다.
• 수치 결과는 비선형 운동 방정식, 예를 들어 BGK 방정식에서도 동일한 일관된 정확도 성질이 유지됨을 확인한다. 이는 스킴의 비선형 은밀한 성격에도 불구하고 성립한다.
• 분석 결과, 1차 스킴는 ε가 작을 경우 해에서 O(ε) 항을 포착하지 못함을 보여주며, 중간 척도의 역학을 해석하기 위해 고차 스킴의 필요성을 강조한다.
• q=4인 경우, 저자들은 계수를 10^-31 정밀도로 계산하였으며, 관련 행렬 G와 A의 최소 고유값은 각각 약 5.03×10^-6과 0.0765로 확인되어 수치적 안정성이 확인된다.