嵌套网格教程.docVIP

CFD-FASTRAN中Chimera帮助 鉴于群里面有一些成员里面有一些嵌套网格初学者，因此应该写一些帮助来提供参考。对于已经掌握嵌套网格的人，可以继续补充和修正一下我的文章。 1 理论基础 使用嵌套方法主要有两个目的：（1）考虑运动体的计算；（2）减少复杂外形网格生成时间。嵌套网格可以相对其它网格独立运动，而不会干扰其它区域的网格，对非定常流场仿真时可以得到较精确的结果。而且嵌套区域可独立生成网格，将复杂的结构网格生成过程简单化。 嵌套方法的基本步骤如下： (1) 在每个区域搜寻确认嵌套区域的边界单元。嵌套边界既可以是区域外边界也可以是内部集合体相邻边界。 (2) 在每个嵌套区域进行挖洞。是否挖洞取决于该区域是否与边界或者其它区域的壁面嵌套。 (3) 对嵌套区域边界单元进行插值。将背景区域的边界单元变量信息插值到嵌套区域的边界单元。 1.1 搜寻确认过程 采用嵌套方法需要寻找到每个网格嵌套区域单元，这个过程一直进行，直到仿真结束。交替数字树（Alternating Digital Tree，简称ADT）算法可以优化这个搜寻过程。ADT以树形结构方式传播的网格单元信息，如网格单元中心坐标，网格边框坐标等。当ADT搜寻到网格与所给的嵌套网格一致时，ADT算法能快速将信息传递到目标单元。当维数为3时，ADT算法与八叉树特别相似。CFD-FASTRAN中的ADT将每个区域的边框坐标信息储存到6维的矩阵中。 1.2 挖洞过程 挖洞过程就是将每个区域的壁面边界部分挖空，用嵌套区域网格部分来代替。在这个过程中同时进行网格信息的获取，并将每个流体变量从一个区域插值到另一个区域中。 要进行单元挖空的区域是主要区域，有壁面边界穿过单元的区域是次要区域。如果在嵌套区域中有穿过壁面边界的单元，则这个单元既属于主要区域又属于次要区域，见图1.1。ADT算法首先寻找并确认这样区域的单元。如果在壁面边界中有网格边穿过，ADT算法就在主要区域中找到那些网格边的两个端点，见图1.2。接下来再确认网格边上的网格节点是在壁面边界内部还是外部。如果节点有一部分在内部，区域单元就标记成为切除单元(cut-cell)，见图1.3。随后把在壁面边界内部的节点也标记出来，见图1.4。挖空单元和切除单元区域合起来就是嵌套网格洞，见图1.5。洞周围的主要区域中的单元就是嵌套边界单元，见图1.6。然后就可以将主要区域单元信息向次要区域单元进行插值。 在嵌套网格洞的外层通常是缓冲层(Buffer Layer)，用来进行信息插值，这样可以避免主要区域的流体变量离壁面太进而无法进行插值。如果用了不止一层的缓冲层，插值就离次要区域的壁面边界相应单元层数开始进行，同时为此采用了边缘层。 1.3 插值过程 以上确认嵌套网格边界和嵌套网格挖洞过程中运用到了ADT算法。在找到嵌套网格边界后，由信息单元节点向边界相邻的单元进行信息传递。这些信息单元节点是次要区域单元中挖洞后新生成的网格点。这个信息传递过程进行三次线性插值，将主要区域中的信息传给嵌套网格边界单元。 在嵌套边界区域内确认孤立网格。孤立网格是那些经过ADT计算后没有信息的嵌套边界单元。这主要是因为在进行挖洞时两个网格体的表面很近导致在进行ATD运算时找不到提供信息的单元。而这些孤立网格单元的信息通常由邻近网格的信息来代替。为了减少或避免孤立网格出现，有相互重叠网格区域网格层数至少为5层。 图1.1 挖洞算法的例子 实线网格区域是主要区域，虚线网格区域是嵌套区域，黑色实体是次要区域的壁面边界。 图1.2 穿过壁面边界网格边的端点状态 图1.3 主要区域中的切除单元 图1.4 次要区域内壁面边界内部的节点 图1.5 在次要区域内部挖空单元或主要区域边上的区域单元生成洞 图1.6 主要区域的嵌套边界单元 2 实际运用 下面介绍一个机翼下挂载弹体的例子，来看看在CFD-FSATRAN中是如何运用嵌套网格的。 在CFD-FSATRAN的程序中会进行，对嵌套网格的搜索、差值和挖洞三个步骤，因此我们要做的就是建立出一套合适的嵌套网格来让CFD-FSATRAN进行以上的三个步骤。 2.1建立一个背景区域 机翼周围生成“H”型结构网格。整个背景计算流场区域为矩形。整个流场一共划分了12个区域，总的网格数量为843500个，并且做了近壁边界处理，同时满足嵌套网格要求。机翼周围网格情况如图2.1所示。如果要进行投弹过程模拟时，在运动的区域内要增加适当的网格数量以便插值时流场信息传递准确。 图2.1 背景区域中机翼根部界面的网格 2.2 建立背景区域内的另一网格 在弹体圆柱型周围采用“O”型网格，弹体尾翼周围采用“H”型网格，弹体两端部分也采用“H”型网格。整个嵌套区域一共划分了26个区域，总的网格数量为244500个。图2.2为嵌套网格