传热学第四章课件.ppt

(1) 理论分析方法 在理论分析的基础上，直接对微分方程在给定的定解条件下进行积分，这样获得的解称之为分析解(解析解)，或叫理论解；,第四章导热问题的数值解法,1. 求解导热问题的三种基本方法,(1) 理论分析法；(2) 数值计算法；(3) 实验法,2. 三种方法的基本求解过程,4-1 导热问题数值求解的基本思想,3. 三种方法的优缺点,(1) 分析法能获得所研究问题的精确解，可以为实验和数值计算提供比较依据；分析解具有普遍性，各种情况的影响清晰可见；局限性很大，对复杂的问题无法求解；,(2) 数值计算法 把原来在时间和空间连续的物理量的场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，从而获得离散点上被求物理量的值。这些离散点上被求物理量值的集合称为该物理量的数值解。(3) 实验法 就是在传热学基本理论的指导下，对所研究对象的传热过程进行观察、测量,(2) 数值法 在很大程度上弥补了分析法的缺点，适应性强，特别对于复杂问题更显其优越性；与实验法相比成本低(3) 实验法 是传热学的基本研究方法适应性不好；费用昂贵数值解法：有限差分法（finite-difference） 有限元法（finite-element） 边界元法（boundary- element）,4. 导热问题的数值求解基本步骤,二维矩形域内稳态无内热源，常物性的导热问题。,(1) 建立控制方程及定解条件,控制方程：,边界条件：,(2) 区域离散化,网格线：与坐标轴平行的线节 点：网格线的交点步 长：两相临节点间的距离控制容积：节点所代表的区域界面线：控制容积的边界线,(3) 区域离散化,离散方程：节点上物理量的代数方程，如tm,n,(4) 设立迭代初场,对各点物理量设置初始值,(5) 求解代数方程组,采用迭代法求解方程组,(6) 解的分析,根据温度分布,求热流,1. 建立离散方程的常用方法,(1) Taylor(泰勒)级数展开法(2) 控制容积平衡法(也称为热平衡法)(3) 控制容积积分法(4) 多项式拟合法,4-2 内节点离散方程的建立方法,2. 泰勒级数展开法,根据泰勒级数展开式，用节点(i，j)的温度ti，j来表示节点(i+1，j)的温度ti+1，j,用节点(i，j)的温度ti，j来表示节点(i-1，j)的温度ti-1，j,若取上面式右边的前三项，并将两式相加，移项整理即得二阶导数的中心差分：,截断误差 未明确写出的级数余项中的x的最低阶数为2,同样可以写出：,根据,如果x=y ，则,对于无内热源，且x=y,有内热源,无内热源,一阶、二阶导数的常用差分表达式,例：有人对一阶导数 提出了如下表达式,试判断该表达式是否正确，为什么？,将 代入 可得,答：该表达式不正确，正确式应为,(2) 控制容积平衡法(热平衡法),基本思想：对每个有限大小的控制容积应用能量守恒，从而获得温度场的代数方程组 它从基本物理现象和基本定律出发，不必事先建立控制方程，依据能量守恒和Fourier定律即可。,流入控制体的总热流量控制体内热源生成热流出控制体的总热流量控制体内能的增量,能量守恒：,即：,从所有方向流入控制体的总热流量控制体内热源生成热控制体内能的增量,稳态、无内热源时：,内部节点(m，n)：,从所有方向流入控制体的总热流量0,根据傅立叶定律,如果x=y，则,如果有内热源，且x=y ，则,4-3 边界节点离散方程的建立及代数方程的求解,1. 边界节点离散方程的建立：,(1) 平直边界上的节点,(2) 外部角点,(3) 内部角点,(1) 第二类边界条件：将qw带入边界点离散方程即可 绝 热 qw=0 非绝热 qw=const,(3) 辐射边界条件：将qw=(Tf4 T 4m,n)带入边界点离散方程,边界热流密度qw的情况,(2) 第三类边界条件：将qw=h(tftm,n)带入边界点离散方程,平直边界,外部角点,内部角点,2.节点方程组的求解,写出所有内节点和边界节点的温度差分方程，n个未知节点温度，n个代数方程式：,代数方程组的求解方法：直接解法、迭代解法,直接解法：通过有限次运算获得代数方程精确解; 矩阵求逆、高斯消元法,迭代解法：先对要计算的场作出假设、在迭代计算过程中不断予以改进、直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。,缺点：所需内存较大、方程数目多时不便、不适用于非线性问题（若物性为温度的函数，节点温度差分方程中的系数不再是常数，而是温度的函数。这些系数在计算过程中要相应地不断更新）,迭代解法有多种：简单迭代（Jacobi迭代）、高斯-赛德尔迭代、块迭代、交替方向迭代等,高斯-赛德尔迭代的特点：每次迭代时总是使用节点温度的最新值,在计算后面的节点温度时应按下式（采用最新值）,例如：根据第 k 次迭代的数值,可以求得节点温度,判断迭代是否收敛的准则：,K、k+1迭代次数；,第k次迭代得到的最大值,当有接近于零的t 时，第三个较好,允许的偏差；相对偏差一般取10-310-6,判断迭代能否收敛的准则：,系数矩阵主对角占优，即主对角的系数大于或等于其他变量系数绝对值之和。,例 如图有一矩形薄板，已知a=2b，在边界x=0和y=0处是绝热的，在x=a处给出第三类边界条件，即给定h和tf，边界y=b处壁面保持恒定温度tb。试写出各节点的离散方程。,解：对于无内热源二维稳态导热,采用均匀网格，x=y=b/2,用矩阵的形式可写为,例 有一各向同性材料的方形物体，其导热系数为常数。已知各边界的温度如图所示，求其内部网格节点的温度。,解：各节点的温度表达式为,设各节点的初始温度为,4-4 非稳态导热问题的数值解法,稳态导热,非稳态导热,本节重点： (1) 非稳态项离散的方法； (2) 扩散项离散时所取时间层的不同对计算带来的影响。,1. 一维非稳态导热时间空间区域的离散化,时间步长：从一个时间层到下一个时间层的间隔，,。,(1) 基本概念,(2) 非稳态项的离散,非稳态项的离散有三种不同的格式： 向前差分 向后差分 中心差分,向前差分 将函数t在节点(n，i+1)对点(n，i)作泰勒展开,向后差分 将函数t在节点(n，i-1)对点(n，i)作泰勒展开,中心差分 将函数t在节点(n，i+1)及(n，i-1)对点(n，i)的泰勒展开式相减,2. 一维非稳态导热微分方程的离散方法,(1) 泰勒级数展开法,一维非稳态导热微分方程中的扩散项离散与稳态导热微分方程中的方法相同，对一维非稳态导热微分方程中 扩散项中心差分 非稳态项向前差分,i时层差分,i+1时层差分,显式差分格式,隐式差分格式,显式差分格式：若已知i时层上各节点的温度值，根据该差分格式即可算出(i+1)时层上各内点的温度，而不必求解联立方程。即 是前一时刻i的第n节点及当前时刻(i+1)与n相邻两节点温度的显函数。,优点：计算工作量小；缺点：受时间及空间步长的限制。,隐式差分格式：即使已知i 时层上各节点的温度值，根据差分格式也不能直接算出(i+1)时层上各节点的温度，而必须求解(i+1)时层上的一个联立方程组，才能算出(i+1)时层各节点的温度。即 是前一时刻i的第n节点及当前时刻(i+1)与n相邻两节点温度的隐函数。,优点：不受时间及空间的步长影响；缺点：计算工作量大。,(2) 控制容积平衡法(热平衡法),优点：不受网格是否均匀限制；不受物体是否为常数限制。,如图所示，一无限大平板，右侧面受周围流体的冷却，表面传热系数为h ，对于边界节点N代表了宽为x/2的微元体,在i 时层上，从节点N-1传导给节点N的热流量,平板右侧被冷却时，在i时层上其单位面积损失的热流量,从i时刻到(i+1)时刻微元体的内能增量,根据能量守恒,整理得,整理得外边界节点的迭代式,内节点的显式迭代式,对于对称加热的中心节点，可以假想有一节点(-1)存在，则节点(1)就成为内节点，就可采用内节点迭代公式对其进行计算,3. 解的稳定性判断,要求 (含 及 )前的系数0。对显式差分格式,对流边界节点,内节点,内节点,对流边界节点,对隐式差分格式,内节点,边界节点,例 厚2=60mm的无限大平板受对称的冷却，初始温度t0=100。突然将其置于t=0的流体中。已知平板的=40w/mK，h=1000w/m2K，试其导热系数为常数。已知各边界的温度如图所示，求Fo=0.1、0.5时各节点随时间的变化。,解 取x=10mm=0.01m，则,可知,解才会稳定,各节点的温度表达式为,Fo=0.1时各节点的温度变化,Fo=0.5时各节点的温度变化,Fo=0.1时各节点的温度变化,如图，将厚度为0.05mm，导热系数为40w/mK的金属薄膜敷设到厚度为10cm，导热系数为0.04w/mK的干木板上，当给金属薄膜通电时，其发热率为2000w/m2，此时，测得tw1=50，tw2=20，tf=20，求金属薄膜表面与流过该表面空气之间的对流换热系数。设金属薄膜表面的发射率为0.2。,两层平壁内为一维稳态导热，温度分布如图所示，已知导热系数1、2均为常数，试确定：(1)热流密度q1、q2、q3的相对大小；(2) 导热系数1和2的相对大小。,从宇宙飞船伸出一根细长的散热棒，以辐射换热形式将热量散发到温度为绝对零度的外部空间，已知棒的表面发射率为，导热系数为，长度为l，横截面积为A，截面周长为P，根部温度为T0，试写出导热微分方程及边界条件。,一初始温度为t0的固体，被置于室温为tf的房间中。物体表面的发射率为，表面与空气间的表面传热系数为h。物体的体积为V，参与换热的面积为A，比热容和密度分别为c和。物体的内热阻可忽略不计，试列出物体温度随时间变化的微分方程。,如图，将厚度为0.05mm，导热系数为40w/mK的金属薄膜敷设到厚度为10cm，导热系数为0.04w/mK的干木板上，当给金属薄膜通电时，其发热率为2000w/m2，此时，测得tw1=50，tw2=20，tf=20，求金属薄膜表面与流过该表面空气之间的对流换热系数。设金属薄膜表面的发射率为0.2。,解：对金属薄膜进行能量平衡分析可知，薄膜产生的电热量q应等于薄膜上表面向外界的辐射换热量qr、与空气之间的对流换热量qc以及通过模板的导热量qd三者之和。,据题有,解得,两层平壁内为一维稳态导热，温度分布如图所示，已知导热系数1、2均为常数，试确定：(1)热流密度q1、q2、q3的相对大小；(2) 导热系数1和2的相对大小。,解：(1) 对于一维平板壁，其稳态导热微分方程的通式为,对第一层平板壁，即 区间内，温度分布曲线上凸，则,可得 即有内热源,对于稳态导热问题，根据能量守恒定律有,由此可得,对第二层平板壁，即 区间内，温度分布是直线，则,可得 即无内热源,对于稳态导热问题，根据能量守恒定律有,(2) 对于通过x=x2面上的热流密度有两种表示方法,由图可知,故,从宇宙飞船伸出一根细长的散热棒，以辐射换热形式将热量散发到温度为绝对零度的外部空间，已知棒的表面发射率为，导热系数为，长度为l，横截面积为A，截面周长为P，根部温度为T0，试写出导热微分方程及边界条件。,解：对于细长圆棒，假设温度只在杆长方向变化，则此题属于一维稳态导热问题。分析厚度为dx的微元段的导热：,微元段的净导热：,微元段的辐射散热：,根据能量守恒有：,导热微分方程为：,边界条件：,无内热源，常物性二维导热物体在某一瞬时的温度分布为 。试说明该导热物体在x=0，y=1处的温度随时间的增加是逐渐增高还是逐渐降低？,解：常物性无内热源二维物体导热的微分方程为,由某一瞬时的温度分布 可得,在x=0，y=1处的温度随时间的增加逐渐增高,一初始温度为t0的固体，被置于室温为tf的房间中。物体表面的发射率为，表面与空气间的表面传热系数为h。物体的体积为V，参与换热的面积为A，比热容和密度分别为c和。物体的内热阻可忽略不计，试列出物体温度随时间变化的微分方程。,解：该物体在加热或冷却过程中，与周围流体间存在对流换热，与周围环境间存在辐射换热。设环境温度与流体温度相等，则在单位时间内对流换热c，辐射换热r和内能增量e分别为,根据能量守恒有,微分方程为,本章小结,1. 基于有限差分的数值解法是用差商代替导数，用求解代数方程组直接得到的离散的节点温度值，代替求解导热微分方程式及相应的单值条件得到的温度与坐标和时间的函数关系。它是一种工程上非常有用的得到近似解的方法。2. 导数的有限差分表达式常用的有向前差分、向后差分和中心差分三种。一般温度对空间坐标的二阶导数采用中心差分表达式，温度对时间的一阶导数采用向前差分或向后差分表达式。因此。对于非稳态导热问题有显式差分格式和隐式差分格式。3. 节点温度的离散方程可以用泰勒级数法和热平衡法导得，两者的结果是完全一致的。对于物性随温度变化的问题，用热平衡法更为方便。应掌握用热平衡法写出内节点和边界节点的温度离散方程,4. 稳态导热问题的节点离散方程组都可以表示为At=C。求解这一线性代数方程组的常用方法是高斯赛德尔迭代法。5. 非稳态导热问题，由于应用温度对时间的导数的差分格式不同，可分为显式与隐式格式。显式格式可以直接求解；隐式格式需联立求解代数方程组，可用迭代法求解。6. 对非稳态导热问题的显式格式，其数值解的稳定性要受稳定性条件的限制，所以空间步长和时间步长的选择不是任意的。必须满足稳定性条件。隐式格式是无条件稳定的。,作业：4-3，4-5，4-9，4-10，4-15，4-16，4-24，4-25,思考题：1.节点的概念.2.向前差分, 先后差分, 中心差分的概念.3.利用能量守恒定律和傅立叶定律, 推导内点和边界点离散方程的基本方法.4.两个导热系数不同的物体紧紧贴在一起, 不计接触热阻, 如何推导接触面上节点离散方程.5.显式差分方程及稳定性判据. 6.显式差分方程和隐式差分方程在求解时的差别.,