2 research outputs found
A finite element method with mesh adaptivity for computing vortex states in fast-rotating Bose-Einstein condensates
Numerical computations of stationary states of fast-rotating Bose-Einstein
condensates require high spatial resolution due to the presence of a large
number of quantized vortices. In this paper we propose a low-order finite
element method with mesh adaptivity by metric control, as an alternative
approach to the commonly used high order (finite difference or spectral)
approximation methods. The mesh adaptivity is used with two different numerical
algorithms to compute stationary vortex states: an imaginary time propagation
method and a Sobolev gradient descent method. We first address the basic issue
of the choice of the variable used to compute new metrics for the mesh
adaptivity and show that simultaneously refinement using the real and imaginary
part of the solution is successful. Mesh refinement using only the modulus of
the solution as adaptivity variable fails for complicated test cases. Then we
suggest an optimized algorithm for adapting the mesh during the evolution of
the solution towards the equilibrium state. Considerable computational time
saving is obtained compared to uniform mesh computations. The new method is
applied to compute difficult cases relevant for physical experiments (large
nonlinear interaction constant and high rotation rates).Comment: to appear in J. Computational Physic
A finite-element toolbox for the stationary Gross-Pitaevskii equation with rotation
We present a new numerical system using classical finite elements with mesh
adaptivity for computing stationary solutions of the Gross-Pitaevskii equation.
The programs are written as a toolbox for FreeFem++ (www.freefem.org), a free
finite-element software available for all existing operating systems. This
offers the advantage to hide all technical issues related to the implementation
of the finite element method, allowing to easily implement various numerical
algorithms.Two robust and optimised numerical methods were implemented to
minimize the Gross-Pitaevskii energy: a steepest descent method based on
Sobolev gradients and a minimization algorithm based on the state-of-the-art
optimization library Ipopt. For both methods, mesh adaptivity strategies are
implemented to reduce the computational time and increase the local spatial
accuracy when vortices are present. Different run cases are made available for
2D and 3D configurations of Bose-Einstein condensates in rotation. An optional
graphical user interface is also provided, allowing to easily run predefined
cases or with user-defined parameter files. We also provide several
post-processing tools (like the identification of quantized vortices) that
could help in extracting physical features from the simulations. The toolbox is
extremely versatile and can be easily adapted to deal with different physical
models