DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: On differentiable local bounds preserving stabilization for Euler equations

Abstract

This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.

Authors:
 [1];  [2];  [3];  [4]
  1. Monash Univ., Clayton, VIC (Australia); Centre Internacional de Mètodes Numèrics en Enginyeria (CIMNE), Castelldefels (Spain)
  2. Centre Internacional de Mètodes Numèrics en Enginyeria (CIMNE), Castelldefels (Spain); Univ. Politecnica de Catalunya, Barcelona (Spain)
  3. Clemson Univ., SC (United States)
  4. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Center for Computing Research; Univ. of New Mexico, Albuquerque, NM (United States)
Publication Date:
Research Org.:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); La Caixa Foundation
OSTI Identifier:
1639088
Alternate Identifier(s):
OSTI ID: 1637576
Report Number(s):
SAND-2020-6831J
Journal ID: ISSN 0045-7825; 687134
Grant/Contract Number:  
AC04-94AL85000; NA0003525; LCF/BQ/DE15/10360010
Resource Type:
Accepted Manuscript
Journal Name:
Computer Methods in Applied Mechanics and Engineering
Additional Journal Information:
Journal Volume: 370; Journal ID: ISSN 0045-7825
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
42 ENGINEERING; 97 MATHEMATICS AND COMPUTING; Shock capturing; Euler equations; hyperbolic systems; positivity preservation

Citation Formats

Badia, Santiago, Bonilla, Jesús, Mabuza, Sibusiso, and Shadid, John N. On differentiable local bounds preserving stabilization for Euler equations. United States: N. p., 2020. Web. doi:10.1016/j.cma.2020.113267.
Badia, Santiago, Bonilla, Jesús, Mabuza, Sibusiso, & Shadid, John N. On differentiable local bounds preserving stabilization for Euler equations. United States. https://doi.org/10.1016/j.cma.2020.113267
Badia, Santiago, Bonilla, Jesús, Mabuza, Sibusiso, and Shadid, John N. Thu . "On differentiable local bounds preserving stabilization for Euler equations". United States. https://doi.org/10.1016/j.cma.2020.113267. https://www.osti.gov/servlets/purl/1639088.
@article{osti_1639088,
title = {On differentiable local bounds preserving stabilization for Euler equations},
author = {Badia, Santiago and Bonilla, Jesús and Mabuza, Sibusiso and Shadid, John N.},
abstractNote = {This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.},
doi = {10.1016/j.cma.2020.113267},
journal = {Computer Methods in Applied Mechanics and Engineering},
number = ,
volume = 370,
place = {United States},
year = {Thu Oct 01 00:00:00 EDT 2020},
month = {Thu Oct 01 00:00:00 EDT 2020}
}

Works referenced in this record:

Runge-Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems
journal, September 2001

  • Cockburn, Bernardo; Shu, Chi-Wang
  • Journal of Scientific Computing, Vol. 16, Issue 3, p. 173-261
  • DOI: 10.1023/A:1012873910884

Fusion Simulation Project Workshop Report
journal, June 2008


A finite difference scheme for a system of two conservation laws with artificial viscosity
journal, January 1979


Invariant regions for systems of conservation laws
journal, February 1985


Maps of Convex Sets and Invariant Regions¶for Finite-Difference Systems¶of Conservation Laws
journal, November 2001

  • Frid, Hermano
  • Archive for Rational Mechanics and Analysis, Vol. 160, Issue 3
  • DOI: 10.1007/s002050100166

Synchronized flux limiting for gas dynamics variables
journal, December 2016


Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems
journal, May 2018

  • Mabuza, Sibusiso; Shadid, John N.; Kuzmin, Dmitri
  • Journal of Computational Physics, Vol. 361
  • DOI: 10.1016/j.jcp.2018.01.048

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations
journal, June 2020


Linearity-preserving flux correction and convergence acceleration for constrained Galerkin schemes
journal, March 2012


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
  • DOI: 10.1016/j.cma.2016.09.035

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
  • DOI: 10.1016/j.cma.2017.03.032

Maximum-principle preserving space–time isogeometric analysis
journal, September 2019

  • Bonilla, Jesús; Badia, Santiago
  • Computer Methods in Applied Mechanics and Engineering, Vol. 354
  • DOI: 10.1016/j.cma.2019.05.042

FEMPAR: An Object-Oriented Parallel Finite Element Framework
journal, October 2017

  • Badia, Santiago; Martín, Alberto F.; Principe, Javier
  • Archives of Computational Methods in Engineering, Vol. 25, Issue 2
  • DOI: 10.1007/s11831-017-9244-1

A tutorial-driven introduction to the parallel finite element library FEMPAR v1.0.0
journal, March 2020


The group finite element formulation
journal, April 1983


Analysis of a group finite element formulation
journal, August 2017


Strong Stability-Preserving High-Order Time Discretization Methods
journal, January 2001


A discontinuity-capturing crosswind-dissipation for the finite element solution of the convection-diffusion equation
journal, December 1993


Flux Correction Tools for Finite Elements
journal, January 2002


Approximate Riemann solvers, parameter vectors, and difference schemes
journal, October 1981


Multidimensional FEM-FCT schemes for arbitrary time stepping: MULTIDIMENSIONAL FEM-FCT SCHEMES
journal, May 2003

  • Kuzmin, D.; Möller, M.; Turek, S.
  • International Journal for Numerical Methods in Fluids, Vol. 42, Issue 3
  • DOI: 10.1002/fld.493

A new finite element formulation for computational fluid dynamics: X. The compressible Euler and Navier-Stokes equations
journal, August 1991

  • Shakib, Farzin; Hughes, Thomas J. R.; Johan, Zdeněk
  • Computer Methods in Applied Mechanics and Engineering, Vol. 89, Issue 1-3
  • DOI: 10.1016/0045-7825(91)90041-4

Stabilization and shock-capturing parameters in SUPG formulation of compressible flows
journal, February 2006

  • Tezduyar, Tayfun E.; Senga, Masayoshi
  • Computer Methods in Applied Mechanics and Engineering, Vol. 195, Issue 13-16
  • DOI: 10.1016/j.cma.2005.05.032

Convergence Analysis of Pseudo-Transient Continuation
journal, April 1998


Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems
journal, May 2018

  • Mabuza, Sibusiso; Shadid, John N.; Kuzmin, Dmitri
  • Journal of Computational Physics, Vol. 361
  • DOI: 10.1016/j.jcp.2018.01.048

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations
journal, June 2020