| 查看: 925 | 回复: 6 | ||
[求助]
关于流场中极小直径颗粒所受曳力的问题 已有1人参与
|
|
请教各位，我最近在做颗粒流的模拟，遇到一个问题 根据斯托克斯公式，当雷诺数小于0.1时，球形颗粒所受曳力为Fd=6*PI*μ*r*Δu，对于极小直径颗粒来说由于其直径较小，基本都是雷诺数小于0.1的 这样一来在下个时间步，曳力反馈给颗粒速度的值即为： Δv=aΔt=F/m*Δt=(6*PI*μ*r*Δu)/(ρ*4/3*PI*r^3)*Δt ; 由于颗粒半径r极小，上述公式化完之后分母上的r^2会造成反馈给颗粒的速度值极大，甚至超过了本来的颗粒-基液之间的速度差。 请教大家，这种情况该如何处理？ （我曾经想增大网格尺度使颗粒直径大于1，但是网格太大计算机算不了了。。。@wuming524 |
回复此楼» 猜你喜欢
急需调剂
已经有6人回复
-
RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一
已经有18人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有7人回复
-
291求调剂
已经有8人回复
-
271求调剂
已经有36人回复
-
0854调剂
已经有9人回复
-
322求调剂
已经有15人回复
-
085404 298分求调剂
已经有13人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有24人回复
-
327求调剂
已经有29人回复
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
【答案】应助回帖
感谢参与，应助指数 +1
|
你忘记了delta t。这个问题如果你想用显式差分求解收敛必然要求一个极小的时间增量。这就能保证delta u是逐渐降为0的了。另外当r极小时，这个问题的尺度也将发生变化。动能与热涨落之间的对比可能发生反转，从而导致颗粒在强制对流中作有偏随机行走（布朗运动） 发自小木虫IOS客户端 |
2楼2018-07-10 23:05:26
|
谢谢帮助。在计算的时候，一开始我是进行了量纲转化的，计算时公式里的量也是用的转化后的量。比如，我的模型里一个时间步长（dt）等于现实中的1.2×10∧-5次方秒。您是指在计算时，这一步需要乘上1.2×10∧-5次方吗？还是就乘上转换后的dt(1秒)？再次感谢帮助 发自小木虫IOS客户端 |
3楼2018-07-10 23:28:42
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ウンガTA: 金币+15, ★有帮助 2018-07-13 16:37:08
ウンガTA: 金币+15, ★有帮助 2018-07-13 16:37:08
|
你看，你的这个delta t是对加速度求导的差分格式得到的，所以它和你公式的其他部分没关系，是可以随意取的一个小值。我不太明白你所说的1.2e-5s是怎么来的。跟时间相关的差分问题都对时间步长有要求，特征长度越小，时间步长就需要越小才能收敛 发自小木虫IOS客户端 |
4楼2018-07-11 00:48:04
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
5楼2018-07-11 00:54:33
|
是这样的，我的计算就是按每个时间步长一步步进行的，每个时间步长为dt。我算的是一个时间步长内，这个力对速度的改变。在一个时间步长内，力对速度的作用时间为dt，所以在计算过程中 delta t就是1。 发自小木虫IOS客户端 |
6楼2018-07-11 16:06:39
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
7楼2018-07-11 16:53:27
