Dynamical simulations of magnetically channelled line-driven stellar winds - II. The effects of field-aligned rotation

Building upon our previous magnetohydrodynamics (MHD) simulation study of magnetic channelling in radiatively driven stellar winds, we examine here the additional dynamical effects of stellar rotation in the (still) 2D axisymmetric case of an aligned dipole surface field. In addition to the magnetic confinement parameter η_* introduced in Paper I, we characterize the stellar rotation in terms of a parameter W ≡ V_rot/V_orb (the ratio of the equatorial surface rotation speed to orbital speed), examining specifically models with moderately strong rotation W = 0.25 and 0.5, and comparing these to analogous non-rotating cases. Defining the associated Alfvén radius R_A ~ η^1/4_* R_* and Kepler corotation radius R_K ~ W^-2/3 R_*, we find rotation effects are weak for models with R_A < R_K, but can be substantial and even dominant for models with R_A >~ R_K. In particular, by extending our simulations to magnetic confinement parameters (up to η_* = 1000) that are well above those (η_* = 10) considered in Paper I, we are able to study cases with R_A >> R_K; we find that these do indeed show clear formation of the rigid body disc predicted in previous analytic models, with however a rather complex, dynamic behaviour characterized by both episodes of downward infall and outward breakout that limit the build-up of disc mass. Overall, the results provide an intriguing glimpse into the complex interplay between rotation and magnetic confinement, and form the basis for a full MHD description of the rigid body discs expected in strongly magnetic Bp stars like σ Ori E.