[論文レビュー] Semi-Lagrangian Finite-Element Exterior Calculus for Incompressible Flows

Wouter Tonnon, Ralf Hiptmair|arXiv (Cornell University)|Jan 12, 2023

Computational Fluid Dynamics and Aerodynamics被引用数 1

ひとこと要約

本稿では、単体メッシュ上での非圧縮性ナビエ＝ストークス方程式のための半ラグランジュ型有限要素外微分計算（FEEC）スキームを提示する。運動量方程式は速度1形式に対する非線形輸送問題に再定式化される。本手法は空間的・時間的2次精度を達成し、粘性係数がゼロに近づく極限でも優れた安定性を示し、ラグランジュ乗数を用いてエネルギーを保存する。これにより、構造保存性を有する特性を持つ非粘性エーラー流の高精度なシミュレーションが可能となる。

ABSTRACT

We develop a mesh-based semi-Lagrangian discretization of the time-dependent incompressible Navier-Stokes equations with free boundary conditions recast as a non-linear transport problem for a momentum 1-form. A linearly implicit fully discrete version of the scheme enjoys excellent stability properties in the vanishing viscosity limit and is applicable to inviscid incompressible Euler flows. Conservation of energy and helicity are enforced separately.

研究の動機と目的

• 任意の単体メッシュ上での有限要素外微分計算（FEEC）を用いた非圧縮性ナビエ＝ストークス方程式の構造保存型離散化を開発すること。
• 高い安定性とエネルギー保存性を備えた非粘性非圧縮性流れのシミュレーションという課題に取り組むこと。
• 半ラグランジュ法を微分形式、特に運動量1形式に拡張し、幾何的に一貫した輸送を実現すること。
• 空間的・時間的両方で2次精度を達成するとともに、非圧縮性およびエネルギーといった重要な物理的不変量を保存すること。

提案手法

• 微分形式および外微分計算を用いて、非圧縮性ナビエ＝ストークス方程式を運動量1形式 ω に対する非線形輸送問題に再定式化する。
• 流れに起因する引き戻しに基づく、運搬された1形式スナップショットの後退差分商を用いて、物質微分 Duω を離散化する。
• 固定された単体メッシュ上にガラーキン型有限要素法を適用し、空間離散化にホッジ分解および離散的外微分計算を用いる。
• 非粘性（ϵ=0）の場合にエネルギー保存を強制するためにラグランジュ乗数を導入し、構造保存性を確保する。
• 自由度を豊かにするために「小さな辺（small edges）」を有限要素空間に導入し、近似の安定性および精度を向上させる。
• 粘性係数がゼロに近づく極限でも安定性を保証する線形的陰解法を用いた時間積分スキームを実装する。

実験結果

リサーチクエスチョン

• RQ1 一般の単体メッシュ上での半ラグランジュ型有限要素外微分計算スキームは、非圧縮性流れに対して2次収束を達成できるか？
• RQ2 本手法は、粘性係数がゼロに近づく極限、特に非粘性エーラー流れに対してどのように性能を発揮するか？
• RQ3 非圧縮性ナビエ＝ストークス方程式の半ラグランジュ設定において、エネルギー保存を効果的に強制できるか？
• RQ4 自由度として「小さな辺（small edges）」を用いることで、スキームの安定性および精度にどのような影響を与えるか？
• RQ5 計算コストおよび保存性の観点から、標準的なエーラリアン手法と比べて本手法はどのように差をつけるか？

主な発見

• 線形的陰解法を用いた完全離散スキームは、粘性係数がゼロに近づく極限において優れた安定性を示し、非粘性非圧縮性エーラー流のロバストなシミュレーションを可能にする。
• 数値実験により、単体メッシュ上での空間的・時間的両方で2次収束が確認された。
• ラグランジュ乗数を用いてエネルギー保存が成功裏に強制され、離散エネルギーが連続エネルギー関係を適切に追跡し、有界性を保った。
• 境界上でのトレース作用素および外微分計算作用素を用いることで、非圧縮性および接線成分境界条件が維持される。
• 数値結果から、保存型バージョンのスキームは長時間にわたるシミュレーションでもエネルギーを保存するのに対し、非保存型バージョンではエネルギーのずれ（エネルギードリフト）が生じることが示された。
• 本手法は複雑な領域に対しても適用可能であり、粒子や辺が領域外に部分的に運搬されても精度を保ち、境界データの取り扱いが一貫的である。