空气动力学仿真技术：格子玻尔兹曼方法(LBM)：LBM的数学基础.pdf

空气动力学仿真技术：格子玻尔兹曼方法(LBM)：LBM的数

学基础

1空气动力学仿真技术：格子玻尔兹曼方法(LBM)：LBM的

数学基础

1.1格子玻尔兹曼方法简介

1.1.1LBM的发展历史

格子玻尔兹曼方法（LatticeBoltzmannMethod,LBM）起源于20世纪80年

代末，最初由FrancescoHiguera和Rapoport在1986年提出，称为Higuera-

Rapoport模型。随后，Shan和Chen在1993年引入了流体动力学相互作用，使

得LBM在流体动力学领域得到了广泛应用。LBM的发展历程中，它逐渐从一个

理论模型转变为一种强大的数值模拟工具，特别是在处理复杂流体动力学问题

时，如多相流、微流体和高雷诺数流动。

1.1.2LBM的基本原理

LBM基于玻尔兹曼方程，但将其简化并离散化到一个格子系统中。这种方

法将流体视为由大量粒子组成的集合，每个粒子在格子上遵循特定的运动规则。

LBM的核心在于其离散化的速度空间和时间空间，以及粒子的碰撞和传输过程。

离散化速度空间

在LBM中，速度空间被离散化为有限个方向，每个方向对应一个速度向量。

最常见的离散化模型是D2Q9模型，其中D2表示二维空间，Q9表示九个速度

方向。这些方向覆盖了所有可能的粒子运动方向，包括沿格子的轴向和对角线

方向。

离散化时间空间

时间也被离散化，系统在每个时间步长内更新粒子的分布函数。分布函数

描述了在特定位置和速度方向上粒子的密度。

碰撞和传输过程

在每个时间步长内，粒子在格子上进行碰撞和传输。碰撞过程通过玻尔兹

曼碰撞算子来描述，但在LBM中，这个过程被简化为一个局部的平衡态恢复过

1

程。传输过程则是粒子沿着其速度方向移动到相邻的格子点。

数学描述

LBM的基本方程是离散化的玻尔兹曼方程，可以表示为：

+,+−=

其中，是粒子在位置、时间和速度方向上的分布函数，是碰撞算子，

描述了粒子在碰撞过程中的变化。

简化碰撞算子

在LBM中，碰撞算子通常被简化为Bhatnagar-Gross-Krook(BGK)算子，

可以表示为：

1

=−

其中，是松弛时间，是局部平衡态分布函数，它依赖于流体的宏观属

性，如密度和速度。

宏观属性的计算

LBM中的宏观属性，如密度和速度，可以通过分布函数的加权平均来计算：

​

,=,

​

1

,=,

,

边界条件处理

LBM在处理边界条件时，通常采用反射或滑移边界条件。反射边界条件意

味着粒子在遇到固体边界时会反弹回来，而滑移边界条件则允许粒子在边界上

滑