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

求助,进出口速度不一样

[复制链接]
发表于 2011-3-11 23:26:06 | 显示全部楼层 |阅读模式

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

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

x
最近编了一个计算程序,想模拟矩形通道内的流动(square  duct),入口和出口条件都是压力边界条件,但发现了一个问题,在迭代第一次时,进出口速度基本相同,但迭代50次后进出口速度差了1000倍以上,进口速度稳定在10E-2的量级,出口速度却是10E-6量级...我用的D3Q15模型,1.3的入口密度,1.0的出口密度,相对动力粘度u=0.2,tua值随着密度的变化而计算...网格80X30X30...请问各位大虾...有没有什么好的建议解决这个问题?我菜鸟,勿笑...
发表于 2011-3-14 21:24:29 | 显示全部楼层
Have you checked mass and momentum conservation locally and globally? If not, why?!
 楼主| 发表于 2011-3-15 22:13:34 | 显示全部楼层

回复 2# luo@odu.edu 的帖子

I think this is a privilege I am lucky to have, thank you very much  Prof. Luo.

At the beginning of the implementation, I do not think mass and momentum should be in conservation, because it is the pressure boundary condition. I initialize the density in square duct with the inlet density , after streaming, the mass cannot be conservated in duct,I think.

If I initialize the density in duct with (inlet density+outlet density)/2.0, can mass and momentum be in conservation?

Moreover,in my code, I have checked every distribution and if it appeare negative density ,the implementation will stop.

My English is poor,sorry...

[ 本帖最后由 血色蔷薇 于 2011-3-16 08:10 编辑 ]
发表于 2011-3-16 20:34:29 | 显示全部楼层
You can respond to my comments in Chinese if that's more convenient to you. But I cannot reply in Chinese because it's too slow for me to type Chinese on my laptop.

You should check the mass influx and outflux, which must be balanced. The mass flux MUST be equal at each and every cross-section of the duct. Also, when the pressure BCs are imposed, only the averaged pressure, but NOT the profile, should be imposed. This is discussed in 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).

I believe some paper(s) by Prof. Guo ZL et al. also discuss this issue.

Good luck.
发表于 2011-3-16 20:39:14 | 显示全部楼层
PS the mass flux is equal only when the "steady" of final state is attained. You should ALWAYS check mass and momentum conservations no matter what. Why should the pressure BCs allow violation of conservation laws?
 楼主| 发表于 2011-3-17 18:30:51 | 显示全部楼层

回复 4# luo@odu.edu 的帖子

非常感谢 Prof. Luo,我会学习您所提及的文献的,有问题的话还希望您能解答
 楼主| 发表于 2011-3-18 13:09:07 | 显示全部楼层

回复 5# luo@odu.edu 的帖子

哦...Prof. Luo,我在学校的数据库中没有找到您提及的文章Lattice Boltzmann modeling of microchannel flow in slip flow regime,能麻烦您发一份给我吗?不胜感激!
wizard_0_0_0_0@163.com
发表于 2011-3-18 23:21:33 | 显示全部楼层
You can find ALL the paper with my name on my website:

http://www.lions.odu.edu/~lluo/journal-pubs.html

I also upload the paper in this message.

Verhaeghe_JCP-2009.pdf

396.08 KB, 下载次数: 524

 楼主| 发表于 2011-3-19 21:30:03 | 显示全部楼层

回复 8# luo@odu.edu 的帖子

谢谢,真的很感谢!
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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