| 查看: 3923 | 回复: 25 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[交流]
LBM方法演化控制方程是否应该考虑外力项？已有5人参与
|
|||||
|
查阅了一些常见的LBM代码，其中演化方程： f=f+tau*(feq-f) 其中tau为弛豫时间，feq为稳态分布函数 然而在一个poiseuille流的模拟代码中，采用的演化方程考虑了外力项： f=f+tau*(feq-f)+dt*F 其中F为外力 考虑外力项对最终模拟计算的结果有怎样的影响呢？何时要考虑呢？ 本人LBM新手，希望与这方面的朋友相互交流讨论（QQ 2241790403） |
回复此楼» 收录本帖的淘贴专辑推荐
 资源收集 | 科研路上必备 |
» 猜你喜欢
Error in routine lambda (100): wrong # or too many modes
已经有0人回复
-
祈福国自然面上
已经有61人回复
-
物理学I论文润色/翻译怎么收费?
已经有243人回复
-
请问castep计算出的能带图横轴数据缩放因子怎么计算
已经有1人回复
-
Result in Physics 从见刊开始 一般多久可以被SCI收录
已经有26人回复
-
面青地会评时间？？？
已经有8人回复
-
科研求助
已经有1人回复
-
深圳大学高等研究院暑期夏令营
已经有10人回复
-
深圳大学高等研究院暑期夏令营
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 中国石油大学出来了
已经有33人回复
3楼2014-10-20 15:03:37
孤独的漫步
木虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 4461.4
- 散金: 677
- 红花: 2
- 帖子: 83
- 在线: 328小时
- 虫号: 2591375
- 注册: 2013-08-11
- 性别: GG
- 专业: 传热传质学
2楼2014-10-16 21:35:12
孤独的漫步
木虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 4461.4
- 散金: 677
- 红花: 2
- 帖子: 83
- 在线: 328小时
- 虫号: 2591375
- 注册: 2013-08-11
- 性别: GG
- 专业: 传热传质学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
你说的那是Zou/He速度压力边界吧，那是用于处理已知进(出)口密度(即压力)或者速度而采取的边界处理形式。 按我的理解 ，压力边界和外力项是统一的，即知道进出口压差，可以转化为外力加到演化方程里面去。 另，Poiseuille流是很好的基础检验程序，基本能把各个边界条件都涉及到，建议露珠再把非平衡外推，反弹等都试一遍，D2Q9模型就能满足。 对于这个Zou/He边界，个人建议，仅仅是个人建议！了解是怎么回事就好，一般能用这个边界方式处理的，非平衡外推都能处理，且非平衡外推在数学方程及其程序语言上都较Zou/He边界简单。 当然，不知露珠是临时涉及LBM还是要长久研究，很多边界都可以编程封装好，在不同case里套用即可。 |
» 本帖已获得的红花（最新10朵）
4楼2014-10-25 22:30:58
送红花一朵|
非常感谢您的回复 lz是新手，估计要长期使用LBM，主要希望通过学习这个边界条件加深LBM的了解。 您的理解是否是以下两种二选一： （1）采用外力项模型直接修正离散控制方程，采用普通的边界条件 （2）采用无外力项的离散控制方程，但边界条件采用Zou/He边界 a）目前依然在艰难摸索（2）中，已知流体进出口密度进行计算，并且在进出口均采用了Zou/He边界，同时假设进出口Y方向的速度均为零，不知在出口这样处理是否合适？ b）初始化时需要计算弛豫时间tau，通常tau是利用动力粘度niu=L*U0/Re得到的，但对于Zou/He密度（压力）边界，U0未知，是否是需要利用进出口密度差计算压力差G，再利用G求出动力粘度niu？ c）关于LBM计算流程的困惑，之前看过的几个算例都是遵循：迁移并碰撞→边界处理→计算宏观量 的步骤，也有：迁移→边界处理→碰撞→计算宏观量 ，其中计算宏观量这一步是否需要在边界上进行？对于Zou/He密度（压力）边界，在边界处理的过程中就能得到进出口边界上X方向的速度，那应该采用哪种计算流程呢？目前始终无法得到结果，怀疑在这一块我的代码也有问题。 困惑较多，望不吝赐教！ |
5楼2014-10-26 10:41:22
10