A Simple and Efficient Numerical Method for Computing the Dynamics of Rotating Bose--Einstein Condensates via Rotating Lagrangian Coordinates
KAUST Grant NumberKUK-I1-007-43
MetadataShow full item record
AbstractWe propose a simple, efficient, and accurate numerical method for simulating the dynamics of rotating Bose-Einstein condensates (BECs) in a rotational frame with or without longrange dipole-dipole interaction (DDI). We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range DDI, state the twodimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential, and review some dynamical laws related to the 2D and 3D GPEs. By introducing a rotating Lagrangian coordinate system, the original GPEs are reformulated to GPEs without the angular momentum rotation, which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable, and very efficient in computation. It is spectral-order accurate in space and second-order accurate in time and conserves the mass on the discrete level. We compare our method with some representative methods in the literature to demonstrate its efficiency and accuracy. In addition, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation, and center of mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without or with the long-range DDI.Copyright © by SIAM.
CitationBao W, Marahrens D, Tang Q, Zhang Y (2013) A Simple and Efficient Numerical Method for Computing the Dynamics of Rotating Bose--Einstein Condensates via Rotating Lagrangian Coordinates. SIAM Journal on Scientific Computing 35: A2671–A2695. Available: http://dx.doi.org/10.1137/130911111.
SponsorsThis research was supported by the Singapore A*STAR SERC "Complex Systems" Research Programme grant 1224504056.This author's work was supported by research award KUK-I1-007-43 from King Abdullah University of Science and Technology (KAUST).This author's work was supported by the University of Missouri Research Board, Simons Foundation award 210138, and NSF grant DMS-1217000.
CollectionsPublications Acknowledging KAUST Support
Showing items related by title, author, creator and subject.
Development of Pre-set Counter-rotating Streamwise Vortices in Wavy ChannelBudiman, A.C.; Mitsudharmadi, Hatsari; Bouremel, Y.; Winoto, S.H.; Low, H.T. (Elsevier BV, 2015-10-23)Development of counter-rotating streamwise vortices in a rectangular channel with one-sided wavy surface has been experimentally quantified using hot-wire anemometry. The wavy surface has fixed amplitude of 3.75 mm. The counter-rotating vortices are pre-set by means of a sawtooth pattern cut at the leading edge of the wavy surface. Variations of the central streamwise velocity Uc with a channel gap H = 35 mm and 50 mm (corresponding to a Reynolds number from 1600 to 4400) change the instability of the flow which can be distinguished from the velocity contours at a certain spanwise plane. The streamwise velocity contours and turbulence intensity for Reynolds number Re = 3100 and H = 35 mm show the disappearance of the mushroom-like vortices prior to turbulence near the second peak of the wavy surface, while for higher Re, this phenomenon occurs earlier. Under certain conditions, for example, for Re = 4400 and H = 50 mm, the splitting of the vortices can also be observed.
Pseudo-spectral method using rotated staggered grid for elastic wave propagation in 3D arbitrary anisotropic mediaZou, Peng; Cheng, Jiubing (Wiley-Blackwell, 2017-05-10)Staggering grid is a very effective way to reduce the Nyquist errors and to suppress the non-causal ringing artefacts in the pseudo-spectral solution of first-order elastic wave equations. However, the straightforward use of a staggered-grid pseudo-spectral method is problematic for simulating wave propagation when the anisotropy level is greater than orthorhombic or when the anisotropic symmetries are not aligned with the computational grids. Inspired by the idea of rotated staggered-grid finite-difference method, we propose a modified pseudo-spectral method for wave propagation in arbitrary anisotropic media. Compared with an existing remedy of staggered-grid pseudo-spectral method based on stiffness matrix decomposition and a possible alternative using the Lebedev grids, the rotated staggered-grid-based pseudo-spectral method possesses the best balance between the mitigation of artefacts and efficiency. A 2D example on a transversely isotropic model with tilted symmetry axis verifies its effectiveness to suppress the ringing artefacts. Two 3D examples of increasing anisotropy levels demonstrate that the rotated staggered-grid-based pseudo-spectral method can successfully simulate complex wavefields in such anisotropic formations.
A Stefan model for mass transfer in a rotating disk reaction vesselBOHUN, C. S. (Cambridge University Press (CUP), 2015-05-04)Copyright © Cambridge University Press 2015. In this paper, we focus on the process of mass transfer in the rotating disk apparatus formulated as a Stefan problem with consideration given to both the hydrodynamics of the process and the specific chemical reactions occurring in the bulk. The wide range in the reaction rates of the underlying chemistry allows for a natural decoupling of the problem into a simplified set of weakly coupled convective-reaction-diffusion equations for the slowly reacting chemical species and a set of algebraic relations for the species that react rapidly. An analysis of the chemical equilibrium conditions identifies an expansion parameter and a reduced model that remains valid for arbitrarily large times. Numerical solutions of the model are compared to an asymptotic analysis revealing three distinct time scales and chemical diffusion boundary layer that lies completely inside the hydrodynamic layer. Formulated as a Stefan problem, the model generalizes the work of Levich (Levich and Spalding (1962) Physicochemical hydrodynamics, vol. 689, Prentice-Hall Englewood Cliffs, NJ) and will help better understand the natural limitations of the rotating disk reaction vessel when consideration is made for the reacting chemical species.