Discontinuous finite element solution of the radiation diffusion equation on arbitrary polygonal meshes and locally adapted quadrilateral grids
Abstract
In this paper, we propose a piecewise linear discontinuous (PWLD) finite element discretization of the diffusion equation for arbitrary polygonal meshes. It is based on the standard diffusion form and uses the symmetric interior penalty technique, which yields a symmetric positive definite linear system matrix. A preconditioned conjugate gradient algorithm is employed to solve the linear system. Piecewise linear approximations also allow a straightforward implementation of local mesh adaptation by allowing unrefined cells to be interpreted as polygons with an increased number of vertices. Several test cases, taken from the literature on the discretization of the radiation diffusion equation, are presented: random, sinusoidal, Shestakov, and Z meshes are used. The last numerical example demonstrates the application of the PWLD discretization to adaptive mesh refinement.
 Authors:
 Publication Date:
 OSTI Identifier:
 22382158
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 280; Other Information: Copyright (c) 2014 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; APPROXIMATIONS; DIFFUSION; DIFFUSION EQUATIONS; FINITE ELEMENT METHOD; IMPLEMENTATION; MATRICES; RANDOMNESS; SYMMETRY
Citation Formats
Ragusa, Jean C., Email: jean.ragusa@tamu.edu. Discontinuous finite element solution of the radiation diffusion equation on arbitrary polygonal meshes and locally adapted quadrilateral grids. United States: N. p., 2015.
Web. doi:10.1016/J.JCP.2014.09.013.
Ragusa, Jean C., Email: jean.ragusa@tamu.edu. Discontinuous finite element solution of the radiation diffusion equation on arbitrary polygonal meshes and locally adapted quadrilateral grids. United States. doi:10.1016/J.JCP.2014.09.013.
Ragusa, Jean C., Email: jean.ragusa@tamu.edu. Thu .
"Discontinuous finite element solution of the radiation diffusion equation on arbitrary polygonal meshes and locally adapted quadrilateral grids". United States.
doi:10.1016/J.JCP.2014.09.013.
@article{osti_22382158,
title = {Discontinuous finite element solution of the radiation diffusion equation on arbitrary polygonal meshes and locally adapted quadrilateral grids},
author = {Ragusa, Jean C., Email: jean.ragusa@tamu.edu},
abstractNote = {In this paper, we propose a piecewise linear discontinuous (PWLD) finite element discretization of the diffusion equation for arbitrary polygonal meshes. It is based on the standard diffusion form and uses the symmetric interior penalty technique, which yields a symmetric positive definite linear system matrix. A preconditioned conjugate gradient algorithm is employed to solve the linear system. Piecewise linear approximations also allow a straightforward implementation of local mesh adaptation by allowing unrefined cells to be interpreted as polygons with an increased number of vertices. Several test cases, taken from the literature on the discretization of the radiation diffusion equation, are presented: random, sinusoidal, Shestakov, and Z meshes are used. The last numerical example demonstrates the application of the PWLD discretization to adaptive mesh refinement.},
doi = {10.1016/J.JCP.2014.09.013},
journal = {Journal of Computational Physics},
number = ,
volume = 280,
place = {United States},
year = {Thu Jan 01 00:00:00 EST 2015},
month = {Thu Jan 01 00:00:00 EST 2015}
}

Diffusion synthetic acceleration (DSA) schemes compatible with adaptive mesh refinement (AMR) grids are derived for the SN transport equations discretized using highorder discontinuous finite elements. These schemes are directly obtained from the discretized transport equations by assuming a linear dependence in angle of the angular flux along with an exact Fick's law and, therefore, are categorized as partially consistent. These schemes are akin to the symmetric interior penalty technique applied to elliptic problems and are all based on a secondorder discontinuous finite element discretization of a diffusion equation (as opposed to a mixed or P1 formulation). Therefore, they only havemore »

Discontinuous diffusion synthetic acceleration for S{sub n} transport on 2D arbitrary polygonal meshes
In this paper, a Diffusion Synthetic Acceleration (DSA) technique applied to the S{sub n} radiation transport equation is developed using PieceWise Linear Discontinuous (PWLD) finite elements on arbitrary polygonal grids. The discretization of the DSA equations employs an Interior Penalty technique, as is classically done for the stabilization of the diffusion equation using discontinuous finite element approximations. The penalty method yields a system of linear equations that is Symmetric Positive Definite (SPD). Thus, solution techniques such as Preconditioned Conjugate Gradient (PCG) can be effectively employed. Algebraic MultiGrid (AMG) and Symmetric Gauss–Seidel (SGS) are employed as conjugate gradient preconditioners for themore » 
A piecewise linear finite element discretization of the diffusion equation for arbitrary polyhedral grids
We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multidimensional radiation diffusion equation. It uses recently introduced piecewise linear weight and basis functions in the finite element approximation and it can be applied on arbitrary polygonal (2D) or polyhedral (3D) grids. We first demonstrate some analytical properties of the PWL method and perform a simple mode analysis to compare the PWL method with Palmer's vertexcentered finitevolume method and with a bilinear continuous finite element method. We then show that this new PWL method gives solutions comparable to those from Palmer's. However, since the PWL method producesmore » 
Discontinuous finite element method for solving the twodimensional neutron transport equation on Lagrangian meshes
It is often desirable to solve the twodimensional multigroup transport equation for (rz) geometries directly given by hydrodynamic calculations. Usually, only Monte Carlo codes can compute kappa or k eigenvalues on such geometries. Most deterministic codes use an orthogonal mesh or restrict the mesh to a regular triangular grid. Other methods were developed at Commissariat a l'Energie Atomique and Los Alamos National Laboratory but do not solve the problem of sliding between two Lagrangian blocks. Thus, a production code has developed that solves these problems and is able to obtain kappa or k eigenvalues with a good degree of accuracymore »