Project 1 - Parsing NASA thermodynamic data

• First we write the main code and then we write the calling functions for the properties of individual species from THERMO.dat file such as Specific Heat (Cp), Enthalpy (H), Entropy (S) and Molecular Weight.
• In the main code we extract global temperatures, each species and their polynomial coeffiecients and their local temperatures from the NASA file.
• Then we pass them to calling functions.
• Here we are using two functions one for Cp,H and S collectively and the other for Molecular Weight.
• Cp, H and S are calculated using the given formulae in the function code.
• Molecular Weight function is calculated using for loops and if statements and the atomic weight of the main elements in the species such as Hydrogen (H), Carbon (C), Oxygen (O), Nitrogen (N) and Argon (A).
• We print the Molecular weight in the command window for each species in the main code.
• Then we plot the required figures for each species, save them in individual folders and place them in a parent folder called ' NASA_Species '.

MAIN CODE:

% MATLAB PROJECT 1 - PARSING NASA THERMODYNAMIC DATA

% Main Program 

% clearing all previous stored values and clearing command window
clear all
close all
clc

% Creating a new Folder to store the output
mkdir Nasa_Species 
R=8.314; % Molar Gas Constant

% Opening and reading THERMO.dat file
f1= fopen('THERMO.dat','r'); % Open and reading THERMO.dat file

% Skipping the first line
fgetl(f1); % Line 1

line2=fgetl(f1); % Line 2 assigned to a variable to extract further data
line2_split=strsplit(line2,' '); % using strsplit to split each value by separating from the spaces 
                                 % this variable is now a cell array

% Accessing each value from the cell array and storing into separate variables  
% Global Temperatures
global_low_temp=str2num(line2_split{2}); 
global_mid_temp=str2num(line2_split{3});
global_high_temp=str2num(line2_split{4});

% Temperature range
T=linspace(global_low_temp,global_high_temp,1000);

% Skipping the next 3 lines
fgetl(f1); % Line 3
fgetl(f1); % Line 4
fgetl(f1); % Line 5

%% Reading Species, Calcuating Cp,H,S,MW and Plotting

for i=1:53 % for 53 species

    tline=fgetl(f1); 

    A=strsplit(tline,' ');  % Splitting the line to cell array
    name_of_species=(A{1}); % Getting the name of the species
   
    lowtemp=str2num(A{end-3}); % Getting the low tempeature
    midtemp=str2num(A{end-2}); % Getting the low tempeature
    hightemp=str2num(A{end-1}); %Getting the low tempeature
 
    % Finding 7 higher and 7 lower coefficients  

    tline1=fgetl(f1); 
    a=findstr(tline1,'E'); % Finding the indices/indexes of character E

    a1=tline1(1:a(1)+3); % First value
    a1=str2num(a1); % Changing the string value to numerical value
    a2=tline1(a(1)+4:a(2)+3);
    a2=str2num(a2);
    a3=tline1(a(2)+4:a(3)+3);
    a3=str2num(a3);
    a4=tline1(a(3)+4:a(4)+3);
    a4=str2num(a4);
    a5=tline1(a(4)+4:a(5)+3);
    a5=str2num(a5);

    tline2=fgetl(f1); 
    a=findstr(tline2,'E'); 

    a6=tline2(1:a(1)+3); 
    a6=str2num(a6); 
    a7=tline2(a(1)+4:a(2)+3);
    a7=str2num(a7);
    a8=tline2(a(2)+4:a(3)+3);
    a8=str2num(a8);
    a9=tline2(a(3)+4:a(4)+3);
    a9=str2num(a9);
    a10=tline2(a(4)+4:a(5)+3);
    a10=str2num(a10);

    tline3=fgetl(f1); 
    a=findstr(tline3,'E'); 

    a11=tline3(1:a(1)+3); 
    a11=str2num(a11); 
    a12=tline3(a(1)+4:a(2)+3);
    a12=str2num(a12);
    a13=tline3(a(2)+4:a(3)+3);
    a13=str2num(a13);
    a14=tline3(a(3)+4:a(4)+3);
    a14=str2num(a14);
   
    % Function Calling  
    % Calculating Specific Heat(Cp), Enthalpy(H), Entropy(S)
    [Cp,H,S] = Properties(a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13,a14,R,T,global_mid_temp);
    
    % Calculating Molecular Weight and printing
    MW = Molecular_Weight(name_of_species);
    fprintf('Molecular Weight of %s= %f\n',name_of_species,MW);

    % Changing the directory to Sub Folder Nasa_Species
    cd Nasa_Species 
    mkdir(name_of_species); % Creating subfolder for each species with names
    cd(name_of_species); % Changing current folder to species folder


    % Plotting CP,H,T

    % Plotting Cp vs T
    figure(1)
    plot(T,Cp,'linewidth',3,'Color','b');
    xlabel('Temperature (T)');
    ylabel('Specific heat (Cp)');
    title(name_of_species, 'Specific Heat vs Temperature')
    grid on
    saveas(figure(1),'Specific Heat.jpg') % Saving the plot as JPG file

    % Plotting H vs T
    figure(2)
    plot(T,H,'linewidth',3,'Color','r');
    xlabel('Temperature (T)');
    ylabel('Enthalpy (H)');
    title(name_of_species, 'Enthalpy vs Temperature')
    grid on
    saveas(figure(2),'Enthalpy.jpg')
    
    % Plotting S vs T
    figure(3)
    plot(T,S,'linewidth',3,'Color','y');
    xlabel('Temperature (T)');
    ylabel('Entropy (S)');
    title(name_of_species, 'Entropy vs Temperature')
    grid on
    saveas(figure(3),'Entropy.jpg')
   
    % Changing back to Parent Folder
    cd('C:\Users\Rinil\Desktop\SKILL LYNC\Matlab Challenges\Project1-File_Parsing')

