导热问题数值解法基础.pptVIP

二、内部节点的隐式差分格式常物性、无内热源、一维、非稳态导热内节点i、在k??时刻热平衡关系式：温度对时间的一阶导数改为向后差分第30页，共50页，星期日，2025年，2月5日内节点i、k??时刻：第31页，共50页，星期日，2025年，2月5日改写为：隐差分格式：(k+1)??时刻各节点的温度不能直接利用前一个时刻（k??）各节点的温度以显函数的形式表示（含有未知的）第32页，共50页，星期日，2025年，2月5日隐差分格式无条件稳定需将所有的内节点和边界节点方程式联立求解第33页，共50页，星期日，2025年，2月5日三、边界节点温度差分方程第一类边界条件：问题简单；边界节点温度给定第二类、第三类边界条件：根据给定的具体条件，针对边界节点所在的网格单元，写出热平衡关系式，建立边界节点温度差分方程显式或隐式差分格式假设边界条件为第三类：边界节点0，k??时刻热平衡：1、边界节点的显式差分格式第34页，共50页，星期日，2025年，2月5日第35页，共50页，星期日，2025年，2月5日第1页，共50页，星期日，2025年，2月5日§4-1建立离散方程的方法一、区域和时间的离散化—有限差分法的基本原理把物体分割为有限个离散的单元体，用有限差商代替导数，从而将微分方程转化为差分方程。通过数值计算求取各网格单元节点的温度二维导热问题；网格线；沿x、y方向的间距为?x、?y；网格单元节点：网格线的交点；p(i,j)p(i,j)边界节点（Nodalpoint）第2页，共50页，星期日，2025年，2月5日每个节点温度就代表了它所在网格单元的温度此方法求得的温度场在空间上不连续网格越细密、节点越多，结果越接近分析解p(i,j)对于非稳态导热：除了在空间上把物体分割成网格单元，时间上也要分割成许多间隔??网格越细密，计算所花时间越长——在时间上也是不连续的建立有限差分离散方程的常用方法：(1)Taylor（泰勒）级数展开法；(2)热平衡法；(3)多项式拟合法；(4)控制容积积分法（Energybalancemethod）第3页，共50页，星期日，2025年，2月5日p(i,j)Taylor（泰勒）级数展开法和多项式拟合法偏重于从数学的角度进行推导；而热平衡法和控制容积积分法则着重于从物理的观点来分析建立差分方程的方法(1)Taylor（泰勒）级数展开法★(2)多项式拟合法(3)热平衡法★(4)控制容积积分法前一类方法优点：便于对离散方程进行数学特性分析缺点：变步长网格的离散方程形式复杂、导出过程的物理概念不清晰、不能保证差分方程具有守恒特性后一类方法优点：推导过程物理概念清晰、离散方程系数具有一定物理意义、保证差分方程具有守恒特性缺点：不便于对离散方程进行数学特性分析第4页，共50页，星期日，2025年，2月5日1、Taylor（泰勒）级数展开法将函数t(x,y,z)在某点（i+1,j）对点（i,j）作Taylor展开二、建立差分方程的方法二阶导数和更高阶导数项之和——截断误差第5页，共50页，星期日，2025年，2月5日二阶导数和更高阶导数项之和——截断误差节点（i,j）一阶导数的向前差分表达式；一阶截差公式第6页，共50页，星期日，2025年，2月5日将函数t(x,y,z)在均匀网格点（i-1,j）对点（i,j）作Taylor展开：节点（i,j）一阶导数的向后差分表达式；一阶截差公式第7页，共50页，星期日，2025年，2月5日节点（i,j）一阶导数的中心差分表达式；二阶截差公式二方程相减，得：节点（i,j）一阶导数的三种表达式中，中心差分表达式的截断误差最小；尽可能采用中心差分表达式第8页，共50页，星期日，2025年，2月5日节点（i,j）二阶导数的中心差分表达式；二阶截差公式二方程相加，得：在表示温度对时间的一阶导数时一般采用向前或向后差分表达式——温度对时间的中心差分表达式不稳定第9页，共50页，星期日，2025年，2月5日节点P(i,j)的温度离散方程在写出导数差分表达式后，可以很容易建立导热微分方程的离散方程。以常物性、无内热源、二维稳态导热为例：第10页，共50页，星期日，2025年，2月5日2、热平衡法（Energybalancemethod）(i,j)oyx(i-1,j)W(i+1,j)E(i,j-1)S(i,j+1)NP?x?x?y?y热平衡法的