混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

高压燃气涡轮径向内冷叶片气热耦合的数值分析

采用气热耦合方法对采用径向内冷方式的MarkⅡ型跨声速高压燃气涡轮金属导叶进行数值模拟,通过分析叶片通道内的传热和流动过程发现叶片表面附面层内流动非常复杂,包含层流流动、转捩和湍流流动状态,所以只有使用转捩模型计算的叶片附面层内流动与实际情况相符,叶片壁面温度和换热系数分布与实验结果吻合的较好,使用其他湍流模型由于不能准确描述附面层内流动而使得计算结果误差相对较大,但是所有的湍流模型都能很好的模拟附面层以外流动.      

CE/SE法模拟爆震波点火及传播过程

运用时空守恒元和解元(CE/SE)法模拟一维爆震波的形成及传播过程,分析产生爆震波的不同机理和DDT(爆燃-爆震转捩)过程的影响因素.采取激波点火和热点点火两种点火方式,分别对应SDT(激波-爆震转捩)和DDT过程,热点点火初始场采用线性温度分布和阶梯性温度分布,整个管道均为常压.用隐式梯形方法实现刚性源项积分.结果表明,SDT和DDT对应于不同的爆震波机理,并且DDT受点火区不同参数的影响,为研究爆震波点火和DDT过程提供了理论依据.     

基于零控脱靶量的大气层外拦截中制导研究

针对大气层外拦截中段采用耗尽关机固体燃料发动机的情况,设计了基于零控脱靶量的鲁棒中制导律。推导了系统状态转移矩阵,通过数值方法求解得到考虑地球扁率J2摄动项的拦截器零控脱靶状态,在此基础上将推力分解到视线方向和垂直视线方向,并根据垂直视线方向需要速度增量大小设计了推力分配比例系数,使得推力方向平稳变化且能快速减小系统的零控脱靶量。仿真表明,所设计的中制导律简单、可靠,便于工程实现,能有效提高拦截精度,并对发动机工作参数具有较强的鲁棒性。     

前体边界层状态对高超声速进气道流动结构及性能的影响

采用CFD方法研究了前体边界层状态对轴对称高超声速进气道的波系结构、激波/边界层干扰特性以及流量捕获能力、总压恢复系数等的影响,获得了不同马赫数下进气道性能参数随一级锥边界层强制转捩位置的变化规律,并对采用二维粗糙体作为强制转捩措施给流场所带来的干扰进行了分析,研究中还对比了具体进气道设计方案间的差别.      

雷诺数变化对翼型边界层发展及失速特性的影响

翼型失速及其边界层发展是飞行器设计中的基础科学问题,而雷诺数变化对其影响很大。针对后缘失速翼型,采用Menter k-ω SST模型及耦合扰动放大因子输运方程的转捩模型,进行雷诺数变化对层流-湍流转捩边界层特性和失速特性的影响分析。结果表明:雷诺数增大时,对于转捩边界层,当地涡量雷诺数增大,转捩前移且分离泡减小,流动能量耗散减小,翼型整体表面剪切效应增强,动能更充沛,流动自持能力增强,压力分布可以维持较长距离的梯度抵抗分离能力增强;因此雷诺数增大使翼型失速迎角提高、升力系数增加。     

基于双eN方法的翼身组合体流动转捩自动判断

发展翼身组合体复杂外形流动转捩自动判断方法,对高亚声速民机自然层流（NLF）机翼设计具有重要意义。使用多块结构化网格和三维雷诺平均Navier-Stokes（RANS）方程求解器,耦合边界层方程求解和基于线性稳定性理论（LST）的完全双e^N方法,发展了一套可同时计及Tollmien-Schlichting波和横流不稳定性扰动诱导转捩的翼身组合体流动转捩自动判断方法。对DLR-F4翼身组合体绕流进行了转捩自动判断,将得到的转捩位置与试验结果进行比较,验证了所发展方法的正确性。使用上述方法对配置自然层流机翼的中短程民机翼身组合体外形进行了数值模拟,并将结果与单独机翼的转捩位置进行了对比,结果表明机身三维位移效应增强了自然层流后掠机翼边界层的横流不稳定性强度,导致翼根转捩位置提前至前缘区。     

结冰风洞中液滴相变效应数值模拟

为明晰结冰风洞中液滴相变效应,发展了基于Euler法的气液两相传质传热耦合流动计算方法,探索了3 m×2 m结冰风洞主试验段构型相变效应对液滴传热过程的影响,开展了参数影响研究,评估了试验段内液滴过冷状态。结果表明：构型内液滴经历了先蒸发后凝结两个阶段,蒸发效应促进了准一维传热阶段液滴温度的下降趋势,使液滴温度趋于湿球温度,而凝结效应则抑制了三维收缩阶段液滴温度的下降趋势,进而增大了试验段液气温差,影响液滴过冷状态;增大初始相对湿度和试验段气流速度,会导致蒸发效应减弱而凝结效应增强,进而增强了相对湿度和气流速度对液滴过冷状态的影响程度,与此相反,初始液滴尺寸的增加,则会导致蒸发效应和凝结效应均减弱,进而减弱了液滴尺寸对液滴过冷状态的影响程度;在典型工况下,小尺寸液滴（液滴直径范围为40~100 μm）在高风速时（试验段气流速度大于164 m/s）将偏离过冷状态（液气温差超过3℃）。     

考虑维修力量影响及载荷动态分配的k/n系统模糊可靠性分析

工程实际中,维修活动开展前往往存在一定时长的准备期,且由于环境时变性、系统长期运行后的劣化累积等因素导致部件状态性能水平存在不确定性,使得系统可靠性建模较为困难。对此,运用模糊数表征系统部件的失效转移率、修复转移率及修理工维修准备率的同时,以Power Law规则刻画部件间的故障相关关系,认为部件承担载荷超过某阈值时才会引发故障相关现象,并考虑了修理工数量与故障件数量之间关系对系统可靠性的影响,研究了载荷动态分配条件下带维修准备期的多修理工n中取k模糊多状态系统模型,建立了状态转移微分方程组,提出用逆向逐层分析的思路建立系统稳态概率系数的递推关系,应用α水平截集及Zadeh扩张原理确定了模糊状态概率的截集区间,得到了系统模糊稳态指标,最后通过算例给出了修理工数量及部件参数模糊程度对系统稳态指标的影响,验证了模型的适用性。     

基于混合卡尔曼滤波器组故障传感器定位方法

为了在发动机性能蜕化与传感器故障并存的情况下实现故障传感器的定位与部件蜕化情况的估计,并实现故障诊断基准数据的修正,构建了1种包含了机载模型与线性卡尔曼滤波器的组合结构混合卡尔曼滤波器组。该卡尔曼滤波器组能够在之前所描述的故障/蜕化耦合情况下定位故障传感器,并得到较为准确的部件蜕化估计结果。为了验证了混合卡尔曼滤波器组的有效性,进行了相关仿真。仿真结果表明:混合卡尔曼滤波器组能够在发动机动态过程中遭遇传感器故障与部件蜕化并存的情况下完成故障定位与蜕化估计。