基于Khokhlov-Zabolotskaya-Kuznetsov方程的声爆频域预测法

建立了一种快速预测声爆传播特性的频域方法,基于传统Khokhlov-Zabolotskaya-Kuznetsov(KZK)方程,描述声爆沿激波波阵面法线方向的传播.为了验证模型正确性,以NASA TD N-161中试件C为对象,首先基于不同网格,利用计算流体动力学(CFD)方法得到超声速流场;将CFD结果作傅里叶变换后代入模型方程,快速求解声爆传播特性.预测结果与NASA实验结果符合很好,研究表明:预测方法能够捕捉声爆的非线性传播,声衍射项使得声场在远场趋于轴对称分布,远距离传播后声能量集中于低频分量.      

基于广义Burgers方程的超声速客机远场声爆高精度预测方法

声爆高精度预测技术是新一代环保型超声速客机设计的核心关键技术之一。基于广义Burgers方程发展了可考虑"大气风"效应的远场高精度预测方法,开发了声爆预测程序"bBoom",并研究了近场声爆信号提取位置和"大气风"对远场声爆计算结果的影响。首先,给出了广义Burgers方程和声爆传播射线的计算方法,重点讨论了方程中计算分子弛豫效应和热黏吸收效应的关键参数。其次,通过简单轴对称构型标模算例、NASA的C25D构型和洛马的LM1021构型等复杂超声速客机算例对所发展的方法进行了验证,表明本文发展的方法在预测远场声爆时具有较高可信度。最后,基于所发展的方法对比了由不同近场提取位置传播到地面的波形,研究了飞机向不同方向飞行时,"大气风"对地面声爆强度和地面影响域的影响。结果表明:对于类C25D标模构型,为了确保远场声爆预测结果具有较高精度,应取机身下方约3倍机身长度位置处的近场压强信号作为传播方程的输入;另外,"大气风"会影响地面声爆强度及地面影响域,在预测时有必要加以考虑。     

飞机超声速机动飞行条件声爆预测方法

针对飞机超声速机动飞行条件下的超声爆问题，基于广义Tricomi方程建立了表征焦散线附近压力波强非线性演化过程的空间模型，采用伪时间步推进和算子分裂算法，结合时域和频域方法开发了可求解超声爆问题的计算程序ARI＿Superboom。构建的求解方法可分为3个步骤：首先，应用航空工业气动院自主研制的CFD软件ARI＿Overset在三维空间求解N-S方程，得到作为声爆远场传播初始值的近场空间压力分布；其次，基于自主研发的ARI＿Boom声爆预测程序开展射线追踪和声爆远场传播计算，获得焦散线及指定高度处声爆信号；最后，通过ARI＿Superboom预测程序获得焦聚区的地面超声爆信息。通过飞机超声速匀加速状态下地面产生的超声爆预测算例，验证了基于广义Tricomi方程建立的超声爆预测方法的可行性和准确性。     

基于流场/声爆耦合伴随方程的超声速公务机声爆优化

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

基于增广Burgers方程的超声速客机远场声爆预测研究

声爆问题是限制超声速客机投入使用的巨大阻碍,精确的声爆预测方法是解决超声速客机声爆问题的关键。声爆远场预测方法主要有波形参数法和增广 Burgers 方程法,相较于波形参数法,增广 Burgers 方程法有更完善的模型和更高的预测精度。本文首先对增广 Burgers 方程求解,提出一种有效的非线性效应求解方法;然后基于算子分裂法构建声爆远场预测方法,利用第二届国际声爆研讨会的标准算例对该方法进行验证;最后分析时间、空间网格密度和缓冲信号对预测结果的影响。结果表明：基于 Burgers 方程的声爆求解方法不需要时间、空间网格达到收敛所需的网格密度,而是按照实际精度需求选取适当网格密度;缓冲信号的引入对于大型超声速飞机的声爆预测更加准确。     

基于广义Burgers方程的声爆传播特性大气湍流影响

声爆是发展超声速民机不可回避的关键问题之一。目前流行的声爆预测技术主要针对静止大气，对大气中的湍流扰动效应考虑不足，尚未建立高效高逼真度的预测方法。基于广义Burgers方程的远场声爆预测方法，通过与射线法相结合，建立了一套可考虑热黏性吸收、分子弛豫等物理效应的大气湍流声爆影响快速预测方法，并采用该方法开展了大气湍流强度和大气边界层厚度对典型远程超声速民机的声爆特性影响规律研究。计算结果表明：建立的预测方法能够合理表征热黏性吸收、分子弛豫等大气物理效应，相比前期基于波形参数法框架的预测方法，能够更加真实地反映大气湍流对声爆传播特性的影响；相比于前期的典型超声速公务机，采用的远程超声速民机声爆波形更加复杂，该预测方法仍能给出复杂波系的大气湍流影响规律；随着湍流强度和边界层厚度的增加，大气湍流效应对声爆特性产生的随机性影响呈增强趋势；同时，声爆在地面到达点的位置也呈现出更加分散的趋势，其可能会改变声爆毯对地面的影响范围，应在飞行轨迹规划中予以考虑。     

