请教用郭老师的非平衡分布方法做poiseuille流遇到的问题

lakershzh981

用郭老师的非平衡态外推方法来处理压力边界，在做poiseuille流时出问题。下面是流场中的u/umax速度曲线图，不同的taub时，中心轴线上的最大速度会变化，而且中心轴线最大速度与理论值也有差异，请问会是什么原因呢？（实线为理论值，其他点为不同taub时对应的数值值）

lakershzh981

哎呀，这图发得有点对不起观众了。绿色格点taub=0.75，蓝色格点taub=1.0，粉色taub=1.5，红色实线为理论值，Re=10.0，请问可能是哪里出了问题呢，想了好久。

magicmadmiao

如何用非平衡态外推方法来处理压力边界啊，我看书上还是不怎么了解，
能用C++代码的形式举个例子吗？比如流场的入口。谢谢了

luo@odu.edu

You can take a look at the following paper:

F. Verhaeghe, L.-S. Luo, and B. Blanpain.
"Lattice Boltzmann modeling of microchannel flow in slip flow regime."
Journal of Computational Physics 228(1):147-157 (January, 2009).

Also the paper by Prof. ZL Guo (which was cited in the above paper).

luo@odu.edu

This is the problem due to the LBGK model, assuming you used the bounce-back
boundary conditions. You can take a look at the following paper:

L.-S. Luo, W. Liao, X. Chen, Y. Peng, and W. Zhang.
"Numerics of the lattice Boltzmann method: Effects of collision models on the lattice Boltzmann simulations."
Physical Review E 83(5):056710 (May 2011).

wdlxmzd

你这应该是程序有问题，在你提到的几个松弛因子下，不会有这么大的误差。