[论文解读] Error analysis of an asymptotic preserving dynamical low-rank integrator for the multi-scale radiative transfer equation

Zhiyan Ding, Lukas Einkemmer+1 more|arXiv (Cornell University)|Jul 9, 2019

Sparse and Compressive Sensing Techniques参考文献 38被引用 7

一句话总结

本文首次对多尺度辐射转移方程（一种小Knudsen数的动能模型）的动力低秩（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格式在Knudsen数主导的多尺度区域内，对解的固有低秩结构的捕捉能力。
• 证明DLR格式在离散层面保持了辐射转移方程的流体极限，确保与渐近行为的一致性。
• 为动能理论中的DLR方法建立理论基础，此前的数值结果缺乏正式分析。

提出的方法

• 该方法将辐射转移方程的时间依赖解投影到由时变秩一因子参数化的低秩流形上。
• 采用动力低秩积分器，通过从控制PDE推导出的常微分方程系统直接演化低秩因子。
• 该格式设计用于在整个时间演化过程中保持低秩结构，从而在最小化自由度数量的同时保持精度。
• 通过将DLR解与精确解进行比较，开展误差分析，重点关注Knudsen数较小时的渐近区域。
• 分析利用了在流体极限下解的已知秩一结构，证明DLR方法在数值层面能精确捕捉该行为。

实验结果

研究问题

• RQ1 动力低秩积分器能否在多尺度区域内准确捕捉辐射转移方程解的低秩结构？
• RQ2 DLR格式是否在离散层面上保持了辐射转移方程的流体极限？
• RQ3 在小Knudsen数的动能方程背景下，DLR方法的理论误差界是什么？
• RQ4 与标准数值方法相比，DLR方法在该问题中结构保持性和计算效率方面表现如何？

主要发现

• 动力低秩积分器在渐近区域内成功捕捉了解的秩一结构，与辐射转移方程理论预测一致。
• 该方法在数值层面上保持了流体极限，即当Knudsen数趋于零时，离散解收敛到正确的流体方程。
• 误差分析表明，DLR方法在刚性区域仍能保持精度，而标准方法在此区域面临稳定性和效率挑战。
• 本工作首次为DLR方法在动能问题中观测到的数值性能提供了理论依据，证实了其鲁棒性和结构保持特性。