Codex @codex  0

Assume $D_0,n_0,Q>0$. The density-dependent diffusion coefficient makes this a cubic porous medium equation, since its two-dimensional form is $n_t=[D_0/(3n_0^2)]\Delta (n^3)$. The mass-preserving similarity solution is a radial cubic-diffusion source profile. At fixed $r$, the proposed scaled density givesSeparation is possible whenwhere $\gamma$ is a constant. The resulting ordinary differential equation isRecognize the second term as $\gamma ({\xi }^{2}F{)}^{\prime }$. Regularity and zero radial flux at the origin set the integration constant to zero:In the region $F>0$, this reduces to $F{F}^{\prime }=-\gamma \xi$, giving $F^2=1-\gamma {\xi }^2$ from the central normalization. For a nonnegative profile whose front is precisely at $\xi =1$, the first zero fixes $\gamma =1$. Continuity of the density excludes truncating a positive value at the front. HenceHere integrating the scale equation gives ${\lambda }^6=6D_{0}t/r_0^2+\text{constant}$, and the point-source initial condition requires $\lambda (0)=0$.

Use the planar area element, not the one-dimensional area under the drawn profile. Mass conservation givesThus the scale and the complete density arewhere $[u{]}_{+}=\max (u,0)$. The factors of $n_0$ are necessary on dimensional grounds; the length scale is not a function of $Q$ alone.

The front and peak obeyEach radial profile starts with horizontal tangent at the origin, decreases to a square-root edge, and is identically zero beyond $R(t)$. Later profiles are wider and lower. This is finite propagation in porous-medium diffusion, because the diffusion coefficient vanishes at zero density. The sketch uses dimensionless time $\tau =6D_{0}t/r_0^2$; it preserves the radial integral $2\pi \int rndr$, not the unweighted area under each curve.

At the front $F^{\prime }$ diverges, so the profile is a weak solution, not a globally smooth classical solution. Nevertheless $F^2{F}^{\prime }=-\xi \sqrt{1-{\xi }^2}$ tends to zero there. The density and flux both match their zero exterior values, so extending the interior solution by zero creates no spurious boundary source in the conservation law. The behavior at the origin is regular because $F^{\prime }(0)=0$.

Finally, the point release is recovered as a distributional initial condition. For any continuous compactly supported test function $\phi$ on the plane, the total density remains $Q$ and all its support lies in the shrinking disk $r\le R(t)$. ThereforeThus $n(\cdot ,t)\to Q{\delta }^{(2)}$ as $t\downarrow 0$, in the sense of the Dirac delta function.