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

Title: High-order local maximum principle preserving (MPP) discontinuous Galerkin finite element method for the transport equation

Authors:
; ; ; ; ; ;
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
1396716
Grant/Contract Number:
LLNL-JRNL-684083
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 334; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-10-04 15:18:03; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English

Citation Formats

Anderson, R., Dobrev, V., Kolev, Tz., Kuzmin, D., Quezada de Luna, M., Rieben, R., and Tomov, V. High-order local maximum principle preserving (MPP) discontinuous Galerkin finite element method for the transport equation. United States: N. p., 2017. Web. doi:10.1016/j.jcp.2016.12.031.
Anderson, R., Dobrev, V., Kolev, Tz., Kuzmin, D., Quezada de Luna, M., Rieben, R., & Tomov, V. High-order local maximum principle preserving (MPP) discontinuous Galerkin finite element method for the transport equation. United States. doi:10.1016/j.jcp.2016.12.031.
Anderson, R., Dobrev, V., Kolev, Tz., Kuzmin, D., Quezada de Luna, M., Rieben, R., and Tomov, V. Sat . "High-order local maximum principle preserving (MPP) discontinuous Galerkin finite element method for the transport equation". United States. doi:10.1016/j.jcp.2016.12.031.
@article{osti_1396716,
title = {High-order local maximum principle preserving (MPP) discontinuous Galerkin finite element method for the transport equation},
author = {Anderson, R. and Dobrev, V. and Kolev, Tz. and Kuzmin, D. and Quezada de Luna, M. and Rieben, R. and Tomov, V.},
abstractNote = {},
doi = {10.1016/j.jcp.2016.12.031},
journal = {Journal of Computational Physics},
number = C,
volume = 334,
place = {United States},
year = {Sat Apr 01 00:00:00 EDT 2017},
month = {Sat Apr 01 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.jcp.2016.12.031

Citation Metrics:
Cited by: 3works
Citation information provided by
Web of Science

Save / Share:
  • No abstract prepared.
  • Discontinuous Galerkin finite element methods are used to estimate solutions to the non-scattering 1-D spherical neutron transport equation. Various trial and test spaces are compared in the context of a few sample problems whose exact solution is known. Certain trial spaces avoid unphysical behaviors that seem to plague other methods. Comparisons with diamond differencing and simple corner-balancing are presented to highlight these improvements.
  • We develop 3rd order maximum-principle-satisfying direct discontinuous Galerkin methods [8], [9], [19] and [21] for convection diffusion equations on unstructured triangular mesh. We carefully calculate the normal derivative numerical flux across element edges and prove that, with proper choice of parameter pair (β 01) in the numerical flux formula, the quadratic polynomial solution satisfies strict maximum principle. The polynomial solution is bounded within the given range and third order accuracy is maintained. There is no geometric restriction on the meshes and obtuse triangles are allowed in the partition. As a result, a sequence of numerical examples are carried outmore » to demonstrate the accuracy and capability of the maximum-principle-satisfying limiter.« less
  • Shallow water equations with a non-flat bottom topography have been widely used to model flows in rivers and coastal areas. An important difficulty arising in these simulations is the appearance of dry areas where no water is present, as standard numerical methods may fail in the presence of these areas. These equations also have still water steady state solutions in which the flux gradients are nonzero but exactly balanced by the source term. In this paper we propose a high order discontinuous Galerkin method which can maintain the still water steady state exactly, and at the same time can preservemore » the non-negativity of the water height without loss of mass conservation. A simple positivity-preserving limiter, valid under suitable CFL condition, will be introduced in one dimension and then extended to two dimensions with rectangular meshes. Numerical tests are performed to verify the positivity-preserving property, well-balanced property, high order accuracy, and good resolution for smooth and discontinuous solutions.« less
  • Purpose: The advent of magnetic resonance imaging (MRI) guided radiotherapy systems demands the incorporation of the magnetic field into dose calculation algorithms of treatment planning systems. This is due to the fact that the Lorentz force of the magnetic field perturbs the path of the relativistic electrons, hence altering the dose deposited by them. Building on the previous work, the authors have developed a discontinuous finite element space-angle treatment of the linear Boltzmann transport equation to accurately account for the effects of magnetic fields on radiotherapy doses. Methods: The authors present a detailed description of their new formalism and comparemore » its accuracy to GEANT4 Monte Carlo calculations for magnetic fields parallel and perpendicular to the radiation beam at field strengths of 0.5 and 3 T for an inhomogeneous 3D slab geometry phantom comprising water, bone, and air or lung. The accuracy of the authors’ new formalism was determined using a gamma analysis with a 2%/2 mm criterion. Results: Greater than 98.9% of all points analyzed passed the 2%/2 mm gamma criterion for the field strengths and orientations tested. The authors have benchmarked their new formalism against Monte Carlo in a challenging radiation transport problem with a high density material (bone) directly adjacent to a very low density material (dry air at STP) where the effects of the magnetic field dominate collisions. Conclusions: A discontinuous finite element space-angle approach has been proven to be an accurate method for solving the linear Boltzmann transport equation with magnetic fields for cases relevant to MRI guided radiotherapy. The authors have validated the accuracy of this novel technique against GEANT4, even in cases of strong magnetic field strengths and low density air.« less