2025-11-14T22:25:11.017657

Lattice Boltzmann model for non-ideal compressible fluid dynamics

Hosseini, Feinberg, Karlin
We present a lattice Boltzmann formulation for the simulation of compressible, non-ideal fluid flows. The method employs first-neighbor lattices and introduces a consistent set of correction terms through quasi-equilibrium attractors, ensuring positive-definite and Galilean-invariant Navier-Stokes dissipation rates. This construction circumvents the need for extended stencils or ad hoc regularization, while maintaining numerical stability and thermodynamic consistency across a broad range of flow regimes. The resulting model accurately reproduces both Euler- and Navier-Stokes-level hydrodynamics. As a stringent validation, we demonstrate, for the first time within a lattice Boltzmann framework, quantitatively accurate simulations of drop-shock interactions at Mach numbers up to 1.47. The proposed approach thus extends the applicability of lattice Boltzmann methods to high-speed, non-ideal compressible flows with a minimal kinetic stencil.
academic

非理想圧縮性流体力学のための格子ボルツマンモデル

基本情報

  • 論文ID: 2510.14712
  • 題目: Lattice Boltzmann model for non-ideal compressible fluid dynamics
  • 著者: S. A. Hosseini, M. Feinberg, I. V. Karlin (ETH Zurich)
  • 分類: physics.flu-dyn
  • 掲載誌: Journal of Fluid Mechanics
  • 論文リンク: https://arxiv.org/abs/2510.14712

要旨

本論文は、圧縮性非理想流体流動をシミュレートするための格子ボルツマン法を提案している。本手法は第一近傍格子を採用し、準平衡吸引子を通じて一貫した修正項集合を導入することで、正定値性とガリレイ不変なNavierStokes散逸率を確保する。この構成は拡張テンプレートまたは一時的な正則化の必要性を回避しながら、広範な流動状態において数値安定性と熱力学的一貫性を維持する。得られたモデルはオイラーおよびNavierStokes階層の流体力学を正確に再現する。厳密な検証として、著者らは格子ボルツマンフレームワーク内でマッハ数1.47までの液滴-衝撃波相互作用の定量的に正確なシミュレーションを初めて実証した。

研究背景と動機

問題定義

非理想圧縮性流体力学は、近臨界、跨臨界および超臨界状態で動作する方法および技術の出現により、急速に発展している流体力学の新興分野である。このような流体には以下が含まれる：

  1. 超臨界流体：臨界点以上の圧力および温度条件下の流体
  2. 稠密蒸気：圧縮係数が理想気体の挙動から逸脱する気体
  3. BZT流体：負の基本気体動力学導関数を有する高分子複雑性流体

重要性

  1. 産業応用の需要：有機ランキンサイクル、超臨界CO₂タービンなどのエネルギー技術
  2. 実験データの不足：実験設定が複雑で、データ取得が困難
  3. 数値ツールの需要：複雑な物理現象を理解するための一貫性のある効率的な数値ツールが必要

既存手法の限界

  1. 従来のCFD手法：熱力学的に不安定なスピノーダル領域内では、断熱音速の二乗が負となり、発展方程式が双曲型から放物型に変わる
  2. 既存のLBMモデル：主に沸騰応用を対象としており、圧縮性非理想流動に対する体系的なモデルが不足している
  3. 数値安定性の問題：高マッハ数および強い非理想効果下での安定性の課題

核心的貢献

  1. 新しい二分布関数LBMモデルの提案：古典的第一近傍格子（D3Q27）を使用して非理想圧縮性流動シミュレーションを実現
  2. 準平衡吸引子機構の導入：fᵢおよびgᵢ修正項を通じて正定値性とガリレイ不変な散逸率を確保
  3. 独立体粘度制御の実現：BGK構造における体粘度が負になる可能性の問題を解決
  4. 高マッハ数液滴-衝撃波相互作用シミュレーションの初実現：LBMフレームワーク内でMa=1.47の定量的に正確なシミュレーションを達成
  5. 熱力学的一貫性の維持：広範な流動状態において数値安定性を保持

方法の詳細

タスク定義

