Lattice Boltzmann framework for multiphase flows by Eulerian-Eulerian Navier-Stokes equations

本論文は、有限差分補正を一切必要とせず、非常に大きな密度比や現実的な抗力係数を扱える、6 つの結合された格子ボルツマン法（LBM）スキームからなる新しいマルチフェーズ流シミュレーションフレームワークを初めて提案し、その有効性を従来の数値解法との比較で実証したものである。

要約

1. 背景：なぜこの研究が必要なのか？

私たちが普段見ている「泡立つお風呂」や「工場の大きなタンクの中」では、空気（気体）と水（液体）が激しく混ざり合っています。これを「多相流（マルチフェーズフロー）」と呼びます。

この動きをコンピュータでシミュレーションするのは、従来の方法だと**「重くて、遅くて、計算が複雑」**でした。

• 従来の方法（有限差分法など）： 地図のマス目一つ一つで、隣り合うマスとの関係性を手計算のように丁寧に計算していく方法。正確ですが、計算量が膨大で、大規模なスーパーコンピュータを使っても時間がかかります。
• 新しい方法（格子ボルツマン法・LBM）： 流体を「小さな粒子」に見立て、マス目の上を「跳ね回る」ように計算する方法。これは**「並列処理（何人もの人が同時に作業する）」**が得意で、現代の超高性能コンピュータ（HPC）に非常に合っています。

しかし、これまでの LBM は、液体と気体が混ざった複雑な計算をするには「足し算や引き算の補正（有限差分補正）」が必要で、そのせいで「並列処理の利点」が半減してしまっていました。

2. この論文のすごいところ：6 つのチームで協力する

この論文の著者たちは、**「補正なしで、純粋な LBM だけで、液体と気体の複雑な動きを再現できる」**という新しい枠組み（フレームワーク）を提案しました。

これを理解するために、**「巨大な工場」**を想像してください。

• 従来の工場： 液体と気体の動きを計算するには、特殊な「補正係数」というマニュアルが必要で、作業員がそれを参照しながら計算していたため、スピードが出なかった。
• 新しい工場（この論文の手法）：
  1. 6 つの専門チームが同じ工場（同じ格子）で働きます。
  2. チーム A & B（気体と液体の動き）： 空気と水そのものがどう動くかを計算。
  3. チーム C & D（泡の量）： 「どこにどれくらい泡があるか（体積分率）」を計算。
  4. チーム E & F（源流）： 「なぜ泡が動いたのか？」という原因（圧力や摩擦）を計算。

これら 6 つのチームが**「同じリズムで、互いに協力しながら」**働きます。

• メリット： 誰かが「補正マニュアル」を参照する必要がありません。全員が自分の役割を素直にこなすだけで、全体として完璧な結果が出ます。
• 結果： スーパーコンピュータの力を最大限に引き出し、これまで難しかった「空気と水の激しい混ざり合い」や「密度が極端に違う（水と空気のような）状態」も、安定して計算できるようになりました。

3. 具体的な工夫：2 つの「魔法の道具」

この新しい工場がうまく回るために、著者たちは 2 つの工夫をしました。

① 「人工的な圧縮」のアイデア

本来、液体は「圧縮できない（硬い）」ものですが、計算をスムーズにするために、**「少しだけ柔らかく（圧縮可能に）見せる」**というトリックを使います。

• 例え話： 硬いゴムボールを、少しだけ風船のように見せて計算するイメージです。
• 工夫： 液体と気体で「柔らかさの度合い」を少し変えて調整し、最終的には「硬い液体」の動きとして正しい結果を導き出します。これにより、複雑な圧力計算を、LBM が得意な「粒子の動き」だけで解決しました。

② 「密度の差」が極端な場合の対策

空気と水のように、重さが 1000 倍も違う場合、計算が不安定になりがちです（小さな計算ミスが大きな暴走に繋がる）。

• 工夫： 計算の「振動」を抑えるために、**「少しだけ damping（制動）」**という仕組みを入れました。
• 例え話： 揺れやすい橋に、揺れを吸収するダンパー（緩衝材）を取り付けるようなものです。これにより、密度差が極端な場合でも、計算が暴走せず、安定してシミュレーションが進みます。

