一维非稳态导热问题的数值解答辩.doc

计算传热学程序报告 题目：一维非稳态导热问题的数值解 姓名： 学号： 学院：能源与动力工程学院 专业：工程热物理 日期：2014年5月25日一维非稳态导热问题数值解 求解下列热传导问题： 1.方程离散化 对方程进行控制体积分得到： 非稳态项：选取T随x阶梯式变化，有 扩散项：选取一阶导数随时间做显示变化，有 进一步取T随x呈分段线性变化，有 ， 整理可以得到总的离散方程为： 2.计算空间和时间步长 取空间步长为： h=L/N 网格Fourier数为： （小于0.5时稳定） 时间步长为： 建立温度矩阵与边界条件