广义Burgers方程声爆传播模型高阶格式离散

声爆预测技术是研制新一代环保型超声速民机的关键技术之一。针对当前基于广义Burgers方程的声爆远场传播模型的数值求解格式精度较低的情况，开展了高阶格式离散方法研究。通过分析该模型中各效应项的物理性质，应用合适的高阶精度格式对各项分别进行离散求解，并采用经典几何声学射线法计算声爆传播路径，实现了对地面声爆波形的精准预测。通过对美国F-5E机型的声爆飞行试验和第二届国际声爆预测研讨会的典型算例进行数值模拟，验证了本文高阶离散方法的可靠性。进一步的结果分析表明，采用高阶离散格式的预测方法相比传统二阶精度预测方法具有明显的低耗散特性，在同等网格数目上可获得更高分辨率的计算结果，且其网格收敛性远优于传统二阶精度格式，同时在计算效率上也有一定优势；在Burgers方程的各效应项中，非线性项的影响更为明显，因而采用高阶离散格式时其优势更为突出；热黏性吸收项对数值结果影响较小，采用高阶离散格式对计算精度的提升效果也不明显，实际计算中仍可以采用传统离散格式，甚至忽略该项的贡献。     

超声速条件下旋成体声爆典型影响因素分析

声爆问题是制约未来超声速客机进入真正商业运营的关键技术之一。航空工业气动院发展了基于近场CFD模拟和远场传播模型相结合的混合声爆数值预测平台ARI_Boom,基于典型旋成体模型。在Ma=1.41～6.0范围内,研究了飞行参数和外形参数对地面声爆特征(包括最大过压和特征正冲量)的影响。首先给出了ARI_Boom预测平台计算方法介绍,并通过典型算例对其进行验证;其次,基于所建立平台研究了飞行参数和外形参数对地面声爆特征影响。结果表明:针对所研究构型,减小半锥角、减小长细比、降低飞行马赫数和提高飞行高度均能减弱地面声爆特征;随着钝度系数的增加,最大过压先减小后增加,并存在最佳钝度。     

基于声扰动方程的气动噪声传播积分计算方法

传统声比拟思想假设声传播介质均匀,所获得的远场噪声积分方法无法考虑近场流动非均匀性引起的声折射效应,声场预测结果往往存在较大的数值误差。为解决上述问题,本文采用声扰动方程描述声波运动,并利用声学互易原理,在频域内构建关联近场声源与远场声压的伴随格林函数,从而将非均匀介质折射效应纳入伴随格林函数。基于伴随格林函数和声扰动方程源项的数值结果,利用积分方程计算气动噪声的空间传播。结合上述方法对二维圆柱和NACA0012翼型绕流噪声开展了数值预测研究,远场声压预测结果与CFD直接计算值吻合。数值研究表明,本文提出的基于声扰动方程的气动噪声传播积分计算方法可以考虑非均匀介质对声传播的影响,能够提高气动噪声的预测精度。     

均匀各向同性大气湍流对声爆传播特性的影响

基于离散Fourier模态有限和生成的随机大气湍流场，采用修正波形参数方法，开展了均匀各向同性大气湍流对典型超声速客机声爆传播特性的影响分析。计算采用的超声速客机模型为自行生成的简化超声速公务机模型。首先，应用航空工业空气动力研究院自主研制的CFD软件ARI_Overset在三维空间求解Navier-Stokes方程，得到作为声爆远场传播初始值的近场空间压力分布；其次，基于Von Karman能量谱，采用离散Fourier模态有限和形式生成随机均匀各向同性大气湍流场；最后，采用修正波形参数方法模拟了声爆信号在随机速度湍流场中的传播过程。数值结果表明：各向同性大气湍流对于地面声爆特征有重要影响。与无湍流状态相比，均匀各向同性大气湍流使得地面声爆特征增强的概率约为55%，使得地面声爆特征减弱的概率约为45%，故总体而言大气湍流效应更倾向于增强地面声爆特征；均匀各向同性大气湍流对于声爆传播路径影响相对较小，但是这一变化仍会导致地面信号接收点的不确定性。     

基于增广Burgers方程的音爆远场计算及应用

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

基于代理模型的高效全局低音爆优化设计方法

研究发展高效实用的低音爆优化设计方法,对于新一代低音爆超声速客机的研制具有重要的理论意义和应用价值。目前国内外发展的低音爆优化方法主要包括遗传算法（GA）和基于Adjoint的梯度优化。遗传算法虽然具有较强的全局优化能力,但其优化效率较低,无法很好满足实际应用的需要;而梯度优化虽然优化效率高,但易陷入局部最优。将最新发展的代理优化算法与音爆预测方法相结合,发展了一种具有全局优化能力的高效低音爆优化设计方法。首先,概述了所采用的线性音爆预测方法,并用NASA超声速圆锥体模型进行验证,表明其计算效率高、预测精度可满足飞行器初步设计的需要。其次,对所采用的代理优化（SBO）方法进行了概述,包括试验设计、代理模型建模、优化加点准则和收敛标准等。再次,运用所发展的方法开展了NASA多段圆锥体模型的低音爆优化设计算例研究,并与遗传算法和梯度优化的结果进行了比较,表明其优化效率比遗传算法提高了2个量级以上,且优化结果优于梯度方法。最后,将所发展的方法应用于AIAA音爆预测大会提供的翼身组合体外形（69°后掠三角翼）的低音爆优化设计,将远场音爆N型波峰值减少了27.4%,表明所发展的方法在复杂外形低音爆优化设计中具有很好的应用潜力。     

