[论文解读] Structure preserving nodal continuous Finite Elements via Global Flux quadrature

Wasilij Barsukow, Mario Ricchiuto+1 more|arXiv (Cornell University)|Jul 15, 2024

Advanced Numerical Methods in Computational Mathematics被引用 1

一句话总结

本论文提出一种新型的保持结构的有限元方法，用于双曲型PDE，在结构化笛卡尔网格上采用全局通量积分，实现相容性稳定化，且在散度为零的状态下该稳定化项消失。该方法证明了其具有保持涡度的特性，对静止状态具有超收敛性，并能准确刻画离散核，克服了经典SUPG与OSS方法在保持约束（如∇·v = 0）方面的局限性。

ABSTRACT

Numerical methods for hyperbolic PDEs require stabilization. For linear acoustics, divergence-free vector fields should remain stationary, but classical Finite Difference methods add incompatible diffusion that dramatically restricts the set of discrete stationary states of the numerical method. Compatible diffusion should vanish on stationary states, e.g. should be a gradient of the divergence. Some Finite Element methods allow to naturally embed this grad-div structure, e.g. the SUPG method or OSS. We prove here that the particular discretization associated to them still fails to be constraint preserving. We then introduce a new framework on Cartesian grids based on surface (volume in 3D) integrated operators inspired by Global Flux quadrature and related to mimetic approaches. We are able to construct constraint-compatible stabilization operators (e.g. of SUPG-type) and show that the resulting methods are vorticity-preserving. We show that the Global Flux approach is even super-convergent on stationary states, we characterize the kernels of the discrete operators and we provide projections onto them.

研究动机与目标

• 解决经典稳定化有限元方法（如SUPG与OSS）在保持离散约束（如∇·v = 0）方面的失败问题。
• 开发一种在线性声学及相关双曲系统中保持静止性与对合性（involution）的方法。
• 构建与梯度-散度结构相容的稳定化算子，确保在散度为零状态下数值耗散消失。
• 在笛卡尔网格上，通过表面积分算子建立约束相容稳定化的系统性框架。
• 刻画所得算子的离散核，并证明在静止解下具有超收敛性。

提出的方法

• 利用全局通量积分定义表面积分（三维中为体积分）算子，模拟拟麦考利有限差分的性质。
• 采用基于高斯–勒让德-洛巴托积分点与权重的节点连续有限元格式，构建一致且稳定的离散算子。
• 提出一种基于散度梯度（∇(∇·v)）的稳定化机制，确保当∇·v = 0时耗散消失。
• 推导离散算子Dx、Dx^x与Zx为表示离散导数与混合算子的矩阵，对边界与界面条件进行细致处理。
• 通过矩阵分解与质量矩阵求逆，定义稳定化算子Zx = Dx^x − Dx M^{-1} Dx，确保对称性与半正定性。
• 利用矩阵转置可逆性与子式分解进行核分析，刻画离散算子的零空间。

实验结果

研究问题

• RQ1 能否构造一种连续有限元方法，使得数值耗散在散度为零的速度场下消失，从而在线性声学中保持静止性？
• RQ2 为何经典SUPG与OSS方法尽管有理论基础，却仍无法保持梯度-散度结构？
• RQ3 如何利用全局通量积分在笛卡尔网格上定义一致、稳定且保持结构的稳定化方法？
• RQ4 稳定化算子的离散核具有何种结构？其与连续核∇·v = 0的关系如何？
• RQ5 所提方法在静止状态下是否具有超收敛性？若成立，其条件为何？

主要发现

• 该方法通过确保速度的离散旋度在时间上保持恒定，实现了涡度保持性，与连续对合关系∂t(∇×v) = 0一致。
• 该方法在静止状态下具有超收敛性，即当解为散度为零时，误差衰减速度超过一般精度阶次。
• 离散算子Dx^x的核为一维，由常数函数1张成，对应于∇·v = 0的连续核。
• 稳定化算子Zx = Dx^x − Dx M^{-1} Dx的核包含仿射函数⟨1, x⟩，其维数随多项式阶次与单元数量增加而增大，表明其结构比经典方法更丰富。
• 矩阵D^T可逆，意味着离散系统适定，且稳定化算子唯一定义。
• 与Zx相关的双线性形式对称且半正定，确保了离散格式的稳定性和物理一致性。