9-10导热问题的数值解法（第四章）.pptVIP

传热学 Heat Transfer 华北电力大学 第四章 导热问题的数值解法 内容简介： 本章介绍了数值求解物理问题的基本思想，以二维稳态导热问题和一维非稳态导热为例，介绍了采用有限差分法数值求解导热问题的基本过程。 4-1 数值求解的基本思想 一、 数值解法的本质 数值解法是用物理问题所涉及的空间和时间区域内有限个离散点(称为节点)的物理量近似值来代替物体内实际连续的物理量分布，将连续物理量分布函数的求解问题转化为各节点物理量值的求解问题。 二、数值解法的基本步骤 建立控制方程及定解条件 确定节点（区域离散化） 建立节点物理量的代数方程 设立物理量的初值 求解代数方程 是否收敛 解的分析 改进初场 是 否 将求解区域按照一定规则划分为许多小的小区域，这个过程称作区域的离散。每个小的区域的物理量值由一个点—节点来表示。 1.区域的离散 x y m n 每一个节点都与它周围相邻的节点存在一定的关系，通过相应的物理定律，可建立它们之间的关系式（属于代数方程式），此关系式又称作节点的离散方程。 2.建立节点物理量的代数方程 (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) 把所有节点的离散方程联立起来，会组成一个封闭的方程组，对代数方程组的求解可采用直接解法或迭代求解，更多的是采用迭代解法。 3.求解代数方程组 (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) 对解的结果进行进一步的分析，获得定性或定量上的一些新的结论。 4.进行解的分析 三、数值解法的应用 数值解法是目前各行各业解决工程问题和进行科学研究的一种重要方法。如：航空航天、冶金、电力、电子、化工、建筑、能源动力、机械等。 四、大型的数值计算商用软件 目前市场上有许多针对某一领域的商用软件，在我们专业应用的有： FLUENT——是目前国际上比较流行的商用CFD软件包，在美国的市场占有率为60%。凡跟流体，热传递及化学反应等有关的工业均可使用。 CFX——是由英国AEA公司开发，是一种实用流体工程分析工具，用于模拟流体流动、传热、多相流、化学反应、燃烧问题。 此外还有：Phoenics和 STAR-CD 等。 4-2 二维稳态导热问题的数值求解 一、物理问题 接下来以一个二维稳态导热问题为例，介绍数值求解导热问题的的具体过程，重点是节点离散方程的建立和代数方程组的迭代求解。 如图所示：二维矩形区域内的稳态、含有内热源、常物性的导热问题。 y x 0 H W h,tf q=0 h,tf q=0 二、数学描述 三、区域的离散 将二维矩形区域按照下图的格式进行离散化，相邻两个节点间的距离称为步长，记为Δx、Δy 。区域内部的节点称为内节点，边界上的节点称为外节点。 x y m n 四、建立节点的离散方程 1．内部节点 （1）微商的差分表示法 如对于图4-3中的节点(m，n)，导热微分方程表示为： 由微商的定义： 当 较小时，可以用差商近似代替微商 或 同理可得 若Δx=Δy （2）热平衡法 热平衡方法的思想是，对每个节点所代表的控制体列能量守恒方程式，从而得出该点与其它节点的关系式。 假定节点（m,n）所代表的控制区域四个边界分别为n、s、e、w ，节点（m,n）所代表的控制体的能量守恒方程为： (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) n e s w 节点（m,n）所代表的控制体的能量守恒方程为： 从四个边界导入控制体的能量可应用傅立叶定律得到，如界面w ： (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) n e s w 代入能量守恒方程式，并整理得 (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) 2．边界节点离散方程的建立 边界节点离散方程只能采用热平衡的方法建立。 对于平直边界（如右边）上的节点 热平衡法的优点 （1）物理意义明确； （2）适用于内节点和边界节点。 对于其它外节点，如外部角点和内部角点等可按同样的方法推导，区别是其所代表的控制体体积不同。 写出所有内节点和边界节点的温度差分方程n个未知节点温度，n个代数方程式： 代数方程组的求解方法：直接解法、迭代解法 五、节点方程组的求解 直接解法：通过有限次运算获得代数方程精确解; 如：矩阵求逆、高斯消元法等。 迭代解法：先对要计算的场作出假设、在迭代计算过程中不断予以改进、直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。 缺点：所需内存较大、方程数目多时不便