Computations of non-reacting and reacting two-fluid interfaces
|
|
- Reynard Todd
- 5 years ago
- Views:
Transcription
1 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 Flows (LEGI), Grenoble INP LaMSID, Sep 6, 0
2 Context & motivation Investigation of the effect of premixed combustion of hydrogen-air mixtures in case of nuclear accident in a reactor containment. Difficulties: Characteristic dimensions: EPR (European Pressurized Reactor): (40 m) 3 ; Reaction zone in a laminar deflagration at atmospheric condition: from about mm (X H 0.96) to 0 mm (X H 0.). Direct Numerical Simulation (DNS) for flame acceleration is impossible in such large geometries. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 / 9
3 Context & motivation Simplifications: Surface reaction (infinite reaction rate, discontinuous interface); Diffusion phenomena (thermal and species) are neglected, and fundamental flame speed is imposed. Challenges: Accurate and robust computation of interfaces on coarse grids. Extension to multi-d unstructured meshes for complex geometries. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 / 9
4 D Chapman-Jouguet deflagration (500 cells) 4 0 State-of-the-art Exact solution Pressure (bar) Position (m) Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 3 / 9
5 Numerical strategies: existing tools Interface modelling: Two-fluid model (two pressures & two velocities): [Delhaye, 968, Ishii, 975, Drew and Passman, 998]; Method for interface propagation: Discrete Equations Method (DEM): [Abgrall and Saurel, 003]; Reactive DEM (RDEM): [Le Métayer et al., 005, Beccantini and Studer, 00]; Accurate discretization of interface: Anti-diffusive method: [Després and Lagoutière, 00, Kokh and Lagoutière, 00]. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 4 / 9
6 Present contribution Combine anti-diffusive strategy with DEM/RDEM. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 5 / 9
7 Contents Preliminaries for flow computation Interface & two-fluid model 3 Discrete Equations Method (DEM) 4 Upwind downwind-controlled splitting (UDCS) 5 Numerical results 6 Conclusions & Perspectives Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 6 / 9
8 Compressible Euler equations Euler equations of gas dynamics: ρ t + ρv x = 0, ρv t + (ρv + p) = 0, x ρe t t + (ρet + pv) x = 0. The system is hyperbolic. The solution contains three characteristic waves associated to the eigenvalues below: λ = v c, λ = v, λ 3 = v + c. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 7 / 9
9 Equations of state (EOS) Non-reacting flow: calorically perfect gas, p = (γ )ρe. stiffened gas (liquids under very high pressures), p = (γ )ρe γp. Reacting flow: thermally perfect gas, p = ρrt; ẽ = h 0 + T 0 c v (τ)dτ. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 8 / 9
10 Finite volume method Integrating the Euler equations We note the integral form as U n+ i u t + f(u) x = 0, = U n i t x (F i+ F i ). U U i F i U i F i+ U i+ i i i i + i + Numerical flux F i± to be defined. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 9 / 9
11 Riemann problem (RP) for D Euler equations red: fluid, blue: fluid U l U l U F(U l ) U r U r x rarefaction v c x i+ t v contact shock v + c U l U r U l = U i U r = U i+ x i+ x Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 0 / 9
12 Interface model j Σ X = X = 0 Σ X = X = 0 P X D I C i The interface is modelled as a discontinuity of characteristic function X k. Interface evolution [Drew and Passman, 998]: dx k = X k + D I X k = 0. dt t Euler equations inside each phase Σ k : u k t + f(u) k = 0. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 / 9
13 Numerical interface model [Abgrall and Saurel, 003] Σ Σ.0 α(volume fraction) Σ Σ Σ Σ Σ x i x i x i+ x Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 / 9
14 Fully conservative Godunov s Method One-fluid method for two-fluid problem. Σ + Σ F i Σ + Σ i i i x t n+ U l U l * t U r * U r F Averaging here i Σ + Σ Σ + Σ t n Ui = U l Ui+ i i = U i r x Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 3 / 9
15 Discrete Equations Method (DEM) [Abgrall and Saurel, 003] α(t = t n ) Σ RP Σ RP Σ α(t = t n + t) Σ F Σ F RP Σ D H C C Σ Σ G B B F x Two-phase Riemann problem 0 i F A i Σ i E x Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 4 / 9
16 Local Riemann problem Non-reacting case (we do not need a model for D I ). t Lgnl CD Rgnl U,L U,R i U,L U,R i Reacting case (D I needs to be modelled). Lgnl CD U,L t Taylor RS U,R U,R { i x U,R Rgnl i U,L i U,R x i Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 5 / 9
17 Upwind & downwind-controlled [Lagoutière, 000] Linear advection (D > 0). α t + D α = 0 αn+ i = α n i x D t x (α i+ α i ). n = 0 upwind everywhere... i i i + i +... n = upwind at i downwind at i +... i i i + i +... i i + exact solution n = downwind controlled upwind... i i i + i +... x x x Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 6 / 9
18 Wave propagation point of view α Σ Σ Σ Σ Σ α i α i,i+ 0 Σ i i i i + i + x Detonation is OK, since precursor shock wave travels together with reacting shock. But in general, shock wave can violate the stability condition. Thus, straightforward coupling of DEM/RDEM with downwind-controlled is unstable. Proposition: Upwind downwind-controlled splitting. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 7 / 9
19 New: Upwind downwind-controlled splitting Do Upwind DEM method, let αu = (α,αρ,αρv,αρe t ) T ; α(t n ) (αu) n k,i Upwind method (αu) n+,up k,i. (a) up n i α n i up n i+ (b) α n i+ 0 i 3 x i i + α(t n + t) α n+,up i α n+,up i+ 0 x i 3 i i + Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 8 / 9
20 New: Upwind downwind-controlled splitting Do non-linear downwind phase splitting ; (b) α(t n + t) α n+,up i α n+,up i+ (c) 0 i 3 x i i + α(t n + t) α n+ i re-arrange the volume α n+ i+ 0 x i 3 i i + Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 9 / 9
21 D reacting flow Before combustion, we have stoichiometric mixture of hydrogen-air; Left chamber: the stoichiometric mixture has been completely burnt (red); H + 0.5O H O. Right chamber: the unburnt (blue) stoichiometric mixture of hydrogen-air; T = 800 K v = 0 m/s p =.03 bar ξ = (burnt) T = 90 K v = 0 m/s p =.03 bar ξ = 0 (unburnt) 0 0 m 0 m Fundamental flame speed K 0 = 00 m/s. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 0 / 9
22 Exact solution of reactive Riemann problem Chapman-Jouguet deflagration [Beccantini and Studer, 00] 4 Exact solution 0 Pressure (bar) t contact Position (m) U r Taylor C-J deflagration shock U l U l U r U r precursor shock x Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 / 9
23 D result of present work (00 cells) Exact solution -nd minmod UDCS anti α (burnt).5.5 ρ (kg/m 3 ) 0. Exact solution -nd minmod UDCS anti x (m) Exact solution -nd minmod UDCS anti p (bar) x (m) x (m) Exact solution -nd minmod UDCS anti v (m/s) x (m) Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 / 9
24 D non-reacting shock bubble interaction 75 mm front side r = 5 mm x = 5 mm y = 44.5 mm R air shock front back side 89 mm 445 mm Location ρ p u x u y γ Air (back) Air (front) R Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 3 / 9
25 Numerical results (DEM+UDCS) Experiment Minmod Minmod Anti-diffusive Anti-diffusive 0 5 quad tri. 0 5 quad tri. (000 00) Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 4 / 9
26 Numerical results (DEM+UDCS) 0.9 -nd 500x00 -nd 000x00 Anti 500x α (air) ρ (kg/m 3 ) nd 500x00 -nd 000x00 Anti 500x x (mm) x (mm) Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 5 / 9
27 D computation of D cylindrical combustion Cartesian grid. First-order (top) v.s. second-order (bottom)..9.8 p(bar).4. ρ(si) Reference x-axis diag r(m) Reference x-axis diag r(m). p(bar).6.4 ρ(si) Reference x-axis diag r(m) Reference x-axis diag Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 6 / 9 r(m)
28 D computation of D cylindrical combustion Cartesian grid. UDCS method. Computed normal (top) v.s. analytical normal (bottom)..4. p(bar).6.4 ρ(si) Reference x-axis diag r(m) Reference x-axis diag r(m).9.8 p(bar).4. ρ(si) Reference x-axis diag r(m) Reference x-axis diag Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 7 / 9 r(m)
29 Conclusions Accurate and robust results obtained using DEM/RDEM+UDCS for reacting and non-reacting flows. Extra benefit: CPU cost reduction. Internal Riemann problem is no longer required; Less diffused fronts leads to a reduced number of (expensive) two-phase Riemann problems to solve..0 α(volume fraction) α(volume fraction) Σ Σ RP Σ Σ Σ Σ Σ RP RP Σ RP Σ x x i x i x i+ x i x i x i+ x UDCS implemented in fast dynamic fluid-structure interaction code EUROPLEXUS (with multi-d unstructured meshes). Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 8 / 9
30 Perspectives Time discretization: first-order forward Euler method for UDCS anti-diffusive. Extension to -level Runge-Kutta stepping. New approach to compute the normal at the flame interface. Extension of UDCS to simplified multiphase models (e.g. one-pressure and one-velocity two-fluid model). Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 9 / 9
31 Bibliography I [Abgrall and Saurel, 003] Abgrall, R. and Saurel, R. (003). Discrete equations for physical and numerical compressible multiphase mixtures. Journal of Computational Physics, 86: [Beccantini and Studer, 00] Beccantini, A. and Studer, E. (00). The reactive Riemann problem for thermally perfect gases at all combustion regimes. International Journal for Numerical Methods in Fluids, 64: [Delhaye, 968] Delhaye, J.-M. (968). Équations fondamentales des écoulements diphasiques. Technical report, Commissariat à l énergie atomique, CEA-R-349. [Després and Lagoutière, 00] Després, B. and Lagoutière, F. (00). Contact Discontinuity Capturing Schemes for Linear Advection and Compressible Gas Dynamics. Journal of Scientific Computing, 6(4): [Drew and Passman, 998] Drew, D. A. and Passman, S. L. (998). Theory of Multicomponent Fluids, Applied Mathematical Sciences, Vol. 35. Springer, New York. [Ishii, 975] Ishii, M. (975). Thermo-fluid Dynamic Theory of Two-phase Flow. Eyrolles, Paris. [Kokh and Lagoutière, 00] Kokh, S. and Lagoutière, F. (00). An anti-diffusive numerical scheme for the simulation of interfaces between compressible fluids by means of a five-equation model. Journal of Computational Physics, 9(8): Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 30 / 9
32 Bibliography II [Lagoutière, 000] Lagoutière, F. (000). Modélisation mathématique et résolution numérique de problèmes de fluides compressibles à plusieurs constituants. PhD thesis, Université Pierre et Marie CURIE. [Le Métayer et al., 005] Le Métayer, O., Massoni, J., and Saurel, R. (005). Modelling evaporation fronts with reactive Riemann solvers. Journal of Computational Physics, 05: Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 3 / 9
33 Averaged equations for two-fluid model Volume fraction: α k def = X k ; α k t + D I α k = 0; Mass α k ρ k t + α k ρ k v k = ρ k,i ( v k,i D I ) α k ; α k ρ k v k t Momentum ) + α k (ρ k v k v k + p k I = Energy The equation is similar. ) ] [ρ k,i v k,i (v k,i D I + p k,i I α k ; Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 3 / 9
34 Closure laws of two-fluid model Rankine-Hugoniot jump conditions: (f,i u,i D I ) α +(f,i u,i D I ) α = 0. (7.) Reacting interface velocity [Beccantini and Studer, 00]: D flame I = v unburnt + K 0 n. (7.) Equations of state (EOS): Non-reacting problem: stiffened gas EOS, p k = (γ k )ρ k (ẽ k h 0,k ) γ k p,k. Reacting problem: thermally perfect gas, p k = ρ k R k T k ; ẽ k = h 0,k + Tk 0 c v,k (τ)dτ. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 33 / 9
35 Multi-dimensional UDCS algorithm Divide the interfaces as inlet and outlet ones: C i,in, n i,in, C i,in,3 C i,in, C i,in, D n i,in, C i C i,out,4 n i,out,4 C i,out,4 D n i,out,4 C i,out,5 Generated volume at j-th inlet/outlet interface: up n i,in,j = ( α n i,in,j α n ) i Ci,in,j D n i,in,j n i,in,j t, up n i,out,j = ( α n i α n i,out,j) Ci,out,j D n i,out,j n i,out,j t. Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 34 / 9
36 Multi-dimensional UDCS algorithm Upwind scheme: α n+ i = α n i + j up n i,in,j. C i (7.3) Downwind scheme (unconditionally unstable): α n+ i = α n i + j up n i,out,j. C i (7.4) Designed high order scheme: α n+ i = α n i + j ( λ n i,in,j ) upn i,in,j C i + j λ n i,out,j 0 λ ; λ = 0 upwind; λ = downwind. up n i,out,j. C i (7.5) Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 35 / 9
37 Two approaches Second order accuracy (use limiter to determine the reconstructed value of α n f,out,j ): λ n i,out,j = αn i αn f,out,j α n i, when α n αn i α n i,out,j. (7.6) i,out,j Downwind-controlled (anti-diffusive) type approach: For simplicity, take a unique value for all outlet interfaces: No Local Extremum of α n+ i : λ n i,out,j = λ n i,out, for any j. m n i αn i + j ( λ n i,in,j ) upn i,in,j C i + λ n up n i,out i,out Mi n C i. (7.7) Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 36 / 9
38 Restrictions 0 λ n i,out ; 0 λ n i,in,j, j; λ n i,out up n i,out if up n i,out < 0; up n i,out [ C i (α n i m n i ) + j [ C i (M n i α n i ) j ( λ n i,in,j) up n i,in,j ( λ n i,in,j) up n i,in,j if up n i,out > 0. (7.8) ], ], Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 37 / 9
39 Sub-region of solution 8 >< >: 0 λ n i,out ; 8 h C up n i (α n i m n i ) + X i sign( up ni,in,j ) up n i,in,j, i,out j >< λ n if up n i,out < 0; i,out h C n i (Mi up n α n i ) X i + sign( up ni,in,j i,out ) up n i,in,j, j >: if up n i,out > 0. (7.9) Tang et al. CEA/Grenoble-INP Two-fluid computation Sep 6, 0 38 / 9
Pour obtenir le grade de. Arrêté ministériel : 7 août 2006
THÈSE Pour obtenir le grade de DOCTEUR DE L UNIVERSITÉ DE GRENOBLE Spécialité : Mécanique des fluides, Énergétique, Procédés Arrêté ministériel : 7 août 6 Présentée par Kunkun TANG tel-8984, version -
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 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 informationModélisation et simulation d écoulements diphasiques avec contrainte granulaire
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 Context Flow
More informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
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 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 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 informationA General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations
A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations Ronald P. Fedkiw Xu-Dong Liu Stanley Osher September, 2 Abstract Standard
More informationTAU Extensions for High Enthalpy Flows. Sebastian Karl AS-RF
TAU Extensions for High Enthalpy Flows Sebastian Karl AS-RF Contents Motivation Extensions available in the current release: Numerical schemes for super- and hypersonic flow fields Models for gas mixtures,
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 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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
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 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 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 informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
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 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 informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationA wave propagation method for compressible multicomponent problems on quadrilateral grids
A wave propagation method for compressible multicomponent problems on quadrilateral grids K.-M. Shyue Department of Mathematics, National Taiwan University, Taipei 16, Taiwan Abstract We describe a simple
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 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 informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
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 informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
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 Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationLecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory
Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory 8.-1 Systems, where fuel and oxidizer enter separately into the combustion chamber. Mixing takes place by convection and diffusion. Only where
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationImplementation of a symmetry-preserving discretization in Gerris
Implementation of a symmetry-preserving discretization in Gerris Daniel Fuster Cols: Pierre Sagaut, Stephane Popinet Université Pierre et Marie Curie, Institut Jean Le Rond D Alembert Introduction 10/11:
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 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 informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
More informationThe multi-stage centred-scheme approach applied to a drift-flux two-phase flow model
The multi-stage centred-scheme approach applied to a drift-flux two-phase flow model Svend Tollak Munkejord 1 Steinar Evje 2 Tore Flåtten 2 1 Norwegian University of Science and Technology (NTNU), Department
More informationLecture 7 Detonation Waves
Lecture 7 etonation Waves p strong detonation weak detonation weak deflagration strong deflagration / 0 v =/ University of Illinois at Urbana- Champaign eflagrations produce heat Thermal di usivity th
More informationAvailable online at Energy Procedia 00 (2010) GHGT-10
Available online at www.sciencedirect.com Energy Procedia 00 (2010) 000 000 Energy Procedia www.elsevier.com/locate/xxx GHGT-10 Depressurization of carbon dioxide in pipelines models and methods Halvor
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 informationAn Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations. Thesis by Cheng Chi
An Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations Thesis by Cheng Chi In Partial Fulfillment of the Requirements For the Degree of Masters of Science King Abdullah
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 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 informationCFD Analysis of Vented Lean Hydrogen Deflagrations in an ISO Container
35 th UKELG Meeting, Spadeadam, 10-12 Oct. 2017 CFD Analysis of Vented Lean Hydrogen Deflagrations in an ISO Container Vendra C. Madhav Rao & Jennifer X. Wen Warwick FIRE, School of Engineering University
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
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 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 informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationA Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with Mie Grüneisen Equation of State
Journal of Computational Physics 171, 678 707 (2001 doi:10.1006/jcph.2001.6801, available online at http://www.idealibrary.com on A Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with
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 informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
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 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 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 informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
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 informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationA WIMF Scheme for the Drift-Flux Two-Phase Flow Model
A WIMF Scheme for the Drift-Flux Two-Phase Flow Model Steinar Eve A,B,D, Tore Flåtten B,A and Svend Tollak Munkeord C A Centre of Mathematics for Applications (CMA), 153 Blindern, NO-316 Oslo, Norway B
More informationDetonation Simulations with a Fifth-Order TENO Scheme
Commun. Comput. Phys. doi: 10.4208/cicp.OA-2018-0008 Vol. x, No. x, pp. 1-37 xxx 201x Detonation Simulations with a Fifth-Order TENO Scheme Haibo Dong 1, Lin Fu 2, Fan Zhang 3,, Yu Liu 4 and Jun Liu 1
More informationNUMERICAL SIMULATION OF HYDROGEN COMBUSTION. Jan-patrice SIMONEAU, FRAMATOME - FRANCE
FR0200503 NUMERICAL SIMULATION OF HYDROGEN COMBUSTION Jan-patrice SIMONEAU, FRAMATOME - FRANCE Novatome - 10, rue Juliette Recamier- F 69456 LYON cedexo6 - France Ph : +33 4 72 74 73 75 - Facs : +33 4
More informationThe Dynamics of Detonation with WENO and Navier-Stokes Shock-Capturing
The Dynamics of Detonation with WENO and Navier-Stokes Shock-Capturing Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM and
More informationccsd , version 1-9 Oct 2005
THE EXPLICIT SIMPLIFIED INTERFACE METHOD FOR COMPRESSIBLE MULTICOMPONENT FLOWS. BRUNO LOMBARD AND ROSA DONAT Abstract. This paper concerns the numerical approximation of the Euler equations for multicomponent
More informationDepartment of Mathematics University of California Santa Barbara, Santa Barbara, California,
The Ghost Fluid Method for Viscous Flows 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 9435 Email:fedkiw@cs.stanford.edu Xu-Dong Liu Department of Mathematics
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationMACH REFLECTION INDUCED DETONATION IN A REACTIVE FLOW. Supervising Professor Name Frank Lu. Donald Wilson. Albert Tong
MACH REFLECTION INDUCED DETONATION IN A REACTIVE FLOW The members of the Committee approve the master s thesis of Walid Cederbond Supervising Professor Name Frank Lu Donald Wilson Albert Tong Copyright
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 informationA Computational Study of High-Speed Droplet Impact
Copyright 2011 Tech Science Press FDMP, vol.7, no.4, pp.329-340, 2011 A Computational Study of High-Speed Droplet Impact T. Sanada 1, K. Ando 2 and T. Colonius 2 Abstract: When a droplet impacts a solid
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationMulticomponent transport impact on turbulent premixed H 2 /O 2 flames
Multicomponent transport impact on turbulent premixed H /O flames Julien de Charentenay and Alexandre Ern ONERA, F-93 Châtillon, France and EC, ECP, F-995 Châtenay-Malabry cedex, France CERMICS, Ecole
More informationDirect numerical simulation of compressible multi-phase ow with a pressure-based method
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2902 Direct numerical simulation of compressible multi-phase ow with a pressure-based
More informationCarbon Science and Technology
ASI RESEARCH ARTICLE Carbon Science and Technology Received:10/03/2016, Accepted:15/04/2016 ------------------------------------------------------------------------------------------------------------------------------
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 informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationA Higher-Order Generalized Ghost Fluid Method for the Poor for the Three-Dimensional Two-Phase Flow Computation of Underwater Implosions
A Higher-Order Generalized Ghost Fluid Method for the Poor for the Three-Dimensional Two-Phase Flow Computation of Underwater Implosions Charbel Farhat, Arthur Rallu, Sriram Shankaran Department of Mechanical
More informationNUMERICAL METHOD FOR THREE DIMENSIONAL STEADY-STATE TWO-PHASE FLOW CALCULATIONS
' ( '- /A NUMERCAL METHOD FOR THREE DMENSONAL STEADY-STATE TWO-PHASE FLOW CALCULATONS P. Raymond,. Toumi (CEA) CE-Saclay DMT/SERMA 91191 Gif sur Yvette, FRANCE Tel: (331) 69 08 26 21 / Fax: 69 08 23 81
More informationMarie Billaud Friess 1 and Samuel Kokh 2
ESAIM: PROCEEDINGS, March 2012, Vol. 35, p. 203-209 Fédération Denis Poisson Orléans-Tours et E. Trélat UPMC, Editors AN ANTI-DIFFUSIVE LAGRANGE-REMAP SCHEME FOR MULTI-MATERIAL COMPRESSIBLE FLOWS WITH
More informationK. Ambika and R. Radha
Indian J. Pure Appl. Math., 473: 501-521, September 2016 c Indian National Science Academy DOI: 10.1007/s13226-016-0200-9 RIEMANN PROBLEM IN NON-IDEAL GAS DYNAMICS K. Ambika and R. Radha School of Mathematics
More informationA G-equation formulation for large-eddy simulation of premixed turbulent combustion
Center for Turbulence Research Annual Research Briefs 2002 3 A G-equation formulation for large-eddy simulation of premixed turbulent combustion By H. Pitsch 1. Motivation and objectives Premixed turbulent
More informationVerified Calculation of Nonlinear Dynamics of Viscous Detonation
Verified Calculation of Nonlinear Dynamics of Viscous Detonation Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM and Joseph
More informationMarie Billaud Friess 1, Benjamin Boutin 2, Filipa Caetano 3, Gloria Faccanoni 4, Samuel Kokh 5, Frédéric Lagoutière 3, 6 and Laurent Navoret 7
ESAIM: PROCEEDINGS, Vol.?, 20, -0 Editors: Will be set by the publisher A SECOND ORDER ANTI-DIFFUSIVE LAGRANGE-REMAP SCHEME FOR TWO-COMPONENT FLOWS Marie Billaud Friess, Benjamin Boutin 2, Filipa Caetano
More informationAnalysis for Steady Propagation of a Generic Ram Accelerator/ Oblique Detonation Wave Engine Configuration
Analysis for Steady Propagation of a Generic Ram Accelerator/ Oblique Detonation Wave Engine Configuration J.M. Powers1, D.R. Fulton2, K.A. Gonthier3, and M.J. Grismer4 Department of Aerospace and Mechanical
More informationCombustion of Kerosene-Air Mixtures in a Closed Vessel
Excerpt from the Proceedings of the COMSOL Conference 010 Paris Combustion of Kerosene-Air Mixtures in a Closed Vessel C. Strozzi *, J-M Pascaud, P. Gillard Institut PRISME, Université d Orléans. *Corresponding
More informationViscous non-linear theory of Richtmyer-Meshkov Instability. Abstract
Viscous non-linear theory of Richtmyer-Meshkov Instability Pierre Carles and Stéphane Popinet Laboratoire de Modélisation en Mécanique, Université Pierre et Marie Curie, Case 162, 4 place Jussieu, 75252
More informationAn interface treating technique for compressible multi-medium flow with. Runge-Kutta discontinuous Galerkin method. Abstract
An interface treating technique for compressible multi-medium flow with Runge-Kutta discontinuous Galerkin method Chunwu Wang 1 and Chi-Wang Shu 2 Abstract The high-order accurate Runge-Kutta discontinuous
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationMATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS
K CHAPMAN & HALL/CRC Monographs and Surveys in Pure and Applied Mathematics I 18 MATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS ANDREI G. KULIKOVSKII NIKOLAI V. POGORELOV ANDREI YU. SEMENOV
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
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 informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationAdvanced Numerical Methods for non-premixed Flames
0 Advanced Numerical Methods for non-premixed Flames Annarita Viggiano Department of Environmental Engineering and Physics University of Basilicata Italy 1. Introduction Engine designers are under increasing
More informationRecently, a novel concept of a multi-mode pulse detonation wave based propulsion system for hypersonic
Detonation Wave Propagation in an Ejector-Augmented Pulse Detonation Rocket Tae-Hyeong Yi, Donald R. Wilson and Frank K. Lu Aerodynamics Research Center, University of Texas at Arlington, Arlington, TX
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
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 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 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 information