Nucleide

Depletion

Load a depletion chain and run one CRAM burnup step in the browser.

Depletion solves the Bateman equation for atom densities n(t)\vec{n}(t):

dndt=An,n(t)=eAtn0\frac{d\vec{n}}{dt} = A \vec{n}, \qquad \vec{n}(t) = e^{A t} \vec{n}_0

where AA is the transmutation matrix built from decay, capture, and fission reaction rates. Nucleide uses the Chebyshev Rational Approximation Method (CRAM) with published 16- and 48-order coefficients.

Depletion chain

A chain XML lists nuclides, half-lives, decay branches, reaction targets, and fission yields. The demo below uses a tiny three-nuclide chain:

I-135βXe-135βCs-135\text{I-135} \xrightarrow{\beta} \text{Xe-135} \xrightarrow{\beta} \text{Cs-135}

One-step solve

deplete(chain, n0, dt, rates, order) integrates over one time step dt. In production you would loop this call, update reaction rates from a transport solve, and accumulate burnup. The demo also offers a Burnup curve run that repeats the step and plots nuclide densities versus cumulative time.

Loading Nucleide WASM…

See also