机翼电加热防冰表面内外传热的耦合计算

分析了机翼防冰表面的传热传质现象,对电加热防冰表面内外传热耦合计算进行了研究.采用附面层积分法对机翼表面传热系数进行计算.根据传热传质比拟得到了防冰表面蒸发量.建立表面微元的能量与质量守恒差分方程组,利用高斯-赛德尔迭代法求解,将外部防冰载荷、蒙皮弦向导热和内部加热热流三者耦合,采用亚松弛系数稳定迭代过程,最后得到平衡表面温度、溢流水和溢流结冰的范围.计算结果和文献中试验以及其他数值仿真结果比较,吻合较好,表明了模型及计算方法的正确性.      

基于内外传热耦合的热气防冰系统仿真计算

基于内外传热耦合原理,建立了热气防冰系统的AMESim仿真模型,研究热气防冰系统在不同引气状态下,管路流量、压力及蒙皮温度的变化规律.采用所建立的仿真方法,计算飞机机翼热气防冰系统的内部流动特性和内外耦合传热特性,将计算结果与热气防冰系统流量分配试验结果进行对比.结果表明:管路流量、压力和蒙皮温度的仿真计算结果与试验结果最大误差为9%,验证了该仿真模型的正确性.在此基础上,对飞机短舱热气防冰系统进行了仿真,分析了飞行包线下系统内外的换热、温度变化、加热效率等关键参数的瞬态特性.仿真结果为热气防冰系统的设计、分析与优化提供依据.      

液态水含量对防冰表面水膜流动换热的影响

为了研究液态水含量对防冰表面水膜流动换热的影响,基于机翼防冰表面水膜及空气相互作用机理,并考虑水膜表面传热传质过程,建立了水膜与空气的流动换热模型,得到溢流水膜及空气边界层流动换热的积分控制方程,通过对比文献试验结果验证了模型的准确性.在此基础上,比较了不同液态水含量条件下防冰表面水膜厚度及主要热流量的分布情况.结果表明:液态水含量对水膜沿表面厚度分布有明显影响,而对换热过程中各项热流的影响主要集中在水滴撞击区域,加热热流与散热热流随液态水含量的增加呈现相反的变化趋势.     

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

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

迎风面三维积冰过程中水膜流动的计算方法

对影响迎风面三维积冰过程的冰层表面水膜流动进行分析,建立描述其流动的数学模型,引入在流动前锋下游预设薄水膜的方法处理其在表面干区的流动推进过程,提出超过临界厚度的水膜被吹除的方法处理角状冰下游局部水膜被吹离表面的现象.通过对翼型-平板结构中翼型上的冰形对比和水膜流动规律分析,验证了该水膜流动数学模型的合理性和计算方法的可行性.研究表明:与Messinger模型相比,积冰模型中引入水膜流动后,可以较好地模拟三维明冰积冰,同时明冰极限的模拟结果也更接近试验值;冰层表面的水膜厚度约为10-5m,流动速度约为10-2m/s.      

基于传热与压力损失的气-水换热器耦合设计

针对换热器设计中如何合理兼顾传热效率和压力损失的难题,以矩形流通截面的燃气-水列管式换热器为对象,从理论上分析了设计参数之间的内在联系,揭示了换热面积与燃气流通面积、横向管间距、换热管直径之间的匹配关系,建立了基于压力损失和传热计算的耦合设计方法。这种方法已成功应用于某大型燃气冷却器的设计,也可供其它换热器设计时参考。     

来流速度对防冰表面溢流水流动换热的影响

为研究来流速度对防冰表面溢流水流动形态及换热的影响,基于空气-水两层相互作用的质量、动量和能量守恒,建立防冰表面溢流水水膜流动换热及破裂的数学模型,分析了防冰表面溢流水在不同来流条件下的流动形态和表面换热情况.计算分析表明:来流速度增加时,防冰表面相同位置处的连续水膜厚度减小,水膜破裂位置随之延后;较高来流速度条件下,破裂处水膜厚度稍有增加,使得破裂后形成的溪流厚度和宽度增大;作为主要的表面散热项,连续水膜表面蒸发及对流换热热流均随来流速度的增加而增大.此外,由水膜破裂引起的表面溢流水流态变化对防冰表面蒸发热流有一定影响.     

基于Lagrange液膜方法的小推力液体火箭发动机耦合传热策略研究

为了研究小推力液体火箭发动机边区冷却液膜的冷却效果,提出了一种综合考虑燃烧流动、气体辐射、液膜冷却、壁面传热和辐射冷却的耦合传热计算方法,构建了基于Lagrange方法的液膜计算模型。采用该耦合传热方法对490N液体火箭发动机进行了数值仿真研究。结果表明:液膜最大厚度60μm,最长30mm,覆盖区域的壁面温度都被限制在420K左右,起到了较好的冷却作用。喷注角度较小射流形成的液膜更集中,生存时间更长,可减小向壁面的热流;喷注角度较大射流形成的液膜周向范围更大,能缓解壁面热流向头部的传递。     

基于TLBM-FVM耦合的飞行器舱内热环境跨尺度预测方法

舱内热环境的有效预测是优化飞行器热控与防热设计、减小系统冗余并保障飞行器热安全的重要基础。在国家数值风洞（NNW）工程支持下,针对目前舱内热环境多尺度、多机制复合传热特点及其数值预测面临的精度与效率提升难题,发展了多区域协同推进的时空耦合模型及流/固界面的自适应分辨率识别算法,建立了基于热格子Boltzmann方法（TLBM）与有限体积法（FVM）相互耦合的舱内复合传热跨尺度预测方法,开展了典型飞行器仪器舱的综合热分析,验证了耦合方法的计算精度及效率。研究表明,相关方法可实现舱内热环境的局部精细化与整体大规模的协同模拟,用于开展整体自然对流与设备局部热量传递的多尺度数值模拟,掌握不同环境参数对舱内热质传递过程的影响规律,从而为飞行器热防护/热管理一体化设计提供重要技术支撑。     

