?Biv??h
1 越小,表示内热阻越小,外部热阻越大。此时采用集总参数法求解更为合适。
2 )物理意义: 的大小反映了物体在非稳态导热条件下,物体内温度场的分布规律。 2 、 1 )定义: 表征两个时间间隔相比所得的无量纲时间。
Fov??l2a
分子 τ 是从边界上开始发生热扰动的时刻起到所计时刻为止的时间间隔。分母可视为边界上发生的有限大小的热扰动穿过一定厚度的固体层扩散到 的面积上所需的时间。
2 )物理意义:表示非稳态导热过程进行的程度, 越大,热扰动就越深入地传播到物体内部,因而物体内各点的温度越接近周围介质的温度。
§3 — 3 一维非稳态导热的分析解
本节介绍第三类边界条件下:无限大平板、无限长圆柱、球的分析解及应 用。如何理解无限大物体,如:当一块平板的长度、宽度 >> 厚度时,平板 的长度和宽度的边缘向四周的散热对平板内的温度分布影响很少,以至于可 以把平板内各点的温度看作仅是厚度的函数时,该平板就是一块“无限大” 平板。若平板的长度、宽度、厚度相差较小,但平板四周绝热良好,则热量 交换仅发生在平板两侧面,从传热的角度分析,可简化成一维导热问题。 一、 无限大平板的分析解
已知:厚度 的无限大平板,初温 ,初始瞬间将其放于温度为 的 流体中,而且 > ,流体与板面间的表面传热系数为一常数。 试确定在非稳态过程中板内的温度分布。
解:如图 3-5 所示,平板两面对称受热,所以其内温度分布以其中心截面为对称
面。
对于 x 0 的半块平板,其导热微分方程:
(0 ) ( 3-8 ) 定解条件:t(x,0)= (0 x ) (边界条件) (边界条件) 引入过余温度: 则 ( 0 (x,0)= (0 x ) (初始条件) (边界条件) (边界条件) 对偏微分方程 分离变量求解得: ( 3-10 ) 其中离散值 是下列超越方程的根,称为特征值。( 3-11 ) 其中 Bi 是以特征长度为 的毕渥数。 由此可见:平板中的无量纲过余温度 ?? 与三个无量纲数有关:以平板厚度一半 为特征长度的傅立叶数、毕渥数及 即: ( 3-12 ) 二、非稳态导热的正规状况阶段 1 、平板中任一点的过余温度与平板中心的过余温度的关系 前述得到的分析解是一个无穷级数,计算工作量大,但对比计算表明,当 Fo>0.2 时,采用该级数的第一项与采用完整的级数计算平板中心温度的误差小于 1% ,因此,当 Fo>0.2 时,采用以下简化结果: ( 3-13 ) 其中特征值 之值与 Bi 有关。 由上式( 3-13 )可知: Fo>0.2 以后平板中任一点的过余温度 (x , τ) 与 平板中心的过余温度 (0 , τ)= ( τ )之比为: ( 3-14 ) 此式反映了非稳态导热过程中一种很重要的物理现象:即当 Fo>0.2 以后,虽然 (x , τ) 与 ( τ )各自均与 τ 有关,但其比值则与 τ 无关,而仅取决于几何位置(x?)及边界条件( Bi )。也就是说,初始条件的影响已经消失,无 ?m(?)之值是一个常数,也就是无量纲论初始条件分布如何,只要 Fo>0.2 , 的温度分布是一样的。 由此可见,当 Fo>0.2 时,非稳态导热过程进入正规状况阶段。 2 、在一个时间间隔内非稳态导热过程中传递的热量 1 ) 从物体初始时刻平板与周围介质处于热平衡,这一过程中传递的热量: ?(x,?)( 3-15 ) 此值为非稳态导热过程中传递的最大热量。 2 ) 从初始时刻到某一时间 τ ,这段时间内所传递的热量 ( 3-16 ) 3 ) 之比: (3-17) : 其中: 是时刻 τ 物体的平均过余温度, 。 对于无限大平板,当 Fo>0.2 ,将式( 3-13 )代入 的定义式,可得: ( 3-18 ) 对圆柱体、球体 及 可表示为: >0.2 时,无穷级数的解也可用第一项近似代替,并且 ( 3-19 ) ( 3-30 ) 其中: η 为无量纲几何位置,对平板 ,对柱体及球体 , R 为外表面半径,系数 A 、 B 及函数 的表达式取决于几何形状,见教材表 3-2 所示。 三、正规阶段状况的实用计算方法 当 Fo>0.2 时,可采用上述计算公式求得非稳态导热物体的温度场及交换的热量,也可采用简化的拟合公式和诺模图求得。 1 、诺模图:工程技术中,为便于计算,采用按分析解的级数第一项绘制的一些图线,叫诺模图。 2 、海斯勒图:诺模图中用以确定温度分布的图线,称海斯勒图。 首先根据( 3—13 )式给出 随 Fo 及 Bi 变化的曲线(此时 x/δ=0 ), 然后根据( 3 — 14 )式确定 的值,于是平板中任意一点的 值便为: ( 3-21 ) 同样,从初始时刻到时刻 τ 物体与环境间所交换的热量,可采用( 3 — 15 )、( 3 — 17 )作出 曲线。 3 、诺模图法评述 优点:简洁方便。 缺点:准确度有限,误差较大。 目前,随着计算技术的发展,直接应用分析解及简化拟合公式计算的方法受到重视。 四、分析解应用范围的推广及讨论 1 、推广范围 1 )对物体被冷却的情况也适用; 2 )也适于一侧绝热,另一侧为第三类边界条件的厚为 δ 的平板; 3 )当固体表面与流体间的表面传热系数 h 时,即表面换热热阻 0 时,所以 时分析解就是固体表面温度发生一突然变化然后保持不变时的解,即第一类边界条件的解。 2 、讨论 Bi 与 Fo 对温度场的影响: 1 )傅立叶数 Fo : 由 (3-10) 、 (3-13) 式及诺模图可知:物体中各点的过余温度随时间 τ 的增加而减小;而 Fo 与 成正比,所以物体中各点过余温度亦随 Fo 的增大而减小。 2 )毕渥数 Bi Bi 对温度的影响从以下两方面分析: 一方面,从教材图 3 — 6 可知, Fo 相同时, Bi 越大, 越小。因为, Bi 越大,意味着固体表面的换热条件越强,导致物体的中心温度越迅速地接近周围介质的温度;当 Bi 时,意味着在过程开始瞬间物体表面温度就达到介质温度,物体中心温度变化最快,所以在诺模图中 1/Bi=0 时的线就是壁面温度保持恒定的第一类边界条件的解。 另一方面 Bi 的大小决定于物体内部温度的扯平程度。如:对于平板,从诺模图 3 — 7 中可知: 当 >10 (即 Bi<0.1 )时,截面上的过余温度差小于 5 % 当 Bi 下限一直推到 0.01 时,其分析解与集总参数法的解相差极微。 综上可得如下结论:介质温度恒定的第三类边界条件下的分析解;当 Bi 转化为第一类边界条件下的解, Bi 0 时,则与集总参数法的解相同。 时, §3 — 4 二维及三维非稳态导热问题的求解 一、求解方法 对于典型的几何形状的物体,可利用一维非稳态导热问题分析解的组合求得。 如图 3-9 所示:无限长方柱体的非稳态导热问题,属二维导热问题。截面 尺寸为: 的方柱体可视为两块厚度分别为 及 的无限大平板 垂直相交所截出的物体。 讨论的目的:找出二维温度场与两块无限大平板的温度场之间的关系。 已知:方柱体初温为 to ,初始时放于 t 流体中,表面传热系数为 h. 试求:温度场分布。 解:如图 3-9 所示,建立坐标系,由于其对称性,只研究其 截面的温度分布,截面上的温度分布由下列导热微分方程和定解条件确定:( 3-22 ) ( I ) ( II ) ( III ) ( IV ) 式中: 为无量纲过余温度。 如果无量纲过余温度 与 分别是处于与方柱体同样定解条件下的厚度分别为 及 的无限大平板的分析解,则它们必须满足各自的导热微分方程及定解条件,即: ( 3-23 )