Construction and Analysis of Structure-Preserving Numerical Algorithm for Two-Dimensional Damped Nonlinear Space Fractional Schrödinger equation

In this paper, we present a novel high-order structure-preserving numerical scheme for solving the damped nonlinear space fractional Schrödinger equation (DNSFSE) in two spatial dimensions. The main idea of constructing new algorithm consists of two parts. Firstly, we introduce an auxiliary exponential variable to transform the original DNSFSE into a modified one. The modified DNSFSE subjects to the conservation of mass and energy, which is crucial to develop structure-preserving numerical schemes. Secondly, we construct a high-order numerical differential formula to approximate the Riesz derivative in space, which contributes to a semi-discrete difference scheme for the modified DNSFSE. Combining the semi-discrete scheme with the variant Crank–Nicolson method in time, we can obtain the fully-discrete difference scheme for solving the modified DNSFSE. The advantage of the proposed scheme is that a fourth-order convergent accuracy can be achieved in space while maintaining the conservation of mass and energy. Subsequently, we conduct a detailed study on the boundedness, uniqueness, and convergence of solution for fully-discrete scheme. Furthermore, an improved efficient iterative algorithm is proposed for the fully-discrete scheme, which has the advantage of maintaining the same convergence order as the original difference scheme. Finally, extensive numerical results are reported to further verify the correctness of theoretical analysis and the effectiveness of the proposed numerical algorithm.