Numerical Integrators for Continuous Disordered Nonlinear Schrödinger Equation

In this paper, we consider the numerical solution of the continuous disordered nonlinear Schrödinger equation, which contains a spatial random potential. We address the finite time accuracy order reduction issue of the usual numerical integrators on this problem, which is due to the presence of the random/rough potential. By using the recently proposed low-regularity integrator (LRI) from (33, SIAM J Numer Anal, 2019), we show how to integrate the potential term by losing two spatial derivatives. Convergence analysis is done to show that LRI has the second order accuracy in L2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^2$$\end{document}-norm for potentials in H2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H^2$$\end{document}. Numerical experiments are done to verify this theoretical result. More numerical results are presented to investigate the accuracy of LRI compared with classical methods under rougher random potentials from applications.