计算流体力学有限体积法.pptxVIP

数值方法、计算机科学与计算流体力学的结合--CFD;一、CFD求解流体力学问题的过程;4.建立离散方程;5. 离散初始条件和边界条件; 如线性方程组可以用高斯消去法或高斯-赛德尔迭代法。非线性方程组采用牛顿牛顿-辛普森方法等。;二、CFD数值计算方法和分类;（二）有限元法（FEM）;基本思路和步骤; 有限体积法又称为控制体积法。其基本思路是:将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。 其中的未知数是网格点上的因变量的数值。为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值分段分布的剖面。 从积分区域的选取方法看来，有限体积法属于加权剩余法中的子区域法；从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。简而言之，子区域法属于有限体积法的基本方法。 有限体积法的基本思路易于理解，并能得出直接的物理解释。离散方程的物理意义就是因变量在有限大小的控制体积中的守恒原理，如同微分方程表示因変量在无限小的控制体积都得到满足，在整个计算区域，自然也就得到满足一样，这是有限体积法吸引人的优点。; 某些离散方法，如有限差分法，仅当网格极其细密时，离散方程才满足积分守恒；而有限体积法即使在粗网格情况下，也显示出准确的积分守恒。 就离散方法而言，有限体积法可视为有限单元法和有限差分法的中间物，有限单元法必须假定值符号网格点之间的变化规律(即插值函数)，并将其作为近似解；有限差分法只考虑网格点上的数值而不考虑其在网格点之间如何变化；有限体积法只寻求节点值，这与有限差分法相类似，但有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布，这又与有限单元法相类似。 在有限体积法中，插值函数只用于计算控制体积的积分，得出离散方程后，便可忘掉插值函数；如果需要的话，可以对微分方程中不同的项采取不同的插值函数。 ;三、有限体积法计算区域的离散; 一般把节点看成是控制体积的代表。在离散过程中，将一个控制体积上的物理量定义并存储在该节点处。一维问题的有限体积法计算网格如图1示，二维问题的有限体积法计算网格如图2所示。 ; 结构化网格的节点排列有序，即给出了一个节点的编号后，立即可以得出其相邻节点的编号，所有内部节点周围的网格数目相同。 结构化网格具有实现容易、生成速度快，网格质量好、数据结构简单化的优点，但不能实现复杂边界区域的离散。 非结构化网格的内部节点以一种不规则的方式布置在流场中，各节点周围的网格数目不尽相同。这种网格且然生成比较复杂却有极大的适应性，对复杂算问题特别有效。;?; 有限体积法的关键一步是在控制体积上积分控制方程，在控制体积节点上产生离散的方程。对一维模型方程2，在图1所示的控制体积P上进行积分，有： 式3中，△是控制体积的体积值。当控制体积很微小时，△可以表示为△V?A，这里A是控制体积界面的面积。从而有： 从式4可以看到，对流项和扩散项均已转化为控制体积界面上的值。有限体积法最显著的特点之一就是离散方程中具有明确的物理插值，即界面的物理量要通过插值的方式由节点的物理量来表示。;?;?;对于源项S，它通常是时间和物理量的函数。为了简化处理，将S转化为如下线性方式: 式中，Sc是常数，S是随时间和物理量φ变化的项。将式1-6～式8代入式4，有： 整理后得到： ;记为： 式10中：;?;?; 在此基础上，定义一维单元的 Peclet数Pe为 ; 瞬态问题与稳态问题相似，主要是瞬态项的离散。其一维瞬态问题的通用控制方程为： ;将式16在一维计算网格上对时间及控制体积进行积分，有：;?;?;整理后得:;代入方程22，得到:; 式25中：; 进一步将一维问题扩展为二维与三维问题。在二维问题中，计算区域离散见图2。发现只是增加第二坐标y,控制体积増加的上下界面,分别用n(north)和s(south)表示，相应的两个邻点记为N和S。在全隐式时间积分方案下的二维瞬态对流扩散问题的离散方程为: ; 在三维问题中，计算区域离散如图3所示(两个方向的投影)。在二维的基础上增加第三坐标z，控制体积增加的前后界面，分别用t(top)和b(bottom)表示，相应的两个邻点记为T和B。在全隐式时间积分方案下的三维瞬态对流扩散问题的离散方程为;?;图3 三维计算区域离散网格在两个方向上的投影;?;?; 流场计算的基本过程是在空间上用有限体积法(或其他类似方法)将计算区域离散成许多小的体积单元，在每个体积单元上对离散后的控制方程组进行求解。 其