Air Standard Cycle

clear all
close all
clc

%Inputs
p1 = 1;
t1 = 500;
gamma = 1.4;
t3 = 2300;

%Engine geometery 
bore = 0.1;
stroke = 0.1;
cr_length = 0.15;
rk = 12;

%volumes
sv = (pi/4)*((bore)^2)*stroke;
cv = sv/(rk-1);

v1 = sv + cv;
v2 = cv;
v3 = v2;
v4 = v1;

const = p1*v1^gamma;
insvolcomp = volume_calc(bore,stroke,cr_length,rk,0,180); %IMPORTANT
insprcomp = const./(insvolcomp.^gamma);

%point 2 
p2 = ((rk)^gamma)*p1;
t2 = ((rk)^(gamma-1))*t1;

%point 3
p3 = (t3/t2)*p2;

constk = p3*v3^gamma;
insvolexp = volume_calc(bore,stroke,cr_length,rk,180,0);
insprexp = constk./(insvolexp.^gamma);

%point 4 
p4 = p3/((rk)^gamma);
t4 = t3/((rk)^(gamma-1));

efficiency = (1 - 1/((rk)^gamma-1))*100;

figure(1)
hold on
xlabel('Volume')
ylabel('Pressure')
plot(v1, p1, '*', 'color', 'b')
plot(v2, p2, '*', 'color', 'b')
plot([v2 v3], [p2 p3], 'color', 'b')
plot(insvolcomp, insprcomp)
plot(v3, p3, '*', 'color', 'b')
plot(v4, p4, '*', 'color', 'b')
plot([v4 v1], [p4 p1], 'color', 'b')
plot(insvolexp, insprexp)