传热学数值解上机作业.docxVIP

能动A43丁济同2140301156问题提出所谓求解导热问题，就是对导热微分方程在规定的定解条件下的积分求解，此解称为分析解。但分析解只能求解集合形状和边界条件较简单的情况，对于复杂几何条件的导热问题，再加上计算机技术的迅猛发展，数值解的优点彰显了出来。本实验运用数值运算的方法求解一个导热问题。解决思路本实验主要依据了《传热学》第四章第163页的思路：建立控制方程及定解条件：导热方程，定解条件分别为等温边界条件和对流边界条件。节点：正方形节点。对于等温边界条件，做两组实验，第二组实验节点间距离是第一组节点间距离的一半，细化节点。代数方程的建立：运用能量守恒定律和傅立叶导热定律。求解方程组：迭代，即把前一组的运算值代入方程，得出下一组运算值，直到前后两组变量差值小于给定值（0.1）.解的分析：对节点温度值进行格式化输出，直观看出温度场分布。程序运行运行界面导热边界条件程序program djt1996implicit nonereal:: a(12,16)!格点温度数组real:: b(12,16)!迭代差值数组integer i,j,kreal tempreal q1,q2k=0!迭代次数q1=0q2=0a(7:12,7:16)=0!数组多余元素b(7:12,7:16)=0b(1,:)=0!等温边界格点，所以此位置温度恒为定值，不需要计算迭代差值b(:,1)=0b(6,6:16)=