传热讲义2009-第4章精要.ppt

第4章 导热问题的数值解法 4－2 有限差分法的基本原理 4-3 非稳态导热问题的数值解法 一厚度为0.1m的无限大平壁，两侧均为对流换热边界条件，初始时 两侧流体温度与壁内温度一致，tf1=tf2=t0=5 ℃；已知两侧对流换热系数分别为h1=11 W/m2K、h2=23W/m2K, 壁的导热系数?=0.43W/mK，导温系数a=0.3437×10-6 m2/s。如果一侧的环境温度tf1突然升高为23℃并维持不变，计算在其它参数不变的条件下，平壁内温度分布及两侧壁面热流密度随时间的变化规律（用图形表示）。 同内部节点温度方程的显式差分格式的道理一样，上式必须满足显式差分格式的稳定性条件，即 因为 ，所以只要满足上式，自然满足内部节点温度方程显式差分格式的稳定性条件。因此上式是第三类边界条件下一维非稳态导热所有节点温度方程显式差分格式的稳定性条件。 对比 如果节点0的温度对时间的变化率采用向后差分，热平衡方程式可写成 引进网格付里叶数 和网格毕渥数 上式写成显函数的形式 边界节点温度方程的显式差分格式 边界节点温度方程的隐式差分格式 2) 节点温度差分方程的求解 （1）对显式差分格式，在满足稳定性条件的前提下，从初始温度开始,依次求出每个节点下一时刻的温度,不需要不必联立求解方程组; （2）对隐式差分格式，对每一个时刻的所有节点温度，需要同时求解：高斯-塞德尔迭代 第四章小结 （2）掌握有限差分法的原理; 重点掌握以下内容： （3）能够根据导热问题的特点，合理地进行求解域的离散； (1) 理解数值解法的基本思想，熟悉数值解法的基本步骤； （4）重点掌握热平衡法建立节点温度差分方程； （5）会利用计算机求解节点温度差分方程组，并掌握求解结果表达方法。 大作业 要求 将全部计算内容（包括网格的划分、节点方程组、计算框图、程序及计算结果）用A4纸打印。 NCEPU Research Group of Heat Transfer NCEPU Research Group of Heat Transfer 4－1 数值解法的基本思想与步骤 1.基本思想： 用导热问题所涉及的空间和时间区域内有限个离散点(节点)的温度近似值来代替物体内实际连续的温度分布，将连续温度分布函数的求解问题转化为各节点温度值的求解问题，将导热微分方程的求解问题转化为节点温度代数方程的求解问题。 2. 基本步骤: （1）对实际导热问题进行分析，做必要的、合理的简化，建立符合实际的物理模型；。 （2）根据物理模型建立完整的数学模型，即给出导热微分方程和单值性条件； 第(1)、(2)步是导热问题所有求解方法的基础。 （3）求解域离散化：用与坐标轴平行的网络线将所涉及的空间和时间区域划分成有限个子区域，将网络线的交点作为节点, 每个节点就代表以它为中心的子区域（元体或称为控制容积），节点温度就代表子区域的温度； （4）建立节点温度代数方程组； （5）求解节点温度代数方程组，得到所有节点的温度值； （6）对计算结果进行分析，若不符合实际情况，则修正上述步骤，重复进行计算，直到满意为止。 目前常用的数值解法主要有：有限差分法、有限元法、边界元法等。其中有限差分法比较成熟，应用广泛。下面主要介绍有限差分法的基本原理。 建立控制方程及定解条件 确定节点（区域离散化） 建立节点物理量的代数方程 设立温度场的迭代初值 求解代数方程 是否收敛 解的分析 改进初场 是 否 以常物性、无内热源的二维稳态导热为例: 用有限差分近似微分，用有限差商近似微商。 将偏微分方程转化差分代数方程。 1. 求解域的离散化 1)子区域的划分 选择网格宽度?x、?y（步长），划分子区域。步长大小根据问题的需要而定。 2)节点的选择 选择网线和网线及网线与物体边界的交点作为节点，标定节点位置，如( m, n )、( m+1,n ) 等。 2. 节点温度差分方程的建立 两种方法：泰勒级数展开法 控制容积热平衡法。 1) 泰勒级数展开法 对节点（m+1, n）和（m-1, n）分别写出t 在(i, j)节点的泰勒级数展开式： 将上两式相加，略去高阶项，得 中心差分格式 同样可得y方向得二阶偏导数 对于无内热源的二维稳态导热，导热微分方程为 将以上两式代入，得 取?x = ?y，得 2) 热平衡法 基本思想：对每个有限大小的控制容积应用能量守恒，从而获得温度场的代数方程组，它从基本物理现象和基本