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

Title: A new limiting procedure for discontinuous Galerkin methods applied to compressible multiphase flows with shocks and interfaces

Abstract

Although the Discontinuous Galerkin (dg) method has seen widespread use for compressible flow problems in a single fluid with constant material properties, it has yet to be implemented in a consistent fashion for compressible multiphase flows with shocks and interfaces. Specifically, it is challenging to design a scheme that meets the following requirements: conservation, high-order accuracy in smooth regions and non-oscillatory behavior at discontinuities (in particular, material interfaces). Following the interface-capturing approach of Abgrall [1], we model flows of multiple fluid components or phases using a single equation of state with variable material properties; discontinuities in these properties correspond to interfaces. Here, to represent compressible phenomena in solids, liquids, and gases, we present our analysis for equations of state belonging to the Mie–Grüneisen family. Within the dg framework, we propose a conservative, high-order accurate, and non-oscillatory limiting procedure, verified with simple multifluid and multiphase problems. We show analytically that two key elements are required to prevent spurious pressure oscillations at interfaces and maintain conservation: (i) the transport equation(s) describing the material properties must be solved in a non-conservative weak form, and (ii) the suitable variables must be limited (density, momentum, pressure, and appropriate properties entering the equation of state), coupledmore » with a consistent reconstruction of the energy. Further, we introduce a physics-based discontinuity sensor to apply limiting in a solution-adaptive fashion. We verify this approach with one- and two-dimensional problems with shocks and interfaces, including high pressure and density ratios, for fluids obeying different equations of state to illustrate the robustness and versatility of the method. Lastly, the algorithm is implemented on parallel graphics processing units (gpu) to achieve high speedup.« less

Authors:
ORCiD logo [1];  [1];  [1]
  1. University of Michigan, Ann Arbor, MI (United States)
Publication Date:
Research Org.:
Univ. of Michigan, Ann Arbor, MI (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA); US Department of the Navy, Office of Naval Research (ONR); National Science Foundation (NSF)
OSTI Identifier:
1418534
Alternate Identifier(s):
OSTI ID: 1247001
Grant/Contract Number:  
FC52-08NA28616; N00014-12-1-0751; CBET 1253157; DEFC52-08NA28616
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 280; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS

Citation Formats

Henry de Frahan, Marc T., Varadan, Sreenivas, and Johnsen, Eric. A new limiting procedure for discontinuous Galerkin methods applied to compressible multiphase flows with shocks and interfaces. United States: N. p., 2014. Web. doi:10.1016/j.jcp.2014.09.030.
Henry de Frahan, Marc T., Varadan, Sreenivas, & Johnsen, Eric. A new limiting procedure for discontinuous Galerkin methods applied to compressible multiphase flows with shocks and interfaces. United States. https://doi.org/10.1016/j.jcp.2014.09.030
Henry de Frahan, Marc T., Varadan, Sreenivas, and Johnsen, Eric. Thu . "A new limiting procedure for discontinuous Galerkin methods applied to compressible multiphase flows with shocks and interfaces". United States. https://doi.org/10.1016/j.jcp.2014.09.030. https://www.osti.gov/servlets/purl/1418534.
@article{osti_1418534,
title = {A new limiting procedure for discontinuous Galerkin methods applied to compressible multiphase flows with shocks and interfaces},
author = {Henry de Frahan, Marc T. and Varadan, Sreenivas and Johnsen, Eric},
abstractNote = {Although the Discontinuous Galerkin (dg) method has seen widespread use for compressible flow problems in a single fluid with constant material properties, it has yet to be implemented in a consistent fashion for compressible multiphase flows with shocks and interfaces. Specifically, it is challenging to design a scheme that meets the following requirements: conservation, high-order accuracy in smooth regions and non-oscillatory behavior at discontinuities (in particular, material interfaces). Following the interface-capturing approach of Abgrall [1], we model flows of multiple fluid components or phases using a single equation of state with variable material properties; discontinuities in these properties correspond to interfaces. Here, to represent compressible phenomena in solids, liquids, and gases, we present our analysis for equations of state belonging to the Mie–Grüneisen family. Within the dg framework, we propose a conservative, high-order accurate, and non-oscillatory limiting procedure, verified with simple multifluid and multiphase problems. We show analytically that two key elements are required to prevent spurious pressure oscillations at interfaces and maintain conservation: (i) the transport equation(s) describing the material properties must be solved in a non-conservative weak form, and (ii) the suitable variables must be limited (density, momentum, pressure, and appropriate properties entering the equation of state), coupled with a consistent reconstruction of the energy. Further, we introduce a physics-based discontinuity sensor to apply limiting in a solution-adaptive fashion. We verify this approach with one- and two-dimensional problems with shocks and interfaces, including high pressure and density ratios, for fluids obeying different equations of state to illustrate the robustness and versatility of the method. Lastly, the algorithm is implemented on parallel graphics processing units (gpu) to achieve high speedup.},
doi = {10.1016/j.jcp.2014.09.030},
journal = {Journal of Computational Physics},
number = C,
volume = 280,
place = {United States},
year = {Thu Oct 02 00:00:00 EDT 2014},
month = {Thu Oct 02 00:00:00 EDT 2014}
}

