Poiseuille流的LBM模拟问题

孤独的漫步

最近学习LBM模拟，想以Poiseuille流练习练习，从静止到充分发展，现遇到如下问题：

试过两种方法，第一，外力驱动，固体边界为非平衡外推方式处理，出入口边界以周期性边界方式处理，模拟正常；

第二，压力差驱动，（后面的理解就不知道对不对了）确定出入口密度，固体边界还是非平衡外推方式（这里是速度已知，密度以内一层密度代替），出入口也想以非平衡外推方式，求平衡态的密度已知，速度由内一层代替，但是计算不出正确结果。。。

另，提取了以外力驱动时的整个计算场的密度，发现不存在密度梯度，这样我理解的第二种方法，以压力差驱动，即以出入口密度差驱动貌似就有问题了。。。

有木有人可以指导一下啊。。。（目前参考的是郭照立老师的那本书）

[ 本帖最后由 孤独的漫步 于 2014-1-30 14:42 编辑 ]

seanwillian

换D2G9格式试试。按我的理解，若使用D2Q9的话，初始化的时候密度自入口到出口应该有梯度差的。以外加作用力方式驱动的话，密度肯定不会有梯度啊。BGK本身描述的是不可压流体，如果密度变了，就有问题了。

孤独的漫步

之前的问题想通了
现在假设来流已经稳定，且呈抛物线型，用非平衡外推，密度用内一层的（这里貌似有问题），然后发散了，这又是为什么？

lwd1981

可以由宏观矩来给边界条件，速度周期，对应点上压力分布存在压力差，然后可以采用非平衡外推即可。

onesupeng

压力差驱动，驱动结果体现在压力的梯度上

而用压力梯度驱动，相当于把压力分成均匀部分+体积力，压力梯度已经作为外力施加上去，所以你看到的密度是均匀的，即没有压力梯度。在处理周期性边界流动的时候，用施加体积力比较妥当

[ 本帖最后由 onesupeng 于 2014-2-6 22:27 编辑 ]

zhy888

压差驱动肯定是密度已知，就不能用内层的代替啦，是速度未知，需要采用外推

luo@odu.edu

I have some questions:
1) You indicated: "换D2G9格式试试。"  
What is "D2G9" model?

2) You stated: "BGK本身描述的是不可压流体，如果密度变了，就有问题了。"
Could you please explain why "BGK per se describes incompressible flow"?
Any references?

-- LS Luo

u_zwh

luo@odu.edu 发表于 2014-5-1 19:30
I have some questions:
1) You indicated: "换D2G9格式试试。"  
What is "D2G9" model?

D2G9格式是指郭照立老师提出的不可压模型，Journal of Computational Physics 165, 288–306 (2000)，Lattice BGK Model for Incompressible Navier–Stokes Equation。

luo@odu.edu

u_zwh 发表于 2014-10-19 20:47
D2G9格式是指郭照立老师提出的不可压模型，Journal of Computational Physics 165, 288–306 (2000)，Lat ...

LBE is NOT a truly incompressible model, for it does NOT solve the Poisson equation exactly -- the speed of sound in LBE is FINITE; in fact, it is very small.