求助：MRT松弛参数的理解

hellolegend

看郭老师的<<格子Boltzmann方法的原理及应用>>3.5.2节(MRT基本模型)，对其第一段的理解出现疑问，请指点。
黑色的字是照着一般思路写出来的，红色的字是疑问。
具体见下图。

[ 本帖最后由 hellolegend 于 2011-5-22 22:34 编辑 ]

feixiang9

可以这样理解：
在LBM的碰撞项设计的过程中，应该保证对守恒量的碰撞不便性，即，碰撞不改变守恒量。
碰撞过程可以写成：
m_i^+=m_i-s_i*(m_i-m_i^eq)
对于等温流动，质量和动量守恒。对于二维问题，通过M矩阵转化的分布函数矩中矩rho（m_0）, rho*ux（m_3）和rho*uy（m_4）对应的平衡矩与即为（m_0^eq=rho）, （m_3^eq=rho*ux）和（m_4^eq=rho*uy）。这样设计之后，与守恒量相对应的松弛参数就可以选取任意值，为了计算方方便，选取其值为0

hellolegend

谢谢，已理解

luo@odu.edu

Be careful with the relaxation rate for the momentum. When there is an external force, the relaxation rate should be "1".

feixiang9

对！
感谢罗老师的提醒，因为上面解释的是郭老师书中的章节，没有将作用力项考虑进去。