摘要：音爆的精确模拟对于超声速客机的低音爆研究与设计意义重大。由于计算能力的限制,客机巡航高度的音爆全场直接模拟目前还难以实现。现有的音爆预测方法一般分为两步,先通过超声速线化理论或计算流体力学的方法得到音爆近场的过压值(Over-pressure)分布,再通过声学理论将近场声压信号推进至远场,以获得飞行器的地面音爆信号。在远场计算中,传统的波形参数法没有考虑音爆传播过程中的经典吸收和分子驰豫效应所造成的声能损失,得到的激波没有厚度,导致计算得到的远场声压级不准确。基于算子分裂法,开展了非线性声学中的增广Burgers方程的数值解法研究。通过计算第二届音爆预测研讨会(SBPW-2)发布的两个标准算例,验证了该方法可以实现地面音爆波形的精确预测。发现在近场声压信号前加入一段无幅值的缓冲信号可以有效提升"N"波上升时间的模拟精度。网格收敛性研究表明适当加密计算网格是有必要的。在此基础上研究了大气声吸收对地面波形的影响,发现分子驰豫效应的影响要强于经典吸收。最后,研究了不同湿度、温度对地面音爆波形的影响,发现干燥、低温的环境对音爆信号的过压值有抑制作用。

摘要：基于自主研发的大规模并行结构化网格CFD求解器PMB3D以及并行化伴随方程求解器PADJ3D,开展了流场/声爆伴随方程的求解研究。首先采用标准算例,对内部CFD代码PMB3D软件和声爆预测代码进行了声爆计算可信度验证,以及声爆强度对近场声压梯度的校核。针对并行环境下多块对接网格的近场声压提取操作的复杂性,提出了"包围盒"的方法实现并行环境下近场声压装配单元编号、网格块编号以及对应的进程编号确定,基于声爆计算坐标将并行传递的数据进行一维排序,为声爆预测、伴随方程以及梯度求解提供输入条件。通过线性插值雅克比矩阵实现均匀坐标系梯度信息向非均匀坐标转换,并进一步根据结构化网格特征提出了插值原则,简化了近场声压转换雅克比矩阵的变分。通过装配单元记录,实现声爆强度对流场守恒变量的变分结果向各个进程装配,将装配结果作为流场伴随方程的右端项实现流场声爆耦合伴随方程的求解。此外,对小型超声速公务机开展了声爆优化,对比分析了设计前后的声压及其频谱特性。