高温下加热/冷却管内参与性介质的辐射与对流耦合换热研究

本文通过数值模拟方法研究了高温下加热/冷却管内层流充分发展流动的参与性介质热起始段的辐射与对流耦合换热。采用离散坐标法求解圆柱坐标系下的非线性积分微分辐射传递方程,采用有限差分法求解辐射与对流的耦合换热能量方程。通过计算,分析了不同壁温及入口温度下,介质导热系数、壁面发射率、介质散射反照率及散射相函数等因素对局部努谢尔数、壁面热流密度及有效加热/冷却长度的影响,并同时采用努谢尔数和壁面热流密度评价了耦合换热中辐射对对流换热的影响。      

主动冷却通道内正庚烷热裂解与换热过程的数值模拟 *

为了深入研究碳氢燃料在主动冷却通道中的传热特性,采用详细基元反应机理对微小通道内正庚烷热解反应耦合流动换热过程进行了数值模拟。根据燃料的热裂解反应和壁面传热特性可将通道分为三个换热区域。通过对不同的入口温度、运行压力、壁面温度和通道直径等典型工况进行模拟,探索了正庚烷的热裂解反应和流动换热过程耦合的基本规律。研究结果表明,入口温度对通道整体的换热性能影响不大;保持其它参数不变,通道直径从2mm减小到0.5mm,其整体换热量降低约10%,单位面积的换热量提升约4倍。     

考虑发动机冷却通道固壁内耦合导热影响的低温甲烷超临界压力传热研究

 在考虑发动机冷却通道固壁内耦合导热影响的情况下,开展了低温甲烷在矩形冷却通道中的超临界压力湍流换热数值模拟研究;仔细分析了热流密度及管道几何形状对低温甲烷超临界压力下的流动和传热的影响;得到了流体速度、壁面温度、壁面热流密度等参数的详细变化情况以及Nusselt数的变化规律。计算结果表明:在考虑流固耦合作用时,上壁面施加的热流有一部分会通过固体壁面内的热传导,经由侧壁面传入超临界压力流体,并且随着热流密度的增加,经侧壁面传导的热流所占的比例也会随之增大;减小冷却通道内截面的高宽比,可以提高超临界压力下的换热效果,但流动压降会大大增加,因此冷却通道高宽比的选择需综合考虑传热与压力损失的影响,可以引入热性能参数作为参考;修正的Jackson&;Hall对流换热关系式基本可以适用于本文中的各种工况。     

应用Transition k-kl-ω转捩模型对内冷叶片气热耦合的数值模拟

选取NASA-Mark Ⅱ跨声速叶片为算例,研究了Transition k-kl-ω转捩模型在内冷叶片气热耦合计算中的应用,探讨了整场耦合与冷却通道内采用对流换热系数准则耦合的差异。结果表明,该转捩模型相比其它全湍流模型能够更准确预测附面层内的层流和转捩状况;由于Transition k-kl-ω转捩模型转捩前期采用层流动能来描述扰动的发展,避免了使用含有来流湍流度的经验公式,引入了"分裂机制"来描述层流与湍流脉动间的相互作用,并且在旁路转捩和自然转捩源项模化中加入了Tollmien-Schlichting波的影响,对强激波后的温度计算相比常用的间歇因子转捩模型与实验值更吻合;换热系数准则耦合用于冷却通道传热计算,避免了冷却通道边界条件带来的误差,计算结果与实验吻合较好,更易于工程应用。     

超临界压力下正癸烷在水平矩形冷却通道内的流动传热数值模拟

为了研究浮升力效应对超燃冲压发动机燃烧室不同位置冷却通道传热特性的影响,对超临界压力下正癸烷在水平矩形冷却通道内的流动传热-固体导热耦合过程展开了详细的数值模拟研究。重点考察了燃烧室不同位置的冷却通道中浮升力对温度以及热流分配的影响及其机理。结果表明：由浮升力引起的二次流动使燃烧室不同位置的冷却通道温度分布和热流分配呈现显著的差异;浮升力提升了燃烧室不同位置冷却通道的换热效果,其中受热方向和重力作用方向相同的冷却通道换热性能提升得最多;修正的Jackson&Hall经验公式不能预测浮升力对冷却通道壁面对流换热的影响,需要寻找其他的经验公式或使用CFD手段进行计算分析来解决这一问题。     

