3-1 非稳态导热的基本概念1. 定义：物体的温度随时间而变化的导热过程称为非稳态导热， 2. 分类第三章非稳态导热第三章非稳态导热 瞬态非稳态导热：物体的温度随时间的推移逐渐趋近于恒定值；如钢坯在炉内的加热 周期性非稳态导热：物体的温度随时间做周期性的变化；如室式热处理炉炉壁的导热1着重讨论瞬态非稳态导热3. 温度分布：24. 两个不同的阶段 非正规状况阶段(不规则情况阶段) 正规状况阶段(正常情况阶段)温度分布主要取决于边界温度分布主要取决于边界条件及物性条件及物性温度分布主要受初始温度温度分布主要受初始温度分布控制分布控制 非非稳稳态态导导热热过过程程总总会会经经历历：非非稳稳态态导导热热非非正正规规状状况况阶阶段段（起始阶段）、正规状况阶段、新的稳态（起始阶段）、正规状况阶段、新的稳态35. 热量变化1 板左侧导入的热流量2 板右侧导出的热流量120046. 学习非稳态导热的目的(2) 非稳态导热的导热微分方程式(3) 求解方法分 析 解 法：分离变量法、积分变换、拉普拉斯变换近似分析法：集总参数法、积分法数 值 解 法：有限差分法、蒙特卡洛法、有限元法、 分子动力学模拟(1) 温度分布和热流量分布随时间和空间的变化规律高等数学高等数学数值分析数值分析数理方程数理方程57. 毕渥数本章以第三类边界条件为重点。(1) 问题的分析 如图所示，存在两个换热环节： 流体与物体表面的对流换热 物体内部的导热tfhtfhxt 0 tfhxt 0(2) 毕渥数的定义：物理意义：固体内部导热热阻与其界面上换热热阻之比6无量纲数无量纲数当Bi时，r rh ；因此，可以忽略对流换热热阻当Bi0 时，r t,固体与流体间的表面传热系数h，固体的物性参数均保持常数。 求：根据集总参数法确定物体温度随时间的依变关系 9解： 建立非稳态导热数学模型 方法一：椐非稳态有内热源的导热微分方程： 物体内部导热热阻很小，忽略不计。 物体温度在同一瞬间各点温度基本相等，即t仅是的一元函数，而与坐标x、y、z无关，即： 可视为广义热源，而且热交换的边界不是计算边界（零维无任何边界）。 则有：10 界面上交换的热量应折算成整个物体的体积热源，即： tot物体被冷却， 应为负值 由（a），（b）式得： 这就是瞬时时刻导热微分方程式。方法二：根据能量守恒原理，建立物体的热平衡方程，即物体与环境的对流散热量=物体内能的减少量.11 物体温度随时间的依变关系方程式改写为：令= tt，则有初始条件初始条件控制方程控制方程12 积分 其中的指数：13第三章 非稳态导热Fov是傅立叶数物体中的温度呈指数分布方程中指数的量纲：14即与 的量纲相同，当 时，则此时上式表明：当传热时间等于 时，物体的过余温度 已经达到了初始过余温度的36.8。称 为时间常数，用 表示。1536.8%16 如果导热体的热容量(cV)小、换热条件好(h大)，那么单位时间所传递的热量大、导热体的温度变化快，时间常数 ( cV / hA) 小。 对于测温的热电偶节点，时间常数越小、说明热电偶对流体温度变化的响应越快。这是测温技术所需要的(微细热电偶、薄膜热电阻)工程上认为=4 cV / hA时，导热体已达到热平衡状态173. 瞬态热流量导热体在时间 0 内传给流体的总热量： 当物体被加热时当物体被加热时(tt，流体与板面间的表面传热系数h为一常数。 求：在非稳态过程中板内的温度分布。27平板非稳态导热微分方程：导热微分方程( (对称性对称性) )边界条件初始条件28引入变量过余温度令方程可化为：29用分离变量法可得其解析解为：式中n为超越方程 的根2. 无限长圆柱的分析解已知：半径为R的实心圆柱，初温为t0，初始瞬间将其放于温度为t的流体中，而且t0t，流体与板面间的表面传热系数h为一常数。 求：在非稳态过程中圆柱内的温度分布。30圆柱非稳态导热微分方程：导热微分方程边界条件初始条件31式中n为超越方程 的根用分离变量法可得其解析解为：3. 圆球的分析解已知：半径为R的实心球，初温为t0，初始瞬间将其放于温度为t的流体中，而且t0t，流体与板面间的表面传热系数h为一常数。 求：在非稳态过程中圆球的温度分布。第一类 v 阶贝塞尔函数32圆球非稳态导热微分方程：导热微分方程边界条件初始条件33式中n为超越方程 的根用分离变量法可得其解析解为：对于一维非稳态导热问题，其解均可表述为344. 非稳态导热正规状况阶段解的简化 数值计算表明，当Fo数大于0.2时，取无穷级数的首项，其计算误差小于1%对于平板与时间无关与时间无关35与时间无关与时间无关对于圆柱对于圆球与时间无关与时间无关36若令Q为0，内所传递热量考察热量的传递Q0 非稳态导热所能传递的最大热量5. 非稳态导热的热量计算时刻 的平均过余温度对于平板37此处此处的A，B及函数f(1)见P127表3-1 及 可用一通式表达无限大平板长圆柱体及球对无限大平板，长圆柱体及球386. 正规热状况的实用计算方法拟合公式法对上述公式中的A，B， 1 ，J0 可用下式拟合式中常数a，b，c见P128表3-2 a，b，c，d见P128表3-3397. 正规热状况的实用计算方法线算图法诺谟图：工程技术中，为便于计算，采用按分析解的级数第一项绘制的一些图线，叫诺模图。 三个变量，因此，需要分开来画以无限大平板为例，F00.2 时，取其级数首项即可(1) 先画海斯勒图：诺模图中用以确定温度分布的图线，称海斯勒图。 40(2) 再根据公式(3-28) 绘制其线算图(3) 于是，平板中任一点的温度为同理，非稳态换热过程所交换的热量也可以利用(331)和(333)绘制出。解的应用范围 (1) 第三类边界条件 (2) 一侧绝热，另一侧为第三类边界条件 (3) 两侧均为第一类边界条件，且维持相同的温度 F00.241例：一火箭发动机喷管，壁厚为9mm，初始温度为30。在进行静推力试验时，温度为1750的高温燃气送入该喷管，燃气与壁面间的表面传热系数 为 1950W/(m2K)。 喷 管 材 料 的 密 度=8400kg/m3， 导 热 系 数=24.6W/(mK)，c=560J/(kgK)。假设喷管因直径与厚度之比较大而可视为平壁，且外侧可作绝热处理，试确定：(1) 为使喷管的最高温度不超过材料允许温度(1000)而允许的运行时间；(2) 在所允许时间的终了时刻，壁面中的最大温差；(3) 在上述时刻壁面中的平均温度梯度与最大温度梯度。42解：本题可视为为厚度为2=29mm的平板两侧突然受第三类边界条件时的非稳态导热问题。先判断Bi数。不满足集总参数法条件，采用拟合公式计算43根据过余温度的表达式根据壁面中心温度44在上述时刻壁面的平均温度梯度由于外表面绝热，温度梯度为零，故最大温度梯度在内壁面。根据第三类边界条件45例例 一直径为100mm，长为1000mm的圆钢，初始温度为30，今将其置于1300的加热炉中。求加热半小时后圆钢表面温度与中心温度的差。取表面传热系数的平均值为232w/(m2K)， 圆 钢 的 导 热 系 数 为 40.5w/(mK)， 热 扩 散 率 为0.62510-5m2/s。解：解：对于2R=100mm的无限长圆柱，有采用拟合公式计算46过余温度比表面和中心温度为473-4 简单几何形状物体多维非稳态导热的分析解 考察一无限长方柱体(其截面为2122的长方形)，初始温度为t0，突然将其置于温度为t的流体中，表面传热系数为h，求其温度分布。2 2 12 2 2xy048 如果能用两个一维无限大平板的解来表示该二维问题的解，则可使问题变得简单其中其中及无限长方柱体的解可以表示为两个无限大平板的解的乘积，即49假设(x, y, )=x(x, )y(y, )成立，则看控制方程控制方程成立控制方程成立左侧：右侧：一维控制方程50看初始条件看边界条件同理51 证明了x(x，)y(y，)是无限长方柱体导热微分方程的解，即(x, y, )=x(x, )y(y, )。这样便可用一维无限大平壁公式、诺谟图或拟合函数求解二维导热问题初始条件和边界条件成立初始条件和边界条件成立适用条件： (1) 一侧绝热，另一侧三类；(2) 两侧均为一类；(3) 初始温度分布必须为常数；(4) 无内热源52 如果用x(x，)P表示无限大平板的解，用x(x，)C 表示无限长圆柱的解，用x(x，)R表示圆球的解，则对于如下的规则图形可球出其解。(y,)P2(x,)P1(x, y, )(x,)P1 (y,)P2(r,)C(R,)P(r, )(r,)C (R,)R(r,)C(x,)P(r, x, )(r,)C (x,)P(x, y, z, )(x,)P1 (y,)P2 (z,)P353导热量的计算二维问题三维问题54例例4 一直径为600mm，长为1000mm的钢锭，初始温度为30，今将其置于1300的加热炉中。求加热4小时后钢锭中心的温度。取表面传热系数的平均值为232w/(m2K)，钢锭的导热系数为40.5w/(mK)，热扩散率为0.62510-5m2/s。解：解：对厚度为2=1000mm的无限大平板采用拟合公式计算55过余温度比对于2R=600mm的无限长圆柱采用拟合公式计算56短圆柱的中心温度为如果将钢锭看成是无限长圆柱说明短圆柱比无限长圆柱加热的快，为什么？说明短圆柱比无限长圆柱加热的快，为什么？573-5 半无限大的物体定义：几何上是指从x=0的界面开始可以向正的 x方向及其他两个坐标（x，y）方向无限延伸的物体，称为半无限大物体。 实际中不存在该物体，但研究物体中非稳态导热的初始阶段，可把实物看为该物体处理。如：有限厚度的平板，起初有均匀温度，后其侧表面突然受到热扰动，如 (1) 壁温突然升高到一定值并保持不变 (2) 壁面突然受到恒定的热流量密度加热 (3) 壁面受到温度恒定的流体的加热或冷却 当扰动的影响只局限在表面附近，而尚未进入平板内部时，就可视该平板为，“半无限大”物体。58xtt0tw 误差函数误差函数 无量纲变量无量纲变量引入过余温度问题的解为如图所示：已知半无限大物体初始温度均匀为to，当=0时，x=0侧表面温度突然升高到tw，并保持不变，试确定物体内温度随时间的变化和在时间间隔0，内的热流量。第一类边界条件第一类边界条件59第二类边界条件第二类边界条件余误差函数余误差函数第三类边界条件第三类边界条件60误差函数：时，该x处的温度仍等于初始温度61 几何位置 若 则 时刻x处的温度可认为未变化。对一原为2的平板，若 即可作为半无限大物体来处理 时间 若 ，即 则此时x处的温度可认为完全不变。对于有限大的实际物体，半无限大物体的概念只适用于物体的非稳态导热的初始阶段。 两个重要参数:62在0，时间内，物体任一点的热密度：令 即得边界面上的热流通量0,时间内，流过面积A的总热量吸热系数吸热系数63例例5 某地地下1m处埋有自来水管，冬天来临时，地表温度初始为10，后突然受冷空气侵袭，地表温度下降到-15，并维持45天不变，试确定此种条件下45天后自来水管是否会被冻。已知土壤的物性为c=1840J/(kgK)，=2050kg/m3，=0.52w/(mK)。解：解：根据已知条件，计算导温系数对于第一类边界条件下的无限大非稳态导热，有64查书P572附录15则另外还可按书另外还可按书P136例例3-7求解求解答：地下1m处的自来水管不会被冻65例 铜和混凝土在温度为23的房间中放置很久。用手分别接触这两种材料，问哪一种材料使人感觉更凉？已知材料的物性参数为：铜，=401W/(mK)，=8933kg/m3，c=386J/(kgK)；混凝土，=1.4W/(mK)，=2300kg/m3，c=880J/(kgK)解 手感觉凉热与其感受到的热流密度有关，把铜和混凝土看成是半无限大物体，手在接触它们时发生了非稳态导热。根据物体非稳态导热表面热流密度的计算公式66随堂测验：材料相同、初温相同且满足集总参数法条件的金属薄板、细圆柱体和小球，置于同一介质中加热，若薄板厚度与细圆柱体直径和小球直径相等，问当它们被加热到相同温度时所需时间之比。解：假设板厚为2，根据题意有D柱=D球= 2，其温度分布为加热到相同温度时，应有即67作业：作业：3-2，3-5，3-12，3-16，3-21，3-24，3-32，3-42，3-48，3-5468思考题：思考题：非稳态导热的分类及各类型的特点。Bi 准则数, Fo准则数的定义及物理意义。Bi0 和Bi 各代表什么样的换热条件?集总参数法的物理意义及应用条件。使用集总参数法，物体内部温度变化及换热量的计算方法。时间常数的定义及物理意义.非稳态导热的正规状况阶段的物理意义及数学计算上的特。非稳态导热的正规状况阶段的判断条件。无限大平板和半无限大平板的物理概念。半无限大平板的概念如何应用在实际工程问题中。699 如何用查图法计算无限大平板非稳态导热正规状况阶段的换热问题?10如何用近似拟合公式法计算无限大平板非稳态导热问题?11半无限大平板非稳态导热的计算方法。70