计算传热学作业.doc

计算传热课程作业（一） 姓名： 李雪莹 学号： 110488 题目： 对于如下的具有源项的对流扩散问题，采用中心差分格式（控制面上的变量值取两侧节点值的算术平均值）迎风格式（视符号取上游值）求出该问题的数值解，并探讨数值解网格的依赖性。其中取，取。 要求： 先求出解析解； 均匀划分网格，写出离散化代数方程（包括内部节点和边界节点）； 求出不同参数组合（，）时的数值解并与精确解对比； 给出两种不同参数组合时，无网格依赖性的解。 一、解析解 解微分方程得 （1） 由 ； 得 ， 代入到（1）式，得 二、离散方程 在以P点为中心的控制容积上积分 得 即 （2） （1）采用中心差分格式 控制面上的变量值取两侧节点值的算术平均值 ， 将控制面上的导数用节点上的控制变量来表示，（2）式即为 （3） 因为网格采用均匀划分，则，（3）式变为 也即 离散化形式为 其中，，， （） （2）采用迎风格式 假定流动方向为w到e，则取上游值有 ， ，则（2）式可变为： 离散化形式仍为： 其中，，， （） 若假定流动方向为e到w,则-Pe，，， （） 三、不同参数组合时的数值解并与精确解对比 1.利用中心差分格式处理 当,时 当,时 当时 当,时 2.利用迎风格式离散化 其中u的方向假定为w到e 当,时 当, 其中u的方向假定为e到w 当，时 当， 四、给出两种不同参数组合时，无网格依赖性的解 中心差分格式 ，时 迎风格式 ，时 五、结论 从不同的离散格式可以看出，对于中心差分和迎风格式当较小时，数值解与精确解相比是一致的；较大时，其值与精确值之间有较大的差距。