Example_HP_PROPELLANTS.m#

 1% -------------------------------------------------------------------------
 2% EXAMPLE: HP PROPELLANTS
 3%
 4% Compute adiabatic temperature and equilibrium composition at constant
 5% pressure (e.g., 1.01325 bar) for lean to rich LH2-LOX mixtures at
 6% standard conditions, a set of 24 species considered and a set of
 7% equivalence ratios phi contained in (0.5, 5) [-]
 8%   
 9% HYDROGEN_L == {'H','H2O','OH','H2','O','O3','O2','HO2','H2O2',...
10%                'H2bLb','O2bLb'}
11%   
12% See wiki or ListSpecies() for more predefined sets of species
13%
14% @author: Alberto Cuadra Lara
15%          PhD Candidate - Group Fluid Mechanics
16%          Universidad Carlos III de Madrid
17%                 
18% Last update Feb 19 2022
19% -------------------------------------------------------------------------
20
21%% INITIALIZE
22self = App('HYDROGEN_L');
23%% INITIAL CONDITIONS
24self = set_prop(self, 'TR', 300, 'pR', 1 * 1.01325, 'phi', 0.2:0.05:5);
25self.PD.S_Fuel     = {'H2bLb'};
26self.PD.S_Oxidizer = {'O2bLb'};
27%% ADDITIONAL INPUTS (DEPENDS OF THE PROBLEM SELECTED)
28self = set_prop(self, 'pP', self.PD.pR.value); 
29%% SOLVE PROBLEM
30self = SolveProblem(self, 'HP');
31%% DISPLAY RESULTS (PLOTS)
32postResults(self);