Arbitrarily high-order structure-preserving schemes for the Gross-Pitaevskii equation with angular momentum rotation in three dimensions
arXiv:2004.08733 · doi:10.1016/j.cpc.2020.107767
Abstract
In this paper, we design a novel class of arbitrarily high-order structure-preserving numerical schemes for the time-dependent Gross-Pitaevskii equation with angular momentum rotation in three dimensions. Based on the idea of the scalar auxiliary variable approach which is proposed in the recent papers [J. Comput. Phys., 416 (2018) 353-407 and SIAM Rev., 61(2019) 474-506] for developing energy stable schemes for gradient flow systems, we firstly reformulate the Gross-Pitaevskii equation into an equivalent system with a modified energy conservation law. The reformulated system is then discretized by the Gauss collocation method in time and the standard Fourier pseudo-spectral method in space, respectively. We show that the proposed schemes can preserve the discrete mass and modified energy exactly. Numerical results are addressed to verify the efficiency and high-order accuracy of the proposed schemes.
18pages, 40 figures
References in corpus (12)
- Mathematical theory and numerical methods for Bose-Einstein condensation
- Computational methods for the dynamics of the nonlinear Schrodinger/Gross-Pitaevskii equations
- Numerical Approximations for a three components Cahn-Hilliard phase-field Model based on the Invariant Energy Quadratization method
- Kolmogorov spectrum of superfluid turbulence: numerical analysis of the Gross-Pitaevskii equation with the small scale dissipation
- Derivation of the Gross-Pitaevskii Equation for Rotating Bose Gases
- A Rigorous Derivation of the Gross-Pitaevskii Energy Functional for a Two-Dimensional Bose Gas
- Hamiltonian Boundary Value Methods (Energy Conserving Discrete Line Integral Methods)
- An efficient and spectrally accurate numerical method for computing dynamics of rotating Bose-Einstein condensates
- General local energy-preserving integrators for solving multi-symplectic Hamiltonian PDEs
- Global well-posedness for the Gross-Pitaevskii equation with an angular momentum rotational term
- Mass and energy conservative high order diagonally implicit Runge--Kutta schemes for nonlinear Schrödinger equation in one and two dimensions
- High-order mass- and energy-conserving SAV-Gauss collocation finite element methods for the nonlinear Schrödinger equation