Unconditional optimal-order error estimates of linear relaxation compact difference scheme for the coupled nonlinear Schrödinger system

Abstract

This paper presents a linear, decoupled, mass- and energy-conserving numerical scheme for the multi-dimensional coupled nonlinear Schrödinger (CNLS) system. The scheme combines the fourth-order compact difference approximation in space with the relaxation technique in a time-staggered mesh framework, solving the primal unknowns and introduced auxiliary relaxation variables sequentially with high efficiency and high-order accuracy. We show the unique solvability and discrete conservation laws of the developed scheme. In particular, for the first time, leveraging an auxiliary error equation system combined with the cut-off technique, optimal-order error estimates in the discrete H1-norm for the primal variables at the time nodes, and in the discrete L2-norm for the auxiliary relaxation variables at the intermediate time nodes, are rigorously proved without any coupling mesh conditions, which contribute to the primary theoretical contribution of this paper for multi-dimensional CNLS system. Numerical experiments demonstrate convincingly the strong performance of the proposed scheme in long-term simulations, maintaining both physical invariants and high-order accuracy.

Figures (7)

• Figure 1: Evolution of mass with respect to time.
• Figure 2: Evolutions of mass errors $|\mathcal{M}_u^n-\mathcal{M}_u^0|$ (left) and $|\mathcal{M}_v^n-\mathcal{M}_v^0|$ (right) with respect to time.
• Figure 3: Evolution of mass errors $|\mathcal{R}_u^n-\mathcal{R}_u^0|$ (left) and $|\mathcal{R}_v^n-\mathcal{R}_v^0|$ (right) with respect to time.
• Figure 4: Evolutions of energy (left) and absolute error (right) with respect to time.
• Figure 5: Elastic collisions of two solitons.

Theorems & Definitions (36)

• Lemma 2.1
• Lemma 2.2
• Lemma 2.3
• Lemma 2.4
• proof
• Lemma 2.5
• proof
• Lemma 2.6: ZHOU2025
• Lemma 2.7: ZG2023
• Remark 1