A stable 1D multigroup high-order low-order method

Yee, Ben Chung; Wollaber, Allan Benton; Haut, Terry Scot; Park, HyeongKae

doi:10.1080/23324309.2016.1187172

Title: A stable 1D multigroup high-order low-order method

Full Record
Other Related Research

Abstract

The high-order low-order (HOLO) method is a recently developed moment-based acceleration scheme for solving time-dependent thermal radiative transfer problems, and has been shown to exhibit orders of magnitude speedups over traditional time-stepping schemes. However, a linear stability analysis by Haut et al. (2015 Haut, T. S., Lowrie, R. B., Park, H., Rauenzahn, R. M., Wollaber, A. B. (2015). A linear stability analysis of the multigroup High-Order Low-Order (HOLO) method. In Proceedings of the Joint International Conference on Mathematics and Computation (M&C), Supercomputing in Nuclear Applications (SNA) and the Monte Carlo (MC) Method; Nashville, TN, April 19–23, 2015. American Nuclear Society.) revealed that the current formulation of the multigroup HOLO method was unstable in certain parameter regions. Since then, we have replaced the intensity-weighted opacity in the first angular moment equation of the low-order (LO) system with the Rosseland opacity. Furthermore, this results in a modified HOLO method (HOLO-R) that is significantly more stable.

Authors:

Yee, Ben Chung ^[1]; Wollaber, Allan Benton ^[2]; Haut, Terry Scot ^[2]; Park, HyeongKae ^[2]

Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Univ. of Michigan, Ann Arbor, MI (United States)
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

Publication Date:: 2016-07-13

Research Org.:: Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

Sponsoring Org.:: USDOE

OSTI Identifier:: 1296672

Report Number(s):: LA-UR-15-28071
Journal ID: ISSN 2332-4309

Grant/Contract Number:: AC52-06NA25396

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational and Theoretical Transport

Additional Journal Information:: Journal Name: Journal of Computational and Theoretical Transport; Journal ID: ISSN 2332-4309

Publisher:: Taylor and Francis

Country of Publication:: United States

Language:: English

Subject:: 97 MATHEMATICS AND COMPUTING; thermal transfer, stability analysis, discontinuous Galerkin

Citation Formats


                    Yee, Ben Chung, Wollaber, Allan Benton, Haut, Terry Scot, and Park, HyeongKae. A stable 1D multigroup high-order low-order method.  United States: N. p., 2016. 
Web.  doi:10.1080/23324309.2016.1187172.

Copy to clipboard


                    Yee, Ben Chung, Wollaber, Allan Benton, Haut, Terry Scot, & Park, HyeongKae. A stable 1D multigroup high-order low-order method.  United States.  https://doi.org/10.1080/23324309.2016.1187172

Copy to clipboard


                    Yee, Ben Chung, Wollaber, Allan Benton, Haut, Terry Scot, and Park, HyeongKae. Wed .  
"A stable 1D multigroup high-order low-order method".  United States.  https://doi.org/10.1080/23324309.2016.1187172.  https://www.osti.gov/servlets/purl/1296672.

