高等工程热力学与传热学.doc
课程名称： 高等工程热力学与传热学 授课老师： 年级专业： 学生姓名： 学 号： 导 师： 报告时间： 大平板非稳态传热数值解与解析解比照摘要关键词AbstractKey words1问题描述一块大平板0xL，初始温度为常数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-23近似积分法解根据问题可知，该大平板传热可分为二个阶段：一、左侧x=0的边界接受热源的加热，其热量从左往右依次渗透。但当L，即热量只传递到大平板的某一处，假设，而并没有传递至最右侧x=L处的边界时，由于右侧边界的初始温度与环境相同，故其并不对环境放热；二、当>L，即左侧热量已传递至右侧时，大平板的左侧仍接受热源的加热，而其右侧那么对环境放热，直到放热量等于吸热量为止，整个大平板的温度趋于稳定。3.1L，即尚未整体渗透阶段假设此阶段大平板内温度分布为： 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-93.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-204有限差分法解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代入导热微分方程可得： 4-11令，经整理可得： 4-12根据能量守恒定律，对左、右边界条件进行隐式差分，可得： 4-13整理后可得： 4-14因此，式2-2所列的数学模型可离散为下式矩阵形式：4-15上式三角矩阵可通过LU分解，分解为一个上三角矩阵和一个下三角矩阵，再通过追赶法求解出各未知量。5大平板非稳态传热的MATLAB解通过上述分析，已经对大平板非稳态传热解析解以及数值解的求解过程有了一定的了解，此处就通过实例来比拟解析解与数值解之间的差异。5.1参数定义表5.1 大平板参数值(材料选用铝)工程密度比热容长度导热系数环境温度热流密度对流换热系数符号cLktfqw数值270088012203001000100单位K5.2解析解的MATLAB实现参考文献1.李格升, 郭庆平与刘春晓, 大平板第一类边界条件下非稳态传热近似解析解. 武汉理工大学学报交通科学与工程版, 2002: 第432-434页.2.李灿, 高彦栋与黄素逸, 热传导问题的MATLAB数值计算. 华中科技大学学报(自然科学版), 2002(9): 第91-93页.3.王飞与裴永祥, 有限差分方法的MATLAB编程. 新疆师范大学学报(自然科学版), 2003(4): 第22-27页.4.田兵, 用MATLAB解偏微分方程. 阴山学刊(自然科学版), 2006(4): 第12-13页.5.张春路等, 大平板瞬态热传导问题的一种新的近似解法. 工程热物理学报, 1999(1): 第82-85页.6.赵德奎与刘勇, MATLAB在有限差分法数值计算中的应用. 四川理工学院学报(自然科学版), 2005(4): 第61-64页.7.曹钢, 王桂珍与任晓荣, 一维热传导方程的根本解. 山东轻工业学院学报(自然科学版), 2005(4): 第77-80页.
