基于DSMC方法的真空羽流不确定量化研究-空间推进专刊*

真实流动环境下的真空羽流必然存在着多种不确定性,那么确定条件输入的数值模拟会存在偏差,因此需进一步研究不确定性对羽流流动特征的影响规律。本文采用直接模拟蒙特卡洛（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而变化；对于依赖于壁面尾缘膨胀波的情况，声速喉道形成之前流管还需要匹配壁面尾缘膨胀波引起的压力变化；根据反射激波下游与马赫盘下游流管的压力匹配关系，可以求解马赫盘位置。     

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

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

洛仑兹力控制高超声速进气道边界层分离的数值模拟

采用空间HUE格式、时间LU-SGS推进、sst-kw湍流模型、多块结构网格程序,对磁流体动力学(Magnetohydrodynamic:MHD)控制高超声速二维进气道边界层分离进行了数值研究.研究发现,不施加控制时,数值模拟得到的壁面静压和实验结果符合良好,进气道喉道处分离区占据喉道高度的1/3左右.通过施加MHD控制,消除了进气道内部的边界层分离,总压恢复系数从0.502提高到0.56,喉道处流场畸变系数减小了18.6%.     

高超声速钝双楔绕流流动转捩与分离流动的壁温影响

为研究壁温对吸气式高超声速飞行器进气道转捩流动的影响,选取钝双楔这一典型外形,基于德国Aachen工业大学Thomas与Herbert所开展的双楔高超声速风洞试验,分析了一些已有的计算流体力学(CFD)研究内容,并结合本文不同方法的CFD数值模拟结果,讨论了不同壁面温度对该双楔模型高超声速绕流流动转捩与流动分离的影响。对于双楔模型,流动分离一般发生在拐角附近,由于流动分离旋涡的剪切作用会诱发流动转捩,转捩又会改变流动分离强度、分离涡尺寸,若分离流动存在非定常特征则将导致非定常旋涡运动与流动转捩的复杂相互作用。通过比较已有文献的CFD数值模拟结果与本文计算结果,表明只有按照转捩思路开展的数值模拟才能够反映该风洞试验情况。计算结果与试验数据的比较显示,文献中按照第一压缩面层流与第二压缩面湍流状态计算得到的结果能够在一定程度上与风洞试验数据相符,本文使用MUSCL格式、剪切应力输运(SST)湍流模型与γ-Reθ关联转捩模型这种计算方法,得到的结果与试验数据符合较好,正确地反映了风洞试验情况。分析还表明,在分离流动之前的区域,随着壁面温度的升高,壁面热流会下降,近壁区域黏性系数变大,边界层内速度剖面不饱满,速度边界层较厚,厚的速度边界层容易发生流动分离现象。     

解耦N-S/DSMC方法计算推力器真空羽流的边界条件研究

准确模拟羽流流场是合理评估空间推力器真空羽流效应的基础。因羽流流场同时包含连续和稀薄两种流动机理,通常采用解耦方式的Navier-Stokes(N-S)方程和直接模拟蒙特卡罗(DSMC)混合方法对其进行模拟。为保证解耦N-S/DSMC方法应用于羽流计算时的准确性并尽量提高其计算效率,对计算中的DSMC入口和喷管壁面等边界条件设置问题开展了研究。通过与文献中低总压工况的羽流试验数据比较,确定了合理设置DSMC入口边界位置、壁面反射类型、壁温等边界条件的方法。针对将此方法应用于实际推力器工况计算时效率过低的问题,通过多种数值试验,表明从喷管出口处开始沿KnGL为0.05的等值线作为DSMC入口界面可同时保证仿真精度和较高的计算效率;并证明实际工况下壁温设置对羽流场仿真结果影响不大。     

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

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

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

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

翼型模型几何误差对气动性能影响的自动微分分析方法

基于自动微分原理和NS方程有限体积方法建立了一套翼型敏感性导数计算方法和程序,可以一次性获得翼型不同气动力系数、压力分布对模型几何外形误差的敏感性导数和不确定度。计算结果表明,在跨声速范围内,即使翼型的外形误差只有63μm(弦长1m),也可以给翼型压力分布带来0.312(以来流动压为参考)量级的不确定度,而激波处的流动最为敏感。这种自动微分方法对于分析数值模拟结果的分散性、风洞试验结果的分散性或不确定性具有很好的指导意义。     

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

有翼高超声速再入飞行器是近年来的研究热点,气动设计是飞行器设计的关键。为了更清楚地认识有翼高超声速再入飞行器气动设计的难点问题,对有翼高超声速再入飞行器的发展、优势及总体任务剖面进行了介绍,从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.     

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

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

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

 波音公司面临着来自空客公司的巨大挑战,企业战略性创新才是公司成功的关键。为此波音公司的全部战略性研究集中在扩大产品的差异性上,体现在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.      

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

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

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

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

二维翼段颤振的μ控制

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

Experimental Study on Effects of Fuel Injection on Scramjet Combustor Performance

In order to investigate the effects of fuel injection distribution on the scramjet combustor performance, there are conducted three sets of test on a hydrocarbon fueled direct-connect scramjet test facility. The results of Test A, whose fuel injection is carried out with injectors located on the top-wall and the bottom-wall, show that the fuel injection with an appropriate close-front and centralized distribution would be of much help to optimize combustor performances. The results of Test B, whose fuel injection is performed at the optimal injection locations found in Test A, with a given equivalence ratio and different injection proportions for each injector, show that this injection mode is of little benefit to improve combustor performances. The results of Test C with a circumferential fuel injection distribution displaies the possibility of ameliorating combustor performance. By analyzing the effects of injection location parameters on combustor performances on the base of the data of Test C, it is clear that the injector location has strong coupled influences on combus- tor performances. In addition, an inner-force synthesis specific impulse is used to reduce the errors caused by the disturbance of fuel supply and working state of air heater while assessing combustor performances.     

Abnormal Shape Mould Winding

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