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

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

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

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

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

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

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

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

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

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

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

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

超声速飞行器声爆飞行试验与声爆预测评估

声爆问题是超声速民机研制中首要解决的关键问题之一。声爆飞行试验是研究超声速飞行器声爆特性的最直接手段,可为声爆预测方法和低声爆设计技术提供真实可信的验证数据,对新一代低声爆超声速民机设计具有重要意义。中国航空研究院与中国飞行试验研究院发展了基于传感器阵列的地空一体化地面声爆测量技术,开展了国内首次超声速飞机声爆专项测试飞行试验,采集了多组真实大气条件下的声爆实测数据,验证了飞行试验方案的合理性。在飞行航迹正下方测得的声爆波形具有显著相关性,声爆信号头激波峰值相对误差在18%左右,尾激波峰值相对误差在8%左右,声爆持续时间均为0.1 s。对比分析了地面声爆实测数据与数值预测结果,发现：飞行器从测量阵列正上方飞过时,得到的声爆信号基本形态一致、持续时间较接近,声爆信号头激波、机翼前缘激波峰值相对误差小于5%;由于计算模型简化和声爆长距离传播的非线性累积效应等因素,导致声爆信号预测值与实测值在局部特征上有一定差异;后续还需深入研究真实大气环境下的超声速声爆远场传播预测方法。     

低声爆超声速客机声爆预测及不确定度量化分析

声爆精确预测问题是制约超声速客机技术突破的关键瓶颈之一。由于大气来流条件不断发生波动存在不确定性，为得到更可靠的声爆近/远场信号，需考虑来流参数的不确定性对声爆预测结果的影响，并甄别其中影响声爆预测的关键性因素，为工程实际应用提供有价值的参考。本文基于CFD声爆近场信号模拟和增广Burgers方程的声爆远场信号预测方法，对第三届声爆预测研讨会（SBPW-3）的C608低声爆超声速飞行器开展了声爆信号特性分析。首先，采用约5 021万非结构混合网格半模计算了近场过压值，并研究了远场地面波形计算的时间和空间网格收敛性。接着，分析了基准状态下复杂近场流动特征及地面波形特征，通过与C608飞行器公开数据对比，验证了该方法和自研程序的准确度。此外，研究了不同物理模型和大气相对湿度对地面波形的影响。在此基础上，运用基于非嵌入式多项式混沌（NIPC）方法开展了对不同来流参数（来流马赫数、来流攻角和单位雷诺数）的不确定度量化分析和敏感性分析。结果表明，在给定的输入变量不确定度条件下，地面波形波峰与波谷处过压值变化明显。相比而言，来流单位雷诺数对地面波形过压值的影响显著低于来流马赫数和来流攻角的影响。     

基于KZK方程的声爆频域预测法

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

超声速飞行器近场声爆信号反演技术

超声速民机是新一代民机的重要发展方向，其独有的声爆现象是制约其在陆地上空进行超声速飞行的最关键因素。对超声速飞行器的气动外形进行反设计是声爆抑制的有效途径。基于等效面积分布开展反设计，需要远场感知声压级作为直接指导。为此，提出了逆向传播分别与本征正交分解（proper orthogonal decomposition, POD）及伴随方程结合，根据远场声爆信号反演近场声爆信号的方法。采用第二届声爆会议（SBPW）提供的LM1021标模算例，并从远场频域内声压级、响度级和感觉噪声级进行反演可信度评估。结果表明，对于给定的任意远场声爆信号，基于逆向传播结果进行POD反演及伴随方程反演，都可以得到较为准确的近场过压信号，且伴随方程反演方法具有更优的高频信号即局部激波信号反演能力，远场感知声压级更精准。反演结果相应的等效面积分布与参考值高度吻合，表明此方法能够为等效面积指导的低声爆气动优化设计提供基础。     

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

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

兼顾气动和近场声爆特性的伴随优化

发展气动/声爆伴随优化设计方法对于研制新一代超声速客机具有重大意义。典型的气动/声爆伴随优化设计方法分为两种：近场信号的伴随优化和远场声爆预测耦合伴随设计优化。将最新发展的基于离散伴随理论的近场声爆信号反设计与气动特性优化结合,提出兼顾气动和近场声爆特性的伴随优化策略。首先,概述了伴随方程梯度求解和提出的基于信号射线方法的近场过压提取及目标过压装配过程,并验证了优化流程中近场信号及地面波形预测的准确性。其次,针对超声速民机标模,开展了低阻优化、近场低声爆反设计、兼顾气动和近场声爆特性的伴随优化三种策略的对比研究,结果表明兼顾气动和近场声爆特性优化策略的伴随优化在综合提升气动和低声爆特性上表现更优。最后与文献中相同算例的优化结果进行对比分析,进一步证明了该优化策略在超声速民机设计中的应用潜力。     

一种先进超声速民机低声爆高效气动布局设计

