基于直接模拟蒙特卡洛方法的真空羽流不确定量化研究

真实流动环境下的真空羽流必然存在着各种不确定性,那么确定性输入条件的数值模拟必然会存在偏差,因此需进一步研究不确定性对羽流流动特征的影响规律。本文采用直接模拟蒙特卡洛(DSMC)方法,对不确定性输入的羽流流场进行模拟;采用稀疏的概率配置点方法对来流、壁面及模型参数等输入不确定量进行描述,对不确定性的传播和输出目标的平均值、方差及不确定度进行计算。研究表明,流场不确定性沿流线传播至流速最大处之后迅速增强,并在声速线前出现骤减的现象;传播至声速线之后,挡板壁面输入不确定性的影响凸显。其中最为显著的是,压力不确定度在挡板驻点位置达到全场最大值,约为输入不确定量(3.54%)的2.1倍。此外,温度跳跃不确定度受到壁面温度不确定性输入的限制而近似保持为一个恒定值,约为输入不确定度的0.8倍。进而,壁面热流不确定度(5.54%)比壁面正应力不确定度(6.25%)略小,切应力不确定度最小(5.07%)。Sobol’全局敏度分析表明,喉道速度和喉道压力的输入不确定性对气动力/热不确定度的贡献是最大的,且远远超过了壁面温度和模拟分子直径不确定性输入的贡献。     

湍流模型系数不确定度对翼型绕流模拟的影响

使用非嵌入式多项式混沌方法研究了湍流模型系数的不确定度对RAE2822跨声速翼型绕流模拟的影响。计算中关注了数值模拟的积分量（升力系数、阻力系数）和局部量（壁面压力、摩擦系数和空间马赫数分布）的不确定度量化结果。首先，从单输入变量入手，研究卡门常数的不确定度对数值模拟的影响。然后，同时考虑Spalart-Allmaras模型中9个参数的不确定度带来的影响。通过多项式混沌展开，得到系统输出对不确定输入变量的响应，由此可以得到输出的统计特性，包括平均值、方差和极值等信息。最后，在多变量不确定度量化过程中，通过Sobol指标来量化每个输入变量的不确定度对输出不确定度的贡献程度。本文计算只考虑了RAE2822跨声速翼型模拟的单一计算状态，影响规律是否可以推及其他工况和算例需要进一步检验。     

轴对称内锥流动中马赫盘的形成与演化

针对轴对称内收缩直锥流场中的马赫反射问题，采用无黏数值模拟结合理论分析，在来流马赫数6及直锥入口半径R恒定的条件下，以模型壁面前缘角度θ_w和壁面长度w/R为几何参数，研究了内锥形激波汇聚增强过程以及马赫盘下游流动对马赫盘位置的影响。结果表明，壁面前缘角度θw通过改变轴对称内锥形入射激波汇聚增强过程中达到von Neumann强度的位置，限定了马赫盘可能存在的位置范围；马赫盘位置还受其下游流动的影响，在不同的θ_w和w/R条件下，马赫盘下游流管内声速喉道的形成机制可以划分为依赖或不依赖壁面尾缘膨胀波两种类型；对于不依赖于壁面尾缘膨胀波的情况，声速喉道形成之前流管内的压力仅与反射激波的波后非均匀压力相匹配，马赫盘的位置不随w/R而变化；对于依赖于壁面尾缘膨胀波的情况，声速喉道形成之前流管还需要匹配壁面尾缘膨胀波引起的压力变化；根据反射激波下游与马赫盘下游流管的压力匹配关系，可以求解马赫盘位置。     

高超声速喷管气膜冷却流动与传热特性研究

基于雷诺时均方法并结合SST k-ω湍流模型，对空气气膜作用下的马赫数6高超声速喷管流动与传热特性进行了数值研究，分析了气膜流量、气膜狭缝几何参数（狭缝高度、台阶厚度）、狭缝流向位置对喷管壁面尤其是喉道热流密度以及对喷管出口气流品质的影响特性。计算结果表明，引入气膜可以显著降低喷管尤其是喉道的壁面热流，而气膜对喷管出口气流速度分布的影响很小，对出口温度分布有一定的影响。当气膜流量仅为主流的3.53%时，喷管喉道热流的降幅达30.1%；增加这一值至14.25%时，喷管喉道热流的降幅升至87.3%。同时，改变狭缝流向位置以及狭缝高度对气膜冷却效率有一定的影响，而改变台阶厚度对冷却效率的影响很小。     

