"""

普朗特-迈耶膨胀波数值计算方法应用

1. 问题背景

1.1 膨胀波流体力学基本概念

1.1.1 膨胀波的定义

• 膨胀波是流体中扰动区与未扰动区的分界面，流体通过此界面会压力降低。
• 普朗特-迈耶膨胀波是膨胀波的一种典型，即超音速气流绕凸角的定常流动。

1.1.2 膨胀波产生原理

• 定常直匀超音速气流绕直壁流动时，微小毛疵引起的扰动传播形成膨胀波。
• 扰动区与未扰动区的分界处产生压强差，非扰动区气流快速进入扰动区。

1.2 普朗特-迈耶膨胀波的特性

1.2.1 应用背景

• 普朗特-迈耶膨胀波广泛应用于空气动力学，喷气推进和高超声速气动力学等领域。
• 它作为膨胀波，是降低压力密度，增加速度的过程，对飞行器流场形状具有重要影响。

1.2.2 数学模型

• 流场分为上游区、扰动区和下游区。
• 利用微元法对微小扰动进行计算可得出其微分关系式，推导出dθ的表达式。
• 在量热完全气体条件下气体状态以及总温表达式。

2. 数值算法及程序

2.1 数值积分方法的选择

2.1.1 复合Simpson求积算法

• 复合Simpson求积算法适用于解决空气动力学问题。
• 它通过复合求积公式和递归迭代求解积分问题。

2.1.2 Romberg求积算法

• Romberg求积算法是数值积分方法的一种，适用于求解复杂积分问题。
• 它通过逐步增加精度来提高积分结果的准确性。

2.2 代码实现

2.2.1 复合Simpson求积算法源码

• 实现复合Simpson求积算法，完成人徐展。
• 输入Ma和精度，计算θ、迭代次数、误差和运行时间。

2.2.2 Romberg求积算法源码

• 实现Romberg求积算法，完成人徐展。
• 输入Ma和精度，计算θ、迭代次数、误差和运行时间。

3. 结果及分析

3.1 复合Simpson求积结果

3.1.1 计算结果

• 随着Ma的增加，气流扰流凸角θ增大，存在收敛趋势。
• 迭代次数和运行时间也随着Ma的增加而有不同程度的增大。
• 计算结果的误差始终在10-7以下，计算结果的精度比较可靠。

3.1.2 结果图表

• 图1展示了复合Simpson求积结果截图。
• 表1列出了复合Simpson求积结果表。

3.2 Romberg法求积结果

3.2.1 计算结果

• 随着Ma的增加，气流绕流凸角θ增大，计算结果的误差也随之增大。
• Romberg法误差比复合Simpson法误差较大。

3.2.2 结果图表

• 图2展示了Romberg求积结果截图。
• 表2列出了Romberg求积结果分析。

4. 结论

4.1 算法对比

• 两种算法的精确度都符合要求，但复合Simpson求积代码误差较小，准确度更高。
• Romberg求积法误差较大，但迭代次数少，运行时间相对较快。

4.2 实践意义

• 数值计算方法在航空航天工程中具有广泛应用和重要意义。
• 数值计算方法能够帮助工程师更准确地预测和控制飞行器流场。
• 通过数值计算方法，工程师可以优化飞行器设计，提高飞行器性能。

4.3 感谢

• 感谢陈青全教员的辛勤付出和耐心讲解。
• 希望课程越来越好，为航空航天工程培养更多优秀的人才。 """