双马赫反射问题作业.docVIP

计算流体力学作业 ——双马赫反射问题数值解 参考文献： Woodward et al. Journal of Computational Physics, 54,115-173,1984 图1 双马赫问题的初始流动示意图 问题描述： 如图所示，假设初始时刻，Mach 10的正激波向平板运动。激波与平板之间成60°夹角。 试计算流场的演化，画出t=0.2时刻流场的密度及压力分布。 无量纲化： 以激波前的密度、温度及音速无量纲化。 要求： 计算域：0 x 4, 0 y 1； 采用均匀计算网格。 分别采用 480*120 及 960*240 网格数的网格计算。 采用无粘的Euler方程计算。 初值： Ma=10, 与壁面呈60°的正激波。提示： 采用正激波关系式给出 图2 双马赫反射问题的初值及边界条件示意图 如图2所示， 虚线表示激波，左侧部分为波后值?2,u2,v2,T2,p2 (密度，速度，温度，压力)； 虚线右侧为波前值??,u1,v1,T1,p1。 【提示】 波前的绝对速度为0， 无量纲密度和温度均为 1。将坐标系固定在激波上，运用正激波关系式计算波后速度。然后将速度转化到绝对坐标系下。 边界条件： 图中红色部分给定为波后值?2