Past exam of the MPhil in nuclear energy course of the University of Cambridge / 2020

Original paper: Released original PDF, physical pages 28–30; questions on pages 29–30.

Take one-dimensional steady state thermal conduction, constant thermal conductivity $k>0$, and positive volumetric heat generation rate. Write $q_0=q_0^{\prime \prime \prime }$ for brevity, and let $x$ increase away from the illuminated face. Conservation of energy and Fourier's law giveThus $T^{\prime \prime }=-(q_0/k)e^{-\mu x}<0$: the profile is strictly concave and any stationary point is its unique maximum. Two integrations and the equal-surface boundary conditions yield the exponential slab heating solutionIts derivative isDefine $\beta =\mu a$ and $C=(1-e^{-\beta })/\beta$. Setting $T^{\prime }(x_m)=0$ gives the distance from the illuminated surface:For $\beta >0$, $e^{-\beta }<C<1$, so $0<x_m<a$. More specifically $C=e^{-\beta /2}\sinh (\beta /2)/(\beta /2)>e^{-\beta /2}$, hence $x_m<a/2$: the peak is nearer the more strongly heated face.

The apparent singularity at $\mu =0$ is removable. Expansion givesThe uniform-heating limit therefore has its maximum at $a/2$, as symmetry requires.

At the maximum found in part(a), $e^{-\mu x_m}=C$ and $(1-e^{-\beta })x_m/a=C\mu x_m=-C\ln C$. Substitution into the exponential slab heating solution givesFor $0<C<1$, the function $1-C+C\ln C$ is positive: it vanishes at $C=1$ and has derivative $\ln C<0$. Thus the interior temperature exceeds the common surface temperature, consistently with positive internal heating.

A convenient dimensionless form uses $\xi =x/a$:As $\beta \to 0$, $\theta \to \xi (1-\xi )/2$, giving the check

The graph holds the normalization $q_0{a}^2/k$ fixed. Its decreasing peaks with increasing $\beta$ reflect a decreasing source away from the illuminated face under that comparison; they do not assert the same trend at fixed total absorbed power.

The two outward surface normals point in opposite directions. With positive heat flux defined as heat leaving each face, Fourier's law givesUsing the derivative from part(a), the removed fluxes per unit face area areBoth are positive because $e^{-\beta }<C<1$. Their sum satisfies the independent conservation of energy checkThe illuminated face removes more heat because heating is concentrated nearer it. Equivalently, weighting each source slice by the fraction of its heat conducted to either equal-temperature boundary givesA decreasing source weights the left face more strongly. In the uniform-heating limit, each face removes $q_{0}a/2$. If a heat removal rate rather than a flux is wanted, multiply the relevant expression by the face area.

The exponential slab heating solution separates three useful controls: the heat source, the conductive path and the boundary temperature. Reduce deposited power, improve thermal conduction, shorten paths to a cooled surface, or lower the actual surface temperature. Each proposal has to preserve the shield's radiation-protection function.

Reduce incident photon intensity, intercept part of the beam in a separately cooled upstream shield, or distribute the deposited power over a larger cooled region. These reduce the relevant volumetric heat generation rate or local heat load. Intermittent exposure can reduce average heating, but the peak then requires a transient heat equation calculation; a smaller average alone does not establish a safe pulse maximum.

Increase the effective thermal conductivity, for example through a suitable material or a well-bonded conductive path to a heat sink. Reduce contact thermal resistance in any added path. Material substitution must also satisfy shielding, mechanical and environmental constraints.

Shorten the maximum conduction distance by introducing cooled internal surfaces or dividing the shield into thinner, separately cooled layers. Reducing the thickness of a single slab lowers its temperature rise at fixed $q_0$ and $\mu$, but also reduces attenuation; it is not automatically an acceptable shield modification. A layered arrangement can preserve the required total shielding while changing the thermal boundary-value problem.

Lower the cooling-air temperature, improve forced convection by appropriate flow and surface design, enlarge the heat-exchange area, or use a more effective coolant where compatible. These can lower the actual boundary temperature when the external cooling is modeled. In the stipulated model $T_0$ is already prescribed: increasing the external heat transfer coefficient while holding $T_0$ fixed leaves the calculated interior profile unchanged. It helps physically by reducing the surface-to-coolant temperature difference under a heat load.

Do not change the attenuation coefficient independently of the heat-source normalization without explaining the comparison. At fixed incident photon energy flux, the entrance deposition $q_0$ generally changes with absorption properties. A larger $\mu$ with artificially fixed $q_0$ is a different experiment. Similarly, increasing total slab heat capacity delays heating but does not lower the steady state maximum for an unchanged source and unchanged thermal boundary conditions.

