求助：Poiseuille压力边界处理，附代码

zgsydx

求助：单相，二维Poiseuille流，采用D2Q9模型，进口和出口压力恒定。怎么处理进出口比较好呢？
现在尝试：
1）进口和出口的未知量用平衡态处理。
2）用ZOU&He的压力边界条件。
我绘制了ux和uy的图像，但是结果都不怎么好。
其中，用Zou&He边界代码（进口，rho是密度，分布函数都是入口边界上格点的分布函数）：
rho=3*p_in;
u=1-(f9+f2+f4+2*(f3+f7+f6))/rho;
f1=f3+2/3*rho.*u;
f5=f7-1/2*(f2-f4)+1/6*rho.*u;
f8=f6+1/2*(f2-f4)+1/6*rho.*u;
以上代码没有处理4个角点，即使处理过还是有问题。

请各位大神指教！

luo@odu.edu

Two important points:

1) You MUST maintain a constant flow rate, or the total mass in the flow domain will vary.
2) The constant flow rate must be in such a value that the maximum velocity satisfies the "incompressible" condition Ma = U_max / c_s < 0.3 (or even smaller).

-- Li-Shi Luo

luo@odu.edu

Two ways to get it right:

1) For velocity conditions (for ideal incompressible flow), the ANALYTIC velocity profile can be obtained, and that should be used at the entrance.

2) For pressure conditions (for compressible flow), please see 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).

hust_illusion

Dear Luo,

对于压力边界（密度），采用Luo 09 JCP的思想对进出口进行外推时，是否需要始终保持进出口密度恒定为 rho-in和rho-out？

我的计算顺序为：
所有点collision--所有点stream--处理上下壁面边界--进出口边界外推--计算宏观量（所有点）--进出口密度赋值为rho-in和rho-out（覆盖上一步计算出的密度）

但是计算还是有问题，无法收敛，速度会越来越大。

[ 本帖最后由 hust_illusion 于 2014-9-5 09:01 编辑 ]

muzhiyu

请教下最后怎么解决的啊