传热学培训教案.ppt

4-1 导热问题数值求解的基本思想及内部节点离散方程的建立 如图（a）所示二维矩形域内无内热源、稳态、常物性的导热问题采用数值解法的步骤如下： (1) 泰勒级数展开法 根据泰勒级数展开式，用节点(m,n)的温度tm,n 来表示节点(m+1,n)，而温度tm+1,n 用节点(m,n)的温度tm,n来表示节点(m-1,n)的 温度tm-1,n 2 代数方程的求解方法 2） 迭代法：先对要计算的场作出假设（设定初场），在迭代计算中不断予以改进，直到计算前的假定值与计算结果相差小于允许值为止的方法，称迭代计算收敛。 1） 直接解法：通过有限次运算获得精确解的方法，如：矩阵求解，高斯消元法。 迭代法目前应用较多的是： 1）高斯——赛德尔迭代法：每次迭代计算，均是使用节点温度的最新值。 2）用雅可比迭代法：每次迭代计算，均用上一次迭代计算出的值。 设有一三元方程组： 其中 （ i=1,2,3 ； j=1,2,3 ）及 是已知的系数（均不为零）及常数。 采用高斯——赛德尔迭代法的步骤： （1）将三元方程变形为迭式方程： （2）假设一组解（迭代初场），记为： 并代入迭代方程求得第一 次解 每次计算均用最新值代入。 （3）以新的初场重复计算，直到相邻两次迭代值之差小于允许值，则称迭代收敛，计算终止。 判断迭代是否收敛的准则： k及k+1表示迭代次数； —第k次迭代得到的最大值； 当计算区域中有接近于零的t时，第三个较好。 说明： 1 ）对于一个代数方程组，若选用的迭代方式不合适，有可能导致发散，即称迭代过程发散； 2 ）对于常物性导热问题，组成的差分方程组，迭代公式的选择应使一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值的代数和，此时，结果一定收敛。 * * 第四章　导热问题数值解基础 1 、重点内容： ① 掌握导热问题数值解法的基本思路； ② 利用热平衡法和泰勒级数展开法建立节点的离散方程。 2 、掌握内容：数值解法的实质。 3 、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。 求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，这些数值解法主要有以下几种： （1）有限差分法 （2）有限元方法 （3）边界元方法 基本思想是把连续的定解区域用有限个离散点构成的网格来代替， 这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似， 积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组 ， 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 （1）有限差分法 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式 ，借助于变分原理或加权余量法，将微分方程离散求解。 （2）有限元法 边界元法是在有限元法之后发展起来的一种较精确有效的工程数值分析方法 。又称边界积分方程-边界元法。它以定义在边界上的边界积分方程为控制方程，通过对边界分元插值离散，化为代数方程组求解。它与基于偏微分方程的区域解法相比，由于降低了问题的维数，而显著降低了自由度数，边界的离散也比区域的离散方便得多，可用较简单的单元准确地模拟边界形状，最终得到阶数较低的线性代数方程组。 （3）边界元法 分析解法与数值解法的异同点： ?? 相同点：根本目的是相同的，即确定 ① t=f(x,y,z) ；② 。 ?? 不同点：数值解法求解的是区域或时间、空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。 ??数值解法的实质 ???? 对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。该方法称为数值解法。 ???? 这些离散点上被求物理量值的集合称为该物理量的数值解。 建立控制方程及定解条件 确定节点（区域离散化） 建立节点物理量的代数方程 设立温度场的迭代初值 求解代数方程