Copy to clipboard


                    
@article{osti_1296672,

  title        = {A stable 1D multigroup high-order low-order method},

  author       = {Yee, Ben Chung and Wollaber, Allan Benton and Haut, Terry Scot and Park, HyeongKae},

  abstractNote = {The high-order low-order (HOLO) method is a recently developed moment-based acceleration scheme for solving time-dependent thermal radiative transfer problems, and has been shown to exhibit orders of magnitude speedups over traditional time-stepping schemes. However, a linear stability analysis by Haut et al. (2015 Haut, T. S., Lowrie, R. B., Park, H., Rauenzahn, R. M., Wollaber, A. B. (2015). A linear stability analysis of the multigroup High-Order Low-Order (HOLO) method. In Proceedings of the Joint International Conference on Mathematics and Computation (M&C), Supercomputing in Nuclear Applications (SNA) and the Monte Carlo (MC) Method; Nashville, TN, April 19–23, 2015. American Nuclear Society.) revealed that the current formulation of the multigroup HOLO method was unstable in certain parameter regions. Since then, we have replaced the intensity-weighted opacity in the first angular moment equation of the low-order (LO) system with the Rosseland opacity. Furthermore, this results in a modified HOLO method (HOLO-R) that is significantly more stable.},

  doi          = {10.1080/23324309.2016.1187172},

  journal      = {Journal of Computational and Theoretical Transport},

  number       = ,

  volume       = ,

  place        = {United States},

  year         = {2016},

  month        = {7}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1080/23324309.2016.1187172

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 1 work

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

Fourier Analysis of the Multilevel QD Method with {alpha}-Approximation For Time-Dependent Radiative Transfer Problems

Journal Article Anistratov, Dmitriy Y. - Transactions of the American Nuclear Society

Thermal radiation transfer (TRT) problems describe interaction of photon radiation with matter. They are defined by the time-dependent radiative transfer (RT) equation for the specific intensity I coupled with the energy balance (EB) equation. This class of problems is characterised by high dimensionality, multiple scales and strong nonlinearity. In this paper we analyze the multilevel quasi-diffusion (QD) method for solving TRT problems. This is a nonlinear method of moments with exact closures. The QD method enables one to reduce dimensionality of the original radiative transfer problem. It is formulated by a nonlinear system of equations consisting of (i) the high-ordermore »« less
Space-Dependent Wielandt Shift Methods for Multigroup Diffusion Eigenvalue Problems

Journal Article Yee, Ben C. ; Larsen, Edward W. ; Kochunas, Brendan ; ... - Transactions of the American Nuclear Society

In the calculation of neutronic eigenfunctions and eigenvalues, it is often desirable to approximate the neutron transport equation by a diffusion equation that requires significantly less computational resources to solve. Even when a full transport solution is needed, the solution of the diffusion system can be used to accelerate iterative methods for calculating the transport eigenfunction and eigenvalue. The convergence rate of the basic power iteration method for eigenvalue problems is governed by the dominance ratio; realistic reactor core problems often have high dominance ratios and the power iteration method can converge slowly as a result. A common technique formore »« less
Multigroup Radiation-Hydrodynamics with a High-Order, Low-Order Method

Journal Article Wollaber, Allan Benton ; Park, HyeongKae ; Lowrie, Robert Byron ; ... - Nuclear Science and Engineering

Recent efforts at Los Alamos National Laboratory to develop a moment-based, scale-bridging [or high-order (HO)–low-order (LO)] algorithm for solving large varieties of the transport (kinetic) systems have shown promising results. A part of our ongoing effort is incorporating this methodology into the framework of the Eulerian Applications Project to achieve algorithmic acceleration of radiationhydrodynamics simulations in production software. By starting from the thermal radiative transfer equations with a simple material-motion correction, we derive a discretely consistent energy balance equation (LO equation). We demonstrate that the corresponding LO system for the Monte Carlo HO solver is closely related to the originalmore »« less
Cited by 3
https://doi.org/10.13182/NSE16-45

Full Text Available
A Transport Acceleration Scheme for Multigroup Discrete Ordinates with Upscattering

Journal Article Evans, Thomas M ; Clarno, Kevin T ; Morel, Jim E. - Nuclear Science and Engineering

We have developed a modification of the two-grid upscatter acceleration scheme of Adams and Morel. The modified scheme uses a low-angular-order discrete ordinates equation to accelerate Gauss-Seidel multigroup iteration. This modification ensures that the scheme does not suffer from consistency problems that can affect diffusion-accelerated methods in multidimensional, multimaterial problems. The new transport two-grid scheme is very simple to implement for different spatial discretizations because it uses the same transport operator. The scheme has also been demonstrated to be very effective on three-dimensional, multimaterial problems. On simple one-dimensional graphite and heavy-water slabs modeled in three dimensions with reflecting boundary conditions,more »« less
https://doi.org/10.13182/NSE09-28
Space-Dependent Wielandt Shifts for Multigroup Diffusion Eigenvalue Problems

Journal Article Yee, Ben C. ; Kochunas, Brendan ; Larsen, Edward W. ; ... - Nuclear Science and Engineering

We present a new concept—the space-dependent Wielandt shift (SDWS)—for accelerating the convergence of the power iteration (PI) scheme for multigroup diffusion k-eigenvalue problems. The SDWS improves on standard Wielandt shift (WS) techniques, which are empirical in nature and are typically effective only when the current estimate of the solution is reasonably converged. By accounting for the physics of the problem through SDWS, we are able to improve the acceleration for the initial iterates when the current estimate of the solution is not close to convergence. Numerical results from one-dimensional problems suggest that, compared to standard WS techniques, the new SDWSmore »« less
https://doi.org/10.1080/00295639.2017.1350001

Similar Records