A relaxationprojection method for compressible flows. Part I: The numerical equation of state for the Euler equations
Abstract
A new projection method is developed for the Euler equations to determine the thermodynamic state in computational cells. It consists in the resolution of a mechanical relaxation problem between the various subvolumes present in a computational cell. These subvolumes correspond to the ones traveled by the various waves that produce states with different pressures, velocities, densities and temperatures. Contrarily to Godunov type schemes the relaxed state corresponds to mechanical equilibrium only and remains out of thermal equilibrium. The pressure computation with this relaxation process replaces the use of the conventional equation of state (EOS). A simplified relaxation method is also derived and provides a specific EOS (named the Numerical EOS). The use of the Numerical EOS gives a cure to spurious pressure oscillations that appear at contact discontinuities for fluids governed by real gas EOS. It is then extended to the computation of interface problems separating fluids with different EOS (liquidgas interface for example) with the Euler equations. The resulting method is very robust, accurate, oscillation free and conservative. For the sake of simplicity and efficiency the method is developed in a Lagrangeprojection context and is validated over exact solutions. In a companion paper [F. Petitpas, E. Franquet, R. Saurel,more »
 Authors:
 Polytech'Marseille, University Institute of France, Universite de Provence and SMASH Project UMR CNRS 6595  IUSTIINRIA, 5 rue E. Fermi, 13453 Marseille Cedex 13 (France). Email: Richard.Saurel@polytech.univmrs.fr
 Polytech'Marseille, University Institute of France, Universite de Provence and SMASH Project UMR CNRS 6595  IUSTIINRIA, 5 rue E. Fermi, 13453 Marseille Cedex 13 (France)
 Publication Date:
 OSTI Identifier:
 20991578
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 223; Journal Issue: 2; Other Information: DOI: 10.1016/j.jcp.2006.10.004; PII: S00219991(06)004761; Copyright (c) 2006 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; APPROXIMATIONS; COMPRESSIBLE FLOW; EFFICIENCY; EQUATIONS OF STATE; EXACT SOLUTIONS; MIXTURES; MULTIPHASE FLOW; OSCILLATIONS; RELAXATION; THERMAL EQUILIBRIUM
Citation Formats
Saurel, Richard, Franquet, Erwin, Daniel, Eric, and Le Metayer, Olivier. A relaxationprojection method for compressible flows. Part I: The numerical equation of state for the Euler equations. United States: N. p., 2007.
Web. doi:10.1016/j.jcp.2006.10.004.
Saurel, Richard, Franquet, Erwin, Daniel, Eric, & Le Metayer, Olivier. A relaxationprojection method for compressible flows. Part I: The numerical equation of state for the Euler equations. United States. doi:10.1016/j.jcp.2006.10.004.
Saurel, Richard, Franquet, Erwin, Daniel, Eric, and Le Metayer, Olivier. Tue .
"A relaxationprojection method for compressible flows. Part I: The numerical equation of state for the Euler equations". United States.
doi:10.1016/j.jcp.2006.10.004.
@article{osti_20991578,
title = {A relaxationprojection method for compressible flows. Part I: The numerical equation of state for the Euler equations},
author = {Saurel, Richard and Franquet, Erwin and Daniel, Eric and Le Metayer, Olivier},
abstractNote = {A new projection method is developed for the Euler equations to determine the thermodynamic state in computational cells. It consists in the resolution of a mechanical relaxation problem between the various subvolumes present in a computational cell. These subvolumes correspond to the ones traveled by the various waves that produce states with different pressures, velocities, densities and temperatures. Contrarily to Godunov type schemes the relaxed state corresponds to mechanical equilibrium only and remains out of thermal equilibrium. The pressure computation with this relaxation process replaces the use of the conventional equation of state (EOS). A simplified relaxation method is also derived and provides a specific EOS (named the Numerical EOS). The use of the Numerical EOS gives a cure to spurious pressure oscillations that appear at contact discontinuities for fluids governed by real gas EOS. It is then extended to the computation of interface problems separating fluids with different EOS (liquidgas interface for example) with the Euler equations. The resulting method is very robust, accurate, oscillation free and conservative. For the sake of simplicity and efficiency the method is developed in a Lagrangeprojection context and is validated over exact solutions. In a companion paper [F. Petitpas, E. Franquet, R. Saurel, A relaxationprojection method for compressible flows. Part II: computation of interfaces and multiphase mixtures with stiff mechanical relaxation. J. Comput. Phys. (submitted for publication)], the method is extended to the numerical approximation of a nonconservative hyperbolic multiphase flow model for interface computation and shock propagation into mixtures.},
doi = {10.1016/j.jcp.2006.10.004},
journal = {Journal of Computational Physics},
number = 2,
volume = 223,
place = {United States},
year = {Tue May 01 00:00:00 EDT 2007},
month = {Tue May 01 00:00:00 EDT 2007}
}

The relaxationprojection method developed in Saurel et al. [R. Saurel, E. Franquet, E. Daniel, O. Le Metayer, A relaxationprojection method for compressible flows. Part I: The numerical equation of state for the Euler equations, J. Comput. Phys. (2007) 822845] is extended to the nonconservative hyperbolic multiphase flow model of Kapila et al. [A.K. Kapila, Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, Twophase modeling of deflagration to detonation transition in granular materials: reduced equations, Physics of Fluids 13(10) (2001) 30023024]. This model has the ability to treat multitemperatures mixtures evolving with a single pressure and velocity and is particularly interesting formore »

Numerical experiments with the multiresolution scheme for the compressible Euler equations
We investigate the performance of the multiresolution method. The method was recently proposed as a way to speed up computations of compressible flows. The method is based on truncating wavelet coefficients. The method showed good performance on onedimensional test problems. We implement the method in two space dimensions for a more complex compressible flow computation, intended to simulate conditions under which many production CFD codes are running. Our conclusion is that we can in many cases reduce the CPU time, but that the gain in efficiency is not as large as for the onedimensional problems. We furthermore observe that itmore » 
A Quantitative Comparison of Numerical Methods for the Compressible Euler Equations: Fifthorder WENO and PiecewiseLinear Godunov
A numerical study is undertaken comparing a fifthorder version of the weighted essentially nonoscillatory numerical (WENO5) method to a modern piecewiselinear, secondorder, version of Godunov's (PLMDE) method for the compressible Euler Equations. A series of onedimensional test problems are examined beginning with classical linear problems and ending with complex shock interactions. The problems considered are: (1) linear advection of a Gaussian pulse in density, (2) Sod's shock tube problem, (3) the ''peak'' shock tube problem, (4) a version of the Shu and Osher shock entropy wave interaction and (5) the Woodward and Colella interacting shock wave problem. For each problemmore » 
A spectral elementFCT method for the compressible Euler equations
A new algorithm based on spectral element discretizations and fluxcorrected transport concepts is developed for the solution of the Euler equations of inviscid compressible fluid flow. A conservative formulation is proposed based on one and twodimensional cellaveraging and reconstruction procedures, which employ a staggered mesh of GaussChebyshev and GaussLobattoChebyshev collocation points. Particular emphasis is placed on the construction of robust boundary and interfacial conditions in one and twodimensions. It is demonstrated through shocktube problems and twodimensional simulations that the proposed algorithm leads to stable, nonoscillatory solutions of high accuracy. Of particular importance is the fact that dispersion errors are minimal,more » 
A Cartesian grid projection method for the incompressible Euler equations in complex geometries
Many problems in fluid dynamics require the representation of complicated internal or external boundaries of the flow. Here the authors present a method for calculating timedependent incompressible inviscid flow which combines a projection method with a Cartesian grid approach for representing geometry. In this approach, the body is represented as an interface embedded in a regular Cartesian mesh. The advection step is based on a Cartesian grid algorithm for compressible flow, in which the discretization of the body near the flow uses a volumeoffluid representation. A redistribution procedure is used to eliminate timestep restrictions due to small cells where themore »