【教学课件】第四章导热问题的数值解法.ppt

《【教学课件】第四章导热问题的数值解法.ppt》由会员分享，可在线阅读，更多相关《【教学课件】第四章导热问题的数值解法.ppt（44页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第四章 导热问题的数值解法Numerical method for heat conduction1 求解导热问题的三种基本方法：(1)理论分析法；(2)数值计算法；(3)实验法2.三种方法的基本求解过程:(1)理论分析方法:直接对微分方程在给定的定解条件下进行积分，获得解析解(close solution)(2)数值计算法:把在时间和空间连续的物理量的场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些点上的物理量值的代数方程，从而获得离散点上被求物理量的值；并称之为数值解(numerical solution)(3)实验法:在传热学基本理论的指导下，采用实验的方法对所

2、研究对象的传热过程进行实验研究，从而求得所求量的方法3 三种方法的特点三种方法的特点(1)分析法 a 能获得所研究问题的精确解，可以为实验和数值计算提供 比较依据；b 局限性很大，对复杂的问题无法求解；c 分析解具有普遍性，各种情况的影响清晰可见。(2)数值法 在很大程度上弥补了分析法的缺点，适应性 强，特别对于 复杂问题更显其优越性；与实验法相比成本低。(3)实验法 传热学的基本研究方法，a 适应性不好；b 费用昂贵。数值解法：有限差分法(finite-difference method)、有限元法(finite element method)、有限体积法(finite volume met

3、hod)、边界元法(boundary-element method)、离散元法(discrete element method)传热（物理）问题的数值求解过程传热（物理）问题的数值求解过程建立控制方程及定解条件区域离散化建立节点物理量的代数方程求解代数方程组获得数值解并分析结果4.1 有限差分法的基本原理将求解区域离散、以节点网格代替物体，以每个节点的温度作为未知量在节点上用差分代替微分，将微分方程式近似地变成差分方程式线性的代数方程组解此代数方程组，得到节点上温度的近似值1、基本思想：函数在点 的泰勒级数展开形式为：函数 、在点x的泰勒级数展开式分别为：(a)(b)(c)(d)2、函数、函数

4、 f（x）在点）在点 x 的导数的有限差分表达式：的导数的有限差分表达式：由式(a)得：由式(b)得：由式(a)与式(b)相减得：由式(a)和式(c)消去得：由式(b)和式(d)消去 得：(e)(f)(g)(i)(h)由式(a)和式(b)消去得：由(e)式(j)式分别略去以上各项得一阶、二阶导数向前、向后及中心差分公式为：、及(j)、一阶导数向前差分：一阶导数向后差分：一阶导数中心差分：二阶导数向前差分：二阶导数向后差分：二阶导数中心差分：一阶导数的差分及其误差：引入下列关系式：函数 f（x）在点 x 的一、二阶导数的有限差分表达式分别为：图4-2有限差分表达式的几何意义式中：向前和向后差分的

5、误差比中心差分的误差高，中心差分应用较广。4.2 稳态导热问题的差分表达式1。内部节点的差分方程式。内部节点的差分方程式 物理性质参数为常数的具有内热源的二维稳态导热方程：图4-3间距为x、y的矩形网格（a）将整个区域划分步长为 的矩形有限差分网格，节点 的坐标(x，y)为：pi、j为整数节点P的温度t（x，y）和热源 可表示为：在节点P，温度对x和y的二阶导数的有限差分表达式：将上式代入方程(a)中可得二维稳态导热方程的有限差分形式为：如果假定正方形网格为 ,则：物理意义：节点热平衡2 边界上节点的差分方程式边界上节点的差分方程式1.对流边界节点(i，j)边界面上的节点(i，j)满足下面的第

6、三类边界条件：图4-4对流边界节点（e）假想节点那么，在节点(i，j)处的导热方程的有限差分形式为：，再利用中心差分公式，边界条件(e)式的有限差分形式为：（f）(g)联立式(f)和式(g)，并消去 得 如果图中所示边界为绝热边界，则导热方程在节点(i，j)的有限差分形式可直接在上式中令 得到，即 2.两对流边界相交处的节点(i，j)图4-4对流边界节点由于其处于两个边界上，则其边界条件为:在节点(i，j)处导热方程的有限差分形式可写为：边界条件的有限差分形式为：(j)(k)(m)联立式(j)、式(k)和式(m)，并消去 和 得：若两个边界都是绝热的，在上式中令 得：用有限差分法求解二维稳态导

7、热问题的步骤：1 领域划分。将物体分割成很多小方块，以每个小方块的中心为节点，形成有很多节点的网络。2 列方程。根据表4-1中所列出的公式，找出对应节点的节点方程。有多少个温度未知的节点就列出多少个方程，将这些线性方程组成线性方程组。3 求解线性方程组。便得到各节点的温度值。计算精度取决于网格疏密程度。对于传热和流体力学问题的求解，一般认为差分法优于其他数值方法。4.3 线性代数方程组的求解4.3.1 直接法 高斯一约当消元法 对于n阶线性方程组 用矩阵形式表示 消元步骤：(1)首先使第一行主对角线上的元素为主元素绝对值最大的元素。如果主对角线上的元素不为主元素，那么可以利用换行的方法把主元素

8、调到主对角钱上来，使得其绝对值最大。再将第一个方程乘以分别与第i个方程相加（i2，,n）得一个新的n阶线性方程组第一列中主元素以下的各元素都化为零(2)对n-1阶线性方程组进行消元，消元法同上，经过消元以后得另一个新的n阶线性方程组(3)用同样方法进行n1次消元，最后能得到一个三角形矩阵，消元过程结束。三角形矩阵为=回代步骤从第n个方程可得再将解得的 代入第n-1个方程解出 ，再将 ，代人第n-2个方程解出 。如此不断地回代，最后可解得 。以上的消元过程及回代过程都可编成计算机通用程序。4.3.2 迭代法迭代法-逐步逼近求解方程的方法 分类：简单选代法和高斯一赛德尔迭代法 简单迭代法简单迭代法

