激波问题经典数值算例.doc

经典数值算例一维数值算例算例1(Lax激波管问题)算例特点：为左稀疏波，右激波类型，能够检验数值算法捕捉间断的能力；求解区域：x0,1;初始条件： ,v,p=0.445,0.698,3.528 当0x0.5, 0.5,0,0.571 当0.5x0.571,边界条件：左右两边采用紧支边界条件；计算时间：t=0.13.算例2(Sod激波管问题)求解区域：x0,1;初始条件： ,v,p= 1,0,1 当x0.5,0.125,0,0.1 当x0.5,边界条件：左右两边采用紧支边界条件；计算时间：t=0.25.算例3(Shu-Osher问题)算例特点：具有高震荡的性质，用来检验高阶数值格式比低阶数值格式对该类型的解具有更好的逼近效果；求解区域：x0,1;初始条件： ,v,p=3.85714, 2.629369, 10.333333 当0x0.1, 1.0+0.2sin50x-25, 0, 1 当0.1x1,边界条件：左边界采用紧支边界条件，右边界采用流入边界条件；计算时间：t=0.18.算例4(Woodward-Colella激波相互碰撞问题)算例特点：双激波的碰撞，其物理量过程更为复杂，检验数值算法捕捉间断的能力及稳定性；求解区域：x0,1;初始条件： ,v,p=1, 0, 1000 当0x0.1, 1, 0, 0.01 当0.1x0.9 1, 0, 100 当0.9x1,边界条件：左右两边采用反射边界条件；计算时间：t=0.038.二维数值算例算例1(二维Riemann问题)求解区域：x,y0,10,1;初始条件：uvpT= 1.5, 0, 0, 1.5 当x0.8,1, y0.8,1, 0.5323, 1.206, 0, 0.3 当x0,0.8, y0.8,1,0.138, 1.206, 1.206, 0.029 当x0,0.8, y0,0.8, 0.5323, 0, 1.206, 0.3 当x0.8,1, y0,0.8,边界条件：x,y方向均采用紧支边界条件；计算时间：t=0.8.算例2(激波通过气泡的问题)求解区域：x,y:0x325,-44.5y44.5;初始条件：uvp=1.3764, 0.394, 0, 1.5698, 当x225.0, 0.138, 0, 0, 1, 当(x-175)2+y2225.0,1, 0, 0, 1, 其他边界条件：左边界采用紧支边界条件，右边界采用流入边界条件，上下边界均采用反射边界条件；计算时间：t=2.0.算例3(Rayleigh-Taylor不稳定性问题)计算区域：x,y0,140,1;初始条件：uvp=2, 0, -0.025ccos8x, 2y-1, 当x,y0,140,12,1, 0, -0.025ccos8x, y+32, 当x,y0,1412,1,其中c=p/为声速，绝热指数=53；边界条件：左右边界均采用反射边界条件，上下边界均采用Dirichelet边界条件，其中下边界满足, u, v, p=(2, 0, 0, 1)，上边界满足, u, v, p=(1, 0, 0, 2.5)；计算时间：t=1.95.算例4(Double Mach Reflect problem不稳定性问题)问题描述：马赫数为10的强激波，入射与x轴成60度，计算区域：x,y0,40,1;初始条件：在x=16,y=0马赫数为10的激波；边界条件：左边界采用流入边界条件，右边界采用流出边界条件，上边界采用马赫数为10的激波真解，下边界采用exact post-shock条件(0x16);计算时间：t=0.2.算例5(内爆问题)求解区域：x,y:-0.3x0.3, 0.3y0.3;初始条件：uvp= 1, 0, 0, 1, 当x-0.15y-x+0.15,0.125, 0, 0, 0.14, 其他,边界条件：左右边界和下边界均采用反射边界条件，上边界采用流出边界条件；计算时间：t=0.8. 边界处理周期边界条件：密度，速度，压力按周期T取值；反射边界条件：密度，压力取边界对称值，压力取边界对称值的相反数；流入(流出)