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
Permanent link to this recordhttp://hdl.handle.net/10754/597405
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.
Large-eddy simulation of flow over a rotating cylinder: the lift crisis at ReD=6×10^4Cheng, W.; Pullin, D. I.; Samtaney, Ravi (Journal of Fluid Mechanics, Cambridge University Press (CUP), 2018-09-19) [Article]We present wall-resolved large-eddy simulation (LES) of flow with free-stream velocity U∞ over a cylinder of diameter D rotating at constant angular velocity Ω , with the focus on the lift crisis, which takes place at relatively high Reynolds number ReD=U∞D/ν , where ν is the kinematic viscosity of the fluid. Two sets of LES are performed within the ( ReD , α )-plane with α=ΩD/(2U∞) the dimensionless cylinder rotation speed. One set, at ReD=5000 , is used as a reference flow and does not exhibit a lift crisis. Our main LES varies α in 0⩽α⩽2.0 at fixed ReD=6×104 . For α in the range α=0.48−0.6 we find a lift crisis. This range is in agreement with experiment although the LES shows a deeper local minimum in the lift coefficient than the measured value. Diagnostics that include instantaneous surface portraits of the surface skin-friction vector field Cf , spanwise-averaged flow-streamline plots, and a statistical analysis of local, near-surface flow reversal show that, on the leeward-bottom cylinder surface, the flow experiences large-scale reorganization as α increases through the lift crisis. At α=0.48 the primary-flow features comprise a shear layer separating from that side of the cylinder that moves with the free stream and a pattern of oscillatory but largely attached flow zones surrounded by scattered patches of local flow separation/reattachment on the lee and underside of the cylinder surface. Large-scale, unsteady vortex shedding is observed. At α=0.6 the flow has transitioned to a more ordered state where the small-scale separation/reattachment cells concentrate into a relatively narrow zone with largely attached flow elsewhere. This induces a low-pressure region which produces a sudden decrease in lift and hence the lift crisis. Through this process, the boundary layer does not show classical turbulence behaviour. As α is further increased at constant ReD , the localized separation zone dissipates with corresponding attached flow on most of the cylinder surface. The lift coefficient then resumes its increasing trend. A logarithmic region is found within the boundary layer at α=1.0 .
Spontaneous Rotational Inversion in PhycomycesGoriely, Alain; Tabor, Michael (Physical Review Letters, American Physical Society (APS), 2011-03-31) [Article]The filamentary fungus Phycomyces blakesleeanus undergoes a series of remarkable transitions during aerial growth. During what is known as the stagea IV growth phase, the fungus extends while rotating in a counterclockwise manner when viewed from above (stagea IVa) and then, while continuing to grow, spontaneously reverses to a clockwise rotation (stagea IVb). This phase lasts for 24-48Ah and is sometimes followed by yet another reversal (stageAIVc) before the overall growth ends. Here, we propose a continuum mechanical model of this entire process using nonlinear, anisotropic, elasticity and show how helical anisotropy associated with the cell wall structure can induce spontaneous rotation and, under appropriate circumstances, the observed reversal of rotational handedness. © 2011 American Physical Society.