梯形截面微管道内流场的PIV测量

利用micro-PIV系统测量了去离子水在水力直径为237μm,长度为31mm的梯形截面微管道内的流场结构,得到了不同雷诺数下沿流动方向不同位置微管道中间截面上的速度分布,并利用数值模拟方法计算了相应雷诺数下三维梯形截面微管道内的速度分布。实验结果和数值模拟结果显示实验中的梯形截面微管道内从层流到湍流的转捩发生在Re=1500～1800左右,实验中充分发展段长度Le。微管道水力直径Dk和雷诺数之间的关系可由公式Le／Dk=（0．08～0．085）船表示。     

二次流动对气冷涡轮叶栅流场的影响

针对典型跨声速高压涡轮叶型平面叶栅吸力面单排孔气膜冷却,采用数值模拟方法,比较分析了加入气膜冷却前后流场变化。结果表明,由于二次流动的影响,加入气膜冷却以后吸力面后部接近下壁面处没有受到冷气保护而直接暴露于主流高温燃气,在实际高压涡轮中将极大的降低叶片寿命。没有气膜冷却情况下,吸力面接近下壁面处边界层仍有可能因受到二次流动的影响发生转捩;加入气膜冷却情况下,气膜孔中心位置下游边界层由于射流和主流的相互作用将转变为湍流边界层,而由于孔间距的影响,只有射流和主流充分掺混以后才能影响到整个叶片的范围。     

变循环发动机模态转换建模及控制规律设计方法研究

为了实现变循环发动机快速可靠的模态转换,本文发展了变循环发动机模态转换过渡态模型及控制规律设计方法。在变循环发动机动态数值仿真程序的基础上,针对模态选择阀与涵道引射器这两个模态转换过程中的关键变几何部件,建立了高精度的气流突扩局部损失模型。首次提出了可考虑模态选择阀堵塞对发动机性能影响的模态选择阀堵塞模型,消除了由于模型不精确造成的模态转换参数波动。在建立的变循环发动机数学模型的基础上,提出了基于直接推力控制技术的模态转换控制规律设计方法,考虑了变循环发动机的8个可调参数,采用差分进化算法对模态转换过程进行了优化。结果表明,本文提出的模态转换控制规律设计方法可以实现变循环发动机快速平稳的模态转换,双外涵转单外涵的参数变化规律与单外涵转双外涵的参数变化规律基本类似,推力在0.6秒就稳定在目标推力值,其余参数大多在1.4秒之后才趋于稳定。本文提出的变循环发动机模态转换控制规律设计方法还可应用于常规航空发动机的加/减速过渡态控制规律设计中。     

基于蒙特卡洛方法的冰晶撞击特性计算

航空发动机摄入冰晶易发生内部结冰、造成发动机推力损失,严重时引发航空事故。冰晶结冰数值模拟是研究冰晶结冰的重要手段,而撞击收集系数计算是冰晶结冰数值模拟的关键环节。基于NNWICE平台开展冰晶结冰前期的阻力系数、运动轨迹和撞击收集系数的数值模拟与计算研究,实现了基于蒙特卡洛方法的冰晶撞击收集系数的准确计算;在此基础上发展了拉格朗日框架下冰晶粒子运动、传热传质耦合数值计算方法;通过完成二维NACA0012翼型的相关算例,对比分析了不同阻力模型和冰晶相变对撞击收集系数的影响。结果表明阻力模型对冰晶撞击收集系数影响较小,冰晶相变过程会对冰晶运动轨迹及撞击收集系数产生显著影响,为进一步发展冰晶结冰数值模拟计算方法奠定了基础。     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

涡扇发动机低压部件通流耦合计算

为了研究各部件匹配与过渡段几何参数对发动机总体性能的影响,针对某涡扇发动机风扇、低压涡轮、内外涵道、掺混段和尾喷管等部件,基于周向平均Navier-St okes通流控制方程,采用时间推进有限体积法进行了部件耦合计算求解,将沿程各部件性能计算参数与试验做了对比,结果表明:在设计状态下发动机沿程各截面计算特征参数与试验值最大误差不超过5.0%,说明该通流耦合计算方法能够快速预估发动机总体性能,完善各部件匹配关系,在发动机总体设计初期有一定工程应用价值。     

一种带有过渡因子的下压指令设计方法

针对高超声速飞行器由滑翔段转入下压段的过程中运动状态无法快速与制导律进行匹配,从而引起弹道震荡的问题,设计了一种带有过渡因子的下压指令,改善了飞行器的过渡品质。首先,分析了高超声速飞行器平衡滑翔弹道,提出了平滑过渡的基本条件,并引入了一次函数、多项式、正弦、指数形式的过渡因子;然后,设计了3种带有过渡因子的下压指令;最后,通过仿真验证了下压指令的合理性。结果表明,该算法设计出的下压指令计算量小、制导精度高,且具有一定工程应用价值。     

临近空间环境下封闭方腔内耦合换热特性

以临近空间浮空器载荷舱为应用背景,对复杂热边界条件下含热源的三维封闭方腔内自然对流、表面辐射和导热的耦合问题进行了数值模拟。综合考虑对流换热、长波辐射、太阳辐射等因素的影响,建立了临近空间热环境模型。通过Fluent软件用户自定义函数（UDF）引入外部非定常的辐射-对流耦合热边界条件,对腔内换热特性的昼夜变化进行研究,并分析了腔壁厚度、发射率和导热系数对其的影响。数值结果表明,腔内平均温度昼夜变化很小,约为12.9 K,但温度场分布随太阳方位变化而变化;腔内对流换热较弱,同一时刻最大温差约为71.3 K;腔壁热阻和发射率增加会削弱自然对流的强度。      

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

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

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

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