Assessment and improvement of the 2D/1D method stability in DeCART

As part of ongoing work with Consortium for Advanced Simulation of Light Water Reactors (CASL), the 2D/1D code, DeCART, has demonstrated some of the advantages of the 2D/1D method with respect to realistic, full-core analysis, particularly over explicit 3D transport methods, which generally have higher memory and computation requirements. The 2D/1D method performs 2D-radial transport sweeps coupled with ID-axial diffusion calculations to provide a full 3D simulation. DeCART employs the 2D method of characteristics for the radial sweeps and ID one-node nodal diffusion for the axial sweeps, coupling the two methods with transverse leakages to ensure a more consistent representation of the transport equation. It has been observed that refinement of the axial plane thickness leads to instabilities in the calculation scheme. This work assesses the sources of these instabilities and the approaches to improve them, especially with respect to negative scattering cross sections and the tightness of the 2D-radial/ID-axial coupling schemes. Fourier analyses show that the existing iteration scheme is not unconditionally stable, suggesting a tighter coupling scheme is required. For this reason 3D-CMFD has been implemented, among other developments, to ensure more stable calculation. A matrix of test cases has been used to assess the convergence, with the primary parameter being the axial plane thickness, which has been refined down to 1 cm. These cases demonstrate the issues observed and how the modification improve the stability. However, it is apparent that more work is necessary to ensure unconditional stability. (authors)