我为我的作业编写了这些代码，如果您有现成的代码可用于寻找斜激波的声波、冯诺依曼和分离角等，那么解决冲击波动力学问题会更方便。 您还可以使用shockpolar.m 代码轻松找到滑移线角度。

运用Chemkin程序和GRI-Mech 3.0机理,对入射激波诱导下矿井乏风瓦斯点火延迟时间进行数值模拟.首先,根据乏风瓦斯燃烧基元反应敏感性分析,定义CH3峰值出现的时刻为乏风瓦斯的点火延迟时间,然后分别研究瓦斯体积分数、入射激波速度、乏风瓦斯初始温度、初始压力对点火延迟时间的影响.研究结果表明,入射激波速度、乏风瓦斯初始温度、初始压力的增加均会使燃烧温度增加,CH3峰值增加,点火延迟时间缩短,其中入射激波速度的增加缩短点火延迟时间的效果最为显著,乏风瓦斯体积分数变化(0.1%～1.0%)对点火延迟时间的影响较小.

我的程序允许用户选择他们想要解决的变量。 使用在 MATLAB 6.1 中创建的符号工具箱 他们被要求提供另外两个。 这是斜激波的 Beta-Theta-Mach 关系。 注意：求解 theta 将得到精确的答案 而求解马赫数或 Beta 则依赖于数值近似。 我输出近似值的误差 使用的算法：最基本的即 O(N^2)

使用TVD格式进行一维磁激波管或者激波管的数值模拟

8.6 运动正激波在固壁面的反射 有时会遇到静止气体中有固壁边界，运动激波在静止气体中传播后遇到固壁发生反射的 情况。当激波波阵面到达壁面的瞬间，波后气体会突然受到固壁的压缩。实际上，若把参考 系取在波后气体上，波阵面到达固壁时，相当于壁面突然以速度 向左压缩气体。在绝对 坐标系中，波后气体的伴随速度在固壁处瞬间变为零。因此从壁面开始将出现一道正激波， 即反射激波向左在原波后气体中传播。原波后气体现在成为反射激波的波前气体。反射波所 到之处，波后气体速度 ，状态为 ， 。 若已知初始激波波前参数， ， ， 及 ，欲求激波反射后壁面压力 （即 反射波波后压力）及反射激波的运动速度 （向左）。反射激波波前速度不为 0，为此，把 参考系固结在反射激波波前的气体上。在此参考系内反射波波前气体速度为零，激波速度为 （向左），波后速度为 。因此有 2 2 2 2 2 2 2 2 1 1 1 4 2 2 2 N v v v c                 1 1 343.9m sc RT  3 1 1.1991kg m  1 1 15240 44.31 343.9 N M c    2 1 2 2 1 12693.5m s 1 cN v N              2 5 2 1 1 2 1 2320 10 Pa 1 1 N p p c                     31 2 12 1 1 5.985 7.176kg m 1 2 1 1 c N                    2 2 2 112632K p T R   2 2 6728m sc RT  2 2 2 12693.5 1.887 6728 v M c    1 2 5 20 2 2 2 1 1 6.57 15236 10 Pa 2 p p M p              2 20 2 2 2 1 1 1.712 192843K 2 T T M T            2 v 3 0v  3 p 3  1 0v  1 p 1  1 N 3 p 2 N 2 2 N v 2 v 流 体 力 学 （ 国 科 大 教 材 ） 第 八 章

Euler方程组为描述无粘性可压缩流，本程序即为求解sod问题的程序。

对一维Sod激波管问题的分析和数值求解，内附python代码及注释。本算例中初值已经给定，首先通过理论推导和代入数值分析确定了算例中属于左膨胀波右侧激波情形，而后根据已知信息进行编程，用二分法求解未知压强，进而计算所有参数并绘制图像。

二维斜激波反射课件，计算流体力学专用~~~~~~~~~~~~~~

采用L-F分裂，GVC格式对一维激波管求解

使用二阶tvd格式对一维磁激波管/激波管进行数值仿真(matlab)，其良好的抑制非物理震荡优势比roe格式较好，流体激波管只需将磁场值全部设为零。适合磁流体仿真和流体仿真入门