发动机进气帽罩防冰热载荷的数值模拟研究

对某发动机进气帽罩防冰热载荷的计算方法进行研究。利用FLUENT软件对帽罩周围的流场进行模拟,同时采用FLUENT中的离散相模型计算帽罩周围的水滴运动轨迹,得到帽罩表面的水滴撞击效率;利用数值计算结果及流动换热经验关系式分析帽罩表面的各项热流,获得该发动机帽罩的防冰热载荷。计算中还比较了不同湍流模型对防冰热载荷的影响。     

霜状冰结冰数值模拟研究

提出了一种三维霜冰结冰的数值模拟方法。基于欧拉-拉格朗日法计算空气-水滴两相流流场,获得了表面水滴撞击特性以及撞击量,求解质量守恒和能量守恒方程获得网格单元内的结冰量,三维霜冰的结冰计算程序采用用户自定义（UDF）编写,FLUENT提供的动网格功能实现了结冰后的网格重构。计算了4°攻角,液态水含量为1 g/m3,水滴直径20μm,结冰时间6 min,温度分别为-19.4°和-28.3°时NACA0012翼型的结冰,并与相关文献结果进行了对比,计算结果与文献实验结果吻合较好,验证了方法的合理性和可靠性。分析了不同的粒径、液态水含量、来流速度等对撞击特性以及结冰的影响,随着粒径、来流速度的增加,水滴收集系数、撞击区和结冰厚度增大,随着液态水含量的增加,水滴收集系数和撞击区是不变的,但是结冰厚度增加。     

环境参数对飞机防冰热载荷的影响规律

为在飞机热气防冰系统设计过程中准确计算防护表面的防冰热载荷,以某飞机翼段模型为例,对防冰热载荷计算开展了研究,给出了计算流程,研究了来流速度、高度、温度、液态水含量和水滴直径等环境参数变化对翼段表面热载荷分布规律的影响。研究结果表明:①热载荷主要分布在机翼前缘,且越靠近翼尖部分热载荷数值越大,热载荷曲线呈现中间低两边高的分布;②表面热载荷分布范围和量级与来流速度和水滴直径呈正相关关系;飞行高度的变化对表面热载荷分布的影响较小;来流温度降低或液态水含量增加,热载荷数值增大但分布范围不变;水滴直径增加10 μm,上翼面热载荷分布范围增加约50%。      

CJ818客机机翼防冰范围的确定及防冰热载荷计算

通过数值计算的方法,对CJ818客机机翼周围的流场和水滴运动轨迹进行了计算;介绍了水滴撞击特性的计算方法,编写程序计算了机翼的水滴撞击特性,并据此确定出机翼的防冰范围;分析了防冰表面热流的组成,并给出了每种热流的计算方法,在此基础上编写程序计算得到了CJ818客机机翼的防冰热载荷。     

尾部可调的发动机进口支板结冰数值模拟研究

基于模块化思想,在FLUENT软件基础之上开发了发动机进口支板的结冰计算程序,针对尾部可调支板的结冰问题进行了计算研究。通过分析发动机进口支板结冰的物理过程,将支板结冰计算分成四个模块,分别为两相流流场计算、水滴收集计算、热力学分析和表面结冰量计算、动网格计算。首先将各模块单独进行研究,通过FLUENT提供的UDF函数及宏函数逐一实现各模块功能,最后将各模块整合,实现了发动机进口支板的结冰计算。通过计算获得了调整角度为30°时,支板的液态水含量分布、不同时间下支板结冰的区域以及不同温度下支板结冰的特点,给出了可调支板应重点防冰区域。数值模拟方法为工程应用提供了一条可选的途径。     

飞机水平尾翼水滴撞击特性及防冰热载荷计算

建立了利用欧拉法求解水滴撞击特性的方法,并基于Messinger质量和能量守恒方程建立了热载荷计算模型,以某型飞机三维水平尾翼为研究对象,展开了撞击特性和热载荷的计算.结果显示水平尾翼的局部水收集系数极值从翼根到翼尖逐渐增加,不同截面水滴收集系数的分布为设计防冰系统范围提供了依据.驻点附近及机翼上、下表面各存在一个热载荷较大的区域,进而确定该尾翼电加热防冰分为3个区域:中间连续加热区和上、下表面驻点附近的各一个间断加热区.      

二维机翼混合相结冰数值模拟

