Past exam of the mathematics course of the year 2 of the natural sciences course of the University of Cambridge / 2012 / Paper 1 / 2

For a spherically symmetric number density, the radial Laplacian is $r^{-2}{\partial }_r(r^2{\partial }_{r}n)$. Applying separation of variables and choosing the separation constant $-k^2$ givesso the ordinary differential equations areThe radial problem is an eigenvalue problem for $-\Delta$ with Dirichlet boundary conditions, hence $k^2>0$: multiplying the radial equation by $r^2\overline{R}$ and integrating gives $k^2{\int }_0^a{r}^2|R{|}^{2}dr={\int }_0^a{r}^2|R^{\prime }{|}^{2}dr>0$ for a nonzero regular mode vanishing at $a$.

To solve it explicitly, set $R=f/r$, corresponding to the suggested power $p=-1$. The equation reduces to $f^{\prime \prime }+k^{2}f=0$, soFiniteness at the origin forces $B=0$, and the Dirichlet boundary condition at $r=a$ requires $k=m\pi /a$ for an integer $m\ge 1$. The absence of interior zeros excludes $m\ge 2$, selecting the fundamental normal mode. If $n_c>0$ denotes the initial central density, the result iswhere the sine ratio is defined as $1$ at the origin. It is strictly positive in the interior for every finite time. The arbitrary positive scale $n_c$ cannot be fixed without additional initial data.

The spatial factor is time independent, so the density and the total neutron population have the same exponential growth rate. Thus the spherical reaction-diffusion growth threshold isAt equality the fundamental mode is stationary; below it the mode decays. For fixed $\lambda >0$, the equivalent critical-radius condition is $a>\pi /\sqrt{\lambda }$. This is the competition between multiplication in the bulk and diffusion loss at the surface in the stated dimensionless model.