skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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.

Authors:
 [1];  [1];  [2]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  2. Univ. of Michigan, Ann Arbor, MI (United States)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Consortium for Advanced Simulation of LWRs (CASL)
Sponsoring Org.:
USDOE Office of Nuclear Energy (NE)
OSTI Identifier:
1325480
Report Number(s):
ORNL/TM-2016/406; CASL-U-2016-1056-002
NT0304000; NEAF343; CASL-U-2016-1056-002; TRN: US1700079
DOE Contract Number:
AC05-00OR22725
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS; APPROXIMATIONS; CONVERGENCE; EFFICIENCY; SCATTERING; STABILITY; TESTING; ITERATIVE METHODS; M CODES

Citation Formats

Stimpson, Shane, Collins, Benjamin, and Kochunas, Brendan. MOC Efficiency Improvements Using a Jacobi Inscatter Approximation. United States: N. p., 2016. Web. doi:10.2172/1325480.
Stimpson, Shane, Collins, Benjamin, & Kochunas, Brendan. MOC Efficiency Improvements Using a Jacobi Inscatter Approximation. United States. doi:10.2172/1325480.
Stimpson, Shane, Collins, Benjamin, and Kochunas, Brendan. 2016. "MOC Efficiency Improvements Using a Jacobi Inscatter Approximation". United States. doi:10.2172/1325480. https://www.osti.gov/servlets/purl/1325480.
@article{osti_1325480,
title = {MOC Efficiency Improvements Using a Jacobi Inscatter Approximation},
author = {Stimpson, Shane and Collins, Benjamin and Kochunas, Brendan},
abstractNote = {In recent weeks, attention has been given to resolving the convergence issues encountered with TCP0 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 TCP0 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.},
doi = {10.2172/1325480},
journal = {},
number = ,
volume = ,
place = {United States},
year = 2016,
month = 8
}

Technical Report:

Save / Share:
  • The MPACT code, being developed collaboratively by the University of Michigan and Oak Ridge National Laboratory, is the primary deterministic neutron transport solver being deployed within the Virtual Environment for Reactor Applications (VERA) as part of the Consortium for Advanced Simulation of Light Water Reactors (CASL). In many applications of the MPACT code, transport-corrected scattering has proven to be an obstacle in terms of stability, and considerable effort has been made to try to resolve the convergence issues that arise from it. Most of the convergence problems seem related to the transport-corrected cross sections, particularly when used in the 2Dmore » method of characteristics (MOC) solver, which is the focus of this work. Here in this paper, the stability and performance of the 2-D MOC solver in MPACT is evaluated for two iteration schemes: Gauss-Seidel and Jacobi. With the Gauss-Seidel approach, as the MOC solver loops over groups, it uses the flux solution from the previous group to construct the inscatter source for the next group. Alternatively, the Jacobi approach uses only the fluxes from the previous outer iteration to determine the inscatter source for each group. Consequently for the Jacobi iteration, the loop over groups can be moved from the outermost loop$-$as is the case with the Gauss-Seidel sweeper$-$to the innermost loop, allowing for a substantial increase in efficiency by minimizing the overhead of retrieving segment, region, and surface index information from the ray tracing data. Several test problems are assessed: (1) Babcock & Wilcox 1810 Core I, (2) Dimple S01A-Sq, (3) VERA Progression Problem 5a, and (4) VERA Problem 2a. The Jacobi iteration exhibits better stability than Gauss-Seidel, allowing for converged solutions to be obtained over a much wider range of iteration control parameters. Additionally, the MOC solve time with the Jacobi approach is roughly 2.0-2.5× faster per sweep. While the performance and stability of the Jacobi iteration are substantially improved compared to the Gauss-Seidel iteration, it does yield a roughly 8$-$10% increase in the overall memory requirement.« less
  • Potential energy savings, and effect on performance were determined for four furnaces and two boilers when retrofitted with the following efficiency improvement techniques; (a) derate, (b) vent connector restrictor, (c) automatic vent damper, (d) thermal vent damper, (e) IID's, and (f) combinations of the above. Advance technology systems with high seasonal efficiency were studied (pulse combustion boiler, pulse combustion furnaces, power burner furnaces, power boilers) and the results are presented. Also discussed are the energy saving potential of two heat recovery devices for installation in the flue of gas-fired appliances.
  • A simple model was constructed to permit the calculation of the operating characteristics of a thermionic diode. It was assumed that space charge was eliminated by a partially ionized plasma, the sole effect of which on the diode is to introduce an internal electrical resistance and thermal conductance. The efficiency of conversion of diodes of fixed parameters was then optimized with respect to lead resistance and operating voltage. Account was taken of back current and of the voltage drop, thermal conduction, and Joule heating of the leads and plasma. A critenion for the lifetime of a thermionic converter was developed,more » based upon the evaporation rate of the cathode. From this criterion, a relation between the heat of vaporization of the cathode and the maximum temperature at which it can be operated was derived. From the latter quantity, a means of calculating the ultimate efficiency atiainable with a given cathode mnterial was obtained. (auth)« less