Fluent教程--第05章 fluent网格.pdf

网格的读入和使用网格的读入和使用 FLUENT 可以从输入各种类型，各种来源的网格你可以通过各种手段对网格进行修 改，如：转换和调解节点坐标系，对并行处理划分单元，在计算区域内对单元重新排序以减 少带宽以及合并和分割区域等 你也可以获取网格的诊断信息， 其中包括内存的使用与简化， 网格的拓扑结构，解域的信息你可以在网格中确定节点、表面以及单元的个数，并决定计 算区域内单元体积的最大值和最小值， 而且检查每一单元内适当的节点数 以下详细叙述了 FLUENT 关于网格的各种功能 （请参阅网格适应一章以详细了解网格适应的具体内容 ） 网格拓扑结构 FLUENT 是非结构解法器，它使用内部数据结构来为单元和表面网格点分配顺序，以 保持临近网格的接触因此它不需要 i，j，k 指数来确定临近单元的位置解算器不会要求 所有的网格结构和拓扑类型， 这使我们能够灵活使用网格拓扑结构来适应特定的问题 二维 问题，可以使用四边形网格和三角形网格，三维问题，可以使用六面体、四面体，金字塔形 以及楔形单元，具体形状请看下面的图形FLUENT 可以接受单块和多块网格，以及二维 混合网格和三维混合网格另外还接受 FLUENT 有悬挂节点的网格（即并不是所有单元都 共有边和面的顶点） ，有关悬挂节点的详细信息请参阅“节点适应”一节。

非一致边界的网 格也可接受（即具有多重子区域的网格，在这个多重子区域内，内部子区域边界的网格节点 并不是同一的） 详情请参阅非一致网格 Figure 1: 单元类型 可接受网格拓扑结构的例子 正如网格拓扑结构一节所说， FLUENT 可以在很多种网格上解决问题 图 1—11 所示为 FLUENT 的有效网格O 型网格，零厚度壁面网格，C 型网格，一致块结构网格，多块结构 网格， 非一致网格， 非结构三角形， 四边形和六边型网格都是有效的 Note that while FLUENT does not require a cyclic branch cut in an O-type grid， it will accept a grid that contains one. Figure 1: 机翼的四边形结构网格 Figure 2:非结构四边形网格 Figure 3: 多块结构四边形网格 Figure 4: O 型结构四边形网格 Figure 5: 降落伞的零厚度壁面模拟 Figure 6: C 型结构四边形网格 Figure 7:三维多块结构网格 Figure 8: Unstructured Triangular Grid for an Airfoil Figure 9:非结构四面体网格 Figure 10:具有悬挂节点的混合型三角形/四边形网格 Figure 11:非一致混合网格 for a Rotor-Stator Geometry 选择适当的网格类型 FLUENT 在二维问题中可以使用由三角形、四边形或混合单元组成的网格，在三维问 题中可以使用四面体，六面体，金字塔形以及楔形单元，或者两种单元的混合。

网格的选择 依赖于具体的问题，在选择网格的时候，你应该考虑下列问题： ? 初始化的时间 ? 计算花费 ? 数值耗散 后面将会详细讨论各种类型网格的特点 初始化的时间 很多实际问题是具有复杂几何外形的， 对于这些问题采用结构网格或块结构网格可能要 花费大量的时间， 甚至根本无法得到结构网格 复杂几何外形初始化时间的限制刺激了人们 在非结构网格中使用三角形网格和四面体网格然而，如果你的几何外形并不复杂的话，两 种方法所耗费的时间没有明显差别 如果你已经有了结构网格代码如 FLUENT 4 生成的网格，那么在 FLUENT 中使用该网 格会比重新生成网格节约大量的时间这一特点也刺激了人们在 FLUENT 仿真中使用四边 形网格和六面体网格注意：FLUENT 有一个格式转换器允许你从其它程序中读入结构网 格 计算花费 当几何外形太复杂或者流动的长度尺度太大时， 三角形网格和四面体网格所生成的单元 会比等量的包含四边形网格和六面体网格的单元少得多 这是因为三角形网格和四面体网格 允许单元聚集在流域的所选区域， 而四边形网格和六面体网格会在不需要加密的地方产生单 元 非结构的四边形网格和六面体网格为对于一般复杂外形提供了许多三角形和四面体网格 的优点。

