物理参数到LBM格子参数的转换方法--注意雷诺数的L取值与粘度单位换算！

还是老规矩先宣传一下QQ群群： 格子玻尔兹曼救星：293267908。为了早点毕业建的群。

1、Re=U*L/v，这里L取值很有趣。 The definition of the Reynolds number is not unique: on the one hand, some people may choose L for the length, others for the width of the considered system. On the other hand, U may be the average velocity or the maximum velocity—depending on the person defining the Reynolds number. This has to be kept in mind when comparing Reynolds numbers from different sources. --Timm Kruger. 书P269.

   一般来说，槽道流（如泊肃叶流）取槽宽，也就是短边长；圆柱绕流一般取直径（U这时指的是流入面的平均流速），还有一些特殊的，比如粗糙度不同的裂隙，得取等效水力宽度。具体情况具体分析吧...压力驱动和体力驱动，Re另有求法！

2、动力粘度的单位是pa.s，运动粘度标准单位是m^2/s。二者的换算关系是：运动粘度＝动力粘度/密度。
密度单位为Kg/m^3 。推导: pa·s/(kg/m^3)=N/m^2*s*m^3/kg=N/kg*m*s,由F=ma这个公式我们可以知道这个单位换算m/s^2=N/kg,所以有pa·s/(kg/m^3)=m^2/s。角加速度Va的单位是弧度/秒，而弧度的转化参数永远是1.因此C_Va=1/C_time。再举个2D例子求力和压强：

% 2021.2.14-------------------------
tau=0.65;% From: The Immersed Boundary-Lattice Boltzmann Method for Solving Fluid-Particles Interaction Problems
nu=(tau-0.5)/3 ;
P_viscosity=1e-6 ; %unit: mm2/s. 1.00*10^(-6) (m2/s)%  water Kinematic viscosity=1.006*10^(-6) (m2/s). air: 15.5e-6.
C_viscosity=P_viscosity/nu; %粘度的转化参数
C_L=1e-4; %0.001m=10 lattices*C_L %长度转化参数。
C_time=C_L^2/C_viscos