格子玻尔兹曼(LBM)小白的进阶之路

格子玻尔兹曼(LBM)小白的进阶之路

起因

2020年4月16日，距离新型冠状病毒爆发已经五个月的时间了，作为一名科研人员，高校迟迟未发开学通知，作为一名实验人员，第一次那么迫切希望可以熬夜做实验。但无论如何都回不到校园里，做实验，搞科研，写文章了。今天钉钉响了起来，导师发来一段简单的文字：最近的那篇文章好好写写，加点模拟的东西，发个高影响因子的期刊。再加上自己毕业需求，所以开始着手学习LBM模拟。B站一直是我的学习之源，这次很失望，没找到相关资源，网店上也没有找到相关书籍。于是抱着侥幸的心理百度了一下，然后就发现了CSDN这个社区。迫于现实开了一个月VIP， 下载了几篇相关的LBM书籍，目前主要还是先看何雅玲老师他们出版的《格子Boltzmann方法的理论及应用》。然后顺便在此做笔记，记录心得。书籍是在2008年出版的。鄙人能源相关专业，多相流为主要研究方向。
LBM，原来是Lattice Boltzmann Method的缩写，诞生到现在（2020）已经30多年了，算是新发展起来的模拟方法了。为进一步了解流动与传热机理，宏观尺度的实验研究和理论分析远远不够，需要从微观和跨尺度的介观角度进行深刻分析。因此而诞生LBM，该方法基于分子动理论，宏观上是离散方法，微观上是连续方法，因此又被称为介观模拟方法。具有天生的并行特性，边界条件处理简单，程序容易实施等优点，但理论体系尚不完善，国内应用还处于探索阶段（感觉现在很成熟了）。
何雅玲老师他们的书主要分为8个章节，所以我也基本上分成8部分进行笔记记录，心得随写。

第一章绪论

感觉自己基础知识不大牢固，不过感觉绪论应该不需要这些知识，所以就没有先去看附录里面的一些基本概念。先尝试往下走。
LBM是计算流体力学和计算传热学的一个新分支，流体力学最早可追溯到阿基米德时代，而传热学则产生于第一次工业革命时期。三大科学研究手段为实验研究，理论分析，数值模拟，相辅相成，缺一不可。计算流体和传热方法分为宏观方法、介观方法和微观方法。宏观层次上，流体被假设为连续介质，运动满足质量守恒、能量守恒和动量守恒，计算上通过离散将方程组离散成各种代数方程组，如有限差分法、有限容积法、有限元法、有限分析法、边界元法、谱方法，对应的商业软件有PHOENICS, FLUENT, STAR-CD, DFX。介观或微观层次，流体不再被假设为连续介质：微观层面上流体由大量的离散分子组成，如分子动力学模拟；介观层面上流体被离散成一系列的流体粒子（微团），常见的模拟方法有格子气自动机、LBM以及直接模拟蒙特卡罗方法。LBM 方法中除了流体被离散成流体粒子外，物理区域也被离散成一系列的格子，时间被离散成一系列的步长。描述流体粒子运动的方程被称为Boltzmann方程或相应的离散形式。第一篇关于LBM的文章最早出现在1988年。大致的改善过程如下：最早的格子气自动机—采用局部粒子分布函数和Boltzmann方程，形成多粒子碰撞模型（最初的LBM模型）—将碰撞算子线性化，形成矩阵模型—利用单松弛时间法简化平衡态，形成格子BGK(LBGK)模型—为消除降低压缩效应，提出多种等温不可压模型、不可压缩热模型（双分布函数模型），逐渐成熟。专业的商业软件为PowerFlow。多相流方面LBM模型主要有四种：着色模型、伪势模型、自由能模型以及其他模型（因为我主要的研究方向就是多相流）。采用LBM解决物理问题的基本步骤如下：

第二章 流体力学基本方程

大致看了一下第二章，感觉硬货来了，所以怂得我还是先看了附录A：笛卡尔张量的基本知识，把之前丢给流力老师的知识往回捡一捡。流力中物理量按其维数划分为标量、矢量和张量。三维空间中n阶张量可以由3^n个分量组成，所以标量其实是零阶张量，矢量是一阶张量，而二阶张量由9个分量表示。爱因斯坦求和约定：张量计算中同一项中如出现两个相同的下标，意味着对该指标求和，该指标称为哑指标。在方程同一项只出现一次的指标称为自由指标。没大看懂，特意百度了一下，看一下例子很好懂：

一个二阶张量可以看成是两个矢量的并矢。单位张量克罗内克定义及性质如下：

两个基矢量的点积便是单位张量。置换符号的定义和性质如下：

任何一个二阶张量都可以唯一地分解为一个二阶对称张量和一个二阶反对称张量的和。二阶张量的代数运算有相等、加减、点积和双点积。张量的微分、积分运算，主要包括梯度、散度、拉普拉斯算子以及高斯公式（这一块非常重要，需要着重理解）。

1. 梯度 （哈密顿算子）：
   其中标量进行梯度计算如下：
   矢量进行梯度计算如下:
   一个n阶张量的梯度为（n+1）阶张量
2. 散度：相当于哈密顿算子与张量的点积。
   其中矢量的散度为：
   二阶张量的散度为：
   
   一个n阶张量的散度为（n-1）阶张量
3. 拉普拉斯算子：先求梯度，再求散度：
4. 高斯公式：
   自此基础知识就这些，感觉还可以，正式开始第二章的内容。第二章主要是推导了流体力学中的总能方程和内能方程，这两个方程对于LBM非常重要。首先讨论了流体力学中平均尺度的选取问题，进而提出连续介质的概念，引出连续介质假设（对于气体系统是否可以作为连续介质处理根据气体分子间的平均自由程λ与气体所处空间的代表性尺度L决定，即努森数Kn=λ/L。当Kn≤0.001时就可以作为连续介质处理；对于液体近代研究表明一直到微米尺度的液体层厚度，液体依然可以作为连续介质处理）。流体力学中研究的最小物质实体是流体质点，他在微观上充分大，宏观上又充分小。流体运动与能量传递的基本方程是根据质量守恒、动量守恒和能量守恒定律所建立起来的连续方程、动量方程和能量方程。

雷诺输运方程

流体力学中有两种方法描述流体运动：拉格朗日法（跟随流体质点研究流体运动的方法）和欧拉法（着眼于空间坐标研究流体运动的方法）。雷诺运输方程的作用是将一个流体系统的拉格朗日变化率表示为欧拉导数，从而将针对系统的物理量转换到易于研究的欧拉参考系中。通过各种推导和简化得到的雷诺输运方程如下：

连续方程（质量守恒）

根据质量守恒，通过带入雷诺输运方程得到连续方程如下：
或

对于不可压缩流体（密度保持不变为常数）可进一步简化为：