A ROBUST SECOND-ORDER MULTIPLE BALANCE METHOD FOR TIME-DEPENDENT NEUTRON TRANSPORT SIMULATIONS

Abstract

A second-order “Time-Dependent Multiple Balance” (TDMB) method for solving neutron transport problems is introduced and investigated. TDMB consists of solving two coupled equations: (i) the original balance equation (the transport equation integrated over a time step) and (ii) the “balance-like” auxiliary equation (an approximate neutron balance equation). Simple analysis shows that TDMB is second-order accurate and robust (unconditionally free from spurious oscillation). A source iteration (SI) method with diffusion synthetic acceleration (DSA) is formulated to solve these equations. A Fourier analysis reveals that the convergence rates of the proposed iteration schemes for TDMB are similar to those of the common (SI + DSA) schemes for Backward Euler (BE); however, TDMB requires about twice the computational effort per iteration. To demonstrate the theory—accuracy, robustness, and convergence rate—and investigate the efficiency of TDMB, we present results from a discrete ordinates (Sn) research code. Results are discussed, and future work is proposed.

Keywords

• time-dependent neutron transport
• time discretization method
• multiple balance
• diffusion synthetic acceleration (dsa)