有限体积法基本原理教学课件chap10_网格质量.pptxVIP

大连理工大学工业装备结构分析优化与CAE软件全国重点实验室有限体积法基本原理夏广庆,张军军,陈冲,韩亚杰,鹿畅

10.1CFD程序常见网格类型在现代CFD程序中，为了处理各种不同形状的几何模型，常需要将计算域分割为各类不同形状的网格单元。常见的单元类型，在2D情况下，四边形、三角形；在3D情况下，有六面体、棱形、四面体、金字塔型、多面体型，如下图所示。因此，网格质量将决定程序计算结果的准确性与稳定性，采用较高质量的网格有助于仿真程序收敛。在构建CFD网格时，一般认为生成网格中各单元应尽量接近完美正方形和完美（等边）三角形。

10.1CFD程序常见网格类型在现代CFD程序中，为了处理各种不同形状的几何模型，常需要将计算域分割为各类不同形状的网格单元。常见的单元类型，在2D情况下，四边形、三角形；在3D情况下，有六面体、棱形、四面体、金字塔型、多面体型，如下图所示。因此，网格质量将决定程序计算结果的准确性与稳定性，采用较高质量的网格有助于仿真程序收敛。在构建CFD网格时，一般认为生成网格中各单元应尽量接近完美正方形和完美（等边）三角形。

10.2网格质量评价指标网格划分后，通过观察可以初步确定网格划分质量：划分的网格形状与对应正多边形越接近，则网格质量认为越高，如下图所示。较为科学的网格质量评价方法为，通过不同质量评价指标对网格质量进行评价，常用网格质量评价指标有：长宽比（AspectRatio）雅可比行列式（JacobianDeterminant）非正交性（Non-Orthogonality）角度扭曲（EquiangleSkewness）

10.3长宽比对于规则的四边形（所有边都互成90°），长宽比可以直接计算。如图10–4所示，长宽比是最长边和最短边的长度之比，长宽比为1表示完美正方形。对于倾斜四边形和其他类型的单元，如何定义长宽比并不清楚。因此，不同的网格生成器和CFD程序使用不同的长宽比定义。为求解器的稳定性，该比值尽可能接近1。这是因为大长宽比的单元具有比其他单元面积大得多的面，这些面的通量贡献比其他面大很多，这可能会导致系统不稳定。

10.3长宽比——AnsysFluent在ANSYSFluent中，单元中心点到所有面中心点和节点的距离都是单独计算的。图10–5显示了这种计算方法：长宽比取自最大距离（??）与最小距离（??）之比。对于完美正方形，这种方法无法计算出1的长宽比。如图10–6所示，对角线的长度（从中心点到角节点的距离）大于边长（从中心点到面中心的距离），因此计算出的长宽比大于1。尽管存在这种差异，但这种方法是完全通用的，可用于倾斜单元和具有任意数量面的单元。因此，长宽比计算中的细微差别并不重要。

10.3长宽比——OpenFOAMOpenFOAM使用不同的方法计算2D面和3D单元的长宽比。对于2D面，长宽比使用构成该面的最长边和最短边之比来计算。要计算3D单元的长宽比，首先要在单元周围构建一个边界框，这与ANSYSFluent采用的方法不同，边界框是一个长方体，它完全包含单元。因此，长宽比计算是通用的，可用于任何多面体单元格，随后通过最大边界框面积与最小边界框的面积比代表单元长宽比。

10.4雅可比行列式在CFD网格中，需要确保没有一个单元是颠倒（凹陷）的，颠倒单元中的一个或更多节点向内塌陷到单元中心，这会导致内角大于180°。雅可比行列式，通常简称为“Determinant”或“Determinant2×2×2”）用于检查是否有任何单元面倒置。雅可比行列式判定值小于0表示该单元有1个或多个倒置面，而该值大于0则表示没有倒置面。一般建议网格中所有单元的雅可比行列式均大于0.2。有限体积单元是线性的，有限元网格可以使用高阶元素，从而提高求解的精确度。由于多了一个节点，雅可比行列式的计算略有不同。高阶单元的的雅可比行列式计算称为“Determinant3×3×3”（这里的3指的是沿面的节点数）。CFD程序中使用的线性有限体积单元不包含沿面的附加节点。因此，“Determinant2×2×2”和“Determinant3×3×3”的计算结果相同。

10.5非正交性直接影响CFD程序稳定性的主要指标是面的非正交性。图10–11给出了一个网格实例：该网格由一个三角形与一个四边形组成，这两个单元通过公共面相连。?

10.5非正交性得到所有网格面的非正交性后，单元的非正交性即为组成该单元的所有面的非正交性的最大值，一些CFD程序和网格生成器将单元的非正交性表示为非正交质量（??qual