北航计算流体力学大作业创新.docVIP

网格生成方法及网格质量控制 （文献综述） 院 系：能源与动力工程学院 姓 名： 学 号： 指导老师：宁方飞 一、前言 有限元网格生成是工程科学与计算科学相交叉的一个重要研究领域，在经历了30多年发展后的今天依然十分活跃一方面，有限元法己成为一种能够有效地求解各类工程和科学计算问题的通用数值分析方法：另一方面，计算机硬件运算能力的不断提高也容许人们对工程和科学计算的规模、复杂度、效率、精度等方面提出更高的要求。作为有限元走向工程应用的桥梁的有限元网格生成由此获得了源源不断的外在动力。同时，有限元网格生成算法研究中的某些难点问题始终未能获得真正意义上的解决，它们的研究解决对计算几何与计算数学都具有重要的理论价值。 有限元网格生成方法研究领域己取得许多重要成果，形成了独特的方法论体系，提出了许多有效的算法并研制出一些成功的工程化软件产品。近10年来，有限元网格生成方法研究不断地深入、完善和发展，适应性强、应用范围广泛的方法研究重点由二维平面问题转移到三维曲面和三维实体问题，从三角形/四面体网格自动生成转移到四边形/六面体网格自动生成，在并行网格生成、自适应网格生成、贴体坐标网格生成、各向异性网格生成等方面亦取得许多重要进展。 网格质量数值求解效率、收敛性和精度的影响。网格生成后后处理。针对复杂计算区域计算区域获得高质量的计算网格 二、网格生成方法 对不规则区域中的流动与传热问题进行数值计算，首先要解决如何进行区域离散化问题现在有多种对不规则区域进行离散生成计算网格的方法，统称为网格生成技术。本章主要详细介绍结构与非结构网格生成技术。 2.1 概述 积分区域的网格划分直接影响到方程离散的难易，数值计算的快慢和所需计算机内存的大小，也影响到数值解的收敛性和准确性。数值计算中所用网格按网格节点排列是否有序，可分为结构化网格和非结构化网格。在一个区域中，网格的形式可以是单一的，也可以是几种形式的组合。网格生成的详细分类见图。分解和映射法前沿推进法Delaunay三角化 图1 网格生成技术分类 2.2 结构化网格生成方法 所谓结构化网格是指，它的网格所有节点的连接关系是规则的，例如，网格中的所有网格点能够通过索引方便的查找和确定。生成边界固定网格的最常用方法是需要产生一个连续的，并在边界上完全符合的网格。也就是把一组连续毗邻的矩形可计算区域通过曲线边界映射到一个实际的物理区域，如图2。 图2 复杂区域网格的映射 本介绍和分析基于简单区域上的网格生成方法。应用复变函数法构造的网格光滑性较好，在二维翼型计算等方面曾有过广泛的应用，但由于该方法仅限于解决二维问题，且适用的范围较狭小，已逐渐被新的网格生成方法所取代。这里主要讨论的是代数法和微分方程法。 2.2.1代数法 是空间映射方法的一种，在这种方法中我们把物理空间上的坐标作为可计算空间坐标的一个函数，首先在区域的边界开始通过插值方法映射边界上的坐标。由于计算的数据是由区域中无数个点所给出的，所以我们把这种插值算法叫做限插值法。 但是这种函数映射必须满足两个要求(1)这种映射必须是一一映射；(2)计算区域的边界必须映射到物理区域的边界上，并且还存在 二维情况下，是把物理坐标看成是计算坐标的函数，通过物理区域边界上的离散点求该区域内部的点。就是通过两个单向插值的和与它们的积来实现的。两个单向插值函数分别为 这两个映射的积映射为 由布尔和运算定理，最终的TFI函数映射关系为 公式中，矢量表示计算区域边界上的点对应的物理区域边界上的点，两种区域边界上的点的对应关系可预先假定也就是说，它们的对应关系己知，那么通过公式可以确定计算区域内部任意一点处的物理区域的对应点的位置矢量，这样，物理区域的网格便可以计算得到。 图2. 就是采用上面介绍的二维超限插值方法在叶轮机两个直叶片之间所生成的结构化网格示意图。 图2.TFI 方法在叶轮机两个直叶片之间生成的结构化网格 2.2.2微分方程法 下面介绍的这些网格生成方法都是基于一个思想，根据可计算空间坐标系统，为物理空间坐标系统定义一组偏微分方程，然后在可计算空间网格上求解这些方程组，从而得到在物理空间上的对应网格。最常用的偏微分方程是椭圆形偏微分方程，这种方程适用于边界封闭的区域对于边界不封闭的区域，我们会假设一个更大的边界来作为区域的边界进行处理。另外还有两种方程就是双曲线方程和抛物线方程这两种方程主要用于边界不封闭的区域。 椭圆形方程网格生成 这些方法首先对可计算空间坐标和定义一个到物理空间坐标，和的偏微分方程这样的方法可以推广到三维空间。最简单的例子就是拉普拉斯方程 在二维情况下，极值理论保证了这样的映射关系是