摘要
This paper proposes and analyzes an efficient finite difference scheme for the two-dimensional nonlinear Schr?dinger(NLS) equation involving fractional Laplacian. The scheme is based on a weighted and shifted Grünwald-Letnikov difference(WSGD) operator for the spatial fractional Laplacian. We prove that the proposed method preserves the mass and energy conservation laws in semi-discrete formulations. By introducing the differentiation matrices, the semi-discrete fractional nonlinear Schr?dinger(FNLS) equation can be rewritten as a system of nonlinear ordinary differential equations(ODEs) in matrix formulations. Two kinds of time discretization methods are proposed for the semi-discrete formulation. One is based on the Crank-Nicolson(CN) method which can be proved to preserve the fully discrete mass and energy conservation. The other one is the compact implicit integration factor(c IIF) method which demands much less computational effort. It can be shown that the cIIF scheme can approximate CN scheme with the error O(τ~2). Finally numerical results are presented to demonstrate the method’s conservation, accuracy, efficiency and the capability of capturing blow-up.
This paper proposes and analyzes an efficient finite difference scheme for the two-dimensional nonlinear Schr?dinger(NLS) equation involving fractional Laplacian. The scheme is based on a weighted and shifted Grünwald-Letnikov difference(WSGD) operator for the spatial fractional Laplacian. We prove that the proposed method preserves the mass and energy conservation laws in semi-discrete formulations. By introducing the differentiation matrices, the semi-discrete fractional nonlinear Schr?dinger(FNLS) equation can be rewritten as a system of nonlinear ordinary differential equations(ODEs) in matrix formulations. Two kinds of time discretization methods are proposed for the semi-discrete formulation. One is based on the Crank-Nicolson(CN) method which can be proved to preserve the fully discrete mass and energy conservation. The other one is the compact implicit integration factor(c IIF) method which demands much less computational effort. It can be shown that the cIIF scheme can approximate CN scheme with the error O(τ~2). Finally numerical results are presented to demonstrate the method’s conservation, accuracy, efficiency and the capability of capturing blow-up.
基金
supported by National Natural Science Foundation of China(Grant Nos.61573008 and 61703290)
Laboratory of Computational Physics(Grant No.6142A0502020717)
National Science Foundation of USA(Grant No.DMS-1620108)
