| 查看: 1775 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
LBM中压力梯度作用下流体的传输困扰 已有2人参与
|
||
|
各位朋友 你们好，本人正在学习LBM方法，想做一个压力驱动下的流体传输的模拟，从网上找了一些简单的程序，但都是速度边界的，我按着参考书里面给的方法和公式想把其改成压力边界，可程序一直无法运行（由于刚开始学LBM，也不知道怎么修正），有些文献中往往给出左侧边界压力是1，右侧压力边界是0之类的语句，这些都是怎么添加进入的，我在设置初始密度分布函数的时候，又该如何设置呢（很多程序都给定初始密度都为1啊）。希望看到的朋友能帮帮我。谢谢 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有117人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
4楼2017-06-22 14:59:36
2楼2016-10-01 21:00:54
|
左侧为1，右侧为0，这个应该是无量纲化了2如果你有实际模型，首先要进行无量纲化，然后再编程解决。压力边界的问题属于边界条件的设置，参考何雅玲教授的书，里面有很多边界条件设置方法，针对你的情况可以选择。LBM中确实初始平均密度都设为1，主要是为了保证精确度。 发自小木虫IOS客户端 |
3楼2016-10-12 12:50:36
【答案】应助回帖
|
看你选取的平衡态分布函数的格式，如果选取的是一般的伪不可压格式，其密度和压力有显示对应的方程，根据你的密度可以反解出密度边界条件，如果想直接带入的压力边界条件，可以考虑选取郭照立老师提出的D2G9模型，对了最关键的边界选用郭老师提出的非平衡外推格式，该格式可以直接将速度/压力（密度）边界条件带入到平衡太分布函数中，精度二阶且严格保证定速度/压力（密度）边界条件。 参考论文：Non-equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method |
6楼2019-04-23 20:24:20
投票: