3. Nuclear Reactors 6 pts

Modern Physics · Nuclear physics, Elastic collisions, Ideal gas

Neutron moderation in a thermal reactor — speeds, temperatures, number of collisions to thermalise, and gas buildup in spent fuel rods.

Solution and grading scheme by Topi Lind, Melvin Storbacka, Oleg Kosik, Aleksandr Sorokin and Ludmila Belogrudova.

Part i (1 point)

If the speed of the particle is much less than the speed of light, we can use the non-relativistic approach. For non-relativistic particles we know v=2Ek/mv = \sqrt{2E_k/m}. Substituting values gives us vf=2.2×103ms1v_f = 2.2 \times 10^3\,\text{m\,s}^{-1}. This is much less than the speed of light and thus justified. Another way to justify the applicability of the non-relativistic approach would be to say that kinetic energy is significantly less than the rest energy (Efmnc2E_f \ll m_n c^2).

Typo in problem description. Ef=0.025eVE_f = 0.025\,\text{eV} is the mode of the Maxwell–Boltzmann distribution which gives E=kBTE = k_BT. The average gives E=32kBTE = \frac{3}{2}k_BT.

Using the mode of the Maxwell–Boltzmann distribution, E=kBTE = k_BT, and converting from eV to J:

Tf=0.025×1.602×10191.38×1023=290K.T_f = \frac{0.025 \times 1.602 \times 10^{-19}}{1.38 \times 10^{-23}} = 290\,\text{K}.

With E=32kBTE = \frac{3}{2}k_BT we find Tf193KT_f \approx 193\,\text{K}.


Part ii (1 point)

The non-relativistic approach is justified as in the previous task. The same approach gives us v0=2.0×107ms1v_0 = 2.0 \times 10^7\,\text{m\,s}^{-1}.

Grading for parts i and ii (preliminary)

  • Expresses vf=2Ef/mv_f = \sqrt{2E_f/m} and/or v0=2E0/mv_0 = \sqrt{2E_0/m}: 0.3 pts
  • Calculates vf=2.2×103ms1v_f = 2.2\times 10^3\,\text{m\,s}^{-1} and/or v0=2.0×107ms1v_0 = 2.0\times 10^7\,\text{m\,s}^{-1} for two correct numerical values: 0.5 pts; if only one value is correct: 0.3 pts
  • Uses Ef=32kBTE_f = \frac{3}{2}k_BT: 0.3 pts
  • Using this formula calculates T=193KT = 193\,\text{K}: 0.3 pts
  • Justifies the validity of the classical approach in both cases: 0.3 + 0.3 pts

Remark. Using Ef=kBTE_f = k_BT without justification and finding T=290KT = 290\,\text{K} gives 0 pts for the formula and 0.3 pts for the numerical calculation.


Part iii (2.5 points)

In a collision between particle 1 (m1m_1, v1,iv_{1,i} and v1,fv_{1,f}) and particle 2 (m2m_2, v2,iv_{2,i} and v2,fv_{2,f}) momentum is conserved,

m1(v1,fv1,i)=m2(v2,iv2,f),m_1(v_{1,f} - v_{1,i}) = m_2(v_{2,i} - v_{2,f}),

and since the collisions are elastic, kinetic energy is also conserved:

m1(v1,f2v1,i2)=m2(v2,i2v2,f2).m_1(v_{1,f}^2 - v_{1,i}^2) = m_2(v_{2,i}^2 - v_{2,f}^2).

Dividing the latter by the former leads to

v1,i+v1,f=v2,i+v2,f.v_{1,i} + v_{1,f} = v_{2,i} + v_{2,f}.

Substituting this into the conservation of momentum gives for particle 1:

v1,f=m1m2m1+m2v1,i+2m2m1+m2v2,i,v_{1,f} = \frac{m_1 - m_2}{m_1 + m_2}v_{1,i} + \frac{2m_2}{m_1 + m_2}v_{2,i},

and similarly for particle 2:

v2,f=2m1m1+m2v1,i+m1m2m1+m2v2,i.v_{2,f} = \frac{2m_1}{m_1 + m_2}v_{1,i} + \frac{m_1 - m_2}{m_1 + m_2}v_{2,i}.

We see that with m1=m2m_1 = m_2 there is a maximum transfer of momentum. Assuming that particle 1 is the neutron and particle 2 is the target, and that the target is at rest for all intents and purposes, the mass of the moderator atoms should be the same as the neutron’s.

In a single collision with a stationary atom of the moderator, the speed of the neutron decreases by the factor m1m2m1+m2\frac{m_1 - m_2}{m_1 + m_2}, so the speed of the neutron after NN head-on collisions with stationary atoms of moderator will be

vN=v(m1m2m1+m2)N.v_N = v\left(\frac{m_1 - m_2}{m_1 + m_2}\right)^N.

Hence,

N=ln(vf/v0)ln[(m1m2)/(m1+m2)]=12ln(Ef/E0)ln[(m1m2)/(m1+m2)]=614.N = \frac{\ln(v_f/v_0)}{\ln[(m_1 - m_2)/(m_1 + m_2)]} = \frac{1}{2}\cdot\frac{\ln(E_f/E_0)}{\ln[(m_1 - m_2)/(m_1 + m_2)]} = 614.

Grading (preliminary)

  • TTfT \ll T_f, so moderator atoms are essentially at rest: 0.3 pts
  • Justifies that maximum momentum transfer is when mn=Mm_n = M: 0.4 pts
  • Applies energy and momentum conservation: 0.3 + 0.3 pts
  • Expresses u=vm1m2m1+m2u = v\frac{m_1 - m_2}{m_1 + m_2}: 0.4 pts
  • Expresses vf=v0(m1m2m1+m2)Nv_f = v_0\left(\frac{m_1 - m_2}{m_1 + m_2}\right)^N: 0.5 pts
  • Calculates N=614N = 614: 0.3 pts

Part iv (1.5 points)

We can model the gas inside the rod as an ideal gas. Simply due to swelling, the pressure inside the rod would increase from 2.5 MPa to 5 MPa as we know from Boyle’s law P1V1=P2V2P2=P1V1/V2P_1V_1 = P_2V_2 \Rightarrow P_2 = P_1V_1/V_2. Thus, the release of xenon must contribute 1.5 MPa’s worth of pressure due to Dalton’s law ptot=pHe+pXep_\text{tot} = p_\text{He} + p_\text{Xe}. From the ideal gas law we find the amount of xenon moles as

nXe=pXeVRT0=5.5×103mol,n_\text{Xe} = \frac{p_\text{Xe} V}{RT_0} = 5.5\times 10^{-3}\,\text{mol},

where pXe=1.5MPap_\text{Xe} = 1.5\,\text{MPa}, V2=9cm3V_2 = 9\,\text{cm}^3 and T=293KT = 293\,\text{K}. In a similar manner we find that the amount of helium in the beginning was

nHe=pHeV0RT0=1.8×102mol,n_\text{He} = \frac{p_\text{He} V_0}{RT_0} = 1.8\times 10^{-2}\,\text{mol},

where PHe=2.5MPaP_\text{He} = 2.5\,\text{MPa}, V0=18cm3V_0 = 18\,\text{cm}^3 and T0=293KT_0 = 293\,\text{K}. The ratio of the two is

nHenXe3.3.\frac{n_\text{He}}{n_\text{Xe}} \approx 3.3.

Grading (preliminary)

  • Applies Boyle’s law: 0.3 pts
  • Applies Dalton’s law: 0.4 pts
  • Expresses nXe=pXeV/(RT0)n_\text{Xe} = p_\text{Xe}V/(RT_0): 0.2 pts
  • Calculates nXe=5.5×103moln_\text{Xe} = 5.5\times 10^{-3}\,\text{mol}: 0.2 pts
  • Expresses nHe=pHeV0/(RT0)n_\text{He} = p_\text{He}V_0/(RT_0): 0.2 pts
  • Calculates nHe/nXe=3.3n_\text{He}/n_\text{Xe} = 3.3: 0.2 pts