Navier-Stokes-Korteweg (NSK)方程式系を満たす圧縮性非理想流体のシミュレーション：

  • 質量保存：∂ₜρ + ∇·(ρu) = 0
  • 運動量保存：∂ₜ(ρu) + ∇·(ρu⊗u) + ∇P + ∇·T^NS + ∇·T^K = 0
  • エネルギー保存：∂ₜ(ρE) + ∇·((ρE+P)u) + u·(∇·T^K) + ∇·(u·T^NS) + ∇·q = 0

ここでT^Kはkorteweg表面張力テンソル、T^NSは粘性応力テンソルである。

モデルアーキテクチャ

1. 二分布関数設計

二組の分布関数{fᵢ, gᵢ}を使用する（i=1,...,27）：

  • fᵢ：質量および運動量保存を処理
  • gᵢ：エネルギー保存を処理

2. 格子ボルツマン方程式

fᵢ(x+cᵢδt, t+δt) = fᵢ + 2β(fᵢᵉᵠ - fᵢ) + (1-β)(fᵢ* - fᵢᵉᵠ)
gᵢ(x+cᵢδt, t+δt) = gᵢ + 2β(gᵢᵉᵠ - gᵢ) + (1-β)(gᵢ* - gᵢᵉᵠ)

ここでβ = Pδt/(2μ + Pδt)は緩和パラメータである。

3. 平衡分布関数

積の形式を採用する：

fᵢᵉᵠ = ρ ∏_α Ψᵢα(uα, P/ρ + uα²)

ここでΨᵢαは基本関数である：

Ψᵢα(ξα, ζαα) = 1 - cᵢα² + ½[(3cᵢα² - 2)ζαα + cᵢαξα]

4. 準平衡修正項

f分布関数の修正：

  • 力項：Fα = -κρ∂α∇²ρ （Korteweg力）
  • 応力修正：Φααは対角応力テンソル修正と独立体粘度制御を含む

g分布関数の修正：

  • シフト速度：u* = u + δt F/ρ
  • シフト温度：T* = T - δt²F·F/(2ρ²cᵥ)
  • 熱流修正：非理想気体のエンタルピ勾配効果を考慮

技術的革新点

1. 準平衡吸引子機構

fᵢおよびgᵢを導入することで、システムは単純に局所平衡に緩和するのではなく、修正された準平衡状態に緩和し、以下を確保する：

  • 正確なマクロ動力学
  • 数値安定性
  • ガリレイ不変性

2. 独立体粘度制御

従来のBGKモデルでは体粘度η' = ((D+2)/D - ρcs²/P)μが負になる可能性がある。修正項Φ'を通じて独立制御を実現する：

Φ' = P((D+2)/D - ρcs²/P - η/μ)(∇·u)

3. 熱力学的一貫性

注意深く設計されたシフト温度T*およびエネルギー修正を通じて、非線形状態方程式下での熱力学的一貫性を確保する。

実験設定

検証ケース

  1. 基本検証：
    • 音速測定
    • 剪断粘度検証
    • 熱Couette流
    • 法線モード減衰
  2. 多相流検証：
    • 液-気共存密度
    • インターフェース一貫性および収束性
  3. 圧縮性流検証：
    • 一次元非理想衝撃波管
    • 衝撃波-液柱相互作用

物質パラメータ

主に窒素(N₂)を用いた検証：

  • 臨界圧力：Pc = 3.4×10⁶ Pa
  • 臨界密度：ρc = 241.96 kg/m³
  • 臨界温度：Tc = 126.2 K
  • R/cv = 0.4

状態方程式

van der Waals状態方程式を採用：

P(ρ,T) = ρRT/(1-bρ) - aρ²

実験結果

主要結果

1. 基本物理量の検証

  • 音速：飽和液体および蒸気分枝上で、Tr∈0.7,1範囲内で解析解と完全に一致
  • 剪断粘度：異なるマッハ数下でガリレイ不変性を維持し、誤差<1%
  • 熱Couette流：Pr∈{0.6,1.2,4.9}およびMa∈{0.8,1.2,1.6}下で解析解と一致

2. 多相流特性

  • 共存密度：Tr∈0.3,0.99範囲内でMaxwell等面積構成と完全に一致
  • インターフェース収束性：異なるメッシュ分解能下で良好な収束特性を示す

3. 非理想衝撃波管

三種類の異なる非理想衝撃波管構成のシミュレーションに成功：

  • ケースI：Γ<0領域の希薄波
  • ケースII：古典的圧縮波挙動
  • ケースIII：Γ=0線を動的に横切る混合希薄波

