传热与流体流动的数值计算(45章).pptVIP

传热与流体流动的数值计算 [美] S.V. 帕坦卡 著 同济大学机械工程学院 朱 彤 第四章 热传导 4-1 本章的对象 着手构建一个求解通用微分方程的数值方法 构成一个求解通用微分方程的数值方法，略去对流项。 其他一些物理过程也由非常类似于热传导方程的数学方程所控制。 本章完成了随后几章所需要的若干预备性的工作，提出代数方程的求解方法。 4-2 一维稳态热传导 基本方程 稳态一维问题的控制微分方程： 推导出离散化方程 -网格间距 网格点距离(δx)e与(δx)w没有必要相等。 虽然只有在网格相当细时才可能得到精确的解，但是在因变量随x变化相当慢的区域没有必要采用细的网格；在T~x变化较陡的区域则需要细的网格。 误区：不均匀网格的准确度比均匀网格差。 设计一个合适的非均匀网格： 从解的定性预计得到指导。 用初步粗网格的解求得T~x变化形式；然后构成合适的非均匀网格。 先进行预备性的实验或探索性试验，然后应用得到的数据资料确定在最终的实验中所应安装的探头位置和数目。 -界面导热系数ke 最直截了当的方法是假设k在P点和E值间呈线性变化： 讨论这样一种情况：围绕着网格点P的控制容积由具有均匀导热系数kP的材料填满，围绕着E点的控制容积由导热系数kE的材料填满，对于P点和E点之间的组合板，根据稳态无内热源一维导热的分析，有： 应用于系数的定义式，得到aE： -非线性 即便是在热传导问题中我们也经常遇到非线性的情况。如离散化方程中的系数本身与T有关。我们用迭代的方法来处理。过程包括： 一开始在所有各个网格点上，猜测或估计一个T值。 由这些估计的T值，计算出离散化方程中的系数的试探值。 解名义上的线性化方程组，得到一组新的T值。 以这些T值作为较好的估计值，返回到第二步并重复整个过程，直到这种进一步的重复计算（迭代）不再引起T值任何有意义的变化为止。 这种最终不变的状态叫做迭代的收敛。与之相反，迭代永远也不会收敛到一个解的状态称为发散。 -源项的线性化 当源项S与T有关时，用方程（4.4）给出的线性形式表达。 当S是T的一个非线性函数时，必须把它线性化，即规定SC和SP的值。有很多方法可以把给定的S表达式分解成SC和SP。如： 已知S=4-5T3。某些可能的线性化： 1. SC =4-5Tp*3，Sp=0。这种做法不能很好利用已知S~T关系的有利条件。 2. SC =4， Sp= -Tp*3 。看起来像准确的线性化，但已知的曲线比这一关系所反映的曲线要陡。 3. 推荐的方法： 在点Tp* ，所 选择的直线与 S~T曲线相切。 4. SC =4+20Tp*3，Sp= -25Tp*2。这一线性化比已知的S~T曲线陡，使收敛速度降低。 四种可能的线性化与实际曲线比较如图： -边界条件 讨论图中所示网格点组。 在两个边界上各有一个网格 点。其余网格点称为内点。 围绕每个内点有一个控制容 积。对每一个控制容积可以 写一个像方程（4.2）那样的离散化方程，如果看作是关于Tp的方程，那么就有了对所有内网格点上未知温度所必要的方程。其中有两个方程包含着边界网格点上的温度。通过处理这些边界温度，就把已知的边界条件引入到数值解法中。 热传导问题中有三类典型边界条件，对于每一种有： 1. 已知边界温度。此时不需要外加 任何方程。 2. 已知边界热流密度。得到： Boundary conditions 如果边界上的热流密度qB已知，则要求的对TB的方程变成： 3. 通过放热系数和周围流体的温度来规定边界的热流密度。如果热流密度qB是放热系数h以及环境流体温度Tt规定，那么，方程TB方程变为： -线性代数方程的解 一维离散化方程的解可以用标准的高斯消去法得到。 当写这些方程的系数矩阵时，所有的非零系数均排列在矩阵的三角对角线上，这种算法称为TDMA（三对角矩阵算法）。 设网格点标号为1，2，3，…，N，其中1和N代表边界点。有： 边界温度已知时，对边界点的方程只剩下一个无意义的形式。(如c1=0,bN=0) T2可以用T3的一个关系式表示，….， TN 可以由TN+1表示，回代就是TDMA的要点。 TDMA TDMA TDMA TDMA 算法 计算：P1=b1/a1 and Q1=d1/a1 使用迭代关系式，获得 Pi 和 Qi ，i=2,3…N. 设 TN =QN 使用迭代关系式，得到 Ti =Pi Ti+1 +Qi ， i=N-1, N-2…3, 2, 1 从而依次得到 TN-1, TN-2, …T3, T2, T1. 4-3 非稳态一维热传导 通用的离散化方程 时间是一个单向坐标，由一已知的初始温度分布开始，沿着时