numerical simulation of detonation re-initiation following mach reflection
DESCRIPTION
Numerical simulation of detonation re-initiation following Mach reflection. She-Ming Lau-Chapdelaine's Graduate Seminar, October 2012. Department of Mechanical Engineering, University of Ottawa. All rights reserved.TRANSCRIPT
![Page 1: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/1.jpg)
Numerical simulation of detonation re-initiation following Mach reflection
S. She-Ming Lau-ChapdelaineRohit BhattacharjeeMatei I. Radulescu
2012/10/26
![Page 2: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/2.jpg)
Combustion
● Channel filled with premixed gas
● Deflagration● Detonation
![Page 3: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/3.jpg)
Detonation
● Supersonic combustion wave
● Shock wave heats mixture● Mixture reacts● Energy release propels shock wave
![Page 4: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/4.jpg)
Detonation Applications
● Pulse detonation engine
![Page 5: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/5.jpg)
Detonation Applications
● Rotating detonation engine
![Page 6: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/6.jpg)
Detonation Application
● Dust/powder explosion
Lakeland Mills sawmill, Prince George, B.C., April 2012
![Page 7: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/7.jpg)
Detonation Applications
● Nuclear safety - hydrogen explosions
Fukushima, 2011
Chernobyl, 1986
![Page 8: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/8.jpg)
Detonation
● Supersonic combustion wave
● Shock wave heats mixture● Mixture reacts● Energy release propels shock wave
![Page 9: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/9.jpg)
Detonation
● Zeldovich-von Neumann-Doring
![Page 10: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/10.jpg)
Detonation structure
● ZND structure is unstable
![Page 11: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/11.jpg)
Detonation structure
● ZND structure is unstable
plane of symmetry
Shock Reflection
Unburnt Gas
![Page 12: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/12.jpg)
Detonation re-initiation
● Importance of rapid reactions mechanisms in detonations difficult to determine
● Rapid reaction mechanisms– 1) Transverse wave
– 2) Mach-stem
– 3) “Wall” jetting effect
– 4) Kelvin-Helmholtz
– 5) Richtmyer-Meshkov
1
2
34
5
![Page 13: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/13.jpg)
Previous work
Teodorczyk, A., J.H.S. Lee, and R. Knystautas. 1991. Prog. Astronaut. and Aeronaut. 133:223–240.
![Page 14: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/14.jpg)
Previous work
T. Obara, J. Sentanuhady, Y. Tsukada, S. Ohyagi, Reinitiation process of detonation wave behind a slit-plate, Shock Waves 18 (2) (2008) 117–127.
![Page 15: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/15.jpg)
Previous work
R. Bhattacharjee, S.SM. Lau-Chapdelaine, G. Maines, L. Maley, M.I. Radulescu. Detonation re-initiation following the Mach reflection of a quenched detonation. Proceedings of the International Combustion Symposium, (2012).
![Page 16: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/16.jpg)
Previous work
R. Bhattacharjee, S.SM. Lau-Chapdelaine, G. Maines, L. Maley, M.I. Radulescu. Detonation re-initiation following the Mach reflection of a quenched detonation. Proceedings of the International Combustion Symposium, (2012).
![Page 17: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/17.jpg)
Objectives
● Model detonation re-initiation● Isolate re-ignition mechanisms● Gain insight from simulations● Predict detonation re-initiation
![Page 18: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/18.jpg)
Numerical Model
● Reactive Euler equations
● 1-step Arrhenius chemistry
∂ρ
∂ t+
∂ρ u∂ x
+v∂ρ v∂ y
=0
ρu t
+
x(ρu2
+ p )+
y(ρ vu)=0
ρ v t
+
x(ρuv )+
y(ρv2
+ p)=0
t(E+ Q λ)+
x((E+ p+ Q λ)u)+
y((E+ p+ Q λ)v)=0
∂ λ∂ t
+u ∂λ∂ x
+v ∂ λ∂ y
=k (λ−1)eE a
RT
![Page 19: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/19.jpg)
Numerical Model
● Reactive Euler model● 1-step Arrhenius reaction● Non-dimensionalized by half-reaction length
and initial conditions∂ λ∂ t
+u ∂λ∂ x
+v ∂ λ∂ y
=k (λ−1)eE a
RT
Pre-exponential Factor
![Page 20: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/20.jpg)
Numerical Model
● Reactive Euler model● 1-step Arrhenius reaction● Non-dimensionalized● Calibrated for post-shock conditions (CH4 + 2O2)
λ
t+ u
λ
x+ v
λ
y=k (λ−1)e
Ea
RT 0
Activation Energy
![Page 21: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/21.jpg)
Numerical Model: Chemistry (CH4 + 2O2)
τ ig∝eEa
RT
![Page 22: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/22.jpg)
Numerical Model
● Reactive Euler model● 1-step Arrhenius reaction● Non-dimensionalized● Calibrated for post-shock conditions (CH4 + 2O2)
t(E+ Q λ)+
x((E+ P+ Q λ)u)+
y((E+ P+ Q λ)v )=0
Heat Release
![Page 23: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/23.jpg)
Numerical Model
● Reactive Euler model● 1-step Arrhenius reaction● Non-dimensionalized by half-reaction length● Calibrated for post-shock conditions (CH4 + 2O2)
● Domain
– Base grid: 200x24– Adaptive grid refinement technique– Resolution: >32 points per induction length
γ=1.17 ;EaRT 0
=48.3 ;QRT 0
=60.5
![Page 24: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/24.jpg)
Numerical Model: DomainP
ost
-ZN
DE
xtrapolate
Reflect
Reflect
![Page 25: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/25.jpg)
Numerical Model
● Reactive Euler model● 1-step Arrhenius reaction● Non-dimensionalized by half-reaction length● Calibrated for post-shock conditions (CH4 + 2O2)
● Domain
– Base grid: 200x24– Adaptive grid refinement technique– Resolution: >32 points per induction length
● AMRITA CFD
γ=1.17 ;EaRT 0
=48.3 ;QRT 0
=60.5
![Page 26: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/26.jpg)
Results: P0=5.5kPa
![Page 27: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/27.jpg)
Results: P0=5.5kPa
Quenching
![Page 28: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/28.jpg)
Results: P0=10.3kPa
![Page 29: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/29.jpg)
Results: P0=10.3kPa
Hot-spot re-ignition
![Page 30: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/30.jpg)
Results: P0=11.9kPa
![Page 31: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/31.jpg)
Results: P0=11.9kPa
Detonation re-initiation
![Page 32: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/32.jpg)
Results: P0=11.9kPa
Density
Reaction progress
![Page 33: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/33.jpg)
Results: P0=12.5kPa
Experimental results from 11.9kPa
Detonation re-initiation
![Page 34: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/34.jpg)
Results: P0=12.5kPa
![Page 35: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/35.jpg)
Results: P0=17.6kPa
![Page 36: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/36.jpg)
Results: P0=17.6kPa
Detonation transmission
![Page 37: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/37.jpg)
Summary
P0 Experimental Simulation5.5kPa Quenched
10.3kPa Hot-spot ignition behind Mach stem
11.9kPa Detonation re-initiation along Mach stem & transverse wave
Detonation re-initiation along Mach stem
17.6kPa Detonation direct transmission
![Page 38: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/38.jpg)
Summary
Detonation StructureExperiments
Detonation diffraction over obstacle
Numerical simulations
P0 Results
5.5kPa Quenched
10.3kPa Hot-spot
11.9kPa Re-initation
17.6kPa Transmission
Numerical modelof experiment
Reactive Euler
1-step Arrhenius
Post-shock calib.
>32 grid/induc.
![Page 39: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/39.jpg)
Conclusions
● Adiabatic Mach compression appears to play an important role
● Wall jetting could be significant
● Transverse wave not recreated in simulations
● Experiments show rapid combustion in unburnt tongue, strong effect of wall jet
![Page 40: Numerical simulation of detonation re-initiation following Mach reflection](https://reader034.vdocument.in/reader034/viewer/2022050808/55327564550346c3558b4618/html5/thumbnails/40.jpg)
Acknowledgements
Alexander Graham Bell CGS