有限体积法()教学文稿.pptVIP

有限体积法;有限差分法;有限差分法;FVM: Finite Volume Method 有限体积法又称为控制体积法,从物理量守恒这一基本要求出发提出的。 其基本思路是： 将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积； 将待解的微分方程对每一个控制体积积分，便得出一组离散方程。 其中的未知数是网格点上的因变量的数值。 为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值的分段的分布剖面。;发展情况; Patanker与Spalding于1967年发表了求解抛物型流动的P-S方法。 在P-S方法中，把x-y平面上的计算区域（边界层）转换到x-w平面上（w为无量纲流函数），从而不论在边界层的起始段还是在其后的发展段，所设置的计算节点均可落在边界层范围内。 1969年Spalding在英国帝国理工学院（Imperial College）创建了CHAM（Concentration, Heat and Mass, Limited），旨在把他们研究组的成果推广应用到工业界。; 1972年SIMPLE算法问世 于是所谓分离式的求解方法应运而生，即先求解有关一个速度分量，而把其他作为常数，随后再逐一求解其它变量。 于是就产生了这样的问题：就是所谓速度与压力的耦合问题。SIMPLE算法成功地解决了这一问题。 SIMPLE算法的一个基本思想是，在流场迭代求解的任何一个层次上，速度场都必须满足质量守恒方程，这是保证流场迭代计算收敛的一个十分重要的原则。 ;1977年由Spalding及其学生开发的ENMIX程序公开发行。 1979年由Spalding教授及其合作者开发的流动传热计算的大型通用软件PHOENICS第一版问世。 PHOENICS是英语Parabolic, Hyperbolic or Elliptic Numerical Integration Code Series的缩写（意为对抛物型、双曲型、椭圆型方程进行数值积分的系列程序）。 *1980年Patankar教授的名著“Numerical Heat Transfer and Fluid Flow”出版。 这本书内容精炼，说理透彻，注重物理概念的阐述，深受全世界数值传热的研究者与使用者的欢迎。出版后不久，被相继译成俄文、日文、波兰文及中文等，成为数值传热学领域中的一本经典著作;优点: (1)出发点是积分形式的控制方程,同时积分方程表示了特征变量在控制容积内守恒; (2)积分方程中每一项都有明确的物理意义,从而使方程离散时,对各离散项可以给出一定的物理解释; (3) 区域离散的节点网格与进行积分的控制容积分立。 ;方法特点;基本概念;基本概念;控制体积的形成 单元中心方式： 单元位于控制体积的中心，即将单一的网格单元作为控制体积，网格单元互不重叠。 单元顶点方式： 以网格节点为中心，它的形成有多种方式。其常用的构成方式是由以该节点为顶点的格子的形心以及各共顶点的网格线中心点的一系列连接线段所构成一个多边形控制体积。也可以由环绕该节点的一组格子组成控制体。 ;单元顶点方式;对于计算同样多的变量，单元中心方式变量布置简单直观，易于处理边界条件和保持离散的守恒性，而且需要的网格数要比单元顶点方式少得多，可节省计算时间。 ;网格;结构网格 矩形网格最为常用，网格生成方便，但对复杂边界处理过于粗糙； 曲线网格：1）只提供了离散点的变换，而不给出解析函数形式的变换关系,使用不光滑的网格时，对变换关系的差分近似会造成了很大的数值误差，甚至会导致不切实际的值。2)如果网格严重偏离正交性，就会极大损坏原有的迭代方法的收敛速率。3)因变量的选择也须谨慎考虑。在曲线网格中，可取原始笛卡尔坐标系变量或曲线坐标系中沿网格方向的协变量两种作为因变量。 ;非结构网格 非结构网格中单元格分布不再规则一致，其位置很难再凭借行列索引关系确定。非结构网格可以采用任意形状的单元格，单元边的数目也无限制，弥补了结构化网格不能够解决任意形状和任意连通区域的网格剖分的缺欠。 ;非结构网格最重要的一个特征是控制方程离散得到的代数方程的系数矩阵不再是结构网格下有规律的对角结构；若用对角形式存储，其带宽只能通过适当的布置单元编号顺序来减少。非结构网格原则上可应用于任何类型的数值方法，但非结构网格的FVM 算法更成熟，应用更广。 ;●非结构网格最早用于FEM； ●但流体流动是高度非线性问题，而且FEM 计算量较大，这些问题使得基于FEM 的非结构网格技术未能在对流问题为主的地面水流（如浅水流动，水波运动等）计算上得到重视； ●八十年代以来，基于FVM 的非结构网格技术在空气动力学得到了广泛的发展和应用； ●九十年代开始一些专家学者根据浅水流动特征，将这些算法引入到