A flexible nonlinear diffusion acceleration method for the SN transport equations discretized with discontinuous finite elements

Schunert, Sebastian; Wang, Yaqi; Gleicher, Frederick; Ortensi, Javier; Baker, Benjamin; Laboure, Vincent; Wang, Congjian; DeHart, Mark; Martineau, Richard

doi:10.1016/j.jcp.2017.01.070

Title: A flexible nonlinear diffusion acceleration method for the S_N transport equations discretized with discontinuous finite elements

Full Record
Other Related Research

Abstract

This paper presents a flexible nonlinear diffusion acceleration (NDA) method that discretizes both the S_N transport equation and the diffusion equation using the discontinuous finite element method (DFEM). The method is flexible in that the diffusion equation can be discretized on a coarser mesh with the only restriction that it is nested within the transport mesh and the FEM shape function orders of the two equations can be different. The consistency of the transport and diffusion solutions at convergence is defined by using a projection operator mapping the transport into the diffusion FEM space. The diffusion weak form is based on the modified incomplete interior penalty (MIP) diffusion DFEM discretization that is extended by volumetric drift, interior face, and boundary closure terms. In contrast to commonly used coarse mesh finite difference (CMFD) methods, the presented NDA method uses a full FEM discretized diffusion equation for acceleration. Suitable projection and prolongation operators arise naturally from the FEM framework. Via Fourier analysis and numerical experiments for a one-group, fixed source problem the following properties of the NDA method are established for structured quadrilateral meshes: (1) the presented method is unconditionally stable and effective in the presence of mild material heterogeneities if themore »« less

Authors:

Schunert, Sebastian ^[1]; Wang, Yaqi ^[1]; Gleicher, Frederick ^[1]; Ortensi, Javier ^[1]; Baker, Benjamin ^[1]; Laboure, Vincent ^[1]; Wang, Congjian ^[1]; DeHart, Mark ^[1]; Martineau, Richard ^[1]

Idaho National Lab. (INL), Idaho Falls, ID (United States)

Publication Date:: Tue Feb 21 00:00:00 EST 2017

Research Org.:: Idaho National Lab. (INL), Idaho Falls, ID (United States)

Sponsoring Org.:: USDOE Office of Nuclear Energy (NE)

OSTI Identifier:: 1363728

Alternate Identifier(s):: OSTI ID: 1397847

Report Number(s):: INL/JOU-16-37745
Journal ID: ISSN 0021-9991; PII: S0021999117301286

Grant/Contract Number:: AC07-05ID14517

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 338; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 73 NUCLEAR PHYSICS AND RADIATION PHYSICS; neutron transport equation; nonlinear diffusion acceleration; discontinuous finite element method

Citation Formats


                    Schunert, Sebastian, Wang, Yaqi, Gleicher, Frederick, Ortensi, Javier, Baker, Benjamin, Laboure, Vincent, Wang, Congjian, DeHart, Mark, and Martineau, Richard. A flexible nonlinear diffusion acceleration method for the SN transport equations discretized with discontinuous finite elements.  United States: N. p., 2017. 
Web.  doi:10.1016/j.jcp.2017.01.070.

Copy to clipboard


                    Schunert, Sebastian, Wang, Yaqi, Gleicher, Frederick, Ortensi, Javier, Baker, Benjamin, Laboure, Vincent, Wang, Congjian, DeHart, Mark, & Martineau, Richard. A flexible nonlinear diffusion acceleration method for the SN transport equations discretized with discontinuous finite elements.  United States.  https://doi.org/10.1016/j.jcp.2017.01.070

Copy to clipboard


                    Schunert, Sebastian, Wang, Yaqi, Gleicher, Frederick, Ortensi, Javier, Baker, Benjamin, Laboure, Vincent, Wang, Congjian, DeHart, Mark, and Martineau, Richard. Tue .  
"A flexible nonlinear diffusion acceleration method for the SN transport equations discretized with discontinuous finite elements".  United States.  https://doi.org/10.1016/j.jcp.2017.01.070.  https://www.osti.gov/servlets/purl/1363728.

