%otto_cycle_and_efficiency
close all
clear all
clc
% inputs
gamma=1.4
t3=2300
%state variabels
p1=101325
t1=500
%engine parameters
bore=0.1;
stroke=0.1;
cr=12;
con_rod=0.15;
%calculating the swept volume and clearencce volume
v_swept=(pi/4)*bore^2*stroke;
v_clearence = v_swept/(cr-1);
v1= v_swept + v_clearence
v2= v_clearence
%sate variable at stae point2
%p2v2^gamma=p1v1^gamma
%p2=p1*(v1/v2)^gamma=p1*cr^gamma
p2=p1*cr^gamma;
%p1v1/t1=p2v2/t2=p2*v2*t1/(p1*v1)
t2=p2*v2*t1/(p1*v1);
constant_c=p1*v1^gamma;
V_compression= piston_kinematics(bore,stroke,con_rod,cr,180,360);
P_compression=constant_c./V_compression.^gamma
% state varible at a state point3
v3=v2
%p3v3/t3=p2v2/t2 |p3=p2*(t3/t2)
p3=p2*(t3/t2)
constant_c=p3*v3^gamma;
V_expansion= piston_kinematics(bore,stroke,con_rod,cr,0,180);
P_expansion=constant_c./V_expansion.^gamma
%state varible at stae point 4
v4=v1
%p3v3^gama=p4v4^gamma | p4=p3*(v3/v4)^gamma
p4=p3*(v3/v4)^gamma
efficiency=1-(v2/v1)^0.4
figure(1)
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\')
figure(2)
plot([v1 v2 v3 v4],[p1 p2 p3 p4],\'color\',\'r\')