机翼热气防冰系统设计

利用热流体系统仿真分析平台Flowmaster和计算流体力学软件Fluent进行联合仿真,对机翼热气防冰系统进行了设计和性能分析.设计参数包括防冰系统笛形管热气喷口孔径、孔数及孔间距等参数.基于Flowmaster对所设计的供气管路进行供气流量分配计算.基于Fluent模拟防冰腔内流动与换热,计算了防冰腔的热效率.根据防冰热载荷及热效率结果预测每段缝翼防冰所需热气流量,引入参数供气余度来分析防冰系统热性能.所提出的热气防冰系统设计方法可为热气防冰系统的工程设计与优化提供重要的帮助.     

对流换热系数的瞬态测量技术

使用等热流法瞬态实验测量对流换热系数:基于流动稳定及热平衡条件下对测量表面突加一稳定热流,测出测量表面高低两种窄幅热色液晶的显色时间,根据窄幅液晶的显色温度及其时间通过求解一维半无限大平板非稳态导热方程得到对流换热系数.推导出此方法的误差传递公式,在此基础上分析如何得到准确的对流换热系数.最后给出热流法测对流换热系数的实验流程,并通过自由盘实验验证热流法测对流换热系数的可行性.     

基于CFD的水收集系数及防冰表面温度预测

利用计算流体力学(CFD,Computational Fluid Dynamics)软件FLUENT,对二维机翼在结冰气象条件下飞行时的水收集系数,以及防冰系统工作并达到稳定后的防冰表面温度进行了计算.水收集系数计算采用FLUENT中的欧拉两相流模型以及用户自定义函数.特定飞行条件及热气防冰系统结构下的防冰表面平衡温度分布预测采用防冰腔内外热流耦合的方法.所得结果和理论分析一致,表明了利用FLUENT可以有效的进行防冰研究,包括水收集系数及防冰表面温度的预测,为今后进一步研究设计热气或电热防冰系统打下基础.      

振动表面自然对流强化换热特性

基于振动影响流动换热边界层发展的思路,通过实验方法对振动条件下的自然对流换热特性进行研究.实验采用电铃谐振器作为加热膜的激振源,并利用红外测温技术对表面温度场进行了测量.结果表明:振动对自然对流的强化可提高90.7%.在等热流密度条件下,振动能量越大,换热越好;在等振动能量条件下,热流密度越小,换热越强.研究结果为强化表面自然对流换热提供了一种新思路.      

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

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

飞机发动机进气道前缘热气防冰器性能分析

飞机发动机防冰对飞行安全具有十分重要的意义,对于热气防冰系统,前缘防冰器的结构型式对防冰的效果影响很大.本文针对常见的周向及双蒙皮波纹板型弦向飞机发动机进气道热气防冰器进行了热分析,在此基础上对其传热性能及防冰表面温度场进行了比较,结果表明双蒙皮弦向防冰器的防冰效率较高;在通常情况下,周向防冰器也可达到防冰的目的,并且结构简单,因而更为实用.      

相变储热/辐射器式热沉的传热特性研究

对一种用于舱外活动的相变储热/辐射器式热沉进行了实验研究和数值模拟.设计了相变储热器的实验件和实验方案,进行了储热器的冷却实验.建立了具有导热与对流换热耦合边界的相变储热器数学模型,模型中考虑了自然对流及空穴的影响.编制了相变储热器的性能分析计算程序.通过实验值与计算值的比较分析,验证了程序的可靠性.该模型及计算程序为相变储热器的设计和优化提供了有利的分析工具.     

两种层板性能对比

针对常规的141型层板,提出了一种增加50％数量扰流柱的设计,即161型层板.采用实验和数值模拟相结合的办法研究了两种类型层板在流阻和换热方面的差别.换热实验部分采用了辐射式加热设备,测量了层板上下表面平均壁温,以及冷气进、出口截面间的焓差,得到了相应进气雷诺数下的层板体积换热努塞尔数.发现161模型流阻降低20％,换热增强5％.对一些实验工况,采用三维计算流体力学程序进行了流-固耦合传热数值模拟,所得结果在趋势上与实验一致,在数值上,流阻和换热与实验结果分别相差5％和30％.利用数值模拟结果分析并比较了两种模型在流场和表面对流换热系数分布方面的细节差别.      

