A Finite Element Method for Fluctuating Navier--Stokes Equations

技術要約：変動性ナビエ・ストークス方程式のための有限要素法

問題提起
熱ゆらぎは、ブラウン運動からキャビテーション、沸騰、流体不安定性に至るまで、物理現象の重要な駆動源である。マクロなスケールでは無視できることが多いが、近年の研究では、熱ノイズと乱流の結合、特にコルモゴロフ長スケール付近における結合が、システムの挙動を著しく変化させ得ることが示されている。これらの効果をモデル化するために、古典的なナビエ・ストークス方程式に確率的応力テンソルを加え、変動流体力学の枠組みを構築する。この枠組みは揺動散逸定理（FDT）に依拠しており、システムが正しい統計的平衡に達することを保証するために、確率的強制項の共分散が粘性散逸演算子に比例しなければならないことを規定している。

計算流体力学における主要な課題は、これらの方程式を離散化する際に、離散レベルでFDTを維持することである。有限差分法や有限体積法などの標準的な数値手法は、複雑な形状の扱いに苦慮することが多い。一方、有限要素法（FEM）は複雑な領域を扱う上で堅牢な手法を提供するが、標準的なガラーキン離散化は、ゆらぎ場に非物理的な長距離相関を導入してしまう。このアーティファクトは、基底関数のグローバルなサポートと、それに伴う非対角な質量行列に起因し、局所的な平衡統計を歪めてしまう。

手法
著者らは、等温・圧縮性変動ナビエ・ストークス方程式のための、新しい構造保存型有限要素フレームワークを提案している。本手法は、以下の3つの核心的な柱に基づいている。

1. 変分定式化と確率的強制： 本手法は、弱定式化の中で直接、確率的強制項を定義する。離散的な確率的強制の共分散は、離散的な粘性散逸演算子（$a_h(\cdot, \cdot)$）に明示的に比例するように構成される。これにより、離散系が連続体のFDTを模倣することを保証し、非物理的なエネルギー注入や人工的な減衰を防ぐ。
2. 節点求積法（マス・ランピング）： 標準的なFEMに特有の非物理的なメッシュスケールの相関を排除するため、著者らは主要な内積に対して節点求積則（台形則）を採用している。この技術により、対角成分を持つ質量行列（ランピングされた質量行列）が得られる。その結果、離散的な平衡共分散は節点において対角となり、ゆらぎが局所的かつ物理的に現実的なものとなる。これは、有限体積法の統計的特性を模倣しつつ、FEMの変分的利点を保持するものである。
3. 時間積分： 時間離散化には、セミインプリシット（半陰解法）$\theta$-スキームを用いる。著者らは特に、クランク・ニコルソン・スキーム（$\theta = 0.5$）を分析し、検証している。線形化された方程式の数学的解析により、クランク・ニコルソン・スキームは、離散共分散バランス方程式における$\Delta t$依存の誤差項を消失させる唯一の手法であることが示されている。これにより、時間ステップを空間離散化に結合させるような厳しい安定条件を課すことなく、定常的な揺動散逸バランスを満たすことが可能となる。

主な貢献

• 構造保存型離散化： 確率的強制の共分散を離散粘性演算子と一致させることで、離散レベルで揺動散逸バランスを厳密に保持する有限要素法を導入した。
• 人工的相関の排除： 節点求積法を用いることで、標準的な連続ガラーキンFEMアプローチを悩ませる非物理的な平衡共分散行列の非対角相関を排除することに成功した。
• クランク・ニコルソンの検証： 線形化領域において、離散FDTバランスを維持するための最適な選択肢はクランク・ニコルソン・スキームであることを、理論的証明を通じて示した。これは、時間ステップ依存の共分散バランスの偏差を排除するものである。
• 包括的な検証： 三角形、四角形、および四面体メッシュを用いた、1次元、2次元、および3次元空間における数値検証を実施した。

結果
数値実験により、提案手法が様々な離散化パラメータにわたって平衡ゆらぎ統計を正しく再現することが確認された。

• 分散のスケーリング： 密度および速度場の経験的分散は、統計力学によって予測される理論的な平衡値に収束する。結果は、正しいスケーリング挙動（$\langle \delta (\cdot)^2 \rangle \propto |T^i_h|^{-1}$）を示しており、メッシュサイズ（または要素体積）を半分にすると、分散がおよそ2倍になる。
• 堅牢性： 本手法は、異なるメッシュ解像度、時間ステップ、および要素タイプ（2Dにおける三角形 vs 四角形）に対して一貫した性能を示す。
• 質量保存： 連続ガラーキン離散化を用いているにもかかわらず、本スキームは全質量を保存する。これは、揺動散逸バランスを変化させる可能性のあるアップウィンド法やその他の安定化手法によって損なわれやすい重要な特性である。
• 限界： 著者らは、平均密度が非常に小さい場合（例：3Dにおける $\rho \approx 10^{-2}$）、線形化平衡公式の妥当性が低下すること、および本手法が正の密度（$\rho \geq 0$）という物理的制約を明示的に強制しないことを指摘している。しかし、中程度のゆらぎ振幅に対しては、本手法は正の密度を保持する有限差分法に匹敵する結果を再現している。

意義と主張
本論文は、従来の有限差分法や有限体積法が困難に直面する、複雑な幾何学的形状における圧縮性流体の熱ゆらぎをシミュレートするための、堅牢で正確なツールを提供すると主張している。離散スキームが基本的な揺動散逸バランスを尊重し、人工的な相関を排除することを保証することで、本手法は、標準的なFEMの確率的流体力学への適用において一般的なアーティファクトなしに、平衡統計を捉える信頼できる方法を提供する。

著者らは、本研究を、完全に結合された非線形ランダウ・リフシッツ・ナビエ・ストークス方程式へと有限要素法を拡張するための重要な一歩として位置づけている。現在の研究は等温・単成分流体に焦点を当てているが、変分フレームワークは多成分・非等温系、さらには適応メッシュ細分化や複雑な境界条件を伴うアプリケーションへと拡張可能な柔軟性を備えている。結論として、本手法は、将来的な課題として正の密度制約の強制への対処が必要であるものの、現実的な設定における変動流体力学のマルチスケールシミュレーションのための強力なツールであると述べている。