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

    超声速飞行器近场声爆信号反演技术    马 创，黄江涛，*，刘 刚，陈 宪，舒博文，，陈其盛，高正红(1.中国空气动力研究与发展中心，绵阳 621000；2.西北工业大学 航空学院，西安 710072)0 引 言声爆是超声速飞机飞行中特有的气动声学现象[1-2]飞行器在超声速飞行时，其近场产生复杂的激波-膨胀波波系，这些波系传递到地面，形成空间内N形的声压分布[3]，这就是声爆现象声爆现象会对生物体和建筑物产生不良影响，因此早期的超声速客机只能在海洋上空进行超声速飞行20世纪中后期，英法和苏联研制的以“协和”号和“图-144”为代表的第一代超声速客机都因声爆过强而被许多国家禁止在境内飞行[4]，严重影响了其商业运营，最终以失败告终声爆评估与抑制是超声速民机发展必须解决的卡脖子问题[5]声爆抑制涉及空气动力学与声学等多个学科，目前主要途径有流动控制（静音锥[5]、热流注入[6]、矢量发动机[7]等）、新概念低声爆布局（双向飞翼[8-9]布局、可变前掠翼布局[10]、双翼布局[11]等）、气动外形优化（代理优化[12-14]、梯度优化[15-17]、进化算法[18]等）……研究表明，在各类方法中，飞机气动外形设计是声爆抑制的最有效途径之一。

其中，通过等效面积分布（equivalent area, EA）指导气动外形优化[19]是开展低声爆气动优化设计的一项关键技术等效面积分布是沿机身轴线的体积截面积分布与升力分布的叠加，直接决定了超声速飞行器的声爆特性，而飞行器等效面积分布特征较大程度由近场过压决定[20]但是，通过近场过压反设计抑制远场声爆，缺乏直接的远场感知声压级作为指导因此，探索有效的远场声爆信号向近场过压反演方法是开展声爆抑制研究的关键环节传统的低声爆反设计主要基于JSGD（Jones-Seebass-George-Darden）方法开展20世纪70年代末，Jones、Seebass、George、Darden 四人[21-26]首次提出基于超声速线化理论的JSGD方法，为低声爆反设计奠定了理论基础该方法在NASA的高速民用运输计划（high speed civil transport, HSCT）中得到广泛应用[27]国内冯晓强等使用声爆预测算法与遗传算法对JSGD参数进行了优化[28]但JSGD方法存在试凑的问题，且基于线化理论，对于复杂构型的计算误差较大，属于反向等效面积分布设计方法[29]因此，有必要探索一种正向等效面积分布设计方法。

本文拟对近场声爆信号反演方法的“逆向传播+POD（proper orthogonal decomposition）”和“逆向传播+伴随方程”两种实现途径开展研究，以期为等效面积分布指导的超声速低声爆飞行器气动优化设计提供技术支持1 声爆信号的逆向传播以零为初始值开始近场信号反演，无论是POD反演方法还是伴随方程反演方法，都将导致庞大的计算量，尤其对于POD方法，所需的样本空间极大因此获取可靠的初始波形是进一步准确反演的关键远场声爆信号的逆向传播可以得到近场声爆信号大致波形，为进一步细节反演提供有力技术支撑1.1 声信号逆向传播模型本文通过逆向增广Burgers方程来建立声爆信号的逆向传播模型推导过程如下经典增广Burgers方程的无量纲形式为：等号右端项依次对应声信号在大气中传播的非线性效应、经典耗散、非均匀介质和分子弛豫现象式中：τ为无量纲时间，参考时间为1 / ω0，由近场声爆信号的采样频率决定；Γ为无量纲气体耗散系数； θv为无量纲分子松弛时间；Cv为无量纲松弛系数在经典增广Burgers方程中，几何扩散项和大气分层项取决于声管传播路径与面积，与时间相关的项只有非线性效应项、经典耗散项和分子弛豫项。