Effect of longitudinal lift distribution on sonic boom of a canard-wing-stabilator-body configuration

After the last flight of the Concorde in 2003, sonic boom has been one of the obstacles to the return of a supersonic transport aircraft to service. To reduce the sonic boom intensity to an acceptable level, it is of great significance to study the effect of lift distribution on far-field sonic boom, since lift is one of the most important contributors to an intense sonic boom. Existing studies on the longitudinal lift distribution used low-fidelity methods, such as Whitham theory, and in turn,o...     

声爆近场空间压力风洞测量技术

针对暂冲式超声速风洞中的声爆试验，发展了近场空间压力精确测量技术，以航空工业空气动力研究院的FL-60风洞为例，开展了技术验证。FL-60风洞是一座典型的亚跨超三声速下吹式风洞，其试验马赫数范围为0.3~4.2，试验段尺寸为1.2 m×1.2 m，单车次试验时间通常为数十秒。根据暂冲式风洞试验时间短、耗气量大等特点，设计了无反射测压轨以代替传统的静压探针，大幅提高了声爆近场空间压力的测量效率。通过CFD技术对无反射测压轨的流动特性、模型安装位置以及风洞试验段中的波系进行了分析，验证了测压轨设计方案的可行性。采用Seeb-ALR低声爆标模和自行设计的带喷流的旋成体模型进行了验证性试验，采用参考车次方法和空间平均技术获得了高质量的数据，试验测量结果与CFD计算结果一致性较好，验证了声爆近场空间压力测量系统设计的合理性。     

宽频声传播计算中多自由度时域阻抗边界条件的适定性

旨在验证一种多自由度时域阻抗边界条件在宽频声传播问题中的适用性与精度.选用NASA(美国国家航空航天局)流管实验数据进行数值校核;采用高精度计算气动声学方法数值求解线化欧拉方程.计算结果与实验数据及单频结果均吻合很好,表明了该多自由度时域阻抗边界条件具备精确处理宽频声传播问题的能力.      

超音速公务机声爆计算与布局讨论

超音速公务机是航空工业的重要发展方向之一,低声爆设计技术是超音速公务机的关键技术,但国内在该领域几乎没有任何研究基础,无法为超音速公务机提供足够的设计支持。为此,介绍了两种计算声爆的方法：一是从超音速流动的线化方程出发,推导体积和升力产生声爆强度的估算方法,该方法适用于飞机概念设计阶段;二是从非线性声学传播方程出发,使用CFD近场结果作为输入,编程计算声爆强度,该方法适用于飞机初步/详细设计阶段。在此基础上,对影响声爆强度的参数进行初步分析,结果表明：飞机重量和飞行高度对声爆强度影响很大,展弦比、翼载等参数对声爆强度的影响较小;＂细长机身＋鸭式布局＋大后掠三角翼＂布局比较有利于减小声爆强度。     

典型标模音爆的数值预测与分析

精确预测音爆对超声速民机的研制具有重要意义。主流的音爆强度预测方法分为两步,首先通过风洞试验或CFD方法得到近场音爆过压（Over-pressure）分布,再运用修正线化理论或非线性声学理论将近场过压传播至地面,最终获得地面音爆的声压信号。本文运用典型标模对当前音爆数值预测方法的精度进行了验证和确认。在近场音爆过压分布的数值预测方面,分别考察了超声速尖点构型前缘修形尺度、不同空间离散格式和无黏/有黏流动控制方程求解对近场过压计算结果的影响。远场音爆预测方面,以LM1021全机构型近场过压分布为输入,使用基于波形参数法的远场传播工具分别考察了不同离散格式和有/无黏性计算的近场过压分布差异对地面音爆结果的影响。算例结果表明,尖点构型近场音爆预测中进行几何修形是十分必要的,使用相对合理的过渡球半径可以保证近场音爆预测精度,过大的修形尺度会对激波形状、激波和膨胀波的峰值均产生较明显的影响;就近场波形而言,熵相容格式计算得到的结果与试验测量值吻合最好,但不同离散格式导致的近场预测波形差异对传播到远场的波形关键指标（主要是最大过压和上升时间）的影响很小;是否计入黏性对近场波形结果尽管仅有小幅的影响,但将近场信号传播到远场得到地面波形时,这些细微差异会在远场波形的音爆评价关键指标上表现出明显的区别。