すべてのケースで参照データとの一致度は優秀である。

4. 衝撃波-液滴相互作用

  • LBMフレームワーク内でMa=1.47の高マッハ数シミュレーションを初実現
  • 入射衝撃波、透過波、反射衝撃波、Mach棒などの複雑な波構造を正確に捕捉
  • 液滴変形進化は実験データと高度に一致

アブレーション実験

異なるメッシュ分解能のインターフェース試験を通じてモデルの収束性を検証：

  • δxが5μmから0.1μmの四つの分解能レベル
  • 良好なメッシュ収束特性を示す

主要な知見

  1. 数値安定性：高マッハ数および強い非理想効果下で安定性を維持
  2. 計算効率：標準第一近傍格子を使用し、計算コストは制御可能
  3. 物理的忠実度：非理想流体の複雑な物理現象を正確に再現

関連研究

従来のLBM多相流モデル

  1. 色勾配モデル（Gunstensen et al. 1991）
  2. 疑似ポテンシャルモデル（Shan & Chen 1993）
  3. 自由エネルギーモデル（Swift et al. 1996）

圧縮性LBM発展

  • He & Doolen (2002)：圧縮性非理想流動への初の試み
  • 既存モデルは主に沸騰応用を対象
  • 圧縮性非理想流動に対する体系的なモデルが不足

本論文の相対的優位性

  1. 体系的手法：完全な二分布関数フレームワーク
  2. 理論的基礎：厳密なマルチスケール解析
  3. 検証の完全性：基本物理量から複雑流動現象までの包括的検証

結論と考察

主要な結論

  1. 圧縮性非理想流体に適用可能なLBMモデルの開発に成功
  2. 準平衡吸引子機構を通じて数値安定性と物理的一貫性を実現
  3. LBMフレームワーク内で高マッハ数液滴-衝撃波相互作用シミュレーションを初実現
  4. モデルは広範な流動状態において熱力学的一貫性を維持

限界

  1. 単一緩和時間モデル：現在は単純なBGK衝突モデルを使用しており、乱流シミュレーションには十分でない可能性がある
  2. 状態方程式の制限：主にvan der Waals方程式を検証しており、他の複雑な状態方程式はさらなる検証が必要
  3. 三次元検証の不足：ほとんどの検証は一次元および二次元構成に集中

今後の方向性

  1. より堅牢な衝突モデルへの拡張：乱流構成のシミュレーション用
  2. より複雑な状態方程式：Peng-Robinsonなど
  3. 大規模三次元応用：実際の工学問題への応用

深い評価

利点

  1. 理論的厳密性：完全なマルチスケール解析と数学的導出
  2. 方法の革新性：準平衡吸引子機構は重要な革新
  3. 検証の包括性：基本物理から複雑応用までの体系的検証
  4. 実用的価値：重要な工学応用問題を解決
  5. 数値効率：標準格子を使用して計算複雑性を回避

不足

  1. 衝突モデルの単純性：BGKモデルの複雑流動への適用性は限定的
  2. パラメータ依存性：特定の修正パラメータの選択に体系的指導が不足
  3. 三次元検証の不足：より多くの三次元複雑幾何の検証が必要

影響力

  1. 学術的貢献：LBMを圧縮性非理想流動分野に新たな方向性を開く
  2. 工学的価値：超臨界流体技術に重要な数値ツールを提供
  3. 方法論的意義：準平衡吸引子の概念は他の数値手法に影響を与える可能性がある

適用シナリオ

  1. 超臨界流体技術：CO₂タービン、有機ランキンサイクル
  2. 多相流動：液滴破裂、相変化過程
  3. 衝撃波動力学：高速流における複雑な波相互作用
  4. 直接数値シミュレーション：高精度が必要な基礎物理研究

参考文献

  1. Guardone, A. et al. (2024). Nonideal compressible fluid dynamics of dense vapors and supercritical fluids. Annual Review of Fluid Mechanics, 56, 241-269.
  2. Hosseini, S.A. & Karlin, I.V. (2023). Lattice Boltzmann for non-ideal fluids: Fundamentals and practice. Physics Reports, 1030, 1-137.
  3. Shan, X. & Chen, H. (1993). Lattice Boltzmann model for simulating flows with multiple phases and components. Physical Review E, 47(3), 1815.