Comparison of Split-Step and Hamiltonian Integration Methods for Simulation of the Nonlinear Schrödinger Equation
arXiv:2008.03938 · doi:10.1016/j.jcp.2020.110061
Abstract
We provide a systematic comparison of two numerical methods to solve the widely used nonlinear Schrödinger equation. The first one is the standard second order split-step (SS2) method based on operator splitting approach. The second one is the Hamiltonian integration method (HIM). It allows the exact conservation of the Hamiltonian at the cost of requiring the implicit time stepping. We found that numerical error for HIM method is systematically smaller than the SS2 solution for the same time step. At the same time, one can take orders of magnitude larger time steps in HIM compared with SS2 still ensuring numerical stability. In contrast, SS2 time step is limited by the numerical stability threshold.
References in corpus (7)
- Numerical simulation of surface waves instability on a discrete grid
- Strong Collapse Turbulence in Quintic Nonlinear Schrödinger Equation
- Fully parallel algorithm for simulating dispersion-managed wavelength-division-multiplexed optical fiber systems
- Cascades in nonlocal turbulence
- Proof of concept implementation of the massively parallel algorithm for simulation of dispersion-managed WDM optical fiber systems
- Langmuir wave filamentation in the kinetic regime. I. Filamentation instability of Bernstein-Greene-Kruskal modes in multidimensional Vlasov simulations
- Intermittency in generalized NLS equation with focusing six-wave interactions