y-chapter5.pptVIP

Advanced Heat Transfer North China Electric Power University 第五章 导热问题的数值解法 内容简介： 本章介绍了数值求解物理问题的基本思想，以二维稳态导热问题和一维非稳态导热为例，介绍了采用有限差分法数值求解导热问题的基本过程。 5-0 数值求解的基本思想 一、 数值解法的本质 数值解法是用物理问题所涉及的空间和时间区域内有限个离散点(称为节点)的物理量近似值来代替物体内实际连续的物理量分布，将连续物理量分布函数的求解问题转化为各节点物理量值的求解问题。 二、主要的数值解法 有限差分法（FDM），主要介绍的方法 有限元法（FEM） 边界元法（BEM） 三、数值解法的基本步骤 给出问题的数学描述 确定节点（区域离散化） 建立节点物理量的代数方程 求解代数方程 解的分析 1.给出问题的数学描述 对问题进行分析，确定问题的类型，并进行合理的简化，给出正确的数学描述。比如下面的问题 y x 0 H W h,tf q=0 h,tf q=0 将求解区域按照一定规则划分为许多小的小区域，这个过程称作区域的离散。每个小的区域的物理量值由一个点—节点来表示。若是非稳态问题，还需将时间离散化。 2.离散化（区域和时间的离散） x y m n 每一个节点都与它周围相邻的节点存在一定的关系，通过相应的物理定律，可建立它们之间的关系式（属于代数方程式），此关系式又称作节点的离散方程。 3.建立节点物理量的代数方程 (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) 把所有节点的离散方程联立起来，会组成一个封闭的方程组，对代数方程组的求解可采用直接解法或迭代求解，更多的是采用迭代解法。 4.求解代数方程组 对解的结果进行进一步的分析，获得定性或定量上的一些新的结论。 5.进行解的分析 四、数值解法的应用 数值解法是目前各行各业解决工程问题和进行科学研究的一种重要方法。如：航空航天、冶金、电力、电子、化工、建筑、能源动力、机械等。 五、大型的数值计算商用软件 目前市场上有许多针对某一领域的商用软件，在我们专业应用的有： FLUENT——是目前国际上比较流行的商用CFD软件包，在美国的市场占有率为60%。凡跟流体，热传递及化学反应等有关的工业均可使用。 CFX——是由英国AEA公司开发，是一种实用流体工程分析工具，用于模拟流体流动、传热、多相流、化学反应、燃烧问题。 此外还有：Phoenics和 STAR-CD 等。 5-1 导数的有限差分近似表达式 有限差分的数学基础是用差商代替微商，即用有限差分代替导数。 利用泰勒级数可以得到更多形式的微分的差商表示，也可以看到其误差。 改写，可得 式中 是用差商代替微商的误差。 此种格式被称为向前差分 同样： 改写，可得 此种格式被称为向后差分 可得一阶导数的中心差分表达式 如果将下面两式相减， 可得二阶导数的中心差分表达式 如果将下面两式相加， 5-2 稳态导数的数值分析 下面以一个有内热源的二维稳态导热问题为例。 y x 0 H W h,tf q=0 h,tf q=0 一、给出问题的数学描述 二、区域的离散 将二维矩形区域按照下图的格式进行离散化，相邻两个节点间的距离称为步长，记为Δx、Δy 。区域内部的节点称为内节点，边界上的节点称为外节点。 x y m n 三、建立节点的离散方程 1．内部节点 （1）微商的差分表示法 如对于图中的节点(m，n)，导热微分方程表示为： (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) 利用二阶导数的中心差分表达式 可得 (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) （2）热平衡法 热平衡方法的思想是，对每个节点所代表的控制体列能量守恒方程式，从而得出该点与其它节点的关系式。 假定节点（m,n）所代表的控制区域四个边界分别为n、s、e、w ，节点（m,n）所代表的控制体的能量守恒方程为： n e s w (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) 节点（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．边界节点离散方程的建立 边界节点离散方程采用热平衡的方法建立更方便。 对于平直边界（如右边）上的