Example_EV.m#

 1% -------------------------------------------------------------------------
 2% EXAMPLE: EV
 3%
 4% Compute equilibrium composition at adiabatic temperature and constant
 5% volume for lean to rich CH4-air mixtures at standard conditions, a set
 6% of 26 species considered and a set of equivalence ratios (phi) contained
 7% 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 ListSpecies() 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 Oct 22 2021
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 = {'O2'};
28self.PD.S_Inert    = {'N2', 'Ar', 'CO2'};
29self.PD.proportion_inerts_O2 = [78.084, 0.9365, 0.0319] ./ 20.9476;
30%% ADDITIONAL INPUTS (DEPENDS OF THE PROBLEM SELECTED)
31% No additional data required. The internal energy and volume are constants.
32%% SOLVE PROBLEM
33self = SolveProblem(self, 'EV');
34%% DISPLAY RESULTS (PLOTS)
35postResults(self);