end

Function for Specific Heat, Enthalpy and Entropy:

% Specific Heat, Enthalpy and Entropy function
% for lower or higher temperarures

function [Cp,H,S] = Properties(a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13,a14,R,T,global_mid_temp)

if(T>=global_mid_temp)
    Cp=R*(a1+(a2*T)+(a3*(T).^2)+(a4*(T).^3)+(a5*(T).^4));
    H=R*((a1*T)+((a2*(T).^2)/2)+((a3*(T).^3)/3)+((a4*(T).^4)/4)+((a5*(T).^5)/5)+a6);
    S=R*(a1.*log(T)+(a2*T)+((a3*(T).^2)/2)+((a4*(T).^3)/3)+((a5*(T).^4)/5)+a7);
else
    Cp=R*(a8+(a9*T)+(a10*(T).^2)+(a11*(T).^3)+(a12*(T).^4));
    H=R*((a8*T)+((a9*(T).^2)/2)+((a10*(T).^3)/3)+((a11*(T).^4)/4)+((a12*(T).^5)/5)+a13);
    S=R*(a8.*log(T)+(a9*T)+((a10*(T).^2)/2)+((a11*(T).^3)/3)+((a12*(T).^4)/4)+a14);
end
end

% Molecular Weight Function

function W = Molecular_Weight(name_of_species)

    % Defining Atomic name of the main elements in the species list
    atomic_name = ['H','C','O','N','Ar'];
    
    % Defining Atomic masses of the main elements
    atomic_weight = [1.008 12.011 15.999 14.007 39.948];

    % Molecular weight starting range
    W = 0;  
    species = name_of_species;

    for i=1:length(species)   % Loop runs to the length of the species
        for j=1:length(atomic_name)  % Loop runs to the length of the atomic names
            if strcmp(species(i), atomic_name(j))  % Function to compare 2 strings
                W = W + atomic_weight(j);
                position = j;
            end
        end

        % Finding if there is more elements in the species 
        n = str2num(species); 
        
        if n>1
            W = W + atomic_weight(position)*(n-1);  % Incrementing the atomic weight of the extra element
        end
    end
   
end