高等工程热力学与传热学

《高等工程热力学与传热学》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、课程名称：高等工程热力学与传热学授课老师：年级专业：学生姓名：学号：导师：报告时间：大平板非稳态传热数值解与解析解对比摘要关键词AbstractKeywords1问题描述一块大平板（0≤x≤L），初始温度为常数ti，且等于环境温度tf，τ>0时，x=0的边界接受恒热流qw，x=L的边界向环境放热，对流换热系数为α，试用有限差分法和解析法确定t(x,τ)，并将数值解和解析解对比之。2数学模型由上述所描述的问题可知，该大平板传热属于无内热源、常物性一维温度场的瞬态导热问题，故其在直角坐标系下的导热微分方程可表示如下：（2-1）其中，t——温度，是时间与空间的函数，单位：K；τ——

2、时间，单位：s；k——大平板的热传导系数，单位：ρ——大平板的密度，单位：c——大平板的比热容，单位：x——大平板上的位置，单位：m。又根据已知条件可得以下三个初始条件，即：（1）当τ=0时，整个大平板温度均匀，即；（2）当τ>0，x=0时，；（3）当τ>0，x=L时，。引入过于温度，，可得数学模型为：（2-2）3近似积分法解根据已知问题可知，该大平板传热可分为二个阶段：一、左侧x=0的边界接受热源的加热，其热量从左往右依次渗透。但当τ﹤τL，即热量只传递到大平板的某一处，假设，而并没有传递至最右侧x=L处的边界时，由于右侧边界的初始温度与环境相同，故其并不对环境放热；二、当

3、τ>τL，即左侧热量已传递至右侧时，大平板的左侧仍接受热源的加热，而其右侧则对环境放热，直到放热量等于吸热量为止，整个大平板的温度趋于稳定。3.1τ﹤τL，即尚未整体渗透阶段假设此阶段大平板内温度分布为：（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τ>τ

4、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）

5、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非稳态传热分析针对式（2-2）所列的数学模型，此处采用隐式差分法（即朝后差分格式）对其进行分析，可得：（4-9）（4-10）把式（4-9）和式（4-10）代入导热微分方程可得

6、：（4-11）令，经整理可得：（4-12）根据能量守恒定律，对左、右边界条件进行隐式差分，可得：（4-13）整理后可得：（4-14）因此，式（2-2）所列的数学模型可离散为下式矩阵形式：（4-15）上式三角矩阵可通过LU分解，分解为一个上三角矩阵和一个下三角矩阵，再通过追赶法求解出各未知量。5大平板非稳态传热的MATLAB解通过上述分析，已经对大平板非稳态传热解析解以及数值解的求解过程有了一定的了解，此处就通过实例来比较解析解与数值解之间的差异。5.1参数定义表5.1大平板参数值(材料选用铝)项目密度比热容长度导热系数环境温度热流密度对流换热系数符号ρcLktfqwα数值27

7、0088012203001000100单位K5.2解析解的MATLAB实现参考文献[1].李格升,郭庆平与刘春晓,大平板第一类边界条件下非稳态传热近似解析解.武汉理工大学学报（交通科学与工程版）,2002:第432-434页.[2].李灿,高彦栋与黄素逸,热传导问题的MATLAB数值计算.华中科技大学学报(自然科学版),2002(9):第91-93页.[3].王飞与裴永祥,有限差分方法的MATLAB编程.新疆师范大学学报(自然科学版),2003(4):第22-27页.[4].田兵,用MATLAB解偏微分