计算传热学第3讲数学模型与求解区域的离散化.pptVIP

计算传热学第3讲 数学模型与求解区域的离散化 Discretization of Mathematical Models and Computation Domain 作业与阅读要求 阅读：陶文铨《数值传热学》第2章 作业：P44 题2－3，题2－4 作业：P46 题2－11 本讲主要内容 求解区域的离散化 Taylor级数展开法 控制方程的离散化－Taylor级数法 控制方程的离散化－控制容积法 控制方程的离散化－变物性的情况 控制容积法 Taylor级数法 交界面参数的计算 四个基本原则 源项的线性化 三个关键环节 建立恰当的数学模型 Proper Mathematical Modelling 对求解区域进行离散化处理 Discretization of Computational Domain 对数学模型进行离散化处理 Discretization of Mathematical Model 3-1求解区域的离散化 3.1.1求解区域的界定： 有限区域（finite domain）： 求解区域（Computational domain）＝实际区域 无限区域（infinite domain）： 求解区域?实际区域 界定原则：计算结果不敏感原则，亦即，求解区域的大小对计算结果没有明显的影响。 例子： 求解区域的界定：例子 流动问题的出口界面： 求解区域的界定：例子 无穷大区域的“无穷远界面” 半无限大介质中的稳态导热 求解区域的界定：例子 无穷大区域的“无穷远界面” 无限大介质中的非稳态导热 求解区域的界定 对称区域：对称问题的求解区域 3.1.2 求解区域的离散化 什么是求解区域的离散化 将求解区域划分为若干个互不重合的子区域（CV) 不重合 子区域（sub-region) 控制容积（control volume) 确定节点在子区域中的位置 给出节点位置坐标 节点所代表的区域及其大小 方法： 用一组正交的网格线（可以是曲线）将求解区域进行分割 求解区域的离散化：方法一 外节点法或节点-控制容积法 网格线的交点作为节点 节点所代表的求解区域（控制容积） 由两节点间中心位置的对称界面围成的区域。 例子：二维矩形区域 求解区域的离散化 节（结）点：网格线的交点 控制容积（节点所代表的求解区域）：两节点中间界面所围成的区域。 节点的分类： 相邻接点：坐标轴方向上相差一个步长的节点 内部节点：所有相邻节点都属于求解区域的节点 边界节点：至少有一个相邻节点不属于求解区域 节点的命名 研究对象点：P（i,j) 相邻节点：按方位关系或位置坐标 求解区域的离散化 确定区域离散化的要素 节点位置坐标 控制界面位置 节点间距 控制容积的大小 求解区域的离散化：方法二 内节点法或控制容积-节点法 先划定控制容积（节点所代表的求解区域） 节点：控制容积的几何中心 例子：二维矩形区域 3.1.3求解区域的离散化： 方法比较 边界节点所代表的求解区域（控制容积）不同： 外节点法：半个控制容积 内节点法：容积为0的控制容积 节点在控制容积中的位置不同 外节点法： 控制界面始终位于两节点中间位置上：导数计算准确 不能保证节点始终位于CV的几何中心上 内节点法： 节点始终位于CV的几何中心上：非稳态项计算准确 不能保证控制界面始终位于两节点中间位置上 求解区域的离散化：方法比较 当网格划分足够细时，两者没有本质区别 内节点法： 边界节点处理较简单 边界相邻节点：要特别注意处理方法，与其它内部节点有所不同 历史及习惯的原因：内节点应用较广泛 内节点法在边界相邻节点处始终是非均匀网格 可能会产生较大的误差 求解区域的离散化：网格参数 一维为例 网格参数：名称与定义 (?x)w＝(?x)+w＋(?x)-w 节点W－P之间的距离 (?x)e＝(?x)+e＋(?x)-e 节点P－E之间的距离 (?x)+w 控制界面w－节点P之间的距离 (?x)-e 节点P－控制界面e之间的距离 ?x ＝ (?x)+w ＋(?x)-e 控制容积 w ， e 左、右控制面 网格参数：各参数之间的关系 外节点法 (?x)+w＝? (?x)w ； (?x)-e ＝? (?x)e ?x ＝? [(?x)w ＋(?x)e ] 内节点法 (?x)+w＝ (?x)-e ＝? ?x 3.2 Taylor级数展开法 控制方程离散化的方法： Taylor级数法 多项式拟合法 控制容积法 。。。。 Taylor级数法和控制容积法最为重要 Taylor级数法的基本思路 借助Taylor级数展开给出各阶导数的差商表达式 将方程中的各阶导数用相应的差商表达式代替 整理化简 3.2.1 Taylor级数展开法－等步