Subfield: Modern Physics

Uranium occurs in nature as $\mathrm{UO}_{2}$ with only $0.720 \%$ of the uranium atoms being ${ }^{235} \mathrm{U}$. Neutron induced fission occurs readily in ${ }^{235} \mathrm{U}$ with the emission of 2-3 fission neutrons having high kinetic energy. This fission probability will increase if the neutrons inducing fission have low kinetic energy. So by reducing the kinetic energy of the fission neutrons, one can induce a chain of fissions in other ${ }^{235} \mathrm{U}$ nuclei. This forms the basis of the power generating nuclear reactor (NR).

A typical NR consists of a cylindrical tank of height $H$ and radius $R$ filled with a material called moderator. Cylindrical tubes, called fuel channels, each containing a cluster of cylindrical fuel pins of natural $\mathrm{UO}_{2}$ in solid form of height $H$, are kept axially in a square array. Fission neutrons, coming outward from a fuel channel, collide with the moderator, losing energy and reach the surrounding fuel channels with low enough energy to cause fission (Figs I-III). Heat generated from fission in the pin is transmitted to a coolant fluid flowing along its length. In the current problem we shall study some of the physics behind the (A) Fuel Pin, (B) Moderator and (C) NR of cylindrical geometry.

<img_4559>

# Schematic sketch of the Nuclear Reactor (NR) 

Fig-I: Enlarged view of a fuel channel (1-Fuel Pins)

Fig-II: A view of the NR (2-Fuel Channels)

Fig-III: Top view of NR (3-Square Arrangement of Fuel Channels and 4-Typical Neutron Paths).

Only components relevant to the problem are shown (e.g. control rods and coolant are not shown).

## A Fuel Pin

| Data <br> for $\mathrm{UO}_{2}$ | 1. Molecular weight $M_{w}=0.270 \mathrm{~kg} \mathrm{~mol}^{-1}$ | 2. | Density $\rho=1.060 \times 10^{4} \mathrm{~kg} \mathrm{~m}^{-3}$ |
| :---: | :--- | :--- | :--- |
|  | 3. Melting point $T_{m}=3.138 \times 10^{3} \mathrm{~K}$ | 4. | Thermal conductivity $\lambda=3.280 \mathrm{~W} \mathrm{~m}^{-1} \mathrm{~K}^{-1}$ |
Context question:
A.1 Consider the following fission reaction of a stationary ${ }^{235} \mathrm{U}$ after it absorbs a neutron of negligible kinetic <br> energy. <br> $\qquad{ }^{235} \mathrm{U}+{ }^{1} \mathrm{n} \rightarrow{ }^{94} \mathrm{Zr}+{ }^{140} \mathrm{Ce}+2{ }^{1} \mathrm{n}+\Delta E$ <br> $m\left({ }^{94} \mathrm{Zr}\right)=93.9063 \mathrm{u} ; m\left({ }^{140} \mathrm{Ce}\right)=139.905 \mathrm{u} ; m\left({ }^{1} \mathrm{n}\right)=1.00867$ u and $1 \mathrm{u}=931.502 \mathrm{MeV} \mathrm{c}{ }^{-2}$. Ignore charge <br> imbalance.
Context answer:
\boxed{$\Delta E=208.684 $}


Context question:
A.2 Estimate $N$ the number of ${ }^{235} \mathrm{U}$ atoms per unit volume in natural $\mathrm{UO}_{2}$.
Context answer:
\boxed{$=1.702 \times 10^{26} $}


Context question:
A.3 Assume that the neutron flux $\phi=2.000 \times 10^{18} \mathrm{~m}^{-2} \mathrm{~s}^{-1}$ on the fuel is uniform. The fission cross-section (effective area of the target nucleus) of a ${ }^{235} \mathrm{U}$ nucleus is $\sigma_{f}=5.400 \times 10^{-26}$ $\mathrm{m}^{2}$. If $80.00 \%$ of the fission energy is available as heat, estimate $Q$ (in $\mathrm{W} \mathrm{m}^{-3}$ ) the rate of heat production in the pin per unit volume. $1 \mathrm{MeV}=1.602 \times 10^{-13} \mathrm{~J}$.
Context answer:
\boxed{$4.917\times 10^8$}


