skip to main content
OSTI.GOV 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

Journal Article · · Journal of Computational Physics

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.

Research Organization:
Univ. of Michigan, Ann Arbor, MI (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA); US Department of the Navy, Office of Naval Research (ONR); National Science Foundation (NSF)
Grant/Contract Number:
FC52-08NA28616; N00014-12-1-0751; CBET 1253157; DEFC52-08NA28616
OSTI ID:
1418534
Alternate ID(s):
OSTI ID: 1247001
Journal Information:
Journal of Computational Physics, Vol. 280, Issue C; ISSN 0021-9991
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 31 works
Citation information provided by
Web of Science

References (66)

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
Superconvergence of discontinuous Galerkin solutions for a nonlinear scalar hyperbolic problem journal May 2006
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
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
Hierarchical reconstruction for discontinuous Galerkin methods on unstructured grids with a WENO-type linear reconstruction and partial neighboring cells journal April 2009
A Simple Method for Compressible Multifluid Flows journal January 1999
Computations of Compressible Multifluids journal May 2001
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
A high-resolution scheme for compressible multicomponent flows with shock waves journal March 2010
A solution-adaptive method for efficient compressible multifluid simulations, with application to the Richtmyer–Meshkov instability journal April 2013
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
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
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
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
A Five-Equation Model for the Simulation of Interfaces between Compressible Fluids journal September 2002
Spherical Solutions of an Underwater Explosion Bubble journal January 1998
Numerical methods for nonconservative hyperbolic systems: a theoretical framework. journal January 2006
Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes journal September 2008
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
Space–time discontinuous Galerkin finite element method for two-fluid flows journal February 2011
On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws journal January 1983
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
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
Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves journal February 2010
A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws journal December 2003
Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities journal September 2009
Numerical simulations of a shock interacting with successive interfaces using the Discontinuous Galerkin method: the multilayered Richtmyer–Meshkov and Rayleigh–Taylor instabilities journal November 2014
Observation and modeling of mixing-layer development in high-energy-density, blast-wave-driven shear flow journal May 2014
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
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
On Godunov-type methods near low densities journal February 1991
A contribution to the great Riemann solver debate journal March 1994
Preventing numerical errors generated by interface-capturing schemes in compressible multi-material flows journal July 2012

Cited By (3)


Figures / Tables (11)