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

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

考虑配平特性的超声速客机低声爆气动布局优化研究

降低声爆水平是下一代超声速运输机研制需要解决的关键问题之一。低声爆优化通常使飞行器布局向着机翼后掠角增大、机翼沿机身方向分布范围增大的趋势发展,给飞行器的配平和低速特性带来不利影响。以某超声速客机基本构型为研究对象,建立基于类别/形状函数的翼身组合体参数化建模方法;基于超声速线化理论分析外形几何参数对声爆水平的影响。在此基础上,分别针对机身轮廓、机翼平面形状以及扭转角分布对该构型进行低声爆优化和俯仰力矩特性优化,并采用CFD 方法对优化结果进行校核。结果表明：与基准构型相比,在不显著增加俯仰力矩的基础上,优化构型的阻力降低了19 cts,近场过压显著降低,地面声爆响度降低5.1 PLdB。     

超声速客机低音爆布局反设计技术研究

音爆已经成为限制民用飞机在陆地上空进行超声速飞行最关键的因素.降低超声速客机的音爆水平,使其能够在陆地上空超声速飞行,将会给超声速客机带来巨大的潜在市场.基于SGD(Seebass-George-Darden)方法,构建了相关的设计分析环境,对超声速客机低音爆布局的反设计技术作了研究分析,在此基础上首次提出了一种“梭式...     

超声速低声爆布局分层优化方法

声爆抑制是发展新一代超声速民机必须突破的关键技术。总体布局参数的合理设计可以使飞行器具有良好的声爆特性。为了提高全局进化算法在布局设计中的优化效率,提出一种基于数据挖掘的分层优化方法,利用数据挖掘中的决策树算法提取设计知识,获得设计变量分层信息,指导低声爆布局分层优化;针对某超声速低声爆飞行器,选取后掠角、上反角、展弦比、梢根比、长细比五个总体布局参数作为设计变量,开展分层优化数值实验,并与一体化优化形成对比验证。结果表明：分层优化方法能够搜索到与一体化优化高度吻合的最优解,分层优化的收敛速度显著高于一体化优化,且对随机寻优历程的表现更稳健。     

低声爆静音锥设计方法研究

 声爆精确预测及低声爆设计方法已成为新一代军民用超声速飞机研制过程中必须解决的关键难题之一。基于计算流体力学(CFD)、波形参数法以及MARK-VII方法构建了高精度声爆预测方法,利用该方法对低声爆静音锥的设计展开研究。研究结果表明,静音锥的设计存在临界长度,静音锥长度小于临界长度时,静音锥产生的激波在传播过程中与机头弓形激波合并,静音锥无法起到降低声爆的作用;静音锥长度大于临界长度时,声爆水平也会略有上升。静音锥临界长度随飞行高度和飞行马赫数的变化而变化,可以根据实际飞行状态采用可伸缩设计,达到最佳的低声爆效果。多级静音锥利用多道弱激波取代机头强弓形激波,其声爆水平较单级静音锥也更低,同样,多级静音锥每一级的长度也要达到临界长度。不同静音锥头部形状产生的脱体激波形状不同,脱体距离也不同,导致阻力系数以及静音锥壁面温度有所不同,但静音锥头部形状对远场声爆信号的影响并不明显。采用静音锥的低声爆方案与原始方案比较,声爆水平得到大幅降低,阻力系数略有上升。     

Inverse design of low boom configurations using proper orthogonal decomposition and augmented Burgers equation

Mitigation of sonic boom to an acceptable stage is a key point for the next generation of supersonic transports. Meanwhile, designing a supersonic aircraft with an ideal ground signature is always the focus of research on sonic boom reduction. This paper presents an inverse design approach to optimize the near-field signature of an aircraft, making it close to the shaped ideal ground signature after the propagation in the atmosphere. Using the Proper Orthogonal Decomposition(POD) method, a guessed input of augmented Burgers equation is inversely achieved. By multiple POD iterations, the guessed ground signatures successively approach the target ground signature until the convergence criteria is reached. Finally, the corresponding equivalent area distribution is calculated from the optimal near-field signature through the classical Whitham F-function theory. To validate this method, an optimization example of Lockheed Martin 1021 is demonstrated. The modified configuration has a fully shaped ground signature and achieves a drop of perceived loudness by 7.94 PLdB. This improvement is achieved via shaping the original near-field signature into wiggles and damping it by atmospheric attenuation. At last, a nonphysical ground signature is set as the target to test the robustness of this inverse design method and shows that this method is robust enough for various inputs.     

