润滑系统管路流动与换热耦合计算研究

为研究沿程管路换热与滑油流动之间的耦合问题,通过分析主要元件的流动和换热特性,引入元件换热的能量方程,结合管路流动的连续性方程和动量方程,建立了润滑系统管路热流耦合计算模型。针对一段发动机供油管路,利用MATLAB/Simulink建模仿真工具包,建立其热流耦合网络计算模型,获得滑油流量、温度以及管壁温度等参数,采用商业软件Flowmaster对其进行验证。模型计算结果与Flowmaster计算结果非常接近,滑油流量和温度的最大误差均不超过3%,管道壁温的计算误差不超过2.4%,说明建立的热流耦合计算模型具有良好的仿真精度,适用于考虑沿程管路换热的发动机润滑系统计算分析。     

复燃对液氧煤油发动机尾焰冲击特性影响研究

为了研究复燃对液氧煤油发动机尾焰冲击特性的影响,建立了液氧煤油发动机尾焰冲击数值计算模型,并基于模型研究了喷管出口距离平板3m,5m两种工况下复燃对尾焰冲击特性的影响。结果表明:模型考虑了发动机内部燃烧对尾焰冲击特性的影响,计算得到了主射流区的激波结构;复燃增大了尾焰自由射流区和壁面射流区的高温区域,改变了自由射流区和滞止区的形状结构;平板壁面压力随着径向距离增大而逐渐减小,并且3m工况时在1.8m和2.5m处分别出现2.5倍环境压力和1.5倍环境压力的波动,5m工况时在2m处出现1.5~2倍环境压力的波动,在波动之后平板壁面上压力很快降为环境压力,复燃对5m工况的波动较3m工况影响大。     

幂律型流体撞击射流破碎特性直接数值模拟

凝胶推进剂的双股射流撞击雾化广泛应用于液体火箭发动机的燃烧室中,其破碎特征及雾化效果直接影响燃烧效率。为探究雾化特性的发展规律,采用直接数值模拟DNS方法,对射流速度为100m/s的剪切稀化非牛顿液体正交撞击产生的雾化特征、液体表面积、表面波、涡特性以及非牛顿特性开展研究。结果表明,射流下形成的雾化流场迅速扩张形成液膜,液膜两侧边缘破碎成大量的液丝与液滴,核心部分产生撞击波后在气体力的作用下逐步发展为带有凸起和褶皱的不稳定表面波,其撞击波波长最大可达2.46倍射流直径。液体表面积不断增长,但无量纲表面积总体呈现先下降再上升的趋势。气体中的涡量分布则分为有序附着区和无序爆炸区两类,并且涡量主要集中分布于气相区域。此外,射流撞击时产生强剪切使该液体内部的粘性系数下降,最低仅为初始粘性系数的0.3倍。     

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

低压涡轮动叶叶顶间隙泄漏流动特性数值计算

为了获得低压涡轮动叶叶顶间隙泄漏流动特性,采用数值方法研究带冠动叶叶顶间隙对低压涡轮级气动性能的影响。 对试验测量的扇形带冠叶顶叶栅气动性能进行数值计算,得到的叶栅总压恢复系数与试验数据吻合良好,验证了数值方法的可靠 性。对没有间隙和3 种动叶叶顶间隙条件下低压涡轮级的气动性能和动叶叶顶间隙泄漏流动特性进行对比分析。结果表明：带冠动 叶叶顶间隙增大导致转子叶尖的出口气流角增大和叶顶间隙泄漏量增加;涡轮级气动效率降低和叶顶间隙泄漏量增加与带冠动叶 叶顶间隙增大呈近似线性变化;在相同叶顶间隙条件下,出口马赫数对叶顶间隙泄漏量的影响较小;随着涡轮级出口马赫数的提 高,出口气流角度减小。     

基于欧拉方法的2 维翼型冰晶结冰数值计算

针对现有航空发动机冰晶结冰试验方法存在数据采集困难,且难以揭示冰晶结冰机理问题,对某楔形翼型进行冰晶结 冰数值计算。采用S-A 湍流模型获得流场结果,应用欧拉方法获得冰晶和液滴轨迹结果,基于Messinger 模型获得冰形,并以 NASA-NRC 的第139 号试验结果验证了数值方法的可行性。结果表明：在其他工况相同时,压力越低（飞行高度越高）,冰晶/ 液滴蒸 发越强烈,结冰情况越严重;分析不同液态水含量与总水含量比例（CLW/CTW）对结冰的影响可知,液态水过少或过多都不利于冰层形成。     

壁面催化效应对高超声速气动热影响研究

高超声速飞行器壁面催化效应会导致激波层中原子在壁面处复合释热,加剧周围气动热环境。针对高超声速流动壁面催化特性,选择不同飞行马赫数及高度条件,采用完全催化和非催化两种条件对球锥模型壁面热流率进行数值模拟计算,研究壁面催化效应对气动热的影响规律。结果表明,固定飞行高度时,壁面催化效应对气动热的影响随马赫数增加而加强,Ma=25条件下驻点处完全催化与非催化热流比值高达1.92;固定飞行马赫数时,在50km高空以上壁面催化效应对气动热的影响随高度增加而减弱;壁面催化效应不仅会影响壁面附近的流场特性及组分分布状态,而且对整个激波层都有一定的影响作用。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。