传热学课件：第四章 数值解法.pptVIP

一、建立节点方程组 原理：把有限差商代替微商 原理：对于稳态热传导，单位时间内进入每一个节点的能量和必为零。 二、解方程组 松弛法（余数调节法） ⑤重复步骤③ ④ ，直到全部余数为零。 例题：有一正方形截面，边界长为1m，边界上的温度已知，求t1，t2，t3，t4。 迭代法 松弛法 三、利用“导热形状因子S”计算导热量 导热系数为常数，无内热源稳态导热体内，两壁温度为定值，即有Q=λS(t1－t2)。 稳定热传导的热源法（虚拟热源法） 定义：如果一个物体有内热源作用时，我们可以通过导热微分方程式和相应的单值性条件求温度分布。但如果知道温度分布，我们反过来找导致这种温度分布的原因—实际存在的热源或假想的热源。这种方法称虚拟热源法或称映象法。 因管道很长，从而可以看作是二维稳定导热 由于微分方程式和边界条件均是线性的，二个热源各自造成的温度场互不干涉。两个温度场叠加就得合成。 求y0，考虑管道表面温度均匀，所以是一等温圆。 由于是稳定导热，从而流过每一等温面的热流量是相同的 * * 第四章 热传导问题的数值解法 2. 求 解 方 程 组 。 数值解法的步骤 1. 建立节点方程组； 多维稳态导热 m,n m+1,n m,n-1 m,n+1 m-1,n m+1/2,n m-1/2,n （1）有限差分法 内 部 结 点 tm-1,n tm,n Φ1 Φ4 Φ3 Φ2 2. 能量平衡法 tm-1,n tm,n Φ1 Φ3 Φ2 tm,n＋1 tm,n－1 (a) m,n h,tf tm-1,n m,n m,n-1 m-1,n h,tf m,n h,tf (b) (c) (d) 1. 直接法：计算一次就可得出结果，适用于中小型问题。 高斯消元法 逆矩阵法 AX=B → X=A-1B 2. 间接法（迭代法）经过有限次的迭代，求出近似解，对于计算机来说，存储量较少。 高斯—赛德尔迭代法 （1）松弛法 ①设初值； ②求R1，R2，…，Rn，找Rmax；(余数) ③如设R4为最大，改变t4，使R4 ≈0，t4=t4+R4/4： ④重新计算有关节点的余数； ②一次次的直接计算t1，t2，…，tn ，注意计算tn时， tn前面的温度全部用新值代替。如知道t1后，求t2时，用t1代替原设的初值。 （2）高斯—赛德尔迭代法 ①选初值； 500℃ 100℃ 100℃ 100℃ 1 3 2 4 解（1）列节点方程式 150.04 150.08 250.08 250.16 5 150.15 150.31 250.31 250.61 4 150.61 151.18 251.26 252.35 3 152.35 154.69 254.69 259.38 2 159.38 168.75 268.75 275 1 200 200 300 300 0 t4 t3 t2 t1 n 0 150 0 150 0 250 0 250 11 5 150 0 150 -20 255 5 250 10 5 150 5 150 -15 255 -15 255 9 15 150 -35 160 -5 255 -5 255 8 -25 160 -25 160 -5 255 -5 255 7 -10 160 -25 160 -65 270 10 255 6 -10 160 -5 160 -45 270 -70 275 5 30 160 -165 200 -45 270 -30 275 4 -130 200 -125 200 -5 270 -30 275 3 -100 200 -125 200 -125 300 0 275 2 -100 200 -100 200 -100 300 -100 300 1 R4 t4 R3 t3 R2 t2 R1 t1 No t1=250 ℃ ，t2=250 ℃ ，t3=150 ℃ ，t4=150 ℃ 补充内容 有一热力管道，外径d=2r，埋于地平面下h米深处。土壤为均质且导热系数λ为常数。管子表面温度及地表面温度也是均匀的常量，为tW和tF，设管道很长，求单位管道的热损失。 p(x,y) x y r r’ r” h y0 M N tF 下面用此方法求地下埋管与土壤间的导热量 ②假设在x=0，y=y0有一线热源作用（源），管道的温度是由于此线热源作用而造成的，而在对称于地表面的位置上x=0，y=－y0有一负热源－ql W/m作用称汇。 这是一个半无限大物体的热传热问题。想法变成一个无限大物体： ①映象系统将地面往上无限延伸。 p(x,y) r r’ r” h y0 M N tF 热源单独作用时 二个热源作用的