Journal Article:

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

Figures / Tables:

Figure 1 Figure 1: Pressure profile at t = 2 for the advection of a sharp material interface (P = 2, Δ$\mathcal{x}$ = 1/128). Solid red: conservative transport equation and limiting of the conserved variables (fully conservative approach). Dashed green: non-conservative equation and limiting of the conserved variables. Dash-dotted blue: non-conservative equationmore » and modified limiting (our approach).« less

Save / Share:

Works referenced in this record:

How to Prevent Pressure Oscillations in Multicomponent Flow Calculations: A Quasi Conservative Approach
journal, April 1996


A posteriori error estimation for discontinuous Galerkin solutions of hyperbolic problems
journal, January 2002

  • Adjerid, Slimane; Devine, Karen D.; Flaherty, Joseph E.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 191, Issue 11-12
  • DOI: 10.1016/S0045-7825(01)00318-8

Superconvergence of discontinuous Galerkin solutions for a nonlinear scalar hyperbolic problem
journal, May 2006

  • Adjerid, Slimane; Massey, Thomas C.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 195, Issue 25-28
  • DOI: 10.1016/j.cma.2005.06.017

TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: One-dimensional systems
journal, September 1989


The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems
journal, December 1998


The Runge–Kutta Discontinuous Galerkin Method for Conservation Laws V
journal, April 1998

  • Cockburn, Bernardo; Shu, Chi-Wang
  • Journal of Computational Physics, Vol. 141, Issue 2
  • DOI: 10.1006/jcph.1998.5892

Parallel, adaptive finite element methods for conservation laws
journal, April 1994


Limiters for high-order discontinuous Galerkin methods
journal, September 2007


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


Central Discontinuous Galerkin Methods on Overlapping Cells with a Nonoscillatory Hierarchical Reconstruction
journal, January 2007

  • Liu, Yingjie; Shu, Chi-Wang; Tadmor, Eitan
  • SIAM Journal on Numerical Analysis, Vol. 45, Issue 6
  • DOI: 10.1137/060666974

Hierarchical reconstruction for discontinuous Galerkin methods on unstructured grids with a WENO-type linear reconstruction and partial neighboring cells
journal, April 2009

  • Xu, Zhiliang; Liu, Yingjie; Shu, Chi-Wang
  • Journal of Computational Physics, Vol. 228, Issue 6
  • DOI: 10.1016/j.jcp.2008.11.025

A Simple Method for Compressible Multifluid Flows
journal, January 1999


Computations of Compressible Multifluids
journal, May 2001

  • Abgrall, Rémi; Karni, Smadar
  • Journal of Computational Physics, Vol. 169, Issue 2
  • DOI: 10.1006/jcph.2000.6685

Implementation of WENO schemes in compressible multicomponent flow problems
journal, December 2006


WENO Scheme with Subcell Resolution for Computing Nonconservative Euler Equations with Applications to One-Dimensional Compressible Two-Medium Flows
journal, February 2012

  • Xiong, Tao; Shu, Chi-Wang; Zhang, Mengping
  • Journal of Scientific Computing, Vol. 53, Issue 1
  • DOI: 10.1007/s10915-012-9578-7

