Numerical Simulation of Complex and Multiphase Flows Porquerolles, 18-22 April 2005

1. Numerical Simulation of Complex and Multiphase FlowsPorquerolles, 18-22 April 2005 • thanks to: ERCOFTAC, Conseil Général Var, Région PACA, USTV and Stana. • copy of the students passports (ERCOFTAC grants). • extra nights have to be paid today. • lectures (morning) and advanced communications (afternoon).

2. Spring school « numerical simulations of multiphase and complex flows », 18-22 April 2005, Porquerolles. Applications of the finite volumes method for complex flows: from the theory to the practice Philippe HELLUY, ISITV, Université de Toulon, France.

3. I) Introduction to finite volumes for hyperbolic systems of conservation laws. II) An industrial application. III) Mixtures thermodynamics and numerical schemes for phase transition flows.

4. IIntroduction to finite volumes for hyperbolic systems of conservation laws

5. Characteristic curve A characteristic is a curve along which a regular solution of a first order Partial Differential Equation (PDE) is constant Example: let u(x,t) be a solution to u(x(t),t)=Cst  The characteristic has thus the equation That implies that u is an arbitrary function of x-ct.

6. Other example: Burgers equation A regular solution satisfies The characteristic equation is thus or Consider now an initial condition of the form If u were regular, it would imply u(x,t)=1=0 if x>1 and t>1

7. Admissible shock wave A shock wave is a discontinuous solution of a first order PDE, wt+f(w)x=0, emerging from an intersection of characteristics. If the shock curve is parametrized by (x(t),t), the normal vector to the shock is n=(1,-s), with s=x’(t). We denote by [w]=wR -wL the jump of the quantity w in the shock. According to distribution theory, we must have the Rankine-Hugoniot relation (For Burgers: s=1/2(uL+uR) ) The characteristic intersection condition gives Lax’s condition (for Burgers: uL>uR)

8. Entropy condition More practical criterion: often, a supplementary conservation law can be deduced from the conservation law wt+f(w)x=0 If U is convex (Lax entropy), we require that In the sense of distributions. Example: Burgers

9. The Riemann problem Example of the Burgers equation + entropy condition The solution is noted R(x/t,uL,uR) and is Rarefaction wave Shock wave

10. t x/t=l1+ x/t=l2 – x/t=l2+ w1 x/t=l1– w0=wL wm=wR x Hyperbolic systems In order to generalize to a System of Conservation Laws (SCL), wt+f(w)x=0, w in Rm, we must suppose that f’(w) has only real eigenvalues. The solution of Riemann problem has then the form (for m=2) (m+1 constant states separated by shock or rarefaction waves in the (x,t) plane)

11. In one dimension: let t be a time step, h a space step, xi=ih, tn=nt, The cell (or finite volume) Ci is the interval ]xi-1/2,xi+1/2[ Cj Ci Finite volume schemes Godunov flux "Exact" in some sense, and satisfies a discrete entropy inequality In two dimensions the SCL reads wt+f(w)x+g(w)y=0, and one solves a Riemann problem for each edge between two finite volumes in the normal direction to the edge.

12. IIIndustrial application

13. Multifluid model Gas generator: industrial tool to eject device • Two phases (air and water) • Compressible (pressure up to 300 atm) • Duration: 50 ms, thus evaporation and viscosity neglected.

14. The unknowns are the density r(x,t), the velocity u(x,t), the internal energy e(x,t) the pressure p(x,t) and the mass fraction of gas y(x,t) Euler equations (in 1D in order to simplify) Pressure law:

15. Hyperbolicity In order to prove hyperbolicity, we consider the change of variables w(r,u,s,y), where s is the physical entropy satisfying Tds=de+pd(1/r). One gets The eigenvalues are u-c, u, u, u+c with With a good pressure law, the system is thus hyperbolic. We have also that U(w)= –rs is a Lax entropy

16. Practical pressure law We use a pressure law that allows the exact resolution of the Riemann problem and acceptable precision: the stiffened gas equation of state (EOS) then For numerical reasons, it appears that the last conservation law should be replaced by a non-conservative transport equation (Abgrall-Saurel 1996)

17. Numerical results Rouy, 2000

18. IIIMixtures thermodynamics and numerical schemes for phase transition flows.

19. boiling Cavitation

20. Demonstration

21. Bubble collapse near a rigid wall Liquid area heated at the center by a laser pulse Wall 2.0 mm, 70 cells Heated liquid (1500 atm) 1.4 mm 0.15 mm 0.45 mm Ambient liquid (1atm) 2.4 mm, 70 cells

22. Bubble close to a rigid wall Mixture pressure (from 0 to 2ns)

23. Bubble close to a rigid wall Volume Fraction of Vapor (from 0 to 66ns)

24. Mixtures thermodynamics and numerical schemes for phase transition flows • 1) Thermodynamics of a single fluid and of an immiscible mixture of two fluids; • 2) Relaxation scheme for flow with phase transition (coupling the hydrodynamics and the thermodynamics); • 3) Miscible mixtures and super-critical fluids. • T. Barberon and Helluy. Finite volume simulation of cavitating flows. Computers and Fluids, 2004. • P. Helluy, N. Seguin. A simple model for super-critical fluids, 2005, preprint.