Use a common comparison for all of question2: retain the same number and height of fuel rods and approximately the same normalized rod-power shape; reduce fuel lattice pitch by making the active core cross-section smaller. Compare otherwise similar operating states of an initially under-moderated reactor, typical of a pressurized water reactor. Rod number, core height, reflector and control settings are not all fixed explicitly by the question, so these are stated assumptions. Keeping the external core envelope fixed while inserting more rods would be a different comparison.

Reactor reactivity normally decreases at unchanged controls and reference temperatures. Less water per fuel atom reduces neutron moderation and increases resonance absorption during slowing down. Reduced absorption in the water partly compensates, but on the under-moderated side the net lattice tendency is lower multiplication. This is the usual under-moderation argument, rather than a universal rule for arbitrary fuel lattices.^[ref]

Neutron leakage would usually tend to increase. The smaller transverse core extent increases the geometric tendency to leak, and reduced moderation generally increases the slowing-down distance. There are competing changes in fuel absorption, transport lengths and reflection. Thus a definite leakage fraction needs the actual core and reflector arrangement; the usual increase is a qualitative tendency, not a geometric theorem for every redesign.

The Doppler coefficient of reactivity remains negative under the ordinary uranium-fuel assumptions. Increasing fuel temperature broadens absorption resonances, increasing effective uranium-238 resonance capture and reducing reactor reactivity. This stabilizing fuel-temperature effect is distinct from moderator-density feedback.^[ref] Greater importance of resonance absorption can make this feedback more important in a tighter lattice, but its change in magnitude is not fixed by pitch alone: spectrum, resonance self-shielding, fuel composition, burnup and the reference fuel temperature all enter the derivative.

The fraction of fissions in uranium-238 tends to increase. The harder neutron spectrum gives more opportunity for fast fission before neutrons slow below the uranium-238 fission threshold. This concerns the fraction of the total nuclear fission rate, not an assertion that the total thermal power increases; that power is held fixed.

Continue the same-rod-number and same-height comparison from part(a), and assume unchanged pump speed, similar inlet conditions and single-phase forced flow for the normal-operation hydraulic argument.

Core power density increases. The active core volume decreases while total thermal power stays fixed. The fuel's own average volumetric power and the mean fuel linear heat rate do not increase merely from bringing the same rods closer together: fuel volume and rod count are unchanged.

Coolant mass flow rate normally decreases. Tighter passages reduce flow area and hydraulic diameter, increasing the core's hydraulic resistance. The unchanged centrifugal pumps do not enforce constant flow. The pump operating point moves to the intersection of the same pump curve with the steeper system-resistance curve, normally giving smaller flow at greater pump head.^[ref]

The core pressure drop normally increases for the usual pump curve. Its resistance at any given flow is higher. At the new, smaller flow, unchanged external piping and component losses are smaller, whereas a conventional falling pump head-versus-flow curve gives greater available head. The core takes the difference. This conclusion uses the stated pump-curve assumptions; holding flow fixed or holding pump head fixed would be different models. The exact new flow and pressure drop need both system and pump characteristics.

Fuel temperature would normally tend to rise, especially downstream, if the loss of flow increases coolant heating without a compensating improvement in film heat transfer. At fixed inlet state,shows that smaller total mass flow rate raises the bulk coolant temperature rise. At unchanged local fuel linear heat rate, the fuel-to-surface conductive contribution is approximately unchanged if material properties and gap conditions are held fixed. Its boundary temperature and the coolant film drop can nevertheless change.

The local coolant mass flux is $G=\dot{m}/A_f$. Both total flow and flow area decrease, so the sign of the change in $G$, the heat transfer coefficient and the fuel temperature cannot be deduced from total flow alone. Changes in power peaking and inlet conditions add further qualifications to the usual hotter-fuel tendency.

The minimum departure from nucleate boiling ratio would commonly decrease, implying reduced boiling margin, but its sign is conditional. The defining local ratio isThe numerator is the critical heat flux for the local conditions, and the denominator is actual rod-surface heat flux.^[ref] The mean surface heat loading remains approximately unchanged for the same rods and total power. Greater coolant enthalpy rise tends to reduce margin, while changes in $G$, geometry, pressure and local power determine the actual critical heat flux. A geometry-appropriate correlation or experiment is needed for a definite MDNBR change; unchanged hardware does not supply it.