Copy to clipboard


                    
@article{osti_1363728,

  title        = {A flexible nonlinear diffusion acceleration method for the SN transport equations discretized with discontinuous finite elements},

  author       = {Schunert, Sebastian and Wang, Yaqi and Gleicher, Frederick and Ortensi, Javier and Baker, Benjamin and Laboure, Vincent and Wang, Congjian and DeHart, Mark and Martineau, Richard},

  abstractNote = {This paper presents a flexible nonlinear diffusion acceleration (NDA) method that discretizes both the SN transport equation and the diffusion equation using the discontinuous finite element method (DFEM). The method is flexible in that the diffusion equation can be discretized on a coarser mesh with the only restriction that it is nested within the transport mesh and the FEM shape function orders of the two equations can be different. The consistency of the transport and diffusion solutions at convergence is defined by using a projection operator mapping the transport into the diffusion FEM space. The diffusion weak form is based on the modified incomplete interior penalty (MIP) diffusion DFEM discretization that is extended by volumetric drift, interior face, and boundary closure terms. In contrast to commonly used coarse mesh finite difference (CMFD) methods, the presented NDA method uses a full FEM discretized diffusion equation for acceleration. Suitable projection and prolongation operators arise naturally from the FEM framework. Via Fourier analysis and numerical experiments for a one-group, fixed source problem the following properties of the NDA method are established for structured quadrilateral meshes: (1) the presented method is unconditionally stable and effective in the presence of mild material heterogeneities if the same mesh and identical shape functions either of the bilinear or biquadratic type are used, (2) the NDA method remains unconditionally stable in the presence of strong heterogeneities, (3) the NDA method with bilinear elements extends the range of effectiveness and stability by a factor of two when compared to CMFD if a coarser diffusion mesh is selected. In addition, the method is tested for solving the C5G7 multigroup, eigenvalue problem using coarse and fine mesh acceleration. Finally, while NDA does not offer an advantage over CMFD for fine mesh acceleration, it reduces the iteration count required for convergence by almost a factor of two in the case of coarse mesh acceleration.},

  doi          = {10.1016/j.jcp.2017.01.070},

  journal      = {Journal of Computational Physics},

  number       = ,

  volume       = 338,

  place        = {United States},

  year         = {Tue Feb 21 00:00:00 EST 2017},

  month        = {Tue Feb 21 00:00:00 EST 2017}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.jcp.2017.01.070

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 12 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

Crystal plasticity modeling of irradiation effects on flow stress in pure-iron and iron-copper alloys

Journal Article Chakraborty, Pritam ; Bulent Biner, S. - Mechanics of Materials

This work presents a flexible Nonlinear diffusion acceleration (NDA) method that discretizes both themore »« less
https://doi.org/10.1016/j.mechmat.2016.07.013
Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)

Technical Report Wang, Dean ; Downar, Thomas ; Xu, Yunlin ; ...

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to develop and implement advanced acceleration schemes for the DOE NEAMS neutronics code PROTEUSSN. The project team included the University of Massachusetts Lowell, The Ohio State University, University of Michigan, and Argonne National Laboratory. The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transportmore »« less
https://doi.org/10.2172/1511575

Full Text Available
Implementation of CMFD Acceleration Scheme in PROTEUS-MOC

Journal Article Hsieh, Albert ; Yang, Won Sik - Transactions of the American Nuclear Society

The PROTEUS-MOC code is a three-dimensional (3D) transport code, which is being developed at Argonne National Laboratory under the DOE-NE's Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. It provides a fully consistent 3D discretization of the transport equation by combining the two-dimensional (2D) method of characteristics (MOC) for the planar direction with the discontinuous Galerkin method for the treatment of the axial variable. PROTEUS-MOC uses the unstructured finite element method (FEM) and thus allows the explicit representation of any extruded geometries. It also has a capability of massively parallel computation in energy group, angle, plane, and space using amore »« less
A High-Order Nonlinear Diffusion Acceleration for the SN Equations Discretized with the Discontinuous FEM I: Theory and Numerical Results

Conference Schunert, Sebastian ; Wang, Yaqi ; Ortensi, Javier ; ...

Recently, nonlinear diffusion acceleration (NDA) methods have received widespread attention in the radiation transport community because of their inherent efficiency (Willert 2014, Wang 2014, Park 2012). However, NDA is in general only conditionally stable \cite{Larsen2014}, while diffusion synthetic acceleration methods (DSA) are unconditionally stable if suitable discretization methods for the transport and diffusion systems are chosen \cite{Larsen1982}. The presented work introduces a new generalized NDA method for solving the transport equation discretized withmore »« less
Full Text Available
On the Degradation of the Effectiveness of Nonlinear Diffusion Acceleration with Parallel Block Jacobi Splitting

Journal Article Schunert, Sebastian ; Wang, Yaqi ; Ortensi, Javier ; ... - Transactions of the American Nuclear Society

MOOSE's radiation transport module Rattlesnake comprises a first order S N solver that can be accelerated with a flexible Nonlinear Diffusion Acceleration. At each NDA iteration, the streaming and collision operator of the S{sub N} equations are (approximately) inverted using a transport sweep. When executed on multiple processors parallel block Jacobi (PBJ) splitting is utilized lagging the processor sub-domain's incoming angular fluxes making the implemented method asynchronous, i.e. the convergence properties of NDA become dependent on the number of concurrent tasks and details of the partitioning. Within this work, we study the NDA convergence properties for the two-dimensional C5G7 benchmarkmore »« less

Similar Records

Title: A flexible nonlinear diffusion acceleration method for the SN transport equations discretized with discontinuous finite elements

Abstract

Citation Formats

Title: A flexible nonlinear diffusion acceleration method for the S_N transport equations discretized with discontinuous finite elements