非结构网格上迎风有限元格式及其迎风类差分格式在非结构网格上的推广.pdfVIP

非结构网格上迎风有限元格式 及迎风类差分格式在非结构网格上的推广 蔡庆东1 吴望一2 (1北京应用物理与计算数学研究所计算物理实验室：北京大学力学与工程科学系) 摘要： 本文在非结构网格上给出一种基于修正积分区域的迎风有限元格式，它与～阶迎风差分 格式相当，可应用于构造各种不同的非结构网格上的迎风型数值格式。 1引 言 纵观差分方法的发展，从以Lax—Wendroff格式为代表的中心型格式到基于迎风型 式等，差分法中发展了一大批基于迎风格式的数值方法，它们在计算流体力学领域发挥了 差分方法中那样，构造出一批基于此的迎风型有限元格式，有限元法中从未出现相当于 察其原因，有限元中的迎风是通过修正权函数实现的，在一维情况，可以选择合适的参数， {．．；； ：牌一侧的贡献消去，从而与差分法中的迎风格式等同，但在多维情况，尤其是非结构网格 却很难做到这一点，它总是使用了与中心型格式相同的基架点集。因此，多维非结构网格 上的迎风有限元格式实际上并不具有差分法中迎风格式的特点，也就不能用来构造有限元 中的迎风型高阶格式。 本文提出了基于修正积分区域的迎风有限元格式，通过对积分区域的修正，达到使 用迎风型的基架点集的目的。这种迎风型有限元格式与差分法中的一阶迎风格式一样，可 用来构造高阶的迎风型有限元格式。本文构造了Warming—Beam型有限元格式，二阶迎 风型有限元格式及NND型的有限元格式，依照本文的方法，可以构造出各种非结构网格 H≮瑚雹l。 上的有限元格式。 2一阶导数的单侧迎风型有限元格式 差分法中的一阶迎风格式，实际上就是计算一点的单侧导数。在结构性的差分网格 中，实现这一点是容易的，但在非结构网格上，目前尚无好的方法。为了在非结构网格上 实现单侧导数的计算，我们首先将其转换为求积分的问题，然后构造一种基于修正积分区 域的迎风有限元格式，即积分仅在一点的迎风区域上进行。 我们考虑二维三角形网格，对每个单元的每个节点，引入迎风系数a。‘．它与单元 e，及局部节点编号有关。以图l所示区域为例，不失一般性，假设有四个单元 e，(f=1,2，3，4)，中心点总体编号为l，四个单元的局部节点号按逆时针排序，且总以中 心点为第一点，那么对于正x轴方向的通量，有如下表的迎风系数： 134 P 单元编号 巳 P2 3 已4 局部节点号 l 2 3 l 2 3 l 2 3 1 2 3 1 0 1 0 0 1 0 1 ％‘ Al、／Ael A2x／Ae2 A3x／Ad A4x／A。4 其中A。为单元P，的面积，爿。为单元e，相对于 中心点的迎风区域的面积。我们欲计算定义于节 矗F 点上的函数F的单测导数G=j≥，采用线性插 m 值，插值基函数为Ⅳ，，那么，F=(Ⅳf． G=G．Ⅳ，它的积分形式的有限元方程为(仅在 羹嗣 单元e．上考虑问题) 么 A孙 。。N GNjdr2 l 2l秘艘 (1)2图．一一图3 q