计算传热学第讲扩散方程的数值解.pptVIP

边界条件的处理说明：尽管它与一阶Taylor级数展开法的结果形式上相同，但它却是二阶精度的！请大家证明这一结论。采用内节点法划分网格时，即使在均匀网络的前提下，第1个近边界节点也不是等步长的。qBx=0(?x)2(?x)3123(?x)1(?x)2从图中可以清楚地看出这一点即使(?x)2=(?x)3(?x)1也不等于(?x)2所以要对第一个内部节点给予特别注意。第30页，共83页，星期日，2025年，2月5日边界条件的处理例如，对于直角坐标系，对于节点2，qBx=0(?x)2(?x)3123(?x)1(?x)2注意：CV2的左控制面w与节点W（节点1）重合，即与左边界重合！控制面e！第31页，共83页，星期日，2025年，2月5日边界条件的处理或写成，第32页，共83页，星期日，2025年，2月5日边界条件的处理附加源项法（Additionalsourcetermmethod)以内节点法为例由方程（19）解出边界节点上的待求变量T1，代入与第1个近边界节点的差分方程（21），第33页，共83页，星期日，2025年，2月5日边界条件的处理代入与第1个近边界节点的差分方程（21），整理后得到，第34页，共83页，星期日，2025年，2月5日边界条件的处理整理后得到，或者写成，其中，第35页，共83页，星期日，2025年，2月5日边界条件的处理其中Additionalsourceterm!第36页，共83页，星期日，2025年，2月5日边界条件的处理对于第3类边界条件，也可以做类似的处理，但是这时，请大家证明，第37页，共83页，星期日，2025年，2月5日边界条件的处理附加源项法的实质边界节点消去法不仅能用于内节点网格，也能用于外节点网格实施方法：计算附加源项：Sc,ad，Sp,ad把附加源项计入该控制容积中的源项中令与边界节点对应的系数（aW）等于0第38页，共83页，星期日，2025年，2月5日特别提示边界条件的处理是传热问题数值计算最重要的环节之一元体能量平衡法的基础地位尽可能采用外节点法划分网格边界节点消去法广泛应用提高收敛速度尽可能采用均匀网格边界节点的离散化方程在形式上与内部节点的相同第39页，共83页，星期日，2025年，2月5日4.1.7差分方程的求解将前面得到的差分方程（8）改写为，或者简单的写成矩阵的形式，其中，第40页，共83页，星期日，2025年，2月5日4.1.7差分方程的求解第41页，共83页，星期日，2025年，2月5日4.1.7差分方程的求解与方程（33）对比，知，由方程（33）系数阵[A]的特殊性，通常称之为三对角方程（Tri-diagonalequation)三对角方程可以采用非常高效的追赶法（TDMA法）求解基于矩阵分解属于必须掌握的内容第42页，共83页，星期日，2025年，2月5日4.1.7差分方程的求解TDMA法Fortran源程序第43页，共83页，星期日，2025年，2月5日4.1.8计算机实现：算例求解下面的一维稳态导热问题：第44页，共83页，星期日，2025年，2月5日4.1.8计算机实现：算例求解区域的离散化：内节点法：先划分控制容积，在确定节点均匀网格：?x=?x将整个求解区域划分为（N-2）个控制容积，N个节点（包括2个边界节点）内部节点的差分方程第45页，共83页，星期日，2025年，2月5日4.1.8计算机实现：算例其中，第46页，共83页，星期日，2025年，2月5日4.1.8计算机实现：算例注意：采用内节点法划分网格时，近边界节点与其它内部节点不尽相同，所以必须单独考虑。123NN-1(?x)w(?x)e(?x)e当i＝2时(?x)w=??x?w=?W=?1所以，当i＝2时第47页，共83页，星期日，2025年，2月5日4.1.8计算机实现：算例所以，当i＝2时，第48页，共83页，星期日，2025年，2月5日4.1.8计算机实现：算例当i＝3，4，。。。，N-2时，第49页，共83页，星期日，2025年，2月5日4.1.8计算机实现：算例同样，当i＝N-1时(?x)e=??x?e=?E=?N所以，当i＝N-1时123NN-1(?x)w(?x)e(?x)e第50页，共83页，星期日，2025年