《热传导与~热辐射大作业报告》由会员分享，可在线阅读，更多相关《热传导与~热辐射大作业报告（28页珍藏版）》请在金锄头文库上搜索。
1、. 热传导与热辐射大作业报告 .目录一、作业题目- 1 -二、作业解答- 2 -个人感想- 17 -附件.计算中所用程序- 18 - 一、作业题目一矩形平板, ，内有均匀恒定热源，在及处绝热，在及处保持温度，初始时刻温度为，如右图1所示：1、求时，矩形区域内的温度分布的解析表达式；2、若，热传导系数，热扩散系数。请根据1中所求温度分布用MATLAB软件绘出下列结果，加以详细物理比较和分析：(a) 300s内，在同一图中画出点、（单位：m）温度随时间的变化；(b) 200s内，画出点、（单位：m）处，分别沿x、y方向热流密度值随时间的变化；(c) 画出时刻区域内的等温线；(d) 300s内，在同
2、一图中画出点（单位：m）在分别等于，情况下的温度变化；(e) 300s内，比较点(9,6) （单位：m）在其它参数不变情况下热导率分别为、和的温度、热流密度变化；(f) 300s内，比较点(9,6) （单位：m）在其它参数不变情况下热扩散系数分别为、和的温度、热流密度变化；3、运用有限差分法计算2中(b)、(d)和(e)，并与解析解结果进行比较，且需将数值解与解析解的相对误差减小到1以下；4、附上源程序和个人体会；以报告形式整理上述结果，用A4纸打印上交。二、作业解答1、求时，矩形区域内的温度分布的解析表达式；解答：我们令，则可以得到一个方程和边界条件： （1-1） 将上式分解为一个的稳态问题
3、： （1-2） 和一个的其次问题： （1-3） 其中则原问题的解根据下式求得： （1-4）发热强度为常数的特解可从表2-4中查的，则新变量可定义为： （1-5）将（1-5）带入（1-2）整理得到： （1-6） 若令常数，则上式可以变为： （1-7） 其中假定可以分离出如下形式： （1-8）对应于的分离方程为： （1-9） （1-10） 在中特征值问题的解可以直接从表2-2第6条中得到，只需要用a代替L， （1-11） （1-12）是下面方程的正根： （1-13）方程（1-10）的解可以取为 （1-14）的完全解由下式组成： （1-15）此式满足热传导问题（1-7）及三个齐次边界条件，其中，系数
4、可以根据方程的解还应满足非齐次的边界条件来决定。利用的边界条件可得： （1-16）利用函数的正交性可以求得系数， （1-17）式中：将这个表达式带入式（1-15），其中范数在前面已经给出，解得结果为 （1-18）则： （1-19）假定分离成如下表达式 （1-20）对应于函数和的分离方程为： （1-21） ： （1-22） 的解为： （1-23）上述问题的完全解为： （1-24）其中0xa,0yb。当t=0时，上式变为： （1-25）其中0xa,0yb。确定未知系数的方法是，在上式两边逐项用如下算子作运算： 及并利用这些函数的正交性，得到： （1-26）最终得到问题的解为： （1-27）式中出现
5、的特征函数，特征值及范数可以从表2-2中直接查得： （1-28） （1-29）且为如下方程的正根： （1-30）满足特征值问题的函数对应于表2-2中的第6条，得到： （1-31） （1-32）且是如下方程的正根： 最后得到： （1-33）令 其中令令 根据余弦函数的正交性，只有当m=n时积分才不为0，故上式可以化为：再令 所以令所以由（1-4）、（1-19）及（1-33）可知 。以上是解析解的全过程，具体值的计算采用MATLAB编程计算求取。2、若，热传导系数，热扩散系数。请根据1中所求温度分布用MATLAB软件绘出下列结果，加以详细物理比较和分析：300s内，在同一图中画出点、（单位：m）温
6、度随时间的变化；图1.不同点温度随时间变化曲线图分析：开始时刻通过右、上边界向内部导热，这时候尽管有内热源，但谁相对离右、上边界越近，温度曲线越陡。即开始时刻（0,8）点比（0,4）点温度曲线陡，（12,0）点比（6,0）点温度曲线陡，一定时间后由于有内热源，内部温度逐渐高于边界温度，这时内部开始向边界导热。这时谁离两个绝热边交点越近，谁的温度会越高，这就是为什么最后（0,4）点比（0,8）点温度高，（6,0）点温度比（12,0）点温度高。（b）.200s内，画出点、（单位：m）处，分别沿x、y方向热流密度值随时间的变化； 图2.200s内x方向不同点的热流密度曲线（解析解） 图3.200s内
7、y方向不同点的热流密度曲线（解析解）图4.200s内x方向不同点的热流密度曲线（数值解） 图5.200s内y方向不同点的热流密度曲线（数值解） 图7.不同点x方向热流密度数值解与解析解相对误差 图6.不同点x方向热流密度数值解与解析解相对误差分析：右边界（18,4）和（18,8）这两点开始时X方向两侧温差较大，故热流密度也会特别大，开始时内部温度较边界温度低，向内部导热，热流密度为负值。后来内部温度大于边界温度，向外散热，热流密度为正值。而上边界点温度相同，故在X方向不存在热传导，故导热系数为零。而中间点开始时从右向左导热，热流密度为负，随着边界层温度影响的深入，热流密度绝对值越来越大，但由于
8、有内热源，会使此影响逐渐减弱，故在一段时间后待热流密度达到一个顶峰以后会逐渐减小，后来由于内热源的作用，导热由内向外进行，热流密度也由负值变为正值。Y方向分析类似。由于（9,6）离上边界更近，故沿Y方向达到的下边界峰值更大。（c）.画出时刻区域内的等温线； 图8.50s时区域内的等温线 图9.75s时区域内的等温线 图10.100s时区域内的等温线 图11.125s时区域内的等温线 图12.150s时区域内的等温线分析：开始时刻，尽管有内热源的存在，但边界温度比内部温度高，此时边界向内部传热，故开始时靠近边界的温度比内部高，这就是为什么50、75、100s时等温线呈现由坐下到右上温度逐渐升高。
9、过一段时间后，中间部分由于内热源和边界热传导的共同作用，而坐下边界此时收到的内热源和边界热传导的作用小于中间部分，故造成了中间部分温度反而比其他部分高。一段时间后，内热源起主导作用，向外散热，这事等温线上的温度由左下到右上逐渐降低。（d）.300s内，在同一图中画出点（单位：m）在分别等于，情况下的温度变化；图13.不同内热源下温度变化曲线（解析解） 图14.不同内热源下温度变化曲线（数值解）图15.不同内热源下数值解与解析解相对误差分析：内热源越大，单位时间内内部产生的能量越多，节点温度升高的越快。在其它条件相同的情况下，内热源越大，最终内部温度也越高。开始时，由于温度变化剧烈，此时解析解和
10、数值解的误差也相对较大，一段时间以后温度趋于稳定，这个时候相对误差也趋于一个较小的稳定值。（e）.300s内，比较点(9,6) （单位：m）在其它参数不变情况下热导率分别为、和的温度、热流密度变化； 图16.不同导热数下温度变化曲线（解析解） 图17.不同导热数下温度变化曲线（数值解）图18.不同导热系数下数值解和解析解的相对误差 图20.不同导热系数下X方向热流密度曲线（数值解）图19.不同导热系数下X方向热流密度曲线（解析解） 图20.不同导热系数下X方向热流密数值解与解析解相对误差图22。不同导热系数下Y方向热流密度曲线（数值解）图21.不同导热系数下Y方向热流密度曲线（解析解） 图23.(9,6)点不同导热系数下y方向数值解和解析解的相对误差分析：导热系数K越大，内部温度越能快速的传递给外界，这就是问什么导热系数越大，节点最终温度低。根据热流密度方程，可知。K越大，热流密度越大，这就是为什么K越大，热流密度最低点峰值越大。而最后由于内热源相同，根据能量守恒，最后导热系数也必然趋近于一个定值。开始时由于温度变化剧烈，在不同的导热系数下同一点温度随之间变化值得数值解和解析解的相对误差较大，一段时间后温度趋于稳定，此时数值解和解析解的相对温差是一个较小