A tighter core is likely to have a more demanding cooling response, but the specification does not determine a unique peak cladding temperature or reflood time in a loss-of-coolant accident. Several effects compete.

At the same pre-accident power and comparable irradiation history, the total decay heat source remains of broadly comparable scale; increasing core power density does not multiply the total decay power by the inverse change in core volume. The core nevertheless has less water within its tighter rod lattice. Its local thermal buffer is smaller, and the same heat must be removed through narrower passages. If the normal-operation fuel temperature is higher, more stored fuel thermal energy must also be removed after shutdown.

Reduced passage size and increased hydraulic resistance can impede emergency-coolant penetration and steam discharge during reflood. Local two-phase flow, entrainment and boiling heat transfer can change, so single-phase pump-curve reasoning is insufficient during blowdown and reflood. Less local water can favor earlier uncovering, but blowdown flow and the much larger inventory elsewhere in the unchanged primary circuit also matter. It would be unjustified to infer a proportionate reduction in the entire primary-circuit inventory.

Loss of water also removes neutron moderation. In the assumed under-moderated reactor, this normally contributes negative reactor reactivity, helping suppress fission power. It does not eliminate decay heat after shutdown. The thermal challenge must therefore be assessed separately from the sign of this feedback.

A quantitative comparison needs break characteristics, initial states, emergency injection performance, bypass paths and a transient two-phase calculation. The defensible conclusion is potentially reduced thermal margin and more difficult reflood, requiring reassessment of the emergency core cooling system, rather than an unconditional prediction that every large-break transient becomes worse in every metric.

The conversion ratio usually tends to increase as moderation is reduced. A harder spectrum gives greater importance to capture in fertile fuel, while less water can reduce some parasitic neutron absorption. These tendencies can increase replacement of fissile material during irradiation. Competing leakage, fast-fission and spectral effects mean the exact change is not fixed by the hydrogen-to-heavy-metal ratio alone.

Higher conversion can support greater discharge burnup and a longer nuclear fuel cycle length, but neither follows automatically. The improvement tends to slow the loss of useful fissile inventory during irradiation. The same change of lattice, however, lowers beginning-of-cycle reactor reactivity under the usual under-moderated assumptions. With enrichment unchanged, the available initial reactivity margin may become inadequate before any later benefit from improved conversion can be realized.

Under the same-rod-number and same-height assumptions, the initial heavy-metal mass of nuclear fuel $M_{\mathrm{HM}}$ and thermal power $P$ are unchanged. Thus the rate of accumulating energy-per-mass burnup is unchanged:For a constant-power irradiation interval $t$, $B=Pt/M_{\mathrm{HM}}$. A longer achievable interval would consequently increase the burnup accumulated in that interval. This relation alone says nothing about whether that interval remains possible.

If sufficient initial reactor reactivity, shutdown capability and thermal margins are retained, slower depletion can allow increased discharge burnup and a longer cycle. If initial reactivity is too low, or the reduced boiling margin or fuel limits become controlling, either can instead decrease. Discharge burnup also depends on batch residence and refuelling strategy, so it is not generally identical to burnup accumulated in one operating cycle.

The appropriate qualitative answer is therefore conversion tends upward; achievable discharge burnup and cycle length are conditional on the full reactivity history and operating limits. An increased conversion ratio by itself does not determine a longer cycle.

Original paper: Released original PDF, physical pages 32–35; questions on pages 33–34.

The characteristic effect is an irradiation-induced upward shift of the ductile-to-brittle transition in ferritic steel. A component previously used above its transition can then fail by cleavage fracture at the same service temperature, with much less warning from plastic deformation.

Ferrite has a body-centered cubic lattice. At low temperature, dislocation glide becomes difficult enough that brittle crack propagation can precede substantial plastic relaxation. Irradiation hardening raises the yield strength by introducing obstacles to dislocation glide; the competition with cleavage fracture consequently shifts toward brittle failure. Austenite, with its face-centered cubic lattice, does not normally show this same sharp cleavage transition. This does not make austenitic steel immune to loss of ductility, helium embrittlement or other radiation damage in materials.

Examine the fracture surface for cleavage fracture facets and river patterns, and compare Charpy impact test energy-versus-temperature curves from exposed material and an unirradiated reference of the same composition and heat treatment. The reference is important: a low impact energy from one failed specimen cannot establish an irradiation-induced shift. Transition curves can also use lateral expansion or the proportion of ductile fracture appearance.^[ref]

