数值传热学(课件)-2

《数值传热学(课件)-2》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第12章求解二维椭圆型流动与换热问题的通用程序应用举例PresentedbyProfessorWQTao2003-12-24,西安目录12－1例题1直角坐标中二维稳态无内热源的导热12－2例题2空心圆柱内的稳态热传导12－3例题3正方形管道内的充分发展对流换热12－4例题4内壁上有直肋的环形通道内充分发展对流换热12－5例题5二维突扩通道中的流动与换热12－6例题6方形通道内的复杂充分流动12－7例题7旋转圆盘上的冲击流动12－8块修正技术12－9例题10带中心射流的通道内的紊流换热12－1直角坐

2、标中二维稳态无内热源的导热－熟悉USER的结构12-1－1物理问题及其数学描写已已知：稳态常物性、无内热源的导热问题，求解区域如图1所示，四个边界上的温度由下式决定：T=x+y+xy求：该区域中的温度分布。与通用方程相对比，相当于＝1，＝0.的稳态扩散问题。yx21图1******************************************************************--------------------EXAMPLEONE--------------------

3、-----------****************************************************************DIMENSIONT(22,22)EQUIVALENCE(F(1,1,4),T(1,1))*ENTRYGRIDLAST=10LSOLVE(4)=.TRUE.LPRINT(4)=.TRUE.TITLE(4)='.TEMP.'XL=1.YL=2.L1=7M1=7CALLUGRIDRETURN12-1－2程序阅读ENTRYSTARTDO100J=1,M1D

4、O100I=1,L1T(I,J)=0.IF(I.EQ.1.OR.I.EQ.L1)GOTO101IF(J.EQ.1.OR.J.EQ.M1)GOTO101GOTO100101T(I,J)=(X(I)+Y(J)+X(I)*Y(J))!不变的边界条件可在START中规定100CONTINUERETURN*ENTRYDENSERETURN！空块*ENTRYBOUNDRETURN！空块ENTRYOUTPUTIF(ITER.NE.0)GOTO400WRITE(8,401)401FORMAT(1X,'ITER',

5、13X,'T(4,4)',14X,'T(5,3)')400WRITE(8,403)ITER,T(4,4),T(5,3)！每次迭代输出两个温度，！观察收敛情况403FORMAT(1X,I5,2F20.6)IF(ITER.EQ.LAST)CALLPRINT！收敛后输出二维温度场RETURN*ENTRYGAMSORIF(ITER.GT.0)GOTO501！常物性一次赋值即可DO500J=1,M1DO500I=1,L1500GAM(I,J)=1.501CONTINUERETURNENDCOMPUTATIO

6、NINCARTESIANCOORDINATES****************************************ITERT(4,4)T(5,3)00.0000000.00000011.9999781.72036422.0000001.72000132.0000001.72000042.0000001.72000052.0000001.72000062.0000001.72000072.0000001.72000082.0000001.72000092.0000001.72000010

7、2.0000001.72000012-1－3计算结果分析**************************.TEMP.**************************I=1234567J72.00E+002.30E+002.90E+003.50E+004.10E+004.70E+005.00E+0061.80E+002.08E+002.64E+003.20E+003.76E+004.32E+004.60E+0051.40E+001.64E+002.12E+002.60E+003.08E+00

8、3.56E+003.80E+0041.00E+001.20E+001.60E+002.00E+002.40E+002.80E+003.00E+0036.00E-017.60E-011.08E+001.40E+001.72E+002.04E+002.20E+0022.00E-013.20E-015.60E-018.00E-011.04E+001.28E+001.40E+0010.00E+001.00E-013.00E-015.00E-017.00E-019.00E-011.00E+0