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.