A direct fracture toughness test over the transition region complements the Charpy impact test; the two measurements are not interchangeable. Combine these observations with exposure history, local temperature, microstructure and crack location. Brittle appearance supports the failure mechanism, while a measured transition shift and suitable controls support its attribution to irradiation.

Neutron capture by boron-10 can generate helium, leading to helium embrittlement, especially at elevated temperature. The relevant nuclear reaction isThe emitted alpha particle eventually becomes a helium atom in the metal. The product branch can also involve gamma emission, without changing this helium-producing balance. The production depends on boron isotopic content and the neutron energy spectrum; a thermal component is particularly relevant to this reaction.

Helium has very low solubility in steel and can be trapped at vacancies, interfaces and grain boundaries. Boron already segregated near grain boundaries can therefore create a locally unfavorable distribution of helium. Stabilized helium–vacancy clusters develop into bubbles or cavities. Under tensile stress, boundary cavities can grow and link, reducing ductility, fracture toughness and creep-rupture life through intergranular fracture. Early thin-foil observations and mechanical tests associated high-temperature ductility loss with helium-containing boundary bubbles and their coalescence.^[ref]

Irradiation temperature controls the transport and clustering of both helium and vacancies. At lower temperature, reduced long-range transport can retain small trapped clusters in the matrix. These can contribute to irradiation hardening, but extensive boundary-cavity growth is kinetically less favorable. At higher temperature, migration, cavity growth and material creep become faster, making intergranular fracture and loss of high-temperature ductility more likely. Trap distribution and prior microstructure can redistribute the gas, so there is no composition-independent critical temperature or universal monotonic law.

The mechanical-test temperature also matters. Material irradiated cold can become embrittled during subsequent hot testing or service when trapped gas and defects become mobile. Distinguish gas production during irradiation from the later conditions under which its mechanical consequences are expressed.

The key imbalance is preferential loss of self-interstitials to biased defect sinks, leaving a vacancy supersaturation that can feed void swelling. The initial displacement cascade creates both vacancies and self-interstitials, often in equal net numbers as Frenkel pairs; it need not create more vacancies than interstitials at birth.

Many nearby Frenkel pairs recombine. Others separate, become trapped or form dislocation loops. Surviving self-interstitials are usually much more mobile than vacancies over the relevant range. Dislocations also have an interstitial bias: their elastic interaction with defects makes them more efficient sinks for self-interstitials. Preferential interstitial removal to the dislocation structure can leave the mobile vacancy population above its thermal equilibrium value.

A minimal radiation damage rate theory makes the distinction explicit:Here $c_i,c_v$ are defect number densities, $G$ is the surviving mobile-pair production rate per volume, $K$ the recombination coefficient, $D_i,D_v$ the diffusion coefficients, and $k_i^2,k_v^2$ the respective sink strengths. The vacancy loss term includes thermal emission in an effective form. Equal production terms do not require equal steady populations or equal partitioning among sinks. This simplified model omits an explicit cluster-size distribution.

A stable cavity receiving a greater vacancy influx than interstitial influx and vacancy emission grows. Vacancies therefore condense into voids, producing void swelling. Helium can stabilize small cavity nuclei, but a growing vacancy void need not be a high-pressure gas bubble. Dislocation loops, precipitates and changes in the dislocation density evolve alongside the voids and modify the competing sink strengths.

There is no violation of atom conservation: atoms represented by removed self-interstitials have joined extended defects or surfaces, while the corresponding unfilled lattice volume can accumulate in cavities. The classic electron-microscope observations of cavities in irradiated type316 steel provide direct evidence for this microscopic basis of swelling.^[ref]

Void swelling is generally favored by an intermediate irradiation temperature window. The balance involves defect mobility, stable nucleation, absorption at defect sinks and thermal emission; a larger vacancy population alone does not guarantee a larger void.

At low temperature, vacancy migration is slow. Self-interstitials may migrate and form dislocation loops, and small clusters can cause considerable irradiation hardening, while vacancies cannot readily feed large growing cavities. Defects retained in the matrix can evolve during later heating.

At intermediate temperature, vacancies can migrate far enough to reach stable cavities. Interstitial bias and the competition among defect sinks sustain a net vacancy flux into these cavities. Once the incubation period has produced suitable nuclei and sink populations, substantial void swelling can develop.

At high temperature, the thermal equilibrium vacancy population and vacancy emission from small cavities rise. A small void's curvature further favors emission, so maintaining stable nuclei becomes harder. Recovery changes the dislocation population, and recombination, gas transport and evolving precipitates also change the flux balance. The irradiation-driven vacancy supersaturation can fall, suppressing void nucleation and net swelling. Existing cavities may still coarsen, and gas bubbles, material creep and other damage can remain important.

