Nonradial and nonpolytropic astrophysical outflows XI. Simulations of the circumstellar environment of RY Tauri

Context. Recent observational evidence has shown that RY Tau may present two different outflow stages, a quiescent and a more active stage. We try to model this phenomenon. Aims. We have performed new 2.5D magnetohydrodynamical simulations of the possible accretion-outflow environment of RY Tau based on analytical solutions with the aim to reduce the relaxation time. Methods. We used the analytical self-similar solution that we used to model the RY Tau microjet as initial conditions. In the closed field line region of the magnetosphere, we reversed the direction of the flow and increased the accretion rate by increasing the density and velocity. We also implemented the heating rate and adjusted it according to the velocity of the flow. The accretion disk was treated as a boundary condition. Results. The simulations show that the stellar jet and the accreting magnetosphere attain a steady state in only a few stellar rotations. This confirms the robustness and stability of self-similar solutions. Additionally, two types of behavior were observed that are similar to the behavior observed in RY Tau. Either the steady stellar outflow and magnetospheric inflow are separated by a low static force-free region or the interaction between the stellar jet and the magnetospheric accretion creates episodic coronal mass ejections that originate from the disk and bounce back onto the star. Conclusions. The ratio of mass-loss rate to mass-accretion rate that coincides with the change in behavior observed in RY Tau lies within the range of ratios that have been measured during the period in which the initial microjet was analyzed.