Past exam of the mathematics course of the University of Cambridge / 2016 / iii / Paper 321 / 1 / b

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.

Use characteristic, vertically averaged quantities at fixed radius. Let $H$ be a disk thickness, $\rho \sim \Sigma /H$, and let the alpha disk prescription give $\mu \sim \alpha p/\Omega$, up to a convention-dependent constant. From the preceding radiation pressure and constant-opacity balance, $\mu \sim c/\kappa$ is independent of $\Sigma$; thusVertical hydrostatic equilibrium gives $p\sim \rho {\Omega }^2{H}^2\sim \Sigma {\Omega }^{2}H$, henceThe vertically averaged kinematic viscosity satisfies $\bar{\nu }\Sigma =\int \mu dz\sim \mu H$. ConsequentlyAt fixed radius the omitted factors are independent of $\Sigma$, so this is the requested inverse-density scaling of the viscous torque in an accretion disk.

This is an order-of-magnitude vertical model, as requested. Imposing both exact pointwise $\mu =\text{constant}$ and exact pointwise $\mu \propto p(z)$ would force $p$ to be constant, contrary to nontrivial vertical hydrostatic equilibrium. The characteristic-pressure prescription supplies a scaling, not an exact simultaneous vertical solution to all those idealized equations.

Write $W(\Sigma )=\bar{\nu }\Sigma$. In a small radial patch $r=r_c+\chi$ with disturbance scale much smaller than $r_c$, the leading disk diffusion equation is ${\partial }_t\Sigma \simeq 3{\partial }_{\chi }^{2}W$. Taking the linearization around a uniform surface density ${\Sigma }_c$ givesFor a local Fourier mode $\delta \Sigma \propto e^{st+ik\chi }$,because $W\propto {\Sigma }^{-1}$ has $W^{\prime }<0$. Thus the effective density diffusion coefficient is negative and the disk has a viscous instability of an accretion disk.

A localized overdensity is a maximum of $\Sigma$ but a minimum of $W$. At its center ${\partial }_{\chi }^{2}W>0$, so ${\partial }_t\Sigma >0$: it grows. The local radial mass flux is proportional to $-{\partial }_{\chi }W$, directing matter toward that minimum from both sides. This is the physical feedback: an overdense patch has a smaller viscous torque in an accretion disk and attracts further accumulation. The arbitrarily rapid $k\to \infty$ growth in this ideal negative-diffusion model is not a resolved prediction below the validity scale of the thin-disk and local viscous closure.