[论文解读] Comment on the paper Li-Shi Luo, Wei Liao, Xingwang Chen, Yan Peng and Wei Zhang, Numerics of the lattice Boltzmann method: Effects of collision models on the lattice Boltzmann simulations, Physical Review E 83, 056710 (2011)

I. V. Karlin, Sauro Succi|arXiv (Cornell University)|Jul 15, 2011

Lattice Boltzmann Simulation StudiesEngineering被引用 16

一句话总结

此评论反驳了一篇2011年发表于《Physical Review E》的论文，该论文声称熵格式格子Boltzmann方程（ELBE）劣于标准LBGK模型。作者指出，文中标记为ELBE的模拟实际上是一种具有恒定松弛时间的标准LBGK格式，而非使用基于熵的自适应松弛参数的真实ELBE。因此，该比较无效且存在循环论证，动摇了该论文关于ELBE稳定性和性能的结论。

ABSTRACT

Critical comments on the entropic lattice Boltzmann equation (ELBE), by Li-Shi Luo, Wei Liao, Xingwang Chen, Yan Peng and Wei Zhang, Physical Review E 83, 056710 (2011), are based on simulations which make use of a model that, despite being called ELBE by the authors, is in fact fully equivalent to the standard lattice Bhatnagar-Gross-Krook equation. As a result, the conclusion of Luo et al on ELBE is circular, hence devoid of scientific bearing.

研究动机与目标

• 纠正罗等人（2011）论文中的科学误述，该文声称ELBE在稳定性和性能上劣于LBGK。
• 阐明ELBE方案与LBGK的根本区别在于其基于熵的自适应松弛参数，该参数受H定理约束。
• 证明罗等人并未模拟真正的ELBE，而是使用了恒定的松弛时间，导致其比较无效。
• 强调ELBE的真正价值在于其通过在未充分解析的模拟中引入可变粘度实现亚网格稳定化，而非直接数值模拟。
• 反驳认为ELBE中可变粘度不物理的观点，指出此类有效粘度在计算流体动力学中极为常见，且对高雷诺数流动至关重要。

提出的方法

• 指出罗等人在其模拟中使用了恒定的松弛时间τ，这违反了ELBE的核心机制，即通过熵条件动态确定松弛参数α。
• 重申正确的ELBE公式：分布函数演化满足 f_i(x+v_i,t+1) - f_i(x,t) = αβ(f^eq_i - f_i)，其中α通过求解熵条件 H(αf^eq + (1-α)f) = H(f) 确定，以确保H定理成立。
• 强调ELBE与LBGK的唯一区别在于平衡分布函数和松弛参数，其中LBGK的松弛参数固定，而ELBE的松弛参数自适应。
• 证明当分布函数接近平衡（如在解析良好的模拟中），ELBE退化为α = 2的LBGK，因此在直接数值模拟（DNS）范围内与LBGK完全相同。
• 指出罗等人采用的熵平衡分布的多项式近似无法恢复真正的ELBE行为，因为其未强制执行熵条件（2）。
• 引用凯蒂ng等人、斯帕索夫等人、格里丁克与霍伊斯特拉、米田与佐藤的独立ELBE实现，均证实其正确应用了熵条件和自适应松弛机制。

实验结果

研究问题

• RQ1 具有恒定松弛时间的熵格式格子Boltzmann方程（ELBE）是否与标准LBGK模型具有相同的数值稳定性？
• RQ2 当松弛参数固定时，ELBE方案与LBGK模型是否在本质上不同？
• RQ3 若未强制执行确保H定理成立的熵条件（2），能否将标记为ELBE的模拟视为科学上有效的？
• RQ4 为何声称ELBE是LB模型中‘最差’的论断在科学上站不住脚？
• RQ5 ELBE中的可变粘度是否不物理？其是否以牺牲计算效率的方式使该方法失效？

主要发现

• 罗等人（2011）标记为ELBE的模拟并非ELBE，而是一种具有恒定松弛时间τ的标准LBGK方案，导致其比较无效。
• 罗等人在自身论文中承认未测试ELBE的可变松弛时间，证实其实施并非真正的ELBE。
• 熵条件（2）对ELBE的H定理和亚网格稳定化至关重要；若缺失，该模拟不能被视为ELBE。
• 熵平衡与其实用多项式近似之间的差异为O(u^4)量级，处于低马赫数流动下格子Boltzmann方法的截断误差范围内。
• ELBE的真正优势在于其在未充分解析模拟中的自适应稳定能力，仅在必要时和必要位置部署可变粘度，而LBGK在类似条件下会失效。
• 求解非线性熵条件的计算开销，被高雷诺数下所需减少的网格分辨率所抵消，使ELBE在实际应用中总体更高效。