 Solution

The four equations describe, respectively, vertical hydrostatic equilibrium in the central object's gravity, the vertical divergence of radiative flux balancing viscous dissipation, optically thick radiative diffusion, and the equation of state combining ideal gas gas pressure with radiation pressure. Here $\mu$ is dynamic viscosity, not the mean molecular weight ${\mu }_m$; $\kappa$ is the opacity, $\sigma$ the Stefan-Boltzmann constant, and $\Omega$ the local orbital frequency.

In the radiation-pressure limit, $p=4\sigma T^4/(3c)$, so radiative diffusion givesCompare this with vertical hydrostatic equilibrium. Wherever $\rho \ne 0$,At fixed radius, constant opacity makes $dF/dz=c{\Omega }^2/\kappa$. The viscous dissipation equation therefore yieldsindependent of height within the optically thick radiation-supported layer. The result assumes the supplied gravity, heating and diffusion equations; optically thin surface corrections are outside that approximation.