hyacinth 发表于 2006-5-22 16:59

[转帖]网格划分的几种基本处理方法

<P><STRONG><FONT color=#0000ff>贴体坐标法:</FONT></STRONG></P>
<P>      贴体坐标是利用曲线坐标,并使其坐标线与燃烧室外形或复杂计算区域边界重合,这样所有边界点能够用网格点来表示,不需要任何插值。一旦贴体坐标生成通过变换,偏微分方程求解可以不在任意形状的物理平面上,而在矩形或矩形的组合(空间问题求解域为长方体或它们的组合)转换平面上进行。这样计算与燃烧室外形无关,也与在物理平面上网格间隔无关。</P>
<P>      而是把边界条件复杂的问题转换成一个边界条件简单的问题;这样不仅可避免因燃烧室外形与坐标网格线不一致带来计算误差,而且还可节省计算时间和内存,使流场计算较准确,同时方便求解,较好地解决了复杂形状流动区域的计算,在工程上比较广泛应用。</P>
<P><FONT color=#0000ff><STRONG>区域法:</STRONG></FONT></P>
<P>      虽然贴体坐标系可以使坐标线与燃烧室外形相重合,从而解决复杂流动区域计算问题。但有时实际流场是一个复杂的多通道区域,很难用一种网格来模拟,生成单域贴体网格,即使生成了也不能保证网格质量,影响流场数值求解的效果。因此,目前常采用区域法或分区网格,其基本思想是,根据外形特点把复杂的物理域或复杂拓扑结构的网格,分成若干个区域,分别对每个子区域生成拓扑结构简单的网格。由这些子区域组合而成的网格,或结构块网格。对区域进行分区时,若相邻两个子域分离边界是协调对接,称为对接网格;若相邻两子域有相互重叠部分,则此分区网格称为重叠网格。根据实际数值模拟计算的需要,把整个区域(燃烧室)分成几个不同的子区域,并分别生成网格。这样不仅可提高计算精度,而且还可节省计算机内存,提高收敛精度。但是计算时,必须考虑各区域连接边界处耦合以及变量信息及时、准确地传递问题。处理各个区域连接有多种方法,其中一个办法是在求解各变量时各区域可以单独求解若干次而对压力校正方程.设压力校正值在最初迭代时为零,为了保证流量连续各个区域应同时求解,然后对各个速度和压力进行校正。或者采用在两个区域交界处有一个重叠区,两个区域都对重叠区进行计算,重叠区一边区域内的值,要供重叠区另一边区域求解时用。或通过在重叠内建立两个区域坐标对应关系,实现数据在重叠区内及时传递。如果两个区采用网格疏密分布不相同,要求重叠区二边流量相等。区域法能合理解决网格生成问题,已被大量用来计算复杂形状区域流动。</P>
<P><STRONG><FONT color=#0000ff>区域分解法:</FONT></STRONG></P>
<P>       对于复杂几何形状的实际燃烧装置,为了保证数值求解流场质量,目前常采用区域分解法。该法基本要点是:根据燃烧室形状特点和流场计算需要,把计算区域分成一个主区域和若干个子区域,对各个区域(块)分别建立网格,并对各个区域分别进行数值求解。区域分解原则是尽量使每个子区域边界简便以便于网格建立,各个子区域大小也尽可能相同,使计算负载平衡有利于平行计算。各区域的网格间距数学模型以及计算方法都可以不同,通常在变量变化梯度大的区域,可以布置较细网格,并采用高阶紊流模型和描述复杂反应的紊流燃烧模型,以便更合理模拟实际流场。对于变量变化不太大区域,可采用较疏的网格和较简单的数学模型,这样可节省计算时间。各子区域的解在相邻子区域边界处通过耦合条件来实现光滑,相邻子区域连接重叠网格或对接网格来实现,在各子区域交界处通过插值法提供各子域求解变量的信息传递,满足各子域流场计算要求通量和动量守恒条件以便实现在交界面处各子域流场解的匹配和耦合,从而取得全流场解。</P>
<P><STRONG><FONT color=#0000ff>非结构网格法:</FONT></STRONG></P>
<P>       上述各方法所生成的网格均属于结构化网格,其共同特点是网格中各节点排列有序,每个节点与邻点之间关系是固定的,在计算区域内网格线和平面保持连续。特别是其中分区结构网格生成方法已积累了较多经验,计算技术也较成熟,目前被广泛用来构造复杂外形区域内网格。但是,若复杂外形稍有改变,则将需要重新划分区域和构造网格,耗费较多人力和时间。为此,近年来又发展了另一类网格——非结构网格。此类网格的基本特点是:任何空间区域都被以四面体为单元的网格所划分,网格节点不受结构性质限制,能较好地处理边界,每个节点的邻点个数也可不固定,因此易于控制网格单元的大小、形状及网格的位置。与结构网格相比,此类网格具有更大灵活性和对复杂外形适应性。在20世纪80年代末和90年代初,非结构网格得到了迅速发展。生成非结构网格方法主要有三角化方法和推进阵面法两种。虽然非结构网格容易适合复杂外形,但与结构网格相比还存在一些缺点:(1)需要较大内存记忆单元节点之间关联信息;(2)需要更多CPU时间,这不仅是因为网格结构不规则而增加寻址时间,而且因网格不具备方向性,导致计算工作量增大;(3)结构网格中成熟流场计算方法不能简单地用于非结构网格,离散时所形成代数方程求解过程收敛性差。</P>
<P><STRONG><FONT color=#0000ff>多重网格法:</FONT></STRONG></P>
<P>       多重网格法是一种具有快速收敛特点的计算技术。该法在求解偏微分方程时用一系列逐步加密或减疏的网格去离散求解区域,不同粗细网格可以消除不同波长的误差,从而加快收敛。该法的基本思想是在粗细不同的网格上用迭代法求解差分方程,在每层网格上求出的解包括两部分:一是上一层的解在该层网格上的插值,另一是该层网格消除的误差。将该层所消的误差(上一层未能消除)插值到上一层网格上,作为对上一层原有解的修正,从而得出差分方程在该层上的解。</P>
<P>       多重网格法可把现有计算程序的计算速度提高l~2个数量级,因此近年来得到迅速发展,已推广应用于可压缩反应流和非结构网格、贴体网格系统等。随着计算问题越来越复杂,需求解方程数目越来越多,形式也越来越复杂,为了提高计算速度把多重网格应用于各种工程技术问题显得格外重要。</P>
页: [1]
查看完整版本: [转帖]网格划分的几种基本处理方法