静音锥对超声速客机声爆水平的影响分析研究

静音锥低声爆的基本原理是通过在超声速飞机头部加装静音锥将机头产生的强激波转化为一系列互不叠加的弱激波,从而达到降低声爆的目的。论文以一种“梭式”布局的超声速客机为基本模型,采用CFD和波形参数法相结合的方法,研究了不同参数的单级和多级静音锥对超声速客机声爆水平的影响。结果表明,静音锥的长度可以调节激波的干涉程度;静音锥的直径和圆锥顶角能改变静音锥的初始超压值;静音锥级数对上升时间影响显著。     

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

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

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

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

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

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

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

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

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

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

新一代环保型超声速客机气动相关关键技术与研究进展

更快的旅行速度是人类永恒的追求。虽然以"协和"号和"图144"为代表的第一代超声速客机商业运营失败,但之后人类从来就没有停止过对新一代更加环保的超声速客机的探索与研究。本文首先梳理总结了其中首要突破的四大关键技术(声爆预测及其抑制技术、超声速减阻技术、变循环发动机技术、低声爆低阻布局与综合优化设计技术),并对其国内外研究进展情况进行了文献综述,对研究现状进行了分析。其次,介绍了西北工业大学超声速客机研究中心在声爆预测理论与方法、声爆抑制技术、低声爆低阻布局与综合优化技术、超声速层流减阻技术等方面的研究进展。最后,针对发展新一代环保型超声速客机当前急需突破的关键科学与技术问题,探讨了未来需重点研究的方向。     

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

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

发动机喷管羽流对近场声爆特性影响的风洞试验技术

超声速飞行所引发的声爆问题是困扰新一代环保型超声速客机发展的关键技术难题，发动机喷管羽流对全机声爆特性尤其是后激波特性具有重要影响。设计了单喷管喷流试验模型及声爆试验装置，评估了风洞试验段洞壁反射激波对模型近场压力测量的影响，重点针对通气支臂对喷管羽流的支撑干扰问题进行了分析与优化。基于中国航空工业空气动力研究院FL-60风洞，开展了发动机喷管羽流对旋成体单喷管模型近场声爆特性影响试验技术研究，试验来流马赫数2.0、落压比（NPR）范围1～20.39。研究结果表明，通过对来流马赫数、通气支臂外形、喷流模型长度、通气支臂与模型的相对位置等参数的综合优化，消除了通气支臂带来的支撑干扰对喷管羽流的影响，确保在风洞试验段受限空间内模型近场压力测量不受洞壁反射和通气支臂波系的影响；喷管羽流主要对模型尾部的近场压力特征产生影响，在来流马赫数一定的条件下，提高喷管NPR使喷流状态从过膨胀到欠膨胀，喷管唇口激波逐渐增强、位置逐渐向上游移动，抑制了喷管船尾膨胀波的发展。     

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

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

不同密度梯度的多层流体界面上的Richtmyer-Meshkov不稳定性研究

在气-液界面Richtmyer-Meshkov（R-M）不稳定性研究的基础上,利用不同流体间密度的差异,构造了空气-硅油-水、空气-酒精-硅油两种流体粘度方向迥异的气-液-液三相界面,实验中以氮气作为高压驱动气体,在不同激波马赫数下对这两种界面R-M不稳定性后期尖钉与气泡区发展进行了测量与统计,通过对实验数据的分析,得出了相关规律,并用这些规律与已有的单层气-液界面R-M不稳定性研究成果作比较,得出了异同点。同时,还研究了两种三相界面在R-M不稳定发展中的差异。实验结果表明：当流体的粘度梯度方向（从小到大的方向）与激波方向一致时,界面失稳更加明显,湍流混合更为显著。     

激波抛撒水膜成雾参数实验研究

使用激光相位多普勒粒子测速系统对激波抛撒水膜形成的云团进行了观测,得到激波抛撒水膜形成的云雾中液滴速率和尺寸参数。结果表明：不同厚度水膜在同等强度激波的作用下抛撒,破碎产生的液滴轴向速率和尺寸以及形成云团的高度都随水膜厚度的增加呈“W”形变化。不同厚度水膜其破碎液滴轴向运动速率存在显著差异。而径向运动速率差距不大。结果说明,在激波作用初期,激波驱动力起着关键作用．随着抛撒距离的增加．液滴的蘑力及在气流中所受摩擦阻力的影响越大。液膜厚度越小,透射激波越强,液膜被抛撒的范围越大。形成云团的直径随着液滴的径向速率增加而变大。     

冷气喷射法控制激波强度的数值研究

对冷气喷射时激波受到的影响进行了多方案数值研究.在高压级静叶吸力面反射激波生成点前后5个不同位置上采用相同总压及相同温度的冷气喷射,分析了不同位置的冷气喷射对激波强度和方向影响.结果表明:在吸力面激波折射点附近喷射冷气是一种行之有效的控制激波强度的方法.冷气喷射位置位于激波折射点附近时能够对激波的强度产生影响,在接近激波折射点前部位置注入冷气对流动有积极作用,减弱了激波.      

波转子非定常泄漏流动机理

