Skip to content
M. S. M. HusseinReactor & Radiation PhysicistRev. 2026-10-01

§ Research · doctoral design study

A CANDU that burns warheads, and cannot run on half of itself.

The Multispectrum CANDU Reactor puts a small helium-cooled fast-neutron core in the middle of a CANDU 6 thermal core. The fast core destroys surplus weapons-grade uranium or plutonium; the thermal core makes power from natural uranium. Neither core is critical on its own.

Thesis
RMC, May 2017, 363 pp.
Supervisors
Hugues Bonin · Brent Lewis
Codes
Serpent · MCNP5/6 · WIMS-AECL
Funding
NSERC · UNENE · HPCVL

§1 Problem

Disarmament left a stockpile that has to go somewhere.

Post-Cold-War treaties dismantled thousands of warheads and released large quantities of highly enriched fissile material. Storing it is a proliferation risk; burning it in reactors turns it into power and into isotopes that are far harder to weaponise.

Thermal reactors burn it slowly and breed new plutonium as they go. Fast reactors burn it efficiently but Canada does not have one. The question the thesis asks: can the reactor Canada already knows how to build be given a fast core?

§2 Concept

Two cores, one calandria.

Select a model. Geometry is to scale from the thesis tables.

Cross-section of the Multispectrum CANDU Reactor, Model IA circular CANDU 6 core of 320 thermal fuel channels on a 28.575 cm square lattice surrounding an inner fast-spectrum core of 32 channels within a 108.86 cm radius, separated by a partition wall.1 m
histories 0collisions 0in fast core 0 %
Fast-core radius
108.9 cm
Thermal channels
320
Fast channels
32
Fast-core pitch
28.6 cm
Fast-core fuel
U-235 · 19.9 %

Same pitch as the thermal core. Reached subcriticality after 135 days.

  • Thermal core · nat. UO₂ in D₂O
  • Fast core · U or Pu MOX in He
  • Partition wall · SS316
  • D₂O reflector
  • Neutron history · schematic
Multispectrum CANDU Reactor, six models from the PhD thesis (RMC, 2017), Table 4-1. Geometry to scale: CANDU 6 lattice pitch 28.575 cm, core radius 335 cm. Colours indicate the neutron-spectrum region, not flux. Fast-core channel positions are inferred from the tabulated counts. The tracks are an illustrative random walk with a longer mean free path in the helium-cooled fast core, not a transport calculation.
43.3 mm1 + 6 + 12 + 18 = 37
The standard CANDU 37-element bundle, used in both cores. Thermal core: Zircaloy-4 sheath, heavy-water coolant. Fast core: stainless steel 316L(N) sheath, helium coolant. Drawn from published AECL/IAEA bundle geometry; the thesis keeps standard dimensions.

Why the safety argument holds

Each core is designed subcritical alone. Only the neutrons each sends the other push the combined system to keff = 1. Drain the moderator, void the coolant, or disturb either core and the coupling breaks and the whole reactor shuts down. That is the defining property of a coupled reactor, and it is why the design starts from coupled-reactor theory rather than from a single-core model.

Numbers that fix the design

  • CANDU 6 reference: 380 channels on a 28.575 cm pitch, core radius 335 cm, 594 cm channels, 12 bundles per channel.
  • Three fast-core radii: 137.32, 108.86 and 80.39 cm, measured from the partition wall.
  • Fast-core pitch 28.575 cm (same as CANDU) or 14.3–14.6 cm (twice the channel density).
  • Fast fuel: 19.9 % enriched uranium, or MOX with 13.89 % Pu-239 in depleted uranium.

§3 Method

Trust nothing until two codes agree.

Three stages, each a paper before the thesis.

  1. 01

    Verify the theory

    Two-point coupled-reactor theory was checked numerically on the Japanese Deuterium Critical Assembly, a two-region core whose inner and outer water levels can be varied. MCNP5 and Serpent gave the same multiplication factors and coupling coefficients as the analytical model.

    CANDU Fuel 2013 ↗

  2. 02

    Validate the tools

    Full-core CANDU 6 (Gentilly-2) models in MCNP6 and Serpent were run through burnup and compared with WIMS-AECL 3.1.2.1, the lattice code the Canadian industry uses. k-effective against burnup agreed; Serpent was the faster of the two.

    PBNC 2014 ↗

  3. 03

    Design the reactor

    With the codes trusted, six MSCR models were built in Serpent: fast-core radius, lattice pitch, number of fast channels and fuel (19.9 % U-235 or 13.9 % Pu-239 MOX) were varied and criticality, flux, power, form factor, burnup and actinide inventories were tracked over years of operation.

    ICMSNSE 2015 ↗

Eq. Two-point coupled-reactor multiplication factor

keff=k11+k222+(k11−k222)2+k12 k21k_{\text{eff}} = \frac{k_{11}+k_{22}}{2} + \sqrt{\left(\frac{k_{11}-k_{22}}{2}\right)^{2} + k_{12}\,k_{21}}
kij=fission neutrons produced in core i by neutrons born in core jneutrons born in core jk_{ij} = \frac{\text{fission neutrons produced in core } i \text{ by neutrons born in core } j}{\text{neutrons born in core } j}

k11, k22 are the isolated cores' factors; k12, k21 the coupling coefficients, computed from track-length flux tallies. Thesis Ch. 5, verified on the DCA with MCNP5 and Serpent.

§4 Result

What each model did.

From the thesis abstract and Chapters 7–8. Models without a headline number are marked as such.

Six MSCR models
ModelFast radiusChannels th / fastFast pitchFast fuelOutcome
I108.86 cm320 / 3228.575 cmU-235 · 19.9 %Same pitch as the thermal core. Reached subcriticality after 135 days.
II108.86 cm320 / 14814.575 cmU-235 · 19.9 %Destroyed 10.7 % of fissile isotopes over 1000 days without refuelling.
III137.32 cm292 / 5228.575 cmPu-239 MOX · 13.9 %Largest fast core at the thermal pitch. Detailed results in Chapter 8.
IV137.32 cm292 / 24014.287 cmPu-239 MOX · 13.9 %Highest fissile destruction of the plutonium models: 15.3 % over 2600 days.
V108.86 cm320 / 14414.287 cmPu-239 MOX · 13.9 %Mid-size plutonium core. Detailed results in Chapter 8.
VI80.39 cm348 / 7614.287 cmPu-239 MOX · 13.9 %Smallest fast core, highest consumption rate: 420 days to subcriticality.

The thesis concludes that destruction of fissile isotopes from dismantled weapons can be achieved in an MSCR, with Model II (uranium) and Model IV (plutonium) the strongest performers, and flags that Models II, IV, V and VI carry form factors above the desired safety margin, so power flattening is the next design task.

§5 Sources

Read the originals.