Project 1 - Parsing NASA thermodynamic data

AIM: To write a program in matlab to parse the NASA thermodynamic data file and then calculate t various gas species.

OBJECTIVE:

To write a function that extracts the 14 co-efficient and calculates the enthalpy ,entropy and s data file.

Calculating the molecular weight of each species and displaying it in command window. Plotting cp, enthalpy and entropy for local temperature range specific for each species.

Save the plot as images with appropriate names and dump them in separate folders for each

INTRODUCTION:

PARSING DATA:

Parsing means to make something understandable. For programming this means to convert i form into another form that's easier to work with. This is done by partially analysing the data, u structure (by making some assumptions based on what you're expecting to see), and then ex in the code.

Parsing a file means reading in a data stream of some sort and building an in memory model data. This aims at facilitating perfoming some kind of transformation on the data.

EQUATIONS:

The equations that we are going to use in the Functions to calculate cp, H,S are:

cp                                   2                 3                          4

= a₁ + a₂T + a₃T  + a₄T  + a₅T

R

H                   a2T            T ²           T ³           T ⁴       a6

= a₁ +           + a₃                           + a₄                           + a₅                           + RT                   2               3              4              5         T

S                                         T ²           T ³           T ⁴

= a₁ ln T + a₂T + a₃                                              + a₄        + a₅        + a7 R            2                            3                                     4

Where, Cp=specific heat

R =molar gas constant

T=temperature S=Entropy H=Enthalpy a=coefficients

BODY OF THE CONTENT:

Code for calling function Specific_heat( ):

function Cp = Specific_heat(a1,a2,a3,a4,a5,a8,a9,a10,a11,a12,T,R,global_midtemp)

%calculating specific_heat for higher and lower temperature respectively if(T>= global_midtemp)

Cp=R*(a1+(a2*T)+(a3*(T).^2)+(a4*(T).^3)+(a5*(T).^4));

else

Cp=R*(a8+(a9*T)+(a10*(T).^2)+(a11*(T).^3)+(a12*(T).^4));

end end

Code for calling function Enthalpy( ):

function H = Enthalpy(a1,a2,a3,a4,a5,a6,a8,a9,a10,a11,a12,a13,T,R,global_midtemp)

%calculating enthalpy for higher and lower temperature respectively if(T>=global_midtemp)

H=R*((a1*T)+((a2*(T).^2)/2)+((a3*(T).^3)/3)+((a4*(T).^4)/4)+((a5*(T).^5)/5)+a6);

else

H=R*((a8*T)+((a9*(T).^2)/2)+((a10*(T).^3)/3)+((a11*(T).^4)/4)+((a12*(T).^5)/5)+a13);

end end

Code for calling function Entropy( ):

function S=Entropy(a1,a2,a3,a4,a5,a7,a8,a9,a10,a11,a12,a14,T,R,global_midtemp)

%calculating entropy for higher and lower temperature respectively if(T>= global_midtemp)

S=R*(a1.*log(T)+(a2*T)+((a3*(T).^2)/2)+((a4*(T).^3)/3)+((a5*(T).^4)/5)+a7);

else

S=R*(a8.*log(T)+(a9*T)+((a10*(T).^2)/2)+((a11*(T).^3)/3)+((a12*(T).^4)/4));

end end

function Molecular_weight=wt(species_name)

Elements=['H' 'C' 'O' 'N' 'A' 'S'];

Atomic_weight =[1 12 16 14 40 32]; species=upper(species_name)

Molecular_weight=0;

for i=1:length(species)

for    j=1:length(Elements)

if strcmp(species(i),Elements(j))

Molecular_weight=Molecular_weight+Atomic_weight(j); k=j;

else

n=str2double(species(i)); if (n>1)

Molecular_weight=Molecular_weight+Atomic_weight(k)*(n-1); break

end

end end

end

end

MainCode-