许爱国：从物理学角度评述LBM方法.pdfVIP

物理学进展 Progress in Physics 2014年 （已录用稿件） 非平衡与多相复杂系统模拟研究 ——Lattice Boltzmann 动理学理论与应用 1* 1 2 1 许爱国 ，张广财 ，李英骏 ，李华 1 北京应用物理与计算数学研究所计算物理重点实验室，北京 100088 2 中国矿业大学（北京）深部岩土力学与地下工程国家重点实验室，北京 100083 * Corresponding author. Email: Xu_Aiguo@ 收稿日期：********* 摘 要：在自然界和工程物理领域存在大量的非平衡、多相等复杂系统和复杂行为。Lattice Boltzmann （LB ） 方法起源于复杂系统复杂行为研究的格子气或元胞自动机模型；其中，现代版的 Lattice Boltzmann Kinetic Model （LBKM ）植根于非平衡统计物理学的基本方程——Boltzmann 方程。本文从物理学视角评述 LB 方法， 给出单松弛因子和多松弛因子 LBKM 构建的统一理论，介绍其在非平衡与多相复杂系统研究方面的应用。简 单列举 LB 在多相流、可压流、材料动理学等方面的进展，重点介绍使用 LB 研究流体界面不稳定性、燃烧等 问题的工作。本文所重点传递的信息为：可以通过宏观量研究系统的非平衡行为、可以提供系统偏离热力学 平衡引发的宏观效应是 LBKM 建模优越于宏观连续介质建模的地方；除了可以从更基本的层面理解相应物理 系统的特征、机制和规律外，这类研究结果可以为现有程序或软件中宏观模型的改进（例如修正项的构造） 提供物理参考。 关键词：格子玻尔兹曼；动理学模型；非平衡效应；复杂系统 中图分类号：O41 1 1. 引言 （Richtmyer-Meshkov ，瑞奇迈尔-莫西科夫）不 稳定性的增长，而 RM 不稳定性的发展又会带来 除了日常生活中常见的各类复杂流动过程 其它不稳定性的增长： 在内爆加速阶段，轻材 外，工业领域存在大量的种类各异、起源不同的 料 包 层 加 速 重 材 料 （ 金 属 ） 会 引 发 RT 非平衡与多相复杂流体系统。首先，强冲击与爆 （Rayleigh-Taylor ，瑞利-泰勒）不稳定性；在反 轰产生的压强可以达到数十万个大气压，是多数 弹阶段，气芯的轻材料加速重材料壳层，同样引 固体材料屈服应力的几十倍。当如此强的冲击波 发 RT 不稳定性。在 RT 不稳定性发展的后期， 在固体材料（包括炸药）内进行传播时，波后物 轻重流体两侧的大变形切向流又会导致 KH 质表现出流体行为，相当一部分动力学、动理学 （Kelvin-Helmholtz ，开尔文-亥姆霍兹）不稳定 特征可采用流体模型来描述。由于实际材料内往 性的增长，出现卷曲、翻滚的现象，