FILE PARSING

-

 NASA THERMODYNAMIC DATA:

Thermodynamics is a branch of physics that deals with heat, work, and temperature, and their relation to energy, radiation, and physical properties of matter.  Thermodynamics applies to a wide variety of topics in science and engineering, especially physical chemistry, biochemistry, chemical engineering and mechanical engineering.

Some of thermodynamic data were given by NASA,that has been shown here…

In this senerio, we just extract the values,which is already written in that particular file and that values used to find some other datas like specific heat,entropy and enthalpy with the help of matlab.

Specific heat:

 The specific heat capacity of a substance is the heat capacity of a sample of the substance divided by the mass of the sample.

Entropy:

Entropy is a scientific concept as well as a measurable physical property that is most commonly associated with a state of disorder, randomness, or uncertainty.

Enthalpy:

Enthalpy, a property of a thermodynamic system, is the sum of the system's internal energy and the product of its pressure and volume. It is a state function used in many measurements in chemical, biological, and physical systems at a constant pressure.

Formulla to calculate specific heat(cp),entropy(h),enthalpy(s):

Write a function, to extarct the values form data file and to calculate specific heat(cp),entropy(h),enthalpy(s):

function [cp h s]=co_eff(a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13,a14,T)
R=8.314;
for i=1:length(T)
    if T(i)<300;
        cp = (a1+a2*T+a3*T.^2+a4*T.^3+a5*T.^4)*R
        h=(a1+(a2*T)./2+(a3*T.^2)./3+a4*T.^3+a5*T.^4+a6./T).*(R.*T)
        s=(a1*log(T)+a2*T+(a3*T.^2)./2+(a4*T.^3)./3+(a5*T.^4)./4+a7)*R

    else
          cp=(a8+a9*T+a10*T.^2+a11*T.^3+a12*T.^4)*R
        h=(a8+(a9*T)./2+(a10*T.^2)./3+a11*T.^3+a12*T.^4+a13./T).*(R.*T)
        s=(a8*log(T)+a9*T+(a10*T.^2)./2+(a11*T.^3)./3+(a12*T.^4)./4+a14)*R
    end
end
end

Already,we know how to write a function and all the stuff....

write a function to calculate molecular weight:

formulla:

molecular_wt=atomic_wt*mole

function molecular_weight = molecular_func(species_name)

elements = ['H','O','C','N','A'];
molecular_weight = 0;
atomic_weight = [1,12,14,16,40];

for i = 1:length(species_name)

   for j = 1:length(elements)
 
          if(strcmp(species_name(i),elements(j)))

             molecular_weight = molecular_weight + atomic_weight(j);
             position = j;
          end

   end

   n = str2double(species_name(i));
       if(n>1)
        molecular_weight = molecular_weight + atomic_weight(position)*(n-1);
       end
   
end

disp(molecular_weight)
end

In this function we find the element and compare with atmic wt throungh the strcmp.finally, we calculate molecular weight.

 

MAIN PROGRAM:

 

clear all
close all
clc
%FILE PARSING
%INPUTS WHICH TAKEN BY THE FILE..
f1=fopen('THERMO.dat','r');
%fgetl func is used to pick first line of the file...
s1=fgetl(f1);
s2=fgetl(f1);
%strcmp func split the numbers
S=strsplit(s2,' ');
%str2num func is used to convrt string to number
temp1=str2num(S{2});
temp2=str2num(S{3});
temp3=str2num(S{4});
%avoiding header lines
for i=1:3
s3=fgetl(f1);
end
%extract the co-eff and calc cp,h,s,molecular weight
for j=1:53
tline=fgetl(f1);
species_name=tline(1:6);
%strfind func is used to find string in that line
line1=fgetl(f1);
l=strfind(line1,'E');
a1=str2num(line1(1:l(1)+3));
a2=str2num(line1(l(1)+4:l(2)+3));
a3=str2num(line1(l(2)+4:l(3)+3));
a4=str2num(line1(l(3)+4:l(4)+3));
a5=str2num(line1(l(4)+4:l(5)+3))

line2=fgetl(f1);
m=strfind(line2,'E');
a6=str2num(line2(1:m(1)+3));
a7=str2num(line2(m(1)+4:m(2)+3));
a8=str2num(line2(m(2)+4:m(3)+3));
a9=str2num(line2(m(3)+4:m(4)+3));
a10=str2num(line2(m(4)+4:m(5)+3));


