Shock & Detonation Toolbox - Matlab

• Shock & Detonation Toolbox Matlab Demos

Constant Volume Explosion
Initial Conditions: Shock Heated Mixture
(CJ Speed)

• demo_cvCJ.m
  Matlab program used to generate the sample output described below
  
  
• Numerical Solution Methods for Control Volume Explosions and ZND Detonation Structure
  Documentation discussing algorithm used for this demo

Input Data
Mixture	Stoichiometric H2/O2 'H2:2 O2:1 N2:3.76'
Mechanism	h2air_highT.cti
Initial Pressure	1 bar = 100000 Pa
Initial Temperature	300 K
Output Data
CJ speed	1967.8451 m/s
gas object at equilibrium
Output file

The CJ velocity is the minimum of the plot	Curve fit to find the minimum of the data

Constant Volume Temperature & Pressure Profiles
& Induction Time (tau) calculated as the time to the maximum temperature gradient.