A hybrid scheme to compute contact discontinuities in one-dimensional Euler systems

Thierry Gallouët; Jean-Marc Hérard; Nicolas Seguin[1]

  • [1] Laboratoire J.-L. Lions, Université Paris VI, France.

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique (2002)

  • Volume: 36, Issue: 6, page 1133-1159
  • ISSN: 0764-583X

Abstract

top
The present paper is devoted to the computation of single phase or two phase flows using the single-fluid approach. Governing equations rely on Euler equations which may be supplemented by conservation laws for mass species. Emphasis is given on numerical modelling with help of Godunov scheme or an approximate form of Godunov scheme called VFRoe-ncv based on velocity and pressure variables. Three distinct classes of closure laws to express the internal energy in terms of pressure, density and additional variables are exhibited. It is shown first that a standard conservative formulation of above mentioned schemes enables to predict “perfectly” unsteady contact discontinuities on coarse meshes, when the equation of state (EOS) belongs to the first class. On the basis of previous work issuing from literature, an almost conservative though modified version of the scheme is proposed to deal with EOS in the second or third class. Numerical evidence shows that the accuracy of approximations of discontinuous solutions of standard Riemann problems is strengthened on coarse meshes, but that convergence towards the right shock solution may be lost in some cases involving complex EOS in the third class. Hence, a blend scheme is eventually proposed to benefit from both properties (“perfect” representation of contact discontinuities on coarse meshes, and correct convergence on finer meshes). Computational results based on an approximate Godunov scheme are provided and discussed.

How to cite

top

Gallouët, Thierry, Hérard, Jean-Marc, and Seguin, Nicolas. "A hybrid scheme to compute contact discontinuities in one-dimensional Euler systems." ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique 36.6 (2002): 1133-1159. <http://eudml.org/doc/246043>.

@article{Gallouët2002,
abstract = {The present paper is devoted to the computation of single phase or two phase flows using the single-fluid approach. Governing equations rely on Euler equations which may be supplemented by conservation laws for mass species. Emphasis is given on numerical modelling with help of Godunov scheme or an approximate form of Godunov scheme called VFRoe-ncv based on velocity and pressure variables. Three distinct classes of closure laws to express the internal energy in terms of pressure, density and additional variables are exhibited. It is shown first that a standard conservative formulation of above mentioned schemes enables to predict “perfectly” unsteady contact discontinuities on coarse meshes, when the equation of state (EOS) belongs to the first class. On the basis of previous work issuing from literature, an almost conservative though modified version of the scheme is proposed to deal with EOS in the second or third class. Numerical evidence shows that the accuracy of approximations of discontinuous solutions of standard Riemann problems is strengthened on coarse meshes, but that convergence towards the right shock solution may be lost in some cases involving complex EOS in the third class. Hence, a blend scheme is eventually proposed to benefit from both properties (“perfect” representation of contact discontinuities on coarse meshes, and correct convergence on finer meshes). Computational results based on an approximate Godunov scheme are provided and discussed.},
affiliation = {Laboratoire J.-L. Lions, Université Paris VI, France.},
author = {Gallouët, Thierry, Hérard, Jean-Marc, Seguin, Nicolas},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique},
keywords = {Godunov scheme; Euler system; contact discontinuities; thermodynamics; conservative schemes},
language = {eng},
number = {6},
pages = {1133-1159},
publisher = {EDP-Sciences},
title = {A hybrid scheme to compute contact discontinuities in one-dimensional Euler systems},
url = {http://eudml.org/doc/246043},
volume = {36},
year = {2002},
}

TY - JOUR
AU - Gallouët, Thierry
AU - Hérard, Jean-Marc
AU - Seguin, Nicolas
TI - A hybrid scheme to compute contact discontinuities in one-dimensional Euler systems
JO - ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique
PY - 2002
PB - EDP-Sciences
VL - 36
IS - 6
SP - 1133
EP - 1159
AB - The present paper is devoted to the computation of single phase or two phase flows using the single-fluid approach. Governing equations rely on Euler equations which may be supplemented by conservation laws for mass species. Emphasis is given on numerical modelling with help of Godunov scheme or an approximate form of Godunov scheme called VFRoe-ncv based on velocity and pressure variables. Three distinct classes of closure laws to express the internal energy in terms of pressure, density and additional variables are exhibited. It is shown first that a standard conservative formulation of above mentioned schemes enables to predict “perfectly” unsteady contact discontinuities on coarse meshes, when the equation of state (EOS) belongs to the first class. On the basis of previous work issuing from literature, an almost conservative though modified version of the scheme is proposed to deal with EOS in the second or third class. Numerical evidence shows that the accuracy of approximations of discontinuous solutions of standard Riemann problems is strengthened on coarse meshes, but that convergence towards the right shock solution may be lost in some cases involving complex EOS in the third class. Hence, a blend scheme is eventually proposed to benefit from both properties (“perfect” representation of contact discontinuities on coarse meshes, and correct convergence on finer meshes). Computational results based on an approximate Godunov scheme are provided and discussed.
LA - eng
KW - Godunov scheme; Euler system; contact discontinuities; thermodynamics; conservative schemes
UR - http://eudml.org/doc/246043
ER -

