2. 西北工业大学 航空学院, 西安 710072
2. School of Astronautics, Northwestern Polytechnical University, Xi'an 710072, China
声爆是下一代超声速民机必须解决的卡脖子问题,从目前开展的研究工作来看,在众多的声爆抑制技术研究中,气动外形优化是超声速民用飞机声爆抑制最有效的手段,因此,开展超声速低声爆气动外形综合设计技术研究具有重要意义。
国外在气动/声爆优化方面的起步较早,开展了较为系统性的研究,大部分研究工作基于伴随方程优化方法,基于伴随方程的优化进而也分为两个方向:近场声压伴随[1]与流场/声爆伴随方程,由于伴随系统具有计算代价小,梯度计算量与各个学科设计变量个数均无关等优点,在气动/声爆综合优化领域具有不可替代的优势,受到低声爆设计研究人员的积极关注。国内在2018年以来,工业部门、研究机构、高校开始普遍重视这方面的研究:2018年8月18日,首届全国空气动力学大会(绵阳)设置了“超声速飞行器声爆数值模拟与低声爆布局设计技术”专题研讨会;2018年11月19日,首届全国性的“超声速商用飞机技术研讨会”在北京召开,来自全国近19家单位参会;2018年12月5日《空气动力学学报》组织出版了首个声爆专栏;2019年“绿色超声速民机设计技术”入选中国科协发布的重大科学问题和工程技术难题之中,核心攻关内容就是声爆预测与抑制等,标志着国内在超声速声爆预测与抑制方面开始了较为系统化的研究。
在超声速民机优化设计方面,西北工业大学[4-10]、中国航空研究院[11]、中国航天空气动力技术研究院[12]、航空工业气动院[13]、北京航空航天大学[14-15]、中国空气动力研究与发展中心[16]等单位在声爆预测、代理优化、布局设计以及耦合伴随优化等关键技术方面开展了系列研究,取得了一定的进展。
超声速民机低声爆设计中,在飞行高度、飞行速度、重量这三类重要影响因素确定的前提下,气动外形综合设计则成为超声速飞行器减阻/低声爆抑制的核心工作。
本文研究工作主要分两个主要层次:第一,飞行器的气动布局形式选型设计;第二,气动外形精细化综合设计。两个层次的工作对设计技术的需求不同。本文基于笛卡尔求解器以及大型并行数值优化软件AMDEsign[17],开展了低声爆气动外形综合设计技术研究,系统总结典型的布局设计参数对近场过压分布的影响,对比分析详细设计前后气动外形细节变化以及近场过压形态,为设计人员提供技术参考。
1 声爆信号非线性预测文中基于Burgers方程进行声爆预测[18],基于该方程发展了声爆预测代码,包含非线性、吸收影响、分子松弛、射线管扩张以及大气分层等项,可以模拟声爆在分层有损耗的大气中传播:
$ \begin{array}{l} \frac{{\partial P}}{{{\rm{ }}\partial \sigma }}{\rm{ }} = {\rm{ }}P{\rm{ }}\frac{{\partial P}}{{\partial \tau }}{\rm{ }} + {\rm{ }}\frac{1}{{{\rm{ }}\mathit{\Gamma } }}{\rm{ }}\frac{{{\partial ^2}P}}{{\partial {\tau ^2}}}{\rm{ }} + \sum\nolimits_v {\frac{{{C_v}\frac{{{\partial ^2}}}{{\partial {\tau ^2}}}}}{{1 + {\theta _v}\frac{\partial }{{\partial \tau }}{\rm{ }}}}{\rm{ }}P - } \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{1}{{2G}}{\rm{ }}\frac{{\partial G}}{{\partial \sigma }}{\rm{ }}P + \frac{{{\rm{ }}1}}{{2{\rho _0}{c_0}}}{\rm{ }}\frac{{\partial {\rho _0}{c_0}}}{{\partial \sigma }}{\rm{ }}P \end{array} $ | (1) |
式(1)可以采用通过算子分裂法进行快速求解,各项的具体含义参考文献[18]。图 1为LM1021标准模型示意图,图 2给出了周向角为0°时的声爆计算结果与NASA代码sBOOM的对比,可以看出,文中采用的非线性声爆预测技术精度较高,能够为选型以及详细设计提供可靠平台。
2 流场声爆耦合伴随优化方法
流场、声爆耦合伴随方程是实现地面声爆信号反设计的关键,可以克服近场变分无法进行地面信号设计的缺点,因此文中采用流场/声爆耦合伴随方程开展综合优化。其表达式为:
$ \begin{array}{*{20}{l}} {\frac{{\partial {I_b}}}{{\partial {\rm{ }}\mathit{\boldsymbol{p}}{_0}}} + {\rm{ }}\mathit{\boldsymbol{\lambda }}_b^{\rm{T}} = 0}\\ {\mathit{\boldsymbol{\lambda }}_f^{\rm{T}}\frac{{\partial \mathit{\boldsymbol{R}}}}{{\partial \mathit{\boldsymbol{W}}}}{\rm{ }} - \mathit{\boldsymbol{\lambda }}_b^{\rm{T}}\frac{{\partial \mathit{\boldsymbol{T}}}}{{\partial \mathit{\boldsymbol{W}}}} = 0} \end{array} $ | (2) |
该方程组仍然由声爆伴随方程和流场伴随方程构成,通过声爆伴随变量实现方程组耦合。因此,核心仍然是构造声爆伴随方程和流场伴随方程,耦合伴随方程的各项表达式含义以及详细推导可以参考文献[16],这里不再赘述。图 3给出了声爆伴随方程的精度校核,其伴随变量的物理含义是远场声爆信号目标函数对近场声压的梯度,与有限差分结果的平均误差低于0.5%,完全满足设计需要。
3 低声爆气动布局设计探讨布局选型基于笛卡尔网格无黏求解器进行,笛卡尔网格见图 4。需要指出的是,由于采用无黏计算,该部分并不是用来进行远场声爆计算,而是进行近场波系结构定性分析,进行布局参数研究;基于影响声爆强度的几个主要因素,设计选型基本原则是增加激波上升时间、减小波系过压峰值,尽可能实现弱激波设计,对机身轴线弯曲度、机翼平面形状、翼面布置等影响进行了系统研究,选取一种近场过压分布符合上述特征的布局,开展进一步综合优化。
![]() |
图 4 笛卡尔网格 Fig.4 Cartesian grid |
通过对机翼前缘后掠角进行分段描述,实现机翼平面参数化建模。图 5给出了两种机翼方案,图 6给出了不同机翼方案条件下的波系分布与近场过压对比,可以看出连续变后掠前缘在提高翼身结合处激波上升时间的同时,对激波强度有明显抑制作用,同时机翼后缘位置带来相应膨胀程度的弱化,对尾激波强度起到积极的弱化作用。
![]() |
图 5 两种机翼方案 Fig.5 Two different wings |
![]() |
图 6 机翼方案A和B的波系分布与近场过压对比 Fig.6 Comparison of wave distribution, near-field overpressure for different wing plane shapes |
该小节主要研究了有无鸭翼以及鸭翼的位置、机翼位置对波系分布的影响。鸭翼方案A相对于鸭翼方案B,前移1.6 m;机翼位置A相对于机翼位置B前移2.5 m。图 7~图 9给出了不同机翼方案条件下的波系分布与近场过压对比,从过压分布上看鸭翼后移导致第一道激波后的膨胀区有所弱化,从而对翼身结合处的激波起到一定的抑制作用;机翼后移对尾激波起到一定的减缓作用。
![]() |
图 7 不同鸭翼方案的波系分布对比 Fig.7 Comparison of wave distribution for different canards |
![]() |
图 8 鸭翼不同位置过压对比 Fig.8 Comparison of overpressure for different >positions of canard wings |
![]() |
图 9 机翼不同位置过压对比 Fig.9 Comparison of overpressure for different positions of wing |
对机身横截面椭圆建模,通过改变z向位置、椭圆纵轴长度实现对机身轴线弯曲度的描述,图 10给出了轴线零弯曲机身方案A与轴线弯曲机身方案B的波系分布以及近场过压分布,从波系分布上看,弯曲轴线机身近场第一道激波更弱,且利于弱激波系形成,弱激波传播中迅速耗散,且激波之间的合并能够得到抑制;从过压对比上看,弯曲轴线机身同时使得头激波强度有所减缓。
![]() |
图 10 轴线零弯曲机身方案A与轴线弯曲机身方案B的波系分布以及近场过压分布对比 Fig.10 Comparison of wave distributions, near-field overpressure distribution for the fuselage axis with or without bending |
通过对机身横截面采用椭圆建模,改变纵向面积分布,进行近场过压分布的对比,图 11给出了机身面积分布方案A与方案C的近场波系与过压分布对比,其中方案A最大截面积靠前,方案C最大截面积后移。从波系分布上看,最大横截面积位置后移的C方案头部激波基本消除,整个过压幅值明显降低,且激波上升时间增加。
![]() |
图 11 机身面积分布方案A与方案C的近场波系与过压分布对比 Fig.11 Comparison of near-field wave system distribution, overvoltage distribution for different airframe area distributions |
通过上述特征参数影响研究,提出一种连续变后掠机/尾翼、弯曲轴线机身、翼上动力超声速民机布局,用于低声爆技术研究。设计状态:Ma=1.5,H=16.7 km,重量41 t,如图 12所示。由于缺少发动机工作参数,因此,在全机气动优化过程中采用通气短舱构型。
![]() |
图 12 超声速民机气动布局 Fig.12 Supersonic civil aircraft aerodynamic configuration |
采用多区域自由式变形技术(FFD)对全机参数化建模,设计变量控制顶点如图 13所示。机身控制顶点为最上和最下方的节点,FFD技术的基函数为三次NURBS基函数,对机身形状、机翼、平尾进行综合优化,共116个设计变量,几何约束为机身容积不低于初始的95%;空间网格单元共5900万左右,为准确捕捉近场波系,空间网格拓扑按照激波形态进行x方向拉伸。采用序列二次规划方法进行大规模设计变量低声爆气动外形优化设计,采用128核进行并行计算,SQP优化搜索收敛判断标准为连续三次搜索结果相差量1.0×10-9。
![]() |
图 13 FFD参数化控制框 Fig.13 FFD Lattice |
采用多目标加权优化,优化数学模型如下:
$ \begin{array}{l} {\rm{min}}J = {\rm{ }}{\omega _1}\sum\limits_{{\rm{ }}i = 1}^N {0.5{\alpha _i}{{\left( {{P_i} - {P_{{\rm{target}},i}}} \right)}^2} + } \\ \;\;\;\;\;\;\;\;\;\;\;\;{\omega _2}{(1 - {\rm{ }}\frac{{{C_L}}}{{C_L^*}})^2} + {\omega _3}{(1 - \frac{{{C_D}}}{{C_D^*}})^2}\\ {\rm{st}}{\rm{.Vo}}{{\rm{l}}_{{\rm{body}}}} \ge {\rm{0}}{\rm{.95Vo}}{{\rm{l}}_{{\rm{body,initial}}}}\\ \;\;\;\;{\rm{Thickness}} \ge {\rm{Thicknes}}{{\rm{s}}_{{\rm{initial}}}} \end{array} $ | (3) |
其中,目标函数第一项作用为声爆信号与目标特征的差量最小化,第二、三项是对升力、阻力系数的目标函数,使得在优化过程中气动性能不过于下降,权系数的选择为ω1=1.0、ω2=0.1、ω3=0.1, 均为经验性参数,将升力、阻力进行低权重加权平均,较大程度上是对声爆强度的优化。为简化伴随方程右端项装配过程,将上述目标函数J的变分分解为两部分进行:
1) 声爆伴随变量求解完毕后单独装配到耦合伴随方程的右端项;
2) 升阻力等气动力变分后单独装配到流场伴随方程的右端项。
相当于求解两次伴随方程,分别求解梯度,最后进行加权平均,其中目标函数的第一部分
图 14、图 15分别给出了马赫数、耦合伴随方程第一伴随变量对称面分布云图,从图 14中可以定性分析,该气动外形对声爆目标函数产生主要影响的区域主要集中在机翼前沿翼身结合处、机翼50%根弦长翼身结合处、平尾安装位置以及机身收尾等附近;图 16给出了优化收敛历程,可以看出经过35代优化,基本趋于收敛;图 17(a)为优化前后气动外形的对比,主要变化量集中在机身中段/后段/末尾部分、机翼上下表面,这与上述定性分析一致。图 17(b)是优化前后地面声爆信号、目标信号的对比,经过优化后过压峰值明显降低,声爆信号呈多弱波系状态,整体趋近于目标特征分布,由于声爆信号优化形态依赖于指定的地面过压目标和设计变量分布,因此,下一步将对上述现象影响因素和机理进行系统研究分析。
![]() |
图 14 机体附近空间波系形态 Fig.14 Wave distribution near the body |
![]() |
图 15 耦合伴随方程第一伴随变量云图 Fig.15 The first adjoint variable contour of coupled adjoint equations |
![]() |
图 16 优化设计历程 Fig.16 Optimization design process |
![]() |
图 17 初始与优化外形初气动/声爆特性比较 Fig.17 Aerodynamic/sonic boom feature of the initial and optimized configuration |
近场呈现若干个弱激波/膨胀波的多波系形态,该系列弱波系在传播中得以迅速耗散,图 17(c)给出了优化前后近场相同位置附近波系形态对比,能够清晰看到优化前后的波系演化过程,优化外形波系强度较弱且迅速耗散。图 17(d)为伴随优化前后H/L=0.6处过压分布,优化过后近场过压被明显抑制,结合图 17(a),可以看到机身轴线进一步弯曲、机身末端形状以及平尾的扭转,对第一道激波和尾激波起到了明显的抑制作用,从而远场声爆强度显著下降。
图 17(e)给出了优化前后声压级频谱特性,整个频段的声压级均明显降低;图 17(f)给出了响度级频域特性,本文采用史蒂文斯响度级方法[19],首先将地面声爆时域信号进行快速傅里叶变换,得到每个频段上的声压级(三分之一倍频程),根据各频段中心频率的声压级,采用线性插值确定各频段响度。最后将各个频段的宋值按以下求和方式获得总强度:
$ {S_t} = {S_m} + F(\sum {{S_i} - {S_m}} ) $ | (4) |
其中,St是总强度,Sm是响度最大值,
$ \begin{array}{l} {\rm{(Perceived}}\;{\rm{ noise}}\;{\rm{Level, PLdB)}};\\ \;\;\;\;\;\;\;P = 32 + 9\;{\rm{lo}}{{\rm{g}}_2}{S_t} \end{array} $ | (5) |
经过优化后,声爆感知声压级从初始的90.6 PLdB,降低到79.8 PLdB,效果较为明显。
表 1给出了指定权重系数的选择下,优化前后升阻力、声爆感觉噪声级的对比,优化后外形升力系数仅降低0.005,同时阻力系数降低0.000 34,验证了气动力/声爆一体化设计思路的可行性。
表 1 初始与优化外形气动/声爆特性 Table 1 Aerodynamic/sonic boom feature of the initial and optimized configuration |
![]() |
图 17(g~i)分别为优化前后等效面积分布、横截面积分布以及升力分布对比,其中各项分布无量纲化参考量均采用初始外形对应的最大值,优化后等效面积分布变化更为缓和,最大横截面积位置后移,尾部升力分布变大,带来尾激波强度大幅减小的效应。
5 结论本文开展了全机气动/声爆综合优化工作,验证了考虑升阻力特性的综合设计技术的可行性, 得到以下结论:
1) 连续变后掠机翼外形、弯曲机身轴线、最大横截面后移的面积分布是低声爆气动布局设计的关键。
2) 基于流场/声爆耦合伴随优化全机一体化声爆/气动力优化,感觉噪声降低10.8 PLdB,达到79.8 PLdB,优化效果明显。
3) 从优化前后近场、远场过压分布对比上看,尾部机身形状、平尾扭转的变化对尾部激波具有明显的抑制作用;机身轴线进一步弯曲使第一道激波得到抑制,从而远场声爆强度显著下降。
4) 优化外形呈现多弱波系形态,平行弱波系在传播中得以迅速耗散,减小了声爆强度,在整个频段上响度、声压均得到有效抑制。
5) 综合优化后低声爆外形等效面积分布变化率更为缓和,最大横截面积位置后移,尾部升力分布变大。
文中将升力、阻力进行低权重加权进行,较大程度上是对声爆强度的优化,而低阻、低声爆设计对“面积分布”要求存在差异,不同“面积分布”与阻力/声爆之间相互作用机制与影响规律,是非常值得关注的一项研究内容,后续将结合声爆反演及其耦合伴随方程开展详细研究。
[1] |
NADARAJAH S, JAMESON A, ALONSO J. Sonic boom reduction using an adjoint method for wing-body configurations in supersonic flow[C]//9th AIAA/ISSMO Symposium on Multidisciplinary Analysis and Optimization, Atlanta, Georgia: 2002. doi: 10.2514/6.2002-5547
|
[2] |
RALLABHANDI S K. Sonic boom adjoint methodology and its applications[C]//29th AIAA Applied Aerodynamics Conference, Honolulu, Hawaii, 2011. doi: 10.2514/6.2011-3497
|
[3] |
MUNGUIA B C, ECONOMON T D, ALONS J J. A discrete adjoint framework for low-boom supersonic aircraft shape optimization[C]//18th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference, Denver, Colorado: 2017. doi: 10.2514/6.2017-3326
|
[4] |
张绎典, 黄江涛, 高正红. 基于增广Burgers方程的音爆远场计算及应用[J]. 航空学报, 2018, 39(7): 96-107. ZHANG Y D, HUANG J T, GAO Z H. Far field simulation and applications of sonic boom based on augmented Burgers equation[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(7): 96-107. DOI:10.7527/S1000-6893.2018.22039 (in Chinese) |
[5] |
ZHANG Y D, HUANG J T, GAO Z H, et al. Inverse design of low boom configurations using proper orthogonal decomposition and augmented Burgers equation[J]. Chinese Journal of Aeronautics, 2019, 32(6): 1380-1389. DOI:10.1016/j.cja.2019.03.017 |
[6] |
乔建领, 韩忠华, 丁玉临, 等. 基于广义Burgers方程的超声速客机远场声爆高精度预测方法[J]. 空气动力学学报, 2019, 37(4): 663-674. QIAO J L, HAN Z H, DING Y L, et al. Sonic boom prediction method for supersonic transports based on augmented Burgers equation[J]. Acta Aerodynamica Sinica, 2019, 37(4): 663-674. DOI:10.7638/kqdlxxb-2018.0267 (in Chinese) |
[7] |
乔建领, 韩忠华, 宋文萍. 基于代理模型的高效全局低音爆优化设计方法[J]. 航空学报, 2018, 39(5): 62-75. QIAO J L, HAN Z H, SONG W P. An efficient surrogate-based global optimization for low sonic boom design[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(5): 62-75. DOI:10.7527/S1000-6893.2017.21736 (in Chinese) |
[8] |
冯晓强, 李占科, 宋笔锋. 超声速客机低音爆布局反设计技术研究[J]. 航空学报, 2011, 32(11): 1980-1986. FENG X Q, LI Z K, SONG B F. A research on inverse design method of a lower sonic boom supersonic aircraft configuration[J]. Acta Aeronautica et Astronautica Sinica, 2011, 32(11): 1980-1986. (in Chinese) |
[9] |
冯晓强, 宋笔锋, 李占科, 等. 超声速飞机低声爆布局混合优化方法研究[J]. 航空学报, 2013, 34(8): 1768-1777. FENG X Q, SONG B F, LI Z K, et al. Hybrid optimization approach research for low sonic boom supersonic aircraft config-uration[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(8): 1768-1777. DOI:10.7527/S1000-6893.2013.0057 (in Chinese) |
[10] |
王刚, 马博平, 雷知锦, 等. 典型标模音爆的数值预测与分析[J]. 航空学报, 2018, 39(1): 164-176. WANG G, MA B P, LEI Z J, et al. Simulation and analysis for sonic boom on several benchmark cases[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(1): 164-176. DOI:10.7527/S1000-6893.2017.21458 (in Chinese) |
[11] |
王宇航, 徐悦. 史蒂文斯响度在超声速民机低声爆设计中的应用[J]. 空气动力学学报, 2019, 37(4): 683-689. WANG Y H, XU Y. Application ofStevens'loudness in low﹣boom design of supersonic civil aircraft[J]. Acta Aerodynamica Sinica, 2019, 37(4): 683-689. DOI:10.7638/kqdlxxb-2018.0285 (in Chinese) |
[12] |
郝璇, 苏诚, 刘芳, 等. 超声速飞行器低声爆气动布局优化设计研究[J]. 空气动力学学报, 2018, 36(2): 327-333. HAO X, SU C, LIU F, et al. Optimization design research on low sonic boom configuration for supersonic transport[J]. Acta Aerodynamica Sinica, 2018, 36(2): 327-333. DOI:10.7638/kqdlxxb-2016.0047 (in Chinese) |
[13] |
冷岩, 钱战森, 刘中臣. 超声速条件下旋成体声爆典型影响因素分析[J]. 空气动力学学报, 2019, 37(4): 655-662, 689. LENG Y, QIAN Z S, LIU Z C. Analysis on typical parameters of bodies of revolution affecting the sonic boom[J]. Acta Aerodynamica Sinica, 2019, 37(4): 655-662, 689. DOI:10.7638/kqdlxxb-2018.0291 (in Chinese) |
[14] |
朱自强, 兰世隆. 超声速民机和降低音爆研究[J]. 航空学报, 2015, 36(8): 2507-2528. ZHU Z Q, LAN S L. Study of supersonic commercial transport and reduction of sonic boom[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(8): 2507-2528. DOI:10.7527/S1000-6893.2014.0207 (in Chinese) |
[15] |
兰世隆. 超声速民机声爆理论、预测和最小化方法概述[J]. 空气动力学学报, 2019, 37(4): 646-654, 645. LAN S L. Overview of sonic boom theory, prediction and minimization methods for supersonic civil aircraft[J]. Acta Aerodynamica Sinica, 2019, 37(4): 646-654, 645. DOI:10.7638/kqdlxxb-2018.0294 (in Chinese) |
[16] |
黄江涛, 张绎典, 高正红, 等. 基于流场/声爆耦合伴随方程的超声速公务机声爆优化[J]. 航空学报, 2019, 40(5): 51-61. HUANG J T, ZHANG Y D, GAO Z H, et al. Sonic boom optimization of supersonic jet based on flow/sonic boom coupledadjoint equations[J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(5): 51-61. (in Chinese) |
[17] |
黄江涛, 高正红, 余婧, 等. 大型民用飞机气动外形典型综合设计方法[J]. 航空学报, 2019, 40(2): 52-62. HUANG J T, GAO Z H, YU J, et al. A typical integrated design method for aerodynamic shape optimization of large civil aircraft[J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(2): 52-62. (in Chinese) |
[18] |
RALLABHANDI S. Advanced sonic boom prediction using augmented Burger's equation[C]//49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Orlando, Florida: 2011. doi: 10.2514/6.2011-1278
|
[19] |
STEVENS SS. Perceived level of noise by mark VⅡ and decibels (E)[J]. Journal of the Acoustical Society of America, 1972, 51: 575-601. DOI:10.1121/1.1912880 |