建立逆向Burgers方程即对这几项进行修改，令时间反向流动[30]描述声信号逆向传播的逆向增广Burgers方程为：1.2 逆向增广Burgers方程求解逆向增广Burgers方程数值求解时，对过压信号进行时间离散，对声管传播路径进行空间离散，并对相关参量进行无量纲化处理使用算子分裂法[31]对式（2）进行求解，可分裂成以下五项：求解过程中，在当前高度依次更新无量纲压力P，空间推进直至目标高度，完成求解[32]式（3）可整理为如下形式：式（9）～式（11）给出了式（8）的离散方法：式（4）经典耗散项和式（5）分子弛豫项均为病态方程为稳定求解，引入伪抛物线型正则化方法进行修正，修正形式如下：依照工程经验选择 ε =1×10-4引入式（12）所示的伪抛物型方程改写经典耗散项和分子弛豫项经典耗散项可以写成：分子弛豫项可以写成：其中，α取0.5改写后的经典耗散项和分子弛豫项均为三对角方程使用TDMA方法[33]进行求解式（6）和式（7）为守恒方程，求解方法如下：2 基于POD理论的声爆信号细节反演由于正向传播与逆向传播均存在一定耗散，导致逆向传播结果中局部激波信号丢失，需要进一步对声爆信号进行精细化反演计算，以尽可能还原真实波形的细节。

本征正交分解是一种基于快照来描述高维空间的方法，其衍生的Gappy POD方法可进一步修复此高维空间内的任一残缺快照，因而可以用于反演近场声爆信号细节将设计好的目标远场波形作为残缺快照中的已知部分，待求的近场声爆信号作为残缺快照中的未知部分，通过对残缺快照的修复来实现对远场声爆信号的反演2.1 本征正交分解（POD）方法POD本质上是一种降维方法对于近远场样本组成的快照集 {Yi，i=(1,2,···,N)}(Yi∈Rn)，POD可以从中寻找主要模态作为基模态使所有快照在各阶基模态组成的子空间上正交投影最大，表达式为：然后，计算快照的脉动矩阵P：采用奇异值分解（singular value decomposition, SVD）方法进行特征值分解奇异值分解计算效率高于传统特征值分解方法，且在计算高阶模态时更精确对矩阵P进行奇异值分解可得：其中：U∈RN×N，V∈RN×N Λ ∈RN×n且仅有对角线元素Λii= σi(σ1≥ σ2≥ ···≥ σr≥0)非零，并得到各阶模态：任意一个快照（如第i个）都可以由各阶模态的线性叠加得到，定义各模态的能量为其对应的奇异值平方，按照能量原则来选取基模态，定义能量下限为99.9%，即满足式（27）的前p阶模态被选择作为基模态。

2.2 Gappy POD方法对于一个残缺样本，在POD的基础上发展出的Gappy POD方法可对其进行修复对于一个如下形式的残缺样本其中：K为样本中的已知元素，在本文中为远场声爆信号值；U为样本中的未知元素，在本文中为待求的相应近场声爆信号值POD方法提取出前p阶模态作为主模态，按以下形式给出：基于最小二乘法，使用主模态对残缺样本中的残缺部分进行修复：其中， Γ =(γ1,γ2···γL)T为残缺样本I在主模态张成的新的p维空间中的坐标修复得到的残缺样本中的缺失数据为：2.3 声爆信号的“逆向传播+POD”反演方法前文对POD理论的介绍也表明POD方法对样本空间的依赖性由于超声速飞行器近场激波-膨胀波波系的存在，近场声爆信号波形大多呈现“N”形起伏直接使用POD方法进行反演，对主模态的要求过于严苛，难以得到理想结果本文提出了逆向传播与POD结合的反演方法首先进行声爆信号的逆向传播，得到近场声爆信号的波形主特征以此为初始值进行POD计算，大大减少了POD的计算成本，仅需构造波形细节所需的样本空间，再对逆向传播的近场波形进行扰动取样即可采用“扰动取样—POD”迭代计算的方法进行反演，尽可能准确地还原声爆信号的局部激波。

