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

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

高温圆管内多孔介质耦合换热的辐射作用

通过数值模拟圆管内填充多孔介质的辐射对流耦合换热,研究高温固体骨架辐射效应对温度分布及换热的影响.基于局部非热平衡模型分别建立流体和固相能量方程,采用蒙特卡罗法求解固体骨架的辐射换热;对不考虑热辐射引起的温度场偏差和管壁发射率以及多孔结构参数的影响进行讨论.结果表明:多孔固体骨架的辐射效应对入口段温度场的影响明显,不考虑辐射将导致较大偏差,壁温为1500K时最大偏差为16%.管壁发射率对温度场影响较小,壁温为1500K时影响小于3%.孔隙率或孔径增大,壁面辐射热流密度比例增加,辐射效应体现明显.      

飞机发动机冷气道与隔热层的耦合传热分析

数值研究了某飞机发动机外侧冷气道与隔热层的耦合传热过程。采用低雷诺数k-ε模型与SIMPLEC算法计算通道内可压缩变物性气流的湍流对流换热,采用蒙特卡罗法求解通道壁面间的辐射换热。通道内气流湍流对流换热、壁面间辐射换热与隔热层内导热耦合求解。通过模拟计算,分析了通道与隔热层的耦合传热机制,考察了相关参数的影响。结果表明,在所考虑的通道结构与空气流条件下,冷气道外环壁面的温度高于气流温度,气流对内外环壁面均起冷却作用;在隔热层参数不变条件下,壁面间的辐射换热与气流的对流冷却是该传热过程的控制机制,增大冷气流量、降低壁面发射率均可显著降低隔热层的外壁面温度。     

矢量喷管内燃气辐射与壁面温度的耦合计算

选取内、外调节片和隔热屏建立几何模型,基于封闭腔净辐射模型和壁面热平衡模型建立了燃气辐射与喷管壁面温度的耦合算法.波段为1~5μm的气体辐射采用窄波带模型计算,其他波段不考虑气体辐射,建立辐射净热流密度-有效辐射亮度-壁面温度的关联式求解燃气与壁面的辐射换热,采用牛顿-拉斐尔森迭代法求解壁面热平衡方程计算其温度.对某轴对称矢量喷管(偏转20°),计算了喷管壁面的红外光谱辐射和辐射净热流密度,以及各部分结构的温度.作为验证,还计算了文献中某液体火箭发动机轴对称矢量喷管壁面的辐射净热流密度,与文献的结果进行对比一致性较好.研究表明:轴对称矢量喷管偏转段沿周向的辐射热流密度和温度差异很大,沿偏转方向部位壁面的温度和辐射热流密度都较低,偏转方向壁面的温度比相反方向大约低10%,辐射热流密度大约低50%.      

三维数值模拟再生冷却喷管的换热

为了解液体火箭发动机喷管再生冷却的换热特点，采用数值模拟的方法，对内喷管燃气、壁面和冷却液建立不同的三维控制方程，进行流动和传热的耦合计算。在计算中，假定喷管流动为冻结流动，考虑燃气向壁面的对流换热和辐射换热；采用二阶迎风格式离散控制方程，采用DO模型离散求解辐射换热方程，水蒸气的吸收系数根据Leckner公式计算。计算模型采用缩比热试车发动机，数值计算结果与实验结果吻合较好，较准确地模拟出了喷管的壁面热流密度，得到了喷管燃气和冷却液的流场和温度场，对高压再生冷却喷管的设计具有指导意义。     

不同热边界对超音速机翼气动热数值分析的影响

飞行器在超音速飞行时受到的气动加热效应给结构强度及热防护设计带来极大影响，且真实状态下的气动热环境需要考虑外流场与结构的耦合及内壁面边界条件的影响。采用S-A湍流模型求解Navier-Stokes方程，通过流场与固体壁面交界处的信息传递，实现外流场与结构场的耦合数值分析。针对三种不同翼型的超音速绕流气动加热进行耦合数值研究，对比翼型内壁面在不同热边界条件下的气动热效应，结果表明：不同翼型具有与气动力相似的的气动热效应；内壁面考虑对流换热的边界条件最接近真实；考虑机翼燃油箱满油时，三维机翼前缘驻点处热流密度最高可达4200w/m2。     

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

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