用GAO-YONG湍流模型计算激波/边界层干扰

采用GAO-YONG可压缩湍流方程组,模拟了平板激波/湍流边界层干扰现象.运用SIMPLE算法求解方程组,并分别采用三阶精度的QUICK格式和中心格式离散对流项和扩散项.计算结果较好预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征: 分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波.对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布和速度型与实验值比较吻合很好.结果表明GAO-YONG可压缩湍流方程组能够高精度模拟平板激波/湍流边界层干扰流动.      

基于地面实验的蒙皮换热器高空换热性能分析

面向未来我国大功率机载电子设备冷却需求,需要发展新型机载冷却方式,减少对传统环控制冷系统的依赖,为此研究了一种气/气蒙皮换热器用于电子设备热管理.该蒙皮换热器是一种无能耗被动换热方式,在不增加飞机代偿损失的情况下实现大功率电子设备热管理.为了考核所设计的蒙皮换热器机载飞机环境换热性能,设计加工了实验样件,在地面开展了换热性能实验研究工作,获得热边流动传热准则和对流换热系数.在此基础上推演出高空巡航状态下的蒙皮换热器换热性能,获得传热最不利安装位置和改进位置时的传热性能.研究工作预期为未来我国基于蒙皮换热器的飞机综合热管理提供设计及机载应用技术支持.     

Boundary-layer analysis for the convection/diffusion transition in dendritic growth

Earlier ground-based studies have shown that a sharp transition occurs in the dominant mode of heat transport during dendritic solidifcation. Specifically, convective heat transport dominates below the transition supercooling and diffusive heat transport dominates above. Thus the influence of gravitational acceleration on dendritic crystallization depends critically on the overall thermodynamic driving force. Since many practical solidification processes occur at relatively small supercoolings, understanding the nature of the transition in transport mode is essential for both terrestrial and space processing. A boundary-layer analysis has been carried out wherein the diffusion boundary-layer thickness (Stefan length) is compared with the thermal boundary layer thickness determined from the fluid mechanics of natural convection. Our analysis shows that the Stefan length δ_S decreases rapidly with supercooling (ΔT^?2.5) whereas the thermal boundary layer thickness δ_T decreases slowly (ΔT^?.25). The transition supercooling at which the transport mode switches from convective to diffusive occurs where δ_T = δ_S. The influence of gravitational acceleration, and the thermo-physical properties of the molten phase will be discussed, along with avenues for conducting future research.     

火星表面热环境对航天器热控影响分析

火星大气对太阳辐射产生吸收和散射作用,同时还将与火星表面航天器发生对流换热.热设计时难以直接评估对流、辐射和导热三种换热对航天器的影响,从而确定主要的控温途径.在调研火星表面辐射、大气等热环境的基础上,从线性化传热系数和对流辐射比的角度对比分析了辐射、对流和导热对航天器的影响.器表辐射传热系数随光学属性和温度的变化范围...     

矩形通道中亚尺度肋片的流动换热数值分析

对装有不同结构亚尺度肋片矩形通道的流动和换热进行了数值模拟,获得了通道流场、温度场分布以及平均努塞尔数的基本特征,对各种肋片通道的换热特性进行了对比分析.计算结果表明:亚尺度肋片距主肋片根部越近越有利于散热;在远离热面区域,亚尺度肋片的比表面积越大传热效果越好.当扩展表面积相同时,亚尺度肋片的长宽比越大换热效果越好;计算结果及分析揭示了控制肋片最优几何形状的统一原则——广义温度梯度均匀化原则,并依照此原则确定了可以强化换热的较优的肋片结构.计算分析表明广义温度梯度均匀化原则可以控制对流换热过程的传热强化.      

矩形微槽内FC-72的单相流动和换热实验研究

微槽换热器是解决航空电子设备元器件热控制问题的一条有效途径.以新型电子设备冷却液FC-72为工质,利用4种不同结构尺寸的矩形微槽,进行了微槽内单相流动和强迫对流换热性能实验研究;分析了工质流速、过冷度以及微槽结构等对换热性能的影响,给出了适用于层流和紊流的单相流动阻力特性实验关联式和单相强迫对流换热实验关联式,并对实验结果进行了不确定度分析.结果表明,在矩形微槽换热器中,FC-72的单相强迫对流换热能力能够满足中高强度的电子设备散热要求.      