Thus low-temperature defect retention, intermediate-temperature cavity growth and high-temperature recovery/emission are different regimes. The boundaries depend on alloy composition, initial microstructure, gas production, accumulated displacements per atom and displacement dose rate. No alloy-independent peak temperature can be inferred from the paper's data.

Irradiation softening occurs when irradiation removes pre-existing strengthening obstacles faster than it creates new ones. It is particularly plausible in an initially hardened alloy, rather than an initially soft annealed reference.

In a material strengthened by cold working, the starting dislocation density is high. Irradiation-generated point defects can assist dislocation climb, rearrangement and annihilation, reducing the contribution of work hardening. A schematic estimate iswhere $\mu$ is the shear modulus, $b$ the magnitude of the Burgers vector, ${\rho }_d$ the dislocation density, and the last two terms represent precipitate and irradiation-obstacle strengthening. It is an illustrative decomposition, not an exact additive law for every alloy. A sufficiently large decrease in the pre-existing contributions can exceed the increase in $\Delta {\sigma }_{\mathrm{irr}}$.

In a material strengthened by precipitation hardening, ballistic mixing can dissolve or disorder the original obstacles. Radiation-enhanced diffusion can change their composition or accelerate coarsening; increasing obstacle spacing can reduce strength. In an ordered alloy, irradiation-induced loss of order can also remove a strengthening contribution. None of these outcomes is inevitable: irradiation can instead create new strengthening precipitates.

The required conditions therefore include an appropriate initial microstructure, enough displacement dose to modify it, and a temperature and displacement dose rate at which obstacle removal or recovery can compete with new defect accumulation. Compare with an unirradiated specimen held for the same time at the same temperature to separate ordinary thermal aging from an irradiation effect.

A drop in yield strength or hardness need not restore ductility: cavities, helium embrittlement or localization of dislocation glide can coexist with softening. Conversely, an initially annealed steel may show strong irradiation hardening under the same exposure conditions.

Irradiation creep adds stress-directed, time-dependent strain driven by the irradiation-generated defect population. It can be appreciable where ordinary thermal material creep would be slow, because irradiation maintains mobile point defects far above their equilibrium populations.

Absorption or emission of vacancies and self-interstitials permits dislocation climb. Applied stress can make differently oriented dislocations or dislocation loops absorb defects at different rates: stress-induced preferential absorption then produces a directional strain. Climb can also allow a pinned dislocation to bypass an obstacle and continue dislocation glide, providing a second route to strain.

A useful small-strain empirical form, over a limited range of conditions, isfor constant stress $\sigma$ and approximately constant irradiation-creep compliance $B$. Here $d$ is the accumulated displacements per atom and $\dot{d}$ the displacement dose rate; $B$ is measured per unit stress and per unit dpa. This relation explains why creep can scale approximately with both stress and neutron flux when the damage spectrum is fixed. It is not a constitutive law valid at every dose or temperature.

The total strain can also include thermal material creep, void swelling, irradiation growth and changes due to precipitation or recovery. Stress-free swelling is not itself irradiation creep. Experiments on pressurized type316 tubes separated these contributions and found that correcting stress-enhanced swelling materially affected comparisons between annealed and specimens strengthened by cold working.^[ref]

At sufficiently high temperature, thermal material creep can dominate. Changes in microstructure caused by irradiation hardening, softening, bubbles and cavities can either alter the creep rate or shorten rupture life. Thus enhanced defect-driven creep is the central effect, while the net rate and lifetime remain material- and history-dependent.

Assume a homogeneous Fe–Cr–Ni alloy, normal incidence, independent scattering events and the stated effective scattering cross-sections. Neglect absorption, other constituents and changes in neutron energy within the thin sheet. The atomic fractions are not supplied; take $x_{\mathrm{Fe}}=0.70$, $x_{\mathrm{Cr}}=0.20$, $x_{\mathrm{Ni}}=0.10$ as an explicit representative assumption, not a composition measured or specified in the paper.

