1% -------------------------------------------------------------------------
2% EXAMPLE: DET REFLECTED
3%
4% Compute pre-shock and post-shock state for a reflected planar detonation
5% considering Chapman-Jouguet (CJ) theory for lean to rich CH4-air mixtures
6% at standard conditions (T = 300 K, p = 1 atm), and a set of equivalence
7% ratios (phi) contained in (0.5, 5) [-]
8%
9% See wiki or setListspecies method from ChemicalSystem class for predefined
10% sets of species
11%
12% @author: Alberto Cuadra Lara
13%
14% Last update October 06 2025
15% -------------------------------------------------------------------------
16
17% Import packages
18import combustiontoolbox.databases.NasaDatabase
19import combustiontoolbox.core.*
20import combustiontoolbox.shockdetonation.*
21
22% Get Nasa database
23DB = NasaDatabase();
24
25% Define chemical system
26system = ChemicalSystem(DB);
27
28% Initialize mixture
29mix = Mixture(system);
30
31% Define chemical state
32set(mix, {'CH4'}, 'fuel', 1);
33set(mix, {'N2', 'O2', 'Ar', 'CO2'}, 'oxidizer', [78.084, 20.9476, 0.9365, 0.0319] / 20.9476);
34
35% Define properties
36mixArray1 = setProperties(mix, 'temperature', 300, 'pressure', 1.01325, 'equivalenceRatio', 0.5:0.01:5);
37
38% Initialize solver
39solver = DetonationSolver('problemType', 'DET_R');
40
41% Solve problem
42[mixArray1, mixArray2, mixArray3] = solver.solveArray(mixArray1);
43
44% Generate report
45report(solver, mixArray1, mixArray2, mixArray3);