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

计算传热学程序报告 题目：一维非稳态导热问题的数值解 姓名： 学号： 学院：能源与动力工程学院 专业：工程热物理 日期：2014 年 5 月 25 日 一维非稳态导热问题数值解 求解下列热传导问题： 2 T 2 x 1 T t T ( 0 x,0) (0 0 x L) T (0, t) 1,T(L,t) 0 L 1, 1 1.方程离散化 对方程进行控制体积分得到： t t t e w 2 T 1 t t d x d t 2 x t e w T t d x d t t t t [( T x ) ( e T x ) w ] dt 1 e w (T t t t T )dx 非稳态项：选取 T 随 x 阶梯式变化，有 e w (T t ) ( ) t t t t t T dx T T p p x 扩散项：选取一阶导数随时间做显示变化，有 t t t [( T T T T t t ) ( ) ]dt [( ) ( ) e w e w x x x x ] t 进一步取 T 随 x 呈分段线性变化，有 T T T E P ( ) ， e x x ( ) e ( T x ) T T P W w x ( ) w 整理可以得到总的离散方程为： 1 t T E t t T P t T E t 2T P t T W t 2 x 2.计算空间和时间步长 取空间步长为： h=L/N 网格 Fourier 数为： t t F （小于 0.5 时稳定） 0 2 2 x x 时间步长为： n F 0 2 h 3.建立温度矩阵与边界条件 T=ones(N+1,M+1) T(:,1)=Ti (初始条件温度都为 0） T(1,:)=To (边界条件 x=0 处温度为 1） T(N+1,:)=Te (边界条件 x=L 处温度为 0） 4.差分法求解温度 由离散方程可得到： t t t t t t TE F0(T 2T T ) T E P W P 转化为相应的温度矩阵形式： T (m, k 1) F0 [T(m 1,k) T (m 1, k) 2 T(m,k)] T(m, k) 5.输入界面 考虑到方程的变量，采用 inputdlg 函数设置 5 个输入变量，对这 5 个变量设置了 默认值，如图 1 所示。在计算中可以改变不同的数值，得到不同的结果，特别注 意稳定条件的临界值是 0.5。根据设置的默认值，得到的计算结果如图 2 所示。 图 1 matlab 变量输入界面 图 2 默认值的计算结果 6.结果分析 根据上面的分析， 给出了程序的输入界面， 以及默认值状态下的数值解。 可以通 过改变不同的输入值，得到需要的分析结果，总结出了下面 4 点结论： （1）取 F0=0.48，得到一维非稳态导热结果如下图所示 图 2 F0=0.48 时一维非稳态导热 从图中可以看出，对于长度 L=1 的细杆，初始时刻 t=0 时温度为 0，边界条 件 x=0 时，T=1,边界条件 x=1 时，T=0。随着时间的增加，温度从 x=0 通过导热 的形式传递到 x=1，不同时刻不同位置杆的温度都不同，并且随着时间的增加， 杆的温度也逐渐增加。 图 2 默认值的计算结果 6.结果分析 根据上面的分析， 给出了程序的输入界面， 以及默认值状态下的数值解。 可以通 过改变不同的输入值，得到需要的分析结果，总结出了下面 4 点结论： （1）取 F0=0.48，得到一维非稳态导热结果如下图所示 图 2 F0=0.48 时一维非稳态导热 从图中可以看出，对于长度 L=1 的细杆，初始时刻 t=0 时温度为 0，边界条 件 x=0 时，T=1,边界条件 x=1 时，T=0。随着时间的增加，温度从 x=0 通过导热 的形式传递到 x=1，不同时刻不同位置杆的温度都不同，并且随着时间的增加， 杆的温度也逐渐增加。 图 2 默认值的计算结果 6.结果分析 根据上面的分析， 给出了程序的输入界面， 以及默认值状态下的数值解。 可以通 过改变不同的输入值，得到需要的分析结果，总结出了下面 4 点结论： （1）取 F0=0.48，得到一维非稳态导热结果如下图所示 图 2 F0=0.48 时一维非稳态导热 从图中可以看出，对于长度 L=1 的细杆，初始时刻 t=0 时温度为 0，边界条 件 x=0 时，T=1,边界条件 x=1 时，T=0。随着时间的增加，温度从 x=0 通过导热 的形式传递到 x=1，不同时刻不同位置杆的温度都不同，并且随着时间的增加， 杆的温度也逐渐增加。 图 2 默认值的计算结果 6.结果分析 根据上面的分析， 给出了程序的输入界面， 以及默认值状态下的数值解。 可以通 过改变不同的输入值，得到需要的分析结果，总结出了下面 4 点结论： （1）取 F0=0.48，得到一维非稳态导热结果如下图所示