########学院计算流体力学与传热学学号：专业：学生姓名：任课教师：教授2013年12月目录第一章验证显式格式的稳定性 (4)1.1 概述 (4)1.2 数学推导 (4)1.3 问题描述 (4)1.4 数值模拟 (4)1.5 结果及分析 (5)第二章判断肋片可以按一维问题处理的主要依据 (6)2.1 概述 (6)2.2 问题描述及算法 (6)2.3 数值模拟 (7)2.4 结果及分析 (8)第三章三层墙导热 (9)3.1 概述 (9)3.2 问题描述 (9)3.3 TDMA算法 (9)3.4 结果 (10)第四章一维无源稳态对流扩散问题 (11)4.1 公式及初值 (11)4.2 情况一 (11)4.3 情况二 (12)4.4 情况三 (13)第五章用ADI算法计算长方肋内的温度分布 (14)5.1 问题描述 (14)5.2 初始参数 (14)5.3 情况一，一列列扫 (14)5.4 情况二，一行行扫 (14)5.5 情况三，采用ADI算法 (15)5.6 结果分析 (15)参考文献 (16)第一章 验证显式格式的稳定性1.1 概述将一维非稳态热传导方程用显式格式差分化为代数方程，在求解的迭代过程中必须满足一定的条件，才能使方程收敛且结果正确。

此处即验证β≤½。

1.2 数学推导方程： 22T t T x α∂∂=∂∂（1）显式离散格式： 此处时间向前差分，空间中心差分11122n n n n ni i i i i T T T T T t x α+-+--+=∆∆1112(2)n n n n ni i i i i t T T T T T xα+-+∆-=-+∆ 令β=2tx α∆∆则： 111(2)n n n n ni i i i i T T T T T β+-+-=-+ （2）误差也应该满足上式，故：()()1()()()2()()iiiiiIkx Ikx Ik x x Ikx Ik x x n n n n nT e T e T e T e T e ψψβψψψ----∆--+∆+⎡⎤-=-+⎣⎦()()()1()12()()()iiiiIkx Ikx Ik x x Ik x x n n n nT e T e T e T e ψβψβψψ----∆-+∆+⎡⎤=-++⎣⎦()()1()12()()iiiIkx Ikx Ikxn n Ik x Ik x n T e T e e e T e ψβψβψ---+-∆∆=-++()()1()121()n Ik x Ik x nT e e T ψββψ+-∆∆=-++≤ 因此 β≤½。

即当β≤½ 时方程（2）才会有收敛的解。

1.3 问题描述在验证过程中同时可模拟一个实际问题，即冬季里墙壁中的温度分布。

此时室内壁温设为Tl=30.0℃，室外壁温Tr=-25.0℃，墙壁以11号楼为例，L=1m ，热扩散系数ɑ=alfa=1.33e-6m 2/s然后分别取β=0.4，n=10和β=0.6，n=10两种情况，看最后的结果是否收敛和正确。

1.4 数值模拟程序流程如下：1.5 结果及分析㈠当beta=0.4、n=10时，运行程序可收敛，得结果如下：30.00 24.47 18.94 13.42 7.91 2.41 -3.09 -8.58 -14.06 -19.53 -25.00 30.00 24.47 18.95 13.43 7.91 2.41 -3.09 -8.57 -14.05 -19.53 -25.00 30.00 24.47 18.95 13.43 7.92 2.41 -3.08 -8.57 -14.05 -19.53 -25.00 30.00 24.47 18.95 13.43 7.92 2.42 -3.08 -8.57 -14.05 -19.53 -25.00 30.00 24.48 18.95 13.44 7.92 2.42 -3.08 -8.56 -14.05 -19.52 -25.00可看出最后已经稳定。

用图表显示可得：可看出最后温度在墙壁内基本成线性分布。

输入ɑ、β、L 、n 值 设定初值：T(0)=Tl T(i)=Tr(i=1，n)迭代计算：TN(i)=T(i)+β×(T(i+1)+T(i-1)-2×T(i))求解输出 判断是否收结束是 否T(i)=TN(i)开始㈡当beta=0.6、n=10时，运行程序不收敛，得到如下结果：30.00 -184.31 395.52 -498.83 562.01 -562.79 483.59 -411.54 253.44 -160.76 -25.00 30.00 292.17 -488.99 674.28 -749.37 739.92 -681.32 524.52 -394.07 169.21 -25.00 30.00 -333.83 677.67 -877.87 998.39 ******* 894.93 -750.14 495.06 -285.28 -25.00 30.00 491.37 -862.55 1181.21 ******* 1337.27 ******* 984.02 -720.26 339.09 -25.00 30.00 -597.80 1176.06 ******* 1777.14 ******* 1639.36 ******* 937.92 -514.98 -25.00 30.00 843.19 ******* 2082.30 ******* 2410.96 ******* 1820.11 ******* 650.75 -25.00 30.00 ******* 2060.31 ******* 3169.92 ******* 2985.10 ******* 1746.07 -935.83 -25.00 30.00 1467.32 ******* 3688.80 ******* 4341.70 ******* 3337.52 ******* 1219.81 -25.00 综上所述，只有当beta≦½时，结果才会稳定收敛。

