1 E. Romenski , D. Drikakis Fluid Mechanics & Computational Science Group, Cranfield University, UK CONSERVATIVE FORMULATION AND NUMERICAL METHODS FOR MULTIPHASE COMPRESSIBLE MEDIA The financial support from the EU Marie Curie Incoming International Fellowship Programme (contract MIF1-CT-2005-021368) is acknowledged.
24
Embed
1 E. Romenski, D. Drikakis Fluid Mechanics & Computational Science Group, Cranfield University, UK CONSERVATIVE FORMULATION AND NUMERICAL METHODS FOR MULTIPHASE.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
1
E. Romenski , D. Drikakis Fluid Mechanics & Computational Science Group,
Cranfield University, UK
CONSERVATIVE FORMULATION AND NUMERICAL METHODS FOR MULTIPHASE COMPRESSIBLE MEDIA
The financial support from the EU Marie Curie Incoming International Fellowship Programme(contract MIF1-CT-2005-021368) is acknowledged.
2
The problem of multiphase flow modeling lies in the mathematical and numerical formulation of the problem -
there is not yet a widely accepted formulation for the governing equations of multiphase flows
The challenge is associated with the development of a mathematical model that satisfies three important properties:
. hyperbolicity (symmetric hyperbolic system in particular)
. fully conservative form of the governing equations
consistency of the mathematical model with thermodynamic laws.-------------------------------------------------------------------------
3
Two-phase compressible flow models with different velocities, pressures and temperatures
1. Bayer-Nunziato-type model (Baer&Nunziato,1986; Saurel&Abgrall,1999). Governing equations are based on the mass, momentum, and energy balance laws for each phase in which interfacial exchange terms (differential and algebraic) included. Equations are hyperbolic (non-symmetric), but non-conservative
2. We propose extended thermodynamics approach and thermodynamically compatible systems formalism (Godunov-Romenski) to develop multiphase model. Governing equations are written in terms of parameters of state for the mixture and are taking into account a two-phase character of a flow.
Equations are hyperbolic (symmetric) and conservative
4
Class of thermodynamically compatible systems of hyperbolic conservation laws
Thermodynamically compatible system is formulated in terms of generating potential and variables
All equations of the system are written in a conservative form and the system can be transformed to a symmetric hyperbolic form
Many well-posed systems of mathematical physics and continuum mechanics can be written in the form of thermodynamically compatible system.
Examples: gas dynamics, magneto-hydrodynamics, nonlinear and linear elasticity, electrodynamics of moving media, etc.
Flux terms are formed by invariant operatorsgrad, div, curl
6
Thermodynamically compatible system in Eulerian coordinates can be obtained by passing to the new coordinates and corresponding transformation of generating potential and variables
The system can be transformed to the symmetric oneand hyperbolic if L is a convex function
7
Development of the two-phase flow model
1. Introduction a new physical variables in addition to the classical variables (velocity, density, entropy) characterizing two-phase flow.
2. Formulation of new conservation laws for these new variables
in addition to the mass, momentum and energy conservation laws.
3. Introduction a source terms modelling phase interaction and dissipation.
4. Formulation of closing relationships, such as Equation of State for the mixture and functional dependence of source terms on the parameters of state
Instantaneous pressure relaxation is assumed, drag force is neglected
GFORCE is comparable with the 1st order linearized solver
19
WATER/AIR SHOCK TUBETwo-phase flow test case (Saurel&Abgrall, 1999)
Numerical solution with the use of single temperature model (GFORCE flux)
Riemann problem, initial discontinuity at x=0.7 m
Initial data:
Black – exact solutionRed -- 200 cellsBlue --- 800 cells, gives a verygood agreement with the exact solution
Modelling of the moving water/air interface – instantaneous pressure relaxation, infinite drag coefficient
left: right:
water +small amount of
air
air +small amount of
water
20
Shock – bubble interaction (single temperature model) Interaction of shock wave propagating in air with a cylindrical bubble
Shock wave with the Mach number 1.23
Mesh: 300 x 100 cellsRegion: 225 x 44.5 (in millimeter), Bubble radius: 25 (in millimeter)
Light (Helium) and Heavy (Freon R22) bubbles have been considered
[experiments: Haas&Sturtevant (1987)]
2D Test case
21
Shock –bubble interaction (single temperature model)
Light (Helium) bubble
Mixture density
Perfect gas EOS with γ =1.4 for air and γ =1.648 for Helium
Both gases are initially at atmospheric pressure
Shock wave with the Mach number 1.23. The pressure behind the wave is 1.68 atmosphere
Instantaneous pressure relaxation is assumed.
Drag coefficient is
___________________________________________
___________________________________________
22
Shock –bubble interaction (single temperature model)
Heavy (Freon R22) bubble
Mixture density
Perfect gas EOS with γ =1.4 for air and γ =1.249 for R22Both gases are initially at atmospheric pressure
Shock wave with the Mach number 1.23. The pressure behind the wave is 1.68 atmosphere
___________________________________________
___________________________________________
Instantaneous pressure relaxation is assumed.
Drag coefficient is
23
Further Developments
1. Implement a phase to phase transition kinetics2. Include dispersed phase coalescence and breakdown3. Develop high-accuracy numerical methods for 3D flows
Conclusions
A new approach in multiphase flow modelling based on thermodynamicallycompatible systems theory is proposed. A hierarchy of conservative hyperbolic models for two-phase compressible flow is presented and robust numerical method for solving equations of the models is developed.
24
REFERENCES1. Godunov S.K., Romenski E., Elements of continuum mechanics and conservation laws. Kluwer Academic/Plenum Publishers, New York (2003).
2. Godunov S.K., Romenski E., Thermodynamics, conservation laws and symmetric forms of differential equations in mechanics of continuous media, in Computational Fluid Dynamics Review 1995. John Wiley & Sons, New York, (1995), 19--31.
3. Godunov S.K. , Mikhailova T.Yu., Romenski E.I., Systems of thermodynamically coordinated laws of conservation invariant under rotations. Siberian Math. J., V. 37 (1996), 690--705.
4. Romensky E., Thermodynamics and hyperbolic systems of balance laws in continuum mechanics, in Godunov Methods: theory and applications, Kluwer Academic/Plenum publishers, (2001), 745--761.5. Romenski E., Toro E.F., Compressible two-phase flow models: two-pressure models and numerical methods, Computational Fluid Dynamics Journal, V.13 (2004), 403--416.6. Romenski E., Resnyansky A.D., Toro E.F., Conservative hyperbolic formulation for compressible two-phase flow with different phase pressures and temperatures, Quarterly of Applied Mathematics V. 65 (2007), 259-279.
7. Romenski E., Drikakis D., Compressible two-phase flow modelling based on thermodynamically compatible systems of hyperbolic conservation laws, submitted toInt. J. for Numerical Methods in Fluids.