第4章导热问题的数值解法课件.pptVIP

第四章　热传导问题的数值解法 计算传热学商业软件简介（CFD/NHT） 网格生成质量对计算精度与稳定性影响极大。 在几何形状复杂的区域上要生成好的网格相当困难 超过90％的精力要用在生成合适的网格上 任何CFD计算结果都要进行检验 应用计算流体力学理论 应用实验数据 其他CFD数据 高质量的CFD模拟结果是建立在对流体机理的正确理解 导热部分小结 x y qw 边界节点 (m,n) 只代表半个元体，若边界上有向该元体传递的热流密度为qw ，据能量守恒定律： 4.3.1 边界节点离散方程的建立 (1) 平直边界上的节点 (2) 外部角点 如图所示，二维墙角计算区域中，该节点外角点仅代表 1/4 个以 为边长的元体。假设边界上有向该元体传递的热流密度为 ，则据能量守恒定律得其热平衡式为： (3) 内部角点 内部角点代表了 3/4 个元体，在同样的假设条件下 x y qw 讨论关于边界热流密度的三种情况： （1）绝热边界 即令上式 即可。 （2） 值不为零 （3）对流边界 此时 ，将此表达式代入上述方程，并将此项中的 与等号前的 合并。 对于 的情形有： 流入元体， 取正，流出元体， 取负 （a）平直边界 （b）外部角点 （c）内部角点 4.3.2 处理不规则区域的阶梯型逼进法 当计算区域出现曲线边界或倾斜边界时，常常采用阶梯形的折线来模拟真实边界，然后用上述方法建立边界节点的离散方程。 4.3.3 代数方程的求解方法 2）迭代法：先对要计算的场作出假设（设定初场），在迭代计算中不断予以改进，直到计算前的假定值与计算结果相差小于允许值为止的方法，称迭代计算收敛。 1）直接解法：通过有限次运算获得精确解的方法，如：矩阵求解，高斯消元法。 2 迭代法目前应用较多的是： 1 ）雅可比迭代法（简单迭代）：每次迭代计算，均用上一次迭代计算出的值。 2 ）高斯——赛德尔迭代法：每次迭代计算，均是使用节点温度的最新值。 在计算后面的节点温度时应按下式（采用最新值） 例如：根据第 k 次迭代的数值 可以求得节点温度： 设有一三元方程组： 其中 （ i=1,2,3 ； j=1,2,3 ）及 是已知的系数（均不为零）及常数。 采用高斯——赛德尔迭代法的步骤： （1）将三元方程变形为迭式方程： （2）假设一组解（迭代初场），记为： 并代入迭代方程求得第一 次 解 每次计算均用最新值代入。 （3）以新的初场重复计算，直到相邻两次迭代值之差小于允许值，则称迭代收敛，计算终止。 判断迭代是否收敛的准则： k及k+1表示迭代次数； —第k次迭代得到的最大值 当有接近于零的t 时，第三个较好 迭代能否收敛的判据 1 ）对于一个代数方程组，若选用的迭代方式不合适，有可能导致发散，即称迭代过程发散； 2 ）对于常物性导热问题，组成的差分方程组，迭代公式的选择应使一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值的代数和，此时，结果一定收敛。 3 ）采用热平衡法导出差分方程时，若每一个方程都选用导出该方程中心节点的温度作为迭代变量，则上述条件必满足，迭代一定收敛。 这一条件数学上称主对角线占优（对角占优）； §4-4 非稳态导热问题的数值求解 在非稳态导热问题中，不但需要对空间区域进行离散，还需要对时间变量进行离散，接下来以一个一维非稳态导热问题为例，重点介绍对非稳态项的离散方法，以及不同离散方法对计算带来的影响等。 一、物理问题和数学描述 厚度 2? 的无限大平壁，?、a为已知常数，?=0时温度为 t0，突然将其放置于两侧温度为 t?并保持不变的流体中，两侧表面与流体之间的表面传热系数为h。 2δ h, t∞ h, t∞ 二、离散方程的建立 0 x t n i n,i n-1,i n+1,i n,i+1 n,i-1 1.热平衡法 对于节点 n 在 i 时刻的所代表的控制容积建立热平衡关系。 从左右两侧导入的热量等于单位时间容积的内能增量。 ①内节点 0 x t n i n,i n-1,i n+1,i n,i+1 n,i-1 ②右端边界节点 ③左端对称 2.直接用差分代替微分 ①向前差分(forward difference) 0 x t n i n,i n-1,i n+1,i n,i+1 n,i-1 ②向后差分(backwar