非正交界面下的插值耦合传热方法

针对复杂热防护结构,开发了一种求解非正交性流/固界面的耦合传热程序。流体和固体区域采用同一积分、守恒型的RANS(Reynolds-Averaged Navier-Stock)方程,通过有限体积法进行离散求解。为了保证非正交界面上的温度和热流密度连续,提出了一种结合网格周边信息计算界面热流密度的插值方法。利用该插值耦合方法模拟了双层复合平板和喷管的热传导。数值结果表明:该插值方法在空间上具有2阶精度;喷管壁面上的对流换热系数沿轴向先增大后减小,在喉部上游达到最大值,当喷管入口压力增加3.38倍时,对流换热系数的最大值相应增加了3.13倍;喉衬与壳体界面上温度的计算值和试验结果存在一定差异,这是由于本文数值计算未考虑接触热阻引起的。     

火箭发动机喷管辐射及壁温的计算研究

研究了火箭发动机复合喷管的耦合换热及瞬态壁温的计算方法,考虑了扩散控制的C/C喉衬的稳态烧蚀,建立了封闭腔-净辐射模型计算喷管内的辐射换热,改进了喷管辐射换热的计算方法.计算了两种复合喷管的壁温随时间的变化,并与有关文献的结果进行对比.研究表明:在火箭发动机喷管的入口段,高温燃气对壁面的辐射换热热流比对流换热热流强;采用封闭腔-净辐射法计算火箭发动机喷管的辐射精度较高.      

TBCC燃烧室/喷管一体化壁面温度计算

基于Navier Stokes(N-S)方程组对包括隔热屏、隔热屏内外流、大气外流在内的涡轮基组合动力(TBCC)发动机燃烧室/喷管进行了一体化的气/热耦合数值模拟,考虑了燃气组分输运、辐射换热等影响,研究了其在某典型飞行状态下TBCC冲压发动机燃烧室/喷管筒体及隔热屏内外壁壁面温度、辐射换热热流及对流换热热流分布.结果表明:燃烧室/喷管筒体与对称面上下交线的壁面温度在轴向距离为0.5~2.6m内变化较小,在轴向距离为2.6~3.1m内急剧增加,在轴向距离为3.1~3.5m内急剧下降.之后,上交线筒体壁面温度沿流向减小,下交线筒体壁面温度先升高后降低.筒体壁面温度最高点在喷管下调节板收缩段,为1577K.隔热屏内壁面辐射热流在370~500kW/m2变化,上下交线处的辐射热流较外壁面的辐射热流约高300kW/m2,辐射热流沿流向先减小后增加.隔热屏外壁面辐射热流在50~200kW/m2范围内分布.      

涡轮发动机高温隔热涂层内的传热研究

高温环境中的金属部件, 通常采用表面涂层技术以降低金属材料的温度和热力梯度。氧化锆、陶瓷等涂层材料对热辐射是半透明的, 因此, 高温下涡轮发动机的热分析需要考虑涂层内辐射换热的影响。本文考虑了各种界面辐射特性: 不透明界面和半透明界面, 镜反射和漫反射。采用光线踪迹 /节点分析法计算了吸收、各向同性散射性介质内的辐射传递系数。结合控制容积法数值模拟了涡轮叶片和燃烧室衬垫涂层内的辐射与导热瞬态耦合换热。      

金属蜂窝夹芯面板有效导热系数的数值计算

