找回密码
 注册
查看: 5845|回复: 8

Poiseuille流的LBM模拟问题

[复制链接]
发表于 2014-1-30 14:39:45 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

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

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

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

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

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

[ 本帖最后由 孤独的漫步 于 2014-1-30 14:42 编辑 ]
发表于 2014-1-30 19:00:25 | 显示全部楼层

回复 1# 孤独的漫步 的帖子

换D2G9格式试试。按我的理解,若使用D2Q9的话,初始化的时候密度自入口到出口应该有梯度差的。以外加作用力方式驱动的话,密度肯定不会有梯度啊。BGK本身描述的是不可压流体,如果密度变了,就有问题了。
 楼主| 发表于 2014-2-6 16:53:07 | 显示全部楼层
之前的问题想通了
现在假设来流已经稳定,且呈抛物线型,用非平衡外推,密度用内一层的(这里貌似有问题),然后发散了,这又是为什么?
发表于 2014-2-6 20:03:34 | 显示全部楼层
可以由宏观矩来给边界条件,速度周期,对应点上压力分布存在压力差,然后可以采用非平衡外推即可。
发表于 2014-2-7 06:24:00 | 显示全部楼层
压力差驱动,驱动结果体现在压力的梯度上

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

[ 本帖最后由 onesupeng 于 2014-2-6 22:27 编辑 ]
发表于 2014-5-1 16:57:51 | 显示全部楼层
压差驱动肯定是密度已知,就不能用内层的代替啦,是速度未知,需要采用外推
发表于 2014-5-1 19:30:49 | 显示全部楼层

回复 2# seanwillian 的帖子

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
发表于 2014-10-19 20:47:19 | 显示全部楼层
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。

点评

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.  详情 回复 发表于 2014-10-24 20:54
发表于 2014-10-24 20:54:38 | 显示全部楼层
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.
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表