% Shock and Detonation Toolboox
% http://www.galcit.caltech.edu/EDL/public/cantera/html/SD_Toolbox/
%
% Calculate the CJSpeed using the Minimum Wave Speed Method.

clear;clc;
display('demo_CJ')

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

gas = importPhase(mech);
set(gas,'Temperature',T1,'Pressure',P1,'MoleFractions',q);

[cj_speed, curve, fit_goodness, dratio_cj] = CJspeed(P1, T1, q, mech, fig_num);

display(['CJ computation for ',mech,' with composition ',q])
display(['CJ speed ',num2str(cj_speed),' (m/s)']);
display(['Use demo_CJstate to find thermodynamic properties at CJ point']);