颗粒无序分布的堆积床内部流动与传热分析

为了详细分析单组元发动机催化床内部的流动与传热,针对颗粒无序分布的堆积床,提出了一种网格生成途径,适合于对大量颗粒自然堆积的管道流动进行数值模拟。根据该方法,对不同颗粒直径、不同进口速度的管道堆积床进行了液体介质的流动与非定常传热模拟,并与实验结果和Ergun方程计算进行了对比。结果表明,本文提出的方法可以有效模拟颗粒无序分布的堆积床内部流动与传热,适合于大量颗粒的堆积床仿真;网格尺度小于1/20颗粒直径时,模拟结果与实验结果符合很好;在相同雷诺数下,摩擦系数随颗粒直径增大而减小,Ergun方程在高雷诺数下计算的流阻偏大;非定常传热时,壁面效应会使壁面附近温度先接近流体温度。     

一维和二维电热除冰相变传热特性的参数影响分析

基于飞机电热除冰过程冰层融化的特征,开展了一维和二维电热除冰相变传热特性数值计算研究和参数影响分析,重点考虑了加热模式、冷却时间、加热功率和加热单元间隔等参数对冰层相变传热的影响。采用了基于焓-多孔介质方法的热焓模型,将计算区域看作是包括多层材料和冰、水及其混合区的多孔介质,采用了结构化网格拓扑结构对计算区域进行划分,采用了有限体积方法对控制方程组进行离散,采用线性插值的方法获得混合区的物性参数,耦合能量方程和液态水体积分数公式,迭代求解了计算域的温度分布,获得了不同材料间界面温度的变化,重点分析了冰-保护层界面的温度变化。基于二维电热除冰模型的冰-保护层界面温度不均匀特征,提出了耦合考虑冰-保护层界面热点和冷点温度的冰脱落温度判断准则。研究表明:高功率的周期性加热模式要优于低功率的连续性加热模式,采用合理的冷却时间和加热功率,可获得更低的能量消耗和更好的除冰效果。合理设置加热单元间隔可以提高周期性的除冰效率,但也会形成冷点和热点,造成冰-保护层界面温度分布的不均匀。冰-保护层界面的冷点类似于锚点,即使此刻热点的冰已经融化,整个冰层也无法脱落和剥离。因此,冰脱落的判断要耦合考虑界面热点和冷点的温度特征。      

直升机旋翼桨叶防/除冰系统防护范围研究

直升机主旋翼桨叶在结冰气象条件下有结冰的危险, 为了满足直升机全天候仪表飞行的要求, 必须对桨叶采取结冰防护.本文针对某翼型直升机桨叶, 通过分析表面各项热流, 计算表面温度为零时极限结冰状态的临界速度, 确定了桨叶展向结冰防护范围为沿展向从翼根到翼尖0～95%;采用分析水滴在气流场中运动的拉格朗日法计算水滴撞击特性, 确定桨叶弦向结冰防护范围为上表面12%, 下表面33%.在此基础上对桨叶表面加热形式和结构进行设计, 可为直升机旋翼桨叶防/除冰系统设计奠定基础, 此研究具有重要的实用价值.      

冰脱落对电热除冰传热特性的影响研究

冰脱落对电热除冰传热特性的影响是电热除冰过程中的重要特征.采用有限体积法和混合网格技术,基于热焓方法将传统的相变传热模型应用到电热除冰计算中,研究了冰脱落对电热除冰传热特性的影响,主要研究了加热模式和加热单元间隔对除冰结构传热特性的影响.研究发现:无间隔除冰单元是有间隔除冰单元的特殊情况,其蒙皮表面温度历程在本质上是一致的.冰脱落后,无间隔除冰单元蒙皮表面温度开始剧烈上升然后逐渐缓慢上升.而对于有间隔除冰单元,随着间隔上加热功率的降低,温升过程逐渐缓和.单元间隔的加热功率降低到一定程度时,表面温度出现温降过程.当传递到蒙皮表面的热量和对流换热带走的热量差逐渐趋于稳定时,表面温度又开始逐渐稳定的上升.     

多股自耦合射流相干流动和换热特性数值研究

采用数值模拟的方法对多股自耦合射流在相干条件下的流动和换热特性进行了研究.研究结果表明:三股自耦合射流都在出口附近耦合成一股射流向下游运动,两侧的射流受到中间一股射流的卷吸融合在一起,融合后仍然是一股沿中心线对称的自耦合射流.在短轴中心截面自耦合射流型面较宽,射流喷出后在下游的扩展角度大,而长轴中心截面射流的型面较窄,扩展角度较小.在冲击靶板上,整个冲击区域呈现出中心温度最低的近似圆形的对称性分布,对流换热系数曲线呈现出中心对称的波峰分布形状,表明自耦合射流驻点区换热最强.对流换热系数沿靶板不同方向的分布趋势基本保持一致.      

克努森数对气膜冷却流动换热相似特性的影响

数值模拟研究克努森数Kn对不同尺度气膜冷却流动换热相似特性的影响,选取孔径为0.3,3,10mm平板单排圆柱形气膜孔模型,计算得到模型在Kn不同与相同条件下冷却效率和孔流量系数.结果表明:Kn相同时,小孔与大孔模型流动换热特性相似,冷却效率及流量系数均吻合很好;Kn不同时,两模型流动换热特性存在一定差异,冷却效率及流量系数吻合较差.因此为在气膜冷却相似放大研究中得到准确的相似结果,需要保证小孔及大孔模型的Kn相同.