Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations
Abstract
In this work, we introduce a new residual distribution (RD) framework for the design of bound-preserving high-resolution finite element schemes. The continuous and discontinuous Galerkin discretizations of the linear advection equation are modified to construct local extremum diminishing (LED) approximations. To that end, we perform mass lumping and redistribute the element residuals in a manner which guarantees the LED property. The hierarchical correction procedure for high-order Bernstein finite element discretizations involves localization to subcells and definition of bound-preserving weights for subcell contributions. Using strong stability preserving (SSP) Runge–Kutta methods for time integration, we prove the validity of discrete maximum principles under CFL-like time step restrictions. The low-order version of our method has roughly the same accuracy as the one derived from a piecewise (multi)-linear approximation on a submesh with the same nodal points. In high-order extensions, we use an element-based flux-corrected transport (FCT) algorithm which can be interpreted as a nonlinear RD scheme. The proposed LED corrections are tailor-made for matrix-free implementations which avoid the rapidly growing cost of matrix assembly for high-order Bernstein elements. The results for 1D, 2D, and 3D test problems compare favorably to those obtained with the best matrix-based approaches.
- Authors:
-
- Technische Univ. of Dortmund (Germany). Inst. of Applied Mathematics
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Center for Applied Scientific Computing
- Univ. Zurich (Switzerland). Inst. of Computational Sciences
- Publication Date:
- Research Org.:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Org.:
- USDOE National Nuclear Security Administration (NNSA); German Research Foundation (DFG); SNF
- OSTI Identifier:
- 1836210
- Alternate Identifier(s):
- OSTI ID: 1580794
- Report Number(s):
- LLNL-JRNL-768125
Journal ID: ISSN 0045-7825; 958809
- Grant/Contract Number:
- AC52-07NA27344; KU 1530/23-1; 200020 175784
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Computer Methods in Applied Mechanics and Engineering
- Additional Journal Information:
- Journal Volume: 359; Journal ID: ISSN 0045-7825
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; advection equation; discrete maximum principles; Bernstein finite elements; matrix-free methods; residual distribution; flux-corrected transport
Citation Formats
Hajduk, Hennes, Kuzmin, Dmitri, Kolev, Tzanio, and Abgrall, Remi. Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations. United States: N. p., 2019.
Web. doi:10.1016/j.cma.2019.112658.
Hajduk, Hennes, Kuzmin, Dmitri, Kolev, Tzanio, & Abgrall, Remi. Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations. United States. https://doi.org/10.1016/j.cma.2019.112658
Hajduk, Hennes, Kuzmin, Dmitri, Kolev, Tzanio, and Abgrall, Remi. Wed .
"Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations". United States. https://doi.org/10.1016/j.cma.2019.112658. https://www.osti.gov/servlets/purl/1836210.
@article{osti_1836210,
title = {Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations},
author = {Hajduk, Hennes and Kuzmin, Dmitri and Kolev, Tzanio and Abgrall, Remi},
abstractNote = {In this work, we introduce a new residual distribution (RD) framework for the design of bound-preserving high-resolution finite element schemes. The continuous and discontinuous Galerkin discretizations of the linear advection equation are modified to construct local extremum diminishing (LED) approximations. To that end, we perform mass lumping and redistribute the element residuals in a manner which guarantees the LED property. The hierarchical correction procedure for high-order Bernstein finite element discretizations involves localization to subcells and definition of bound-preserving weights for subcell contributions. Using strong stability preserving (SSP) Runge–Kutta methods for time integration, we prove the validity of discrete maximum principles under CFL-like time step restrictions. The low-order version of our method has roughly the same accuracy as the one derived from a piecewise (multi)-linear approximation on a submesh with the same nodal points. In high-order extensions, we use an element-based flux-corrected transport (FCT) algorithm which can be interpreted as a nonlinear RD scheme. The proposed LED corrections are tailor-made for matrix-free implementations which avoid the rapidly growing cost of matrix assembly for high-order Bernstein elements. The results for 1D, 2D, and 3D test problems compare favorably to those obtained with the best matrix-based approaches.},
doi = {10.1016/j.cma.2019.112658},
journal = {Computer Methods in Applied Mechanics and Engineering},
number = ,
volume = 359,
place = {United States},
year = {Wed Oct 09 00:00:00 EDT 2019},
month = {Wed Oct 09 00:00:00 EDT 2019}
}
Web of Science
Works referenced in this record:
A Posteriori Analysis of a Positive Streamwise Invariant Discretization of a Convection-Diffusion Equation
journal, August 2011
- Bernardi, Christine; Chacón Rebollo, Tomás; Restelli, Marco
- Journal of Scientific Computing, Vol. 51, Issue 2
High-order local maximum principle preserving (MPP) discontinuous Galerkin finite element method for the transport equation
journal, April 2017
- Anderson, R.; Dobrev, V.; Kolev, Tz.
- Journal of Computational Physics, Vol. 334
Flux-corrected transport algorithms for continuous Galerkin methods based on high order Bernstein finite elements
journal, September 2017
- Lohmann, Christoph; Kuzmin, Dmitri; Shadid, John N.
- Journal of Computational Physics, Vol. 344
A Second-Order Maximum Principle Preserving Lagrange Finite Element Technique for Nonlinear Scalar Conservation Equations
journal, January 2014
- Guermond, Jean-Luc; Nazarov, Murtazo; Popov, Bojan
- SIAM Journal on Numerical Analysis, Vol. 52, Issue 4
Flux Correction Tools for Finite Elements
journal, January 2002
- Kuzmin, D.; Turek, S.
- Journal of Computational Physics, Vol. 175, Issue 2
A partition of unity approach to adaptivity and limiting in continuous finite element methods
journal, August 2019
- Kuzmin, Dmitri; Quezada de Luna, Manuel; Kees, Christopher E.
- Computers & Mathematics with Applications, Vol. 78, Issue 3
An Example of High Order Residual Distribution Scheme Using non-Lagrange Elements
journal, July 2010
- Abgrall, R.; Trefilík, J.
- Journal of Scientific Computing, Vol. 45, Issue 1-3
On maximum-principle-satisfying high order schemes for scalar conservation laws
journal, May 2010
- Zhang, Xiangxiong; Shu, Chi-Wang
- Journal of Computational Physics, Vol. 229, Issue 9
High-Order Multi-Material ALE Hydrodynamics
journal, January 2018
- Anderson, Robert W.; Dobrev, Veselin A.; Kolev, Tzanio V.
- SIAM Journal on Scientific Computing, Vol. 40, Issue 1
Embedded discontinuous Galerkin transport schemes with localised limiters
journal, April 2016
- Cotter, C. J.; Kuzmin, D.
- Journal of Computational Physics, Vol. 311
Construction of a p-Adaptive Continuous Residual Distribution Scheme
journal, March 2017
- Abgrall, R.; Viville, Q.; Beaugendre, H.
- Journal of Scientific Computing, Vol. 72, Issue 3
Improved detection criteria for the Multi-dimensional Optimal Order Detection (MOOD) on unstructured meshes with very high-order polynomials
journal, July 2012
- Diot, S.; Clain, S.; Loubère, R.
- Computers & Fluids, Vol. 64
High-Resolution Conservative Algorithms for Advection in Incompressible Flow
journal, April 1996
- LeVeque, Randall J.
- SIAM Journal on Numerical Analysis, Vol. 33, Issue 2
Essentially non-oscillatory Residual Distribution schemes for hyperbolic problems
journal, May 2006
- Abgrall, R.
- Journal of Computational Physics, Vol. 214, Issue 2
A Review of Residual Distribution Schemes for Hyperbolic and Parabolic Problems: The July 2010 State of the Art
journal, April 2012
- Abgrall, Remi
- Communications in Computational Physics, Vol. 11, Issue 4
Multidimensional upwinding: Its relation to finite elements
journal, April 1995
- Carette, J. -C.; Deconinck, H.; Paillere, H.
- International Journal for Numerical Methods in Fluids, Vol. 20, Issue 8-9
Differentiable monotonicity-preserving schemes for discontinuous Galerkin methods on arbitrary meshes
journal, June 2017
- Badia, Santiago; Bonilla, Jesús; Hierro, Alba
- Computer Methods in Applied Mechanics and Engineering, Vol. 320
Further variation diminishing properties of Bernstein polynomials on triangles
journal, December 1987
- Goodman, T. N. T.
- Constructive Approximation, Vol. 3, Issue 1
Monotonicity-preserving finite element schemes based on differentiable nonlinear stabilization
journal, January 2017
- Badia, Santiago; Bonilla, Jesús
- Computer Methods in Applied Mechanics and Engineering, Vol. 313
Limiters for high-order discontinuous Galerkin methods
journal, September 2007
- Krivodonova, Lilia
- Journal of Computational Physics, Vol. 226, Issue 1
Strong Stability-Preserving High-Order Time Discretization Methods
journal, January 2001
- Gottlieb, Sigal; Shu, Chi-Wang; Tadmor, Eitan
- SIAM Review, Vol. 43, Issue 1
A posteriori subcell limiting of the discontinuous Galerkin finite element method for hyperbolic conservation laws
journal, December 2014
- Dumbser, Michael; Zanotti, Olindo; Loubère, Raphaël
- Journal of Computational Physics, Vol. 278
Multidimensional upwind schemes based on fluctuation-splitting for systems of conservation laws
journal, January 1993
- Deconinck, H.; Paill�re, H.; Struijs, R.
- Computational Mechanics, Vol. 11, Issue 5-6
A posteriori correction of high-order discontinuous Galerkin scheme through subcell finite volume formulation and flux reconstruction
journal, June 2019
- Vilar, François
- Journal of Computational Physics, Vol. 387
Staggered Grid Residual Distribution Scheme for Lagrangian Hydrodynamics
journal, January 2017
- Abgrall, Rémi; Tokareva, Svetlana
- SIAM Journal on Scientific Computing, Vol. 39, Issue 5
A vertex-based hierarchical slope limiter for <mml:math altimg="si27.gif" display="inline" overflow="scroll" xmlns:xocs="http://www.elsevier.com/xml/xocs/dtd" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.elsevier.com/xml/ja/dtd" xmlns:ja="http://www.elsevier.com/xml/ja/dtd" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:tb="http://www.elsevier.com/xml/common/table/dtd" xmlns:sb="http://www.elsevier.com/xml/common/struct-bib/dtd" xmlns:ce="http://www.elsevier.com/xml/common/dtd" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:cals="http://www.elsevier.com/xml/common/cals/dtd"><mml:mi>p</mml:mi></mml:math>-adaptive discontinuous Galerkin methods
journal, April 2010
- Kuzmin, Dmitri
- Journal of Computational and Applied Mathematics, Vol. 233, Issue 12
Fast inversion of the simplicial Bernstein mass matrix
journal, February 2016
- Kirby, Robert C.
- Numerische Mathematik, Vol. 135, Issue 1
Analysis of Algebraic Flux Correction Schemes
journal, January 2016
- Barrenechea, Gabriel R.; John, Volker; Knobloch, Petr
- SIAM Journal on Numerical Analysis, Vol. 54, Issue 4
Works referencing / citing this record:
A subcell-enriched Galerkin method for advection problems
journal, July 2021
- Rupp, Andreas; Hauck, Moritz; Aizinger, Vadym
- Computers & Mathematics with Applications, Vol. 93
Entropy stabilization and property-preserving limiters for discontinuous Galerkin discretizations of nonlinear hyperbolic equations
preprint, January 2020
- Kuzmin, Dmitri
- arXiv