fluent网格检查2009.docVIP

fluent网格检查2009-05-05 18:371. Fluent检查网格质量的方法，网格导入Fluent中之后，grid-check，可以看看网格大致情况，有无负体积，等等；在Fluent窗口输入，grid quality然后回车，Fluent会显示最主要的几个网格质量。 在这里我谈一下，Fluent计算对网格质量的几个主要要求： 1）网格质量参数： Skewness （不能高于0.95，最好在0.90以下；越小越好） Change in Cell-Size （也是Growth Rate，最好在1.20以内，最高不能超过1.40） Aspect Ratio （一般控制在5：1以内，边界层网格可以适当放宽） Alignment with the Flow（就是估计一下网格线与流动方向是否一致，要求尽量一致，以减少假扩散） 2）网格质量对于计算收敛的影响： 高Skewness的单元对计算收敛影响很大，很多时候计算发散的原因就是网格中的仅仅几个高Skewness的单元。 举个例子：共有112,000个单元，仅有7个单元的Skewness超过了0.95，在进行到73步迭代时计算就发散了！ 高长宽比的单元使离散方程刚性增加，使迭代收敛减慢，甚至困难。也就是说，Aspect Ratio尽量控制在推荐值之内。 3）网格质量对精度的影响： 相邻网格单元尺寸变化较大，会大大降低计算精度，这也是为什么连续方程高残差的原因。 网格线与流动是否一致也会影响计算精度。 4）网格单元形状的影响： 非结构网格比结构网格的截断误差大，因此，为提高计算精度计，请大家尽量使用结构网格，对于复杂几何，在近壁这些对流动影响较大的地方尽量使用结构网格，在其他次要区域使用非结构网格。 2. 不要使用那些书上写的y+与yp的计算公式，那个公式一般只能提供数量级上的参考。推荐大家使用NASA的粘性网格间距计算器，设定你想要的y+值，它就能给你计算出第一层网格高度，与计算结果的y+很接近。 3. 关于边界层网格高度与长度的比例，有本CFD书上说，大概在1/sqrt(Re)就可以；另外，也有这种说法，在做粘性计算时，这个比值可以在100-1000之间，无粘有激波计算时，这个比值要相应小点儿，在10-100之间，因为要考虑激波捕捉精度问题 对于一般的cfd程序，结构化网格要求正交性和光滑性（就是你说的最大最小比率相差不大，想不出一个名次就用这个了）要比较好但是对于fluent这样基于非结构网格的，尤其是其中程序中加入了很多加快收敛速度的方法的软件，后者要求就不要太高。我觉得真正需要考虑网格影响的，一般应该在结构网格上才需要，基于非结构网格的有限体积法，计算通量的时候存在相邻节点的通量计算本身就可能存在计算误差，所以精度不会高到那儿，。顺便说一下，对于fluent，顶多二阶格式就够了，而且绰绰有余，一般我都用一阶因为完全没有必要，其在计算中的误差远远达不到二阶的精度。网格质量本身与具体问题的具体几何特性、流动特性及流场求解算法有关。因此，网格质量最终要由计算结果来评判，但是误差分析以及经验表明，CFD计算对计算网格有一些一般性的要求，例如光滑性、正交性、网格单元的正则性以及在流动变化剧烈的区域分布足够多的网格点等。对于复杂几何外形的网格生成，这些要求往往并不可能同时完全满足。例如，给定边界网格点分布，采用Laplace方程生成的网格是最光滑的，但是最光滑的网格不一定满足物面边界正交性条件，其网格点分布也很有可能不能捕捉流动特征，因此，最光滑的网格不一定是最好的网格。??对计算网格的一个最基本的要求当然是所有网格点的Jacobian必须为正值，即网格体积必须为正，其他一些最常用的网格质量度量参数包括扭角（skew angle）、纵横比（aspect ratio、Laplacian、以及弧长（arc length）等。通过计算、检查这些参数，可以定性的甚至从某种程度上定量的对网格质量进行评判。Parmley等给出了更多的基于网格元素和网格节点的网格质量度量参数。有限元素法关于插值逼近误差估计的理论，实际上也对网格单元的品质给出了基本的规定：即每个单元的内切球半径与外切球半径之比，应该是一个适当的，与网格疏密无关的常数