二阶半线性椭圆问题的无稳定化 todo serendipity 虚拟元方法

《Communications in Nonlinear Science and Numerical Simulation》：Stabilization-free serendipity virtual element method for second-order semilinear elliptic problems

　　

本研究提出并分析了一种用于求解二维二阶半线性椭圆问题的无稳定化serendipity虚拟元方法（Stabilization-Free Serendipity Virtual Element Method, S-VEM）。该方法的核心特征在于定义了新型多项式投影算子，这些投影使得数值格式能够保持问题结构并实现自由度的有效约减。同时，研究人员提出了一种针对非线性载荷项的近似方法，该近似可利用离散解的自由度以较低计算代价进行求值，从而提升了方法的计算效率。有限维格式的构造及其实现细节被详细阐述，并在H¹-半范数下严格推导了最优误差估计。此外，数值结果表明所提方法在L²-范数和H¹-半范数下均达到了最优收敛率。

**研究背景与动机**

虚拟元方法（Virtual Element Method, VEM）作为经典有限元方法（Finite Element Method, FEM）的扩展形式，最初由文献[1]提出。其最显著的特征在于放宽了对单元形状和凸性的限制，单元形状可为任意凸或凹多边形或多面体，这一优势使其在网格加密和处理悬挂节点等场景中具有独特价值。VEM通过引入非多项式函数作为基函数来解决不规则几何单元上形函数构造的困难，并通过从VEM空间到多项式空间的投影算子来解决刚度矩阵和载荷项在单元上的可计算性问题。

然而，VEM在内部自由度方面比传统有限元具有更高的计算代价。为此，serendipity虚拟元方法（Serendipity VEM, S-VEM）被引入，该方法在保持VEM收敛率和稳定性的同时削减了内部自由度数量，这一特性对三维VEM尤为重要。半线性椭圆问题因广泛应用于相变建模、非线性扩散、燃烧理论以及流体与固体力学等领域而受到广泛关注。尽管有限元方法在这一领域已有丰富成果，但VEM在求解半线性椭圆问题方面的应用仍处于发展阶段。现有VEM处理半线性椭圆问题时，通常依赖标准稳定化机制离散扩散算子，而稳定化项的各向同性特征可能在某些情况下产生不利影响，且会使自适应程序的设计复杂化。基于此，研究人员采用无稳定化虚拟元方法求解半线性椭圆问题，以消除显式稳定化项，同时结合serendipity构造进一步降低计算成本。

**技术方法**

研究人员为开展本项研究，主要采用了以下几项关键技术方法：首先，基于文献[5]中引入的serendipity构造，建立serendipity虚拟元空间V^S_h,k，通过引入新定义的投影算子实现内部自由度的系统约减；其次，设计无稳定化离散双线性形式a_h(·,·)，仅通过投影构造离散扩散算子，避免使用任何显式稳定化项；再次，提出基于自由度的非线性载荷项近似方法，利用投影算子Π^?,k_E的离散解自由度以低计算代价近似非线性项；最后，在星形多边形网格假设下，应用标准椭圆投影技巧和插值理论推导H¹-半范数最优误差估计。

**研究结果**

**离散格式构造**：基于新型投影算子的定义，研究人员建立了无稳定化的离散双线性形式和载荷项离散形式，提出了半线性椭圆问题的serendipity虚拟元离散格式。通过引入serendipity虚拟元空间，在保证最优收敛性的前提下显著降低了自由度数量，增强了方法的计算效率。

**H¹-半范数误差估计**：在标准正则性假设和网格正则性假设（Assumption 3.1）下，研究人员严格证明了离散格式的适定性，并推导了H¹-半范数下的最优误差估计。该估计表明，当精确解u∈H^k+1(Ω)时，数值解u_h满足|u-u_h|_1,Ω≤Ch^k|u|_k+1,Ω，其中C为与网格尺寸h无关的正常数。

**数值实验**：研究人员设计了数值实验以验证理论分析的正确性并展示S-VEM的优势。在凸和非凸多边形网格上的数值结果显示，所提方法在L²-范数和H¹-半范数下均达到了最优收敛率。通过与标准VEM的自由度数量对比，证实了serendipity构造在降低计算成本方面的有效性。

**讨论与结论**

研究人员在本工作中构建了面向半线性椭圆问题的serendipity虚拟元空间。通过引入新型投影算子，构建了不含外在稳定化项的离散双线性形式，并在H¹-半范数下推导了最优阶误差估计。理论收敛率通过在凸和非凸多边形网格上的数值实验得到验证。当前工作限于二维情形的分析，将其扩展至三维情形是未来研究的方向。该论文发表在《Communications in Nonlinear Science and Numerical Simulation》期刊上。