4. 結果：従来の方法と比べてどう？

著者たちは、この新しい方法をテストしました。

• テスト 1（基本）： 従来の方法（有限差分法）と見比べて、**「ほぼ同じ結果」**が出ました。
• テスト 2（高難易度）： 密度差が極端な場合や、現実的な「抵抗（ドラッグ）」の計算を含めた場合でも、**「安定して、かつ正確」**に計算できました。

特に重要なのは、**「補正（Finite Difference Correction）を一切使っていない」**ことです。これにより、LBM が本来持っている「並列処理の速さ」と「計算の単純さ」を損なわずに、複雑な多相流を扱えるようになりました。

5. まとめ：これがなぜ重要なのか？

この研究は、**「流体シミュレーションの未来」**を開くものです。

• エネルギー分野への貢献： 石油・ガスの精製、化学反応、原子炉の冷却など、液体と気体が絡み合う複雑なプロセスを、より安く、早く、正確に設計できるようになります。
• スーパーコンピュータの活用： 最新の並列コンピュータの性能を 100% 発揮できるため、より巨大で詳細なシミュレーションが可能になります。

一言で言えば：
「これまで、複雑な流体シミュレーションをするには『重たい荷物を背負って（補正を付けて）』歩かなければならなかった。でも、この新しい方法なら、**『荷物を下ろして、軽やかに、何人もの仲間と手分けして走れる』**ようになったのです。」

この「新しい歩き方（フレームワーク）」は、将来のエネルギー技術や環境問題の解決に大きく貢献することが期待されています。

技術概要

以下は、arXiv:2409.10399v4 に掲載された論文「Lattice Boltzmann framework for multiphase flows by Eulerian-Eulerian Navier-Stokes equations」の技術的な要約です。

1. 研究の背景と課題 (Problem)

多相流（気液二相流など）のシミュレーションは、石油・ガス産業や化学反応器（バブルカラムなど）の設計・最適化において極めて重要です。従来の数値流体力学（CFD）では、有限差分法（FD）や有限体積法（FVM）を用いてオイラー - オイラー（二流体）アプローチの支配方程式を解くことが一般的です。しかし、このアプローチには以下の課題があります。

• 数値的複雑性: 各相の運動量方程式における特異点（体積分率が 0 に近づく際の発散）、相間結合による数値的不安定性、体積分率の物理的範囲（0〜1）の保証、チェッカーボード不安定性の回避などが困難です。
• LBM の限界: ラティスボルツマン法（LBM）は並列計算に優れていますが、従来のオイラー - オイラー方程式を直接適用するには、有限差分による補正が必要になるなど、LBM の「局所性」という最大の利点を損なう修正が求められていました。特に、非常に大きな密度比（気体と液体の密度差が大きい場合）や現実的な抗力モデルを扱う際の安定性が課題となっていました。

2. 提案手法 (Methodology)

本研究では、有限差分補正を一切用いずに、オイラー - オイラー多相流方程式を解くための新規 LBM フレームワークを提案しています。この手法は、同じ格子上で実行される6 つの結合された LBM スキームで構成されています。

2.1. 人工圧縮性アプローチと 6 つのスキーム

非圧縮性の支配方程式を直接解くのではなく、各相に対して「人工圧縮性（Artificial Compressibility）」概念を導入し、擬似圧縮性の連続の式と運動量方程式を解くことで、共通の圧力場を導出します。

1. 気相・液相の運動量・連続の式用 (fg, fl):
   • 各相の速度場と圧力を計算します。
   • 気相と液相で異なる状態方程式（Equation of State）を適用し、同じ圧力勾配（空間的結合）と時間的な圧力変化（時間的結合）を両方の相に正しく伝播させるよう設計されています。
   • 非圧縮性平衡分布関数を用いて、密度比が大きい場合でも安定性を保つように調整されています。
