高波数且具角奇性的非线性 Helmholtz（NLH）方程的适应有限元与连续内罚有限元（CIPFEM）方法：后验误差估计、自适应收敛性及光学双稳态模拟

《European Journal of Applied Mathematics》：Adaptive FEM and CIPFEM for a nonlinear Helmholtz equation with high wave number

　　

具有高频率和角奇异性（corner singularities）的非线性 Helmholtz（NLH）方程被线性有限元方法（FEM）和连续内罚有限元方法（CIPFEM）离散化。在推导出若干波数显式（wave-number-explicit）稳定性估计和 NLH 问题的奇异性分解（singularity decomposition）之后，针对形状正则网格（shape regular meshes），包括局部加密网格（locally refined meshes）情形，建立了 (CIP)FEM 的先验（a priori）稳定性和误差估计。随后，针对 NLH 问题的 (CIP)FE 解，利用一种新的残量型误差估计子（residual-type error estimator）推导了后验（a posteriori）上界和下界，该估计子与标准估计子等价。这些后验估计证实了一个重要事实，该事实对 NLH 问题同样成立，即残量型估计子在预渐近区域（preasymptotic regime）严重低估 FE 解的误差；这一现象最初由 Babu?ka et al. [Int J Numer Methods Eng 40 (1997)] 针对一维线性问题观察到。基于该后验误差估计子，当初始网格尺寸位于预渐近区域时，首次针对 NLH 问题证明了所得自适应有限元算法的收敛性和准最优性（quasi-optimality）。最后，给出了数值算例以验证理论发现，并表明采用适当罚参数（penalty parameters）的 CIP 技术能够有效减少污染误差（pollution errors）。特别地，具有 Gaussian 入射波的光学双稳态（optical bistability）非线性现象被自适应 CIPFEM 成功模拟。

本文研究背景为高波数非线性 Helmholtz（nonlinear Helmholtz, NLH）方程在非线性光学介质波散射中的应用，模型包含 Kerr 介质（Kerr medium）、阻抗边界条件（impedance boundary condition）和齐次 Dirichlet 边界条件（homogeneous Dirichlet boundary condition）。高波数下有限元方法（finite element method, FEM）存在污染效应（pollution effect）；精确解还受重入角奇异性、非线性项 |u|²u 引起的局部快速变化以及界面 ?Ω₀ 不连续性影响，导致均匀加密网格效率低。线性 Helmholtz 方程已有自适应有限元研究，但 NLH 问题的非线性项不满足既有工作中常假设的单调性或线性增长条件，因此需要新的理论工具。研究人员开展的研究是构造并分析自适应 FEM 与连续内罚有限元方法（continuous interior penalty finite element method, CIPFEM），即 A(CIP)FEM，建立波数显式先验估计、后验误差估计以及自适应算法的收敛性和准最优性，并通过数值算例验证。

关键技术方法概括如下：研究人员采用线性 FEM 和 CIPFEM 离散 NLH 问题，引入辅助线性化问题、Newton 迭代和椭圆投影（elliptic projection）论证；建立波数显式稳定性估计与解奇异性分解；构造修改的残量型误差估计子并证明其与标准估计子等价；使用 D?rfler 标记和最新顶点二分法生成局部加密网格；在预渐近区域建立后验上、下界、估计子等价关系、收缩性和准最优分析；数值上以 Gaussian 入射波模拟光学双稳态。无样本队列来源。

研究结果如下。1. 引言：研究人员指出 NLH 模型刻画高频率、大入射波能量下非线性光学介质中的波散射；高波数 FEM 受污染效应制约，低正则和高振荡使均匀加密效率低下，因此需要自适应与 CIP 技术。2. 预备知识：研究人员给出几何假设、能量空间和变分形式，为后续分析奠定基础。2.1 辅助线性化问题：通过带 ψ 的辅助线性问题，得到波数显式稳定性估计以及解分解为正则部分和角奇异函数部分，并给出相应范数估计。2.2 NLH 问题估计：通过 Newton 迭代证明在低能量条件下 NLH 问题存在唯一解，并满足能量范数和 L^∞(Ω₀) 稳定性估计及奇异性分解。2.3 FEM 和 CIPFEM 的先验估计：在形状正则网格和局部加密网格上，研究人员建立离散问题的稳定性、内点 L^∞ 估计以及先验预渐近误差估计，条件涉及 ε|ln h_??,min|M(f,g)² 和 k³h_??^1+α 足够小。3. 后验误差估计：研究人员建立 (CIP)FEM 与椭圆投影的后验误差分析框架。3.1 误差估计子：定义单元残量、边跳残量和振荡项，并引入修改的残量型估计子 η_T 及辅助估计子 η?_T。3.2 上界和下界：证明后验上界，如 ‖?(u-u_h)‖₀?(1+k²h_??)η_??，以及 L² 和边界范数上界；同时证明 η_??² 可由 ‖u-P_hu‖₁² 和振荡项控制，并给出局部下界。研究人员确认残量型估计子严重低估预渐近区域 FE 解误差，而 η_?? 能较好拟合 P_hu 的误差。3.3 两个误差估计子之间的等价关系：证明 η_?? 与 η?_?? 在整体网格及子网格上的等价关系，并导出 D?rfler 标记策略之间的等价性。4. AFEM 和 ACIPFEM 的收敛性与准最优性：研究人员利用椭圆投影论证将 CIPFE 解分析转化为椭圆投影分析。4.1 自适应算法：算法包括求解、计算 η_T、D?rfler 标记、最新顶点二分加密和循环迭代。4.2 收敛性：通过引入辅助等价估计子和估计子缩减性质，证明当初始网格尺寸处于预渐近区域时 A(CIP)FEM 的收缩性，即误差以 C_convⁿ 衰减。4.3 准最优性：在适当标记参数和充分小条件下，证明误差被 (1+k²h_n)kM?(f,g)(#??_n-#??₀)^-1/2 控制，表明所得自适应算法具有准最优复杂度。5. 数值结果：数值验证理论。5.1 精度：在含重入角、Kerr 介质和局部高能区的算例中，自适应网格在奇点和高能区加密；AFEM 误差先以 O(N^-1) 快速下降，后接近 O(N^-1/2)；ACIPFEM 误差更小，η_?? 与 P_hu 误差接近，效应指数显示 FE 解误差被低估而 CIPFE 解不明显低估；小波数算例验证奇异性插值误差率。5.2 光学双稳态：在纯阻抗边界、双六边形区域、Gaussian 入射波和 ε=10^-12 等设置下，ACIPFEM 模拟出散射场能量随入射强度的三支曲线；在 I=349000 附近从下支跳至上支，在 I=130000 附近回落，并计算中间支点 B；表明 1300000 和界面 ?Ω₀ 附近加密，AFEM 估计子初期下降快且最终达到 O(N^-1/2)。

讨论与结论可总结并翻译如下：研究人员首次为 NLH 问题发展 A(CIP)FEM，并系统分析其收敛性和准最优性；后验估计证实残量型估计子在预渐近区域严重低估 FE 解误差，这一事实对 NLH 问题同样成立；基于新估计子的自适应算法在初始网格尺寸位于预渐近区域时具有收敛性和准最优性；数值结果表明，采用适当罚参数的 CIP 技术能有效减少污染误差，且自适应 CIPFEM 成功模拟了 Gaussian 入射波下的光学双稳态非线性现象。研究结论为：高波数 NLH 方程可通过 A(CIP)FEM 在预渐近区域获得可靠的自适应误差控制和准最优收敛，CIPFEM 在减少污染误差与模拟非线性光学现象方面具有优势。