| 查看: 1794 | 回复: 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次达到准稳态再计算下一个时间间隔，以此类推。希望大家能指点我一二， 。  1.jpg  2.jpg |
回复此楼» 猜你喜欢
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Lattice boltzmann modeling 格子boltzmann模型 英文版
已经有96人回复
- 流场计算方法-格子Boltzmann方法LBM的源代码
已经有14人回复
- 防线菌的对数期前期如何确定？急急急！！！
已经有4人回复
JBXiao1th007
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 107.9
- 帖子: 152
- 在线: 28.1小时
- 虫号: 3356333
- 注册: 2014-08-07
- 专业: 工程热物理相关交叉领域
2楼2015-01-23 15:49:55
feilin0315
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 3
- 帖子: 4
- 在线: 2.6小时
- 虫号: 3658655
- 注册: 2015-01-22
- 性别: GG
- 专业: 工程热物理相关交叉领域
3楼2015-08-10 08:41:36
4楼2016-01-10 00:04:38
haoxinqing15
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3179.7
- 散金: 10
- 红花: 1
- 帖子: 950
- 在线: 132.9小时
- 虫号: 4106204
- 注册: 2015-09-27
- 性别: GG
- 专业: 等离子体物理
5楼2016-01-10 09:19:42
100