四边形和六边形单元的一个特点就是它们在某些情况下可以允许比三角形/四面体单元 更大的比率三角形/四面体单元的大比率总会影响单元的歪斜因此，如果你有相对简单 的几何外形，而且流动和几何外形很符合，比如长管，你就可以使用大比率的四边形和六边 形单元这种网格可能会比三角形/四面体网格少很多单元 数值耗散 多维条件下主要的误差来源就是数值耗散又被称为虚假耗散 （之所以被称为虚假的， 是 因为耗散并不是真实现象，而是它和真实耗散系数影响流动的方式很类似） 关于数值耗散有如下几点： ? 当真实耗散很小时，即对流占主导地位时，数值耗散是显而易见的 ? 所有的解决流体问题的数值格式都会有数值耗散，这是因为数值耗散来源于截断误差， 截断误差是描述流体流动的离散方程导致的 ? FLUENT 中所用的二阶离散格式可以帮助减少解的数值耗散的影响 ? 数值耗散量的大小与网格的分辨率成反比 因此解决数值耗散问题的一个方法就是精化 网格 ? 当流动和网格成一条直线时数值耗散最小（所以我们才要使用结构网格来计算啊） 最后一点和网格选择最有关系很明显，使用三角形/四面体网格流动永远不会和网格 成一条直线， 而如果几何外形不是很复杂时， 四边形网格和六面体网格可能就会实现流动和 网格成一条线。

只有在简单的流动，如长管流动中，你才可以使用四边形和六面体网格来减 少数值耗散，而且在这种情况下使用四边形和流面体网格有很多优点，因为与三角形/四面 体网格相比你可以用更少的单元得到更好的解 网格所需条件和所要考虑的问题 本节讨论了特殊几何图形和网格的必要条件以及网格质量的一般评价方法 几何图形和网格的必要条件 在计划解决你的问题的开始，应该注意下面的几何图形设定以及网格结构的必要条件 ? 对于轴对称图形来说，必须定义笛卡尔坐标系的 x 轴为旋转轴 (Figure 1). Figure 1:轴对称图形必须以 x 轴为中线 ? 周期性边界条件要具有周期性网格，虽然 GAMBIT 和 TGrid 能够产生真正的周期性边 界， 但是 GeoMesh 和大多数 CAD 软件包是无法产生周期性边界条件的 如果下面的条 件需要满足的话，TGrid 提供了 GeoMesh 和大多数 CAD 软件产生的三角形表面网格生 成周期性边界的功能 1. 周期及其内部在它们的边界曲线上有相同的节点分布 2. 周期及其内部的节点与常数平动因子和转动因子有关 详情请见 GAMBIT 和 TGrid 的帮助文件 如果你用 GeoMesh 和大多数 CAD 软件产生四边形网格和六面体网格， 你必须保证在周 期性区域内的网格是相同的。

然后便可以在 FLUENT 中使用 make-periodic 命令建立周期性 边界详细内容请参阅“创建周期性区域”一节你能够在解算器中对三角形或四面体网 格创建周期性边界条件而不用上面所述的 TGrid 来创建) 网格质量 网格质量对计算精度和稳定性有很大的影响网格质量包括：节点分布，光滑性，以及 歪斜的角度（skewness） 节点密度和聚集度 连续性区域被离散化使得流动的特征解 （剪切层， 分离区域， 激波， 边界层和混合区域） 与网格上节点的密度和分布直接相关 在很多情况下， 关键区域的弱解反倒戏剧化的成了流 动的主要特征比如：由逆压梯度造成的分离流强烈的依靠边界层上游分离点的解 边界层解（即网格近壁面间距）在计算壁面剪切应力和热传导系数的精度时有重要意义这 一结论在层流流动中尤其准确，网格接近壁面需要满足： 1≤ ∞ x p v u y 其中 p y= 从临近单元中心到壁面的距离； ∞ u= 自由流速度 ；v= 流体的动力学粘性系数 ；X= 从边界层起始点开始沿壁面的距离上面的方程基于零攻角层流流动的 Blasius 解[139] 网格的分辨率对于湍流也十分重要 由于平均流动和湍流的强烈作用， 湍流的数值计算 结果往往比层流更容易受到网格的影响。