9、：将方程改写为假定初值为 将其代入方程上式得此式为解的第一次近似，把第一次近似得到的解再代入前式，得到解的第二次近似。在已得到解的第k次近似后，代入前式的右端，得*为解的第kl次近似。只要方程组存在唯一解，且采用适当的迭代，则不论初值如何选取，当k充分大时，所得到的解必然收敛，且收敛于方程组的解。实际计算中 就结束迭代过程，而取 ()作为方程组的近似解。高斯一赛德尔迭代法高斯一赛德尔迭代法 将*式改写为 计算机只需一套单元存放或 ，节省了工作单元。同时，由于在迭代过程中尽可能地采用了新值，加快了收敛速度。例4-1 如图所示为物体的一部分。上表面绝热，左侧面为对流边界，下表面与右侧面温度给定，物

10、体的导热系数为 ，试用高斯一赛德尔迭代法求节点1、2、4、5的温度。解解 列出各节点的差分方程，并注意区域无内热源，使 。节点1 节点2 节点4节点530cm30cm1123456789其中，把各已知量代入上面各节点方程，整理得：假定各节点的初始值为 进行迭代计算，各次迭代计算结果列于表4-2。节点温度节点温度迭代次数迭代次数05.005.0020.0020.0017.759.0021.4419.6129.449.1521.6719.7139.559.2021.7419.7449.609.2221.7719.754.4 非稳态导热问题的有限差分法稳态导热：边界值问题非稳态导热：初值问题、边界值

11、问题非稳态导热过程的有限差分法：也是将导热微分方程变成差分方程。节点温度不仅和周围节点温度有关，而且还和时间有关。将时间用等间隔划分为若干等分，则当需要确定（kl）时刻的温度分布时必须知道前一时间间隔k时刻的温度分布。空间离散和时间离散二维非稳态导热微分方程式的有限差分表达式二维非稳态导热微分方程式的有限差分表达式常物性、无热源二维非稳态导热问题的微分方程式：将求解区域离散成有限个网格，从过程开始起计算，以 划分时间间隔，以k表示其序号。节点 在 时刻的温度可表示为：(4-23)ooo节点 时刻的温度对坐标的二阶导数的中心差分公式为(a)(b)对时间的一阶偏导数的向前差分可表示为(c)取 ，再

12、将式(a)、(b)、(c)代入微分方程式(4-23)得(4-24)(4-25)以准则 代入上式，得二维非稳态导热的节点差分方程 已知初始时刻各节点温度时,可由上式计算出下一时刻各节点的温度。对各个节点列式，不需解联立方程，式（4-24）称为显式格式。对于一维非稳态导热问题可用同样的方法得到其节点差分方程。(4-26)差分方程求解非稳态导热问题时收敛与稳定的必要条件：式（4-24）及式（4-26）右端第二项可能出现负值而导致温度计算值波动，出现不稳定。为避免这种异常情况出现，在选择 、和 时，必须满足下面的条件：对于一维非稳态导热问题对于二维非稳态导热问题 Why?Because:对于一维非稳态

13、导热问题若右端第二项出现负值，即 不合理边界条件的有限差分表达式边界条件的有限差分表达式 已知对流换热系数 和流体的温度 恒定，边界条件为 对于一维问题如图所示的边界点1，假设在边界点1以外再虚设一节点0，根据式（4-26）得到其差分方程：边界条件式（a）的中心差分公式为（a）（b）(c)120联立两式消去虚设节点的温度 ，并令 ，得(4-29)式（4-29）计算值稳定的条件为 即对于第二类边界条件，当q=0时，绝热边界条件的差分方程可在式（4-29）中令 ，即 得到(4-31)(4-30)(4-29)例4-2 保温炉内隔墙由耐火粘土砖构成。停炉检修时此隔墙由 800开始降温，问壁面温度降至1

14、50以下需要多少时间？隔墙厚160mm，材料的导温系数 ，导热系数 。隔墙表面由于冷空气的自然对流而冷却，对流换热系数 ，空气温度 。请列出降温2、4、6、8、10h时的温度分布。为简化计算可把隔墙视作无限大平板。分析：由于降温过程开始时墙内温度均匀，且其两侧表面的对流换热系数和空气温度均相同，所以降温过程中不同时刻墙内的温度分布均对称于中心线。一维非稳态导热问题解：选取 ,网格划分如图，则11160 mm3 2由式（4-31）、（4-26）和式（4-29）可求得不同时刻各节点的温度，即将初始值 代入上式可求得不同时刻各节点的温度，列于表4-3。由表4-3可见，10h后表面温度可冷却至150以

15、下 1点：2点：3点：表4-3计算结果08008008002.0322.6 665.6 759.60.4540.18008004.0258.0 538.2 639.50.8429.8 769.38006.0214.7 441.8 526.61.2375.8 732.9 792.88.0.181.0 366.3 434.51.6344.3 697.8 778.610.0153.6 303.4 358.0从初始时刻开始，逐步求出各点的值，不需要迭代。本章作业：4、7、8 FDM始于1910年，随着计算机的发展，它的应用范围不断扩大。对于传热学和流体力学问题，差分法依然是最有力的数值计算方法。当然，网格划分越细，计算精度越高。数值风洞差分法在空气动力学中的应用。