针对制约波转子性能的泄漏问题,通过提取波转子中与泄漏相关的主要流动现象并建立波转子单通道泄漏模型,对非定常泄漏流动机制进行了详细的数值研究,对本文给出的非定常泄漏损失预测模型进行了数值验证。结果表明:转子通道中存在不同程度的压力波动,波动幅值与间隙宽度有关;连续反射膨胀波、周期性出现的弓形激波及其反射激波是压力波动的根本原因;间隙内部泄漏过程存在3个主要的流动阶段;泄漏过程中通道激波传播速度不变、波后时均压力不变;在一定间隙宽度范围内,激波马赫数、激波静增压比与无量纲间隙宽度均呈线性关系,当间隙宽度从0增大到0.08时,激波马赫数衰减7.3%,激波静增压比衰减10.1%;泄漏流动通过泄漏产生的主膨胀波对激波传播过程施加影响,通道激波衰减本质上是理想激波与主膨胀波叠加效应的结果;泄漏损失预测模型与数值结果吻合良好。     

水下强声波脉冲引发的空化气泡运动和声辐射

通过建立由等离子体声源、PCB压力传感器和高速相机等构成的实验测量系统,研究了强声波脉冲引发的空化气泡运动和二次声辐射。通过实验,观察到了强声波脉冲特有的波形结构,空化气泡的成长、膨胀和坍缩过程以及气泡坍缩时辐射的声信号。其次,利用Gilmore方程和Bernoulli方程,从理论上对空化气泡的运动和声辐射进行了数值模拟和分析。研究结果表明：（1）在强声波脉冲正压区的作用下,空化气泡将受到压缩并围绕＂准平衡态半径＂振荡;在强声波脉冲负压区的作用下,空化气泡将出现膨胀、坍缩和回弹的物理过程。（2）空化气泡坍缩时周期性辐射的声脉冲持续时间极短,具有＂冲击波＂的特点。     

旋转冲压压缩转子二维进气流道数值研究

借鉴二维超声速进气道的设计方法,设计了一种旋转冲压压缩转子的二维进气流道,并采用二维雷诺平均N-S方程和Spalart-Allmaras湍流模型对其流场进行了数值仿真,研究了转速、背压对二维进气流道中波系结构、内部流动特性和性能的影响.计算结果表明:所提出的旋转冲压压缩转子二维进气流道设计方法是可行的;二维进气流道中产生的激波系与二维超声速进气道中产生的激波系相类似,所不同的是二维进气流道中产生的激波为弯曲激波,而二维超声速进气道中产生的激波为平直激波;转速和背压对二维进气流道的性能有较大影响.      

Shock wave configurations and reflection hysteresis outside a planar Laval nozzle

When the pressure ratio increases from the perfectly expanded condition to the third limited condition in which a normal shock is located on the exit plane, shock wave configurations outside the nozzle can be further assorted as no shock wave on the perfectly expanded condition, weak oblique shock reflection in the regular reflection(RR) pressure ratio condition, shock reflection hysteresis in the dual-solution domain of pressure ratio condition, Mach disk configurations in the Mach reflection(MR) pressure ratio condition, the strong oblique shock wave configurations in the corresponding condition, and a normal shock forms on the exit plane in the third limited condition. Every critical pressure ratio, especially under regular reflection and Mach reflection pressure ratio conditions, is deduced in the paper according to shock wave reflection theory. A hysteresis phenomenon is also theoretically possible in the dual-solution domain. For a planar Laval nozzle with the cross-section area ratio being 5, different critical pressure ratios are counted in these conditions, and numerical simulations are made to demonstrate these various shock wave configurations outside the nozzle. Theoretical analysis and numerical simulations are made to get a more detailed understanding about the shock wave structures outside a Laval nozzle and the RRMMR transition in the dual-solution domain.     

水下爆轰燃气泡形态与激波传播过程研究*(爆震专刊)

针对脉冲爆轰发动机在水下工作过程中形成的燃气射流问题,搭建了水下爆轰燃气实验系统,研究了第一个爆轰循环在水下的燃气泡发展变化过程。建立了基于气液两相双流体模型的脉冲爆轰发动机水下喷射模型,采用时-空守恒元和求解元方法,模拟了爆轰波退化为激波在水中的传播及衰减过程。研究结果表明：燃气泡前期受外界水环境阻滞作用呈现“豌豆状”形态,充分发展阶段气液交界面逐渐失稳,在达到最大尺寸后开始收缩并在中心轴线位置出现凸出的射流；水下爆轰燃气射流发展过程中同时存在脱离燃气泡的水中前导激波和管口燃气泡内的高压区两部分,水中前导激波在传播过程中压力迅速衰减至常压量级,而管口燃气泡内则一直保持较高压力；中心轴线区域气液交界面处反射激波的回传使管口附近出现回击现象,并导致前导激波波阵面上压力峰值逐渐出现在30°方向上。     

旋转爆震发动机喷管非定常流动特性的数值模拟研究

为分析旋转爆震喷管内非定常流动特性,基于特征线方法和旋转爆震燃烧室出口时均参数设计了一种塞式喷管模型,通过改变塞锥式喷管的进出口压比π(π=15,30,45),分析了喷管内非定常流场结构与喷管工作性能之间的相互关系,并探究了喷管进出口压比变化对喷管工作性能的影响.研究结果表明:旋转爆震燃烧室下游斜激波进入喷管后,在喷管...