Project 1 - Parsing NASA thermodynamic data
AIM: To write the code in MATLAB to parse the NASA thermodynamic file.
OBJECTIVE-
1. To write a function that extracts 14 coefficient for a given species from the given information.
2. To write functions that calculate the enthalpy, entropy, and specific heat using NASA polynomials for the given species.
3. To write a function that plots the C, Enthalpy, and Entropy of input species.
4. To write a function that calculates the molecular weight of each species.
INTRODUCTION
File parsing is done to extract data from the file, it should be generic as possible to fetch the data and by using this we calculate enthalpy, entropy & specific heat for different species and condition is that the first 7 coefficients are high-temperature coefficients & the second set of 7 coefficients are low-temperature coefficients. Parsing can be defined as breaking the data block into smaller blocks by following a set of orders. so that it can be more easily interpreted, managed and sent to a computer NASA came up with the polynomials that can be used to evaluate thermodynamic properties such as C, Enthalpy, and Entropy using polynomials. They have also documented the co-efficient that are required to evaluate these polynomials for 1000+ species.
The formulas to calculate specific heat, enthalpy and entropy are given below-
SPECIFIC HEAT :
It is the amount of energy that must be added, in the form of heat, to one unit of mass of the substance in order to cause an increase of one unit in its temperature.
Specific heat is the thermodynamic property, which states the amount of heat required for a single unit of mass of a substance to be raised by one degree of temperature. Varying ranges of specific heat values are seen for substances depending on the extent to which they absorb heat. The term heat capacity can be misleading since heat q is the term given to the addition or removal of energy.
Cp = (a1 + a2 T + a3 T^2 + a4 T^3 + a5 T^4)*R
ENTHALPY :
Enthalpy is a concept used in science and engineering when heat and work need to be calculated. When a substance changes at constant pressure, enthalpy tells how much heat and work was added or removed from the substance. Enthalpy, a property of a thermodynamic system, is equal to the system's internal energy plus the product of its pressure and volume. In a system contained so as to prevent the mass transfer, for processes at constant pressure, the heat absorbed or released equals the change in enthalpy.
H = (a1 + a2 T /2 + a3 T^2 /3 + a4 T^3 /4 + a5 T^4 /5 + a6/T)*R*T
Entropy:
Entropy is an extensive property of a thermodynamic system. It is closely related to the number Ω of microscopic configurations (known as microstates) that are consistent
with the macroscopic quantities that characterize the system (such as its volume, pressure and temperature). The entropy of an object is a measure of the amount of energy which is unavailable to do work. Entropy is also a measure of the number of possible arrangements the atoms in a system can have. In this sense, entropy is a measure of uncertainty or randomness.
S = (a1 lnT + a2 T + a3 T^2 /2 + a4 T^3 /3 + a5 T^4 /4 + a7)*R
MOLECULAR WEIGHT :
The molecular mass is the mass of a given molecule: it is measured in daltons. Different molecules of the same compound may have different molecular masses because they contain different isotopes of an element.
Function Programs -
FUNCTION OF SPECIFIC HEAT CODE :
% specific heat calculation function
function [cp] = Specificheat(a1,a2,a3,a4,a5,a8,a9,a10,a11,a12,temp,R,global_medium_temp)
if temp <= global_medium_temp
cp = R*(a1 + a2*(temp) + a3*((temp).^2) + a4*((temp).^3) + a5*((temp).^4));
else
cp = R*(a8 + a9*(temp) + a10*((temp).^2) + a11*((temp).^3) + a12*((temp).^4));
endFUNCTION OF ENTHALPY CODE :
% enthalpy calculation function
function [H] =Enthalpy(a1,a2,a3,a4,a5,a6,a8,a9,a10,a11,a12,a13,temp,R,global_medium_temp)
if temp <= global_medium_temp
H = (R*(temp)).*(a1 + ((a2*temp)/2) + ((a3*(temp).^2)/3) + ((a4*(temp).^3)/4) + ((a5*(temp).^3)/4) + (a6./(temp)));
else
H = (R*(temp)).*(a8 + ((a9*temp)/2) + ((a10*(temp).^2)/3) + ((a11*(temp).^3)/4) + ((a12*(temp).^3)/4) + (a13./(temp)));
endFUNCTION OF ENTROPY CODE:
% entropy calculation function
function [S] = Entropy_calc(a1,a2,a3,a4,a5,a7,a8,a9,a10,a11,a12,a14,temp,R,global_medium_temp)
if temp <= global_medium_temp
S = R*(a1*(log(temp)) + a2*(temp) + a3*((temp).^(2/2)) + a4*(temp.^(3/3)) + a5*(temp.^(4/4)) + a7);
else
% else maximum coefficients are taken
S = R*(a8*(log(temp)) + a9*(temp) + a10*((temp).^(2/2)) + a11*((temp).^(3/3)) + a12*((temp).^(4/4)) + a14);
endFUNCTION FOR MOLECULAR WEIGHT CODE:
function M_WT = Molecular_weight(Species_SP)
% Standard - Atomic Name
atomic_name = ['H','C','O','N','A'];
% Hydrogen, Carbon, Oxygen, Nitrogen, Argon
fprintf('Molecular Weight of ');
fprintf(Species_SP);
fprintf(' : ');
% Standard Atomic Number
Atomic_Weight = [1.008 12.011 15.999 14.007 39.948];
% Starting Range of Molecule
M_WT = 0;
% Adding Atomic number to obtain particular value
for i = 1:length(Species_SP)
for j = 1:length(atomic_name)
if strcmpi(Species_SP(i),atomic_name(j))
M_WT= M_WT + Atomic_Weight(j);
position = j;
fprintf('%d',M_WT);
disp(' ');
end
end
n = str2double(Species_SP(i));
if n>1
M_WT = M_WT + Atomic_Weight(position)*(n-1);
end
end
fprintf('molecular weight %s :',Species_SP);
fprintf('%d',M_WT);
disp(' ');
endTHE MAIN PROGRAM -
clear all
clear all
clc
% step 1: Load file
f1 = fopen('THERMO.dat','r');
fgetl(f1);
%skip first line
% Obtaing Temperatures - by Reading Header Information
global_temp=fgetl(f1);
% Spliting
t=strsplit(global_temp);
% Converting string cells to numbers using str2num,
global_low_temp=str2double(t{2});
global_medium_temp=str2double(t{3});
global_high_temp=str2double(t{4});
% Temperature range & Gas Constant
% Temperature Range assumed
Temp=linspace(global_low_temp,global_high_temp,1000);
% Universal Gas Constant J/(mol-K)
R = 8.314;
for i=1:3
fgetl(f1);
end
for i = 1:53
line=fgetl(f1);
A=strsplit(line,' ');
Species_SP=(A{1});
line_1=fgetl(f1) ;
a=strfind(line_1,'E');
% Extracting higher and lower temperature coefficients
% Extracting High-temperature coefficients
a1=line_1(1:a(1)+3);
a1=str2double(a1);
a2=line_1(a(1)+4:a(2)+3);
a2=str2double(a2);
a3=line_1(a(2)+4:a(3)+3);
a3=str2double(a3);
a4=line_1(a(3)+4:a(4)+3);
a4=str2double(a4);
a5=line_1(a(4)+4:a(5)+3);
a5=str2double(a5);
line_2=fgetl(f1);
b=strfind(line_1,'E');
a6=line_2(1:b(1)+3);
a6=str2double(a6);
a7=line_2(a(1)+4:b(2)+3);
a7=str2double(a7);
% LOW-temperature coefficients
a8=line_2(b(2)+4:b(3)+3);
a8=str2double(a8);
a9=line_2(b(3)+4:b(4)+3);
a9=str2double(a9);
a10=line_2(b(4)+4:b(5)+3);
a10=str2double(a10);
line_3=fgetl(f1);
c=strfind(line_1,'E');
a11=line_3(1:c(1)+3);
a11=str2double(a11);
a12=line_3(c(1)+4:c(2)+3);
a12=str2double(a12);
a13=line_3(c(2)+4:c(3)+3);
a13=str2double(a13);
a14=line_3(c(3)+4:c(4)+3);
a14=str2double(a14);
% Specific heat calculation
Cp = Specificheat(a1,a2,a3,a4,a5,a8,a9,a10,a11,a12,Temp,R,global_medium_temp);
% Enthalpy calculation
H = Enthalpy(a1,a2,a3,a4,a5,a6,a8,a9,a10,a11,a12,a13,Temp,R,global_medium_temp);
% Entropy calculation
S = Entropy_calc(a1,a2,a3,a4,a5,a7,a8,a9,a10,a11,a12,a14,Temp,R,global_medium_temp);
% Calculating molecular weight of species
M_WT = Molecular_weight(Species_SP);
% Read & Write Molecular list
Molecular_list = fopen('Molecular weight.txt','a+');
fprintf( Molecular_list,'Molecular weight of %s is %d n',Species_SP,M_WT);
fclose(Molecular_list);
% plotting(CP,H,S,Temp,Species_SP)
% Saving as a folder
cur_dir=pwd;
mkdir(Species_SP);
cd(Species_SP);
% Plot 1, Specific heat (CP) vs Temperature Range (Temp)
figure(1)
plot(Temp,Cp,'linewidth',2,'color','r');
xlabel('Temperature [K]');
ylabel('Specific Heat [kJ]');
title(sprintf('Varying Specific Heat with Temperature Range of %s',Species_SP));
grid on
saveas(figure(1),'Specific Heat.jpg');
% Plot 2, Enthalpy (H) vs Temperature Range (Temp)
figure(2)
plot(Temp,H,'linewidth',2,'color','g');
xlabel('Temperature [K]');
ylabel('Enthalpy in [kJ/(mol)]');
title(sprintf('Varying Enthalpy with Temperature Range of %s',Species_SP));
grid on
saveas(figure(2),'Ethalpy.jpg');
% Plot 3, Entropy (S) vs Temperature Range (Temp)
figure(3)
plot(Temp,S,'linewidth',2,'color','b');
xlabel('Temperature [K]');
ylabel('Entropy in [kJ/(mol-K)]');
title(sprintf('Varying Entropy with Temperature Range of %s',Species_SP));
grid on
saveas(figure(3),'Entropy.jpg');
cd(cur_dir);
% Molecular weight
M = Molecular_weight(Species_SP);
f2 = fopen('molecular mass','w');
fprintf('%s n %s n %f n','species','Molecular mass',M)
fclose(f2);
endCODE EXPLANATION-
1. 'clear all' - clear all the variables in the workspace, 'close all' - closes all the plots and Figures and 'clc' - clear the command window.
2. Making a directory using 'mkdir' with location address.
3. After that, reading the text file 'THERMO.dat' file and opening it by using ‘fopen’ and 'r' as reading command.
4. For reading the line in a text file using ‘fgetl’.
5. To obtain the temperature range from the Thermo.dat file, using ‘strsplit’ to split the temperature into a particular string cell.
6. Used ‘str2double’ for converting string cell to numbers.
7. Divided the temperature into a string and then converted to a number so that temperature values can be used separately to solve the NASA
thermodynamics.
8. Using ‘for & end’ loops to skip the next three comment lines in the text file.
9. Reading all the 53 species to plot specific heat(cp), enthalpy(H), and entropy(S), with temperature values, and to find the molecular weights.
10. Used ‘for’ loop to read 53 species.
11. Used ‘strfind' to find string.
12. Using 'linspace' specified the temperature range from low to high.
13. Adding universal gas constant value ‘R’.R VALUE IS 0.287KJ/KG K
14. Extracted higher and lower temperature coefficient by reading line by line from the species.
15. Calling specific heat, enthalpy, entropy function codes.
16. Calculating Molecular weight using functions.
OUTPUT-
Molecular Weight of O : 1.599900e+01
molecular weight O :1.599900e+01
Molecular Weight of O : 1.599900e+01
molecular weight O :1.599900e+01
species n Molecular mass n 15.999000 nMolecular Weight of O2 : 1.599900e+01
molecular weight O2 :3.199800e+01
Molecular Weight of O2 : 1.599900e+01
molecular weight O2 :3.199800e+01
species n Molecular mass n 31.998000 nMolecular Weight of H : 1.008000e+00
molecular weight H :1.008000e+00
Molecular Weight of H : 1.008000e+00
molecular weight H :1.008000e+00
species n Molecular mass n 1.008000 nMolecular Weight of H2 : 1.008000e+00
molecular weight H2 :2.016000e+00
Molecular Weight of H2 : 1.008000e+00
molecular weight H2 :2.016000e+00
species n Molecular mass n 2.016000 nMolecular Weight of OH : 1.599900e+01
1.700700e+01
molecular weight OH :1.700700e+01
Molecular Weight of OH : 1.599900e+01
1.700700e+01
molecular weight OH :1.700700e+01
species n Molecular mass n 17.007000 nMolecular Weight of H2O : 1.008000e+00
1.801500e+01
molecular weight H2O :1.801500e+01
Molecular Weight of H2O : 1.008000e+00
1.801500e+01
molecular weight H2O :1.801500e+01
species n Molecular mass n 18.015000 nMolecular Weight of HO2 : 1.008000e+00
1.700700e+01
molecular weight HO2 :3.300600e+01
Molecular Weight of HO2 : 1.008000e+00
1.700700e+01
molecular weight HO2 :3.300600e+01
species n Molecular mass n 33.006000 nMolecular Weight of H2O2 : 1.008000e+00
1.801500e+01
molecular weight H2O2 :3.401400e+01
Molecular Weight of H2O2 : 1.008000e+00
1.801500e+01
molecular weight H2O2 :3.401400e+01
species n Molecular mass n 34.014000 nMolecular Weight of C : 1.201100e+01
molecular weight C :1.201100e+01
Molecular Weight of C : 1.201100e+01
molecular weight C :1.201100e+01
species n Molecular mass n 12.011000 nMolecular Weight of CH : 1.201100e+01
1.301900e+01
molecular weight CH :1.301900e+01
Molecular Weight of CH : 1.201100e+01
1.301900e+01
molecular weight CH :1.301900e+01
species n Molecular mass n 13.019000 nMolecular Weight of CH2 : 1.201100e+01
1.301900e+01
molecular weight CH2 :1.402700e+01
Molecular Weight of CH2 : 1.201100e+01
1.301900e+01
molecular weight CH2 :1.402700e+01
species n Molecular mass n 14.027000 nMolecular Weight of CH2(S) : 1.201100e+01
1.301900e+01
molecular weight CH2(S) :1.402700e+01
Molecular Weight of CH2(S) : 1.201100e+01
1.301900e+01
molecular weight CH2(S) :1.402700e+01
species n Molecular mass n 14.027000 nMolecular Weight of CH3 : 1.201100e+01
1.301900e+01
molecular weight CH3 :1.503500e+01
Molecular Weight of CH3 : 1.201100e+01
1.301900e+01
molecular weight CH3 :1.503500e+01
species n Molecular mass n 15.035000 nMolecular Weight of CH4 : 1.201100e+01
1.301900e+01
molecular weight CH4 :1.604300e+01
Molecular Weight of CH4 : 1.201100e+01
1.301900e+01
molecular weight CH4 :1.604300e+01
species n Molecular mass n 16.043000 nMolecular Weight of CO : 1.201100e+01
2.801000e+01
molecular weight CO :2.801000e+01
Molecular Weight of CO : 1.201100e+01
2.801000e+01
molecular weight CO :2.801000e+01
species n Molecular mass n 28.010000 nMolecular Weight of CO2 : 1.201100e+01
2.801000e+01
molecular weight CO2 :4.400900e+01
Molecular Weight of CO2 : 1.201100e+01
2.801000e+01
molecular weight CO2 :4.400900e+01
species n Molecular mass n 44.009000 nMolecular Weight of HCO : 1.008000e+00
1.301900e+01
2.901800e+01
molecular weight HCO :2.901800e+01
Molecular Weight of HCO : 1.008000e+00
1.301900e+01
2.901800e+01
molecular weight HCO :2.901800e+01
species n Molecular mass n 29.018000 nMolecular Weight of CH2O : 1.201100e+01
1.301900e+01
3.002600e+01
molecular weight CH2O :3.002600e+01
Molecular Weight of CH2O : 1.201100e+01
1.301900e+01
3.002600e+01
molecular weight CH2O :3.002600e+01
species n Molecular mass n 30.026000 nMolecular Weight of CH2OH : 1.201100e+01
1.301900e+01
3.002600e+01
3.103400e+01
molecular weight CH2OH :3.103400e+01
Molecular Weight of CH2OH : 1.201100e+01
1.301900e+01
3.002600e+01
3.103400e+01
molecular weight CH2OH :3.103400e+01
species n Molecular mass n 31.034000 nMolecular Weight of CH3O : 1.201100e+01
1.301900e+01
3.103400e+01
molecular weight CH3O :3.103400e+01
Molecular Weight of CH3O : 1.201100e+01
1.301900e+01
3.103400e+01
molecular weight CH3O :3.103400e+01
species n Molecular mass n 31.034000 nMolecular Weight of CH3OH : 1.201100e+01
1.301900e+01
3.103400e+01
3.204200e+01
molecular weight CH3OH :3.204200e+01
Molecular Weight of CH3OH : 1.201100e+01
1.301900e+01
3.103400e+01
3.204200e+01
molecular weight CH3OH :3.204200e+01
species n Molecular mass n 32.042000 nMolecular Weight of C2H : 1.201100e+01
2.503000e+01
molecular weight C2H :2.503000e+01
Molecular Weight of C2H : 1.201100e+01
2.503000e+01
molecular weight C2H :2.503000e+01
species n Molecular mass n 25.030000 nMolecular Weight of C2H2 : 1.201100e+01
2.503000e+01
molecular weight C2H2 :2.603800e+01
Molecular Weight of C2H2 : 1.201100e+01
2.503000e+01
molecular weight C2H2 :2.603800e+01
species n Molecular mass n 26.038000 nMolecular Weight of C2H3 : 1.201100e+01
2.503000e+01
molecular weight C2H3 :2.704600e+01
Molecular Weight of C2H3 : 1.201100e+01
2.503000e+01
molecular weight C2H3 :2.704600e+01
species n Molecular mass n 27.046000 nMolecular Weight of C2H4 : 1.201100e+01
2.503000e+01
molecular weight C2H4 :2.805400e+01
Molecular Weight of C2H4 : 1.201100e+01
2.503000e+01
molecular weight C2H4 :2.805400e+01
species n Molecular mass n 28.054000 nMolecular Weight of C2H5 : 1.201100e+01
2.503000e+01
molecular weight C2H5 :2.906200e+01
Molecular Weight of C2H5 : 1.201100e+01
2.503000e+01
molecular weight C2H5 :2.906200e+01
species n Molecular mass n 29.062000 nMolecular Weight of C2H6 : 1.201100e+01
2.503000e+01
molecular weight C2H6 :3.007000e+01
Molecular Weight of C2H6 : 1.201100e+01
2.503000e+01
molecular weight C2H6 :3.007000e+01
species n Molecular mass n 30.070000 nMolecular Weight of CH2CO : 1.201100e+01
1.301900e+01
2.603800e+01
4.203700e+01
molecular weight CH2CO :4.203700e+01
Molecular Weight of CH2CO : 1.201100e+01
1.301900e+01
2.603800e+01
4.203700e+01
molecular weight CH2CO :4.203700e+01
species n Molecular mass n 42.037000 nMolecular Weight of HCCO : 1.008000e+00
1.301900e+01
2.503000e+01
4.102900e+01
molecular weight HCCO :4.102900e+01
Molecular Weight of HCCO : 1.008000e+00
1.301900e+01
2.503000e+01
4.102900e+01
molecular weight HCCO :4.102900e+01
species n Molecular mass n 41.029000 nMolecular Weight of HCCOH : 1.008000e+00
1.301900e+01
2.503000e+01
4.102900e+01
4.203700e+01
molecular weight HCCOH :4.203700e+01
Molecular Weight of HCCOH : 1.008000e+00
1.301900e+01
2.503000e+01
4.102900e+01
4.203700e+01
molecular weight HCCOH :4.203700e+01
species n Molecular mass n 42.037000 nMolecular Weight of H2CN : 1.008000e+00
1.402700e+01
2.803400e+01
molecular weight H2CN :2.803400e+01
Molecular Weight of H2CN : 1.008000e+00
1.402700e+01
2.803400e+01
molecular weight H2CN :2.803400e+01
species n Molecular mass n 28.034000 nMolecular Weight of HCN : 1.008000e+00
1.301900e+01
2.702600e+01
molecular weight HCN :2.702600e+01
Molecular Weight of HCN : 1.008000e+00
1.301900e+01
2.702600e+01
molecular weight HCN :2.702600e+01
species n Molecular mass n 27.026000 nMolecular Weight of HNO : 1.008000e+00
1.501500e+01
3.101400e+01
molecular weight HNO :3.101400e+01
Molecular Weight of HNO : 1.008000e+00
1.501500e+01
3.101400e+01
molecular weight HNO :3.101400e+01
species n Molecular mass n 31.014000 nMolecular Weight of N : 1.400700e+01
molecular weight N :1.400700e+01
Molecular Weight of N : 1.400700e+01
molecular weight N :1.400700e+01
species n Molecular mass n 14.007000 nMolecular Weight of NNH : 1.400700e+01
2.801400e+01
2.902200e+01
molecular weight NNH :2.902200e+01
Molecular Weight of NNH : 1.400700e+01
2.801400e+01
2.902200e+01
molecular weight NNH :2.902200e+01
species n Molecular mass n 29.022000 nMolecular Weight of N2O : 1.400700e+01
4.401300e+01
molecular weight N2O :4.401300e+01
Molecular Weight of N2O : 1.400700e+01
4.401300e+01
molecular weight N2O :4.401300e+01
species n Molecular mass n 44.013000 nMolecular Weight of NH : 1.400700e+01
1.501500e+01
molecular weight NH :1.501500e+01
Molecular Weight of NH : 1.400700e+01
1.501500e+01
molecular weight NH :1.501500e+01
species n Molecular mass n 15.015000 nMolecular Weight of NH2 : 1.400700e+01
1.501500e+01
molecular weight NH2 :1.602300e+01
Molecular Weight of NH2 : 1.400700e+01
1.501500e+01
molecular weight NH2 :1.602300e+01
species n Molecular mass n 16.023000 nMolecular Weight of NH3 : 1.400700e+01
1.501500e+01
molecular weight NH3 :1.703100e+01
Molecular Weight of NH3 : 1.400700e+01
1.501500e+01
molecular weight NH3 :1.703100e+01
species n Molecular mass n 17.031000 nMolecular Weight of NO : 1.400700e+01
3.000600e+01
molecular weight NO :3.000600e+01
Molecular Weight of NO : 1.400700e+01
3.000600e+01
molecular weight NO :3.000600e+01
species n Molecular mass n 30.006000 nMolecular Weight of NO2 : 1.400700e+01
3.000600e+01
molecular weight NO2 :4.600500e+01
Molecular Weight of NO2 : 1.400700e+01
3.000600e+01
molecular weight NO2 :4.600500e+01
species n Molecular mass n 46.005000 nMolecular Weight of HCNO : 1.008000e+00
1.301900e+01
2.702600e+01
4.302500e+01
molecular weight HCNO :4.302500e+01
Molecular Weight of HCNO : 1.008000e+00
1.301900e+01
2.702600e+01
4.302500e+01
molecular weight HCNO :4.302500e+01
species n Molecular mass n 43.025000 nMolecular Weight of HOCN : 1.008000e+00
1.700700e+01
2.901800e+01
4.302500e+01
molecular weight HOCN :4.302500e+01
Molecular Weight of HOCN : 1.008000e+00
1.700700e+01
2.901800e+01
4.302500e+01
molecular weight HOCN :4.302500e+01
species n Molecular mass n 43.025000 nMolecular Weight of HNCO : 1.008000e+00
1.501500e+01
2.702600e+01
4.302500e+01
molecular weight HNCO :4.302500e+01
Molecular Weight of HNCO : 1.008000e+00
1.501500e+01
2.702600e+01
4.302500e+01
molecular weight HNCO :4.302500e+01
species n Molecular mass n 43.025000 nMolecular Weight of NCO : 1.400700e+01
2.601800e+01
4.201700e+01
molecular weight NCO :4.201700e+01
Molecular Weight of NCO : 1.400700e+01
2.601800e+01
4.201700e+01
molecular weight NCO :4.201700e+01
species n Molecular mass n 42.017000 nMolecular Weight of CN : 1.201100e+01
2.601800e+01
molecular weight CN :2.601800e+01
Molecular Weight of CN : 1.201100e+01
2.601800e+01
molecular weight CN :2.601800e+01
species n Molecular mass n 26.018000 nMolecular Weight of HCNN : 1.008000e+00
1.301900e+01
2.702600e+01
4.103300e+01
molecular weight HCNN :4.103300e+01
Molecular Weight of HCNN : 1.008000e+00
1.301900e+01
2.702600e+01
4.103300e+01
molecular weight HCNN :4.103300e+01
species n Molecular mass n 41.033000 nMolecular Weight of N2 : 1.400700e+01
molecular weight N2 :2.801400e+01
Molecular Weight of N2 : 1.400700e+01
molecular weight N2 :2.801400e+01
species n Molecular mass n 28.014000 nMolecular Weight of AR : 3.994800e+01
molecular weight AR :3.994800e+01
Molecular Weight of AR : 3.994800e+01
molecular weight AR :3.994800e+01
species n Molecular mass n 39.948000 nMolecular Weight of C3H8 : 1.201100e+01
3.704100e+01
molecular weight C3H8 :4.409700e+01
Molecular Weight of C3H8 : 1.201100e+01
3.704100e+01
molecular weight C3H8 :4.409700e+01
species n Molecular mass n 44.097000 nMolecular Weight of C3H7 : 1.201100e+01
3.704100e+01
molecular weight C3H7 :4.308900e+01
Molecular Weight of C3H7 : 1.201100e+01
3.704100e+01
molecular weight C3H7 :4.308900e+01
species n Molecular mass n 43.089000 nMolecular Weight of CH3CHO : 1.201100e+01
1.301900e+01
2.704600e+01
2.805400e+01
4.405300e+01
molecular weight CH3CHO :4.405300e+01
Molecular Weight of CH3CHO : 1.201100e+01
1.301900e+01
2.704600e+01
2.805400e+01
4.405300e+01
molecular weight CH3CHO :4.405300e+01
species n Molecular mass n 44.053000 nMolecular Weight of CH2CHO : 1.201100e+01
1.301900e+01
2.603800e+01
2.704600e+01
4.304500e+01
molecular weight CH2CHO :4.304500e+01
Molecular Weight of CH2CHO : 1.201100e+01
1.301900e+01
2.603800e+01
2.704600e+01
4.304500e+01
molecular weight CH2CHO :4.304500e+01
species n Molecular mass n 43.045000 n>>
OUTPUT GRAPHS FOR O2, CO2, N2
SNIP OF THE FILES SAVED
Conclusion-
Using MATLAB code, the necessary coefficients required to evaluate various thermodynamic parameters for different species are extracted from NASA
thermodynamic data file. The specific heat, enthalpy, and entropy for all the species are calculated and plotted against the local temperature range. The
molecular weight is calculated for all the species.
References-
https://www.quora.com/What-does-it-means-to-parse-a-file
https://www.mathworks.com/help/matlab/ref/mkdir.html
google drive link- all files folder link