Subsonic structure and optically thick winds from Wolf–Rayet stars

Mass loss by stellar wind is a key agent in the evolution and spectroscopic appearance of massive main sequence and post-main sequence stars. In Wolf–Rayet stars the winds can be so dense and so optically thick that the photosphere appears in the highly supersonic part of the outflow, veiling the underlying subsonic part of the star, and leaving the initial acceleration of the wind inaccessible to observations. Here we investigate the conditions and the structure of the subsonic part of the outflow of Galactic Wolf–Rayet stars, in particular of the WNE subclass; our focus is on the conditions at the sonic point of their winds. We compute 1D hydrodynamic stellar structure models for massive helium stars adopting outer boundaries at the sonic point. We find that the outflows of our models are accelerated to supersonic velocities by the radiative force from opacity bumps either at temperatures of the order of 200 kK by the iron opacity bump or of the order of 50 kK by the helium-II opacity bump. For a given mass-loss rate, the diffusion approximation for radiative energy transport allows us to define the temperature gradient based purely on the local thermodynamic conditions. For a given mass-loss rate, this implies that the conditions in the subsonic part of the outflow are independent from the detailed physical conditions in the supersonic part. Stellar atmosphere calculations can therefore adopt our hydrodynamic models as ab initio input for the subsonic structure. The close proximity to the Eddington limit at the sonic point allows us to construct a sonic HR diagram, relating the sonic point temperature to the luminosity-to-mass ratio and the stellar mass-loss rate, thereby constraining the sonic point conditions, the subsonic structure, and the stellar wind mass-loss rates of WNE stars from observations. The minimum stellar wind mass-loss rate necessary to have the flow accelerated to supersonic velocities by the iron opacity bump is derived. A comparison of the observed parameters of Galactic WNE stars to this minimum mass-loss rate indicates that these stars have their winds launched to supersonic velocities by the radiation pressure arising from the iron opacity bump. Conversely, stellar models which do not show transonic flows from the iron opacity bump form low-density extended envelopes. We derive an analytic criterion for the appearance of envelope inflation and of a density inversion in the outer sub-photospheric layers.

Accretion onto a noncommutative-inspired Schwarzschild black hole

In this paper, we investigate the problem of ordinary baryonic matter accretion onto the noncommutative (NC) geometry-inspired Schwarzschild black hole. The fundamental equations governing the spherically symmetric steady state matter accretion are deduced. These equations are seen to be modified due to the presence of noncommutativity. The matter accretion rate is computed and is found to increase rapidly with the increase in strength of the NC parameter. The sonic radius reduces while the sound speed at the sonic point increases with the increase in the strength of noncommutativity. The profile of the thermal environment is finally investigated below the sonic radius and at the event horizon and is found to be affected by noncommutativity.

Exact solutions to non-classical steady nozzle flows of Bethe–Zel’dovich–Thompson fluids

Steady nozzle flows of Bethe–Zel’dovich–Thompson fluids – substances exhibiting non-classical gasdynamic behaviour in a finite vapour-phase thermodynamic region in close proximity to the liquid–vapour saturation curve – are examined. Non-classical flow features include rarefaction shock waves, shock waves with either upstream or downstream sonic states and split shocks. Exact solutions for a mono-component single-phase fluid expanding from a reservoir into a stationary atmosphere through a conventional converging–diverging nozzle are determined within the quasi-one-dimensional inviscid flow approximation. The novel analytical approach makes it possible to elucidate the connection between the adiabatic, possibly non-isentropic flow field and the underlying local isentropic-flow features, including the possible qualitative alterations in passing through shock waves. Contrary to previous predictions based on isentropic-flow inspection, shock disintegration is found to occur also from reservoir states corresponding to a single sonic point. The global layout of the flow configurations produced by a monotonic decrease in the ambient pressure, namely the functioning regime, is examined for reservoir conditions resulting in single-phase flows. Accordingly, a classification of steady nozzle flows into 10 different functioning regimes is proposed. Flow conditions determining the transition between the different classes of flow are investigated and each functioning regime is associated with the corresponding thermodynamic region of reservoir states.

A Plane-Parallel Wind Solution for Testing Numerical Simulations of Photoevaporation

Here, we derive a Parker-wind-like solution for a stratified, plane-parallel atmosphere undergoing photoionisation. The difference compared to the standard Parker solar wind is that the sonic point is crossed only at infinity. The simplicity of the analytic solution makes it a convenient test problem for numerical simulations of photoevaporation in protoplanetary discs.

P-mode oscillation on slim discs

We numerically investigate the thermally unstable accretion discs around spinning black holes with different spins. We adopted an additional evolutionary viscosity equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find an interesting oscillation when accretion switches to slim disc mode. The oscillation arises from the sonic point of accretion flow and propagates outwards. We mimic the bolometric light-curve and find a series of harmonics on its power spectrum. The frequency ratio of those harmonics is a regular integer series. The lowest frequency of the harmonics is identical to the prediction of trapped p-mode in QPO theory.