Hans Moritz Günther edited our_equations.tex  over 10 years ago

Commit id: e91067790c0ff73f15ea737459e25337af3bf0dd

deletions | additions      

       

In this section we develop an analytical model for the stellar wind that expands along the jet axis. The wind is confined by the pressure of the surrounding disk wind. Those two flows a re separated by a contact discontinuity, whose exact position is given by pressure equilibrium between the outer, disk wind component and the inner, stellar wind component. We assume that the stellar wind is initially emitted radially. As it encounters the contact discontinuity, the velocity component perpendicular to the discontinuity is shocked. Thus, our model needs to distinguish three zones: (i) the cold pre-shock stellar wind, (ii) the host post-shock stellar wind, and (iii) the disk wind. Our goal is to calculate the geometrical shape of the stellar wind shock, since this determines the velocity jump across the shock front and thus the temperature of the post-shock plasma.  The Rankine-Hugoniot jump conditions relate the density, velocity and pressure on both sides of a strong shock. For ideal gases and non-oblique shocks the conservation of mass, momentum and energy across the shock can be written as follows \citep[][chap.~7, \S~15]{http://adsabs.harvard.edu/abs/1967pswh.book.....Z, \S~15]{http://adsabs.harvard.edu/abs/1967pswh.book.....Z},  where the state before the front of the shock front is marked by the index 0, that behind the shock by index 1: \begin{eqnarray}  \rho_0 v_0 &=& \rho_1 v_1 \label{RH1}\\  P_0+\rho_0 v_0^2 &=& P_1+\rho_1 v_1^2 \label{RH2}\\