02第二章导热型问题的求解.pptVIP

第二章 导热型问题的求解 第一节 空间区域离散方法 第二节 一维导热方程的离散－控制容积积分法 第三节 离散化的四项基本规则和界面物性参数的处理 第四节 边界条件的处理－附加热源法 第五节 TDMA算法 第六节 二维问题的离散化及求解 第七节 导热型问题的计算程序 第一节 空间区域离散方法 一、区域离散化 二、区域离散的四要素 三、区域离散方法的分类 四、区域离散的常用符号 一、区域离散化 在对传热问题进行数值计算时，将所计算的区域划分为许多互不重叠的子区域，确定节点在子区域中的位置及其所代表的容积（控制容积），这一过程称为区域离散化。 二、区域离散的四要素 1.节点：需要求解的未知物理量的几何位置； 2.控制容积：应用控制方程或守恒定律的最小几何单位； 3.界面：规定了与各节点相对应的控制容积的分界面位置； 4.网格线：沿坐标轴方向联结相邻两节点而形成的曲线簇。 三、区域离散方法的分类 1.外节点法（方法A）：先确定节点坐标再计算相应的界面。节点位于子区域的角顶上，划分子区域的曲线簇就是网格线，但子区域不是控制容积。为了确定各节点的控制容积，需在相邻两节点的中间位置作界面线，由这些界面线构成各节点的控制容积。 三、区域离散方法的分类 2.内节点法（方法B）：先规定界面位置而后确定节点的方法。节点位于子区域的中心，这时子区域就是控制容积，划分子区域的曲线簇就是控制体的界面线。 三、区域离散方法的分类 三、区域离散方法的分类 3.两种区域离散方法的区别： 三、区域离散方法的分类 3.两种区域离散方法的区别： 三、区域离散方法的分类 3.两种区域离散方法的区别： 三、区域离散方法的分类 3.两种区域离散方法的区别： (4)采用内节点法可以避免同一控制容积内物性发生突变。 四、区域离散的常用符号 当要对离散方程进行特性分析时，采用i－j－n表示法，即所研究的二维问题的节点位置记为(i, j)，与该节点相邻的界面则分别为i+1/2，i-1/2，j+1/2及j-1/2，n表示非稳态问题的时层。 四、区域离散的常用符号 采用P、N、E、W、S表示所研究的节点及相邻的四个节点，用n、e、w、s表示相应的界面，而用上标0表示非稳态问题中上一时层之值。相邻两节点间的距离，以x方向为例，以 表示，而 则表示相邻两界面间的距离。 四、区域离散的常用符号 第二节 一维导热方程的离散－控制容积积分法 一维导热方程的基本型式 第二节 一维导热方程的离散－控制容积积分法 一、最简化的型式 二、源项的处理 三、不稳态的处理 四、一维导热方程的离散方程 一、最简化的型式 一维稳态无热源 在P点离散有： 即 可整理为： 其中： 二、源项的处理 考虑源项后的型式 一般取线性处理 则方程可离散为： 其中： 三、不稳态的处理 考虑不稳态后的型式： 对时间和空间二重积分有： 即： 三、不稳态的处理 采用 于是离散方程为： 其中： 三、不稳态的处理 采用全显式 离散方程为： 要求 的系数大于0，故要求 采用全隐式 离散方程为： 属于无条件稳定，常被采用。 四、一维导热方程的离散方程 由前面三方面的分析可得到一维导热方程 的离散方程： 其中： 第三节 离散化的四项基本规则和界面物性参数的处理 一、离散化的四项基本规则 二、界面物性参数的处理 一、离散化的四项基本规则 1.控制容积面上的相容性 2.正的系数 3.源项负斜率线性化 4.邻近系数之和 二、界面物性参数的处理 1.采用线性处理： 2.采用调和平均法 第四节 边界条件的处理－附加热源法 一、第二类边界条件的处理 二、第三类边界条件的处理 一、第二类边界条件的处理 离散方程是： 由于 一、第二类边界条件的处理 代入离散方程有： 令 有 与 相比： （1） （2) 源项 相当于增加了源项 （3)校核边界温度 二、第三类边界条件的处理 同第二类边界条件的分析类似，有以下三点： （1） （2) 源项 相当于增加了源项 （3)校核边界温度 第五节 TDMA算法 对 考虑边界条件 i=1时， ， 可用 表示； i=2时，