[論文レビュー] Generalised Lebesgue Stable Flux Reconstruction

Will Trojak|arXiv (Cornell University)|May 31, 2018

Computational Fluid Dynamics and Aerodynamics参考文献 24被引用数 3

ひとこと要約

本稿では、エネルギー安定性をルベーグノルムで保証する一般化されたリーマン安定フラックス再構成（GLSFR）という、フラックス再構成法のための補正関数の新規族を提案する。これにより、任意の構成が可能でありながらも安定なスキームを実現する。バナッハ・ニューマン解析および乱流ターレイ・グリーン渦シミュレーションを通じて、GLSFRは運動エネルギー散逸誤差においてノード的DGより優れた精度を達成でき、最適な補正関数では、メッシュを細分化した際の誤差を最大30%まで低減できることが示された。

ABSTRACT

A unique set of correction functions for Flux Reconstruction is presented, with there derivation stemming from proving the existence of energy stability in the Lebesgue norm. The set is shown to be incredibly arbitrary with the only union to existing correction function sets being show to be for DG. Von Neumann analysis of both advection and coupled advection-diffusion is used to show that once coupled to a temporal integration method, good CFL performance can be achieved and the correction function may have better dispersion and dissipation for application to implicit LES. Lastly, the turbulent Taylor-Green vortex test case is then used to show that correction functions can be found that improve the accuracy of the scheme when compared to the error levels of Discontinuous Galerkin.

研究の動機と目的

• フラックス再構成（FR）法における安定な補正関数の理解が不完全である、特に標準的でない、またはDGでない構成に対して、その問題を解決すること。
• ルベーグノルムにおけるエネルギー安定性を保証する新しい補正関数の族を導出すること。既存の枠組みを超えて拡張すること。
• これらの任意の補正関数が、ルンゲ・クッタ時間積分法と組み合わせた際に質量保存および時間的安定性を維持できるかを検証すること。
• GLSFRの実用的性能を乱流シミュレーションにおいて評価し、特に運動エネルギー散逸の正確な捉え込みを可能にするかを検討すること。
• GLSFRが実世界の流体力学的応用において、標準的なノード的DGをFR法で上回る可能性があるかを特定すること。

提案手法

• ソボレフノルムではなく、ルベーグノルムを用いてFRスキームのエネルギー安定性を導出することで、安定性の証明を簡略化し、より広範な補正関数の族を可能にする。
• 境界条件 $ h_L(-1)=1, h_L(1)=0 $ および $ h_R(-1)=0, h_R(1)=1 $ を満たす補正関数 $ h_L $ と $ h_R $ を定義し、ルベーグ安定性を保証する。
• ルンゲ・クッタ時間積分法を用いた線形移流および移流拡散方程式に対して、バナッハ・ニューマン解析を適用し、時間的安定性および分散・散逸特性を評価する。
• 乱流非圧縮性ターレイ・グリーン渦を $ Re = 1600 $ で用い、時間平均された運動エネルギー散逸率の誤差をDNSデータと比較して評価する。
• 補正関数のパラメータ（例：$ h_{l1} $、$ ilde{h}_{l0} $）を体系的に変化させ、散逸誤差を最小化する最適な設定を同定する。
• GLSFRの結果を、同一のグリッド解像度および時間ステップを用いたノード的DGによるFRとの比較を通じて評価する。

実験結果

リサーチクエスチョン

• RQ1 FR法において、既存の既知の集合を超えて、より広いクラスの補正関数に対して、ルベーグノルムにおけるエネルギー安定性を確立できるか？
• RQ2 導出された補正関数の族は、明示的ルンゲ・クッタ時間積分法と組み合わせた際にも質量保存および時間的安定性を維持できるか？
• RQ3 GLSFR補正関数は、高次のシミュレーションにおいて、ノード的DGによるFRと比較して、乱流運動エネルギー散逸の誤差を低減できるか？
• RQ4 特定の補正関数パラメータは、暗黙のLESに適した分散および散逸特性を改善するものか？
• RQ5 GLSFRの性能向上は、乱流問題における異なるグリッド解像度において一貫しているか？

主な発見

• GLSFRフレームワークは、既存の安定集合を超える、ルベーグ安定性を満たす補正関数の新規族を確立した。唯一の共通点は、FR法によるノード的DGである。
• バナッハ・ニューマン解析により、ルンゲ・クッタ時間積分法を用いた場合、上流化および中央フラックス界面の両方で時間的安定性領域が確認され、実用的応用が可能であることが示された。
• 最適なGLSFR補正関数は、$ 24^3 $ グリッド上でのノード的DGによるFRと比較して、乱流運動エネルギー散逸誤差を最大30%まで低減した。
• $ 16^3 $ および $ 24^3 $ グリッドの両方において、GLSFR補正関数は、誤差最小化においてノード的DGを上回る一貫したパラメータ空間領域に位置していた。
• 最適な補正関数は、$ ilde{h}_{l0} $ の左半平面に位置しており、誤差を最小化するためには、より散逸の少ないスキームが好ましいことが示された。これはCFL安定性の傾向とは逆である。
• 渦構造が形成される臨界遷移段階（$ t \approx 7.5 $）において、GLSFRはDGよりも運動エネルギー散逸率をよりよく再現しており、特に高解像度条件下で顕著であった。