Section: New Results

Scientific Computing for Linear and Nonlinear Wave Problems

Participant : Xavier Antoine.

We consider the Backward Euler SPectral (BESP) scheme that was proposed for computing the stationary states of Bose-Einstein Condensates (BECs) through the Gross-Pitaevskii equation. We show that the fixed point approach introduced earlier fails to converge for fast rotating BECs. A simple alternative approach based on Krylov subspace solvers with a Laplace or Thomas-Fermi preconditioner is given. Numerical simulations (obtained with the associated freely available Matlab toolbox GPELab) for complex configurations show that the method is accurate, fast and robust for 2D/3D problems and multi-components BECs.

This work was published in the journal “Journal of Computational Physics” [9] .