3 基于伴随方程的声爆信号细节反演基于声爆伴随方程求解梯度信息与差分方法相比，避免了大量重复进行Burgers方程求解的正计算过程，是一种高效的梯度求解方法，可基于逆向传播结果反演声爆信号细节将目标远场声爆信号视为目标函数，近场过压信号的每个时间离散点对应的过压值作为设计变量，以逆向传播得到的近场声爆信号为初始值进行梯度寻优，以得到更精细化的反演结果3.1 声爆伴随方程推导与梯度求解将增广Burgers方程式（1）进行算子分裂后写成矩阵形式：式（32）～式（35）分别对应氧气分子（O2）弛豫效应、氮气分子（N2）弛豫效应、经典耗散、非线性扭曲这四个声爆信号传播的物理环节kn为增广Burgers方程式（1）右端第四项与第五项的乘积，表征几何扩散和大气分层对声爆信号幅值的影响P、q、r、t分别表示求解过程中声压的四个中间值由上述矩阵形式，构造式（36）所示的拉格朗日量，推导声爆伴随方程的矩阵形式其中：l为目标函数；N为垂直空间网格数；D表示设计变量，即每个时间离散点对应的过压值； γ0、 γ1、β、λ分别为式（32）～式（35）表示的四个物理环节对应的伴随变量式（36）使用链式求导法则对设计变量D求导，可得：令所有状态变量对设计变量的全导数为0，可以得到声爆伴随方程。

式（39）为声爆伴随方程，可进行迭代求解将迭代结果代入式（38），可求得目标远场声爆信号对当前近场声爆信号设计变量的梯度信息3.2 声爆信号的“逆向传播+伴随方程”反演方法绝大多数气动设计问题都可抽象为优化问题本文的对声爆信号反演的研究，可视为如下所示的单变量最优化问题：其中：Ftarget为目标远场过压信号；Ni为当前近场过压信号，以逆向传播所得的过压信号为初始值；FNi为 相应的远场过压信号；EM为欧氏距离；Ii为优化目标函数，只与当前一轮正计算的远场声爆信号有关通过求解声爆伴随方程更新目标函数对设计变量的梯度信息，基于序列二次规划算法进行梯度寻优，可得近场过压信号的反演结果4 基于LM1021标模的声爆信号反演技术验证采用第二届声爆会议SBPW（the second sonic boom prediction workshop）提供的LM1021标模算例[34]对提出的方法进行可信度验证LM1021外形示意于图14.1 近场声爆信号与远场预测在飞行器流场中，机身周向角0°的近场区域提取压力值，并计算得到过压值近场过压信号在空间内的分布如图2所示图2 近场声爆信号参考值Fig.2 Reference of near-field sonic boom signal图3给出了声爆预测程序得到的远场声爆信号预测结果。

由于预测程序中进行了时空变换，远场声爆信号为地面某一空间位置固定的观测点处的压力值随时间的变化图3 远场声爆信号参考值Fig.3 Reference of far-field sonic boom signal本文提出的反演方法是对图3所示的远场声爆信号进行反演计算，使反演结果具有与图2所示的近场声爆信号等效的可信度4.2 远场声爆信号的逆向传播基于逆向增广Burgers方程求解波形B的逆向传播结果，目标高度为飞行高度H= 16 747 m，声爆信号时间离散点数 ω0=50000为了减少离散格式的数值耗散，垂直方向网格数为 5 0000正则化参数ξ=1×10-6图4给出了逆向传播结果与真实近场过压分布的对比经过逆向传播得到的波形，较为准确地反映了真实近场过压分布的主要特征从图5中可以看出，仅经过逆向传播计算，近场声爆信号仍然存在较大误差，精度未达到设计要求逆向传播误差较大部分主要由病态方程正则化造成，需要进一步对声爆反演信号进行精细化计算，以还原被耗散掉的细节图4 逆向传播结果Fig.4 Result of reverse propagation图5 逆向传播结果相应的远场声爆信号Fig.5 Far-fiel。