Context question:
A.4 The steady-state temperature difference between the center $\left(T_{c}\right)$ and the surface $\left(T_{s}\right)$ of the pin can be expressed as $T_{c}-T_{s}=k F(Q, a, \lambda)$ where $k=1 / 4$ is a dimensionless constant and $a$ is the radius of the pin. Obtain $F(Q, a, \lambda)$ by dimensional analysis.
Context answer:
\boxed{$T_{c}-T_{s}=\frac{Q a^{2}}{4 \lambda}$}


Context question:
A.5 The desired temperature of the coolant is $5.770 \times 10^{2} \mathrm{~K}$. Estimate the upper limit $a_{u}$ on the radius $a$ of the pin.
Context answer:
\boxed{$a_{u}=8.267 \times 10^{-3}$}


Extra Supplementary Reading Materials:

B The Moderator

Consider the two dimensional elastic collision between a neutron of mass $1 \mathrm{u}$ and a moderator atom of mass $A$ u. Before collision all the moderator atoms are considered at rest in the laboratory frame (LF). Let $\overrightarrow{v_{b}}$ and $\overrightarrow{v_{a}}$ be the velocities of the neutron before and after collision respectively in the LF. Let $\overrightarrow{v_{m}}$ be the velocity of the center of mass (CM) frame relative to LF and $\theta$ be the neutron scattering angle in the CM frame. All the particles involved in collisions are moving at nonrelativistic speeds.
Context question:
B.1 The collision in LF is shown schematically with $\theta_{L}$ as the scattering angle (Fig-IV). Sketch the collision schematically in CM frame. Label the particle velocities for 1,2 and 3 in terms of $\overrightarrow{v_{b}}, \overrightarrow{v_{a}}$ and $\overrightarrow{v_{m}}$. Indicate the scattering angle $\theta$.

<img_4396>

Collision in the Laboratory Frame
1-Neutron before collision
2-Neutron after collision
3-Moderator Atom before collision
4-Moderator Atom after collision
Context answer:
<img_4334>


Context question:
B.2 Obtain $v$ and $V$, the speeds of the neutron and the moderator atom in the CM frame after the collision, in terms of $A$ and $v_{b}$.
Context answer:
\boxed{$v=\frac{A v_{b}}{A+1}$ , $V=\frac{v_{b}}{A+1}$}


Context question:
B.3 Derive an expression for $G(\alpha, \theta)=E_{a} / E_{b}$, where $E_{b}$ and $E_{a}$ are the kinetic energies of the neutron, in the LF, before and after the collision respectively, and $\alpha \equiv[(A-1) /(A+1)]^{2}$,
Context answer:
\boxed{$G(\alpha, \theta)=1-\frac{(1-\alpha)(1-\cos \theta)}{2}$}


Question: B.4 Assume that the above expression holds for $\mathrm{D}_{2} \mathrm{O}$ molecule. Calculate the maximum possible fractional energy loss $f_{l} \equiv \frac{E_{b}-E_{a}}{E_{b}}$ of the neutron for the $\mathrm{D}_{2} \mathrm{O}(20 \mathrm{u})$ moderator.

Solution:
The maximum energy loss will be when the collision is head on ie., $E_{a}$ will be minimum for the scattering angle $\theta=\pi$.

So $E_{a}=E_{\min }=\alpha E_{b}$.

For $\mathrm{D}_{2} \mathrm{O}, \alpha=0.819$ and maximum fractional $\operatorname{loss}\left(\frac{E_{b}-E_{\min }}{E_{b}}\right)=1-\alpha=0.181$. [A $\boldsymbol{c}$ ceptable Range (0.170 to 0.190 )]

Final answer: $f_{l}=0.181$