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.
Blind Source Separation Algorithms Using Hyperbolic and Givens Rotations for High-Order QAM ConstellationsShah, Syed Awais Wahab; Abed-Meraim, Karim; Al-Naffouri, Tareq Y. (Institute of Electrical and Electronics Engineers (IEEE), 2017-11-24)This paper addresses the problem of blind demixing of instantaneous mixtures in a multiple-input multiple-output communication system. The main objective is to present efficient blind source separation (BSS) algorithms dedicated to moderate or high-order QAM constellations. Four new iterative batch BSS algorithms are presented dealing with the multimodulus (MM) and alphabet matched (AM) criteria. For the optimization of these cost functions, iterative methods of Givens and hyperbolic rotations are used. A pre-whitening operation is also utilized to reduce the complexity of design problem. It is noticed that the designed algorithms using Givens rotations gives satisfactory performance only for large number of samples. However, for small number of samples, the algorithms designed by combining both Givens and hyperbolic rotations compensate for the ill-whitening that occurs in this case and thus improves the performance. Two algorithms dealing with the MM criterion are presented for moderate order QAM signals such as 16-QAM. The other two dealing with the AM criterion are presented for high-order QAM signals. These methods are finally compared with the state of art batch BSS algorithms in terms of signal-to-interference and noise ratio, symbol error rate and convergence rate. Simulation results show that the proposed methods outperform the contemporary batch BSS algorithms.
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.