Project 1 - Parsing NASA thermodynamic data

clear all
close all
clc

f1 = fopen('THERMO (1).dat','r');
% skipping next 5 line
for i= 1:5
    fgetl(f1);
end

%storting starts

for i = 1:53

    %creating an array species
    tline = fgetl(f1)
    A = strsplit(tline,' ');

    %species
   species = A{1}

    %local temperature 
    tlow = str2num(A{end-3});
    thigh = str2num(A{end-2});
    tmid = str2num(A{end-1});

    T = linspace(tlow,thigh,500)
    %finding coefficent
    %line 1
    tline1 =fgetl(f1)
    pos1 = findstr(tline1,'E');
    
    %high coefficent
    a1 = str2num(tline1(1:pos1(1)+3));
    a2 = str2num(tline1(pos1(1)+4:pos1(2)+3));
    a3 = str2num(tline1(pos1(2)+4:pos1(3)+3));
    a4 = str2num(tline1(pos1(3)+4:pos1(4)+3));
    a5 = str2num(tline1(pos1(4)+4:pos1(5)+3));

    %line2
    tline2 = fgetl(f1)
    pos2 = findstr(tline2,'E');
    a6 = str2num(tline2(1:pos2(1)+3));
    a7 = str2num(tline2(pos2(1)+4:pos2(2)+3));

    %lower coefficent
    a8 = str2num(tline2(pos2(2)+4:pos2(3)+3));
    a9 = str2num(tline2(pos2(3)+4:pos2(4)+3));
    a10 = str2num(tline2(pos2(4)+4:pos2(5)+3));

    % line3
    tline3 = fgetl(f1)
    pos3 = findstr(tline3,'E');
    a11 = str2num(tline3(1:pos3(1)+3));
    a12 = str2num(tline3(pos3(1)+4:pos3(2)+3));
    a13 = str2num(tline3(pos3(2)+4:pos3(3)+3));
    a14 = str2num(tline3(pos3(3)+4:pos3(4)+3));
  R = 8.314
    % specific heat 
       cp = specific_heat(a1,a2,a3,a4,a5,a8,a9,a10,a11,a12,tmid,T,R);
    % entopy
       S = entropy_S(a1,a2,a3,a4,a5,a7,a8,a9,a10,a11,a12,a14,tmid,T,R);
    % enthalpy
       H = enthalpy_H(a1,a2,a3,a4,a5,a6,a8,a9,a10,a11,a12,a13,tmid,T,R);

   mkdir(species)
   cd(species)

       figure(1)
       plot(T,cp)
       xlabel('temperature')
       ylabel('specific_heat')
       title('specific_heat vs temperature of',sprintf('%s',species));
       saveas(figure(1),'cp vs T.jpg')

        figure(2)
       plot(T,S)
       xlabel('temperature')
       ylabel('entropy')
       title('entropy vs temperature of',sprintf('%s',species));
       saveas(figure(2),'S vs T.jpg')

         figure(3)
       plot(T,H)
       xlabel('temperature')
       ylabel('enthalpy')
       title('enthalpy vs temperature of',sprintf('%s',species));
       saveas(figure(3),'H vs T.jpg')
      % comming back to program directory
      
       cd ..
       
       % calculat molicular weight of species 
       M = molecular_weight(species);
end

Explaination :

Aim :

To write a function that extracts the 14 co-efficients and to calculate the enthalpy, entropy and specific heats for all the species in the data file

Equation :

% specific heat

cp = specific_heat(a1,a2,a3,a4,a5,a8,a9,a10,a11,a12,tmid,T,R);
% entopy

S = entropy_S(a1,a2,a3,a4,a5,a7,a8,a9,a10,a11,a12,a14,tmid,T,R);
% enthalpy

H = enthalpy_H(a1,a2,a3,a4,a5,a6,a8,a9,a10,a11,a12,a13,tmid,T,R);

We create the function of enthalpy and entropy and specific heat of the co-efficient.
And to combine the equation on main program.

Objective of the project :

The main objective of the project is to create the temperature to given species
It was use to see the individual temperature of entropy enthalpy and the specific heat vs temperature in each species
For example we can take CO2 , O2 and N2 species and to show the plot.

Body of contant :

In case of programming project :

It uses the Matlab programming language
First we read the file and let we assume the f1
We are going to ommite the first 5 line on termodynamic data table
And it creating array species for using for loop command
We found local temp ,tlow = str2num(A{end-3}); thigh = str2num(A{end-2}); tmid = str2num(A{end-1});
Let we assume tha temp(T) at linspace(tlow,thigh, 500)
Now we finding coefficient of temperature at line1

In case of simulation project :

Now we found local temp on thermodynamic data at a1 to a14

a1 = str2num(tline1(1:pos1(1)+3));

a2 = str2num(tline1(pos1(1)+4:pos1(2)+3));

a3 = str2num(tline1(pos1(2)+4:pos1(3)+3));

a4 = str2num(tline1(pos1(3)+4:pos1(4)+3));

a5 = str2num(tline1(pos1(4)+4:pos1(5)+3));

Meanwhile from a6 to a14 has the same procedure,
R = 8.314

% specific heat
cp = specific_heat(a1,a2,a3,a4,a5,a8,a9,a10,a11,a12,tmid,T,R);

% entopy
S = entropy_S(a1,a2,a3,a4,a5,a7,a8,a9,a10,a11,a12,a14,tmid,T,R);

% enthalpy
H = enthalpy_H(a1,a2,a3,a4,a5,a6,a8,a9,a10,a11,a12,a13,tmid,T,R);

It will be consider the function of entropy and enthalpy and specific heat and finally we calculate the molecular weight.

In case of design project :

We create the Mkdir and cd species of the temp
Let we found figure1 and figure2 and figure3
It will be plotting over all 3 species like cp. and H and S
First plot(Cp,T) second plot (H,T) third plot(S,T) finally we create xlabel and ylabel and title
Finally can create saveas(figure(1),'cp vs T.jpg') and saveas(figure(2),'H vs T.jpg') and saveas(figure(3),'S vs T.jpg'
And attach the function of molicular weight of temp its M = molecular_weight(species);

Conclusion :

Thus the Matlab function to calculate the enthalpy, entropy and specific heats for all the species in the data file is plotted succesfully.