References

top
  1. [1] R. Abgrall, How to prevent pressure oscillations in multicomponent flow calculations: a quasi conservative approach. J. Comput. Phys. 125 (1995) 150–160. Zbl0847.76060
  2. [2] R. Abgrall and S. Karni, Computations of compressible multifluids. J. Comput. Phys. 169 (2001) 594–623. Zbl1033.76029
  3. [3] G. Allaire, S. Clerc and S. Kokh, A five equation model for the simulation of interfaces between compressible fluids. J. Comput. Phys. 181 (2002) 577–616. Zbl1169.76407
  4. [4] G. Allaire, S. Clerc and S. Kokh, A five equation model for the numerical solution of interfaces in two phase flows. C.R. Acad. Sci. Paris Sér. I Math. 331 (2000) 1017–1022. Zbl1010.76055
  5. [5] T. Barberon, P. Helluy and S. Rouy, Finite Volume simulations of cavitating flows, in Proc. of Third Symposium on Finite Volumes for Complex Applications, R. Herbin and D. Kroner Eds., Hermes Penton Science (2002) 455–462. Zbl1177.76207
  6. [6] M. Barret, E. Faucher and J.M. Hérard, Some schemes to compute unsteady flashing flows. AIAA J. 40 (2002) 905–913. 
  7. [7] S. Bilicki and D. Kardas, Approximation of thermodynamic properties for subcooled water and superheated steam. Polish Academy of Sciences (1991). 
  8. [8] S. Bilicki and J. Kestin, Physical aspects of the relaxation model in two phase flows. Proc. Roy. Soc. London A 428 (1990) 379–397. Zbl0701.76108
  9. [9] S. Bilicki, J. Kestin and M.M. Pratt, A reinterpretation of the results of the moby dick experiments in terms of the non equilibrium model. J. Fluid Eng. 112 (1990) 212–217. 
  10. [10] T. Buffard, T. Gallouët and J.M. Hérard, Schéma VFRoe en variables caractéristiques. Principe de base et applications aux gaz réels. EDF-DER Report HE-41/96/041/A (1996) in French. 
  11. [11] T. Buffard, T. Gallouët and J.M. Hérard, A sequel to a rough Godunov scheme. Application to real gas flows. Comput. & Fluids 29 (2000) 813–847. Zbl0961.76048
  12. [12] M. Buffat and A. Page, Extension of Roe’s solver for multi species real gases. LMFA report, École Centrale de Lyon, Lyon, France (1995). 
  13. [13] S. Clerc, Accurate computation of contact discontinuities in flows with general equations of state. Comput. Methods Appl. Mech. Engrg. 178 (1999) 245–255. Zbl0964.76049
  14. [14] S. Clerc, Numerical simulation of the homogeneous equilibrium model for two phase flows. J. Comput. Phys. 161 (2000) 354–375. Zbl0965.76051
  15. [15] F. Coquel and B. Perthame, Relaxation of energy and approximate Riemann solvers for general pressure laws in fluid dynamics equations. SIAM J. Numer. Anal. 35 (1998) 2223–2249 (in Memory of Ami Harten). Zbl0960.76051
  16. [16] E. Faucher, J.M. Hérard, M. Barret and C. Toulemonde, Computation of flashing flows in variable cross-section ducts. Int. J. Comput. Fluid Dyn. 13 (2000) 365–391. Zbl0976.76047
  17. [17] R.P. Fedkiw, T. Aslam, B. Merriman and S. Osher, A non oscillatory eulerian approach to interfaces in multimaterial flows (the ghost fluid approach). J. Comput. Phys. 152 (1999) 457. Zbl0957.76052MR1699710
  18. [18] T. Gallouët, J.M. Hérard and N. Seguin, Some recent Finite Volume schemes to compute Euler equations using real gas EOS. Internat. J. Numer. Methods Fluids 39-12 (2002) 1073–1138. Zbl1053.76044
  19. [19] T. Gallouët, J.M. Hérard and N. Seguin, An hybrid scheme to compute contact discontinuities in Euler systems. LATP Report 01-027, Université de Provence, France (2001). Zbl1137.65419
  20. [20] T. Gallouët, J.M. Hérard and N. Seguin, On the use of some symmetrizing variables to deal with vacuum (submitted). Zbl1072.76042
  21. [21] S. Gavrilyuk and R. Saurel, Mathematical and numerical modelling of two phase compressible flows with inertia. J. Comput. Phys. 175 (2002) 326–360. Zbl1039.76067
  22. [22] E. Godlewski and P.A. Raviart, Numerical approximation for hyperbolic systems of conservation laws. Springer Verlag (1996). Zbl0860.65075MR1410987
  23. [23] S.K. Godunov, A difference method for numerical calculation of discontinous equations of hydrodynamics. Sbornik (1959) 271–300 (in Russian). Zbl0171.46204
  24. [24] X. Hou and P. G. LeFLoch, Why non conservative schemes converge to wrong solutions. Math. Comp. 62 (1994) 497–530. Zbl0809.65102
  25. [25] A. In, Numerical evaluation of an energy relaxation method for inviscid real fluids. SIAM J. Sci. Comput. 21 (1999) 340–365. Zbl0953.76066
  26. [26] A. In, Méthodes numériques pour les équations de la dynamique des gaz complexes et écoulements diphasiques. Ph.D. thesis, Université Paris VI, France (1999). 
  27. [27] M. Ishii, Thermo-fluid dynamic theory of two-phase flows. Collection de la Direction des Etudes et Recherches d’Electicité de France (1975). Zbl0325.76135
  28. [28] A.K. Kapila, S.F. Son, J.B. Bdzil, R. Menikoff and D.S. Stewart, Two-phase modelling of DDT: structure of the velocity relaxation zone. Phys. Fluids 9 (1997) 3885–3897. 
  29. [29] S. Karni, Multicomponent flow calculations by a consistent primitive algorithm. J. Comput. Phys. 112 (1994) 31–43. Zbl0811.76044
  30. [30] S. Karni, Hybrid multifluid algorithms. SIAM J. Sci. Comput. 17 (1996) 1019–1039. Zbl0860.76056
  31. [31] S. Karni and R. Abgrall, Ghost fluid for the poor: a single fluid algorithm for multifluid. Oberwolfach (2001). MR1882900
  32. [32] R. Kee, J. Miller and T. Jefferson, Chemkin: a general purpose, problem independant transportable fortran chemical kinetics code package. SAND Report 80-8003, Sandia National Laboratories. 
  33. [33] S. Kokh, Aspects numériques et théoriques de la modélisation des écoulements diphasiques compressibles par des méthodes de capture d’interface. Ph.D. thesis, Université Paris VI, France (2001). 
  34. [34] F. Lagoutière, Modélisation mathématique et résolution numérique de problèmes de fluides compressibles à plusieurs constituants. Ph.D. thesis, Université Paris VI, France (2000). 
  35. [35] A. Letellier and A. Forestier, Le problème de Riemann en fluide quelconque. CEA-DMT Report 93/451 (1993) in French. 
  36. [36] R. LeVeque, Numerical methods for conservation laws. Birkhauser (1992). Zbl0847.65053MR1153252
  37. [37] R. Pollak, Die thermodynamischen eigenschaften von wasser dargestellt durch eine kanonische zustands gleichung fur die fluiden homogenen und heterogenen zustande bis 1200 Kelvin und 3000 bars. Ph.D. thesis, Ruhr Universitat, Germany (1974). 
  38. [38] P. Rascle and O. Morvant, Interface utilisateur de Thetis - THErmodynamique en Tables d’InterpolationS. EDF-DER Report HT-13/95021B, Clamart, France (1995) in French. 
  39. [39] P.L. Roe, Approximate Riemann solvers, parameter vectors and difference schemes. J. Comput. Phys. 43 (1981) 357–372. Zbl0474.65066
  40. [40] S. Rouy, Modélisation mathématique et numérique d’écoulements diphasiques compressibles. Ph.D. thesis, Université de Toulon et du Var, France (2000). 
  41. [41] R. Saurel and R. Abgrall, A multiphase Godunov method for compressible multifluid and multiphase flows. J. Comput. Phys. 150 (1999) 425–467. Zbl0937.76053
  42. [42] R. Saurel and R. Abgrall, A simple method for compressible multifluid flows. SIAM J. Sci. Comput. 21 (1999) 1115–1145. Zbl0957.76057
  43. [43] J. Sethian, Level set methods. Cambridge University Press (1996). Zbl0859.76004MR1409367
  44. [44] K.M. Shyue, A fluid mixture type algorithm for compressible multicomponent flow with Van der Waals equation of state. J. Comput. Phys. 156 (1999) 43–88. Zbl0957.76039
  45. [45] J. Smoller, Shock waves and reaction diffusion equations. Springer Verlag (1983). Zbl0508.35002MR688146
  46. [46] E.F. Toro, Riemann solvers and numerical methods for fluid dynamics. Springer Verlag (1997). Zbl0801.76062MR1474503
  47. [47] I. Toumi, Contribution à la modélisation numérique des écoulements diphasiques eau-vapeur. Thèse d’habilitation, Université Paris Sud, France (2000). 

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

                
Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.