Flux and current in a spherical fusion chamber

This problem is from NUCL 510, but plays off of the material in (Lamarsh, 2002).

Neutrons are produced uniformly and isotropically throughout a spherical chamber containing a fixture of \(\ce{^3}{H}\) and \(\ce{^{2}H}\) gasses at \(1\times 10^{8} \unit{K}\). Calculate the flux and current in the chamber.

We know that there are fusion reactions

\[\ce{T}\left(d, n\right) \ce{^{4}He}\]

and

\[\ce{D}\left(d, n\right)\ce{^{3}He}\]

.

First we set up the geometry. We have a sphere with its \(z\) axis through point \(A\). \(dV\) is the volume of the skirt \(\overline{AB}\). so we have

\[d\psi = \frac{S}{4 \pi r^{2}} dV\]

We know that \(\phi = \int_{V} d\psi\), so we integrate from \(0\) to \(2R \cos \theta\).

\[\phi = \int_{0}^{\nicefrac{pi}{2}} \int_{0}^{2R\cos \theta} \frac{S}{4\pi r^{2}} \underbrace{r^{2} \sin \theta dr d\theta d\varphi}_{\text{spherical lagrangian}}\]
\[\phi = \frac{S}{2}\int_{0}^{1} \int_{0}^{2R\mu} drd\mu\]
\[\phi = SR\int_{0}^{1}\mu d\mu\]
\[\phi = \frac{SR}{2}\]

For current, we need to apply the solid angle

\[\vec{\Omega} = \left( \sin \theta \cos \varphi, \sin \theta \sin \varphi, \cos \theta \right)\]
\[\vec{J} = \int_{0}{\nicefrac{\pi}{2}} \int_{0}{2R\cos \theta} \int_{0}^{2\pi} \frac{S}{4\pi r^{2}} \vec{\Omega} r^{2}\sin \theta d\varphi dr d\theta\]

and with \(\int_{0}^{2\pi} \cos \varphi d\varphi=\int_{0}^{2\pi} \sin \varphi d\varphi=0\), therefore \(J_{x} = 0\) and \(J_{y}=0\)

\[J_{z} = \int_{0}{\nicefrac{\pi}{2}} \int_{0}{2R\cos \theta} \int_{0}^{2\pi} \frac{S}{4\pi } \cos \theta \sin \theta d\varphi dr d\theta\]
\[J_{z} = \frac{S}{2}\int_{0}^{1} \int_{0}^{2R\mu} \mu dr d\mu\]
\[J_{z} = \frac{SR}{3}\times \vec{e_{z}} = \frac{SR}{3}\times \vec{e_{r}}\]

Bibliography