[논문 리뷰] Error analysis of an asymptotic preserving dynamical low-rank integrator for the multi-scale radiative transfer equation

Zhiyan Ding, Lukas Einkemmer외 1명|arXiv (Cornell University)|2019. 07. 09.

Sparse and Compressive Sensing Techniques참고 문헌 38인용 수 7

한 줄 요약

이 논문은 소수의 쿠드센 수를 가진 다중 척도 방사열전달 방정식, 즉 운동학적 모델에 대해 동적 저질서(DLR) 적분기의 첫 번째 수학적 오차 분석을 제시한다. 이 방법은 시간에 따라 변화하는 질서 인자로 해를 저질서 다양체 위로 투영하며, 해의 질서 일치 구조를 동적으로 포착하고 유체 한계를 유지함으로써 운동학 문제에 대한 이론적 수렴성과 강건성을 입증한다.

ABSTRACT

Dynamical low-rank algorithm are a class of numerical methods that compute low-rank approximations of dynamical systems. This is accomplished by projecting the dynamics onto a low-dimensional manifold and writing the solution directly in terms of the low-rank factors. The approach has been successfully applied to many types of differential equations. Recently, efficient dynamical low-rank algorithms have been applied to treat kinetic equations, including the Vlasov--Poisson and the Boltzmann equation, where it was demonstrated that the methods are able to capture the low rank structure of the solution and significantly reduce the numerical effort, while often maintaining good accuracy. However, no numerical analysis is currently available. In this paper, we perform an error analysis for a dynamical low-rank algorithm applied to a classical model in kinetic theory, namely the radiative transfer equation. The model used here includes a small parameter, the Knudsen number. This setting is particularly interesting since the solution is known to be rank one in certain regimes. We will prove that the scheme dynamically and automatically captures the low-rank structure of the solution, and preserves the fluid limit on the numerical level. This work thus serves as the first mathematical error analysis for a dynamical low rank approximation applied to a kinetic problem.

연구 동기 및 목표

• 운동학 방정식, 특히 방사열전달 방정식에 적용된 동적 저질서(DLR) 방법에 대한 첫 번째 엄밀한 수학적 오차 분석을 제공하는 것.
• 쿠드센 수에 의해 지배되는 다중 척도 영역에서 DLR 스킴이 해의 본질적인 저질서 구조를 얼마나 잘 포착하는지 분석하는 것.
• DLR 스킴이 이산 수준에서 유체 한계를 유지함으로써 점점 감소하는 쿠드센 수에서의 점근적 거동과 일致함을 보여주는 것.
• 이전의 수치 결과가 공식적 분석을 결여한 운동학 이론에서 DLR 방법의 이론적 기초를 구축하는 것.

제안 방법

• 시간에 따라 변화하는 방사열전달 방정식의 해를 시간에 따라 변화하는 질서 일치 요소로 매개변수화된 저질서 다양체 위로 투영하는 방법.
• 지배 PDE로부터 유도된 상미분 방정식 시스템을 통해 직접 저질서 인자를 진화시키는 동적 저질서 적분기를 사용하는 것.
• 시간 진화 전반에 걸쳐 저질서 구조를 유지하도록 설계되어, 자유도 수를 최소화하면서도 정확도를 유지하는 것.
• 쿠드센 수가 작을 때의 점근적 영역에 중점을 두고, DLR 해를 정확한 해와 비교하여 오차 분석을 수행하는 것.
• 유체 한계에서 해의 알려진 질서 일치 구조를 활용하여, DLR 방법이 이 행동을 수치 수준에서 정확히 포착함을 증명하는 것.

실험 결과

연구 질문

• RQ1 동적 저질서 적분기는 다중 척도 영역에서 방사열전달 방정식 해의 저질서 구조를 정확히 포착할 수 있는가?
• RQ2 DLR 스킴은 이산 수준에서 방사열전달 방정식의 유체 한계를 유지하는가?
• RQ3 소수의 쿠드센 수를 가진 운동학 방정식의 맥락에서 DLR 방법의 이론적 오차 한계는 무엇인가?
• RQ4 이 문제에 대해 기존 수치 방법과 비교했을 때 DLR 접근법은 구조 유지 및 계산 효율성 측면에서 어떻게 다른가?

주요 결과

• 동적 저질서 적분기는 이론적으로 예측된 것과 같이 점점 감소하는 영역에서 해의 질서 일치 구조를 정확히 포착한다.
• 이 방법은 수치 수준에서 유체 한계를 유지하므로, 쿠드센 수가 0으로 수렴할 때 이산 해가 올바른 유체 방정식으로 수렴함을 의미한다.
• 오차 분석을 통해 DLR 방법이 표준 방법이 안정성과 효율성 문제를 겪는 강성 영역에서도 정확도를 유지함을 입증한다.
• 이 작업은 DLR 방법이 운동학 문제에서 관찰된 수치 성능에 대한 첫 이론적 근거를 제공하며, 그 강건성과 구조 유지 성격을 확인한다.