不可压格子波尔兹曼模型-D2G9

D2Q9模型是一种弱可压模型，郭老师在此基础上提出了D2G9模型，是一种不可压模型，两个模型主要有以下不同点：

控制方程：

1、平衡态方程定义不同

 对于D2Q9：

，

对于D2G9：

，，，

其中S_i(u)定义为：

 

2、演化方式一个基于密度演化（D2Q9），一个基于压力演化（D2G9），体现在宏观量求解量不同，且定义不同

对于D2Q9：

，

对于D2G9：

，

两者宏观量计算时，下标的开始位置不一样，D2Q9中包含i=0，而D2G9中则不包含，如果编程时，包含了i=0，则计算出现NAN。

参考文献：

Guo Z, Shi B, Wang N. Lattice BGK Model for Incompressible Navier–Stokes Equation[J]. Journal of Computational Physics, 2000, 165(1):288-306.

 

转载于:https://www.cnblogs.com/kljfdsa/p/7978839.html