[論文レビュー] Analysis of a low complexity, time-accurate discretization of the Navier-Stokes equations

Victor DeCaria, William Layton+1 more|arXiv (Cornell University)|Oct 15, 2018

Numerical Methods and Algorithms参考文献 1被引用数 2

ひとこと要約

本稿では、標準的な後退オイラー離散化に2段階線形時間フィルタを適用することにより、ナビエ＝ストークス方程式に対して低複雑性で時間に正確な手法を提案する。この手法は計算コストを増加させることなく精度を向上させ、過剰な減衰を除去する一方で、無条件エネルギー安定性を維持する。数値実験では、抗力および揚力の予測が改善される。

ABSTRACT

This report presents a low computational and cognitive complexity, stable, time accurate and adaptive method for the Navier-Stokes equations. The improved method requires a minimally intrusive modification to an existing program based on the fully implicit / backward Euler time discretization, does not add to the computational complexity, and is conceptually simple. The backward Euler approximation is simply post-processed with a two-step, linear time filter. The time filter additionally removes the overdamping of Backward Euler while remaining unconditionally energy stable, proven herein. Even for constant stepsizes, the method does not reduce to a standard / named time stepping method but is related to a known 2-parameter family of A-stable, two step, second order methods. Numerical tests confirm the predicted convergence rates and the improved predictions of flow quantities such as drag and lift.

研究の動機と目的

• ナビエ＝ストークス方程式に対して、時間に正確で、安定的かつ計算効率の良い時間離散化手法を開発すること。
• 既存の後退オイラー・ソルバを変更する際の認知的および計算的オーバーヘッドを最小限に抑えること。
• 後退オイラー法に内在する過剰な減衰を排除しつつ、無条件エネルギー安定性を保持すること。
• 抗力や揚力などの重要な流れの物理量の予測精度を向上させること。
• 実装において既存のコードベースへの影響を最小限に抑え、シンプルに統合できること。

提案手法

• 2段階線形時間フィルタを、後退オイラー解の後処理として適用する。
• フィルタは数値的減衰を低減しつつ安定性を保持するように時間離散化を変更する。
• 完全に陰的な後退オイラースキームの修正として導出され、システム行列の構築に変更を加える必要がない。
• フィルタは無条件エネルギー安定性であると、論文で証明されている。
• 2パラメータを持つ、A安定性を満たし、2次精度を持つ2段階法の既知の族に関連している。
• 最小限のコード変更で実装可能であり、レガシーソルバへの統合に適している。

実験結果

リサーチクエスチョン

• RQ1 単純な後処理フィルタを用いることで、計算コストを増加させることなく後退オイラー法の精度を向上させることができるか？
• RQ2 フィルタを適用した手法は、無条件エネルギー安定性を維持しつつ、数値的過剰減衰を低減するのか？
• RQ3 この手法は、流れシミュレーションにおける抗力および揚力の予測をどの程度改善するのか？
• RQ4 収束性および安定性の観点から、標準的な2次精度時間積分法と比較して、この手法はどのように異なるのか？
• RQ5 既存の後退オイラー・コードベースに最小限の変更で実装可能か？

主な発見

• フィルタを適用した手法は、定常時間ステップでも時間方向に2次精度を達成するが、標準的な名前が付いた方法に帰着しない。
• 数値実験により、時間離散化の予測された収束率が確認された。
• 標準的な後退オイラー法と比較して、抗力および揚力係数の計算精度が顕著に向上した。
• 時間フィルタは、後退オイラー法に内在する過剰な数値的減衰を効果的に除去したが、無条件エネルギー安定性を維持した。
• 論文で厳密に証明されたように、手法は無条件エネルギー安定性を保ったままである。
• 計算複雑性に増加はなく、既存のコードベースへの影響も最小限に抑えられる。