安全工程数值计算原理及应用.docVIP

评 分 中国矿业大学（北京） 研究生课程考试试卷 考试科目 安全工程数值计算原理与方法 考试时间 学　　号 TSP140102073 姓 名 王静伟 所属学院 资源与安全工程学院 类别（硕士、博士、进修生） 硕士 评 语： 任课教师签名： 一维非稳态导热模型及数值解算 问题的提出 将一块面积较大的高温金属薄板，放入恒温的自然流动的空气室内，高温块体会从内向外传递热量，逐渐冷却，那么金属薄板的表面温度随时间到底是怎么变化的呢？ 问题的简化 由于金属薄板的长（l）和宽（m）远远大于厚度（d），故将金属薄板可以简化为无限大薄壁导热问题，即热量的传递沿厚度方向，为单项导热； 室内空间较大，可将室内气体视为恒温条件； 平板的厚度为d，密度为ρ，初始温度为To，导热系数为λ，比热为c，空气对流换热系数为h，室内温度为Ta； 平板内部导热满足傅里叶定律，表面换热满足牛顿冷却公式。 数理方程描述 由于平板换热关于中心线是对称的，仅对平板一半区域进行计算即可，坐标x的原点选在平板中心线上，因而一般区域的非稳态导热的数学描述为： (1)数理方程 边界条件 数值离散建模 1）求解域的离散 空间步长为(x，时间步长为 ((， (x、((大小的选择需要保证节点温度方程求解的稳定性，如图1。 2）节点温度差分方程的建立 运用热平衡法可以建立非稳态导热物体内部节点和边界节点温度差分方程。表示空间节点i在k((时刻（简称k时刻）的节点温度。 a、内部节点温度差分方程： b、外边界节点温度差分方程： c、中心线边界节点温度差分方程： 程序编写 public class { static int i,k=0; static int N=20; static double h=5;//空气自燃对流换热系数 static double λ=109;//黄铜的导热系数 static double Δx=0.01; static double c=3.9*100;//铜的比热 static double ρ=8900;//铜的密度 static double τ=1; static double a=λ/ρ/c; static double Fo=a*τ/Δx/Δx; static double Bi=h*Δx/λ; static double To=100; static double T1=20; static double []T=new double[N+1]; public static void main(String[] args) { for(i=0;i=N;i++){ T[i]=To; System.out.printf(%7.1f,T[i]); } System.out.println(); do{ k++; for(i=1;iN;i++){ T[i]=Fo*(T[i-1]+T[i+1])+(1-2*Fo)*T[i]; } i=0; T[i]=2*Fo*(T[i+1]-T[i])+T[i]; i=N; T[i]=2*Fo*(T[i-1]+Bi*T1)+(1-2*Bi*Fo-2*Fo)*T[i]; if(k%7200==0){ for(i=0;i=N;i++){ System.out.printf(%7.1f,T[i]); } System.out.println(); } }while(T[N]-T1=0.1); System.out.println(k); } } 结果 图. 黄铜金属薄板表面温度随时间的变化曲线