|
发表于 2010-11-10 16:26:22
|
显示全部楼层
回复 6# ywang 的帖子
有外力作用下,还是需要转换的吧?看了Jonas Latt的 choice of units in lattice boltzmann simulations, 越来越糊涂了。我想算REV尺度下多孔介质的流动,请教一下:
几何模型是 0.3mm( y方向 )*3mm( x方向 )的一个流动区域,流动介质是一种牛顿流体,其rho=910kg/m^3, eta=2480Pa-s,nu=eta/rho=2.72m^2/s,计算poiseuille flow, 压力驱动,进出口压差为2MPa,流动沿X方向,那么根据poisuille的解,
外加压力产生的加速度应为: g=(dp/dx)/rho=7.32e5m/s^2
平均速度 u=g*(y/2)^2/(3*nu)=2.02mm/s
时间 t=3/2.02=1.485s
Re Re=2.23e-7
取无量纲因子 x0=0.03mm, u0=2.02mm/s, t0=0.01485s
无量纲化: xd=0.3/0.03=10, yd=3/0.03=100, td=1.485/0.01485=100, u-d=1,
nu-d=4.48e7
假设:10*100个格子,计算步数取1000,则dx=1,dt=0.1,
根据Choice~~
c(lattice speed)=10,u-lb=u-d*dt/dx=0.1,nu-lb=4.48e6,
因为 Cs^2=c^2/3=100/3 以及 nu-lb=Cs^2(TAU-0.5)*dt
所以TAU=3*nu-lb/10+0.5=1.344e5
同时:g-lb(LB单位的外力加速度)=u-lb*3*nu-lb/(y/2)^2=5.376e5
这样推算正确吗?
[ 本帖最后由 arlu 于 2010-11-10 16:37 编辑 ] |
|