[論文レビュー] Projected dynamics of kinetic equations with energy diffusion in spaces of orthogonal polynomials

Jon Wilkening, Antoine Cerfon+1 more|arXiv (Cornell University)|Feb 12, 2014

Model Reduction and Neural Networks被引用数 4

ひとこと要約

本論文では、速度空間におけるエネルギー拡散を伴う運動方程式の射影に、非古典的重み関数に関して直交する非古典的直交多項式を用いる手法を提案する。著者らは、やや特異な初期条件に対して、同じ計算コストで古典的エルミート多項式と比較して、数個の桁数も高い精度が得られることを示している。また、正確な解は丸め誤差の範囲内で射影部分空間に戻ることを確認しており、これはエルミート多項式を用いた場合とは異なり、大きなずれが生じないことを示している。

ABSTRACT

We examine a partial differential equation describing energy diffusion in velocity space, analyzing its projected dynamics in finite-dimensional subspaces of orthogonal polynomials. Comparison with an expensive but highly accurate spectral transform approach reveals that for a large class of (mildly singular) initial conditions, using polynomials that are orthogonal with respect to a non-classical weight function can be many orders of magnitude more accurate than classical Hermite polynomials for the same computational work. For the new polynomials, the exact solution is found to leave the subspace for a time, but return (up to roundoff accuracy) to the same point evolved to by the projected dynamics in that time. By contrast, using classical polynomials, the exact solution differs significantly from the projected dynamics solution when it arrives in the subspace. We also explore the connection between eigenfunctions of the projected evolution operator and (non-normalizable) eigenfunctions of the full evolution operator, as well as the effect of truncating the computational domain. Techniques for suppressing floating-point arithmetic errors are discussed, including use of the singular value decompositon to solve a better-conditioned version of the stiffness-mass generalized eigenvalue problem.

研究の動機と目的

• 速度空間におけるエネルギー拡散を伴う運動方程式の数値解の精度を、非古典的直交多項式を用いることで向上させること。
• 時間経過とともに顕著な誤差を示す古典的エルミート多項式の限界を是正すること。
• 射影系のダイナミクスを解析し、特に部分空間への再帰性と安定性に注目して、正確な解と比較すること。
• 射影演算子の固有関数と、完全な発展演算子の固有関数との関係を調査すること。
• スティフネス-質量一般固有値問題における浮動小数点演算誤差を特異値分解を用いて低減すること。

提案手法

• 非古典的重み関数から導かれる直交多項式が張る有限次元部分空間へ運動方程式を射影すること。
• 解の精度を評価するために、高精度なスペクトル変換法と射影されたダイナミクスを比較すること。
• スティフネス-質量一般固有値問題をより良好な数値的性質を持つ形に再定式化するために特異値分解を用いることで、数値誤差を低減すること。
• 時間経過に伴い正確な解が射影部分空間から離脱した後、再びその部分空間に戻る様子を解析し、特にその再帰性を評価すること。
• 射影された発展演算子のスペクトル的性質を解析し、それらを完全な演算子の非正規化可能な固有関数と関連付けること。
• 計算領域を切り詰める処理を実施し、その解の精度と安定性への影響を評価すること。

実験結果

リサーチクエスチョン

• RQ1 非古典的直交多項式は、エネルギー拡散方程式を解くにあたり、同じ計算コストで古典的エルミート多項式よりも著しく高い精度を達成できるか？
• RQ2 非古典的多項式を用いた場合、正確な解が一時的に部分空間を離脱した後、丸め誤差の範囲内で元の部分空間に戻るのに対し、エルミート多項式を用いた場合には戻らないのはなぜか？
• RQ3 射影された発展演算子の固有関数は、完全な発展演算子の非正規化可能な固有関数とどのように関係しているか？
• RQ4 計算領域の切り詰めが、射影解の精度と安定性に与える影響は何か？
• RQ5 スティフネス-質量一般固有値問題における浮動小数点演算誤差は、どのようにして効果的に抑制できるか？

主な発見

• やや特異な初期条件を満たす広範なクラスに対して、非古典的直交多項式を用いることで、同じ計算コストで古典的エルミート多項式を用いた場合と比較して、解の精度が数桁も向上することが確認された。
• 非古典的多項式を用いた射影において、正確な解は一時的に部分空間を離脱した後、丸め誤差の範囲内で元の部分空間に戻ることを確認した。これは、射影の忠実性が非常に高いことを示している。
• 対照的に、エルミート多項式を用いた場合、正確な解は部分空間に戻った後も、射影されたダイナミクスの解から著しくずれが生じ、近似の質が低いことが示された。
• 射影された演算子の固有関数は、完全な発展演算子の非正規化可能な固有関数と密接に関係しており、この手法がスペクトル的一致性を有していることが示唆された。
• スティフネス-質量一般固有値問題に対して特異値分解を用いて前処理を行うことで、浮動小数点演算誤差が効果的に抑制され、数値的安定性が向上した。
• 計算領域の切り詰めは誤差を生じるが、非古典的多項式アプローチは古典的代替手法と比較して、依然として頑健で高い精度を維持していた。