低声爆高效气动布局设计是超声速民机研究的重点和关键技术之一。采用基于声爆最小化理论反设计方法、波系有益干扰后体设计方法、参数化近场超压信号的混合可信度反设计方法,提出了一种先进超声速民机低声爆气动布局,对每一步降低声爆的效果进行了分析并研究了该布局的全声爆毯特性;采用CFD数值求解近场声爆信号并通过Burgers方程传播到远场,研究了飞行高度、飞行马赫数等参数对该气动布局地面声爆响度的影响;采用CFD数值模拟方法研究了飞行高度、马赫数等参数对该气动布局气动特性的影响。研究表明,采用基于声爆最小化理论反设计方法降低了基准气动布局的地面声爆响度约6.54 PLdB,采用波系有益干扰后体设计方法进一步使地面声爆响度降低了约0.97 PLdB,采用参数化近场超压信号的混合可信度反设计方法使气动布局地面声爆响度进一步降低了约4.04 PLdB;合理地设计飞行高度、飞行马赫数,可以有效地降低地面声爆响度;合理地选择巡航飞行高度和巡航飞行马赫数,可以有效地提高巡航效率。研究工作对超声速民机气动布局设计具有一定的工程指导价值,对超声速民机总体气动方案设计亦具有一定的工程借鉴意义。     

美国公布低声爆验证机方案

低声爆技术的成功与否，将是民用超声速飞行能否“重现江湖”的关键。飞机在超声速飞行时产生的强压力波传到地面时，会形成如同雷爆的声音，这种声音被称为声爆。声爆会使人们产生不适的感觉，也有可能会损害建筑物，因而美国法律禁止民用飞机在陆地上空进行超声速飞行。     

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

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

分层大气湍流场对远场声爆传播的影响

研究大气湍流对声爆的影响对于新一代超声速民机的低声爆设计具有重要意义。在真实大气环境中，大气湍流是普遍存在的，其会使低声爆波形在由高空向地面传播的过程中发生扭曲、畸变和振荡，并改变波形的能量分布。利用基于单向近似假设的二维HOWARD方程，以类Tu-144飞机的N型波和经过JSGD声爆最小化理论设计的低声爆波形为对象，研究分层大气湍流对N型波和低声爆波形的影响。首先，介绍了考虑大气边界层湍流效应的声爆传播模拟方法，包括改进的二维HOWARD方程及其离散求解方法、分层大气湍流场的建模方法，并采用JAXA抛体试验数据对预测方法的有效性进行了验证，同时与KZK方程预测结果进行了对比。其次，运用本文发展方法研究了N型波和低声爆波形在分层大气边界层湍流场中的传播情况，对地面可观测波形的主观噪声级和超压峰值进行了统计分析，发现大气湍流对增大N型波和低声爆波形主观噪声级（ASEL、PLdB）的概率很小，但很大程度上会增强波形超压峰值。最后，通过改变大气湍流生成过程中3个参数的数值，研究了不同强度的大气湍流场对声爆波形的影响，结果表明：风速脉动幅值和湍流积分尺度的增强对N型波和低声爆波形的感觉声压级影...     

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

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

基于逆向增广Burgers方程的声爆反演技术

声爆问题一直是制约超声速客机发展的障碍之一，因此低声爆设计技术在超声速客机的设计中尤为重要。逆向增广Burgers方程可将中场声压信号反演至近场，能够为低声爆反设计工作提供优化目标。围绕逆向增广Burgers方程展开研究，运用算子分裂法和正则化伪抛物型方程法对逆向方程进行数值求解，并求取对应的反向等效面积。借助3个标准算例，对数值方法的收敛性、求解准确度、不同中场高度及不同周向角下反演技术的计算精度进行了研究。针对逆向增广Burgers方程的特点，探讨了在中场设置目标波形进行反设计方法的可行性，并对声压信号中高频分量的传播特性进行了研究。研究认为上述方法能够较为准确地完成中场至近场的声压信号反演计算，并且使用中场目标波形代替地面目标波形能减轻反演距离对反演计算过程的影响。     

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

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

OS﹣X0试验飞行器声爆特性飞行测量与数值模拟分析

声爆问题是目前制约超声速民机研制的核心关键问题。但是,目前国内外关于声爆的飞行试验测试数据非常匮乏,严重制约了声爆预测与低声爆技术发展。基于零壹空间公司的"重庆两江之星"号OS-X0科学试验飞行器,航空工业气动院进行了声爆的飞行试验研究。采用航空工业气动院研发的飞行过程声爆信号地面测量技术成功获得了OS-X0试验飞行器在真实大气环境下的地面声爆信号,借助ARI_Boom声爆数值模拟平台对数据进行了分析和整理。数值计算结果与真实试验测量结果总体符合较好,但仍存在一定差异,这表明真实飞行环境下影响声爆特征的因素较多,在超声速飞行器声爆理论和预测方法等方面还需更深入的研究,同时飞行试验测量技术也有待进一步提高。