第2章 计算区域剖分与控制方程离散.pptVIP

第二章 区域与方程离散 第二章 计算区域与控制方程离散方法 2.1 数值传热学的基本思想 2.2 区域剖分 2.3 离散方法 2.31Taylor 级数与多项式拟合法 2.32控制容积积分与元体平衡法 2.4 四项基本原则    Taylor 展开（续） 2.32 control volume integration 2.32 control volume integration 2.32控制容积积分control volume integration 2.32控制容积积分control volume integration 2.32元体平衡法 2.4离散方程的四项基本原则 界面通量连续原则 正系数原则  其基本思想： 流入到微元体的+微元体产生的=微元体内的变化量 一旦对界面上的变量值和导数值给定计算方法，上式最终也转化为代数方程。 由此可以看出，无论用什么样的离散技术，最终目标都是把定义在数学点上的偏微分方程离散成定义在控制容积上的代数方程，遍历所有控制容积（或节点）得到整个区域的代数方程组。 1. 控制容积截面上总通量连续原则 微分方程离散结果为： 并不是总能得到合理的解，必须满足适当的原则。 4.相邻节点系数关系准则 3.源项的负斜率原则 2.正系数原则 界面通量连续原则是离散过程中必须保持的，又是非常容易出错的地方。例如，有时采用二次插值，用到三个节点构造。如图所示，虚线截面是P，E控制容积公共界面。但W，P，E三点和P，E，EE三点构造的多项式是不同的。因此界面处的一阶和二阶导数是不同的，故扩散通量就不同。 界面通量不连续，就是在方程中无形地引入了一个源项。其结果已不是原方程的了，并且常导致发散 * 2.1 The idea of numerical heat transfer 描述流体流动与传热现象的是一组非线性方程，它们之间互相影响和制约。除非特别的情况，一般难于得到解的解析函数。因此，为分析物理现象的具体规律，求得离散函数解来代替解析函数解。离散解首先把求解区域离散，再把控制方程离散，建立关于各离散点上变量的值之间的制约关系。这种方法在数学上称为求数值解。求解的第一步是把微分方程的有定义的区域离散或剖分。 2.2 Domain decomposition 区域：就是控制方程有定义的时空区间。 区域离散就是把求解区域划分成不重叠无裂缝的子区域之和。对于一维情况就是分布在求解区间上的离散点（节点），每个离散点代表一定的子区域，称为控制容积。离散点或子区域的划分有两种方法，即内点法和外点法。 外点法就是先布置节点，在两个节点中间位置为控制容积的界面。 内点法则是先布置控制容积的界面，控制容积的中心作为节点。 节点位于界面的中间位置。 内点法 外点法 内点法与外点法的剖分过程 比较这两个点差别 对于二维三维空间也有类似离散方法。例如，下面几个是区域离散的例子。 多边形网格 适体网格 三角形网格 当网格单元的相邻关系和拓扑关系比较简单时，网格称为结构化网格，结构化网格分一般结构化网格和贴体网格（BFC）。当拓扑关系复杂时，称为非结构化网格，例如多边形网格和三角形网格。 网格元素与拓扑关系 为什么？ 区域剖分的任务是生成网格系统，网格系统的元素有：网格线、节点和单元，（line , Node and Cell）。外点法的节点就是网格线的交点。内点法就是把节点布置在单元中心。此两种方法分别称为方法A和B 外点法用网格线的交点作为节点，其控制容积为以格点为中心，左右上下各取网格线间距的一半形成的矩形区域（或四边形）为控制容积，所有节点都有控制容积，位于边界上的节点的控制容积是内部节点的控制容积的一半，角点为1/4，界面上的通量计算体现得较好。 内点法生成的网格以网格线为控制容积的边界（界面），节点位于控制容积的中心，边界节点没有控制容积。源项和不稳态项体现得较好。 内点法与外点法对比： （1）控制容积特性上看，边界点。。。 （2）控制容积代表来看， （3）从材料物性突变来看， （4）从导数计算来看， 区域剖分的结果是什么？ 整个区域剖分成了不重叠、无缝隙的子区域之和。用符号表示就是： 2.3方程的离散方法 1、模型方程 为不使问题复杂化，针对一个有源的一维非稳态对流扩散问题作为建立离散方程的例子。 非守恒型： 守恒型： 2、Taylor 展开法差分方程建立 O（x）表示截断误差（truncation error）。 Taylor 展开（续） 后差分 中心差分 前差分