近年来冰晶引发的飞机结冰问题逐渐引起人们的重视。针对冰晶和过冷水滴同时存在的混合相结冰问题,通过数值模拟的方法实现结冰冰形预测。空气流场和对流换热的计算中采用了转捩剪切应力输运（SST）湍流模型,基于欧拉法计算冰晶和水滴收集系数。在Messinger结冰热力学模型基础上进行扩展,分析了二维结冰部件表面在混合相气象条件下的传热传质过程,建立了适用于混合相的结冰热力学模型,同时考虑冰晶的黏附效应,添加黏附效率经验公式。利用FLUENT的用户自定义函数（UDF）编程求解混合相热力学模型,计算了霜冰和明冰条件下NACA0012翼型表面结冰情况,与国外风洞试验结果进行比较,验证了计算模型和方法的有效性。结果表明,冰晶黏附效应对混合相结冰量及冰形有很大影响,明冰条件结冰形状偏向楔形。     

基于流-固耦合传热的热气防冰系统干空气飞行蒙皮温度场计算研究

以机翼热气防冰系统为研究对象,建立了包含热气防冰系统防冰腔内外流场对流换热和固体结构导热的三维稳态流-固耦合传热物理模型,对整个计算区域生成混合网格,边界条件为第三类边界条件,采用计算流体力学方法以 FLUENT 软件为工具,对干空气飞行状态下流-固耦合传热模型进行了求解,获得防冰腔蒙皮内外表面对流换热系数分布和温度场结果,并对计算结果进行了分析。结果表明：防冰腔铝合金蒙皮沿展向和厚度方向导热显著,温度分布较均匀,防冰引气温度为200℃时,防冰腔蒙皮内外表面上最高温度为101℃,最低温度为21℃,3 mm厚的蒙皮同一点处内外表面最大温差仅为4℃,防冰腔排气口处气体的平均温度为63℃。热气防冰系统蒙皮温度场计算方法和计算结果,能够为热气防冰系统干空气飞行试验设计和测试中温度传感器的选型与布置提供依据。     

发动机进口支板防冰热负荷数值模拟

为研究航空发动机进口防冰问题,应用三维雷诺平均N—S方程和双方程k～ε模型,计算了某发动机进口支板外部空气流场,用欧拉法建立并求解了空气中过冷水滴控制方程;选取某些飞行工况,得到了支板防冰表面的水滴撞击特性,且分别在干、湿蒸发2种情况下得到了表面的防冰热负荷分布。可以此为基础进行支板表面加热形式和结构设计.     

某型飞机发动机短舱防冰系统设计计算

介绍了某型飞机发动机短舱防冰系统的结构形式和防冰原理,按照可利用热量应大于所需热量的原则,通过计算完全蒸发和湿状态下的热载荷,确定了发动机短舱防冰严酷状态及所需热流量.基于机翼水收集系数和数据相似原则,利用欧拉法两相流计算方法,确定了短舱防冰系统水收集系数.按照在湿状态下允许有一定结冰量存在的设计原则,根据发动机的吞冰能力,利用VSAERO/ICE工具对发动机短舱防冰系统进行了计算,计算结果表明,在严酷状态下,防冰系统能够满足允许结冰厚度要求.      

燃烧室壁冲击-逆向对流-气膜冷却特性的数值研究

采用FLUENT软件对冲击-逆向对流-气膜冷却结构进行了流固耦合计算, 得到了对流通道和气膜形成区流场、流体温度场以及固壁温度场分布, 分析了该冷却结构内部复杂的流动和换热情况.计算研究目的是摸索一种火焰筒壁温的计算方法, 以较好地获得火焰筒壁温分布.计算结果和试验结果进行了对比, 两者有较好的吻合性.分析表明采用流固耦合计算获取冲击-逆向对流-气膜冷却结构的壁温分布是可行的.      

小型平板回路热管蒸发器内传热与流动的三维数值模拟

针对小型平板回路热管蒸发器内的流动与传热，建立了多区域耦合的数学物理模型，并应用FLUENT软件进行了三维数值模拟。结果表明:蒸发器传热特性在不同热负荷下呈现出较大的差异，其温度分布不仅取决于热负荷，更依赖于毛细芯表面发生的两种传热机制，即毛细蒸发和热传导。相比高热负荷(Q=120W)和低热负荷(Q=40W)，中等热负荷(Q=80W)下蒸发器各个部位的温度均较低。三种不同热负荷下，毛细芯反向导热均大于侧壁漏热，补偿腔内与毛细芯相邻处易出现高温区。冷凝回流液在补偿腔内的流动形成两个涡，这种流动特点有利于降低毛细芯的温度。当热负荷与系统冷凝能力匹配时，整个系统流动与传热特性最优。      

