波尔兹曼算法lbm单位换算.docx

对于无量纲系统

考虑无量纲系统：假定格子数为N,时间步数为Nt,离散空间和时间尺度苏和员分别给出如下：6x=1/N，流=\/Nt.

此时，Up和ui的关系给出如下：

up=dxdtu\,Cp=dxdtCi.

同时，我们有

Vp=1/Re=3x23tvi.

那么，格子速度和黏性决定如下：

ui=6t6xup,vi=6t6x21Re.

由于无量纲速度为1,所以我们有

Up=1,u]=3t3x.

对于，真实物理系统如法炮制，只需根据物理Up，重新定义Vp，然后得到ui.

真实物理系统

网友问题

物理参数为：Pin=101325PaPout=101225Pa,压强差为100PaLx=25mmLy=1mm,求出特征速度

=13.9m/s,Re=920空气粘度为1.5X(10A-5)，划分为250X10

问题：如何确定格子参数？

1.确定特征尺度，特征速度和Re：

对于压力驱动的流动，对于此问题(Posieuilleflow)，如果为层流，由压力差可以确定特这速度

U(y)=APp2n((L/2)2-y2),y=±L/2(L=Ly)

由此可以计算出入口的平均速度Umean,p.

Umean,p=』L/2—L/2U(y)dy/L.

2.首先确定格子解析度

假定格子解析度为N，那么物理空间网格尺度为dx=L/N.根据苏可以确定其他方向的格子解析度。

3.压力差转换.在巳知了Pp(物理密度)，压力差的转换关系可以下面来定义

APpPpUmean,p/2=Ap]piU2mean,l/2

4.确定物理时间步长.

假定物理时间步长为步，(Cp和C]为物理声速和格子声速)，那么物理速度和格子速度的关系为up=3x3tui.

如果也考虑格子声速和物理声速的关系，也需要下列关系成立

Cp=WCi.

由Re可得物理黏性和格子黏性之间的关系，

Vp=Umean.pLRe=^x2^tvl.

如果关心Ma数或者弱可压缩性，可以根据物理Ma数确定格子速度如下：

ui=upCpCi.

否则，你可以自己定义低Ma对应ui减小可压缩性误差.如此，时间步长可以确定

3t=3xu]up.

实际物理单位：L（长度），U（速度），C（声速），v（黏度），T（时间）

对应的格子单位：Li,Ui,Cs（格子速度），vi，T1

对应关系：Lr=LLl，Ur=UUl=CCs，Vr=vvl，Tr=TTl

因为Re相同，所有：Re=^v=UlLlvl（1），对于一个巳知的格子模型，巳知量有Li（等于N*X，3x一般

取1,所以Ll等于N），vl（可由式子v=Cs（Tp0.5）获得）,v（可查表获得）

且有（1）式可知

vvl=ULUlLl=UrLr（2）

Ur=UUl=CCs（3）该式子的结果是巳知的，因为实际声速可查表获得，而格子声速也是固定的。

将（3）式代人（2）,就可得Lr

Lr=TTl=LULlUl=LrUr（4）

将⑵式左右同除以（3）得平方，就可得到（4）,即Lr

到此为止，Ur，Lr，Tr都巳知道，所以只要知道计算模型的格子长度，格子时间，格子速度就可知道实际物理模型中的值。