对具有漫灰体内壁的金属蜂窝夹芯面板,忽略蜂窝内腔的空气导热,综合考虑蜂窝结构的热传导和热辐射等热传递过程,利用有限元方法求解了周期性分布的蜂窝单胞稳态热传导控制方程,蜂窝内壁的边界条件是由净热量法得到的热辐射换热积分方程,由胞元的温度场分布数据及Fourier定律得到了蜂窝结构的有效导热系数。与现有文献相比,采用了较少近似的模型及较高精度的离散方法,计算结果表明计算模型及方法是可靠有效的。     

气动加热下光学窗口传热特性的数值模拟

针对气动加热条件下光学窗口辐射导热耦合传热开展研究,以了解光学窗口的热辐射发射规律及热防护性能。采用有限体积法、蒙特卡罗法、谱带模型建立了具有镜反射界面的光学窗口辐射传输及辐射导热耦合传热计算模型,对具有两层SiO2玻璃的光学窗口传热特性进行了模拟。结果表明处于高温区玻璃向外发射热辐射,低温区玻璃吸收热辐射;除了加热面附近外的玻璃区域,辐射热流密度远大于导热热流密度,在2.9～4.2μm谱带内的辐射热流密度大于0.4～2.9μm谱带内的;在0.4～2.9μm谱带内,可见光波长内的出射辐射热流密度最小,中红外波长范围内出射辐射热流密度最大。在光学窗口总厚度一定的情况下,降低玻璃层厚度和玻璃夹层气体热导率能够有效降低光学窗口非加热面稳态温度。     

导热-辐射耦合传热的多尺度分析和数值模型

复合材料高温传热特性的准确预测对飞行器热防护结构的设计有重要意义,相关模型也是国家数值风洞工程中多相多介质计算模型的重要组成部分。针对周期性结构复合材料的高温传热问题,利用多尺度渐进分析方法,对包含导热方程和辐射传输方程的导热-辐射耦合传热模型开展了研究。建立了表征单元尺度模型及宏观平均导热-辐射耦合传热模型,获得了材料宏观等效导热系数与表征单元模型间的关系,发现宏观等效辐射吸收和散射等系数可通过表征单元内的体积平均求取。结合有限容积方法与格子Boltzmann方法,建立了复合材料导热-辐射耦合传热多尺度数值模型。采用二维常物性材料传热过程的模拟验证了多尺度模型的有效性,结果表明所建立的多尺度模型能够对温度场给出准确高效的计算结果。该方法有助于为复合材料传热特性的预测提供数值手段。     

纳米隔热材料热辐射特性的理论计算

高温下纳米隔热材料内热辐射的影响将显著增强,其热辐射特性对热辐射传热有很大影响。为了认识高温纳米隔热材料的热辐射特性,采用Mie理论建立了掺杂纤维增韧剂和遮光剂的纳米隔热材料热辐射特性理论计算方法,编写了纳米隔热材料热辐射特性计算程序,对某纳米隔热材料的热辐射特性进行了理论研究,得到了光谱衰减系数、光谱散射反照率以及全光谱平均辐射特性参数及散射相函数。理论模拟结果表明:在3～9.5μm波长范围内,纳米隔热材料对热辐射具有强烈的衰减作用,对3～7μm的热辐射呈现强烈的散射特征,对7～9.5μm的热辐射,随波长增大散射特征逐渐减弱,对9.5μm的热辐射呈现较强的吸收特征。在300～1 300 K,该纳米隔热材料全光谱平均衰减系数>5×104 m-1,平均散射反照率>0.96,具有较强的前向散射特征,这些特征来源于遮光剂粒子,增韧剂影响很小。     

半透明高温涂层的内部辐射效应

为了准确计算涂层和金属基底内的温度和热流，需要考虑涂层内部热辐射的影响。采用光线踪迹／节点分析法对吸收、各向同性散射性非灰涂层内的辐射与导热瞬态耦合换热进行数值模拟。与多个文献的计算结果比较表明，计算方法正确，计算结果精度高。数值模拟结果表明，忽略涂层内的热辐射会使实际运动中的金属部件温度高于设计温度。