Modélisation et simulation d écoulements diphasiques avec contrainte granulaire
|
|
- Robert Richard
- 6 years ago
- Views:
Transcription
1 Modélisation et simulation d écoulements diphasiques avec contrainte granulaire T. Gallouët, P. Helluy, J.-M. Hérard, J. Nussbaum LATP Marseille IRMA Strasbourg EDF Chatou ISL Saint-Louis
2 Context Flow of weakly compressible grains (powder, sand, etc.) inside a compressible gas; importance of the granular stress when the grains are packed; averaged model (the solid phase is represented by an equivalent continuous medium); two-velocity, two-pressure mathematical model; relaxation of the pressure equilibrium in order to gain stability; stable and useful approximation for real applications.
3 Granular stress R
4 Averaging process
5 Outlines Two-phase granular flow models Unstable one-pressure model More stable two-pressure model Entropy and closure Thermodynamically coherent granular stress Shape of the granular stress Entropy and granular stress Hyperbolicity study Stability of the relaxed system Instability of the one-pressure model Numerical approximation: splitting method Convection step Relaxation step Numerical results Stability study 1D Combustion chamber 3D Combustion chamber Conclusion Bibliography
6 Two-phase granular flow models
7 Unstable one-pressure model A gas phase (G) k = 1, a solid phase (S) k = 2. Seven unknowns: partial densities ρ k (x, t), velocities u k (x, t), internal energies e k (x, t), gas volume fraction α 1 (x, t). Pressure law: p k = p k (ρ k, e k ). Other definitions: m k = α k ρ k, α 2 = 1 α 1, E k = e k + u2 k 2. The whole six-equation PDE system reads [3] m k,t + (m k u k ) x = 0, (m k u k ) t + (m k u 2 k + α k p k ) x p 1 α k,x = 0, (m k E k ) t + ((m k E k + α k p k )u k ) x + p 1 α k,t = 0, (1) with k = 1, 2 and with the algebraic relation p 2 = p 1 + R( ).
8 More stable two-pressure model The previous system is generally not hyperbolic. Therefore, we relax the algebraic pressure relation. The seven-equation PDE system reads [1] m k,t + (m k u k ) x = 0, (m k u k ) t + (m k uk 2 + α k p k ) x p 1 α k,x = 0, (m k E k ) t + ((m k E k + α k p k )u k ) x + p 1 α k,t = 0, α 1,t + u 2 α 1,x = P. (2) And the source term P has to contain the granular modeling.
9 Entropy and closure The gas entropy s 1 satisfies the standard relation (T k is the temperature of phase k) T 1 ds 1 = de 1 p 1 ρ 2 dρ 1. 1 For the entropy s 2 of the solid phase, we take into account the granular stress. T 2 ds }{{} 2 = de }{{} 2 energy variation internal energy variation p 2 ρ 2 dρ 2 2 }{{} pressure work R m 2 dα 2 }{{} granular work After some computations, we find the following entropy dissipation equation ( k m k s k ) t + ( k m k u k s k ) x = P T 2 (p 1 + R p 2 ).
10 A natural choice that ensures positive dissipation is P = 1 τ p α 1 α 2 (p 1 + R p 2 ), τ p 0. The choice of the pressure relaxation parameter τ p is very important. We can also write p 1 + R p 2 = τ p α 1 α 2 (α 1,t + u 2 α 1,x ). for τ p = t eq = 0, we recover the instantaneous six-equation equilibrium model p 2 = p 1 + R; τ p = t eq p ref where t eq is the characteristic time of the pressure equilibrium and p ref a reference pressure; for τ p > 0, the system should be more stable.
11 Thermodynamically coherent granular stress
12 Shape of the granular stress Usually, R = R(α 2 ) [4]
13 Entropy and granular stress How to choose the granular stress R for compressible solid phases? Exemple: for a stiffened gas equation of state p 2 = p 2 (ρ 2, e 2 ) = (γ 2 1)ρ 2 e 2 γ 2 π 2, γ 2 > 1. We suppose that R = R(ρ 2, α 2 ). We express that the entropy is a closed differential form (1/T 2 is an integrating factor of the energy form) ds 2 = 1 ( de 2 p 2 dρ 2 R ) dα 2. T 2 m 2 We find Particular choice ρ 2 2 R(ρ 2, α 2 ) = ρ γ 2 2 r(α 2). r(α 2 ) = λα γ 2 2 R = λm γ 2 2.
14 Hyperbolicity study
15 Stability of the relaxed system (τ p > 0) Let Y = (α 1, ρ 1, u 1, s 1, ρ 2, u 2, s 2 ) T. In this set of variables the system becomes B(Y ) = Y t + B(Y )Y x = (P, 0, 0, 0, 0, 0, 0) T u 2 ρ 1(u 1 u 2) α 1 u 1 ρ 1 c 2 p 1 ρ 1 u 1,s1 1 ρ 1 u 1 u 2 ρ 2 p 1 p 2 m 2 c 2 2 ρ 2 u 2 p 2,s2 ρ 2 u 2, c k = γ k (p k + π k ) ρ k. det(b(y ) λi ) = (u 2 λ) 2 (u 1 λ)(u 1 c 1 λ)(u 1 +c 1 λ)(u 2 c 2 λ)(u 2 +c 2 λ). The system is hyperbolic [7, 2].
16 Instability of the one-pressure model (τ p = 0) When τ p 0, formally, we end up with a standard single pressure model p 2 = p 1 + R. We can remove an equation (for example the volume fraction evolution) and we find a six-equation system Z = (ρ 1, u 1, s 1, ρ 2, u 2, s 2 ) T. Z t + C(Z)Z x = 0. In the simple case where R = λm γ 2 2, let = α 1 α 2 + δ(α 1 ρ 2 a α 2 ρ 1 c 2 1 ), and δ = α1 1/γ 2 2 λγ 2 ρ γ. 2 2
17 Then we find The eigenvalues can be computed only numerically. We observe that when R 0, the system is generally not hyperbolic. We observe also that when R, we recover hyperbolicity.
18 Numerical approximation: splitting method The two-pressure model is hyperbolic and more stable. The one-pressure model is commonly used in applications on coarse meshes (helped by the high damping of the numerical viscosity...). We can approximate the one-pressure model by the more general two-pressure model. A time step is made of a convection stage and a relaxation stage. In the relaxation stage, we apply the exact or relaxed pressure equilibrium.
19 Convection step Let w = (α 1, m 1, m 1 u 1, m 1 E 1, m 2, m 2 u 2, m 2 E 2 ) T, The system can be written w t + f (w) x + G(w)w x = (P, 0, 0, 0, 0, 0, 0) T In the first half step the source term is omitted. We use a standard Rusanov scheme w i wi n t f n i+1/2 = f (w n i, w n f (a, b) = + f i+1/2 n f i 1/2 n + G(wi n ) w i+1 n w i 1 n = 0, x 2 x i+1) numerical conservative flux, f (a) + f (b) λ (b a). 2 2 For λ large enough, the scheme is entropy dissipative. Typically, we take λ = max [ ρ ( f (a) + G(a) ), ρ ( f (b) + G(b) )].
20 Relaxation step In the second half step, we have formally to solve α 1,t = P, m k,t = u k,t = 0, (m k e k ) t + p 1 α k,t = 0. Because of mass and momentum conservation we have m k = m k and u k = u k. In each cell we have to compute (α 1, p 1, p 2 ) from the previous state w p 2 = p 1 + λm γ τ α 2 α2 p α 1 α 2 t, m 1 e 1 + m 2 e 2 = m1e 1 + m2e 2, (m 1 e 1 m 1e 1) + p 1 (α 1 α 1) = 0.
21 After some manipulations, we have to solve H(α 2 ) = (π 2 π 1 )(α 1 + (γ 1 1)(α 1 α1))(α 2 + (γ 2 1)(α 2 α2)) +(α 2 γ 2 2 m γ 2 2 θ(α τ p α 2 α2 2) + A 2 )(α 1 + (γ 1 1)(α 1 α (1 α 2 ) t 1)) +A 1 (α 2 + (γ 2 1)(α 2 α 2)) = 0, with A k = α k(p k + π k ) > 0. Theorem: for stiffened gas laws and the simple granular law, H is concave, H(0) < 0 and H(β 2 ) > 0 with 0 < β 2 = (1 α 2 ) + α 2 γ 1 γ 1 < 1. There is a unique α 2 in ]0, β 2 [ ]0, 1[ such that H(α 2 ) = 0.
22 Numerical application
23 Stability study We have constructed an entropy dissipative approximation of a non-hyperbolic system! What happens numerically? Consider a simple Riemann problem in the interval [ 1/2, 1/2]. γ 1 = and γ 2 = We compute the solution at time t = The CFL number is 0.9. Data: (L) (R) ρ u p ρ u p α
24 τ p = 0, R = 0 Figure: Void fraction, 50 cells, τ p = 0, R = 0.
25 τ p = 0, R = 0 Figure: Void fraction, 1, 000 cells, τ p = 0, R = 0.
26 τ p = 0, R = 0 Figure: Void fraction, 10, 000 cells, τ p = 0, R = 0.
27 τ p = 0, R = 0 Figure: Void fraction, 100, 000 cells, τ p = 0, R = 0. Linearly unstable but non-linearly stable...
28 τ p > 0, R = 0 Figure: Pressure, 50, 000 cells, τ p = 0, 2, 10, R = 0. The relaxation clearly stabilizes the results...
29 τ p = 0, R > 0 Figure: Pressures p 1 and p 2, 10, 000 cells, τ p = 0, R = 500m γ2 2.
30 τ p = 0, R > 0 6 Figure: Eigenvalues: imaginary part I = Im (λ i ) 2, 1, 000 cells, R = 0 or R = 500m γ2 2. i=1 The granular stress slightly improves the stability...
31 1D Combustion chamber We consider now a simplified gun [3, 5]. The right boundary of the computational domain is moving. We activate the granular stress and other source terms (chemical reaction and drag), which are all entropy dissipative. The instabilities would occur on much finer grids... Other source terms are added. m k,t + (m k u k ) x = ±M, (m k u k ) t + (m k u 2 k + α k p k ) x p 1 α k,x = ±Q, (m k E k ) t + ((m k E k + α k p k )u k ) x + p 1 α k,t = S k, α k,t + v 2 α k,x = ±P.
32 The mass transfer term is defined by the simplified relations M = α 2 ρ 2 3ṙ r ṙ = m/s (combustion velocity of the grains) (3) r = 10 3 m (radius of the grains) The momentum source term is given by Q = Mu 2 D D = Cα 1 α 2 ρ 2 (u 1 u 2 ) u 1 u 2 (drag force) C = 3 4r (simplified shape factor) (4) The energy source terms are S 1 = u 2 D + MQ ex S 2 = u 2 D Q ex = J/kg (chemical combustion energy) (5)
33 Figure: Pressure evolution at the breech and the shot base during time. Comparison between the Gough and the relaxation model.
34 Figure: Porosity at the final time. Relaxation model with granular stress.
35 Figure: Velocities at the final time. Relaxation model with granular stress.
36 Figure: Pressures at the final time. Relaxation model with granular stress.
37 Figure: Density of the solid phase at the final time. Relaxation model with granular stress.
38 3D Combustion chamber We consider now a 3D virtual gun [6] of length 434 mm. The section is varying from 42 to 40 mm. The igniter has a length of 50 mm and a radius of 2.5 mm. The 3x4 holes of the igniter are at mm, 63.5 mm and mm from the breech. On each row they are distributed at 0, 90, 180 and 270.
39 3D Gun Figure: Gun
40 Mesh Figure: Mesh of the combustion chamber
41 Temperature Figure: t=0.048 ms
42 Temperature Figure: t=0.100 ms
43 Temperature Figure: t=0.372 ms
44 Temperature Figure: t=0.668 ms
45 Temperature Figure: t=0.943 ms
46 Temperature Figure: t=1.204 ms
47 Temperature Figure: t=1.454 ms
48 Temperature Figure: t=1.697 ms
49 Temperature Figure: t=1.933 ms
50 Temperature Figure: t=2.000 ms
51 Temperature Figure: t=0.100 ms
52 Granular stress Figure:
53 Granular stress Figure:
54 Granular stress Figure:
55 Granular stress Figure:
56 Granular stress Figure:
57 Conclusion Good generalization of the one pressure models; Extension of the relaxation approach to flows with granular stress; Rigorous entropy dissipation and maximum principle on the volume fraction; Stability for a finite relaxation time; The instability is (fortunately) preserved by the scheme for fast pressure equilibrium; The model can be used in practical configurations (the solid phase remains almost incompressible).
58 Bibliography M. R. Baer and J. W. Nunziato. A two phase mixture theory for the deflagration to detonation transition (ddt) in reactive granular materials. Int. J. for Multiphase Flow, 16(6): , Thierry Gallouët, Jean-Marc Hérard, and Nicolas Seguin. Numerical modeling of two-phase flows using the two-fluid two-pressure approach. Math. Models Methods Appl. Sci., 14(5): , P. S. Gough. Modeling of two-phase flows in guns. AIAA, 66: , A. K. Kapila, R. Menikoff, J. B. Bdzil, S. F. Son, and D. S. Stewart. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations. Physics of Fluids, 13(10): , Julien Nussbaum, Philippe Helluy, Jean-Marc Hérard, and Alain Carrière. Numerical simulations of gas-particle flows with combustion. Flow, Turbulence and Combustion, 76(4): , J. Nussbaum. Modélisation et simulation numérique d un écoulement diphasique de la balistique intérieure. PhD thesis, Université de Strasbourg, December R. Saurel and R. Abgrall. A multiphase Godunov method for compressible multifluid and multiphase flows. Journal of Computational Physics, 150(2): , 1999.
59 Merci de votre attention.
Numerical simulations of gas-particle flows with combustion
Numerical simulations of gas-particle flows with combustion Julien NUSSBAUM French-german Research Institut of Saint-Louis 5, rue du Général Cassagnou BP 34, 68301 SAINT LOUIS CEDEX Philippe HELLUY Laboratoire
More informationA Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems
A Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems Dia Zeidan,a) and Eric Goncalves 2 School of Basic Sciences and Humanities, German Jordanian University, Amman, Jordan
More informationModeling of two-phase flow > Direct steam generation in solar thermal power plants
Modeling of two-phase flow > Direct steam generation in solar thermal power plants Pascal Richter RWTH Aachen University EGRIN Pirirac-sur-Mer June 3, 2015 Novatec Solar GmbH Direct steam generation Solar
More informationOn the Modeling and Simulation of a Laser-Induced Caviation Bubble
1 On the Modeling and Simulation of a Laser-Induced Caviation Bubble Ali Zein Supervisor: Prof. Dr. Gerald Warnecke Co-supervisor: Dr. Maren Hantke Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationNumerical investigation of shock wave dense particles cloud interaction
25 th ICDERS August 2 7, 215 Leeds, UK Numerical investigation of shock wave dense particles cloud interaction Pavel S. Utkin Moscow Institute of Physics and Technology Dolgoprudny, Moscow region, Russia
More informationHervé Guillard, INRIA Projet Smash, B.P. 93, Sophia-Antipolis Cedex, France,
TRAVELING WAVE ANALYSIS OF TWO-PHASE DISSIPATIVE MODELS Hervé Guillard, INRIA Projet Smash, B.P. 93, 06902 Sophia-Antipolis Cedex, France, Herve.Guillard@sophia.inria.fr Joint work with : Mathieu Labois,
More informationA SPLITTING METHOD FOR THE ISENTROPIC BAER-NUNZIATO TWO-PHASE FLOW MODEL
ESAIM: PROCEEDINGS, Vol.?, 01, 1-10 Editors: Will be set by the publisher A SPLITTING METHOD FOR THE ISENTROPIC BAER-NUNZIATO TWO-PHASE FLOW MODEL Frédéric Coquel 1, Jean-Marc Hérard and Khaled Saleh,
More informationTHE numerical simulation of the creation and evolution
Proceedings of the World Congress on Engineering Vol III WCE, July 4-6,, London, U.K. Numerical Simulation of Compressible Two-phase Flows Using an Eulerian Type Reduced Model A. Ballil, Member, IAENG,
More informationComputations of non-reacting and reacting two-fluid interfaces
Computations of non-reacting and reacting two-fluid interfaces Kunkun Tang Alberto Beccantini Christophe Corre Lab. of Thermal Hydraulics & Fluid Mechanics (LATF), CEA Saclay Lab. of Geophysical & Industrial
More informationNumerical Simulation with Finite Volume of Dynamic Liquid-Vapor Phase Transition
Numerical Simulation with Finite Volume of Dynamic Liquid-Vapor Phase Transition G. Faccanoni * S. Kokh * G. Allaire ** * DEN/DANS/DM2S, Commissariat à l Énergie Atomique, Saclay, 91191 Gif-sur-Yvette
More informationA Simple Model for Cavitation with Non-condensable Gases
A Simple Model for Cavitation with Non-condensable Gases Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène Mathis To cite this version: Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène
More informationThe Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow
The Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow D.W. Schwendeman, C.W. Wahle 2, and A.K. Kapila Department of Mathematical Sciences, Rensselaer Polytechnic
More informationA Godunov-type method for the diphasic Baer-Nunziato model
A Godunov-type method for the diphasic Baer-Nunziato model C. Chalons Joint work with : A. Ambroso, P.-A. Raviart with benefit from discussions with F. Coquel (CNRS, Ecole Polytechnique/CMAP) N. Seguin
More informationFinite Volume simulation of cavitating flows
Finite Volume simulation of cavitating flows Thomas Barberon, Philippe Helluy ISITV, Laboratoire MNC, BP 56, 83162 La Valette, France Abstract We propose a numerical method adapted to the modelling of
More informationWELL-BALANCED ROE-TYPE NUMERICAL SCHEME FOR A MODEL OF TWO-PHASE COMPRESSIBLE FLOWS
J. Korean Math. Soc. 51 2014), No. 1, pp. 163 187 http://dx.doi.org/10.4134/jkms.2014.51.1.163 WELL-BALANCED ROE-TYPE NUMERICAL SCHEME FOR A MODEL OF TWO-PHASE COMPRESSIBLE FLOWS Mai Duc Thanh Abstract.
More informationModel adaptation in hierarchies of hyperbolic systems
Model adaptation in hierarchies of hyperbolic systems Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France February 15th, 2012 DFG-CNRS Workshop Nicolas Seguin (LJLL, UPMC) 1 / 29 Outline of the
More informationAnnalisa Ambroso 1, Christophe Chalons 2, Frédéric Coquel 3 and Thomas. Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher RELAXATION AND NUMERICAL APPROXIMATION OF A TWO-FLUID TWO-PRESSURE DIPHASIC MODEL
More informationOn a simple model of isothermal phase transition
On a simple model of isothermal phase transition Nicolas Seguin Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie Paris 6 France Micro-Macro Modelling and Simulation of Liquid-Vapour Flows
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationMulti-water-bag model and method of moments for Vlasov
and method of moments for the Vlasov equation 1, Philippe Helluy 2, Nicolas Besse 3 FVCA 2011, Praha, June 6. 1 INRIA Nancy - Grand Est & University of Strasbourg - IRMA crestetto@math.unistra.fr (PhD
More informationarxiv: v2 [physics.class-ph] 4 Apr 2009
arxiv:0903.4949v2 [physics.class-ph] 4 Apr 2009 Geometric evolution of the Reynolds stress tensor in three-dimensional turbulence Sergey Gavrilyuk and Henri Gouin Abstract The dynamics of the Reynolds
More informationEulerian Two-Phase Flow CFD Simulation Using a Compressible and Equilibrium Eight- Equation Model. Y. Wang 1 and R. D. Reitz
ILASS Americas 27th Annual Conference on Liquid Atomization and Spray Systems, Raleigh, NC, May 2015 Eulerian Two-Phase Flow CFD Simulation Using a Compressible and Equilibrium Eight- Equation Model Y.
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationEcoulements turbulents des eaux peu profondes Turbulent flows in shallow water
Ecoulements turbulents des eaux peu profondes Turbulent flows in shallow water SERGEY GAVRILYUK a, HENRI GOUIN b a. Université d Aix-Marseille & C.N.R.S. U.M.R. 6595, IUSTI, Project SMASH, 5 rue E. Fermi,
More informationA Study of Detonation Evolution and Structure for a Model of Compressible Two-Phase Reactive Flow
A Study of Detonation Evolution and Structure for a Model of Compressible Two-Phase Reactive Flow D. W. Schwendeman, C. W. Wahle and A. K. Kapila Department of Mathematical Sciences, Rensselaer Polytechnic
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationTCCS Published by Elsevier Ltd. Selection and/or peer-review under responsibility of [name organizer]
Energy Procedia 00 (2011) 000 000 Energy Procedia www.elsevier.com/locate/procedia TCCS-6 Depressurization of CO 2 a numerical benchmark study Sigmund Clausen a,c* and Svend Tollak Munkejord b a Gassco
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationPulsation Amplitude Influence on Free Shear Layer of a Vertical Pulsating Jet in Laminar Regime
Pulsation Amplitude Influence on Free Shear Layer of a Vertical Pulsating Jet in Laminar Regime Nawel Khaldi, Salwa Marzouk, Hatem Mhiri, and Philippe Bournot Abstract In this study, we are proposed to
More informationCommissariat à l Energie Atomique - Saclay Department of Nuclear Energy Fluid Modeling and Simulation. MIT, November 3rd, 2008
ON THE NUMERICAL SIMULATION OF GENERAL MULTI-FIELD TWO-PHASE FLOW MODELS Anela KUMBARO Commissariat à l Energie Atomique - Saclay Department of Nuclear Energy Fluid Modeling and Simulation MIT, November
More informationMicro-scale modelling of internally
Micro-scale modelling of internally unstable soils Dr Tom Shire School of Engineering, University of Glasgow 1 st September 2017 Outline Internal instability Micro-scale modelling Hydromechanical criteria
More informationA unified flow theory for viscous fluids
Laboratoire Jacques-Louis Lions, Paris 27/10/2015 A unified flow theory for viscous fluids ILYA PESHKOV CHLOE, University of Pau, France joint work with EVGENIY ROMENSKI Sobolev Institute of Mathematics,
More informationRecent Developments in Overture
Recent Developments in Overture Bill Henshaw Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA, USA. 11th Symposium on Overset Grids and Solution Technology,
More informationMathematical modeling of impact of two metal plates using two-fluid approach
Journal of Physics: Conference Series PAPER OPEN ACCESS Mathematical modeling of impact of two metal plates using two-fluid approach To cite this article: P S Utkin and S V Fortova 018 J. Phys.: Conf.
More informationNumerical resolution of a two-component compressible fluid model with interfaces
Numerical resolution of a two-component compressible fluid model with interfaces Bruno Després and Frédéric Lagoutière February, 25 Abstract We study a totally conservative algorithm for moving interfaces
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationTwo-fluid compressible simulations on GPU cluster
Two-fluid compressible simulations on GPU cluster Philippe Helluy Jonathan Jung To cite this version: Philippe Helluy Jonathan Jung. Two-fluid compressible simulations on GPU cluster. ESAIM: Proceedings
More informationIntroduction of compressible turbulence
Introduction of compressible turbulence 1 Main topics Derive averaged equations for compressible turbulence Introduce a math. technique to perform averaging in presence of density variation Favre average
More informationThe Discontinuous Galerkin (dg) method is a popular numerical method that has received much attention
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2595 Discontinuous Galerkin method for multifluid Euler equations Marc
More informationEulerian interface-sharpening methods for hyperbolic problems
Eulerian interface-sharpening methods for hyperbolic problems Application to compressible multiphase flow Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan 11:55-12:25, March 05,
More informationA well-balanced approximate Riemann solver for compressible flows in variable cross-section ducts
A well-balanced approximate Riemann solver for compressible flows in variable cross-section ducts Philippe Helluy, Jean-Marc Hérard, Hélène Mathis To cite this version: Philippe Helluy, Jean-Marc Hérard,
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS 1D-time averaged models
A SHORT INTRODUCTION TO TWO-PHASE FLOWS 1D-time averaged models Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +3304 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationConstraints on Equation of State for Cavitating Flows with Thermodynamic Effects
Constraints on Equation of State for Cavitating Flows with Thermodynamic Effects Eric Goncalvès da Silva, Régiane Fortes Patella To cite this version: Eric Goncalvès da Silva, Régiane Fortes Patella. Constraints
More informationComparison of Two Computational Methods for High-Velocity Cavitating Flows around Conical Projectiles : OTi-HULL Hydrocode and Two-Phase Flow Method
CAV200:sessionB3.003 - -Schmitt Page 29/03/0 Comparison of Two Computational Methods for High-Velocity Cavitating Flows around Conical Projectiles : OTi-HULL Hydrocode and Two-Phase Flow Method M. Schaffar
More informationClosure conditions for non-equilibrium multi-component models
Closure conditions for non-equilibrium multi-component models Siegfried Müller, Maren Hante 2 and Pascal Richter 3 Preprint No. 44 December 204 Key words: multi-component flows, entropy, relaxation, phase
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationAn Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term
An Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term Alexander Kurganov Department of Mathematics, Tulane University, 683 Saint Charles Avenue, New Orleans, LA 78,
More informationSCHEDULE. Gonca Aki (Berlin) An incompressible diffuse flow with phase transition
SCHEDULE TUESDAY AFTERNOON 2:00pm Sergey Gavrilyuk (Marseille) Diffuse interface model for compressible fluid - compressible elastic-plastic solid interaction 2:45pm Gonca Aki (Berlin) An incompressible
More informationSteady Deflagration Structure in Two-Phase Granular Propellants
12th ICDERS, Ann Arbor, Michigan July 23-28, 1989 Steady Deflagration Structure in Two-Phase Granular Propellants Joseph M. Powers 1, Mark E. Miller2, D. Scott Stewart3, and Herman Krier4 July 23-28, 1989
More informationPalindromic Discontinuous Galerkin Method
Palindromic Discontinuous Galerkin Method David Coulette, Emmanuel Franck, Philippe Helluy, Michel Mehrenberger, Laurent Navoret To cite this version: David Coulette, Emmanuel Franck, Philippe Helluy,
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More informationTWO-DIMENSIONAL VISCOUS COMPACTION IN COMPRESSIBLE GRANULAR MATERIALS. A Thesis. Submitted to the Graduate School. of the University of Notre Dame
TWO-DIMENSIONAL VISCOUS COMPACTION IN COMPRESSIBLE GRANULAR MATERIALS A Thesis Submitted to the Graduate School of the University of Notre Dame in Partial Fulfillment of the Requirements for the Degree
More informationMinimum fluidization velocity, bubble behaviour and pressure drop in fluidized beds with a range of particle sizes
Computational Methods in Multiphase Flow V 227 Minimum fluidization velocity, bubble behaviour and pressure drop in fluidized beds with a range of particle sizes B. M. Halvorsen 1,2 & B. Arvoh 1 1 Institute
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations
A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationarxiv: v2 [physics.comp-ph] 17 Oct 2014
A Technique for Computing Dense Granular Compressible Flows with Shock Waves Ryan W. Houim, Elaine S. Oran Department of Aerospace Engineering, University of Maryland, College Park, MD 074, USA arxiv:131.190v
More informationTowards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below
Copyright 2008 Tech Science Press FDMP, vol.4, no.4, pp.263-269, 2008 Towards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below G. Accary 1, S. Meradji
More informationRecent advances in numerical methods for compressible two-phase flow with heat & mass transfers
Recent advances in numerical methods for compressible two-phase flow with heat & mass transfers Keh-Ming Shyue Institute of Applied Mathematical Sciences National Taiwan University Joint work with Marica
More informationA second order scheme on staggered grids for detonation in gases
A second order scheme on staggered grids for detonation in gases Chady Zaza IRSN PSN-RES/SA2I/LIE August 1, 2016 HYP 2016, RWTH Aachen Joint work with Raphaèle Herbin, Jean-Claude Latché and Nicolas Therme
More informationUn schéma volumes finis well-balanced pour un modèle hyperbolique de chimiotactisme
Un schéma volumes finis well-balanced pour un modèle hyperbolique de chimiotactisme Christophe Berthon, Anaïs Crestetto et Françoise Foucher LMJL, Université de Nantes ANR GEONUM Séminaire de l équipe
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationInterior ballistic two-phase flow model of guided-projectile gun system utilizing stick propellant charge
Mahmoud Rashad, Xiaobing Interior ballistic two-phase flow model of guided-projectile gun system utilizing stick propellant charge MAHMOUD RASHAD School of Energy and Power Engineering Nanjing University
More informationarxiv: v1 [math.ap] 12 Oct 2017
A THERMODYNAMICALLY CONSISTENT MODEL OF A LIQUID-VAPOR FLUID WITH A GAS arxiv:1710.04525v1 [math.ap] 12 Oct 2017 Abstract. This work is devoted to the consistent modeling of a three-phase mixture of a
More informationTurbulence models and excitation of solar oscillation modes
Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives
More informationESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria
ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr DeCaria References: An Introduction to Dynamic Meteorology, Holton MOMENTUM EQUATIONS The momentum equations governing the ocean or atmosphere
More informationCFD Investigations of Effects of Cohesive Particles Proportion on Fluidization of Binary Particles
Proceedings of the 2 nd World Congress on Momentum, Heat and Mass Transfer (MHMT 17) Barcelona, Spain April 6 8, 2017 Paper No. ICMFHT 122 ISSN: 2371-5316 DOI: 10.11159/icmfht17.122 CFD Investigations
More informationA Reaction Progress Variable Modeling Approach for Non-Ideal Explosives
50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 09-12 January 2012, Nashville, Tennessee AIAA 2012-0980 A Reaction Progress Variable Modeling Approach for
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol II - Computational Methods for Compressible Flow Problems - Remi Abgrall
COMPUTATIONAL METHODS FOR COMPRESSIBLE FLOW PROBLEMS Rémi Abgrall Université Bordeaux I, Talence Cedex, France Keywords: computational methods, numerical schemes, 1- D problems, multidimensional problems,
More informationCalculating Storm Surge and Other Coastal Hazards Using Geoclaw
Calculating Storm Surge and Other Coastal Hazards Using Geoclaw Kyle T. Mandli Department of Applied Mathematics University of Washington Seattle, WA, USA Modeling and Computations of Shallow-Water Coastal
More informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
More informationThe Continuum Navier-Stokes Equations
I (A) The Continuum Navier-Stokes Equations Foundations of Macroscopic Hydrodynamics See also Tennekes & Lumley, Section 2.2 Hamiltonian equations of a simple molecular fluid: Conservation laws: H = p
More informationSimplified Hyperbolic Moment Equations
Simplified Hyperbolic Moment Equations Julian Koellermeier and Manuel Torrilhon Abstract Hyperbolicity is a necessary property of model equations for the solution of the BGK equation to achieve stable
More informationTHEORY OF DETONATION STRUCTURE FOR TWO-PHASE MATERIALS JOSEPH MICHAEL POWERS. B.S., University of Illinois, 1983 M. S., University of Illinois, 1985
THEORY OF DETONATION STRUCTURE FOR TWO-PHASE MATERIALS BY JOSEPH MICHAEL POWERS B.S., University of Illinois, 1983 M. S., University of Illinois, 1985 THESIS Submitted in partial fulfillment of the requirements
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationNumerical Simulation of the Hagemann Entrainment Experiments
CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationNotes. Multi-Dimensional Plasticity. Yielding. Multi-Dimensional Yield criteria. (so rest state includes plastic strain): #=#(!
Notes Multi-Dimensional Plasticity! I am back, but still catching up! Assignment is due today (or next time I!m in the dept following today)! Final project proposals: I haven!t sorted through my email,
More informationNeeds work : define boundary conditions and fluxes before, change slides Useful definitions and conservation equations
Needs work : define boundary conditions and fluxes before, change slides 1-2-3 Useful definitions and conservation equations Turbulent Kinetic energy The fluxes are crucial to define our boundary conditions,
More informationTarget Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu
Muon Collider/Neutrino Factory Collaboration Meeting May 26 28, CERN, Geneva U.S. Department of Energy Target Simulations Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu Center for Data Intensive
More informationShock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids
Shock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids Andrea Mentrelli Department of Mathematics & Research Center of Applied Mathematics (CIRAM) University of Bologna, Italy Summary
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 21, No. 3, pp. 1115 1145 c 1999 Society for Industrial and Applied Mathematics A SIMPLE METHOD FOR COMPRESSIBLE MULTIFLUID FLOWS RICHARD SAUREL AND RÉMI ABGRALL Abstract. A simple
More informationNumerical and mathematical analysis of a five-equation model for two-phase flow
SINTEF Energy Research Master s Thesis PFE-Master Numerical and mathematical analysis of a five-equation model for two-phase flow Author Pedro José Martínez Ferrer Supervisors Svend Tollak Munkejord Tore
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationProspects for High-Speed Flow Simulations
Prospects for High-Speed Flow Simulations Graham V. Candler Aerospace Engineering & Mechanics University of Minnesota Support from AFOSR and ASDR&E Future Directions in CFD Research: A Modeling & Simulation
More informationModified Porosity Distributed Resistance Combined to Flamelet Combustion Model for Numerical Explosion Modelling
Modified Porosity Distributed Resistance Combined to Flamelet Combustion Model for Numerical Explosion Modelling Dr. Sávio Vianna School of Chemical Engineering University of Campinas - UNICAMP University
More informationy =ρw, w=v+p(ρ). (1.2)
TRANSPORT-EQUILIBRIUM SCHEMES FOR COMPUTING CONTACT DISCONTINUITIES IN TRAFFIC FLOW MODELING CHRISTOPHE CHALONS AND PAOLA GOATIN Abstract. We present a very efficient numerical strategy for computing contact
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationJournal of Computational Physics
Journal of Computational Physics 5 (03) 90 309 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics www.elsevier.com/locate/jcp A diffuse interface model with immiscibility
More informationChristoph Altmann 1, Thomas Belat 2, Michael Gutnic, Philippe Helluy, Hélène Mathis, Éric Sonnendrücker 3, Wilfredo Angulo 4 and Jean-Marc Hérard 5
ESAIM: PROCEEDINGS, August 2009, Vol. 28, p. 33-54 M. Ismail, B. Maury & J.-F. Gerbeau, Editors A LOCAL TIME-STEPPING DISCONTINUOUS GALERKIN ALGORITHM FOR THE MHD SYSTEM Christoph Altmann 1, Thomas Belat
More informationTransport equation cavitation models in an unstructured flow solver. Kilian Claramunt, Charles Hirsch
Transport equation cavitation models in an unstructured flow solver Kilian Claramunt, Charles Hirsch SHF Conference on hydraulic machines and cavitation / air in water pipes June 5-6, 2013, Grenoble, France
More informationA variational approach to design a numerical scheme on an arbitrary moving grid for N-fluid flow with thermodynamic consistency
A variational approach to design a numerical scheme on an arbitrary moving grid for N-fluid flow with thermodynamic consistency T. VAZQUEZ-GONZALEZ a, A. LLOR b, C. FOCHESATO c a. CEA, DAM, DIF, 91297
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More information