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

Lattice Boltzmann model for non-ideal compressible fluid dynamics

基本信息

  • 论文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

摘要

本文提出了一种用于模拟可压缩非理想流体流动的格子玻尔兹曼方法。该方法采用第一近邻格子，通过准平衡吸引子引入一致的修正项集合，确保正定和伽利略不变的Navier-Stokes耗散率。这种构造避免了对扩展模板或临时正则化的需求，同时在广泛的流动状态范围内保持数值稳定性和热力学一致性。所得模型准确再现了欧拉和Navier-Stokes层次的流体力学。作为严格验证，作者首次在格子玻尔兹曼框架内展示了马赫数高达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. 非理想激波管

成功模拟三种不同的非理想激波管配置：

  • Case I：Γ<0区域的稀疏激波
  • Case II：经典压缩波行为
  • Case III：动态跨越Γ=0线的混合稀疏波

所有案例与参考数据的吻合度优秀。

4. 激波-液滴相互作用

  • 首次在LBM框架内实现Ma=1.47的高马赫数模拟
  • 准确捕获了入射激波、透射波、反射激波、马赫杆等复杂波结构
  • 液滴变形演化与实验数据高度一致

消融实验

通过不同网格分辨率的界面测试验证了模型的收敛性：

  • δ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.