%A PROGRAM TO SOLVE OTTO CYCLE USING MATLAB
clear all
close all
clc
%Inputs
gamma = 1.4
t3 = 2300
%State variables
p1 = 101325
t1 = 500
%Engine geometric parameters
Bore = 0.1;
Stroke = 0.1;
Con_rod = 0.15;
cr = 12;
%Calculating the swept volume and clearance volume
v_swept = (pi/4)*Bore^2*Stroke;
v_clearance = v_swept/(cr-1);
v1 = v_swept + v_clearance;
v2 = v_clearance
%State variables at point 2
%p2v2^gamma = p1v1^gamma
%p2 = p(v1/v2)^gamma
p2 = p1*cr^gamma
%p1v1/t1 = p2v2/t2
%t2 = p2*v2*t1/(p1*v1)
t2 = p2*v2*t1/(p1*v1)
%State variables at point 3
v3 = v2
%p3v3/t3 = p2v2/t2
p3 = p2*t3/t2
%State variables at point 4
v4 = v1
%p3v3^gamma = p4v4^gamma
p4 = p3*(v3/v4)^gamma
%plot
plot(p1,v1,'*','color','r')
plot(p2,v2,'*','color','r')
plot(p3,v3,'*','color','r')
plot(p4,v4,'*','color','r')