数值传热大作业.docx

数值传热大作业.docx

《数值传热大作业.docx》由会员分享，可在线阅读，更多相关《数值传热大作业.docx（14页珍藏版）》请在冰豆网上搜索。

数值传热大作业

放置竖直孤立平板的二维围场内的空气流动与换热的数值分析

（西安交通大学能源与动力工程学院，710049，西安）

摘要：

针对内部放置孤立平板的二维围场内的空气流动与换热问题，在稳态、常物性和壁面温度以及孤立平板温度恒定的条件下，采用SIMPLER算法，对围场内部的空气进行了流动与换热的数值模拟计算。

在瑞利数Ra=10000时，计算得到了二维围场内的流线、等温线以及热线。

关键词：

SIMPLER算法、孤立平板、流线、等温线、热线

Abstract:

Inordertoinvestigatethefluxionandheattransferofairina2Dsquareenclosurewithanisolatedplate.SIMPLERalgorithmwasadoptedbasedontheReylondsconservationequationsofthesteady-stateconstantpropertylaminarflowandaconstanttemperatureoftheisolatedplateandtheinnerwallsoftheenclosurecondition.SlovefluidvelocityandtemperaturefieldsintheenclosureforRa=10000,anddrawthediagramsofstreamlines,isothermsandheatlines.

Keywords:

SIMPLERalgorithm;isolatedplate;streamlines;isotherms;heatlines

主要符号表

瑞利数

普朗特数

空气运动粘度

m2/s

重力加速度

kg.m/s2

空气导热系数

W/（m℃）

空气体膨胀系数

1/℃

空气比热容

J/kg.℃

空气密度

kg/（m3s）

金属板温度

℃

围场壁面温度

℃

温差

℃

一、引言

封闭空腔内孤立物体自然对流换热是一个重要的研究课题，从某种角度讲，大空间自然对流是封闭腔内孤立物体自然对流的一个特例。

文献[1,2]应用原始变量法求解了圆形空腔内存在同心的孤立圆弧状或八边形物体时的自然对流换热问题，文献[3]研究了封闭腔内竖直孤立物体的自然对流。

作为这一类问题的继续，本文着重研究了二维围场内竖直放置的孤立平板的自然对流，旨在计算分析得出在瑞利数等于10000时，二维围场内的流线、等温线以及热线。

二、物理问题与数值描写

本文研究的二维围场与孤立平板的布置如图1所示，正方形围场的壁面恒定温度为

℃，边长

，内部竖直放置的孤立平板温度恒定为

℃。

围场内其余地方填充着空气。

图1、二维围场平面示意图

数值计算是在如下基本假设的前提下进行的：

1、在整个计算过程中假设空气的物性是恒定不变的；

2、空气的粘性耗散忽略不计；

3、围场外表面假设是绝热的，与外界没有热交换；

4、空气流动以及换热处于稳态、层流状态下。

所研究的问题的控制方程如下：

连续性方程：

动量方程：

X方向：

Y方向：

能量方程：

其中边界条件：

围场壁面：

孤立平板：

温度恒定

四周边界速度

（1.8）

三、数值求解方法

此次数值模拟计算采用SIMPLER算法，并用区域离散方法B以及交错网格。

为保证迭代计算过程中孤立平板处的速度恒为零（或与主流区域相比要小若干数量级的小数），采用了以下方法[4]：

1、在每一层次的迭代计算前令孤立平板的速度为零，以保证孤立平板中的节点对流体区中的节点速度起滞止的影响；

2、在求解速度的代数方程前令孤立平板各速度离散方程主对角元的系数为一很大值（此次计算过程中，该值取为1030），以保证预估值

；

3、计算压力修正值时，使孤立平板各速度修正值计算公式的系数（即

等）取一个近似于零的值，如10-38，以使孤立平板中各速度修正值也为零。

四、计算结果与讨论

在瑞利数

条件下，计算得出速度场和温度场，并利用tecplot对计算数据进行处理得出流线、等温线、热线图如下：

图2、流线

从流线图我们可以看出，在孤立平板的左右两侧的流场中出现了许多封闭的流线，说明空气在上述区域形成了自然循环流动。

图3、等温线

对温度场中的等温线进行观察，我们可以很容易的得出在孤立平板的四个角处的等温线比较密集，说明四个顶角处的温度梯度比较大，而且围场壁面处的等温线也比较密集，温度梯度也比较大。

图4、热线

对热线进行分析，我们可以粗略的看出热流密度的方向是从二维围场墙壁指向孤立平板，这和围场墙壁温度（280℃）大于孤立平板的温度（180℃）的情况是相符合的，也验证了热量是从高温物体向低温物体单向传递的事实。

在孤立平板的右侧，由于空间足够大，热线形成了一个自我闭合的同心曲线簇。

五、结论

1、靠近围板壁面和孤立平板的区域，由于温差比较大，对流效果比较明显，导致速度梯度比较大，流线密集。

2、围板壁面和孤立平板四个顶角等处的温度梯度比较大，等温线比较密集；

3、热流密度的方向是由高温物体指向低温物体，即热量传递方向是由高温物体向低温物体单向传递的。

。

参考文献

[1]杨沫、陶文铨、陈钟欣，“缝隙度对封闭圆内开缝圆环自然对流换热的影响，”工程物理学报，Vol.12,No.2,1991,pp.193-196.

[2]杨沫、陶文铨、陈钟欣，“开缝八边形离相封闭母线自然对流换热的数值计算，”中国电机工程学报，Vol.11,No.4,1991,pp.58-64.

