说明：方法来自外文资料，主要通过无量纲方法进行LBM仿真计算中由实际物理量到格子物理量的转化，保证其一致性
已知量
物理长度尺寸
H
H
H ；
格子长度尺寸 （自定义）
H
~
\tilde{H}
H~；
密度
ρ
\rho
ρ ；
格子密度
ρ
~
\tilde{\rho }
ρ~（自定义，常取为1）；
流体运动粘度
v
v
v；
实际
Re
=
u
m
H
v
\operatorname{Re}=\frac{{{u}_{m}}H}{v}
Re=vumH ；
格子时间步
t
~
\tilde{t}
t~（常取1）；
格子空间步
x
~
\tilde{x}
x~（常去1）；
物理速度量
u
u
u 。
待求量及求解顺序：
1.长度换算因子：
C
H
=
H
H
~
{{C}_{H}}=\frac{H}{{\tilde{H}}}
CH=H~H ；
2.密度换算因子：
C
ρ
=
ρ
ρ
~
{{C}_{\rho }}=\frac{\rho }{{\tilde{\rho }}}
Cρ=ρ~ρ；
3.声速换算因子：
C
s
2
=
1
3
Δ
x
2
Δ
t
2
{{C}_{s}}^{2}=\frac{1}{3}\frac{\Delta {{x}^{2}}}{\Delta {{t}^{2}}}
Cs2=31Δt2Δx2（常取1）；
4.时间换算因子：由
v
=
(
τ
−
1
2
)
C
s
2
Δ
t
=
(
τ
−
1
2
)
1
3
Δ
x
2
Δ
t
=
(
τ
−
1
2
)
3
C
H
2
C
t
v=\left( \tau -\frac{1}{2} \right){{C}_{s}}^{2}\Delta t=\left( \tau -\frac{1}{2} \right)\frac{1}{3}\frac{\Delta {{x}^{2}}}{\Delta t}=\frac{\left( \tau -\frac{1}{2} \right)}{3}\frac{{{C}_{H}}^{2}}{{{C}_{t}}}
v=(τ−21)Cs2Δt=(τ−21)31ΔtΔx2=3(τ−21)CtCH2得出求解公式：
C
t
=
(
τ
−
1
2
)
3
C
H
2
v
{{C}_{t}}\text{=}\frac{\left( \tau -\frac{1}{2} \right)}{3}\frac{{{C}_{H}}^{2}}{v}
Ct=3(τ−21)vCH2 ；
5.速度换算因子：
C
u
=
C
H
C
t
{{C}_{u}}=\frac{{{C}_{H}}}{{{C}_{t}}}
Cu=CtCH ；
6.格子速度：
u
~
=
u
C
u
\tilde{u}\text{=}\frac{u}{{{C}_{u}}}
u~=Cuu ；
计算结束后需要通过前后的Re数等进行一致性验证，另需注意在LBM中格子速度不宜超过0.3，为了满足这一要求，可通过减小长度换算因子
C
H
{{C}_{H}}
CH会导致计算复杂度增加）或者减小松弛因子
τ
\tau
τ会导致LBM变得不稳定），因此需要把握减小的度以获得满足条件的格子速度。
文档来源
《Unit conversion in LBM》
文档链接：
链接：https://pan.baidu.com/s/1_ype63aKxgLZ6Ut9MOm8lg&shfl=sharepset
提取码：1sj3
6551
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