磁场对高超声速弱电离气体流动的影响

对偶极子磁场作用下的三维钝头体高超声速黏性绕流的化学非平衡流动进行了数值模拟。电导率利用组分公式计算,化学模型为7组元、6反应模型。应用诱导磁场方法将磁场与流场耦合,控制方程的空间离散采用有限体积法,扩散项用中心差分格式,对流项采用二阶迎风格式,时间推进采用3步Runge-Kutta法。计算结果表明,外加偶极子磁场使激波脱体距离增加,壁面摩擦系数和表面热流密度增大。与无磁场作用时相比,在0.153T外加磁场作用下,冻结流中的激波脱体距离增加约3倍,局部壁面摩擦系数最大增加63%,局部表面热流密度最大增加61%;而化学非平衡流中的激波脱体距离增加约0.5倍,局部壁面摩擦系数最大增加47%,局部表面热流密度最大增加31%,且化学非平衡流中激波层内温度的最大值约为冻结流中的64%。     

尾迹/势流干扰对跨声速风扇静叶附面层流动影响的数值研究

为揭示尾迹/势流干扰下压气机静叶附面层的非定常流动机理,采用数值方法对单级高负荷跨声速风扇中径处的非定常流场进行模拟,针对尾迹/势流干扰下的静叶附面层非定常流动特征进行研究。基于尾迹/势流与叶片附面层干扰模型,通过叶片壁面摩擦力、近壁面附面层湍动能和壁面静压脉动,详细分析了尾迹和势流干扰下尾迹对高负荷静叶附面层流动状态的影响。研究发现:在非定常条件下,尾迹干扰能够"刺穿"静叶附面层,使得静叶附面层发生跨越转捩现象。尾迹和受其诱导的势流碰撞静叶前缘产生的压力扰动波在压力面附面层内以声速传播,影响压力面整个弦长的静压和摩擦力分布,而在吸力面上,扰动波的传播仅局限在前缘区域部分弦长上,吸力面和压力面传播现象不同与逆压梯度和曲面凸、凹形状相关。     

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

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

二元喉道倾斜矢量喷管调节方法

利用数值模拟方法,对二元喉道倾斜矢量喷管的调节方法进行了研究。研究结果表明,当喉道注气流量不变,随着扩张段注气流量逐渐增大,喷管内出现三种典型流动状态。当扩张段注气流量比较小时,喉道处声速线的形状和位置没有明显变化,在扩张段注气口前面形成亚声速区域;当喷管扩张段注气流量比较大时,喷管中出现了强激波系和上下贯通的亚声速区域,这时喷管内的流动损失最大,推力系数最低;进一步增大扩张段注气流量,声速线会从喷管喉道移动到喉道注气口和扩张段注气口之间,出现典型的喉道倾斜现象。单独控制流量系数可通过保持扩张段注气流量不变,改变喉道注气流量来实现;矢量角的控制可通过保持喉道注气流量不变,改变扩张段注气流量来实现。     

基于改进SST模型的分离流动数值模拟

对SST湍流模型进行了改进,并通过对典型分离流动进行数值模拟,来检验和提高模型预测分离流动的能力.基于亚声速分离流动,提出了提高雷诺应力的模拟精度和分离流修正的改进方法,并进行了对比研究得出结论;在亚声速分离流动分析结论基础上,采用了可压缩性修正方法,开展了跨声速、超声速和高超声速激波/边界层干扰分离流动的数值模拟研究.结果表明:提高雷诺应力的模拟精度和采用分离流修正明显地提高了SST湍流模型对分离流动的模拟能力;适当地可压缩性修正对超声速和高超声速分离流动的计算精度有改善作用.     

均匀磁场中高超声速弱电离气体流动数值模拟

