辐射加热过程中非定常耦合传热数值方法.pdfVIP

辐射加热过程中的非定常耦合传热数值方法 明平剑，龚京风，宣领宽，付丽荣，张文平 (哈尔滨工程大学动力与能源工程学院，哈尔滨，中国，150001 联系电话：0451 E-mail：pingjianming@hrbeu．edu．cn 摘要：对流一结构耦合传热广泛存在于热动力机械、飞行器气动加热、地面结构热试验等 问题中，随着温度升高，辐射传热所占总传热量比重不断增加，而且很多情况下需要考虑 加热过程的时间效应。。在此环境下，热辐射不仅改变温度分布，而且还影响流动分布以 及动量和能量传递¨1。A．Mezrhab等对不同加热程度的空腔内辐射和自然对流进行了数值 模拟，研究发现辐射改变了方腔内均匀的温度分布，而且提高了平均Nusselt数心1。到目 前为止，先后发展了很多预测自然对流问题的数值方法，流函数方法及原始变量法，有限 差分法和有限体积方法等。各方法具有独自的优势，有限体积方法因其守恒性及对任意网 格的通用性而广泛被应用工程计算流体力学软件中。由于有限体积法可以方便地用来求解 对流和热传导，数值传热学问题也主要采用有限体积法。辐射传热本质上不同于对流和导 and 热，在计算方法上也差别较大。为了统一流体力学与传热学计算方法，G．D．Raithby 进和发展。M．Y．Kim等研究了立体角剖分的问题并分析了封闭腔内格板对辐射的影响 。“。C．M．Kim进一步将非结构化网格有限体积法扩展到多边形网格哺1。张吴春等研究了 有限体积法求解辐射问题的计算精度问题瞄1。贺志宏等用有限体积方法研究了炉内辐射 传热问题m3。作者进一步发展了辐射输运方程求解的有限体积法，使其适用于任意混合 网格，从而与CFD求解器相匹配即1。 本文针对流I卉I耦合传热问题进行了研究，提出了一种考虑辐射效应的流一结构耦合传 热数值计算方法，发展了基_丁^非结构化网格有限体积方法的求解热辐射输运方程，然后计 算流体单元中心辐射热源和流川界面上的辐射热流量，最后求解流体运动方程和能量守恒 方程，计算得到速度、压力和温度分布。重点探讨了流吲耦合传热与辐射输运方程的耦合 计算方法对计算结果的影响。本文方法为辐射作用下的自然对流与结构耦合传热问题研究 提供了一种新的计算方法。 关键词：辐射传热；流同耦合传热；非稳态加热过程； 1控制方程 流吲耦合传热控制方程是一组非线性的质量、动量和能量守恒方程。控制方程可写 作如下微分形式： V．U=0 (1) (2) 竽冉(p删)=一跏+V．p[VU+(Vu)厂p 警丹(JDuh)_v．(胛力+薯 采用Boussinesq假设考虑由于温度分布引起的密度差所致浮升力． pb=Poflf厂一ro)g 对于漫灰介质，不考虑介质的散射，辐射输运方程可以写作如下形式 兰：刊+呱 Qs 对于固体边界，辐射强度边界条件满足下面条件： 名=s名+宰磊。∥|睇 (5) 固体区域能量守恒方程如下： 等珊(胛7)+‘ 求解得到辐射强度后，可以通过如下方程计算得到能量方程源项。 ‘=K(40-74_p，棚] 流同耦合界面辐射热流通过下式计算得到：