The Full Flowpath Analysis of a Hypersonic Vehicle

对一种类X43-A吸气式高超飞行器全流道开展了M7一级的三维数值仿真研究,深入探索了进气道处于起动状态和不起动状态时全流道的冷流流场结构和和气动力特性,且部分结果与实验数据进行了对比。研究结果表明:(1)前体横截面上存在显著的展向压强梯度,使得经过预压缩的气流偏离了进气道进口,但同时也减少了进入内通道的边界层气流,提高了进口流场的品质;(2)后体喷流股的膨胀过程受到了周围外流的显著干扰,因而沿流动方向其截面形状不断发生变化,如在喷口附近为近似矩形,而在后体末端附近则演化为近似三角形;(3)当进气道处于不起动状态时,其全流道流动结构发生了显著变化,进气道的外部压缩波系往复振荡,尾喷管出口的喷流股也在不停的膨胀和收缩,具有极强的非定常特征;(4)当进气道处于不起动状态时,全机的气动力特性呈周期性变化,升阻比的变化幅度较大,在最大、最小值分别可达起动状态的2倍和1/4倍。另外,升力、阻力系数的变化曲线之间存在一定的相位差;(5)与实验结果的对照表明,所采用的数值仿真方法具有较高的精度。     

喷流作用下的单边膨胀后体气动载荷

高隐身无人机为提高侧向、后向隐身指标,机身和喷管-后体结构的一体化设计尤为重要,该设计在气动、结构、强度等领域带来了新的挑战,喷流作用下的气动载荷预测是其中一个关键问题。利用CFD方法结合风洞试验数据,对高隐身无人机喷管-后体结构的稳态、动态气动载荷及其影响规律进行研究,结果表明：喷流形成的复杂波系投影到后体壁面上会形成压力、吸力交替的气动载荷,该载荷分布对喷管落压比极为敏感,小幅度的落压比变化即可能导致载荷方向变化;稳态气动载荷分布也会受到来流速度、次/主流流量比等的影响。基于IDDES方法的数值计算对动态气动载荷有不错的预测精度,将流场划分为喷流主流与次流掺混区、喷流核心区、喷流与外流掺混区等几个区域,有助于从流动机理上揭示脉动压力的分布规律,且得到了试验验证。动态载荷极值区位于喷流区边界,正是气流剪切掺混最强的位置,而喷流核心区尽管稳态气动载荷强,但动态气动载荷相对较低。     

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

CFD在螺旋桨飞机滑流影响研究中的应用

通过求解绝对坐标系下的非定常雷诺平均Navier-Stokes方程，计算了非定常滑流对某飞机巡航构型低速状态的气动特性影响。为描述螺旋桨桨叶的相对运动，采用了动态重叠网格技术，并在并行环境下采用全隐式双时间步方法和多重网格技术来保证时间精度和提高计算效率。将计算结果与试验结果相比较，在线性段，计算得到的升力系数与试验值总体偏差不超过2%，阻力系数相对试验值而言整体偏大，基本控制在8%左右，俯仰力矩系数的计算结果与试验值趋势一致。总结了滑流对螺旋桨飞机宏观气动特性的影响规律，由于滑流在空间发展过程中不断与周围空气相混合，使得滑流边界逐渐"模糊"，因此滑流流管在空间发展过程中的巨大畸变对滑流影响的分析工作带来了难度。针对这种现象，本文将有/无滑流的流场进行对比，通过当地动压增量来定义滑流的加速效应边界，以及通过当地气流角增量来定义滑流的洗流效应边界。该方法能较好地捕捉和解释由于滑流对飞机部件干扰而使得飞机方向安定性呈现的非线性现象，初步揭示了螺旋桨滑流复杂尾迹流动的特点，在螺旋桨飞机非定常滑流影响机理研究方面具有一定的参考价值。     

高超声速飞行器表面温度分布与气动热耦合数值研究

