By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov te...By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov tempered difference operator to approximate the Riemann–Liouville tempered fractional partial derivative,the fractional central difference operator to discritize the space Riesz fractional partial derivative,and the classical central difference formula to discretize the advection term,a numerical algorithm is constructed for solving the Caputo tempered fractional advection-diffusion equation.The stability and the convergence analysis of the numerical method are given.Numerical experiments show that the numerical method is effective.展开更多
In this paper,we propose numerical methods for the Riesz space fractional advection-dispersion equations with delay(RFADED).We utilize the fractional backward differential formulas method of second order(FBDF2)and wei...In this paper,we propose numerical methods for the Riesz space fractional advection-dispersion equations with delay(RFADED).We utilize the fractional backward differential formulas method of second order(FBDF2)and weighted shifted Grünwald difference(WSGD)operators to approximate the Riesz fractional derivative and present the finite difference method for the RFADED.Firstly,the FBDF2 and the shifted Grünwald methods are introduced.Secondly,based on the FBDF2 method and the WSGD operators,the finite difference method is applied to the problem.We also show that our numerical schemes are conditionally stable and convergent with the accuracy of O(+h2)and O(2+h2)respectively.Thirdly we find the analytical solution for RFDED in terms Mittag-Leffler type functions.Finally,some numerical examples are given to show the efficacy of the numerical methods and the results are found to be in complete agreement with the analytical solution.展开更多
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ü...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.展开更多
文摘By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov tempered difference operator to approximate the Riemann–Liouville tempered fractional partial derivative,the fractional central difference operator to discritize the space Riesz fractional partial derivative,and the classical central difference formula to discretize the advection term,a numerical algorithm is constructed for solving the Caputo tempered fractional advection-diffusion equation.The stability and the convergence analysis of the numerical method are given.Numerical experiments show that the numerical method is effective.
文摘In this paper,we propose numerical methods for the Riesz space fractional advection-dispersion equations with delay(RFADED).We utilize the fractional backward differential formulas method of second order(FBDF2)and weighted shifted Grünwald difference(WSGD)operators to approximate the Riesz fractional derivative and present the finite difference method for the RFADED.Firstly,the FBDF2 and the shifted Grünwald methods are introduced.Secondly,based on the FBDF2 method and the WSGD operators,the finite difference method is applied to the problem.We also show that our numerical schemes are conditionally stable and convergent with the accuracy of O(+h2)and O(2+h2)respectively.Thirdly we find the analytical solution for RFDED in terms Mittag-Leffler type functions.Finally,some numerical examples are given to show the efficacy of the numerical methods and the results are found to be in complete agreement with the analytical solution.
基金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)
文摘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.
