Flux stabilization in neutron problems with fixed sources
Although critical core calculations are the most common in design and safety analysis, fixed source calculations are needed for specific applications, e.g. to compute ex-core detector response functions, to develop new methodologies for dilution and reload error accidents and more in general for all situations involving sub-critical shut-down states. It is well known that the source problem becomes difficult to be solved with core configuration close to criticality, i.e. with the multiplication factor approaching unity, for the occurrence of numerical ill-conditioning and very high number of iterations, possibly leading to failure in the flux convergence. In this work, the Wielandt Eigen-shift technique used in iterative methods of critical problems is developed for source problems too, in order to stabilize the solution. The mathematical basis and the proof of the convergence are discussed. Compared to the existing methods, this technique allows also more control to avoid singular behavior at inner iterations. Numerical tests with a 1D analytical benchmark are reported to prove the robustness of the technique. (authors)