A high-resolution scheme for compressible multicomponent flows with shock waves
journal, March 2010

  • Kawai, Soshi; Terashima, Hiroshi
  • International Journal for Numerical Methods in Fluids, Vol. 66, Issue 10
  • DOI: 10.1002/fld.2306

An Efficient Shock-Capturing Algorithm for Compressible Multicomponent Problems
journal, May 1998


A Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with van der Waals Equation of State
journal, November 1999


A Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with Mie–Grüneisen Equation of State
journal, August 2001


Numerical simulations of non-spherical bubble collapse
journal, June 2009


A hybrid, center-difference, limiter method for simulations of compressible multicomponent flows with Mie-Grüneisen equation of state
journal, April 2010


A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows
journal, April 1999

  • Saurel, Richard; Abgrall, Rémi
  • Journal of Computational Physics, Vol. 150, Issue 2
  • DOI: 10.1006/jcph.1999.6187

Discrete equations for physical and numerical compressible multiphase mixtures
journal, April 2003


A five equation reduced model for compressible two phase flow problems
journal, January 2005


An interface capturing method for the simulation of multi-phase compressible flows
journal, September 2010

  • Shukla, Ratnesh K.; Pantano, Carlos; Freund, Jonathan B.
  • Journal of Computational Physics, Vol. 229, Issue 19
  • DOI: 10.1016/j.jcp.2010.06.025

HLLC-type Riemann solver for the Baer–Nunziato equations of compressible two-phase flow
journal, May 2010


Runge–Kutta discontinuous Galerkin method for the approximation of Baer and Nunziato type multiphase models
journal, June 2012


A discontinuous Galerkin method for viscous compressible multifluids
journal, March 2010

  • Michoski, C.; Evans, J. A.; Schmitz, P. G.
  • Journal of Computational Physics, Vol. 229, Issue 6
  • DOI: 10.1016/j.jcp.2009.11.033

Runge–Kutta discontinuous Galerkin method for interface flows with a maximum preserving limiter
journal, July 2012


An interface treating technique for compressible multi-medium flow with Runge–Kutta discontinuous Galerkin method
journal, November 2010


A generalized approach for sub- and super-critical flows using the Local Discontinuous Galerkin method
journal, January 2013

  • Gryngarten, Leandro D.; Menon, Suresh
  • Computer Methods in Applied Mechanics and Engineering, Vol. 253
  • DOI: 10.1016/j.cma.2012.08.021

A Five-Equation Model for the Simulation of Interfaces between Compressible Fluids
journal, September 2002

  • Allaire, Grégoire; Clerc, Sébastien; Kokh, Samuel
  • Journal of Computational Physics, Vol. 181, Issue 2
  • DOI: 10.1006/jcph.2002.7143

Spherical Solutions of an Underwater Explosion Bubble
journal, January 1998

  • Wardlaw, Andrew B.; Mair, Hans U.
  • Shock and Vibration, Vol. 5, Issue 2
  • DOI: 10.1155/1998/690105

Numerical methods for nonconservative hyperbolic systems: a theoretical framework.
journal, January 2006

  • Parés, Carlos
  • SIAM Journal on Numerical Analysis, Vol. 44, Issue 1
  • DOI: 10.1137/050628052

Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes
journal, September 2008

  • Castro, Manuel J.; LeFloch, Philippe G.; Muñoz-Ruiz, María Luz
  • Journal of Computational Physics, Vol. 227, Issue 17
  • DOI: 10.1016/j.jcp.2008.05.012

Well-balanced high-order centred schemes for non-conservative hyperbolic systems. Applications to shallow water equations with fixed and mobile bed
journal, June 2009


Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations
journal, January 2008

  • Rhebergen, S.; Bokhove, O.; van der Vegt, J. J. W.
  • Journal of Computational Physics, Vol. 227, Issue 3
  • DOI: 10.1016/j.jcp.2007.10.007

Space–time discontinuous Galerkin finite element method for two-fluid flows
journal, February 2011

  • Sollie, W. E. H.; Bokhove, O.; van der Vegt, J. J. W.
  • Journal of Computational Physics, Vol. 230, Issue 3
  • DOI: 10.1016/j.jcp.2010.10.019

