-
Notifications
You must be signed in to change notification settings - Fork 0
/
specHeatC12H26.m
32 lines (22 loc) · 1 KB
/
specHeatC12H26.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
function [cp, cv, k] = specHeatC12H26(AF, T)
%
% Returns dodecane-air fuel mixture specific heats (see eqn below)
%
% phi*C12H26 + 18.5*(O2 + 3.76*N2) ->
% 12*phi*CO2 + 13*phi*H2O + (3.76*18.5)*N2 + 18.5*(1-phi)*O2
%
% Compute mass fractions
AF_stoich = 18.5*(const.M_O2 + ... % stoichometric air-fuel ratio
3.76*const.M_N2)/const.M_C12H26;
phi = AF_stoich/AF; % equivalence ratio
M_total = 12*phi + 13*phi + 3.76*18.5 + 18.5*(1-phi);
MF_CO2 = 12*phi/M_total; % molar fraction CO2
MF_H2O = 13*phi/M_total; % molar fraction H2O
MF_N2 = 3.76*18.5/M_total; % molar fraction N2
MF_O2 = 18.5*(1-phi)/M_total; % molar fraction O2
MF = [MF_CO2 MF_H2O MF_N2 MF_O2]'; % molar fractions
M = [const.M_CO2; const.M_H2O; const.M_N2; const.M_O2]; % molar mass
ABCD = [const.abcd_CO2; const.abcd_H2O; ...
const.abcd_N2; const.abcd_O2]; % specific heat constants (3-poly)
[cp, cv, k] = specHeatMix(M, MF, ABCD, T);
end