line3=fgetl(f1);
n=strfind(line3,'E');
a11=str2num(line2(1:n(1)+3));
a12=str2num(line2(n(1)+4:n(2)+3));
a13=str2num(line2(n(2)+4:n(3)+3));
a14=str2num(line2(n(3)+4:n(4)+3));
%open a folder and named as that species name,what ever it is...
mkdir(['C:\Users\AMEER\OneDrive\Documents\MATLAB\Examples\R2021b\matlab\FindArrayElementsThatMeetAConditionExample\',species_name])
T=linspace(temp1,temp3,500);
%calling a func(cp,h,s)
[cp,h,s]=co_eff(a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13,a14,T);


%plotting
figure(1)
plot(T,cp,'LineWidth',3)
title(species_name)
xlabel('temperature')
ylabel('specific heat')

figure(2)
plot(T,h,'LineWidth',3)
title(species_name)
xlabel('temperature')
ylabel('enthalpy')

figure(3)
plot(T,s,'LineWidth',3)
title(species_name)
xlabel('temperature')
ylabel('entropy')

directory = (['C:\Users\AMEER\OneDrive\Documents\MATLAB\Examples\R2021b\matlab\FindArrayElementsThatMeetAConditionExample\',species_name]);
cd(directory)
%save that figure in particular folder
saveas(1,'specific_heat.png')
saveas(2,'enthalpy.png')
saveas(2,'entropy.png')
%change tha folder path
cd('C:\Users\AMEER\OneDrive\Documents\MATLAB\Examples\R2021b\matlab\FindArrayElementsThatMeetAConditionExample\')
%calling a func(molecular wt)
molecular_weight = molecular_func(species_name)
end
fclose(f1)

outputs:

  • The spicies name's folder has been created by the code

  • Molecular weight has been displayed in the command window.

o2

N2

 

CO2


Comments

Post a Comment

Popular posts from this blog

simple pendulum using matlab

-
Image
  ORDERLY DEFFRENTIAL EQUATION: In applied engineering and science many physical initial and boundary value problems can be described by ODE. The behavior of many physical problem,particularly those system undergoing time dependent changes can be formulated as ODE. The diff.equ involves one or more physical quantities such as displacement,velocity,temperature etc.. An nth order ODE may be represented as F(x,y,y’,y’’…)=0 SIMPLE PENDULUM: A simple pendulum is a mechanical device,the thin string or massless string is fixed in a supported beam and the other end is fixed to a bob.actually,its act like a newton’s 2 nd law,we give a force to the bob,the pendulum swig due to net force and mass.according to the net force and mass , the bob travels high or low. THE EQU OF A SIMPLE PENDULUM: b-damping co-efficient m-mass of the object g-gravity l-length og the string   First of all,we change the 2 nd order ODE to 1 ST order ODE  According to this formulla,create a ...
Read more

DRAG FORCE

-
Image
DRAG FORCE   DRAG FORCE: Drag force is a force,which resistance offered by the fluid like air or water.if we consider a car travelling a point A to B,we put a hand outside thorugh the window,then we feel the opposive of air.its known as drag force.   DRAG FORCE EQUATION: When the drag takes place in the presence of air, it is known as aerodynamic drag. When it takes place in water, it is known as hydrodynamic drag.   F D  = ½c_dρAv 2   Where, C_d-Drag co-efficient,which is dimensionless with value,depending upon the shape of object.(c_d=0.33-0.41 for mini cars and c_d=0.31-0.40 for sport cars etc...). ρ  -rho is the density of fluid. V – velocity of the object. A – cross sectional area of object( In automonile 'A' consider as frontal area). UNIT OF DRAG FORCE:NEWTON or N Drag force is propotional to the square of the velocity.if the velocity of speed increases two times,drag force increases four times. EXAMPLE: A car travelling on a ...
Read more

2r robotic arm

-
Image
2R ROBOTIC ARM ABOUT: It is a type of mechanical arm,usually it’s perform a similar function like a human arm.in the 2r robotics arm,the two links and manipulator were connected by the joints and the joints allow to perform either rotational motion or linear displacement.the links of the manipulator can be consider to form a kinematic chain To create a animation using matlab: First of all,we should understand the concept of forward kinematics 2r robotics arm like how it’s works?,what are the material has been used?,how many linkages and joints has been used? Then,we should do analyse mathematically, ¨     Calculate length and position of link1.. ¨     Calculate length and position of link2.. For plotting purpose;we pick some co-ordinates,which refers to forward kinematics robotics arm. Starting orgin x0=0;y0=0 At link 1: x1=l1cos⍬;y1=l1sin⍬ At link 2:x2=x1+l2cos⍬;y2=y2+l2sin⍬ [x0 y0][x1 y1][x2 y2] write a code to create animation li...
Read more