在近壁面区域， 不同的近壁面模型需要不同的网格 分辨率 一般说来， 无流动通道应该用少于 5 个单元来描述 大多数情况需要更多的单元来完全 解决 大梯度区域如剪切层或者混合区域， 网格必须被精细化以保证相邻单元的变量变化足 够小不幸的是要提前确定流动特征的位置是很困难的而且在复杂三维流动中，网格是要 受到 CPU 时间和计算机资源的限制的在解运行时和后处理时，网格精度提高，CPU 和内 存的需求量也会随之增加自适应网格技术可用于在流场的发展基础上提高和/或减少网格 密度，并因此而提供了网格使用更为经济的方法 光滑性 临近单元体积的快速变化会导致大的截断误差 截断误差是指控制方程偏导数和离散估 计之间的差值FLUENT 可以改变单元体积或者网格体积梯度来精化网格从而提高网格的 光滑性 单元的形状 单元的形状（包括单元的歪斜和比率）明显的影响了数值解的精度单元的歪斜可以定 义为该单元和具有同等体积的等边单元外形之间的差别 单元的歪斜太大会降低解的精度和 稳定性比方说：四边形网格最好的单元就是顶角为 90 度，三角形网格最好的单元就是顶 角为 60 度比率是表征单元拉伸的度量正如在计算花费一节所讨论的，对于各向异性流 动， 过渡的比率可以用较少的单元产生较为精确的结果。

但是一般说来应该尽量避免比率大 于 5:1 流动流场相关性 分辨率、光滑性、单元外形对于解的精度和稳定性的影响强烈的依赖于所模拟的流场 例如： 在流动开始的区域可以忍受过渡歪斜的网格， 但是在具有大流动梯度的区域这一特点 可能会使得整个计算无功而返 因为大梯度区域是无法预先知道的， 所以我们只能尽量的使 整个流域具有高质量的网格 网格的读入 FLUENT 能够处理大量的具有不同结构的网格拓扑结构因此我们有很多产生网格的 工具， 比如： GAMBIT， TGrid， GeoMesh， preBFC， ICEMCFD， I-DEAS， NASTRAN， PATRAN， ARIES， ANSYS， 以及其它的前处理器， 或者使用FLUENT/UNS， RAMPANT， 以及 FLUENT 4 case 文件中包含的网格， 你也可以准备多个网格文件， 然后把它们结合在一 起创建一个网格 GAMBIT 网格文件 你可以使用 GAMBIT 创建二维和三维结构/非结构/混合网格 详细内容请参阅 GAMBIT 建模向导，并将你的网格输出为 FLUENT 5 格式所有的这样的网格都可以直接读入到 FLUENT，菜单：File/Read/Case. GeoMesh 网格文件 你可以使用 GeoMesh 创建二维四边形网格或三角形网格以及三维六面体网格和三维四 面体网格的三角网格面。

具体请参阅 GeoMesh 用户向导要完成三维四面体网格的创建你 必须把表面网格读入到 TGrid 然后产生体网格 其它的网格都可以直接读入到 FLUENT： 菜 单 File/Read/Case. TGrid 网格文件 你可以用 TGrid 从边界或表面网格产生二维或三维非结构三角形/四面体网格具体方 法请参阅 TGrid 用户向导在 FLUENT 中你可以点击 File/Write/Mesh.菜单保存网格读入 网格请点击 File/Read/Case.菜单，具体内容参阅读入网格文件一节 preBFC 网格文件 你可以用 preBFC 产生两种 FLUENT 所使用的不同类型的网格：结构四边形/六面体网 格和非结构三角形/四面体网格下面详细介绍一下 结构网格文件。