A general formalism to solve the G-group neutron diffusion equation is described. The G-group flux is represented by complementing an ″asymptotic″ mode with (g-1) ″transient″ modes. A particular reduction-to-one-group technique gives a high computational efficiency. MUSIC, a 2-group code g the above formalism, is presented. MUSIC is demonstrated on a fine-mesh calculation and on 2 coarse-mesh core calculations: a heavy-water reactor (HWR) problem and the 2-D lightwater reactor (LWR) IAEA benchmark. Comparison is made to finite-difference results.
Authors
Bonalumi RA; Rouben B; Dastur AR; Dondale CS; Li HYH