展向离散抽吸法控制边界层转捩实验研究

通过实验成功发展了一种生成大振幅稳定条带的有效方法,即展向离散抽吸方法。在此基础上,开展条带控制边界层转捩的研究。实验在水洞中进行,以零压力梯度平板边界层为研究对象,利用氢气泡时间线法观测引入条带前后人工激发转捩边界层中扰动发展变化,分析条带对转捩的控制效果及参数影响。结果表明,展向离散抽吸方法生成的稳定条带最大振幅可达28.4%U（U为自由流速度）;实验中引入的宽度14和28mm的稳定条带都能起到抑制转捩的作用;条带振幅越大、宽度越窄,抑制效果越明显。研究结果为探索降低水/空气中航行器摩阻的新技术提供有价值的参考。     

高超声速飞行器-进气道一体化热流数值计算

采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.      

PSE方法在后部加载叶栅绕流稳定性分析及转捩预测中的应用

采用数值求解与实验测量相结合的方法, 分析后加载叶栅绕流的稳定性并预测转捩.推导出正交曲线坐标系下三维扰动波的抛物化稳定性方程, 在风洞实验中测量叶栅叶片表面的静压分布, 并结合测量结果对设计叶栅的流动稳定性进行特征值分析.分析结果表明:(1)后加载叶型能有效推迟转捩, 提高叶栅的变冲角特性;(2)当冲角变化时叶顶流动稳定性的变化最显著, 负冲角流动稳定, 正冲角流动最不稳定.      

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。     

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

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

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

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

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

采用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℃）。