The conventional cell of a face-centered cubic lattice containsatoms. The number density is thereforeFor a mixture, the scattering probability per unit path length is ${\Sigma }_s={\sum }_j{n}_j{\sigma }_j=n{\bar{\sigma }}_s$, whereThus ${\Sigma }_s=27.751{\mathrm{m}}^{-1}$. For thickness $t=10^{-3}\mathrm{m}$, the unscattered fraction $U$ satisfies $dU/dz=-{\Sigma }_{s}U$ with $U(0)=1$. The fraction undergoing at least one neutron scattering event isThe thin-sheet approximation gives $P_s\simeq {\Sigma }_{s}t=2.78\%$, consistent to the accuracy of the assumed composition. A simple iron-dominated estimate ${\bar{\sigma }}_s\simeq 3.2$ barns gives essentially the same about $3\%$ scattered conclusion.

For any other assumed composition, use the displayed weighted formula; mass fractions must first be converted to atomic fractions. Since the supplied elemental cross-sections span $2.8$–$4.1$ barns, any mixture of these three elements at the stated number density gives approximately $2.39$–$3.48\%$. This is a bound within the idealized mixture model, not a claim that every such composition is stable austenite.

Use the same atomic fractions and collision assumptions as in part(c). Let $A$ be the irradiated area. The number of atoms in the sheet is $N_t=ntA$, while the number of incident neutrons is $\Phi A$, where $\Phi$ is the neutron fluence.

Counting one primary knock-on atom for each neutron that scatters and using $\nu =430$ displacements per recoil gives the first-collision estimateHence the nominal damage is about $55$ displacements per atom under the assumed composition. In the thin-sheet approximation the density and thickness cancel:The $1.4\%$ difference comes from using first-collision attenuation rather than the linear approximation. Treating repeated collisions after deflection would require angular and energy information not provided; this correction is already smaller than the precision warranted by the composition assumption.

There is also a counting convention in “other atoms”: $430$ counts secondary displacements, while including the displaced primary knock-on atom itself gives $\nu =431$. The corresponding first-collision result is $54.58$ rather than $54.45$ dpa, an increase of only $0.23\%$. Both conventions round to about $55$ dpa. The numerical estimate is explicitly conditional on assigning the stated displacement yield to each scattering event.

A value greater than one is physically possible. Displacements per atom counts displacement events relative to the initial atom inventory: atoms can be displaced repeatedly, and many vacancies and self-interstitials recombine. It does not mean $55$ surviving vacant sites per original lattice site, or that every atom has moved exactly $55$ times. This is a large nominal exposure for which evolving dislocation loops, cavities and irradiation hardening matter; the surviving microstructure additionally depends on irradiation temperature, displacement dose rate, recombination and defect sinks.

Original solutions to both printed questions. The University of Cambridge release occupies physical PDF pages 37–40, with the question body on pages 38–39. ^[ref]

The fuel-cycle account uses the stated assays and explicitly identifies the assumed full-power year. The chemical-separation answer has a visible conceptual scope limitation.

Conceptual scope: the safe account below describes the functions of spent nuclear fuel management and the chemical basis of partitioning. It does not supply an ordered fissile-separation flowsheet, reagents, process conditions or purification instructions. In particular, the requested operational eight-step separation scheme is outside this answer's scope.

Eight useful functional objectives are: receiving and identifying the spent nuclear fuel; managing residual decay heat; providing radioactive material confinement and radiation shielding; establishing the chemical and radiological inventory; partitioning reusable constituents from waste at a conceptual level; maintaining accountable control of recovered materials; conditioning solid and liquid wastes; and retaining or treating radioactive effluents. These are management and containment functions, not a claim that the question specifies one universal eight-stage plant design. Several occur throughout a facility's operation rather than sequentially.

The chemical principle of uranium/plutonium partitioning is that an element's oxidation state and complexation affect its chemical speciation, solubility and relative affinity for immiscible phases. solvent extraction can therefore separate elements whose chemistry differs, even though it does not distinguish uranium isotopes by this chemical selectivity. A recovery operation must account for the accompanying fission products, contaminated process materials and secondary wastes; the recovered material is only one output. radioactive material confinement, remote handling and nuclear material accountancy remain necessary because chemical partitioning redistributes both radioactive and fissile inventories.

Use the thermal power when applying nuclear fuel burnup. No capacity factor is supplied, so take a full-power year of $365$ days; if the operating fraction were $u$, the annual quantities in this part would be multiplied by $u$.The corresponding initial heavy-metal nuclear fuel requirement isAll nuclear fuel and concentrate masses below are uranium masses, excluding oxygen. Let $F,P,W$ denote uranium feed, product and enrichment tails, and $x_f,x_p,x_w$ the uranium-235 mass fractions. Conservation of total uranium and uranium-235 givessoWith $P=M$, the uranium oxide concentrate requirement is thereforeThis is not the gross mass of oxide concentrate. Both $F-P-W=0$ and $0.007F-0.050P-0.003W=0$ provide independent balance checks.

