方腔散射辐射与自然对流耦合传热的数值模拟.docVIP

MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h \* MERGEFORMAT SEQ MTSec \r 1 \h \* MERGEFORMAT SEQ MTChap \r 1 \h \* MERGEFORMAT 中国工程热物理学会 传热传质学 学术会议论文 编号：143490 方腔散射辐射与自然对流耦合传热的数值模拟 付丽荣，明平剑，张文平 （哈尔滨工程大学 动力与能源工程学院，哈尔滨，中国，150001） (Tel: 0451 Email:pingjianming@hrbeu.edu.cn) 摘要：为了研究辐射与自然对流非线性耦合问题，本文在自主软件GTEA基础上，针对封闭方腔内参与性介质的辐射与自然对流耦合换热进行模拟，计算结果与文献经典算例及现有商业软件Fluent结果吻合较好。并分析了普朗克数和介质散射反照率的变化对方腔内温度和速度的影响，结果表明：随着普朗克数的减小，增强了辐射换热的影响；考虑各向同性散射介质时，随着散射反照率的减小，中心处的温度有所增大，靠近绝热壁面处，速度U的变化更剧烈。 关键词：辐射传热；混合网格；有限体积方法；各向同性散射； 0 前言 封闭腔的自然对流与辐射换热的耦合作用在工程领域有着非常广泛的应用，比如太阳能集热器、空气制冷设备、电子元件的散热、绕核反应器中心的空气腔，建筑物的绝热、动力电站封闭母线及旋转电机的散热等，因而日益得到重视[1]。由于自然对流内温度与流速的内在的耦合关系，使得辐射换热对自然对流的影响比对强制对流更强烈[2]。已经有多人对封闭腔的自然对流换热问题进行了实验研究以及运用数值方法来对此进行数值计算，或者用它来验证算法与程序[3]。Y. Tian 等通过实验方法研究了竖直方腔内空气的低湍流自然对流，为CFD方法提供了一个标准算例[4]。M.A.Gad等用软件Fluent 6.3数值模拟了二维方腔内的表面辐射与对流，研究发现DOM法比S2S方法更好[5]。A. Mezrhab 等对不同加热程度的空腔内辐射和自然对流进行了数值模拟，研究发现辐射改变了方腔内均匀的温度分布，而且提高了平均Nusselt数[6]。G.Colomer等研究了三维参与性透明介质内的自然对流与辐射换热，用DOM法求解辐射方程，研究发现光学性质对参与性介质辐射换热起着重要的作用[7]。明平剑等用有限体积方法研究了含吸收性非散射性介质的辐射换热与自然对流耦合的问题，分析了数和数对腔内温度的影响[8]。 本文在文献[8]的基础上，进一步研究包含有对辐射具有吸收、发射和各向同性散射的参与性介质的辐射换热与自然对流耦合作用。用Fortran语言编写了一套自然对流与散射辐射耦合换热程序，采用非结构有限体积方法求解辐射输运方程，求解得到辐射强度分布函数后，计算单元中心的热源和辐射热流密度，然后求解质量、动量和能量方程。 1数学模型 流体控制方程 基金项目：国家自然科学基金资助项目 流体控制方程可写作如下微分形式[8]： MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) 采用Boussinesq假设考虑由于温度分布引起的密度差所致浮升力： MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 4) 式中，是流体密度，是流体速度矢量，是压力，为动力粘性系数，为比定压热容，是热传导系数，为温度，为辐射源项，是参考密度，是热