具有发散冷却功能的曲面结构边界层特性实验和数值研究

用实验测量和理论计算相结合的方法,研究了具有发散冷却功能的曲面结构边界层流动和换热特性.半球形实验件由高温合金钢粉末烧结的多孔材料制成,孔隙率和特征尺寸分别为0.37和20μm.实验件外部的热流温度及速度由热风洞的控制系统给定,内部的冷气流量由数字式质量流量控制器给定,外表面温度由NEC TH5104远红外热像系统测量.利用FLUENT软件中的二维、稳态、局部热平衡模型,对实验件进行发散冷却特性数值模拟,并通过与实验数据的比较,验证数学模型和数值方法.      

燃气轮机环形燃烧室燃烧流动的数值研究

由于燃气轮机燃烧室内复杂的物理化学变化，利用数学模拟的方法来研究，对减小燃烧室研制费用，缩短研制周期具有重要意义。对QD128型燃气轮机燃烧室燃烧天然气进行了数值模拟，在模拟过程中采用了雷诺应力模型、EBU—Arrheniue湍流燃烧模型和六通量辐射模型来描述其燃烧流动过程，运用FLUENT软件求解了三维流场和温场分布。计算结果能够很好地反映环形燃烧室燃烧流动特点，对预测环形燃烧室内的燃烧流动有一定参考价值。     

基于螺纹锥形搅拌头的搅拌摩擦焊三维流场仿真研究

本文利用CFD软件包FLUENT建立了带螺纹表面的锥形搅拌头搅拌摩擦焊三维流动模型。通过数值模拟得到了焊缝沿厚度方向的水平切面速度矢量图,截面的速度矢量图以及用于观察材料流动行为的流动线图。同时分析了搅拌头表面螺纹对塑化金属流动行为的影响及焊缝塑化金属的流动规律。模拟结果对于搅拌摩擦焊过程中焊缝处塑化金属流动以及焊缝形成机理研究具有重要的指导意义。     

整体式粒子分离器数值模拟

 为找出分离效率更高、压力损失更小的最优型面设计，以两种型面差异较大的整体式惯性粒子分离器为研究对象，利用FLUENT计算软件对其内二维两相流动进行了数值模拟，得到了分离器对不同粒径的砂粒的分离效率，以及不同清除比对分离效率的影响，并对计算结果进行了比较分析。在此基础上，分析了型面对分离效率的影响因素。计算采用标准k-ε模型，颗粒相采用拉格朗日轨道模型。结果表明：对C规范砂，型面变化剧烈的模型（IPS01）的分离效率比型面变化平缓的模型（IPS02）的平均低1.60%，进气总压损失平均增大0.53%；对AC粗砂，分离效率前者比后者平均高0.25%，进气总压损失平均增大0.58%。相比较而言，IPS02更具实际应用价值。     

单级跨音压气机整圈三维动静叶干涉的数值模拟

本文使用商业软件FLUENT,对某单级跨音压气机的整圈三维非定常动静叶干涉流场进行了数值模拟,计算并比较了拉开动静叶轴向间隙和正常轴向间隙情况下的不同流动状况。非定常整圈三维计算的结果模拟了清晰的动静叶干涉图像,并显示了静叶三个连续通道出口流量的周期性变化和滞后现象,表明可以利用成熟的CFD商业软件研究叶轮机械内部的复杂物理现象。     

二维非对称喷管数值模拟与验证

采用Fluent软件针对二维非对称喷管进行数值模拟,考察了其中几种常用湍流模型及网格密度对模拟精度的影响,并将计算结果与NASA公布的实验数据进行比较,找出了精度较高、求解快捷的模拟方案,结果表明FLUENT软件可以应用于带高超声速外流的非对称超声速喷管内外流场的数值计算。      

稳态热传导下基于多相材料的一体化设计

提出一种基于多相多孔材料的一体化模型。模型以结构散热弱度最小为目标,以结构体积比与材料微结构质量为约束。引入独立的宏观设计变量和微观相设计变量,通过单元相密度建立两类变量间的联系;基于对等混合材料插值模型建立单元相密度与热传导系数间的惩罚关系;推导得到目标函数灵敏度表达式。求解偏微分方程实现单元散热弱度过滤,消除棋盘格及网格依赖性现象。通过二维数值算例讨论并分析了材料特性、热载荷、体积比约束以及质量约束对一体化优化结果的影响。数值实验结果表明,该建模方法在多相材料/结构一体化稳态热传导优化设计中具有可行性和有效性。