On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws
journal, January 1983

  • Harten, Amiram; Lax, Peter D.; Leer, Bram van
  • SIAM Review, Vol. 25, Issue 1
  • DOI: 10.1137/1025002

Weak solutions of nonlinear hyperbolic equations and their numerical computation
journal, February 1954


The calculation of the interaction of non-stationary shock waves and obstacles
journal, January 1962


Restoration of the contact surface in the HLL-Riemann solver
journal, July 1994

  • Toro, E. F.; Spruce, M.; Speares, W.
  • Shock Waves, Vol. 4, Issue 1
  • DOI: 10.1007/BF01414629

A Simple Extension of the Osher Riemann Solver to Non-conservative Hyperbolic Systems
journal, July 2010


A weak formulation of roe's approximate riemann solver
journal, October 1992


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


Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov's method
journal, July 1979


Weighted Essentially Non-oscillatory Schemes
journal, November 1994

  • Liu, Xu-Dong; Osher, Stanley; Chan, Tony
  • Journal of Computational Physics, Vol. 115, Issue 1
  • DOI: 10.1006/jcph.1994.1187

Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves
journal, February 2010

  • Johnsen, Eric; Larsson, Johan; Bhagatwala, Ankit V.
  • Journal of Computational Physics, Vol. 229, Issue 4
  • DOI: 10.1016/j.jcp.2009.10.028

A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws
journal, December 2003

  • Ren, Yu-Xin; Liu, Miao’er; Zhang, Hanxin
  • Journal of Computational Physics, Vol. 192, Issue 2
  • DOI: 10.1016/j.jcp.2003.07.006

Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities
journal, September 2009

  • Geuzaine, Christophe; Remacle, Jean-François
  • International Journal for Numerical Methods in Engineering, Vol. 79, Issue 11
  • DOI: 10.1002/nme.2579

Observation and modeling of mixing-layer development in high-energy-density, blast-wave-driven shear flow
journal, May 2014

  • Di Stefano, C. A.; Malamud, G.; Henry de Frahan, M. T.
  • Physics of Plasmas, Vol. 21, Issue 5
  • DOI: 10.1063/1.4872223

Efficient implementation of essentially non-oscillatory shock-capturing schemes, II
journal, July 1989


Ghost fluid method for strong shock impacting on material interface
journal, September 2003


Analysis of Numerical Errors Generated by Slowly Moving Shock Waves
journal, May 2013


The Effects of Numerical Viscosities
journal, July 1996


PLIF flow visualization and measurements of the Richtmyer–Meshkov instability of an air/SF 6 interface
journal, August 2002


Effects of WENO flux reconstruction order and spatial resolution on reshocked two-dimensional Richtmyer–Meshkov instability
journal, February 2007

  • Latini, Marco; Schilling, Oleg; Don, Wai Sun
  • Journal of Computational Physics, Vol. 221, Issue 2
  • DOI: 10.1016/j.jcp.2006.06.051

Shock-induced collapse of a cylindrical air cavity in water: a Free-Lagrange simulation
journal, September 2000


A conservative interface method for compressible flows
journal, December 2006


High-fidelity interface tracking in compressible flows: Unlimited anchored adaptive level set
journal, June 2007


Treatment of interface problems with Godunov-type schemes
journal, May 1996

  • Cocchi, J. P.; Saurel, R.; Loraud, J. C.
  • Shock Waves, Vol. 5, Issue 6
  • DOI: 10.1007/BF02434010

On Godunov-type methods near low densities
journal, February 1991


A contribution to the great Riemann solver debate
journal, March 1994

  • Quirk, James J.
  • International Journal for Numerical Methods in Fluids, Vol. 18, Issue 6
  • DOI: 10.1002/fld.1650180603

Preventing numerical errors generated by interface-capturing schemes in compressible multi-material flows
journal, July 2012


Works referencing / citing this record:

Entropy stable, robust and high-order DGSEM for the compressible multicomponent Euler equations
journal, November 2021


Central Upwind Scheme for a Compressible Two-Phase Flow Model
journal, June 2015