超音速客机音爆问题初步研究

在超音速客机设计中,音爆已成为决定和衡量方案成败的关键技术指标之一,也是超音速客机设计必须解决的关键技术问题。利用简化音爆预测(SSBP)方法对超音速客机的音爆问题进行了初步研究,实现了音爆特征的定量化计算分析,开发了基于Matlab GUI的音爆特征分析计算软件,通过算例验证了所采用方法的正确性。最后,利用所开发的音爆分析计算软件对一种新概念超音速静音客机的音爆水平进行了预测。     

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...     

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

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

隐身性能约束的多目标气动外形优化设计

从飞机外形设计的总体、气动和隐身设计要求出发,根据多目标优化的基本概念,将Pareto方法与遗传优化搜索相结合,并采用了群体排序、基于共享机制的小生境技术和Pareto解集过滤器等技术使解集具有良好分布特性,在此基础上建立了一套可满足飞行器外形气动/隐身一体化综合优化设计优化模型和优化方法。文中针对飞行器外形优化设计要求,提出了复杂外形参数化和设计变量的选取原则。并根据某飞行器设计要求,进行了在以隐身特性为约束条件下,以亚声速和超声速气动特性为设计目标的飞行器外形综合优化设计,取得了良好的优化设计结果。     

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

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

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

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

基于代理模型的气动外形平面参数多目标匹配设计

将Kriging代理模型和Pareto遗传算法引入气动外形平面参数匹配设计中,提出一种基于代理模型的多目标平面参数匹配设计方法。将拉丁超立方试验设计用于平面参数筛选,确定出参数匹配方案库;基于方案库的计算流体力学(CFD)分析结果构建Kriging气动代理模型;将Kriging模型替代CFD分析,用于气体布局参数匹配优化设计,提高了设计效率并保证了可信度;通过Pareto遗传算法优化解决多点设计要求下气动布局参数匹配问题,一次性给出参数匹配方案的最优解集,从Pareto前沿中根据设计偏向选择气动布局最佳匹配方案。以典型的双后掠布局平面参数多点匹配优化设计问题作为算例,研究结果表明:Kriging气动代理模型与实际CFD分析结果的误差均小于5%,满足精度要求;根据不同设计偏向,选择了3种参数匹配Pareto优化方案,与原样本方案相比超声速阻力减小6.0%~12.8%,跨声速升阻比增加0.01%~3.40%,证明了匹配设计方法的有效性;通过试验设计的Pareto分析与主、交互效应分析,获得了气动布局平面参数对气动性能影响的定量关系,能够为参数匹配设计提供依据。所提出的平面参数匹配设计方法可应用于其他常规与非常规气动布局型式。     

细长杆用于超音速飞行器降噪的数值分析

在超音速飞机头部加装合理设计的细长杆可有效降低飞行噪音.基于Fluent准三维数值模拟手段对多组细长杆设计方案进行模拟分析,对比其对远场压力的影响,以比较各类细长杆的降噪效果,根据气动理论对数值结果作出解释,揭示细长杆的降噪效果与其外形的内在联系;从远场压力最大值、远场压力曲线形状等不同方面分析细长杆对音爆的影响.结果表明:6 m长的12.5°半锥角单级锥形细长杆有明显的降噪效果;合理设计的多级细长杆可明显降低初始过压值,并能延迟远场压力达到最大值的时间;凹凹面、凹锥面、凸凹面和凸锥面的细长杆能够降低远场压力峰值,改善压力曲线的形状,具有较好的降噪作用;细长杆产生的强膨胀波能在激波传播过程中大幅削弱机头激波,从而削弱音爆水平.     

A survey on numerical simulations of drag and heat reduction mechanism in supersonic/hypersonic flows