25. Single fluid of mass M, volume V and energy E. W=(M,V,E). An extensive variable X is a function of W that is Positively Homogeneous of degree 1 (PH1) : X is an intensive variable if it is PH0: We define the temperature T by The pressure p by The chemical potential This gives Single fluid thermodynamics S(W) is the entropy. It is extensive and concave

26. The gradient of a PH1 function is PH0 thus pressure, temperature and chemical potential are intensive. The euler relation for PH1 functions gives Definition of the specific volume, energy and entropy t, e, s.

27. Hyperbolicity Hyperbolicity if The sound speed c of the fluid computed from p=p(t,e) must be real The sign of p is not important… • H. B. Callen. Thermodynamics and an introduction to thermostatistics, second edition. Wiley and Sons, 1985. • J.-P. Croisille. Contribution à l’étude théorique et à l’approximation par éléments finis du système hyperbolique de la dynamique des gaz multidimensionnelle et multiespèces. PhD thesis, Université Paris VI, 1991.

28. Constraints Mixture of two immiscible fluids (1) and (2) M,V,E M1,V1,E1 Mixture entropy out of equilibrium M2,V2,E2 Equilibrium entropy at equilibrium: (Isobaric law)

29. Out-of-equilibrium specific entropy It is more practical to use intensive variables. out-of-equilibrium pressure and temperature:

30. The maximum is under constraints: at equilibrium, is not necessarily zero… Equilibrium specific entropy

31. Simple example (perfect gases mixture) We suppose temperature and pressure equilibrium The fractions a and z can then be eliminated. Out-of-equilibrium specific entropy (before mass transfer)

32. On the other side, Pressure law out of equilibrium and saturation curve Out of equilibrium, we have a perfect gas law The saturation curve is thus a line in the (T,p) plane.

33. Phase 2 is the most stable Phase 1 is the most stable Phases 1 and 2 are at equilibrium

34. Equilibrium pressure law Let We suppose (for a fixed temperature, fluid (2) is heavier than fluid (1))

35. Comparison between the real and simplified models

36. Partial conclusion: • the temperature and pressure of the mixture are obtained from the entropy; • the equilibrium entropy is solution of a constrained convex optimisation problem. • Next step: • coupling with hydrodynamics; • numerical scheme.

37. One-velocity two-fluid models Instantaneous equilibrium No phase transition Other possible models…

38. Formal limit It is natural to study the weak entropy solutions of the formal limit system: • Generally, this system has several Lax solutions • Liu entropy criterion is more adequate • R. Menikoff and B. J. Plohr. The Riemann problem for fluid flow of real materials. Rev. Modern Phys., 61(1):75–130, 1989. • S. Jaouen. Étude mathématique et numérique de stabilité pour des modèles hydrodynamiques avec transition de phase. PhD thesis, Université Paris VI, November 2001.

39. Standard schemes as Rusanov's may converge towards different solutions CFL=0.9418 CFL=0.9419 density numerical entropy production

40. Relaxation scheme based on entropy optimisation When l=0, the previous system can be written in the classical form 1) Finite volumes scheme (relaxation of the pressure law) 2) Projection on the equilibrium pressure law

41. Other works about relaxation schemes • Yann Brenier. Averaged multivalued solutions for scalar conservation laws. SIAM J. Numer. Anal., 1984. • B. Perthame. Boltzmann type schemes for gas dynamics and the entropy property. SIAM J. Numer. Anal., 1990. • F. Coquel and B. Perthame. Relaxation of energy and approximate Riemann solvers for general pressure laws in fluid dynamics. SIAM J. Numer. Anal., 1998. • Saurel, Richard; Abgrall, Rémi A multiphase Godunov method for compressible multifluid and multiphase flows. J. Comput. Phys., 1999. • G. Chanteperdrix, P. Villedieu, and Vila J.-P. A compressible model for separated two-phase flows computations. In ASME Fluids Engineering Division Summer Meeting. ASME, Montreal, Canada, July 2002. • Stéphane Dellacherie. Relaxation schemes for the multicomponent Euler system. M2AN Math. Model. Numer. Anal., 2003. • …

42. Numerical results with CFL=1, comparison with the Liu solution for a simple Riemann problem Density

43. Velocity

44. Pressure

45. Mixture of stiffened gases • We suppose only temperature equilibrium (elimination of z). • Out of equilibrium, the mixture still satisfies a stiffened gas law: exact Riemann solver in the relaxation step. • Equilibrium is obtained after optimization with respect to a and j. • The pressure law is not analytic.

46. Bubble collapse near a rigid wall Liquid area heated at the center by a laser pulse Wall 2.0 mm, 70 cells Heated liquid (1500 atm) 1.4 mm 0.15 mm 0.45 mm Ambient liquid (1atm) 2.4 mm, 70 cells

47. Bubble close to a rigid wall Mixture pressure (from 0 to 2ns)

48. Bubble close to a rigid wall Volume Fraction of Vapor (from 0 to 66ns)

49. Partial conclusion: • the relaxation scheme is based on entropy optimisation • it seems to converge towards the Liu solution • it can be used in practical configurations Next step: super-critical fluids

50. Mixture of two miscible fluids Constraints: Equilibrium entropy At equilibrium: (Dalton's law)