OpenCFD2D RM不稳定性数值模拟教学示例.pdfVIP

RM不稳定性数值模拟

参考文献：

1.，薰，：Richtmyer-Meshkov失稳的数值模拟，计算

18（5）：390-396，2001

2.JacobsWJ,Shock-inducedmixingofalight-gascylinder,J.FluidMech.1992,234:

629-649

3.PiconeJMetal.Vorticitygenerationbyshockpropaionthroughbubbina

gas.J.FluidMech1988,189:23-51

1.网格设置

网格采用直角坐标均匀网格。计算域为

0xxL；0xyL

2.初值设置

如上图所示，柱形界面半径为r0,圆心位于（x0,y0）;

激波位于x=x1处；

计算与及初值设置如图2所示；初始时刻流场设置如下：

激波Mach数1.093，柱型界面内外气体密度比为0.138(柱内为轻气体)，

柱力、速度均相同。

初值设定:

方式1）激波不动，圆柱向右运动

该方式下，激波不动，波前（左侧）量（纲）为：

1=1，u=1,v=0,T=1；

111

根据激波关系式，可知波后（右侧）量为：

2

u22+（-1）M

2

u1（+1）M

2

+

2（1）M

=

2

2+-1

1（）M

-1222

（1+M）（M−1）

T22−1

T2

1（+1）2

M

2（-1）

计算出波后的量值。

波前圆柱内：

，，；

u1,v0ccc=0.138T1/0.138c

方式2）圆柱不动，激波向左运动

与方式1）相比，密度、均不变；速度u在方式1）的基础上叠加

上-1；

方式3）激波-圆柱都运动，激波运动速度为ut

该方式相当于，坐标系相对运动速度为Ut

(1)

u=u−u

t

根据经验Ut设置为激波后的运动速度比较好，这样，激波扫过圆柱后，圆柱

移动的速度为0，这样计算域就无需设置很长，有利于减小计算量。

3.计算步骤

1）生成初值

编译并运行init-RM.f90(可根据情况修改init-RM.in);

ln–sflow2d0.datopencfd2d.dat

注：

Init-RM.in中的x0,y0,r0为圆柱圆心坐标及半径；x1为激波位置；

Ma为激波Mach数，dr为圆柱内气体密度

ut为坐标系相对运动速度（见上文）。

deltr为圆柱界面过渡区的厚度（过渡区，密度从dr过渡到1）；

2)编译opencfd2d-1.5.5,将可执行文件拷贝入工作目录；

3)运行

mpirun–np8./opencfd2d-1.5.5.out

4）后处理

ln–sO