导热问题的数值解法基础-第四章.ppt

第四章导热问题的数值解法上海交通大学第四章导热问题的数值解法本章要点理解导热问题数值解法的基本思路及其实施过程熟练掌握利用热平衡法建立内节点及边界节点离散方程的过程与方法（本章重点）理解非稳态导热问题的显式与隐式差分格式及其稳定性的判别了解数值解法的基本方法及求解代数方程组的常用方法4.1.1基本思想把原来在空间与时间坐标中连续的物理量的场（如温度场、速度场及浓度场等），用一系列有限个离散点上的值的集合来代替，通过一定的原则建立起这些离散点上变量之间关系的代数方程，并进行求解，从而获得离散点上被求物理量的值；并称之为数值解；4.1.2导热问题常用数值解法根据区域离散方式、代数方程离散方式及代数方程求解方法的不同，有：1、有限差分法(FDM):finite-differencemethod2、有限元法(FEM):finite-elementmethod3、有限分析法(FAM):finite-analysismethod4、有限容积法(FVM):finite-volumemethod5、边界元法(BEM):boundary-elementmethod1、物理模型：对实际传热问题进行简化处理，从而构成可用数学语言（公式）来描述的物理问题；2、数学模型：根据所得物理模型建立其控制方程，并确定初始与边界条件；3、区域离散化（网格划分）：将所研究的导热区域沿坐标方向（x、y或z）各自均匀或非等分，得到各个节点，每一个节点可看成是以它为中心的小区域的代表；4、节点方程离散化：根据傅立叶导热定律及热平衡法，建立节点物理量的代数方程（关键）；5、设置初场，求解代数方程组6、解的分析1、有限差分法的基本思想基本思想：把原来在时间和空间坐标中连续变化的物理量（如温度、压力、速度和热流等），用有限个离散点上的数值集合来近似表示数学基础：用差商代替微商（导数）几何意义：用函数在某区域内的平均变化率代替函数的真实变化率差分格式：向后差分格式、向前差分格式、中心差分格式2泰勒级数展开法根据泰勒级数展开式，用节点(i,j)的温度ti,j来表示节点(i+1,j)的温度ti+1,j用节点(i,j)的温度ti,j来表示节点(i-1,j)的温度ti-1,j若取上面式右边的前三项，并将式①和式③相加移项整理即得二阶导数的中心差分：同样可得：

对于二维稳态导热问题，在直角坐标中，其导热微分方程为：其节点方程为：4.3.1时间的三种差分格式1.向前差分2.向后差分3.中心差分§4.3非稳态导热问题的数值解法K-1时刻SWENPSWENPSWENPτxyK时刻K+1时刻ΔxΔy前者为时间上的向前差分，而后者为时间上向后差分。控制体单位时间的内能增加量为：显式差分格式（采用时间的向前差分）隐式差分格式（采用时间的向后差分）显式差分格式：最突出的优点是节点温度表达式的右边只涉及K时刻的节点温度值，那么只要知道这一时刻周围节点的温度值就可以求出该节点的下一时刻的温度值；隐式差分格式：温度表达式的两端都是同一K时刻的节点温度值，这就意味着必须同时计算同一时刻所有节点的温度值，即必须联立求解K时刻所有节点的差分方程组2控制体的边界节点方程SNWPT∞ΔyΔxα1）对流换热边界条件下的节点方程边界节点P，其控制体的热平衡关系式为SNWPT∞ΔyΔxhΔx=ΔySNWPT∞ΔyΔxh显式差分格式Δx=ΔySNWPT∞ΔyΔxh隐式差分格式2）绝热边界条件下的节点方程。SNEΔx绝热边界ΔyP绝热边界节点P，其控制体的热平衡关系式为SNEΔx绝热边界ΔyPΔx=ΔySNEΔx绝热边界ΔyP显式差分格式Δx=ΔySNEΔx绝热边界ΔyP隐式差分格式§4.4离散方程的求解写出所有内节点和边界节点的温度差分方程n个未知节点温度，n个代数方程式：代数方程组的求解方法：直接解法、迭代解法直接解法：通过有限次运算获得代数方程精确解;矩阵求逆、高斯消元法迭代解法：先对要计算的场作出假设、在迭代计算过程中不断予以改进、直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。缺点：所需内存较大、方程数目多时不便、不适用于非线