The recycle feed is the recovered uranium, not the whole initial nuclear fuel charge. After the stated loss it isUsing the same enrichment tails assay and the stipulated higher recycle-product assay, the isotopic mass balance givesTreat the $5.3\%$ recycled product as replacing the same uranium mass of ordinary $5.0\%$ nuclear fuel: the additional uranium-235 compensates for the supplied uranium-236 penalty. It displaces fresh natural-uranium feed in the ratio calculated in part (b), henceThe remaining natural feed is $313.36245\mathrm{t}\mathrm{U}\mathrm{y}{\mathrm{r}}^{-1}$, a 9.885% resource saving relative to the full-power once-through requirement. This credits one uranium recycle only; it adds no second-pass recycling or plutonium credit.

Separative work must be charged for producing the recycled nuclear fuel as well as credited for displaced fresh nuclear fuel. With uranium-235 mass fraction $x$, use the separative value functionMasses in kilograms give work in kg-SWU. For ordinary fresh product the work per product kilogram isThus the fresh work avoided isThe recycled product still requiresThe net separative-work saving is consequentlyAs a whole-cycle check, the fresh route requires $215668.0851$ kg-SWU per year, while the mixed fresh-plus-recycled supply requires $214075.3261$ kg-SWU. The net saving is only 0.7385% of the original work. The recovered stream has a higher uranium-235 assay than natural feed, but the higher required product assay uses most of the work credit. These equations are resource accounting and do not specify an uranium enrichment plant or its operation.

Resource savings do not by themselves establish a financial saving. The one-recycle case avoids purchasing $34374.04$ kg U in natural feed and $1592.76$ kg-SWU per year. Those are the benefits to compare with the incremental costs, not the whole value of the recovered uranium or the gross avoided fresh-product separative work.

Let $p_U$ be the concentrate price per kilogram of uranium, $p_S$ the price per kg-SWU, and $C_{\mathrm{add}}$ the annual net extra cost of recovery, recycle-specific conversion and fabrication, waste management and financing, relative to the once-through alternative. Recycling is cheaper under this simplified boundary only ifEquivalently, if $p_S$ and $C_{\mathrm{add}}$ were supplied, the break-even uranium price would beAn avoided spent-fuel-disposal charge belongs as a credit in $C_{\mathrm{add}}$; a new radioactive waste conditioning or disposal liability belongs as a cost. They must not both be assigned the same sign simply because each concerns waste.

The natural-uranium saving is appreciable, but the fractional separative-work saving is small. Recovery losses and the uranium-236 compensation limit the value of recycling. Additional radiological handling, separate fabrication campaigns, secondary waste streams, capital tied up during storage and recovery, and uncertain future disposal costs can matter. Conversely, a higher natural-uranium price, a lower incremental recovery cost or an independently justified disposal credit can favour reuse. The calculation gives a conditional comparison, not a uniquely determined economic verdict or market-price threshold, because no costs are supplied. plutonium utilisation is a different option and has not been silently credited to this uranium-only recycle.

A parent initially at rest emits an alpha particle and a recoiling heavy daughter. momentum conservation gives equal and opposite momenta, so their nonrelativistic kinetic energies satisfyHere $Q$ is the kinetic energy available to this branch; excitation of the daughter reduces it and may be followed by photon emission. For an illustrative $Q=5$ MeV, $m_{\alpha }=4\mathrm{u}$ and $M_d=236\mathrm{u}$, the energies are $4.917$ MeV in the alpha particle and $83.3$ keV in the daughter. The heavy daughter receives much less energy, but that does not imply much less atomic damage.

The fast alpha particle loses most of its energy through electronic stopping: ionization and electronic excitation along a relatively long path. It also undergoes nuclear collisions, particularly towards the end of its range, and can produce displacement damage. The slower, heavy, highly charged recoil couples strongly to surrounding nuclei and generates a compact displacement cascade. Successive knock-on atoms create vacancies, interstitial defects and clusters. Some Frenkel pairs recombine as the cascade relaxes; the initial collision count is therefore not the number of defects that remain.

A quantitative illustration comes from brannerite rather than an assumed universal ceramic. Simulations in the primary natural-mineral study used a $4.5$ MeV alpha particle and an $80$ keV daughter recoil. The corresponding ranges were about $13$ micrometres and $23$ nanometres. At the study's assumed $50$ eV threshold displacement energy, its tabulated results give roughly $560$–$571$ recoil atomic displacements and $91$–$100$ alpha particle atomic displacements per radioactive decay. The numbers depend strongly on the displacement threshold and material, but the shorter recoil path contains the dominant and much denser ballistic radiation damage.^[ref]

