关于流场中极小直径颗粒所受曳力的问题
请教各位，我最近在做颗粒流的模拟，遇到一个问题
根据斯托克斯公式，当雷诺数小于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
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
你忘记了delta t。这个问题如果你想用显式差分求解收敛必然要求一个极小的时间增量。这就能保证delta u是逐渐降为0的了。另外当r极小时，这个问题的尺度也将发生变化。动能与热涨落之间的对比可能发生反转，从而导致颗粒在强制对流中作有偏随机行走（布朗运动）
，
解的不就是dv/dt=9/2*μ*(u0-v)/r^2这个方程吗？u0是流体速度，v是颗粒速度
你再好好想想delta t是不是要等于1吧，这是个分分钟就能算出来的题，直接解微分方程也可以