【2017年整理】清华大学传热大作业_无限大平壁简单数值计算报告.docxVIP

传热大作业报告 2011010*** 热动** *** 大作业题目 一厚度为0.1m的无限大平壁，两侧均为对流换热边界条件，初始时 两侧流体温度与壁内温度一致，tf1=tf2=t0=5 ℃；已知两侧对流换热系数分别为h1=11 W/m2K、h2=23W/m2K, 平壁材料的导热系数?=0.43W/mK，导温系数a=0.3437×10-6 m2/s。如果一侧的环境温度tf1突然升高为50℃并维持不变，计算在其它参数不变的条件下，平壁内温度分布及两侧壁面热流密度随时间的变化规律（用图形表示）。 要求：将全部计算内容（包括网格的划分、节点方程组、计算框图、程序及计算结果）用A4纸打印。 网格划分 如图，将无限大平板作为一维处理，本题为一维非稳态导热问题，对流换热边界条件。 空间网格划分：平板总厚度为delta=0.1m，定义空间步长为 dx=0.005m，则距离份数为N=delta/dx=20份。定义x{n}为以0为首项，以dx为公差的等差数列，尾项为delta=0.1m，共有N+1项，则x{n}中的每一项即表示一个沿平板厚度方向中的划分点。 时间网格划分：设总时间长度为T= 100000s，定义时间步长为dtao=20s，则时间份数为M=5000份。定义 tao{m}是以0为首项，以dtao为公差的等差数列，尾项为T=100000s，共有M+1项,则tao{m}中每一项即表示一个时刻。 计算框图 程序中的各个变量的名称及意义： 题设中各个常数 lambda=0.43 导热系数；a=0.3437e-6 热扩散率；h1=11 边界对流换热系数；h2=23 边界对流换热系数2；t0=5 初始温度；tf1=50 初始流体温度；tf2=5 初始流体温度2；delta=0.1 总距离长度（无限大平板厚度）； 网格划分所设的变量 T=100000 总时间长度（在T时间内考虑本问题）；dtao=20 定义时间步长；dx=0.005 定义距离步长；M=floor(T/dtao) 时间份数=总时间/时间步长（向下取整）；N=floor(delta/dx 距离份数=总厚度/距离步长(向下取整）；tao=0:dtao:T 定义时间划分单元（以0为首项，以dtao为公差的等差数列，尾项为T），共有M+1项；x=0:dx:delta 定义距离划分单元（以0为首项，以dx为公差的等差数列，尾项为delta），共有N+1项； 判定稳定性的准则数 Bi1=h1*dx/lambda 边界节点网格毕渥数；Bi2=h2*dx/lambda 边界节点网格毕渥数2；Fo=a*dtao/dx^2 傅里叶数； 程序计算框图 开始 开始 输入 输入delta,T,dtao,dx,M,N,tao,x Fo1/(2*Bi1+2)Fo1/(2*Bi2+2) Fo1/(2*Bi1+2)Fo1/(2*Bi2+2) No Yes 建立 建立t(M+1,N+1)温度矩阵，令t(1,:)=t0,令m=1 m=m+ m=m+1 打印“不稳定” 打印“不稳定” t(m,1)=2*Fo*(t(m-1,2)+Bi1*tf1)+(1-2*Bi1*Fo-2*Fo)*t(m-1,1) t(m,1)=2*Fo*(t(m-1,2)+Bi1*tf1)+(1-2*Bi1*Fo-2*Fo)*t(m-1,1) t(m,N+1)=2*Fo*(t(m-1,N)+Bi2*tf2)+(1-2*Bi2*Fo-2*Fo)*t(m-1,N+1) t(m,n)=Fo*(t(m-1,n-1)+t(m-1,n+1))+(1-2*Fo)*t(m-1,n) mM mM+1 No Yes 输出温度矩阵 输出温度矩阵t(M+1,N+1)和相应图象 停机 停机 程序代码 本程序在MATLAB R2008a中运行通过，以下是源代码（%后为注释）： lambda=0.43;%导热系数 a=0.3437e-6;%热扩散率 h1=11;%边界对流换热系数 h2=23;%边界对流换热系数2 t0=5;%初始温度 tf1=50;%初始流体温度 tf2=5;%初始流体温度2 delta=0.1;%总距离长度（无限大平板厚度） T=100000;%总时间长度（在T时间内考虑本问题） dtao=20;%定义时间步长 dx=0.005;%定义距离步长 M=floor(T/dtao);%时间份数=总时间/时间步长（向下取整） N=floor(delta/dx);%距离份数=总厚度/距离步长(向下取整） tao=0:dtao:T;%定义时间划分单元（以0为首项，以dtao为公差的等差数列，尾项为T），共有M+1项 x=0:dx:delta;%定义距离划分单元（以0为首项，以dx为公差的等差数列，尾项为delta），共有N+1项 Bi1=h1*dx/lambda;%边界