[論文レビュー] Dispersion-dissipation condition for finite difference schemes

Xiao Hu, V. K. Tritschler+2 more|arXiv (Cornell University)|Apr 23, 2012

Computational Fluid Dynamics and Aerodynamics参考文献 10被引用数 26

ひとこと要約

本稿では、有限差分スキームにおける不適切な高波数波を減衰させるために必要な最小限の数値的散乱を決定する分散・散乱条件を導入する。修正波数解析を通じて群速度と散乱率を結びつけることで、線形および非線形スキームにおける分散と散乱を最適化する。この条件は、WENO-CU6-M2スキームの改良に応用され、乱流流れの過小解像度シミュレーション、特に自己相似的エネルギー減衰および乱流への遷移の高精度な再現を可能にする。

ABSTRACT

A general dispersion-dissipation condition for finite difference schemes is derived by analyzing the numerical dispersion and dissipation of explicit finite-difference schemes. The proper dissipation required to damp spurious high-wavenumber waves in the solution is determined from a physically motivated relation between group velocity and dissipation rate. The application to a previously developed low-dissipation weighted essentially non-oscillatory scheme (WENO-CU6-M2) [X. Y. Hu and N. A. Adams; Scale separation for implicit large eddy simulation, J. Comput. Phys. 230 (2011) 7240-7249] demonstrates that this condition can serve as general guideline for optimizing the dispersion and dissipation of linear and non-linear finite-difference schemes. Moreover, the improved WENO-CU6-M2 which satisfies the dispersion-dissipation condition can be used for under-resolved simulations. We demonstrate this capability by considering transition to turbulence and self-similar energy decay of the three-dimensional Taylor-Green vortex. Simulations of the inviscid and the viscous Taylor-Green vortex at Reynolds numbers ranging from $Re=400$ to $Re=3000$ show a significant improvement over the classical dynamic Smagorinsky model and demonstrate competitiveness with state-of-the-art implicit LES models, while preserving shock-capturing properties.

研究の動機と目的

• 有限差分スキームにおける分散と散乱の間の定量的関係を確立し、不適切な高波数波を制御する。
• 物理的でない波を減衰させる最適な散乱レベルを特定する長年の課題に取り組むが、解像されたスケールの精度を損なわないようにする。
• 線形および非線形スキームの両方における分散と散乱の最適化の一般的な指針を提供するが、特に過小解像度シミュレーションを想定する。
• WENO-CU6-M2スキームに分散・散乱条件を満たすように改良し、乱流流れの高精度な大渦シミュレーションを可能にする。

提案手法

• 線形移流方程式を用いた線形有限差分スキームの修正波数から分散・散乱条件を導出する。
• 位相誤差（|dξ̃R/dξ − 1|）と散乱率（−ξ̃I）の比を用いて条件を定義し、限界的に解像された波の十分な減衰を保証する。
• WENO-CU6-M2スキームにこの条件を適用するため、再構成手順を変更し、高解像度特性を損なわずに散乱を調整する。
• スペクトル解析を用いて、特にナイキストカットオフ近辺の波数領域における分散および散乱特性を評価する。
• ベンチマーク問題（Sodの衝撃波管、2つの爆風波の衝突、Taylor-Green渦）を用いて、修正されたスキームを検証する。
• 精度の評価として、DNSデータおよび既存のLESモデル（動的Smagorinsky、ALDM）と比較し、エネルギー減衰および散乱率の両面で性能を評価する。

実験結果

リサーチクエスチョン

• RQ1 有限差分スキームにおいて、不適切な高波数波を効果的に減衰させるために必要な最小限の数値的散乱はどれくらいか？
• RQ2 解像されたスケールの精度を保ちながら非物理的振動を抑えるために、分散と散乱をどのようにバランスさせるか？
• RQ3 分散・散乱条件は、WENO-CU6-M2のような非線形で衝撃波捕捉を目的としたスキームに適用可能か？
• RQ4 改良されたWENO-CU6-M2スキームは、古典的なLESモデルと比較して、過小解像度乱流流れのシミュレーションでどのように性能を発揮するか？
• RQ5 分散・散乱条件を満たすことで、Taylor-Green渦における自己相似的エネルギー減衰の予測精度はどの程度向上するか？

主な発見

• 分散・散乱条件を満たすように改良されたWENO-CU6-M2スキームは、三次元非粘性Taylor-Green渦の自己相似的エネルギー減衰を正確に捉えており、理論的t−1.2減衰率と一致する。
• Re=400からRe=3000の粘性Taylor-Green渦シミュレーションにおいて、古典的動的Smagorinskyモデルと比較して、散乱率の予測において顕著な改善を示す。
• Re=400およびRe=800の条件下で、スキームの散乱率は、最先端のimplicit LESモデル（ALDM：適応的局所脱結合モデル）と同等の水準に達する。
• Kolmogorovの慣性領域におけるk−5/3スケーリングがカットオフ波数まで良好に保持されており、解像されたスケールにおける数値的散乱が最小限に抑えられていることを示している。
• Sodの衝撃波管および2つの爆風波の衝突問題を通じて、不適切な波の伝播を低減させつつ、正確な衝撃波捕捉特性を維持していることが確認された。
• 分散・散乱条件により、線形および非線形スキームにおける数値的散乱の体系的最適化が可能となり、乱流流れの高精度シミュレーションのための一般的なフレームワークを提供する。