有限差分法求解Sod激波管问题综述.pptVIP

* 中国石油大学（北京）化工学院 有限差分法求解Sod激波管问题 姓名： 陈晓慧 专业：化工过程机械 2016.05.06 一、 问题描述 Sod激波管问题属于黎曼问题，如图所示左侧为高温高压气体，右侧为低温低压气体，中间用薄膜隔开。t=0时刻，突然撤去薄膜，分析激波管中气体的运动。 初始条件（t=0）： 输出时刻：无量纲时间t=2.0 二、控制方程 一维Euler方程（守恒形式）： 方程拟线性化： 其中， 无量纲压力 无量纲定容热容 三、方程离散 1.Jacibian矩阵特征值的分裂 雅克比系数矩阵A（U）的三个特征值分别对应为λ1=u；λ2=u-c；λ3=u+c，依据如下算法将其分裂成正负特征值： 2.流通矢量分裂（Steger-Warming分裂法） 3.一阶显示格式 其中， 无量纲声速 四、matlab编程 一阶迎风显示格式： 格式稳定性(CFL数)： 其中， 五、结果展示