Volume 633, January 2020
|Number of page(s)||24|
|Section||Numerical methods and codes|
|Published online||23 December 2019|
A 3D short-characteristics method for continuum and line scattering problems in the winds of hot stars
LMU München, Universitätssternwarte, Scheinerstr. 1, 81679 München, Germany
e-mail: firstname.lastname@example.org, email@example.com
2 Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
Accepted: 8 October 2019
Context. Knowledge about hot, massive stars is usually inferred from quantitative spectroscopy. To analyse non-spherical phenomena, the existing 1D codes must be extended to higher dimensions, and corresponding tools need to be developed.
Aims. We present a 3D radiative transfer code that is capable of calculating continuum and line scattering problems in the winds of hot stars. By considering spherically symmetric test models, we discuss potential error sources, and indicate advantages and disadvantages by comparing different solution methods. Further, we analyse the ultra-violet (UV) resonance line formation in the winds of rapidly rotating O stars.
Methods. We consider both a (simplified) continuum model including scattering and thermal sources, and a UV resonance line transition approximated by a two-level-atom. We applied the short-characteristics (SC) method, using linear or monotonic Bézier interpolations, for which monotonicity is of prime importance, to solve the equation of radiative transfer on a non-uniform Cartesian grid. To calculate scattering dominated problems, our solution method is supplemented by an accelerated Λ-iteration scheme using newly developed non-local operators.
Results. For the spherical test models, the mean relative error of the source function is on the 5 − 20% level, depending on the applied interpolation technique and the complexity of the considered model. All calculated line profiles are in excellent agreement with corresponding 1D solutions. Close to the stellar surface, the SC methods generally perform better than a 3D finite-volume-method; however, they display specific problems in searchlight-beam tests at larger distances from the star. The predicted line profiles from fast rotating stars show a distinct behaviour as a function of rotational speed and inclination. This behaviour is tightly coupled to the wind structure and the description of gravity darkening and stellar surface distortion.
Conclusions. Our SC methods are ready to be used for quantitative analyses of UV resonance line profiles. When calculating optically thick continua, both SC methods give reliable results, in contrast to the alternative finite-volume method.
Key words: radiative transfer / methods: numerical / stars: massive / stars: rotation / stars: winds / outflows
© ESO 2019
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.