多孔隔热壁温度场气动传热耦合计算方法

研究了燃烧室多孔隔热壁温度场气动传热藕合计算新方法。特点是综合应用工程热力学、空气动力学、传热学(包括导热、对流换热及辐射换热)的基础理论,结合数值传热学的理论和方法,考虑了压力(压差)、流速(流量)、温度(温差)及热量等诸多参数的相互耦合影响的关系,包括气膜孔内壁与冷气对流换热的影响。所发展程序在给出燃烧室进出口必要的气动参数后,即可直接计算出带环状缝隙及多个气膜冷却孔的燃烧室隔热壁温度场。      

有无内部对流换热对导叶冷却效果影响的数值研究

根据导向叶片的曲面结构,运用剪应力输运方程(SST)湍流模型,在不同冷气入口压力条件下,对涡轮叶栅通道内部的三维流场和圆形气膜孔叶片型面的冷却效果进行了数值模拟.计算针对两种情形,即叶片内部有对流换热的情形和无对流换热的情形,结果表明随着冷气入口压力的增加,叶片有无内部对流换热的冷却效果也在增加;有内部对流换热的冷却效果较纯气膜冷却的高,尤其是在叶片前缘.      

管内高温介质流动入口段的辐射与对流耦合换热数值模拟

将求解辐射传递方程的离散坐标法与求解对流换热的控制容积法相结合 ,数值模拟了高温下圆管流动入口段参与性介质的辐射与对流耦合换热。考察了介质光学厚度及管壁温度对温度分布、壁面热流密度及局部努谢尔数分布的影响。结果表明 ,高温下介质辐射使介质的温度分布和对流换热状况发生较大改变 ,对管内流动入口段的换热起着重要作用     

方形再生冷却通道内超临界正癸烷湍流传热数值研究

为了深入理解再生冷却过程中碳氢燃料的超临界湍流传热特性,以正癸烷为研究对象,对其在非对称受热(上壁面外侧加热)方形通道内的流动传热进行了数值研究。在数值方法可靠性得到充分验证的基础上,详细探讨了进口温度和压力等运行参数对通道上壁面和侧壁面内侧平均壁温和平均努塞尔数分布的影响。计算结果表明:在运行压力接近临界压力且主流温度处于拟临界温度附近的综合条件下,流体热物性剧烈变化导致的类膜态沸腾效应,引起了上壁面内侧显著的传热恶化现象。同时,不平衡压差诱发的通道截面速度场异常分布,进一步影响了上壁面内侧近壁区域流体的换热性能。另外,类膜态沸腾效应致使上壁面热流更多地传递到侧壁面,导致侧壁面平均努塞尔数大幅增大,相比于正常换热的最大增幅约为50%。     

结构参数对圆转矩形喷管换热的影响

为了解结构参数对圆转矩形内喷管再生冷却换热的影响，设计了多个圆转矩形喷管，考虑了三种结构参数：转方位置、出口高宽比和出口圆角大小的影响。采用有限体积法求解三维可压缩的N-S方程对其内部流动和换热进行了数值模拟。湍流模型采用标准的k-ε双方程模型，壁面附近的流动和传热采用壁面函数法处理，速度与压力的耦合采用SIMPLE算法求解。结果表明：在型面一阶导数连续的情况下，转方位置对圆转矩形内喷管的换热影响不大；出口高宽比对圆转矩形内喷管的换热影响较大，出口高宽比不能太小，否则影响内喷管流场和换热；出口圆角大小影响内喷管周向上的温度分布，圆角太小造成周向温度分布不均匀。     

变物性与热辐射对发汗冷却过程的影响规律研究

