Past exam of the mathematics course of the University of Cambridge / 2017 / iii / Paper 321 / 1 / a

Take $0<r_1<r_2$ and positive mean kinematic viscosity $\bar{\nu }$. The steady state of the supplied viscous evolution of an accretion disk equation satisfies ${\mathcal{F}}^{\prime }=-2\pi rS(r)$. The outer no-flux condition therefore givesPositive $S$ produces positive inward accretion rate in this sign convention. Set $Y=r^{1/2}\bar{\nu }\Sigma$. In a Keplerian disk, the viscous torque in an accretion disk is proportional to $Y$, so the zero-torque inner boundary condition gives $Y(r_1)=0$. Since $Y^{\prime }=\mathcal{F}/(6\pi r^{1/2})$, a second integration yieldsThis is the steady viscous disk supplied by a radial source. It does not require $\bar{\nu }$ to be independent of radius: the equation directly determines $\bar{\nu }\Sigma$. If viscosity depends on the local thermodynamic state, the displayed formula is an implicit relation completed by a separate closure. For nonnegative integrable $S$ it gives nonnegative surface density of a disk, and the inner accretion rate equals the total supplied mass per unit time, $2\pi {\int }_{r_1}^{r_2}rS(r)dr$. Differentiating both nested integrals recovers the flux and the differential equation.

Interpret the outer Dirac delta function as a boundary source carrying its full mass, or as the limit of a narrow source lying just inside $r_2$. Then ${\int }_u^{r_2}vS(v)dv=r_2{S}_0$ for $r_1<u<r_2$, andThus the interior is the usual steady Keplerian accretion disk with a constant inward accretion rate. The product $\bar{\nu }\Sigma$ vanishes at the stress-free inner edge, and the surface density of a disk does too if the viscosity has a finite nonzero boundary value; far from that edge $\bar{\nu }\Sigma$ approaches $\dot{M}/(3\pi )$. Its density need not itself be constant because the viscosity may vary.

The outer-edge source of a viscous disk requires care about the boundary trace. The flux just inside the injecting edge is $\mathcal{F}(r_2^{-})=\dot{M}$, while the no-flux value outside the source is $\mathcal{F}(r_2^{+})=0$; the source produces the intervening jump. A simultaneously continuous interior condition $\mathcal{F}(r_2^{-})=0$ would contradict nonzero steady injection. A symmetrically regularized delta integrated only up to its center can instead give half its full weight; that would mean a different supplied mass, not a different interior solution for the same $\dot{M}$. The physically normalized one-sided limiting source removes this convention ambiguity.