On the non-uniqueness of the nodal mathematical adjoint

• We evaluate three CMFD schemes for computing the nodal mathematical adjoint.
• The nodal mathematical adjoint is not unique and can be non-positive (nonphysical).
• Adjoint and forward eigenmodes are compatible if produced by the same CMFD method.
• In nodal applications the excited eigenmodes are purely mathematical entities.

Computation of the neutron adjoint flux within the framework of modern nodal diffusion methods is often facilitated by reducing the nodal equation system for the forward flux into a simpler coarse-mesh finite-difference form and then transposing the resultant matrix equations. The solution to the transposed problem is known as the nodal mathematical adjoint. Since the coarse-mesh finite-difference reduction of a given nodal formulation can be obtained in a number of ways, different nodal mathematical adjoint solutions can be computed. This non-uniqueness of the nodal mathematical adjoint challenges the credibility of the reduction strategy and demands a verdict as to its suitability in practical applications. This is the matter under consideration in this paper. A selected number of coarse-mesh finite-difference reduction schemes are described and compared. Numerical calculations are utilised to illustrate the differences in the adjoint solutions as well as to appraise the impact on such common applications as the computation of core point kinetics parameters. Recommendations are made for the proper application of the coarse-mesh finite-difference reduction approach to the nodal mathematical adjoint problem.