基于流体特征的任意多边形非结构网格方法研讨.pdfVIP

基于流体特征的任意多边形非结构网格方法 研究幸 王瑞利，林忠，葛全文，付汉清，姚彦忠 (北京应用物理与计算数学研究所计算物理实验室北京) 摘要：本文基于非结构网格的思想，发展了任意多边形网格的数值模拟技术。采用结构网格与非结构 网格(非结构三角形、非结构四边形和非结构六边形)联合技术和拉氏有限体积方法的强适应性，以及 非结构网格邻域(拓扑)可随时空灵活变化，研制出基于任意形状网格的拉氏有限体积方法的数值模拟 技术。此方法能较好的适应高温高压多介质大变形流体的计算。几个数值例子也显示了基于流体特征的 网格邻域变化方法对高温高压多介质大变形问题，有很强的数值模拟能力。 关键词：流体特征非结构网格邻域变化有限体积法 1引言 高温高压多介质流体问题的数值模拟也是当今计算流体力学者关心的热门课题。此类问题常需要 数值模拟能精确描述多介质的界面，拉氏方法是此类问题常用的方法。拉氏方法的特点是采用贴体网 格，可以精确描述物质的界面，但网格跟随流体的局部速度的变化而移动，对复杂流体大变形问题常 常由于网格的变形影响计算的精度，甚至计算的成败。为了适应复杂流体大变形的数值模拟，非结构 网格引起国内外计算科学家和学者的喜爱，发展了很多非结构网格的数值模拟方法，但大部分非结构 网格的研究都集中在三角形网格上，且主要应用于初始网格的生成。对于任意多边形非结构网格的研 究工作极少，在计算当中随流体的变化而网格自然跟随流体变化的数值模拟技术几乎没有。本文基于 非结构网格的思想，发展了任意多边形网格的数值模拟技术。采用结构网格与非结构网格(非结构三 角形、非结构四边形和非结构六边形)联合技术和拉氏有限体积方法的强适应性，以及非结构网格邻 域(拓扑)可随时空灵活变化，研制出基于任意形状网格的拉氏有限体积方法的数值模拟技术。此方 法能较好的适应高温高压多介质大变形流体的计算，具有很强的数值模拟能力。 2结构网格和非结构网格 2．1结构网格与非结构嘲格的描述 根据网格点连接方式的不同，网格主要分为两大类。如果内点和它相邻点的连接方式与它的位置 无关，就称为结构网格；如果点点的连接方式不同，就称为非结构网格。对于结构网格，网格点数据 到矩阵元素的存储中就隐含了网格的连接关系。而非结构网格的连接方式必须要使用恰当的数据结构 显式的描述。因此，在非结构网格上构造算法远比在结构网格上复杂得多。但是，非结构网格具有非 常好的几何适应性，当处理复杂几何区域或者要求网格适应流场的复杂特性时，这是至关重要的。 例如，微分算子三u：—Ou_=(x一,y)+—Ov_=(x一,y)。 积 砂 ’中物院基金200320603)和计算物理重点实验室基金资助项目 103 !!兰：!二丝二!!!。兰!!±!二兰：!二! 在结构网格下，离散格式为：(tu)u= Xi+1一Xi—l Yj+1一Yj一1 其中(∥)表示网格编号(一致的邻域数及连接方式，网格间的相邻关系由编号本身即可确定)， 即结构网格用(彬)就可以描述。 在非结构网格下，无论采用什么计算格式，必须建立点与网格之间的邻域关系。通常需建立两个 基本的关系表，这里定义为点格表和格点表。 格点表：如图la，将网格f的全部节点口i，j=l，2，．．．，t，按逆时针方向排列起来，并令 域点数。 点格表：如图lb，付凸啊-114一点口的全部网格tk=-I，2，3，．．．，，％，按逆时针方向排列起来，并令 的邻域网格数。 a 针对上面的微分算子，在非结构网格下，有限体积格式，Lu在网格i上积分 (加)；4≈L三u拟=E(罢+考]姒=t“方一诎≈否l(‰吨一、瓴) (三u)。2去吝(毪咄一、吨) 其中i表示网格编号，二维空间网格一维编号：1，2，…，MC(邻域数不一致，连接方式不一