2. 体積分率の輸送用 (fαg, fαl):
   • 各相の体積分率（$\alpha_g, \alpha_l$）の輸送方程式を解きます。
   • 体積分率が 0〜1 の範囲に収まるよう、Spalding による正規化手法を LBM の衝突・ストリーミングサイクル内で適用しています。
3. 相間連続性のソース項計算用 (fβg, fβl):
   • 人工圧縮性項（$\nabla \cdot u \neq 0$ となる補正項）を計算するために使用されます。
   • これにより、有限差分法を用いずに、体積分率の勾配や応力テンソルなどの微分項を LBM 内部で直接計算可能にしています。

2.2. 数値的安定化と拡張性

• 高密度比への対応: 密度比が極めて大きい場合（例：気体/液体 $\approx 1000$）、運動量方程式の液体項における $1/R$ の項が数値誤差と同程度になり不安定化します。これを防ぐため、連続の式に散逸項（ダッシュポット）を追加し、ソース項の更新頻度を調整する手法を導入しています。
• 現実的な抗力モデル: Clift, Grace, Weber (CGW) モデルなどの、体積分率に依存する複雑な抗力係数を扱えるようにしています。
• 次元独立性: 提案された手法と数式は、物理次元（1 次元、2 次元、3 次元）に依存せず、任意の次元に適用可能です。

3. 主要な貢献 (Key Contributions)

• 初の完全 LBM フレームワーク: 有限差分補正を一切使わず、オイラー - オイラー多相流方程式を解く LBM フレームワークを初めて提案しました。
• 6 つの結合スキーム: 運動量、連続、体積分率、ソース項をそれぞれ独立かつ結合された分布関数で処理する効率的な構造を確立しました。
• 高密度比と現実モデルの統合: 非常に大きな密度比（$R \approx 833$）や、体積分率に依存する剛直な（stiff）抗力モデルを安定してシミュレーションできることを実証しました。
• HPC への適合性: 局所的な計算のみで完結するため、大規模な高性能計算（HPC）施設や新規の並列ハードウェアでのスケーラビリティが極めて高いです。

4. 数値検証と結果 (Results)

1 次元の垂直管（気液二相流）をテストケースとして、提案 LBM 手法と従来の有限差分法（FD）ソルバーによる基準解を比較しました。

• 基本テスト (R=2, R=5): 気相・液相の速度、体積分率、動圧力分布において、LBM と FD の結果は極めて良好な一致を示しました。特に、混合流速が発散自由（$\nabla \cdot (\alpha_g u_g + \alpha_l u_l) = 0$）である一方、個々の相の流速は発散しないという物理的挙動を LBM が正しく再現していることが確認されました。
• 高密度比テスト (R=833.3): 密度比が非常に大きい場合でも、提案された安定化手法（散逸項とソース項の更新制御）により、数値的発散なく安定した解が得られました。
• 現実的抗力モデルテスト (CGW モデル): 体積分率に依存する複雑な抗力モデルを用いた場合でも、定常状態において FD 解と良好な一致を示しました（局所的な不一致は存在するものの、物理的な傾向は正確に捉えられています）。
• 微分量の精度: 有限差分演算子を使わずに、LBM 分布関数から直接計算された体積分率の勾配や応力テンソルが、FD による計算値と高い精度で一致していることが確認されました。

5. 意義と将来展望 (Significance)

• 計算効率の飛躍的向上: 有限差分補正や非局所的な計算を排除したため、LBM の持つ並列計算の利点を最大限に引き出せます。これにより、大規模な多相流シミュレーションが HPC 上で現実的な時間で実行可能になります。
• 産業応用への道筋: バブルカラム反応器や核融合炉の冷却材循環など、気液相互作用が複雑なエネルギー分野の装置設計において、高精度かつ高効率なシミュレーションツールを提供します。
• 汎用性: 手法が次元に依存しないため、今後の 2 次元・3 次元への拡張が容易であり、複雑な幾何形状を持つ実機シミュレーションへの応用が期待されます。

結論として、本研究は多相流シミュレーションにおける LBM の適用範囲を大きく広げ、従来の CFD ソルバーに匹敵する精度を保ちつつ、HPC 環境での大規模計算を可能にする画期的なフレームワークを提示したものです。