In a susceptible ceramic, cascade cores can accumulate into an amorphous fraction; defects and associated relaxation can also change volume, thermal conductivity and mechanical response. molecular dynamics simulations of zirconolite found a strongly disordered core surrounded by more isolated point defects, supporting this distinction between local cascade damage and an entirely transformed specimen.^[ref]

The alpha particle eventually becomes a helium atom after slowing and charge exchange. helium accumulation can have later microstructural effects distinct from the initial knock-on cascade. Material-dependent recombination and annealing limit surviving disorder; a change in diffraction or radiation-induced swelling is not automatically radiation-induced amorphization, and radiation-induced amorphization is not automatically dissolution or release of the contained radionuclide.

Consider a representative bulk volume $\mathcal{V}$ and let $f_a(t)$ be its amorphous fraction. With alpha-decay event rate density $R$, during $dt$ the expected number of independent alpha decay events is $R\mathcal{V}dt$. Each produces effective amorphization volume per decay $V_{\alpha }$, but only the fraction $1-f_a$ of a randomly placed damage region was previously crystalline. The expected new amorphous volume is thereforeDividing by $\mathcal{V}dt$ givesSeparate variables and integrate:Hence, with cumulative alpha-decay event density $D=Rt$,The product $Rt{V}_{\alpha }$ is dimensionless: events per cubic metre accumulated over time multiply cubic metres damaged per event. Here “dose” is an event density, not an absorbed energy per kilogram.

An independent probability interpretation reaches the same result. A particular point is covered by a number of random damage events with a Poisson distribution of mean $D{V}_{\alpha }$. Its probability of remaining crystalline is the zero-hit probability $e^{-D{V}_{\alpha }}$; the complementary fraction has been hit at least once. Overlapping events do not create a second copy of already amorphous material.

At small dose, $f_a\simeq D{V}_{\alpha }$; at larger dose, overlap bends the curve towards one. For example, $D{V}_{\alpha }=1$ gives $f_a=0.6321$, not full radiation-induced amorphization, and $f_a=0.99$ requires $D{V}_{\alpha }=\ln 100=4.6052$.

The derivation assumes a homogeneous, initially crystalline material; independent random event locations; one-event radiation-induced amorphization of a fixed effective volume; and no recovery of that volume. It is the direct-impact model, not a model requiring several prior hits to cross a damage threshold. If $R$ varies while $V_{\alpha }$ stays fixed, replace $Rt$ by ${\int }_0^{t}R(s)ds$. If annealing, material evolution or spatially correlated actinide distributions matter, this simple law need not hold. The subscript on $V_{\alpha }$ indexes an alpha decay event: it does not require the emitted alpha particle itself to produce the damaged volume.

The assumption is a leading approximation for concentrated ballistic radiation damage, not an exact statement that the alpha particle displaces no atoms. The daughter has only about $m_{\alpha }/M_d$ times the alpha particle's kinetic energy, but it is much slower and loses a much larger fraction of that energy in collisions with target nuclei. Those collisions initiate closely packed recoil cascades. The alpha particle deposits most of its much larger energy into the electronic system over a path hundreds of times longer, with a less concentrated displacement contribution.

Thus the dense region susceptible to one-event loss of crystalline order can reasonably be associated chiefly with the daughter recoil. This explains why the effective damaged volume in the direct-impact amorphization model can be estimated or fitted from recoil-dominated cascades. It is the spatial density and stopping channel of the damage, not simply which particle carries more energy, that matters.

The approximation has a measurable limit. Turcotte's comparison of actinide-dioxide self-damage and external alpha irradiation found a nonzero alpha contribution, reporting about $19\%$ of total damage in the studied case. Comparison of diffraction-derived lattice changes with internal-reference measurements helps separate real radiation-induced swelling from instrument shifts, while the distinction between self-decay and external-alpha irradiation probes the omitted recoil channel.^[ref]

Consequently “all ballistic radiation damage” is useful as the stated model assumption but should not be promoted to an exact material-independent law. electronic excitation can also couple to structural change, and recovery alters what survives. The derivation in part (b) remains valid for an effective one-event volume fitted consistently to the relevant regime; it does not establish that all ceramics or spent fuels share the same $V_{\alpha }$ or the same eventual amorphous fraction.