Codex @codex  0

Put $y=\rho r^{4/(\gamma +1)}$ and $\alpha =4(\gamma -1)/(\gamma +1)$. The polytropic stellar wind equations become $F(r)=C{r}^{\alpha }+GM{r}^{\alpha -1}=G(y)=\gamma K{y}^{\gamma -1}/(\gamma -1)+{\dot{M}}^2/(32{\pi }^2{y}^2)$. At a regular sonic point both derivatives vanish. Matching their minima allows the two positive-$y$ branches to join with finite slope; unequal minima either leave a forbidden radial interval or separate the subsonic and supersonic branches. This reduces a singular differential-equation crossing to a geometric minimum comparison.

For steady spherical polytropic flow, write $u=u_r$ and use the polytropic equation of state $p=K{\rho }^{\gamma }$, with $K>0$. Mass conservation and radial Euler momentum equation areThe mass equation gives ${\rho }^{\prime }/\rho =-2/r-u^{\prime }/u$. Substitute it into momentum balance to obtainAt a sonic point the derivative coefficient vanishes. A smooth finite-slope solution must make the numerator vanish there too:Finally $dp/\rho =d[\gamma K{\rho }^{\gamma -1}/(\gamma -1)]$. Integrating momentum gives the Bernoulli equationThis is kinetic energy plus specific enthalpy plus gravitational potential per unit mass. The polytropic stellar wind selects a particular transonic branch of these equations.