1% -------------------------------------------------------------------------
2% EXAMPLE: HP
3%
4% Compute adiabatic temperature and equilibrium composition at constant
5% pressure (e.g., 1.01325 bar) for lean to rich CH4-air mixtures at
6% standard conditions, a set of 26 species considered and a set of
7% equivalence ratios phi contained in (0.5, 5) [-]
8%
9% Soot formation == {'CO2','CO','H2O','H2','O2','N2','He','Ar','Cbgrb',...
10% 'C2','C2H4','CH','CH','CH3','CH4','CN','H',...
11% 'HCN','HCO','N','NH','NH2','NH3','NO','O','OH'}
12%
13% See wiki or list_species() for more predefined sets of species
14%
15% @author: Alberto Cuadra Lara
16% PhD Candidate - Group Fluid Mechanics
17% Universidad Carlos III de Madrid
18%
19% Last update July 22 2022
20% -------------------------------------------------------------------------
21
22%% INITIALIZE
23self = App('Soot formation');
24%% INITIAL CONDITIONS
25self = set_prop(self, 'TR', 300, 'pR', 1 * 1.01325, 'phi', 0.5:0.01:5);
26self.PD.S_Fuel = {'CH4'};
27self.PD.S_Oxidizer = {'N2', 'O2', 'Ar', 'CO2'};
28self.PD.ratio_oxidizers_O2 = [78.084, 20.9476, 0.9365, 0.0319] ./ 20.9476;
29%% ADDITIONAL INPUTS (DEPENDS OF THE PROBLEM SELECTED)
30self = set_prop(self, 'pP', self.PD.pR.value);
31%% SOLVE PROBLEM
32self = solve_problem(self, 'HP');
33%% DISPLAY RESULTS (PLOTS)
34post_results(self);