传热学讲义-2.pptVIP

前 提 上节所讨论的是一维、无内热源、常物性的非稳态导热问题。所谓“无限大”的含义,如果平板沿厚度方向四周绝热，圆柱体端部绝热，均可简化成一维问题（参见图2-5）。 本节的讨论适用于图3-1中（b），（c）两种情形，即物体内部导热热阻不能忽略的情形。 其中图3-1（b）可以看成图3-1（c）的特例。 因此，在本节中，温度是空间坐标（r或x）及时间的函数，即：t=f(t,τ)或t=(r,τ)。 本节重点讨论无限大平板的一维问题，其结论可以类推到无限长圆柱体和球体的导热。 * 物理问题及数学描写 以无限大平板为例，讨论厚为2δ的平板，处于温度为t∞、表面传热系数为h的对流环境中，初始时刻温度为t0。 则描述其温度分布的导热微分方程及定解条件为： 可以看出t=f(a,τ,λ,δ,h,x)，即温度取决于上述6个变量，若引入如下无量纲温度Θ、无量纲坐标X和无量纳时间F0。 则(3—7)式可以化成： 由式(3—8)，当引入上述无量纲量后，无量纲温度只取决于下述3个无量纲量： 即引入上述参数后，变量个数大大减少。 * 求 解 结 果 利用分离变量法，对式(3—7)进行求解的结果为 式中是下列超越方程的根： 上述分析求解结果和无量纲化结果表明，引入准则数Bi，F0及无量纲量Θ和X后，变量个数由原来的7个减少至4个。 对本例而言，意味着用图的形式表示解的结果是有可能的，这就是后面要介绍的海斯勒图。 * 正规状况阶段非稳态导热的计算方法 分析解： 当F00.2可直接利用教材式(3－22)，(3－25)，(3－26)，计算无限大平板内瞬时温度分布或非稳态导热过程中所传递的热量，即： 式中ˉθ＝ˉθ（τ），是时刻τ，物体的平均过余温度。 近似拟合公式 可利用教材式(3－27)，(3－28)等近似计算公式计算无限大平板、无限长圆柱、球的非稳态导热。即： 式中各系数A，B及f(μ1η)的表达式详见教材表3－2。，A，B和J0(x)也可采用近似拟合公式，见教材式(3－29a)～(3－29d)和表3－3、表3－4。 海斯勒线算图 由式(3－10)，Θ＝θ/θ0＝f(Bi,F0,x/δ)，考虑到一张图上只能表示三个参数，所以采用两张图的形式(以无限大平板为例，即教材图3－6，3－7)。教材图3－6中，x/δ ＝0，θm/θ0f(F0,Bi) 。教材图3－7中，F00.2，即正规状况阶段，θ/θm＝cos(μ1x/δ)＝f（x/δ，Bi） 。 利用海斯勒图仍可计算热量，如教材图3－8所示。 * 二维及三维非稳态导热问题的求解 乘积解的理论依据及解的形式 教材已通过方程化简的形式，分别从微分方程及边界条件、初始条件的角度证明了乘积解的存在。满足乘积解的多维非稳态导热问题可以分解为相应的二个或三个一维问题解的乘积形式。 以直角坐标系下的三维问题为例，解的形式为： 或 注意乘积解中温度必须以过余温度或无量纲过余温度的形式出现，而以温度t出现时则不满足乘积解。 乘积解的适用条件 教材所介绍的乘积解是有条件的，即要求初始温度为常数，第一类边界条件中边界温度为定值和第三类边界条件中周围流体温度和h为定值时的情形。 实际求解非稳态导热问题时，常常要判断温度变化最快或最慢点所在的位置。 判断的方法是：温度变化最慢的点位于物体的体心或形心，而变化最快的点则位于离体心或形心最远的点。 * 半无限大物体的非稳态导热 物理问题，数学描写及解的结果 初始温度均匀，第一类边界条件下半无限大物体的非稳态导热问题的数学描写如下： 解的结果用误差函数表示： * 对解的讨论 由于误差函数erf（η）随η的增加而增加，且erf(2)=0．995 3。所以从工程观点来看，可将η≥2作为是否受扰动的标志。 由于η=x/√(4aτ)所以对厚为2δ的平板 从几何位置来看：当δ≥4√（aτ），则在时刻以前可将平板看成半无限大物体。 从时间角度来看：τ≤x2/(16a)时x处的温度仍可认为等于初始温度t∞。其中x2/(16a)称为惰性时间。 半无限大物体吸热量或放热量的计算  其中qw为表面上的热流密度，Q为[0，]时间间隔内通过面积A的总热量。为吸热系数。 对有限大小的物体，半无限大的概念只适用于非稳态导热的初始阶段。 因此，对半无限大物体而言，不存在前面所讲的正规状况阶段。 * 半无限大物体的非稳态导热 非稳态导热的基本概念 集总参数法的简化分析 一维非稳态导热的分析解 二维及三维非稳态导热问题的求解 半无限大物体的非稳态导热 * 导热问题的数值解法 导热问题数值求解的基本思想及内节点离散方程的建立 边界节点离散方程的建立及代数方程的求解 非稳态导热问题的数值解法 导热问题数值计算实例 * 对流换热 对流换热概说