24小时热门版块排行榜    

查看: 1800  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Lionel9527

银虫 (初入文坛)

[求助] 格子Boltzmann单位转化 已有1人参与

如下是我对LBM单位转换的理解,不知正确与否,麻烦虫友能指点一二。
对于格子空间量的取值任意性很大,只要保证主要准则数(如Re)、几何形状比等相似准则不变。
对于一个长为L,直径为D,初始速度为u0的二维管道流动,取特征值为D,则无量纲长度为L/D,无量纲直径D'=1,取格子数为N+1,则被分为N个区间,dx=L/D/n
由于Ma=Re*(tau-0.5)*dx/1.732<0.3(根号3的值,此公式可根据Ma定义结合粘度公式推导得出)且雷诺数、dx,已知可得到tau的取值范围,结合tau的经验范围(0.5,2)可取一个初值tau,任取流体的格子速度Ulb=0.1 密度rho=1,保证雷诺数值不变得到粘度值,再由粘度公式v=C*C(tau-0.5)*dt/3,及C=dx/dt,得出dt的值,至此基本变量dx,dt,Ulb,v等量已知,具体步骤如下图。
    看似合理但又感觉这样的计算存在问题,首先对于非稳态传热问题,需要画出无量纲时间(通常用Fo傅里叶数表达)与温度的关系图,可是怎么联系格子时间与F0,格子单位下一次碰撞就是经历了dt时间,而对于非稳态问题 需要在一个时间间隔Dt内迭代N次达到准稳态再计算下一个时间间隔,以此类推。希望大家能指点我一二, 。

格子Boltzmann单位转化
1.jpg


格子Boltzmann单位转化-1
2.jpg
回复此楼
Doit
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feilin0315

新虫 (初入文坛)

【答案】应助回帖

http://www.opencimulation.org/?p=1517
3楼2015-08-10 08:41:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

JBXiao1th007

新虫 (小有名气)

LBM的问题,你可以去请教一下西交何老师的学生,呵呵
2楼2015-01-23 15:49:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dream070141

捐助贵宾 (小有名气)

LZ对单位转换清楚了吗 刚接触LBM 对这问题一直不明白
4楼2016-01-10 00:04:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haoxinqing15

新虫 (正式写手)

引用回帖:
2楼: Originally posted by JBXiao1th007 at 2015-01-23 15:49:55
LBM的问题,你可以去请教一下西交何老师的学生,呵呵

何老师哪位学生?

发自小木虫Android客户端
5楼2016-01-10 09:19:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见