高等工程热力学与传热学.docVIP

大平板非稳态传热数值解与解析解对比 摘要 关键词 Abstract Key words 1问题描述 一块大平板（0≤x≤L），初始温度为常数ti，且等于环境温度tf，τ0时，x=0的边界接受恒热流qw，x=L的边界向环境放热，对流换热系数为α，试用有限差分法和解析法确定t(x, τ)，并将数值解和解析解对比之。 2数学模型 由上述所描述的问题可知，该大平板传热属于无内热源、常物性一维温度场的瞬态导热问题，故其在直角坐标系下的导热微分方程可表示如下： （2-1） 其中，t——温度，是时间与空间的函数，单位：K； τ——时间，单位：s； k——大平板的热传导系数，单位： ρ——大平板的密度，单位： c——大平板的比热容，单位： x——大平板上的位置，单位：m。 又根据已知条件可得以下三个初始条件，即： （1）当τ=0时，整个大平板温度均匀，即； （2）当τ0，x=0时，； （3）当τ0，x=L时，。 引入过于温度，，可得数学模型为： （2-2） 3近似积分法解 根据已知问题可知，该大平板传热可分为二个阶段：一、左侧x=0的边界接受热源的加热，其热量从左往右依次渗透。但当τ﹤τ，而并没有传递至最右侧x=L处的边界时，由于右侧边界的初始温度与环境相同，故其并不对环境放热；二、当ττL，即左侧热量已传递至右侧时，大平板的左侧仍接受热源的加热，而其右侧则对环境放热，直到放热量等于吸热量为止，整个大平板的温度趋于稳定。 3.1τ﹤τ （3-1） 此时边界条件有： （3-2） 将边界条件带入式（3-1）中，可确定相应的待定系数，于是可得到如下结果： （3-3） 此阶段的温度分布取决于的确定，于是，根据积分方程有： （3-4） 其中： （3-5） （3-6） 将式（3-5）、（3-6）带入式（3-4）可得： （3-7） 当时，可求得： （3-8） 再把式（3-7）代入式（3-3）中，即可得出此阶段平板内温度分布为： （3-9） 3.2ττL，即大平板已整体渗透阶段 假设此阶段大平板内温度分布为： (3-10) 此时边界条件有： (3-11) 将边界条件带入式（3-10）中，可确定相应的待定系数，于是可得到如下结果： (3-12) 此阶段的温度分布取决于的确定，于是，根据积分方程有： （3-13） 其中， （3-14） （3-15） 将式（3-14）、（3-15）带入式（3-13）可得： （3-16） 又因为整个热渗透过程是连续的，故当，时，，即 则可得： （3-17） 将式（3-17）带入式（3-16）中，可得： （3-18） 再将式（3-18）带入式（3-12）中，即可得出此阶段大平板内温度分布为： （3-19） 当时，可得，从而可得： （3-20） 4有限差分法解 4.1有限差分法的基本概念 设函数可以展开为泰勒展开式： （4-1） （4-2） 对以上两式进行改写，可得： （4-3） （4-4） 如果将式（4-3）和式（4-4）相加，即可得到一阶导数的中心差分公式： （4-5） 如果将式（4-3）和式（4-4）相减，即可得到二阶导数的中心差分共识： （4-6） 此外，还可以有朝前差分公式和朝后差分公式： （4-7） （4-8） 有限差分的核心思想就是用各阶差分代替各阶微分。 4.2非稳态传