Optimal $L^2$ error estimates of mass- and energy-conserved FE schemes for a nonlinear Schrödinger-type system
Zhuoyue Zhang, Wentao Cai
TL;DR
This work develops an implicit Crank–Nicolson finite element scheme for nonlinear Schrödinger–type systems, including Schrödinger–Helmholz and Schrödinger–Poisson models, on convex/polyhedral domains. The method preserves discrete mass and energy, with existence and uniqueness of fully discrete solutions shown via Schaefer's fixed point theorem. It provides optimal $L^2$ error estimates for both models by an error-splitting analysis and Ritz projection techniques, validated by numerical experiments. The results demonstrate stable, structure-preserving discretization with convergence rates matching the theoretical predictions, offering a practical approach for simulations of Schrödinger-type phenomena.
Abstract
In this paper, we present an implicit Crank-Nicolson finite element (FE) scheme for solving a nonlinear Schrödinger-type system, which includes Schrödinger-Helmholz system and Schrödinger-Poisson system. In our numerical scheme, we employ an implicit Crank-Nicolson method for time discretization and a conforming FE method for spatial discretization. The proposed method is proved to be well-posedness and ensures mass and energy conservation at the discrete level. Furthermore, we prove optimal $L^2$ error estimates for the fully discrete solutions. Finally, some numerical examples are provided to verify the convergence rate and conservation properties.