[3]王秋旺、王育清、陶文铨、杨沫，“几何位置对封闭方腔内水平孤立平板自然对流换热的影响，”工程物理学报，Vol.15,No.2,1994,pp.195-199.

[4]陶文铨编著，数值传热学（第二版），2001，西安交通大学出版社，244-244.

附录

源程序

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

MODULEUSER_L

C******************************************************************

INTEGER*4I,J,M00,M0,N0,N00

REAL*8TIN,TW,VIN,VOUT,AMU,COND,BETA,TREF,FLOWIN,FL,GRA,TM

C******************************************************************

ENDMODULE

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

SUBROUTINEUSER

C*******************************************************************

USESTART_L

USEUSER_L

IMPLICITNONE

C******************************************************************

C-----------------------PROBLEMSIX-----------------------

CLaminarfluidflowandheattransferinatwo-dimensional

C---------------------suddenexpansion--------------------

C******************************************************************

C

ENTRYGRID

TITLE

（1）='VEL_U'

TITLE

（2）='VEL_V'

TITLE（4）='TEMP.'

TITLE（11）='PRESSURE'

TITLE（8）='EX'

TITLE（9）='EY'

RELAX

（1）=0.8

RELAX

（2）=0.8

LSOLVE

（1）=.TRUE.

LSOLVE（4）=.TRUE.

LPRINT

（1）=.TRUE.

LPRINT

（2）=.TRUE.

LPRINT（4）=.TRUE.

LPRINT（11）=.TRUE.

LPRINT（8）=.TRUE.

LPRINT（9）=.TRUE.

LAST=200

XL=0.017

YL=0.017

L1=63

M1=63

N0=20

N00=28

M0=18

M00=44

CALLUGRID

RETURN

C

ENTRYSTART

TREF=20

DO100J=1,M1

DO101I=1,L1

U（I,J）=0.0

V（I,J）=0.0

T（I,J）=0.0

T（1,J）=280.0

T（L1,J）=280.0

T（I,1）=280.0

T（I,M1）=280.0

U（1,J）=0.0

U（L1,J）=0.0

V（I,1）=0.0

V（I,M1）=0.0

101ENDDO

100ENDDO

DO120J=M0,M00

DO121I=N0,N00

T（I,J）=180.0

U（I,J）=0.0

V（I,J）=0.0

121ENDDO

120ENDDO

AMU=3.472E-5

CPCON=1034.0

COND=3.931E-2

GRA=9.8

BETA=3.676E-3

RHO=0.746

RETURN

C

ENTRYDENSE

RETURN

C

ENTRYBOUND

DO34I=1,L1

DO35J=1,M1

IF（I==1）THEN

DTX（I,J）=（T（I+1,J）-T（I,J））*L1/XL

ELSEIF（I==L1）THEN

DTX（I,J）=（T（I,J）-T（I-1,J））*L1/XL

ELSEIF（N0.LE.I.AND.I.LE.N00）THEN

DTX（I,J）=0.0

ELSE

DTX（I,J）=（T（I+1,j）-T（I-1,J））*L1/XL/2.0

ENDIF

35ENDDO

34ENDDO

DO44J=1,M1

DO45I=1,L1

IF（J==1）THEN

DTY（I,J）=（T（I,J+1）-T（I,J））*M1/YL

ELSEIF（J==M1）THEN

DTY（I,J）=（T（I,J）-T（I,J-1））*M1/YL

ELSEIF（M0.LE.J.AND.J.LE.M00）THEN

DTY（I,J）=0.0

ELSE

DTY（I,J）=（T（I,j+1）-T（I,J-1））*M1/YL/2.0

ENDIF

45ENDDO

44ENDDO

DO24I=1,L1

DO25J=1,M1

EX（I,J）=0.746*U（I,J）*CPCON*T（I,J）-COND*DTX（I,J）

25ENDDO

24ENDDO

DO26J=1,M1

DO27I=1,L1

EY（I,J）=0.746*V（I,J）*CPCON*T（I,J）-COND*DTY（I,J）

27ENDDO

26ENDDO

RETURN

ENTRYOUTPUT

IF（ITER==0）THEN

PRINT401

WRITE（8,401）

401FORMAT（'ITER',7X,'V（4,7）',

19X,'T（4,7）'）

ELSE

PRINT403,ITER,V（4,7）,T（4,7）

WRITE（8,403）ITER,V（4,7）,T（4,7）

403FORMAT（I6,1P4E15.3）

ENDIF

IF（ITER==LAST）THEN

CALLPRINT

ENDIF

RETURN

ENTRYGAMSOR

DO500J=1,M1

DO501I=1,L1

GAM（I,J）=AMU

IF（NF==4）THEN

GAM（I,J）=COND/CPCON

IF（M0.LE.J.AND.J.LE.M00.AND.N0.LE.I.AND.I.LE.N00）THEN

CON（I,J）=180*1E30

AP（I,J）=-1E30

DU（I,J）=1E-25

DV（I,J）=1E-25

ENDIF

ENDIF

501ENDDO

500ENDDO

DO510J=2,M2

DO511I=2,L2

IF（NF==2）THEN

TM=T（I,J）-TREF

CON（I,J）=GRA*BETA*TM*0.746

ENDIF

IF（NF==2.OR.NF==1）THEN

GAM（I,J）=AMU

ENDIF

511ENDDO

510ENDDO

DO552J=M0,M00

DO551I=N0,N00

IF（NF==2.OR.NF==1）THEN

GAM（I,J）=1E38

ENDIF

551ENDDO

552ENDDO

RETURN

END