一维非稳态导热问题的数值解.docxVIP

计算传热学程序报告 题目：一维非稳态导热问题的数值解 学号： 学院：能源与动力工程学院 专业：工程热物理 日期：2014年5月25 R i维非稳态导热问题数值解 a2 a2T =0 (0xL) a dt T(xfi) = 0 7W) = l，7W) = 0 L = \.a = 1 1?方程离散化 对方程进行控制体积分得到: w dt =丄「（厂+4_厂皿 a Jw 非稳态项：选取T随x阶梯式变化，有 ￡ （厂1 - V ）dx = （7；严-T； ）Av 扩散项：选取一阶导数随时间做显示变化，有 )松哙);—( )松哙);—( dT lx ）；」 进一步取T随x呈分段线性变化，有 竺 _Tb：-TP oT _tp-tw dx e ~ (3x)e ，莎 ” (dx)w 整理可以得到总的离散方程为: T；严-T； _军一2耳+邛 △/ 2?计算空间和时间步长 取空间步长为： 网格Fourier数为： h=L/N 汽唱（小于。?5时稳定） 时间步长为： n = F0 — a 建立温度矩阵与边界条件 T=ones(N+1 ,M+1) T(:,l)=Ti (初始条件温度都为0) T(l,:)=To (边界条件x=0处温度为1) T(N+l,:)=Te (边界条件x=L处温度为0) 差分法求解温度 山离散方程可得到： 俨=佗s - 27；+兀)-7； 转化为相应的温度矩阵形式： T{m,k + V)= F() *[T(m+\ik)+T(m — lik)-2*T(m9k)] + T(mik) 5?输入界面 考虑到方程的变量，采用inputdig函数设置5个输入变量，对这5个变量设置了 默认值，如图1所示。在计算中可以改变不同的数值，得到不同的结果，特别注 意稳定条件的临界值是0.5。根据设置的默认值，得到的讣算结果如图2所示。 图1 mntlab变量输入界面  图2默认值的计算结果 6?结果分析 根据上面的分析，给出了程序的输入界面，以及默认值状态下的数值解。可以通 过改变不同的输入值，得到需要的分析结果，总结出了下面4点结论： （1）取F0=0.48,得到一维非稳态导热结果如下图所示 图2 Fo=O?48时一维非稳态导热 从图中可以看岀，对于长度L=1的细机 初始时刻=0时温度为0,边界条 件x=0时，T=l,边界条件x=l时，T=0o随着时间的增加，温度从x=0通过导热 的形式传递到x=l,不同时刻不同位置杆的温度都不同，并且随着时间的增加， 杆的温度也逐渐增加。 （2）取Fo=O.48,可以得到不同位置的温度响应曲线，如下图所示 图3 Fo=O.48时不同x位垃处的温度响应 图中红色曲线代表x=0」位置的温度瞬态响应，黑色曲线代表x=0.2位置的 温度瞬态响应，蓝色曲线代表x=0.4位置的温度瞬态响应。从图中可以看出，随 着x的增加，曲线与x轴的交点值越大，温度开始传递到该位置的所需的时间越 长。随着x的增加，温度响应曲线的变化速率越慢，最终的达到的温度也越低。 （3）取F（）=0.25,得到不同位置的温度响应曲线如下图所示 图4 Fo=0.25时不同x位苣处的温度响应 图中三条曲线分别是x=0.1, x=0.2, x=0.4位置的温度瞬态响应。与图3的 Fo=O.48进行对比，两种情况下的F（）值不同，F（）值越大表明热扩散系数d的值越 大。从图中可以看出热扩散系数对于导热的影响，Fo=0.25时，与Fo=0.48相比较，各位置开始响应时所需的时间较长，而且各位置响应曲线的变化速率较小, 最终的达到的温度也较低，说明了热扩散系数越小，热传导越慢，传递效率越低。 （4）取Fo=O.5H得到非稳定的数值解如图所示 0 100 200 300 400 500 600 700 800 900 1000 时间t 图6 Fo=O.51时不同x位置处的温度响应 从图中可以看出，对于显示格式的离散方程，并不是所有的F（）值都能得到 有意义的解，必须要求F00.5时才能得到稳定的数值解，当F（）0.5时，会出现 物理上不真实的解。 附件：(matlab程序) function heat_conduction() %—维齐次热传导方程 %设置输入界面 options={空间杆长LT空间点数N：时间点数MT扩散系数叫稳定条件的值Fo(临界值0.5)：}; topic』一维非稳态导热；％标题栏显示 lines=l;%输入行为1行 dcf= {r； 100710007 r/0.481}；% 默认值输入 f=inputdlg(options,topicjines,def);% 输入框设置 L=eval(f{l});%设置输入值 N=eval(f{2}); M=eval(f{3}); a=eval(f{4}); Fo=eval(f{5});%Fo的值必须小