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

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

控制方程:

1、平衡态方程定义不同

 对于D2Q9:

对于D2G9:

其中Si(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.

 

posted @ 2017-12-04 20:21  药否  阅读(668)  评论(0编辑  收藏  举报