[論文レビュー] A time-accurate, adaptive discretization for fluid flow problems

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

Advanced Numerical Methods in Computational Mathematics参考文献 30被引用数 36

ひとこと要約

本稿では、ナビエ＝ストークス方程式に対する最小限で適応可能な時間離散化手法を提案する。この手法は、後退オイラー法に二段階線形時間フィルタを適用することで、精度を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次収束を証明すること。
• 特に中程度の時間刻みにおいて、標準の後退オイラー法と比較して、抗力、揚力、圧力損失といった重要な流れの量の精度が向上することを示すこと。
• 円柱の周りの流れといったベンチマーク問題を用いて、後退オイラー法が失敗するのに対し、正しい渦剥離行動を再現できることを検証すること。

提案手法

• 本手法は、後退オイラー解に二段階線形時間フィルタを適用することで、追加の解法や関数評価なしに1次精度の解を2次精度の解に変換する。
• フィルタは $ u^{n+1} = \hat{u}^{n+1} - \frac{1}{3}(\hat{u}^{n+1} - 2u^n + u^{n-1}) $ として適用され、ここで $ \hat{u}^{n+1} $ は後退オイラー解である。
• 2つの選択肢を提供する：オプションAは速度のみをフィルタリング、オプションBは速度と圧力の両方をフィルタリングし、オプションBは圧力の精度を向上させる。
• 本手法は埋め込み型でモジュラー設計であり、後退オイラーコードへの統合が容易で、適応的時間刻みと次数制御のための約20行の追加コードで実現可能である。
• 時間フィルタは無条件エネルギー安定性を保ち、特に大きな時間刻みにおいて典型的な後退オイラーの過剰な分散を回避する。
• 一定時間刻みの下で、速度の安定性と誤差解析を用いて手法を分析し、先行研究ではカバーされていなかった新しい圧力誤差解析を実施した。

実験結果

リサーチクエスチョン

• RQ1 後退オイラー法に対する最小限の修正、具体的には時間フィルタの適用によって、計算複雑性を増加させずに2次精度と改善された動的挙動を達成できるか？
• RQ2 時間フィルタを施した後退オイラー法は、流体流れのシミュレーションにおいて過剰な数値的減衰を低減しながらも、無条件エネルギー安定性を維持できるか？
• RQ3 時間フィルタは、抗力、揚力、圧力損失といった時間に依存する流れの量の予測を、標準の後退オイラー法と比較してどの程度向上させるか？
• RQ4 このフィルタベースのフレームワークを用いて、従来の後退オイラーコードに適応的時間刻みと適応的次数選択を容易に実装できるか？
• RQ5 速度のみをフィルタリングするのと、速度と圧力の両方をフィルタリングするのとでは、ナビエ＝ストークス方程式における圧力近似の精度にどのような差が生じるか？

主な発見

• 時間フィルタを施した後退オイラー法は、速度および圧力の両方で2次収束を達成しており、数値実験により両者とも最適な収束率が確認された。
• 円柱の周りの流れの問題において、後退オイラー法は $ \Delta t = 0.04 $ では渦剥離を捉えられなかったが、フィルタを施した手法は $ \Delta t = 0.02 $ から正しく予測した。
• 最大抗力誤差は、後退オイラー法（$ \Delta t = 0.04 $）の 2.95112558 から、フィルタを施した後退オイラー法に低下し、2.95021463 にまで減少し、基準値 2.950921575 に近づいた。
• フィルタを施した手法は、最大揚力の大きさとタイミングを正しく予測した：$ c_{l,max} = 0.4744 $ が $ t = 5.7 $ で得られ、基準値 0.47795 に近く、後退オイラー法は振幅と位相の両方を低く予測していた。
• 速度と圧力の両方をフィルタリングすると圧力損失の精度がわずかに向上するが、主な効果は流れの力学的誤差の主要因を是正する速度のフィルタリングに起因する。
• 本手法は無条件エネルギー安定性を保ち、時間刻みが十分に小さくない場合でも安定性を維持する。これに対して、標準の後退オイラー法は物理的振動を抑圧してしまう。