[论文解读] Entropic Lattice Boltzmann Method for Large Scale Turbulence Simulation

I. V. Karlin, Santosh Ansumali+3 more|arXiv (Cornell University)|May 30, 2003

Lattice Boltzmann Simulation Studies参考文献 16被引用 7

一句话总结

本文提出熵晶格玻尔兹曼方法（ELBM）作为一种稳定、基于动力学的大型湍流模拟方法，通过强制满足H定理来确保分布函数的正性并抑制数值不稳定性。结果表明，ELBM能准确再现巴托勒的 $k^{-3}$ 能量谱，并表现出间歇性的负有效粘性，表明其具有内在的亚网格尺度建模能力，可捕捉传统涡粘性模型无法描述的反向散射和局部不稳定性。

ABSTRACT

Recently, a minimal kinetic model for fluid flow, known as entropic lattice Boltzmann method, has been proposed for the simulation of isothermal hydrodynamic flows. At variance with previous Lattice Boltzmann methods, the entropic version permits to describe the non-linear dynamics of short scales in a controlled and stable way. In this paper, we provide the first numerical evidence that the entropic lattice Boltzmann scheme provides a quantitatively correct description of the large-scale structures of two-dimensional turbulence, free of numerical instabilities. This indicates that the entropic lattice Boltzmann method might provide a starting basis for the formulation of a new class of turbulence models based on genuinely kinetic principles.

研究动机与目标

• 开发一种用于模拟具有未解析亚网格尺度动力学的大规模湍流的数值稳定晶格玻尔兹曼方法。
• 克服标准晶格玻尔兹曼格式在模拟接近网格尺度湍流时固有的数值不稳定性。
• 检验熵晶格玻尔兹曼方法是否能再现二维湍流的关键统计特征，如巴托勒谱。
• 探究ELBM是否能捕捉非平衡效应（如负有效粘性），这些效应在传统涡粘性模型中缺失。
• 建立ELBM作为一类基于动力学理论和统计力学原理的新湍流模型的基础。

提出的方法

• 熵晶格玻尔兹曼方法通过约束分布函数满足离散熵最大化条件来强制实现H定理，确保其正性与数值稳定性。
• 该方法使用一个可动态调整的弛豫参数 $\alpha\beta$，以维持熵产生，防止近网格尺度模态中的非物理振荡。
• 有效粘性定义为 $\nu_{\rm eff} = c_s^2 \left( \frac{1}{\alpha\beta} - \frac{1}{2} \right)$，从而可直接监测其与分子粘性的偏差。
• 模拟在二维周期性区域内进行，采用粗网格（$\delta \gg l_k$），模拟未解析的湍流以检验亚网格尺度建模能力。
• 在多个时间点计算能量谱和涡量谱，以验证巴托勒谱的 $k^{-3}$ 标度。
• 分析 $\delta\nu/\nu = \nu_{\rm eff}/\nu - 1$ 的概率分布，以检测正、正常和负有效粘性的间歇性事件。

实验结果

研究问题

• RQ1 熵晶格玻尔兹曼方法是否能在未解析模拟中再现二维湍流特征性的 $k^{-3}$ 能量谱？
• RQ2 ELBM 是否表现出间歇性的负有效粘性，表明存在反向散射和标准涡粘性模型无法捕捉的局部不稳定性？
• RQ3 ELBM 内置的亚网格模型在多大程度上源于统计力学原理，而非经验调参？
• RQ4 弛豫参数 $\alpha\beta$ 的分布如何反映湍流中空间间歇性和非平衡动力学？
• RQ5 在无显式涡粘性模型或高阶离散化的情况下，ELBM 是否能保持数值稳定性？

主要发现

• ELBM 模拟在 $t = 15000$ 时再现了巴托勒 $k^{-3}$ 能量谱，最佳拟合常数为 $2 \times 10^{-7}$，转移速率为 $1.7 \times 10^{-11}$，推导出的常数为 4.8，接近理论值 2.626。
• 涡量谱在惯性区间内遵循预期的 $k^{-1}$ 标度，证实了涡量传递的正确动力学。
• 观察到 $\delta\nu/\nu$ 的三模态分布：正常事件（$M_{\rm Normal} \approx 0.84$）、正事件（$M_{\rm Plus} \approx 0.1$）和负事件（$M_{\rm Minus} \approx 0.06$），表明存在间歇性负粘性。
• 负事件（负有效粘性）的比例随时间增加，从 $t=1000$ 时的 1.8% 上升至 $t=10000$ 时的 6.25%，表明局部不稳定性逐渐增强。
• 正常事件的分布与 $\delta\nu/\nu = 0$ 处的狄拉克δ函数一致，表明 $\alpha\beta$ 在大部分时间处于静止状态，仅在罕见、局部的活动爆发中发生扰动。
• 该方法在未解析模拟中成功避免了数值不稳定性，表现出鲁棒性，且无需依赖人工耗散或高阶正则化。