In nuclear reactors, delayed neutron precursors (DNPs) are important for reactor safety and operation. In liquid nuclear fuels, DNPs are transported by the flow, and an advection-reaction balance equation for their concentration must be solved in addition to the neutron balance equation. This research paper applies the method of characteristics to solve the DNP equation using techniques previously developed in petroleum engineering and groundwater analysis. The calculation strategy incorporates a pathline generation algorithm covering all the mesh cells of the geometry for concentration calculations, and pathlines are reconstructed for both step and piecewise linear velocity fields.

The analytical solutions are used on a Cartesian mesh to compute the DNP concentrations, assuming a step source of DNPs. The results obtained on the steady-state phases of a benchmark problem illustrate this new method’s capability to solve advection-dominated problems in liquid-fueled reactors. The method developed in this work is suited to smooth velocity fields (e.g. laminar or Reynolds-average Navier-Stokes flows) where pathlines can be tracked and where DNP diffusivity is negligible.