第二章判断肋片可以按一维问题处理的主要依据2.1 概述在实际问题中，经常将肋片按一维处理，即假设表面上的换热热阻1/h远远大于肋片中的导热热阻δ/λ（即Bi很小），因而在任一截面上肋片温度可认为是均匀的。

这样假设的确便于数学处理，但是否合理，下面就进行验证。

2.2 问题描述及算法粗而短的肋片的分析节点划分计算在下表所示的两种条件下肋片的效率，并与一维分析解的结果比较。

程序设计中各个点的温度为：1 T(i,j)=0.25*(T(i-1,1)+T(i+1,1)+2*T(i,2))2 T(i,j)=0.25*(T(i+1,j)+T(i-1,j)+T(i,j+1)+T(i,j-1)) 3a T(i,j)=(T(i-1,N)+T(i+1,N)+2*T(i,N-1))/(4+2*Bi) 3b T(i,j)=(T(M,j-1)+T(M,j+1)+2*T(M-1,j))/(4+2*Bi) 4a T(i,j)=(T(M-1,N)+T(M,N-1))/(2+2*Bi)4b T(i,j)=(T(M-1,1)+T(M,2))/(2+Bi)2.3 数值模拟程序流程图为：工况 T0/℃ Tf/℃ h/W/(m 2.℃) λ/W/(m.K) δ/m H/m 1 100 20 50 100 0.02 0.04 21002040080.020.08以计算得到的T(i,j)计算临近节点的温度否是结束 判断是否稳定：开始 输入计算参数 h ，λ，δ,L ，M ，N 求Bi设定温度初值: J=1,N T(1,j)=T0-Tf 其余的节点为T(i,j)=0.0 令c(i,j)=T(i,j) 根据初值计算节点温度(共6类) 1 T(i,j)=0.25*(T(i-1,1)+T(i+1,1)+2*T(i,2)) 2 T(i,j)=0.25*(T(i+1,j)+T(i-1,j)+T(i,j+1)+T(i,j-1)) 3a T(i,j)=(T(i-1,N)+T(i+1,N)+2*T(i,N-1))/(4+2*Bi) 3b T(i,j)=(T(M,j-1)+T(M,j+1)+2*T(M-1,j))/(4+2*Bi) 4a T(i,j)=(T(M-1,N)+T(M,N-1))/(2+2*Bi) 4b T(i,j)=(T(M-1,1)+T(M,2))/(2+Bi)2.4 结果及分析肋效率的数值计算结果工况1（Bi=0.01）工况2 （Bi=1）工况 节点 Bi 二维数值计算的η值 一维数值计算的η值 相对偏差 1 9×5 0.01 0.970 0.971 0.1％ 217×510.1850.2008.1％当Bi=0.01时，得出的肋效率与一维公式计算结果的差别完全可以忽略，而对于Bi=1的长肋片，它们之间差别较大。

由此可见，判断肋片中导热可否按一维问题处理的标准应该是Bi 尽量小。

第三章 三层墙导热3.1 概述因为各层材料的导温系数不同，在分界面处的导温系数需用调和平均数计算。

用有限体积法计算时还要注意边界条件的处理。

3.2 问题描述一居民的砖墙结构如图所示：各层参数如下：3.3 TDMA 算法一维稳态不含源项的控制方程为0)(=dx dT k dx d用有限体积法可化为：0)()(=---w x TWTP kw e x TP TE keδδ参数 k （W/㎡.k ） L/㎜ N 塑料板 0.16 12 10 杨木 0.141 25 20 砖块0.72200100室内空气 TL=（25±5）℃ hL=7W/(㎡.k)室外空气 Tr=5℃hr=28W/(㎡.k)塑料板杨木一般建筑用砖块设为：a p *TP=a e *TE+a w *TW其中 a e =e x ke)(δ ；aw=w x kw )(δa p=a e +a w为使方程具一般性，可设为A i *T i =B i *T i-1+C i *T i+1+D i 对于此问题，墙壁两端均为对流换热：当i=1时，A1=hL+ke/xe ；B1=0 ；C1=ke/xe ；D1=hL*TL ；P1=C1/A1；Q1=D1/A1 当i=2--N-1时，Ai=Bi+Ci ；Bi=kw/xw ；Ci=ke/xe ；Di=01*--=i i i i P B A C Pi 11**---+=i i i i i ii P B A Q B D Q当i=N 时，Ai=hr+kw/xw ；Bi=kw/xw ；Ci=0；Di=hr*Tr ； Pi=0；11**---+=i i i i i i i P B A Q B D Q由TDMA 算法知：Ti=Qi然后由公式 T i-1=P i-1*T i +Q i-1 反向回带可依次求出其余温度值。