1. 网格基础
1.1 计算域模型
流体计算域指的是在流体计算过程中,参与积分计算的区域。
根据流动区域的不同,可以将计算域分为内流计算域、外流计算域和混合(内外流混合)计算域。
1.2 计算域生成方法
- 直接建模:直接采用CAD等几何建模软件构建。
- 几何抽取:对于内表面复杂的计算域几何模型,通常采用几何抽取生成计算域模型。
1.3 计算域简化
为了减少计算开销,通常根据问题的特征对几何模型进行简化处理。常见的几何简化包括:利用对称性、3D模型转化为2D模型、利用流动的周期性等。
1.4 多区域计算模型
多区域计算模型指的是计算模型中包含有两个及两个以上计算域。
多区域计算模型主要应用于以下情况:
- 计算模型中涉及运动区域,如旋转机械模拟仿真。
- 计算模型中同时存在固体或流体区域,如共轭传热仿真等。
- 计算模型中存在多孔介质区域或其他需要单独求解的区域。
- 为了网格划分方便,将计算区域切分成多个区域。
1.4.1 Interface——边界
对于多区域计算模型,一般CFD求解器均提供了数据传递方式,用户只需将区域界面进行编组处理,即可完成计算过程中数据的传递。CFD区域分界面数据传递主要采用定义Interface对来完成。Interface主要用于处理多区域计算模型中区域界面间的数据传递。
Interface是边界类型的一种,这意味着Interface是计算域的边界,因此,当计算模型中存在多个计算区域时,若想计算域保持流通,则需要在相互接触的边界上创建Interface。
在CFD计算中,Interface通常都是成对出现的,计算结果数据则通过Interface对进行插值传递。
利用Interface并不要求边界上的网格节点一一对应。
1.4.2 Interior——内部面
Interior指的是内部面,常出现在单计算区域中。在对单计算域进行分区网格划分时,尤其是划分不同类型计算网格时,网格分界面将会被求解器识别为Interior类型。
需要注意的是:与Interface不同,Interior不是计算区域边界,而是计算域内部网格面。Interior面上不能有重合的网格节点。而Interface对上的网格节点既可以重合,也可以不重合。
形象的说法:Interface是两个独立区域边界,是实际存在的边界。而Interior则常常是虚拟形成的。默认情况下,Interior是连通的,而Interface则是非连通的,需要在求解器中设置Interface对才能使计算域保持连通。
1.5 网格基础概念
- 网格:Grid、Cell、Mesh都指网格,网格通常指计算域离散后形成的封闭体积。
- 节点:Node、Vertices,其中固体常用Node,流体常用Vertices,节点指离散计算域的分割线的焦点。
- 网格形状:
2D模型中,常见的网格类型包括三角形网格和四边形网格;
3D模型中,常见的包括四面体和六面体网格、棱柱网格、金字塔网格、多面体网格。 - 结构网格:在拓扑结构上相当于矩形域内的均匀网格,其节点定义在每一层的网格线上,且每一层上节点数都是相等的,这样使复杂外形的贴体网格生成比较困难。
- 非结构网格:没有规则的拓扑结构,也没有层的概念,网格节点的分布是随意的,因此具有灵活性。但是非结构网格计算的时候需要较大的内存。
1.6 网格的度量
1.6.1 网格数量
2D网格由网格节点、网格边及网格面构成;3D网格由网格节点、网格边、网格面、单元体构成。
通常所说的网格数量指网格节点数量以及网格面(2D网格中)或网格体(3D网格体)。网格数量对计算的影响主要体现在以下几个方面。
- 网格数量越多,需要的计算资源(内存、CPU时间、硬盘等)越大。
- 并非网格数量越多,计算越精确。对于物理量变化剧烈区域,采用局部网格加密可以提高该区域计算精度。同时需要进行网格独立性验证。
- 影响计算收敛性的因素是网格质量,而不是网格数量。对于一些瞬态计算,时间步长与网格尺寸有关系。小的网格尺寸意味着需要更加细密的时间步长。
1.6.2 网格质量
- 角度(Angle):用于度量网格边之间的夹角。角度范围0°~90°,0°表示单元退化的网格(质量差),90°为完美网格。角度度量标准比较常用,CFD计算通常要求角度大于18°,但是一些不太敏感区域,14°以上也在可接受范围内。
- 纵横比(Aspect Ratio):主要用于六面体网格,定义为单元最大边长度与最小边长度的比值。其中纵横比为1为完美网格。纵横比最好是限定在20以内,在CFD计算中,只有边界层网格允许较大的纵横比。过大的纵横比会引入较大的计算误差,甚至会导致计算发散。
- 行列式(Determinant 2×2×2):或者说是相对行列式,定义为最大雅克比矩阵行列式与最小雅克比矩阵行列式的比值。正常网格取值范围为0~1。值为1表示为完美网格,值越低表示网格越差。负值表示存在负体积网格,不能被求解器接受。该评判标准应用较多。
- 行列式(Determinant 3×3×3):该评判指标用于六面体网格。与2×2×2不同的是,单元边上的中心点会被增加至雅克比计算中。
- 最小角(Min Angle):计算每一个网格单元的最小内角。值越大表示网格质量越好。
- 质量(Quality): 质量是ICEM CFD中用于标定网格质量的衡量指标,对于不同类型的网格,其采用不同的衡量方式。
- 三角形或四面体网格:计算高度与每一条边的长度比值,取最小值。越接近1,网格质量越好。
- 四边形网格:网格质量利用行列式Determinant 2×2×2进行度量。
- 六面体网格:计算三种度量方式(行列式、最大正交性、最大翘曲度),取最小值作为网格质量评判标准。
- 金字塔网格:采用行列式进行评判。
- 棱柱网格:计算行列式与翘曲度,取最小值作为质量评判指标。