% Shock and Detonation Toolboox
% http://www.galcit.caltech.edu/EDL/public/cantera/html/SD_Toolbox/ 
%
% Generates plots and output files for a constant volume explosion simulation where  
% the initial conditions are shocked reactants and the shock travels at the
% speed of a CJ detonation in the given mixture.

clear;clc;
display('demo_cvCJ')

P1 = 100000; T1 = 300; 
q = 'H2:2 O2:1 N2:3.76';
mech = 'h2air_highT.cti';

fig_num1 = 1; % figure number for Temperature profile
fig_num2 = 0; % for CJspeed plots

[cj_speed,gas] = cv_CJ(fig_num1, fig_num2, P1, T1, q, mech, 'h2air');