DOE OSTI · 1325480
MOC Efficiency Improvements Using a Jacobi Inscatter Approximation
Abstract
In recent weeks, attention has been given to resolving the convergence issues encountered with TCP 0 by trying a Jacobi (J) inscatter approach when group sweeping, where the inscatter source is constructed using the previous iteration flux. This is in contrast to a Gauss-Seidel (GS) approach, which has been the default to-date, where the scattering source uses the most up-to-date flux values. The former is consistent with CASMO, which has no issues with TCP 0 convergence. Testing this out on a variety of problems has demonstrated that the Jacobi approach does indeed provide substantially more stability, though can take more outer iterations to converge. While this is not surprising, there are improvements that can be made to the MOC sweeper to capitalize on the Jacobi approximation and provide substantial speedup. For example, the loop over groups, which has traditionally been the outermost loop in MPACT, can be moved to the interior, avoiding duplicate modular ray trace and coarse ray trace setup (mapping coarse mesh surface indexes), which needs to be performed repeatedly when group is outermost.
Explore related subjects
Keep this discovery
Explore connections, maps & timelines
Stimpson, Shane, Collins, Benjamin, Kochunas, Brendan. 2016-08-31. MOC Efficiency Improvements Using a Jacobi Inscatter Approximation. https://doi.org/10.2172/1325480
Cite the original work for its findings. Save a collection to share your selection of sources.