A second order scheme on staggered grids for detonation in gases
|
|
- Chastity Clark
- 5 years ago
- Views:
Transcription
1 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
2 Introduction 1/23 Fukushima I Nuclear Power Plant, Unit 3, 14 March Hydrogen-air explosion severely damaging the reactor building and leaking radioactive materials in the atmosphere.
3 2/23 Introduction Previous studies at CEA/IRSN in the 80 s concluded that J.-M. Evrard. General description of hydrogen related risks in pressurized water reactor containments. Tech. Rep. SAER/84/392, CEA (1984): If direct initiation of a detonation appears highly improbable, the flame acceleration mechanism, on the hand, is liable to lead to its formation. (...) The turbulences induced by the obstacles in front of the flame cause (...) acceleration of the flame front and may lead to a deflagration detonation transition. A. Forestier and J. Brochard. Tenue d une enceinte de REP 1300 en cas de détonation hydrogène. Tech. Rep. DEMT/87/276, CEA (1987): The chosen scenario leads very likely to a very serious damage of the containment building. Nowadays, the ruin (of the dome in particular) cannot be excluded. Simulating hydrogen-oxygen explosions in ISIS/P2REMICS code: Mixing Low-Mach Low speed Ignition Deflagration Subsonic High speed DDT Detonation Supersonic High speed Projection scheme Explicit scheme Common Staggered Finite-Volume Mesh
4 3/23 Introduction The Zeldovich-von Neumann-Döring (ZND) model for detonation: Steady 1D combustion wave: inert shock wave + chemical reaction (1) with finite rate Self-sustained: shock wave brings reactants to ignition while expansion of products (2) drives shock forward p State B p M S 1 State C State A products (2) (1) reactants ZND states on pressure profile x B C A y = 1 M y = 0 S τ ZND states on Hugoniot curves
5 Introduction The ZND model is a good approximation but detonations are a genuinely multi-dimensional phenomena, driven by transverse instabilities: R T 1 S M R S T 2 I R T 3 S M R S T 4 detonation front at t 1 R S T 1 R S R S M R T 4 detonation S front at t 2 I T 2 T 3 I Front structure: (I)ncident shock (R)eflected shock (M)ach stem (S)lipline (T)riple point Left: soot tracks of a planar detonation in a channel (Strehlow, 1968). Right: open-shutter photograph of a converging cylindrical detonation in a narrow channel (Knystautas Lee, 1971). 4/23
6 Physical Model 5/23 Unimolecular reaction of a perfect gas with reaction rate defined by either Heaviside law or Arrhenius law: ( ) R P EA κ(t ) = κ γ R = γ P = γ 0 exp RT Equation of state: e = p/ρ γ 1 + q 0y Reactive Euler equations: t ρ + div(ρu) = 0 t (ρu) + div(ρu u) + p = 0 t (ρe) + div(ρeu) + div(pu) = 0 t (ρy) + div(ρyu) = ρyκ(t ) Controlling positivity of e : solve an internal energy balance instead?
7 Space Discretization 6/23 Staggered Finite Volume meshes: K P K D σ ε=d σ D σ D σ L u y D σ u y D σ M L u x D σ M Rannacher-Turek/Crouzeix-Raviart grids MAC grid with unknowns: Primal grid ρ K, p K, e K, y K K M Dual grid u Dσ,i σ E S,i, 1 i d
8 Numerical Scheme: Time Algorithm (1) Hydrodynamic step (continuity, mass fraction, internal energy) ρ n+1 K = ρ n K δt div K (ρ n u n ) ȳ n+1 K = 1 ρ n+1 K e n+1 K = 1 ρ n+1 K (2) Chemistry step [ρ n K y n K δt div K (ρ n u n y n )] [ ρ n K en K δt(div K (ρ n e n u n ) + pk n div K (u n ) + 1 ] K Sn K ) y n+1 K = ȳ n+1 K exp( κ(t n K )δt) with T n K = Mpn K Rρ n K (3) Equation of state p n+1 K = (γ 1)ρ n+1 K (en+1 K q 0 y n+1 K ) (4) Hydrodynamic step (momentum) u n+1 σ,i = 1 [ρ n ρ n+1 Dσ uσ,i n δt ( div Dσ (ρ n u n ui n ) + ( ) ] pn+1 ) σ,i D σ 7/23
9 Numerical Scheme: Convection Operators Convection operators on the primal grid discretized as: div K (ρuz) = 1 K σ E(K) F K,σ z σ with F K,σ = σ ρ σ ( d i=1 Convection operators on the dual grid discretized as: div Dσ (ρuz) = 1 F σ,ε z ε D σ ε E(D σ) u σ,i n K,σ,i ) with F σ,ε calculated such that (Ansanay-Alex Babik Latché Vola, 2010): D σ ( ρ n+1 D δt σ ρ n ) D σ + Fσ,ε n = 0 ε E(D σ) Convected quantities: upwinding with respect to the material velocity u n. Improved accuracy at contact waves : MUSCL-like interpolation of ρ σ, e σ with a limitation procedure such that p and u are well preserved at contacts (Therme, 2015) p σ = βp K + (1 β)p L, β [0, 1] 8/23
10 Numerical Scheme: Discrete Regularization Term (1) Hydrodynamic step (continuity, mass fraction, internal energy) ρ n+1 K = ρ n K δt div K (ρ n u n ) ȳ n+1 K = 1 ρ n+1 K e n+1 K = 1 ρ n+1 K (2) Chemistry step [ρ n K y n K δt div K (ρ n u n y n )] [ ρ n K en K δt(div K (e n ρ n u n ) + pk n div K (u n ) + 1 ] K Sn K ) y n+1 K = ȳ n+1 K exp( κ(t n K )δt) with T n K = Mpn K Rρ n K (3) Equation of state p n+1 K = (γ 1)ρ n+1 K (en+1 K q 0 y n+1 K ) (4) Hydrodynamic step (momentum) u n+1 σ,i = 1 [ρ n ρ n+1 Dσ uσ,i n δt ( div Dσ (ρ n u n u n i ) + ( p n+1 ) ] ) σ,i D σ Regularization term S n K compensates upwind diffusion in time and space generated in the momentum balance: S µ x,t u 2 It can be proved: S > 0 e > 0 Under CFL, S > 0 S n K is the counterpart of residual term R n+1 σ from the kinetic energy balance. It is found by deriving a discrete kinetic energy balance from the momentum equation. Then S n K can be designed accordingly. 9/23
11 Numerical Scheme: Discrete Regularization Term (1) Hydrodynamic step (continuity, mass fraction, internal energy) ρ n+1 K = ρ n K δt div K (ρ n u n ) ȳ n+1 K = 1 ρ n+1 K e n+1 K = 1 ρ n+1 K (2) Chemistry step [ρ n K y n K δt div K (ρ n u n y n )] [ ρ n K en K δt(div K (e n ρ n u n ) + pk n div K (u n ) + 1 ] K Sn K ) y n+1 K = ȳ n+1 K exp( κ(t n K )δt) with T n K = Mpn K Rρ n K (3) Equation of state p n+1 K = (γ 1)ρ n+1 K (en+1 K q 0 y n+1 K ) (4) Hydrodynamic step (momentum) u n+1 σ,i = 1 [ρ n ρ n+1 Dσ uσ,i n δt ( div Dσ (ρ n u n u n i ) + ( p n+1 ) ] ) σ,i D σ Regularization term S n K compensates upwind diffusion in time and space generated in the momentum balance: S µ x,t u 2 It can be proved: S > 0 e > 0 Under CFL, S > 0 S n K is the counterpart of residual term R n+1 σ from the kinetic energy balance. It is found by deriving a discrete kinetic energy balance from the momentum equation. Then S n K can be designed accordingly. 9/23
12 Numerical Scheme: Discrete Regularization Term (1) Hydrodynamic step (continuity, mass fraction, internal energy) ρ n+1 K = ρ n K δt div K (ρ n u n ) ȳ n+1 K = 1 ρ n+1 K e n+1 K = 1 ρ n+1 K (2) Chemistry step [ρ n K y n K δt div K (ρ n u n y n )] [ ρ n K en K δt(div K (e n ρ n u n ) + pk n div K (u n ) + 1 ] K Sn K ) y n+1 K = ȳ n+1 K exp( κ(t n K )δt) with T n K = Mpn K Rρ n K (3) Equation of state p n+1 K = (γ 1)ρ n+1 K (en+1 K q 0 y n+1 K ) (4) Hydrodynamic step (momentum) u n+1 σ,i = 1 [ρ n ρ n+1 Dσ uσ,i n δt ( div Dσ (ρ n u n u n i ) + ( p n+1 ) ] ) σ,i D σ Regularization term S n K compensates upwind diffusion in time and space generated in the momentum balance: S µ x,t u 2 It can be proved: S > 0 e > 0 Under CFL, S > 0 S n K is the counterpart of residual term R n+1 σ from the kinetic energy balance. It is found by deriving a discrete kinetic energy balance from the momentum equation. Then S n K can be designed accordingly. 9/23
13 Numerical Scheme: Discrete Regularization Term S K defined to compensate the source term appearing in the following discrete kinetic energy balance (Nguyen, 2012) (Therme, 2015) δt ( ( ρ n+1 D σ u n+1 σ,i ) 2 ρ n Dσ ( u n σ,i ) 2 ) div D σ (ρu u 2 ) + ( p n+1 ) σ,i uσ,i n = 1 D σ Rn+1 σ,i with for all test function φ, when h, δt 0: ( ) d δt SK n φ n+1 K R n+1 σ,i φ n+1 σ 0 σ E S,i n 0 K M i=1 Moreover if the following discrete duality property holds, K M K p K div K uφ K + d D σ u σ,i ( p) σ,i φ σ σ E S,i i=1 Ω pu φ We recover a global total energy balance at the limit. 10/23
14 Numerical Scheme: Discrete Regularization Term 10/23 S K defined to compensate the source term appearing in the following discrete kinetic energy balance (Nguyen, 2012) (Therme, 2015) δt ( ( ρ n+1 D σ u n+1 σ,i ) 2 ρ n Dσ ( u n σ,i ) 2 ) div D σ (ρu u 2 ) + ( p n+1 ) σ,i uσ,i n = 1 D σ Rn+1 σ,i with for all test function φ, when h, δt 0: ( ) d δt SK n φ n+1 K R n+1 σ,i φ n+1 σ 0 σ E S,i n 0 K M The need for the discrete duality property constrains the choice for the discretization of the pressure gradient (centered): K M K p K div K u + d i=1 i=1 σ E S,i D σ u σ,i ( p) σ,i = 0
15 Numerical Scheme: Conserving Entropy Why solve ρ/e/p/u (1) instead of ρ/u/e/p (2)? Formally, the entropy balance can be recovered from the mass and internal energy balances as follows: s = s(ρ, e) = ln(ρ) 1 γ 1 ln e ρ s [ t ρ + div(ρu)] + es ρ [ t(ρe) + div(ρue) + p div u] ( = t s + u s ) p div u γ 1 ρe }{{} =0 With choice (1) we have time-level consistency between div(ρ n u n ) (from mass balance) and p n div u n (from internal energy). On the other hand the entropy balance cannot be recovered with choice (2) : u n+1 σ,i = 1 [ ρ n ρ n+1 Dσ uσ,i n δt (div Dσ (ρ n u n u n i ) + ( p n ) σ,i ) ] D σ e n+1 K = 1 ρ n+1 K [ ρ n K ek n δt(div K (e n ρ n u n ) + pk n div K (u n+1 ) + 1 ] K Sn K) 11/23
16 12/23 Numerical Scheme: Conserving Entropy Choice (1) vs choice (2) : non-conservation of entropy in rarefactions. Both yield a weak solution, but only the one from choice (2) is entropic. 1,000 1, pressure 500 pressure position position ρ u e p ρ e p u
17 Application: CJ Detonation 13/23 Exact ZND solution (q 0 = 25, Arrhenius law with κ 0 = 200, E A = 25) A N A 1/2 y A 1 A 1/2 30 A N A 0 p A x H C R A τ p A 1/2 A 1 A N A x
18 Application: CJ Detonation 13/23 Numerical solution, cells (400 cells/l 1/2 ) A N A 1/2 y A 1 A 1/2 30 A N A 0 p A x 15 C 40 A 1/2 A N 10 5 H R A τ p 20 A 1 A x
19 Application: CJ Detonation 13/23 Numerical solution, 300 cells (4 cells/l 1/2 ). 40 A N 1 A 1 A N 35 y A 0 p A x 15 C 40 A 1 A N 10 R p 20 5 A 0 H τ A x
20 Application: Detonation Diffraction 14/23 Detonation diffraction around a corner (Xu Aslam Stewart, 1997). u CJ unburnt state pressure horizontal axis Kinetic law: Arrhenius κ(t ) = κ 0 exp ( E ) A RT κ 0 = E A /M = 25 Initial data: Given unburnt state Matching ZND solution (CJ)
21 Application: Detonation Diffraction 15/23 Comparison of numerical results on a grid. 120 Isovalues of density at t = 14 (Xu Aslam Stewart, 1997). 120 Isovalues of density at t = 14 with our scheme (15 values)
22 Application: Detonation Diffraction 16/23 Comparison of numerical results on a grid. 120 Shock front at different times (solid line) : our results (red) vs (Xu Aslam Stewart, 1997) (black). 120 Isovalues y {0.1, 0.5, 0.9} (dashed) and y = 1 (solid line) for our results (red) vs (Xu Aslam Stewart, 1997) (black)
23 Application: Detonation in a Channel (Heaviside) Ozone decomposition detonation in a channel (Bao-Jin, 2000) Frame of reference : laboratory frame INFLOW burnt mixture ρ b, p b, u b, y b unburnt mixture ρ u, p u, u u, y u OUTFLOW Kinetic law: Heaviside κ(t ) = κ 0 H(T ) { 0 if T < Tign H(T ) = 1 otherwise T ign = 500 K κ 0 = s 1 Initial data: Given unburnt state Matching burnt state for a CJ detonation but with u b 2u bcj Symmetric deformation of the front into the unburnt state 17/23
24 Application: Detonation in a Channel (Heaviside) 18/23 Transverse displacement of the triple points ( cells) : density contours near the detonation front between t = 0 and t = (Bao Jin, 2000) Our results
25 Application: Detonation in a Channel (Heaviside) 18/23 Transverse displacement of the triple points ( cells) : density contours near the detonation front between t = 0 and t = (Bao Jin, 2000) Our results
26 Application: Detonation in a Channel (Arrhenius) Same problem, with the following changes Frame of reference : attached to shock (u s 1140 m s 1 ) OUTFLOW burnt/reacting mixture ρ b, p b, u b, y b unburnt mixture ρ u, p u, u u, y u INFLOW Kinetic law: Arrhenius κ(t ) = κ 0 H(T ) exp ( E ) A RT T ign = 500 K κ 0 = 10 8 s 1 R = J kg 1 K 1 E A = J kg 1 Initial data: Given unburnt state Matching ZND solution (CJ) Density: random perturbation r [0, 1] in small region after shock as ρ = ( r)ρ 19/23
27 20/23 Application: Detonation in a Channel (Arrhenius) Cellular structure: procedure for tracking triple points (max t0<t<t 1 p(x, t) in the laboratory frame) in parallel. Each rank r performs the following operations: 1 Define the domain of interest D at time t [t 0, t 1 ] y L y Pressure field in shock frame at time t > t 0 D = [x 0 σ(t), L x ] [0, L y ] σ(t) = ε 0 + (t t 0 )u s 0 σ(t) D x 0 L x x 2 Compute partitions of D over self and neighboring ranks y Pressure field in shock frame at time t > t 0 Ω = 0 i nrank Ω i for i = 0,..., n rank : Π i = D Ω i 0 D Π 3 Π 4 Π 5 Ω 1 Ω 2 Ω 3 Ω 4 Ω 5 x
28 Application: Detonation in a Channel (Arrhenius) Cellular structure: procedure for tracking triple points (max t0<t<t 1 p(x, t) in the laboratory frame) in parallel. Each rank r performs the following operations: 3 Determine sub-regions of Π r (local frame) to be sent and to which ranks Π i = Π i shifted by (x 0 σ(t), 0) for i = 0,..., n rank : Ω sendi = Π r Ω i 4 Determine from which ranks to receive Π r and the matching sub-regions (lab frame) for i = 0,..., n rank : Ω recvi = Π i Ω r y 0 Pressure field in laboratory frame at time t > t 0 Π 3 Π4 Π5 x 0 σ(t) Ω 1 Ω 2 Ω 3 Ω 4 Ω 5 Parallel communications at time t 5 Π5 Ω 3 3 Π3 Ω 1 Π4 Ω 3 Π4 Ω1 1 4 Π4 Ω 2 2 x 20/23
29 20/23 Application: Detonation in a Channel (Arrhenius) Cellular structure: procedure for tracking triple points (max t0<t<t 1 p(x, t) in the laboratory frame) in parallel. Each rank r performs the following operations: 5 Non-blocking Send/Recv and update of max t>t0 p(x, t) in the lab frame y Pressure field in shock frame at time t > t 0 σ(t) for i = 0,..., n rank : D if Ω sendi { } : Send p Ωsendi to rank i 0 x 0 x for i = 0,..., n rank : y max t>t0 p(x, t) (laboratory frame) if Ω recvi { } : Receive from rank i update max t 0<t<t 1 p(x, t) on Ω recvi 0 x 0 σ(t) Π 3 Π4 Π5 Ω 1 Ω 2 Ω 3 Ω 4 Ω 5 x
30 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = 0 Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
31 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = t 0 = 19700δt = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
32 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
33 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
34 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
35 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
36 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
37 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
38 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
39 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
40 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = Pressure Mass fraction max p(x, t) t 0 <t<t 1 21/23
41 Application: Detonation in a Channel (Arrhenius) Numerical solution, cells, t = t 1 = 45063δt = Pressure Mass fraction max p(x, t) t 0 <t<t 1 L C = /23
42 Application: Detonation in a Channel (Arrhenius) 22/23 Detonation front structure: (I)ncident shock (weak shock), (R)eflected shock (transverse wave), (M)ach stem and (T)riple point
43 Perspectives 23/23 1 Detailed chemistry 2 Interaction with more complex obstacles 3 3D unstable detonation 4 Coupling with codes for deflagration and low Mach flows 5 Higher order in space and time
On 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 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 informationUNIVERSITY OF CALGARY. Stability of Detonation Wave under Chain Branching Kinetics: Role of the Initiation Step. Michael Anthony Levy A THESIS
UNIVERSITY OF CALGARY Stability of Detonation Wave under Chain Branching Kinetics: Role of the Initiation Step by Michael Anthony Levy A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT
More informationAn adaptive projection scheme for all-mach flows applied to shock diffraction
An adaptive projection scheme for all-mach flows applied to shock diffraction CEA DEN/DM2S/STMF/LMEC AMU I2M UMR 7373/AA Group January 20, 2015 All-Mach solver? Approximate Riemann solvers for M 1 (Guillard-Viozat,
More informationCellular structure of detonation wave in hydrogen-methane-air mixtures
Open Access Journal Journal of Power Technologies 91 (3) (2011) 130 135 journal homepage:papers.itc.pw.edu.pl Cellular structure of detonation wave in hydrogen-methane-air mixtures Rafał Porowski, Andrzej
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 informationInitiation of stabilized detonations by projectiles
Initiation of stabilized detonations by projectiles P. Hung and J. E. Shepherd, Graduate Aeronautical Laboratory, California Institute of Technology, Pasadena, CA 91125 USA Abstract. A high-speed projectile
More informationDETONATION HAZARD CLASSIFICATION BASED ON THE CRITICAL ORIFICE PLATE DIAMETER FOR DETONATION PROPAGATION
DETONATION HAZARD CLASSIFICATION BASED ON THE CRITICAL ORIFICE PLATE DIAMETER FOR DETONATION PROPAGATION by Mitchell Cross A thesis submitted to the Department of Mechanical and Materials Engineering In
More informationTRANSITION TO DETONATION IN NON-UNIFORM H2-AIR: CHEMICAL KINETICS OF SHOCK-INDUCED STRONG IGNITION
TRANSITION TO DETONATION IN NON-UNIFORM H2-AIR: CHEMICAL KINETICS OF SHOCK-INDUCED STRONG IGNITION L.R. BOECK*, J. HASSLBERGER AND T. SATTELMAYER INSTITUTE OF THERMODYNAMICS, TU MUNICH *BOECK@TD.MW.TUM.DE
More informationOn two fractional step finite volume and finite element schemes for reactive low Mach number flows
Fourth International Symposium on Finite Volumes for Complex Applications - Problems and Perspectives - July 4-8, 2005 / Marrakech, Morocco On two fractional step finite volume and finite element schemes
More informationThe Deflagration Phase of Type Ia SNe
The Center for Astrophysical Thermonuclear Flashes The Deflagration Phase of Type Ia SNe Alan Calder ASC FLASH Center Type Ia Supernova Team Type Ia Supernovae and Cosmology August 5, 2004 An Advanced
More informationFundamentals of Rotating Detonation. Toshi Fujiwara (Nagoya University)
Fundamentals of Rotating Detonation Toshi Fujiwara (Nagoya University) New experimental results Cylindical channel D=140/150mm Hydrogen air; p o =1.0bar Professor Piotr Wolanski P [bar] 10 9 8 7 6 5 4
More informationDetonation Structure
Planar Detonations and Detonation Structure Jerry Seitzman. 5 Mole Fraction.5..5 CH4 HO HCO emerature Methane Flame...3 Distance (cm) 5 5 emerature (K) Detonations - Coyright 4-5 by Jerry M. Seitzman.
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 informationAME 513. " Lecture 8 Premixed flames I: Propagation rates
AME 53 Principles of Combustion " Lecture 8 Premixed flames I: Propagation rates Outline" Rankine-Hugoniot relations Hugoniot curves Rayleigh lines Families of solutions Detonations Chapman-Jouget Others
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 informationAn analytical model for direct initiation of gaseous detonations
Issw21 An analytical model for direct initiation of gaseous detonations C.A. Eckett, J.J. Quirk, J.E. Shepherd Graduate Aeronautical Laboratories, California Institute of Technology, Pasadena CA 91125,
More informationUnidirectional Propagation of Gas Detonations in Channels with Sawtooth Walls
Naval Research Laboratory Washington, DC 20375-5320 NRL/MR/6404--10-9255 Unidirectional Propagation of Gas Detonations in Channels with Sawtooth Walls Vadim N. Gamezo Elaine S. Oran Center for Reactive
More informationReaction Rate Closure for Turbulent Detonation Propagation through CLEM-LES
5 th ICDERS August 7, 05 Leeds, UK through CLEM-LES Brian Maxwell, Matei Radulescu Department of Mechanical Engineering, University of Ottawa 6 Louis Pasteur, Ottawa, KN 6N5, Canada Sam Falle School of
More informationINVESTIGATION OF INSTABILITIES AFFECTING DETONATIONS: IMPROVING THE RESOLUTION USING BLOCK-STRUCTURED ADAPTIVE MESH REFINEMENT PRASHAANTH RAVINDRAN
INVESTIGATION OF INSTABILITIES AFFECTING DETONATIONS: IMPROVING THE RESOLUTION USING BLOCK-STRUCTURED ADAPTIVE MESH REFINEMENT by PRASHAANTH RAVINDRAN Presented to the Faculty of the Graduate School of
More informationCombustion Behind Shock Waves
Paper 3F-29 Fall 23 Western States Section/Combustion Institute 1 Abstract Combustion Behind Shock Waves Sandeep Singh, Daniel Lieberman, and Joseph E. Shepherd 1 Graduate Aeronautical Laboratories, California
More informationDetonation Diffraction
Detonation Diffraction E. Schultz, J. Shepherd Detonation Physics Laboratory Pasadena, CA 91125 MURI Mid-Year Pulse Detonation Engine Review Meeting February 10-11, 2000 Super-critical Detonation Diffraction
More informationPeriodic oscillation and fine structure of wedge-induced oblique detonation waves
Acta Mech. Sin. (211) 27(6):922 928 DOI 1.17/s149-11-58-y RESEARCH PAPER Periodic oscillation and fine structure of wedge-induced oblique detonation waves Ming-Yue Gui Bao-Chun Fan Gang Dong Received:
More informationTHERMODYNAMIC ANALYSIS OF COMBUSTION PROCESSES FOR PROPULSION SYSTEMS
2nd AIAA Aerospace Sciences Paper 2-33 Meeting and Exhibit January -8, 2, Reno, NV THERMODYNAMIC ANALYSIS OF COMBUSTION PROCESSES FOR PROPULSION SYSTEMS E. Wintenberger and J. E. Shepherd Graduate Aeronautical
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 informationSIMULATION OF DETONATION AFTER AN ACCIDENTAL HYDROGEN RELEASE IN ENCLOSED ENVIRONMENTS
SIMULATION OF DETONATION AFTER AN ACCIDENTAL HYDROGEN RELEASE IN ENCLOSED ENVIRONMENTS L. Bédard Tremblay, 1 L. Fang, 1 L. Bauwens,1 P.H.E. Finstad,2 Z. Cheng3 and A. V. Tchouvelev3 1 Department of Mechanical
More informationBIFURCATING MACH SHOCK REFLECTIONS WITH APPLICATION TO DETONATION STRUCTURE
BIFURCATING MACH SHOCK REFLECTIONS WITH APPLICATION TO DETONATION STRUCTURE Philip Mach A thesis submitted to the Faculty of Graduate and Postdoctoral Studies in partial fulfillment of the requirements
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 informationPropagation of hydrogen-air detonation in tube with obstacles
Open Access Journal Journal of Power Technologies 91 (3) (2011) 122 129 journal homepage:papers.itc.pw.edu.pl Propagation of hydrogen-air detonation in tube with obstacles Wojciech Rudy, Rafał Porowski,
More informationDetonation of Gas Particle Flow
2 Detonation of Gas Particle Flow F. Zhang 2.1 Introduction Fine organic or metallic particles suspended in an oxidizing or combustible gas form a reactive particle gas mixture. Explosion pressures in
More informationTriple Point Collision and Origin of Unburned Gas Pockets in Irregular Detonations
Triple Point Collision and Origin of Unburned Gas Pockets in Irregular Detonations Yasser Mahmoudi 1*, Kiumars Mazaheri 1 Department of Engineering, University of Cambridge, CB 1PZ Cambridge, United Kingdom
More informationAngular momentum preserving CFD on general grids
B. Després LJLL-Paris VI Thanks CEA and ANR Chrome Angular momentum preserving CFD on general grids collaboration Emmanuel Labourasse (CEA) B. Després LJLL-Paris VI Thanks CEA and ANR Chrome collaboration
More informationAA210A Fundamentals of Compressible Flow. Chapter 13 - Unsteady Waves in Compressible Flow
AA210A Fundamentals of Compressible Flow Chapter 13 - Unsteady Waves in Compressible Flow The Shock Tube - Wave Diagram 13.1 Equations for irrotational, homentropic, unsteady flow ρ t + x k ρ U i t (
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 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 informationDevelopment of One-Step Chemistry Models for Flame and Ignition Simulation
Development of One-Step Chemistry Models for Flame and Ignition Simulation S.P.M. Bane, J.L. Ziegler, and J.E. Shepherd Graduate Aerospace Laboratories California Institute of Technology Pasadena, CA 91125
More informationConsistency analysis of a 1D Finite Volume scheme for barotropic Euler models
Consistency analysis of a 1D Finite Volume scheme for barotropic Euler models F Berthelin 1,3, T Goudon 1,3, and S Mineaud,3 1 Inria, Sophia Antipolis Méditerranée Research Centre, Proect COFFEE Inria,
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 informationNUMERICAL SIMULATION OF UNSTEADY NORMAL DETONATION COMBUSTION AJJAY OMPRAKAS. Presented to the Faculty of the Graduate School of
NUMERICAL SIMULATION OF UNSTEADY NORMAL DETONATION COMBUSTION by AJJAY OMPRAKAS Presented to the Faculty of the Graduate School of The University of Texas at Arlington in Partial Fulfillment of the Requirements
More informationDetonation Characteristics Of Dimethyl Ether, Methanol and Ethanol Air Mixtures
24 th ICDERS July 28 - August 2, 213 Taipei, Taiwan Detonation Characteristics Of Dimethyl Ether, Methanol and Ethanol Air Mixtures Peter Diakow and Gaby Ciccarelli Queen's University Kingston, Ontario,
More informationEVALUATION OF HYDROGEN, PROPANE AND METHANE-AIR DETONATIONS INSTABILITY AND DETONABILITY
EVALUATION OF HYDROGEN, PROPANE AND METHANE-AIR DETONATIONS INSTABILITY AND DETONABILITY Borzou, B. 1 and Radulescu, M.I. 2 1 Mechanical Engineering Department, University of Ottawa, 161 Louis Pasteur,
More informationRayleigh-Taylor Unstable Flames
Rayleigh-Taylor Unstable Flames Elizabeth P. Hicks 1,2 and Robert Rosner 2 CIERA, Northwestern University 1 University of Chicago 2 CIERA Conference: September 2, 2011 1 Type Ia Supernovae Image: NASA
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 informationTriple point shear-layers in gaseous detonation waves
4nd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit 9 - July 6, Sacramento, California AIAA 6-575 Triple point shear-layers in gaseous detonation waves L. Massa, J.M. Austin, T.L. Jackson, University
More informationA numerical study of detonation diffraction
J. Fluid Mech., 529:117-146, 2005. (Preprint - see journal for final version http://dx.doi.org/10.1017/s0022112005003319). Under consideration for publication in J. Fluid Mech. 1 A numerical study of detonation
More informationMild Ignition Phenomena in Rapid Compression Machines
25 th ICDERS August 2 7, 2015 Leeds, UK Kevin P. Grogan a, S. Scott Goldsborough b, Matthias Ihme a a Stanford University, Stanford, CA 94305 b Argonne National Laboratory, Argonne, IL 60439 1 Introduction
More informationNumerical Hydraulics
ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction
More informationCarbon Science and Technology
ASI RESEARCH ARTICLE Carbon Science and Technology Received:10/03/2016, Accepted:15/04/2016 ------------------------------------------------------------------------------------------------------------------------------
More informationPresentation Start. Zero Carbon Energy Solutions 4/06/06 10/3/2013:; 1
Presentation Start 10/3/2013:; 1 4/06/06 What is an Explosion? Keller, J.O. President and CEO,, ISO TC 197, Technical Program Director for the Built Environment and Safety; Gresho, M. President, FP2FIRE,
More informationEffect of Initial Disturbance on The Detonation Front Structure of a Narrow Duct. Hua-Shu Dou *, Boo Cheong Khoo,
Published in: Shock Waves, Vol.20, No.2, March 2010, pp.163-173. Effect of Initial Disturbance on The Detonation Front Structure of a Narrow Duct Hua-Shu Dou *, Boo Cheong Khoo, Department of Mechanical
More informationShock and Expansion Waves
Chapter For the solution of the Euler equations to represent adequately a given large-reynolds-number flow, we need to consider in general the existence of discontinuity surfaces, across which the fluid
More informationQuenching and propagation of combustion fronts in porous media
Quenching and propagation of combustion fronts in porous media Peter Gordon Department of Mathematical Sciences New Jersey Institute of Technology Newark, NJ 72, USA CAMS Report 56-6, Spring 26 Center
More informationarxiv: v1 [physics.flu-dyn] 7 Dec 2015
Chapman-Jouguet deflagrations and their transition to detonation Mohamed Saif a, Wentian Wang a, Andrzej Pekalski b, Marc Levin c, Matei I. Radulescu a a Department of Mechanical Engineering, University
More informationarxiv: v1 [physics.flu-dyn] 9 Aug 2017
arxiv:1708.02682v1 [physics.flu-dyn] 9 Aug 2017 Investigation of Turbulent Mixing and Local Reaction Rates on Deflagration to Detonation Transition in Methane-Oxygen Brian Maxwell a,c,, Andrzej Pekalski
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More informationThe evolution and cellular structure of a detonation subsequent to a head-on interaction with a shock wave
Combustion and Flame 151 (2007) 573 580 www.elsevier.com/locate/combustflame The evolution and cellular structure of a detonation subsequent to a head-on interaction with a shock wave Barbara B. Botros
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 informationSimulation of sympathetic detonation by a CIP Eulerian code
Computational Ballistics II 107 Simulation of sympathetic detonation by a CIP Eulerian code S. Kubota 1, Z. Liu 2, T. Saburi 1, Y. Ogata 1 & M. Yoshida 1 1 Research Center for Explosion Safety, National
More informationAn artificial nonlinear diffusivity method for shock-capturing in supersonic reacting flows
Center for Turbulence Research Annual Research Briefs 005 57 An artificial nonlinear diffusivity method for shock-capturing in supersonic reacting flows By B. Fiorina AND S. K. Lele. Motivation and objectives
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 informationExperimental Study of 2D-Instabilities of Hydrogen Flames in Flat Layers
25 th ICDERS August 2 7, 2015 Leeds, UK Experimental Study of 2D-Instabilities of Hydrogen Flames in Flat Layers M. Kuznetsov 1 *, J. Grune 2, S. Tengah 1, J. Yanez 1 1 Intitute for Energy and Nuclear
More informationNear limit behavior of the detonation velocity
Near limit behavior of the detonation velocity John H.S. Lee 1, Anne Jesuthasan 1 and Hoi Dick Ng 2 1 McGill University Department of Mechanical Engineering Montreal, QC, Canada 2 Concordia University
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 informationAcoustic timescale characterization of hot spot detonability
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 214 Acoustic timescale characterization of hot spot detonability Michael D. Kurtz Iowa State University Follow
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 informationFLAME ACCELERATION AND TRANSITION FROM DEFLAGRATION TO DETONATION IN HYDROGEN EXPLOSIONS
FLAME ACCELERATION AND TRANSITION FROM DEFLAGRATION TO DETONATION IN HYDROGEN EXPLOSIONS A. Heidari and J.X. Wen* Centre for Fire and Explosion Studies, Faculty of Engineering, Kingston University Friars
More informationThe role of diffusion at shear layers in irregular detonations
The role of diffusion at shear layers in irregular detonations Marco Arienti 1 Joseph E. Shepherd 2 1 United Technologies Research Center, 411 Silver Lane, East Hartford, CT 06108 2 California Institute
More informationEQUIVALENCE OF LOW-FREQUENCY STABILITY CONDITIONS FOR MULTIDIMENSIONAL DETONATIONS IN THREE MODELS OF COMBUSTION
EQUIVALENCE OF LOW-FREQUENCY STABILITY CONDITIONS FOR MULTIDIMENSIONAL DETONATIONS IN THREE MODELS OF COMBUSTION HELGE KRISTIAN JENSSEN, GREGORY LYNG, MARK WILLIAMS Abstract. We use the classical normal
More informationDETONATION WAVE PROPAGATION IN SEMI-CONFINED LAYERS OF HYDROGEN-AIR AND HYDROGEN-OXYGEN MIXTURES
DETONATION WAVE PROPAGATION IN SEMI-CONFINED LAYERS OF HYDROGEN-AIR AND HYDROGEN-OXYGEN MIXTURES Grune, J. 1 *, Sempert, K. 1, Friedrich, A. 1, Kuznetsov, M. 2, Jordan, T. 2 1 Pro-Science GmbH, Parkstr.9,
More informationHigh-order finite difference methods with subcell resolution for 2D detonation waves
Center for Turbulence Research Annual Research Briefs 7 High-order finite difference methods with subcell resolution for D detonation waves B W. Wang, C.-W. Shu, H. C. Yee AND B. Sjögreen. Motivation and
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 informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
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 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 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 informationIX. COMPRESSIBLE FLOW. ρ = P
IX. COMPRESSIBLE FLOW Compressible flow is the study of fluids flowing at speeds comparable to the local speed of sound. This occurs when fluid speeds are about 30% or more of the local acoustic velocity.
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 informationDetonation Front Structure and the Competition for Radicals
Detonation Front Structure and the Competition for Radicals Z. Liang, S. Browne, R. Deiterding, and J. E. Shepherd California Institute of Technology, Pasadena, CA 9112 USA Abstract We examine the role
More informationSevere, unconfined petrol vapour explosions
Severe, unconfined petrol vapour explosions Graham Atkinson Health and Safety Laboratory Fire and Process Safety Unit Buncefield - Lessons for whom? Plant managers Safety managers Risk assessors Explosion
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 informationA REDUCED-ORDER METHANE-AIR COMBUSTION MECHANISM THAT SATISFIES THE DIFFERENTIAL ENTROPY INEQUALITY
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Special Issue/2018, pp. 285 290 A REDUCED-ORDER METHANE-AIR COMBUSTION MECHANISM THAT SATISFIES THE DIFFERENTIAL
More informationCFD Analysis of Mixing in Polymerization Reactor. By Haresh Patel Supervisors: Dr. R. Dhib & Dr. F. Ein-Mozaffari IPR 2007
CFD Analysis of Mixing in Polymerization Reactor By Haresh Patel Supervisors: Dr. R. Dhib & Dr. F. Ein-Mozaffari Introduction Model development Simulation Outline Model Setup for Fluent Results and discussion
More informationLaminar Premixed Flames: Flame Structure
Laminar Premixed Flames: Flame Structure Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Course Overview Part I: Fundamentals and Laminar Flames Introduction Fundamentals and mass balances of
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationEQUATION OF STATE FOR MODELING THE DETONATION REACTION ZONE
EQUATION OF STATE FOR MODELING THE DETONATION REACTION ZONE D. Scott Stewart and Sunhee Yoo Department of Theoretical and Applied Mechanics University of Illinois, Urbana, IL 61801, USA William C. Davis
More informationTHRUST CHAMBER DYNAMICS AND PROPULSIVE PERFORMANCE OF AIRBREATHING PULSE DETONATION ENGINES
The Pennsylvania State University The Graduate School Department of Mechanical and Nuclear Engineering THRUST CHAMBER DYNAMICS AND PROPULSIVE PERFORMANCE OF AIRBREATHING PULSE DETONATION ENGINES A Thesis
More informationNumerical simulation of air blast waves
Numerical simulation of air blast waves M. Arrigoni, S. Kerampran, ENSTA Bretagne, France J.-B. Mouillet, Altair Engineering France B. Simoens, M. Lefebvre, S. Tuilard, Ecole Royal Militaire de Bruxelles,
More informationHot Surface Ignition and Flow Separation
25 th ICDERS August 2 7, 215 Leeds, UK Hot Surface Ignition and Flow Separation Melguizo-Gavilanes, J. and Shepherd, J.E. Explosion Dynamics Laboratory Graduate Aerospace Laboratories of the California
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
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 informationChemical reactors. H has thermal contribution, pressure contribution (often negligible) and reaction contribution ( source - like)
Chemical reactors - chemical transformation of reactants into products Classification: a) according to the type of equipment o batch stirred tanks small-scale production, mostly liquids o continuous stirred
More informationGHS: Physical Hazards
A liquid having a flash point of not more than 93 C (199.4 F). 3 4 Extremely Liquid and Vapor Highly Liquid and Vapor Liquid and Vapor Combustible Liquid Liquids Flash point < 23 C (73.4 F) and initial
More informationWell Stirred Reactor Stabilization of flames
Well Stirred Reactor Stabilization of flames Well Stirred Reactor (see books on Combustion ) Stabilization of flames in high speed flows (see books on Combustion ) Stabilization of flames Although the
More informationNumerical Simulation of Premixed V-Flame
Proceedings of the World Congress on Engineering 7 Vol II WCE 7, July -, 7, London, U.K. Numerical Simulation of Premixed V-Flame C.K. Chan, B. Stewart, and C.W. Leung Abstract A Lagrangian front-tracking
More informationNumerical and Analytical Studies of the Dynamics of Gaseous Detonations
Numerical and Analytical Studies of the Dynamics of Gaseous Detonations Thesis by Christopher A. Eckett In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy I A I N S T I T
More informationNumerical simulation of a rotating detonation with a realistic injector designed for separate supply of gaseous hydrogen and oxygen
Numerical simulation of a rotating detonation with a realistic injector designed for separate supply of gaseous hydrogen and oxygen Thomas Gaillard, Dmitry Davidenko, Francis Dupoirieux To cite this version:
More information