Along with the survey on experimental investigations drawing attention to the drag and heat reduction mechanism, the authors simultaneously focus on the recent advances of numerical simulations on the schemes applied to supersonic/hypersonic vehicles. The CFD study has evolved as an irreplaceable method in scheme evaluation and aircraft optimization. Similar to our previous experimental survey, the advances in drag and heat reduction schemes are reviewed by similar kinds of mechanism in this article, namely the forward-facing cavity, the opposing jet, the aerospike, the energy deposition and their combinational configurations. This review article puts an emphatic eye on the flow conditions, numerical methods, novel schemes and analytical conclusions given in the simulations. Further, the multi-objective design optimization concept has also been illustrated due to the observable advantages of using CFD over experimental method, especially those performances conducted in drag reduction and thermal protection practice, and this would possess reference value in the design of aircraft system.     

超声速翼身组合体激波阻力优化的EFCE算法

 超声速飞行器的横截面积分布对其激波阻力的影响十分显著,合理的机翼和机身横截面积分布可以显著降低其激波阻力。使用类别形状函数变换(CST)方法对机身进行基于横截面积分解的CST参数化外形表示,在此基础上提出了扩展的远场组元(EFCE)超声速翼身组合体激波阻力优化算法,并使用该方法对超声速客机翼身组合体进行外形优化,使其激波阻力系数降低了39%。研究结果表明:由于只进行一个方向上的面积分解,机身CST参数化所使用的参数数量和相应优化过程的计算量比机翼大幅降低;经过EFCE激波阻力优化的机身具有较为明显的面积率修形"蜂腰"特征。     

Supersonic biplane—A review

One of the fundamental problems preventing commercial transport aircraft from supersonic flight is the generation of strong sonic booms. Sonic booms are the ground-level manifestation of shock waves created by airplanes flying at supersonic speeds. The strength of the shock waves generated by an aircraft flying at supersonic speed is a direct function of both the aircraft’s weight and its occupying volume; it has been very difficult to sufficiently reduce the shock waves generated by the heavier and larger conventional supersonic transport (SST) configuration to meet acceptable at-ground sonic-boom levels. It is our dream to develop a quiet SST aircraft that can carry more than 100 passengers while meeting acceptable at-ground sonic-boom levels. We have started a supersonic-biplane project at Tohoku University since 2004. We meet the challenge of quiet SST flight by extending the classic two-dimensional (2-D) Busemann biplane concept to a 3-D supersonic-biplane wing that effectively reduces the shock waves generated by the aircraft. A lifted airfoil at supersonic speeds, in general, generates shock waves (therefore, wave drag) through two fundamentally different mechanisms. One is due to the airfoil’s lift, and the other is due to its thickness. Multi-airfoil configurations can reduce wave drag by redistributing the system’s total lift among the individual airfoil elements, knowing that wave drag of an airfoil element is proportional to the square of its lift. Likewise, the wave drag due to airfoil thickness can also be nearly eliminated using the Busemann biplane concept, which promotes favorable wave interactions between two neighboring airfoil elements. One of the main objectives of our supersonic-biplane study is, with the help of modern computational fluid dynamics (CFD) tools, to find biplane configurations that simultaneously exhibit both traits. We first re-analyzed using CFD tools, the classic Busemann biplane configurations to understand its basic wave-cancellation concept. We then designed a 2-D supersonic biplane that exhibits both wave-reduction and cancellation effects simultaneously, utilizing an inverse-design method. The designed supersonic biplane not only showed the desired aerodynamic characteristics at its design condition but also outperformed a zero-thickness flat-plate airfoil. (Zero-thickness flat-plate airfoils are known as the most efficient monoplane airfoil at supersonic speeds.) Also discussed in this paper is how to design 2-D biplanes, not only at their design Mach numbers but also at off-design conditions. Supersonic biplanes have unacceptable characteristics at their off-design conditions such as flow choking and its related hysteresis problems. Flow choking causes rapid increase of wave drag and it continues to be kept up to the Mach numbers greater the cruise (design) Mach numbers due to its hysteresis. Some wing devices such as slats and flaps, which could be used at take-off and landing conditions as high-lift devices, were utilized to overcome these off-design problems. Then supersonic-biplane airfoils were extended to 3-D wings. Because that rectangular-shaped 3-D biplane wings showed undesirable aerodynamic characteristics at their wingtips, a tapered-wing planform was chosen for the study. A 3-D biplane wing having a taper ratio and aspect ratio of 0.25 and 5.12, respectively, was designed utilizing the inverse-design method. Aerodynamic characteristics of the designed biplane wing were further improved by using winglets at its wingtips. Flow choking and its hysteresis problems, however, occurred at their off-design conditions. It was shown that these off-design problems could also be resolved by utilizing slats and flaps. Finally, a study on the aerodynamic characteristics of wing-body configurations was conducted using the tapered biplane wing. In this study a body was chosen in order to generate strong shock waves at its nose region. Preliminary parametric studies on the interference effects between the body and the tapered biplane wing were performed by choosing several different wing locations on the body. From this study, it can be concluded that the aerodynamic characteristics of the tapered biplane wing are minimally affected by the disturbances generated from the body, and that the biplane wing shows promise for quiet commercial supersonic transport.     

