找回密码
 注册
查看: 2510|回复: 6

请教MR(Multireflection)格式中的修正项

[复制链接]
发表于 2011-12-20 15:39:12 | 显示全部楼层 |阅读模式

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

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

x
在Luo教授An evaluation of lattice Boltzmann schemes for porous medium flow simulation,Computers & Fluids 35 (2006) 898–909文章中,MR格式中的修正项f(2)=-M^-1*S*( t - t_eq ),文中给出了t的表达式,请问t_eq的表达式是什么啊?对于DmQn这种形式有没有通用的表达式?

另外我发现这里面的一个小错误,对照文献Multireflection boundary conditions for lattice Boltzmann models,Luo的文章MR中k3=k5=q^2 / (1+q)^2应该是k3= - k5=q^2 / (1+q)^2。

An evaluation of lattice Boltzmann schemes for porous medium flow simulation.pdf

385.93 KB, 下载次数: 249

 楼主| 发表于 2011-12-20 15:55:45 | 显示全部楼层
我觉得应该可以由Multireflection boundary conditions for lattice Boltzmann models里面的式子推导出来,但是这篇文献中的式子太复杂了,看不懂,特请教大家。
 楼主| 发表于 2011-12-23 09:05:54 | 显示全部楼层
再不顶一下这问题就要石沉大海了。
求助求助~~~
 楼主| 发表于 2011-12-24 13:35:44 | 显示全部楼层

结局

经过与Luo教授邮件联系,他非常友善地回答了我的问题,对此问题也感兴趣的同仁可以参考下,以下是他的答复:

t_eq are given in Eqs. (9): they involve only q's and m's.

You may be correct about the formula for k_3 and k_5 -- it may be a mistake of
ours.
发表于 2012-1-3 20:49:36 | 显示全部楼层
I may also add that MR-BCs may not benefit that much. If you can use 2nd-order interpolations
and MRT correctly, it should work well enough.
 楼主| 发表于 2012-1-4 14:26:45 | 显示全部楼层

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

Thank you very much for your kind suggestion!
Happy new year!
发表于 2012-2-8 10:41:17 | 显示全部楼层
Prof. Luo说的对,MR并不像乍看起来那么完美。如果仅仅模拟Poiseuille流,MR非常非常出色。但是在动边界中,MR不如Lallemand&Luo的二次插值稳定。这应该和MR的推导过程有关的。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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