Otto Cycle And Thermal Efficiency

clear all
close all
clc

%Engine Geometrics
bore=0.1;
stroke=0.1;
con_rod=0.15;
cr=12;
gamma=1.4;

v_swept=(pi/4)*(bore^2)*stroke;
v_clearance=v_swept/(cr-1);

v1= v_swept+ v_clearance;
v2=v_clearance;

p1=101325;
t1=500;

%calculating p2
p2=p1*(cr^gamma);

%calculating t2
t2=(p2*v2*t1)/(p1*v1);

con_c=p1*v1^gamma;
v_compression=engine_kinematics(bore,stroke,con_rod,cr,180,0);
p_compression=con_c./v_compression.^gamma;

v3=v2;
%calculating p3;
t3=2300;
p3=((p2*t3)/t2);

con_c=p3*v3^gamma;
v_expansion=engine_kinematics(bore,stroke,con_rod,cr,0,180);
p_expansion=con_c./v_expansion.^gamma;

v4=v1;

%calculating p4
p4=p3*((v3/v4)^gamma);



figure(1)
hold on
xlabel(\"volume\")
ylabel(\"pressure\")
plot(v1,p1,\'*\',\'color\',\'r\')
plot([v2 v3],[p2 p3])
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([v4 v1],[p4 p1])


efficiency=1-(1/(cr^(gamma-1)))