数值传热_传热学上机实验_墙角稳态导热问题数值模拟解析.pdf

热学上机实验报告 墙角稳态导热问题数值模拟 核工 一、 问题描述 有 个用砖砌成的长方形截面的冷空气通道，其截面尺寸如图所示，假设在垂 直于纸面方向上冷空气及砖墙的温度变化很小，可以近似的予以忽略：在下列 两种情况下试计算：（1）砖墙横截面上的温度分布；（2 ）垂直于纸面方向的 每米长度上通过砖墙的导热量。 第 种情况：内外壁分别均匀的维持在0℃及30℃； 第二种情况：内外表面均为第三类边界条 ，且已知： 图一 二、 计算原理 本次上机模拟实验选等温边界条 。墙角是中心对称的，所以取其 1/4 研究， 方便计算机计算。上机模拟选取网格划分方法同实际实验，可根据热平衡法列 出节点方程，各方向导入单元体的热量之和为零。该边界条 下共有四类节点， 内节点、内边界点、外边界点和绝热边界点。 图二 四种节点的节点方程简化如下： 三、 编程求解 使用C++编程如下： #include iostream #include fstream #include iomanip using namespace std; int main() { ofstream out(shuju.txt); out setiosflags(ios::fixed); int i,j ; double temp,q_in,q_out,q; double eps= 1; double A[12][16]; //设置迭代初场 for(i=0;i 16;i++) A[11][i]=30; for(i=0;i 12;i++) A[i][0]=30; for(i=5;i 16;i++) for(j =0;j7;j++) A[j ][i]=0; for(j =5;j 16;j++) for(i=7;i 11;i++) A[i][j ]=6*(i-6); for(i=0;i 11;i++) for(j = 1;j5;j++) A[i][j ]=6*(5-j ); //建立迭代方程组并求解 while(eps 1.0E-4) { for(i=7;i 11;i++) { for(j = 1;j 15;j++) { eps= 1; temp=A[i][j ]; A[i][j ]=(A[i-1][j ]+A[i+ 1][j ]+A[i][j -1]+A[i][j + 1])/4; eps=A[i][j ]-temp; } eps= 1; temp=A[i][15]; A[i][15]=(A[i-1][15]+A[i+ 1][ 15]+2*A[i][14])/4; eps=A[i][ 15]-temp; } for(j = 1;j5;j++) { for(i= 1;i7;i++) eps= 1; temp=A[i][j ]; A[i][j ]=(A[i-1][j ]+A[i+ 1][j ]+A[i][j -1]+A[i][j + 1])/4; eps=A[i][j ]-temp; } eps= 1; temp=A[0][j ]; A[0][j ]=(A[0][j -1]+A[0][j + 1]+2*A[1][j ])/4; eps=A[i][ 15]-temp; } //计算墙体外表面导热量 q_out=0; for(i= 1;i 11;i++) q_out=q_out+A[i][0]-A[i][ 1]; for(j = 1;j 15;j++) q_out=q_out+A[11][j ]-A[10][j ]; q_out=q_out+(A[0][0]-A[0][1]+A[11][ 15]-A[10][15])/2; q_out=q_out*0.53; //计算墙体内表面导热量 q_in=0; for(i= 1;i7;i++) q_in=q_in+A[i][4]-A[i][5]; for(j =5;j 15;j++) q_in=q_in+A[7][j ]-A[6][j ]; q_in=q_in+(A[0][4]-A[0][5]+A[7][15]-A[6][15])/2; q_in=q_in*0.53; //计算平均导热量和相对误差 q=(q_in+q_out)/2; eps=abs(q_in-q_out); } //输出