网格的读入和使用 FLUENT可以从输入各种类型,各种来源的网格。你可以通过各种手段对网格进行修改,如:转换和调解节点坐标系,对并行处理划分单元,在计算区域内对单元重新排序以减少带宽以及合并和分割区域等。你也可以获取网格的诊断信息,其中包括内存的使用与简化,网格的拓扑结构,解域的信息。你可以在网格中确定节点、表面以及单元的个数,并决定计算区域内单元体积的最大值和最小值,而且检查每一单元内适当的节点数。以下详细叙述了FLUENT关于网格的各种功能。(请参阅网格适应一章以详细了解网格适应的具体内容。) 网格拓扑结构 FLUENT是非结构解法器,它使用内部数据结构来为单元和表面网格点分配顺序,以保持临近网格的接触。因此它不需要i,j,k,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软件产生的三角形表面网格生成周期性边界的功能。 周期及其内部在它们的边界曲线上有相同的节点分布。 周期及其内部的节点与常数平动因子和转动因子有关。 详情请见GAMBIT和TGrid的帮助文件。 如果你用GeoMesh和大多数CAD软件产生四边形网格和六面体网格,你必须保证在周期性区域内的网格是相同的。然后便可以在FLUENT中使用make-periodic命令建立周期性边界。详细内容请参阅“创建周期性区域”一节。(你能够在解算器中对三角形或四面体网格创建周期性边界条件而不用上面所述的TGrid来创建) 网格质量 网格质量对计算精度和稳定性有很大的影响。网格质量包括:节点分布,光滑性,以及歪斜的角度(skewness)。 节点密度和聚集度 连续性区域被离散化使得流动的特征解(剪切层,分离区域,激波,边界层和混合区域)与网格上节点的密度和分布直接相关。在很多情况下,关键区域的弱解反倒戏剧化的成了流动的主要特征。比如:由逆压梯度造成的分离流强烈的依靠边界层上游分离点的解。 边界层解(即网格近壁面间距)在计算壁面剪切应力和热传
第05章 fluent网格35.doc
下载此电子书资料需要扣除0点,