采用两层k-ε湍流模型对高温条件下有发汗冷却时的矩形槽道内湍流流动和换热进行了数值模拟。计算结果表明:随着冷却气体流量的增加,在发汗冷却壁面附近主流气体的温度梯度大大减小,壁面温度、局部对流换热系数都迅速下降;变物性条件下同时考虑热辐射时数值计算得到的平均St相对值St/S_t0随修正注入率F/St0的变化比已知的常物性条件下的值小;随着高温气体温度的上升,变物性和热辐射对与高温流体接触的表面温度有较大影响,热辐射对壁面温度的影响比变物性更大,而对流换热系数所受的影响不大;当主流为高温气体(如:2500K)、而壁面温度较低时(如:800K)热辐射的影响不大,发汗冷却效果仍然很明显。      

飞机整体瞬态热状况的数值仿真研究

 在分析飞机内外热环境的耦合传热机制基础上，建立了描述飞机整体瞬态热状况的物理数学模型。引入壁面热流函数，将飞机蒙皮外的气动对流与辐射换热作用转化为舱内热分析的浮动热边界条件，实现了飞机内外耦合热作用的解耦计算，避免了难以实现的直接耦合求解。采用区域分解技术，将飞机整体分解为特性不同的多个求解区域。采用蒙特卡罗法求解舱内设备之间的辐射换热，采用热网络法求解设备内部及各设备之间的辐射-导热-对流耦合换热。利用该方法对某型飞机飞行过程中的瞬态热状况进行了数值模拟，获得了飞机舱内的瞬态温度场，分析了相关因素的影响，为飞机设计及相关研究提供了重要依据。     

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

客机APU舱温度场的数值计算

为研究大型客机APU舱传热过程和进行结构热分析,采用数值方法计算了APU舱温度场。采用商业CFD软件Fluent,选用Realizable k-ε湍流模型和S2S热辐射模型建立了APU舱内流动和传热数值计算模型,得到了流场、温度场和热流信息。APU舱中主要传热过程是辐射传热;APU防火罩和排气管隔热罩起到了热辐射遮热罩的作用;防火罩内的APU冷却空气对于降低防火罩温度有重要作用;排气管内的引射冷却空气有效降低了壁面温度。     

液体推进剂火箭发动机推力室再生冷却通道三维流动与传热数值计算

应用湍流模型对液体推进剂火箭发动机再生冷却推力室通道的流动与传热进行了三维数值模拟,冷却工质为氢气,其密度、导热系数、动力粘度随着温度和压力而变化,冷却剂比热容及金属固体物性随着温度而变化。计算采用标准k-ε双方程湍流模型及气-固耦合算法。结果表明:推力室燃气侧壁面的温度和热流密度的最高点均发生在喉部附近,喉部横截面固体区域最大温度梯度靠近燃气,喉部附近氢气在垂直主流方向的截面上产生了二次流。气固耦合面最大热流密度及最大对流换热系数同样位于推力室喉部附近。      

液体火箭发动机不同室压下冷却方案适用范围

为了精准评估不同冷却方案对高压液氧烃火箭发动机推力室传热特性的影响,建立了一套再生通道-液膜屏蔽-隔热镀层-辐射换热的整机模型,采用Ievlev半经验模型计算燃气侧壁面的对流换热过程,引入Shruvik安全裕度评估准则,计算推力室径向的分区温度和热流密度。基于某型大推力液氧煤油火箭发动机,研究了不同冷却结构组合的换热能力上限,分析了不同推力室压力对冷却设计方案的影响。结果表明:推力室压力在12 MPa及以下时,可主要依靠再生冷却技术满足冷却需求;在16 MPa及以下时需要配合内冷却环带满足冷却需求;在18 MPa及以下时需进一步设置隔热镀层提高热防护能力;室压在20 MPa甚至更高时,必须采用其他强化换热措施。      

考虑辐射效应的多孔材料传热有限元算法

应用有限元法求解一维半透明多孔介质辐射对流传热问题。通过对辐射传递方程与边界条件中积分项的离散 ,推导出有限元求解列式。对于多孔介质传热控制方程和气体传热控制方程的求解 ,采用时间积分技术对两相温度统一进行求解。而对于边界方程处理 ,则通过传热方程与边界方程的循环求解 ,求出任意时刻温度场的瞬态解。给出数值算例 ,得出了放置在通风口、同时受到伴随辐射作用的半透明多孔介质传热的瞬态解 ,讨论了部分参数对瞬态温度场以及换热效果的影响