Abstract:As a basic principle, benefits of adaptive discretisations are an improved balance between required accuracy and efficiency as well as an enhancement of the reliability of numerical computations. In this work, the capacity of locally adaptive space and time discretisations for the numerical solution of low-dimensional nonlinear Schrödinger equations is investigated. The considered model equation is related to the time-dependent Gross-Pitaevskii equation arising in the description of Bose-Einstein condensates in dilute gases. The performance of the Fourier-pseudo spectral method constrained to uniform meshes versus the locally adaptive finite element method and of higher-order exponential operator splitting methods with variable time stepsizes is studied. Numerical experiments confirm that a local time stepsize control based on a posteriori local error estimators or embedded splitting pairs, respectively, is effective in different situations with an enhancement either in efficiency or reliability. As expected, adaptive time-splitting schemes combined with fast Fourier transform techniques are favourable regarding accuracy and efficiency when applied to Gross-Pitaevskii equations with a defocusing nonlinearity and a mildly varying regular solution. However, the numerical solution of nonlinear Schrödinger equations in the semi-classical regime becomes a demanding task. Due to the highly oscillatory and nonlinear nature of the problem, the spatial mesh size and the time increments need to be of the size of the decisive parameter [Formula: see text], especially when it is desired to capture correctly the quantitative behaviour of the wave function itself. The required high resolution in space constricts the feasibility of numerical computations for both, the Fourier pseudo-spectral and the finite element method. Nevertheless, for smaller parameter values locally adaptive time discretisations facilitate to determine the time stepsizes sufficiently small in order that the numerical approximation captures correctly the behaviour of the analytical solution. Further illustrations for Gross-Pitaevskii equations with a focusing nonlinearity or a sharp Gaussian as initial condition, respectively, complement the numerical study.

A Frequency-Dependent P-Adaptive Technique for Spectral Methods

Adaptive Hermite Spectral Methods in Unbounded Domains

SYMPLECTIC PSEUDOSPECTRAL TIME-DOMAIN SCHEME FOR SOLVING TIME-DEPENDENT SCHRODINGER EQUATION

EFFICIENT SCALING AND MOVING TECHNIQUES FOR SPECTRAL METHODS IN UNBOUNDED DOMAINS

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

A Novel Pseudospectral Time Domain Method for Complex EM Simulations

An Adaptive Spectral Method for Oscillatory Second-Order Linear ODEs with Frequency-Independent Cost

Development of Time-Space Adaptive Smoothed Particle Hydrodynamics Method with Runge-Kutta Chebyshev Scheme

Adaptive mesh strategies for the spectral element method

Spectral measures and iterative bounds for effective diffusivity of steady and space-time periodic flows

Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations

A numerical study of adaptive space and time discretisations for Gross-Pitaevskii equations

A novel spectral method for the subdiffusion equation

Solving time-dependent PDEs with the ultraspherical spectral method

Spectrally Adapted Physics-Informed Neural Networks for Solving Unbounded Domain Problems

Spectral Collocation Technique for Absorbing Boundary Conditions with Increasingly High Order Approximation

Convergence study of the h-adaptive PUM and the hp-adaptive FEM applied to eigenvalue problems in quantum mechanics

An adaptive finite element method for high-frequency scattering problems with smoothly varying coefficients

A Highly Efficient Numerical Method for the Time-Fractional Diffusion Equation on Unbounded Domains

Spectral Methods for Multiscale Stochastic Differential Equations