Week 4.1 - Genetic Algorithm

clear all
close all
clc

%input
gamma = 1.4;
t3 = 2300;    %in kelvin

%state variable 
p1 = 1;    %in pa
t1 = 500;       %in kelvin

% engine geometeric parameters
bore = 0.1;          %in meters
stroke = 0.1;        %in meters
con_rod = 0.15;      %in meters
cr = 12;

%calcuating swept volume and clearance volume 
v_s = (pi/4)*bore^2*stroke;
v_c  = v_s/(cr-1);
 v1 = v_s+v_c 
 v2 = v_c
 
 %state variable at state point 2
 %p2v2^gamma = p1v1^gamma
 %p2 = p1*(v1/v2)^gamma = p1*cr^gamma
 p2 = p1*cr^gamma;
 
 
 %p1v1/t1 = p2v2/t2 | t2 = p2*v2*t1/(p1*v1)
 t2 = p2*v2*t1/(p1*v1);

 constant_c=p1*v1^gamma;
 V_compression=engine_kinematics(bore,stroke,con_rod,cr,0,180);
 P_compression=constant_c./V_compression.^gamma;
 
 
 
 %state variable at state point 3
  v3 = v2;
 
 %p3v3/t3 = p2v2/t2 | p3 = p2*t3/t2
 p3 = p2*t3/t2;
 
 constant_c=p3*v3^gamma;
 V_expansion=engine_kinematics(bore,stroke,con_rod,cr,0,180);
 P_expansion=constant_c./V_expansion.^gamma;
 

 %state variable at state point 4
 v4 =v1
 
 %p3v3^gamma = p4v4^gamma | p4 = p3(v3/v4)^gamma
 p4 = p3*(v3/v4)^gamma;
 t4=p4*v4*t3/(p3*v3);
 
 %thermal efficiency
 thermal_efficiency = 1-((cr)^(1-gamma))
 

 hold on
 plot(v1,p1,'*','color','r')
 plot(V_compression,P_compression)
 plot(v2,p2,'*','color','r')
 plot(v3,p3,'*','color','r')
 plot(V_expansion,P_expansion)
 plot(v4,p4,'*','color','r')
 plot([v2 v3],[p2 p3],'color','b')
 plot([v4 v1],[p4 p1],'color','b')