直升机旋翼防/除冰电加热控制律仿真

电热防除/冰系统的控制涉及到电加热与外流场的传热耦合,计算较为复杂,可利用的数据资料较为稀少。为探索电热防/除冰系统工作时与外流场的耦合传热规律,建立了二维电热除冰的数学模型,该模型在Messinger模型和改进的焓法模型基础上耦合了外表面与环境的复杂换热以及融冰和重新结冰过程的相变换热;采用控制容积法对控制微分方程进行离散后,使用TDMA（Tri-Diagonal Matrix Algorithm）和ADI（Alternating Direction Implicit）方法对离散得到的线性方程组进行求解,进而得到了除冰表面温度分布,同时揭示了电热防/除冰系统的耦合传热规律;分析了不同结冰条件下,加热时间控制律和加热热流密度对除冰表面温度的影响。计算发现合理设计加热热流密度大小及分布和加热时间控制律,可实现电热除冰系统能源的高效利用,进而确保飞行安全。     

贮箱内低温推进剂汽化过程的CFD数值仿真

为研究贮箱内低温推进剂相变对推进剂温度和贮箱压力的影响,对贮箱内的传热传质过程进行了仿真.仿真涉及的物理过程包括贮箱与外界环境的换热、推进剂的自然对流、推进剂与贮箱内壁面的换热以及低温推进剂的相变过程等.根据热力学平衡原理建立了低温推进剂相变模型,使用CFD(Computational Fluid Dynamic)方法对处于地面常压停放状态的液氢贮箱进行了450 s的仿真.研究表明随着贮箱壁面传热过程的稳定,推进剂的温度分布、流动状态以及相变情况会趋于稳定;通过仿真获得了推进剂单位时间的汽化量;影响相变的主要因素是贮箱壁面漏热以及推进剂自身的对流运动.     

粗糙度对涡轮叶片流动转捩及传热特性的影响

为研究表面粗糙度对涡轮叶片流动转捩以及传热特性的影响,在自行开发的CFD程序平台上提出了对γ-Re_θ转捩模型的粗糙度修正方法,并参考平板绕流和涡轮叶栅的实验数据对该方法进行验证。考虑粗糙度效应的影响,对Mark Ⅱ涡轮导叶5411工况进行数值模拟,得到如下结论:表面粗糙度对层流边界层换热系数影响不大,而对湍流边界层则有较大影响,进而显著改变壁面温度分布;与光滑壁面相比,5μm的等效沙粒粗糙度使吸力面湍流区域壁面温度升高约5.7K,100μm粗糙度使壁面温度升高28.4 K,增幅达5%左右;当壁面粗糙度较低时,激波干涉对吸力面边界层的转捩起主导作用,而当粗糙度大于某临界值时,其作用会使转捩位置突然变化,本算例中该临界值近似为150μm。      

含内热源的多孔方腔自然对流非正交MRT-LBM模拟

为了增强多孔方腔内流体流动与传热效果,采用非正交多松弛格子Boltzmann方法（MRT-LBM）对含有内热源的多孔方腔自然对流传热现象进行了数值模拟。研究了不同冷源布置方案（Scheme A~Scheme F）、内热源结构形式（Case 1、Case 2、Case 3）、内热源位置（a、b）、Darcy数、Rayleigh数等对多孔方腔内流体流动与传热的影响。计算结果表明:冷源布置方案对腔内流体流动与传热具有重要影响,当冷源左右对称布置时,腔内温度场及流场亦对称分布;在高Rayleigh数下采用Scheme A的双上部冷源布置方案能明显提高腔内的传热强度;内热源的形状对腔内对流传热影响很大,高Rayleigh数下,Case 3的布置方式更好。内热源的位置a和b对腔内的传热影响明显,提出了热壁面平均Nusselt数与位置a的拟合关系式,存在最佳的位置a（a=0.25）,使得腔内的对流传热最强;热壁面平均Nusselt数亦随b值变化表现出特定的变化规律。随着b值的增加,热壁面平均Nusselt数呈现先增后减再增的趋势。