针对高超声速飞行器热防护设计中的高温气体非平衡效应问题和气动热环境精确预测问题,基于流场的非平衡Navier-Stokes方程、表面的能量守恒方程和内部的热传导方程,考虑流场的非平衡效应、表面的热辐射效应、催化效应和烧蚀效应以及热防护层内部的热传导效应,建立了初步的表面温度分布与气动热的耦合计算方法,完善了高超声速飞行器气动物理流场计算软件(AEROPH_Flow)。在表面材料为碳-碳(C-C)条件下,对飞行高度为65km和飞行速度为8,10km/s的半球以及飞行高度为50km和飞行速度为8km/s的球锥模型,开展了表面温度分布与气动热的耦合计算,验证了计算方法和计算软件,分析了表面温度分布对气动热环境的影响。研究结果表明:表面温度分布对气动热的计算结果有较大影响,在气动热环境的预测中,不仅要考虑热化学非平衡效应和表面催化效应的影响,还要考虑表面温度分布的影响,最好是采用表面温度分布与气动热耦合计算的方法,以减小表面温度分布对气动热计算结果的影响。为此,需要发展完善非平衡流场/表面催化和烧蚀/热传导温度场(气/表/固)的计算模型、耦合求解技术和计算软件,实现对高超声速飞行器的真实飞行条件下高温气体非平衡效应和气动热环境的精确模拟。     

考虑气穴影响的2D压力伺服阀稳定性

为研究双自由度（2D）压力伺服阀工作时的气穴现象对阀芯稳定性的影响，利用AMESim中的基本库和HCD库建立2D压力伺服阀各组件仿真模型。理论分析得出气穴现象会降低阀体中含气油液的有效体积弹性模量，得出考虑气穴影响的阀芯稳定性条件。仿真结果表明气穴现象会导致主阀芯到达稳定位置后出现振荡现象，影响主阀芯的穿越频率从21 Hz降低至6 Hz，限制主阀芯的频宽，适当增加主阀芯阻尼比能提高主阀芯工作稳定性。提出一种新型阻尼活塞结构，建立模型分析得出其能在不影响伺服阀频响、阀芯位移和抗污染能力情况下提高主阀芯阻尼比以提高阀芯工作稳定性，并且设计实验进行验证。实验结果表明，应用阻尼活塞将2D压力伺服阀压力输出波动从9%降低至2%，能够提高飞行器制动时刹车阀工作稳定性。     

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.     

航天器返回地球的气动特性综述

航天器返回地球的飞行过程中,气动特性是实现将宇宙飞行速度减到落地前速度、保证再入飞行得到有效控制以及再入防热安全可靠的关键因素。针对简单旋成体气动外形、半弹道式再入控制、烧蚀防热类返回航天器,综述了返回地球过程中变化的空气流域特性、航天器周围的气体绕流环境、空气与航天器作用产生的动力学与热效应等。系统地给出了该类航天器的再入气动特性参数与飞行性能的共性规律,包括:气动阻力与再入减速、气动升力与再入轨迹控制、配平攻角与飞行稳定性、气动加热与防热,以及再入过程中不同气动特性航天器、气象条件变化等对再入飞行性能的影响规律。为航天器开展返回飞行过程的跨流域气动性能工程研制提供设计参考。     

OPTIMAL CONTROL OF BRAKING MOMENT FOR A WHEEL AND TYRE

ＯＰＴＩＭＡＬＣＯＮＴＲＯＬＯＦＢＲＡＫＩＮＧＭＯＭＥＮＴＦＯＲＡＷＨＥＥＬＡＮＤＴＹＲＥＺｈａｎｇＬｉｎｇ（Ｘｉ＇ａｎＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｘｉ＇ａｎ，Ｃｈｉｎａ，７１００４９）ＺｈｕＤｅｐｅｉ；ＺｈａｎｇＹｕ（Ｎｏｒｔｈｗｅｓｔｅ...     

锥导乘波体构型的气动特性不确定度分析

为研究锥导乘波体偏离设计条件下气动特性变化情况,采用稀疏的非嵌入式混沌多项式方法,对乘波体气动特性进行了不确定性量化及全局非线性灵敏度分析。首先,采用CATIA二次开发技术对锥导乘波体进行参数化建模;其次,在来流速度、温度、密度和迎角满足特定扰动的条件下,通过拉丁超立方试验设计生成样本,并采用CFD进行计算;最后,根据试验设计样本建立响应面,通过混沌多项式分析得到了乘波体气动力系数的不确定度。灵敏度分析结果表明,迎角在锥导乘波体的气动特性变化中起主导作用。对马赫数和压强的流场不确定性分析结果表明,气动特性变化主因是乘波体前缘处的压力泄漏,影响了上表面压力分布,导致了气动性能的改变。