Abstract
Currently, the theory of radiation-driven winds of massive stars possess three known solutions for the velocity and density profiles of the stellar winds, namely: the fast, Ω -slow and δ -slow solutions. In order to confirm their stability we use a time-dependent numerical hydrodynamic code called ZEUS-3D, and then we compare their results with the stationary solutions from our numerical hydrodynamic code. ZEUS-3D needs an initial trial solution to start to integrate, for this we use the stationary solution (from our code) or a β-law for the velocity field. In both cases we obtain the same results. Fast and both slow stationary solutions are attained in ZEUS-3D and are all stable. Furthermore, there is a very good agreement with the velocity and density fields from ZEUS-3D and our code, having differences between the terminal velocities lower than 3%. In addition, we found that ZEUS-3D is very sensitive to the boundary conditions (base density and velocity profile), in some cases we obtain kinks in the velocity profiles, similar to the ones obtained by Madura et al. (2007) for stars with high rotation. Such kinks are most likely the result of the wind being mass overloaded, but further investigation is needed to understand its nature better. Currently, we are exploring the effects of small perturbation at the base of the wind in order to study possible transitions or oscillations between δ-slow and fast solutions.
