[論文レビュー] On the stability of robust dynamical low-rank approximations for hyperbolic problems

Jonas Kusch, Lukas Einkemmer+1 more|arXiv (Cornell University)|Jan 1, 2021

Model Reduction and Neural Networks参考文献 35被引用数 2

ひとこと要約

本稿では、双曲型問題の空間的・時間的離散化の前に行われる動的低ランク近似（DLRA）を適用する、新たなプロジェクタースプリット統合手法を提案する。これにより古典的CFL条件が回復され、安定性が向上する。一次モーメントにおいて優れた安定性を示し、運動論的輸送における散乱項の効率的かつ安定な更新が可能になることが、数値実験で検証された。

ABSTRACT

The dynamical low-rank approximation (DLRA) is used to treat high-dimensional problems that arise in such diverse fields as kinetic transport and uncertainty quantification. Even though it is well known that certain spatial and temporal discretizations when combined with the DLRA approach can result in numerical instability, this phenomenon is poorly understood. In this paper we perform a L2 stability analysis for the corresponding nonlinear equations of motion. This reveals the source of the instability for the projector splitting integrator when first discretizing the equations and then applying the DLRA. Based on this we propose a projector splitting integrator, based on applying DLRA to the continuous system before performing the discretization, that recovers the classic CFL condition. We also show that the unconventional integrator has more favorable stability properties and explain why the projector splitting integrator performs better when approximating higher moments, while the unconventional integrator is generally superior for first order moments. Furthermore, an efficient and stable dynamical low-rank update for the scattering term in kinetic transport is proposed. Numerical experiments for kinetic transport and uncertainty quantification, which confirm the results of the stability analysis, are presented.

研究の動機と目的

• 双曲型問題の空間的・時間的離散化後にDLRAを適用した場合の数値不安定性の原因を解明すること。
• 空間的・時間的離散化の前に行われるDLRAを適用することで、古典的CFL条件を回復する安定化されたプロジェクタースプリット統合手法を開発すること。
• 新規統合手法と従来のプロジェクタースプリット法との安定性および精度を、一次モーメントおよび高次モーメントにおいて比較すること。
• 運動論的輸送方程式における散乱項に対して、効率的かつ安定な動的低ランク更新を設計すること。
• 理論的安定性の予測を、運動論的輸送および不確実性定量化の分野における数値実験で検証すること。

提案手法

• 離散化後にDLRAを適用した場合に生じる非線形運動方程式のL2安定性解析を行い、不安定性の原因を同定する。
• 空間的・時間的離散化の前に行われる連続系へのDLRAを適用する新しいプロジェクタースプリット統合手法を提案する。
• 早期の低ランク射影によって安定性を確保することで、古典的CFL条件を回復する修正された統合手法を導出する。
• 一次モーメントと高次モーメントの両者において、従来手法と非従来手法の安定性特性を分析・比較する。
• 運動論的輸送における散乱項に対して、安定性と効率性を維持する専用の動的低ランク更新スキームを設計する。
• 運動論的輸送および不確実性定量化問題における数値実験で、提案された統合手法を実装・テストする。

実験結果

リサーチクエスチョン

• RQ1 DLRAを空間的・時間的離散化後に適用した場合、プロジェクタースプリット統合手法に数値不安定性が生じる原因は何か？
• RQ2 離散化の前に行われる低ランク近似を適用することで、古典的CFL条件を回復するDLRAベースの統合手法を設計できるか？
• RQ3 新規統合手法と従来のプロジェクタースプリット法との間で、一次モーメントと高次モーメントの両方において安定性および精度はどのように比較されるか？
• RQ4 運動論的輸送方程式の散乱項に対して、効率的かつ安定な動的低ランク更新を構築できるか？
• RQ5 数値実験は、提案された統合手法の理論的安定性予測を確認できるか？

主な発見

• 従来のプロジェクタースプリット統合手法の不安定性は、離散化後にDLRAを適用することによって、安定性に必要な元の構造が破壊されることが原因である。
• 離散化の前に行われるDLRAを適用する提案された統合手法は、古典的CFL条件を正常に回復し、安定性を確保する。
• 非従来の統合手法は、全体的により優れた安定性特性を示し、特に一次モーメントにおいて顕著である。
• 高次モーメントの近似においては、従来のプロジェクタースプリット統合手法が、非線形運動方程式における特定の構造のおかげで優れた性能を示す。
• 運動論的輸送における散乱項に対して、効率的かつ安定な動的低ランク更新が正しく導出され、数値的に検証された。
• 運動論的輸送および不確実性定量化における数値実験により、理論的安定性解析の妥当性が確認され、提案された統合手法の優位性が示された。