Aerodynamic multi-objective integrated optimization based on principal component analysis

Based on improved multi-objective particle swarm optimization (MOPSO) algorithm with principal component analysis (PCA) methodology,an efficient high-dimension multiobjective optimization method is proposed,which,as the purpose of this paper,aims to improve the convergence of Pareto front in multi-objective optimization design.The mathematical efficiency,the physical reasonableness and the reliability in dealing with redundant objectives of PCA are verified by typical DTLZ5 test function and multi-objective correlation analysis of supercritical airfoil,and the proposed method is integrated into aircraft multi-disciplinary design (AMDEsign) platform,which contains aerodynamics,stealth and structure weight analysis and optimization module.Then the proposed method is used for the multi-point integrated aerodynamic optimization of a wide-body passenger aircraft,in which the redundant objectives identified by PCA are transformed to optimization constraints,and several design methods are compared.The design results illustrate that the strategy used in this paper is sufficient and multi-point design requirements of the passenger aircraft are reached.The visualization level of non-dominant Pareto set is improved by effectively reducing the dimension without losing the primary feature of the problem.     

伸缩机翼变体飞机气动特性研究

伸缩机翼变体飞机通过机翼伸缩调整机翼展长,从而改变机翼面积和展弦比,改变飞机的气动布局和机翼的气动特性,满足多任务点的设计要求。简要介绍伸缩机翼变体飞机的发展历史,重点研究一种采用伸缩机翼设计的超音速飞机的气动特性变化。研究结果表明：亚音速时机翼展长伸长,展弦比增大,飞机诱导阻力降低,升阻比提高,可以明显提高飞机的航程;超音速时机翼展长缩短,展弦比减小,飞机的波阻降低,升阻比增大,提高了超音速飞行性能。伸缩机翼概念用于超音速飞机设计时能很好地兼顾亚音速巡航和超音速冲刺。     

Far-field sonic boom prediction considering atmospheric turbulence effects: An improved approach

Accurate prediction of sonic boom is one of key challenges for the design of a low-boom supersonic aircraft. For most of available far-field prediction methods, the effect of atmospheric turbulence appearing in the planetary boundary layer cannot be considered, which results in remarkable inaccuracy of predicting ground-level sonic boom waveform. Although some efforts have been made to overcome the shortcoming, the turbulence effects are not yet well described so far. This article proposes an improved method by extending the two-dimensional Heterogeneous One-Way Approximation for the Resolution of Diffraction (HOWARD) equation to account for the axial and transverse convections of wind fluctuation as well as the effect of temperature fluctuation. The proposed method is validated by comparing the predictions with the flight-test data of JAXA D-SEND#1 LBM, which shows that the result of the proposed method is in better agreement with the flight-test data than that of the method without considering atmospheric turbulence effects. Then, distortion mechanism of sonic boom waveforms caused by atmospheric turbulence is analyzed by using the proposed method. It is indicated that the effect of turbulent convection makes uniform sonic-boom wavefronts irregular, which creates the condition of diffraction effect to perturb waveforms. Finally, the proposed method is applied to investigate the behavior of two types of waveforms given by the sonic boom minimization theory. Results show that a far-field waveform with a weaker initial shock is more beneficial for low-boom design of a supersonic aircraft.