重力加速度g

you1840

1) 实际系统当中的重力加速度g_p=9.8m/s2方向向下，请问这个转化到LB系统中时候，g_lb是取正还是取负？
2) 带有外力项的LB方程中，求宏观速度是否要加上外力项？如下：
p[x][y] =f[x][y][0]+f[x][y][1]+......f[x][y][8];
u[x][y] =(f[x][y][0]*e[0][0]+f[x][y][1]*e[1][0]+......f[x][y][8]*e[8][0])/p[x][y]+0.5*G[0];
v[x][y] =(f[x][y][0]*e[0][1]+f[x][y][1]*e[1][1]+......f[x][y][8]*e[8][1])/p[x][y]+0.5*G[1];
G[0]，G[1]表示重力加速度x和y方向的分量
谢谢

onesupeng

1、这个和你建立的坐标相关
2、方法不止一种，可参考Guo_PRE_2002的文章，我更喜欢这种方法

you1840

谢谢

tsfed

是不是应该 u[x][y] =(f[x][y][0]*e[0][0]+f[x][y][1]*e[1][0]+......f[x][y][8]*e[8][0]+0.5*dt*G[0]])/p[x][y];

DelvonLiu

重力加速度怎么 转换啊

luo@odu.edu

the acceleration g = a constant, then du = g dt is a constant in each time step. The key is to set the limit in the (max.) Mach number so the code won't blow up.

xiaozhi

还是不怎么明白，这里只是有重力项的影响，若有别的外力（不止重力一种力的作用），这种情况又该如何处理？是否将实际的g转化为格子单位下的Glb?

luo@odu.edu

g can be a sum of ALL forces. In LBE, forcing is local.

DelvonLiu

在多孔介质中流动 ，重力驱动，介质阻力不是固定值与速度有关，我怎么确定最大的mach数呢？

luo@odu.edu

I assume you are dealing Brinkmann flow. In such case, the max. velocity can be estimate. You can always do some estimation. If you don't know how to do that, then use the approach of trial and error to find out.

绿藻心lv

若考虑重力的话，g==9.8？？？是值得怀疑的。。。。具体大小可以通过无量纲数推导出来。。。基本不是9.8.。。。。。因为是加速度的量纲，转换到格子单位大小肯定是会变的

绿藻心lv

在考虑外力作用的模型时候，，，宏观速度大小的求法，需根据力的模型来考虑，，，参照guo的那本书，有详细介绍。