采用7组元化学模型并应用组分公式计算电导率,通过求解黏性MHD(magneto-hydro-dy-namics)方程组,研究了不同强度均匀磁场对三维钝头体高超声速绕流化学非平衡流动的影响.结果表明,随着磁场的增强,激波脱体距离逐渐增加;总阻力系数和壁面温度逐渐减小.在B_y=0.03 T磁场作用下,与无磁场的结果相比,化学非平衡流中的激波脱体距离增加约7%,总阻力系数减小约5%,局部壁面温度最大降低74%;而冻结流中的激波脱体距离增加约43%,总阻力系数减小约6.9%,局部壁面温度最大降低18%.      

有翼高超声速再入飞行器气动设计难点问题

有翼高超声速再入飞行器是近年来的研究热点,气动设计是飞行器设计的关键。为了更清楚地认识有翼高超声速再入飞行器气动设计的难点问题,对有翼高超声速再入飞行器的发展、优势及总体任务剖面进行了介绍,从5个方面详细介绍了该类飞行器气动设计的难点问题,包括多约束复杂面对称气动布局设计、高温真实气体效应对气动特性影响、天地差异与天地换算方法、反作用控制系统(RCS)喷流干扰对气动特性的影响以及气动数据不确定度等,简要阐明了这些难点问题对总体设计的重要性以及初步的解决思路,为有翼高超声速再入飞行器气动设计提供了一些参考。     

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

Simulation and Analysis of Crashworthiness of Fuel Tank for Helicopters

Crashworthiness requirement of fuel tanks is one of the important requirements in helicopter designs. The relations among the protection frame, textile layer and rubber layer of the fuel tank are introduced. Two appropriate FE models are established, one is for an uncovered helicopter fuel tank without protection frame, and the other is for fuel tank with protection frame. The dynamic responses of the two types of fuel tanks impinging on the ground with velocities of 17.3 m/s are numerically simulated for the purpose of analyzing energy-absorbing capabilities of the textile layer and protection frame. The feasibility of the current crashworthiness design of the fuel tank is examined though comparing the dynamic response behaviors of the two fuel tanks.     

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

突破技术趋同，波音再现竞争优势

 波音公司面临着来自空客公司的巨大挑战,企业战略性创新才是公司成功的关键。为此波音公司的全部战略性研究集中在扩大产品的差异性上,体现在3个方面：电子化（e-Enabled）运营环境、整体复合材料机身部件的制造技术和支持波音787客机的全球协同环境（GCE）。     

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。     

基于弯曲激波压缩系统的高超声速进气道反设计研究进展

总结了近十年来弯曲激波压缩研究的主要成果。提出了弯曲激波压缩系统的新概念,即利用特殊设计的楔形弯曲压缩面或空间弯曲压缩面,产生一系列与前缘弱激波相互交汇或叠加的压缩波系,从而使前缘激波弯曲,形成特殊的弯曲激波,它与波后的等熵压缩波来共同完成对气流的压缩。在此基础上,实现了由给定出口气动参数的超声速内流道反设计,实现了由给定压缩面压力分布和给定压缩面马赫数分布要求的型面反设计,实现了由给定激波波面的压缩型面反设计。研究证明,弯曲压缩面-弯曲激波压缩系统具有良好的综合气动性能,为高性能高超声速进气系统的气动设计提供了一种全新的设计方法。     

临近空间飞行器测控与信息传输系统频段选择

 临近空间飞行器是高性能信息化武器平台,测控（TT&;C）与信息传输系统是其信息保障的核心,而选择合理、可行的频段是展开系统设计的前提和基础。频段选择影响到整个技术方案的制定,是一个需综合考虑、影响深远并具有战略意义的关键问题,从国际电联（ITU）国际标准、高速数传、接收信噪比（SNR）、“三抗”、超视距中继、黑障、雨衰以及设备研制成熟度8个方面全面、细致论证了近空间平台测控系统的频段选择问题,最终得出在视距链路中以Ka频段为宜,在超视距链路中以Ku/Ka双频段为宜的结论。     

基于H∞性能指标的质量矩拦截弹鲁棒控制

 以所建立的质量矩拦截弹数学模型为基础，利用微分几何的反馈线性化理论，得到一个解耦线性系统，考虑到拦截弹的鲁棒性要求和3个滑块的协调控制问题，提出采用双回路的设计方法，内回路采用线性二次调节器（LQR），外回路采用考虑混合灵敏度问题的H_∞控制设计。仿真结果证明了该方法动态性能满足设计要求，同时对系统参数的不确定性具有较强的鲁棒性。