数值传热大作业.docx

放置竖直孤立平板的二维围场内的空气流动与换热的数值分析（西安交通大学能源与动力工程学院，710049，西安）摘要：针对内部放置孤立平板的二维围场内的空气流动与换热问题，在稳态、常物性和壁面温度以及孤立平板温度恒定的条件下，采用SIMPLER算法，对围场内部的空气进行了流动与换热的数值模拟计算。在瑞利数Ra=10000时，计算得到了二维围场内的流线、等温线以及热线。关键词：SIMPLER算法、孤立平板、流线、等温线、热线Abstract: Inorder to investigate the fluxion and heat transfer of air in a 2D square enclosure with an isolated plate. SIMPLER algorithm was adopted based on the Reylonds conservation equations of the steady-state constant property laminar flow and a constant temperature of the isolated plate and the inner walls of the enclosure condition. Slove fluid velocity and temperature fields inthe enclosure for Ra=10000,and draw the diagrams of stream lines ,isotherms and heat lines.Key words:SIMPLER algorithm; isolated plate; stream lines; isotherms; heat lines主要符号表瑞利数普朗特数空气运动粘度m2/s重力加速度kg.m/s2空气导热系数W/(m℃)空气体膨胀系数1/℃空气比热容J/kg. ℃空气密度kg/(m3s)金属板温度℃围场壁面温度℃温差℃一、引言封闭空腔内孤立物体自然对流换热是一个重要的研究课题，从某种角度讲，大空间自然对流是封闭腔内孤立物体自然对流的一个特例。文献[1,2]应用原始变量法求解了圆形空腔内存在同心的孤立圆弧状或八边形物体时的自然对流换热问题，文献[3]研究了封闭腔内竖直孤立物体的自然对流。作为这一类问题的继续，本文着重研究了二维围场内竖直放置的孤立平板的自然对流，旨在计算分析得出在瑞利数等于10000时，二维围场内的流线、等温线以及热线。二、物理问题与数值描写本文研究的二维围场与孤立平板的布置如图1所示，正方形围场的壁面恒定温度为℃，边长，内部竖直放置的孤立平板温度恒定为℃。围场内其余地方填充着空气。图1、二维围场平面示意图数值计算是在如下基本假设的前提下进行的：在整个计算过程中假设空气的物性是恒定不变的；空气的粘性耗散忽略不计；围场外表面假设是绝热的，与外界没有热交换；空气流动以及换热处于稳态、层流状态下。所研究的问题的控制方程如下：连续性方程：动量方程：X方向： Y方向：能量方程：其中边界条件：围场壁面：孤立平板：温度恒定四周边界速度(1.8)三、数值求解方法此次数值模拟计算采用SIMPLER算法，并用区域离散方法B以及交错网格。为保证迭代计算过程中孤立平板处的速度恒为零（或与主流区域相比要小若干数量级的小数），采用了以下方法[4]：在每一层次的迭代计算前令孤立平板的速度为零，以保证孤立平板中的节点对流体区中的节点速度起滞止的影响；在求解速度的代数方程前令孤立平板各速度离散方程主对角元的系数为一很大值（此次计算过程中，该值取为1030），以保证预估值；计算压力修正值时，使孤立平板各速度修正值计算公式的系数（即等）取一个近似于零的值，如10-38，以使孤立平板中各速度修正值也为零。四、计算结果与讨论在瑞利数条件下，计算得出速度场和温度场，并利用tecplot对计算数据进行处理得出流线、等温线、热线图如下：图2、流线从流线图我们可以看出，在孤立平板的左右两侧的流场中出现了许多封闭的流线，说明空气在上述区域形成了自然循环流动。图3、等温线对温度场中的等温线进行观察，我们可以很容易的得出在孤立平板的四个角处的等温线比较密集，说明四个顶角处的温度梯度比较大，而且围场壁面处的等温线也比较密集，温度梯度也比较大。图4、热线对热线进行分析，我们可以粗略的看出热流密度的方向是从二维围场墙壁指向孤立平板，这和围场墙壁温度（280℃）大于孤立平板的温度（180℃）的情况是相符合的，也验证了热量是从高温物体向低温物体单向传递的事实。在孤立平板的右侧，由于空间足够大，热线形成了一个自我闭合的同心曲线簇。五、结论1、靠近围板壁面和孤立平板的区域，由于温差