An interface treating technique for compressible multi-medium flow with Runge-Kutta discontinuous Galerkin method

The high-order accurate Runge-Kutta discontinuous Galerkin (RKDG) method is applied to the simulation of compressible multi-medium flow, generalizing the interface treating method given in Chertock et al. (2008) [9]. In mixed cells, where the interface is located, Riemann problems are solved to define the states on both sides of the interface. The input states to the Riemann problem are obtained by extrapolation to the cell boundary from solution polynomials in the neighbors of the mixed cell. The level set equation is solved by using a high-order accurate RKDG method for Hamilton-Jacobi equations, resulting in a unified DG solver for the coupled problem. The method is conservative if we include the states in the mixed cells, which are however not used in the updating of the numerical solution in other cells. The states in the mixed cells are plotted to better evaluate the conservation errors, manifested by overshoots/undershoots when compared with states in neighboring cells. These overshoots/undershoots in mixed cells are problem dependent and change with time. Numerical examples show that the results of our scheme compare well with other methods for one and two-dimensional problems. In particular, the algorithm can capture well complex flow features of the one-dimensional shock entropy wave interaction problem and two-dimensional shock-bubble interaction problem.

[1]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[2]  Boo Cheong Khoo,et al.  A Real Ghost Fluid Method for the Simulation of Multimedium Compressible Flow , 2006, SIAM J. Sci. Comput..

[3]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[4]  Boo Cheong Khoo,et al.  Ghost fluid method for strong shock impacting on material interface , 2003 .

[5]  Alexander Kurganov,et al.  Interface tracking method for compressible multifluids , 2008 .

[6]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[7]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[8]  Boo Cheong Khoo,et al.  The ghost fluid method for compressible gas-water simulation , 2005 .

[9]  Ning Zhao,et al.  Conservative front tracking and level set algorithms , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[10]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[11]  Tiegang Liu,et al.  Runge-Kutta discontinuous Galerkin methods for compressible two-medium flow simulations: One-dimensional case , 2007, J. Comput. Phys..

[12]  Qiang Zhang,et al.  Three-Dimensional Front Tracking , 1998, SIAM J. Sci. Comput..

[13]  R. Fedkiw,et al.  Coupling an Eulerian fluid calculation to a Lagrangian solid calculation with the ghost fluid method , 2002 .

[14]  R. Fedkiw,et al.  A numerical method for two-phase flow consisting of separate compressible and incompressible regions , 2000 .

[15]  Zhiliang Xu,et al.  Conservative Front Tracking with Improved Accuracy , 2003, SIAM J. Numer. Anal..

[16]  S. F. Davis,et al.  An interface tracking method for hyperbolic systems of conservation laws , 1992 .

[17]  Bernardo Cockburn,et al.  The P(1-) RKDG Method for Two-Dimensional Euler Equations of Gas Dynamics , 1991 .

[18]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[19]  Rémi Abgrall,et al.  Computations of compressible multifluids , 2001 .

[20]  R. Abgrall How to Prevent Pressure Oscillations in Multicomponent Flow Calculations , 1996 .

[21]  Chi-Wang Shu TVB uniformly high-order schemes for conservation laws , 1987 .

[22]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

[23]  Timothy G. Leighton,et al.  Shock-induced collapse of a cylindrical air cavity in water: a Free-Lagrange simulation , 2000 .

[24]  James J. Quirk,et al.  On the dynamics of a shock–bubble interaction , 1994, Journal of Fluid Mechanics.

[25]  Chi-Wang Shu,et al.  The Runge-Kutta local projection $P^1$-discontinuous-Galerkin finite element method for scalar conservation laws , 1988, ESAIM: Mathematical Modelling and Numerical Analysis.

[26]  Tiegang Liu,et al.  An adaptive ghost fluid finite volume method for compressible gas-water simulations , 2008, J. Comput. Phys..

[27]  J. Sethian,et al.  A Fast Level Set Method for Propagating Interfaces , 1995 .

[28]  Khoon Seng Yeo,et al.  The simulation of compressible multi-medium flow. I. A new methodology with test applications to 1D gas–gas and gas–water cases , 2001 .

[29]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[30]  Chi-Wang Shu,et al.  A discontinuous Galerkin finite element method for directly solving the Hamilton-Jacobi equations , 2007, Journal of Computational Physics.

[31]  John E. Field,et al.  Shock-induced collapse of single cavities in liquids , 1992, Journal of Fluid Mechanics.

[32]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[33]  R. Fedkiw,et al.  The Ghost Fluid Method for de agration and detonation discontinuities , 1998 .

[34]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[35]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[36]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[37]  R. Fedkiw,et al.  A boundary condition capturing method for incompressible flame discontinuities , 2001 .

[38]  Jianxian Qiu,et al.  Simulations of compressible two-medium flow by Runge-Kutta discontinuous Galerkin methods with the ghost fluid method , 2008 .