Air standard Cycle Matlab

AIM:- Write MATLAB code that can solve or check the engine kinematics of an engine operating on the OTTO cycle and make plots for it. Create a PV diagram and also output the thermal efficiency of the engine.

DESCRIPTION:- An Otto cycle is an idealized thermodynamic cycle that describes the functioning of a typical spark ignition piston engine. It is the thermodynamic cycle most commonly found in automobile engines.

The Otto cycle is a description of what happens to a mass of gas as it is subjected to changes of pressure, temperature, volume, addition of heat, and removal of heat. The mass of gas that is subjected to those changes is called the system. The system, in this case, is defined to be the fluid (gas) within the cylinder. By describing the changes that take place within the system, it will also describe in inverse, the system's effect on the environment. In the case of the Otto cycle, the effect will be to produce enough net work from the system so as to propel an automobile and its occupants in the environment.

PROCESS:-

• Process 0–1 a mass of air is drawn into piston/cylinder arrangement at constant pressure.
• Process 1–2 is an adiabatic (isentropic) compression of the charge as the piston moves from bottom dead centre (BDC) to top dead centre (TDC).
• Process 2–3 is a constant-volume heat transfer to the working gas from an external source while the piston is at top dead centre. This process is intended to represent the ignition of the fuel-air mixture and the subsequent rapid burning.
• Process 3–4 is an adiabatic (isentropic) expansion (power stroke).
• Process 4–1 completes the cycle by a constant-volume process in which heat is rejected from the air while the piston is at bottom dead centre.
• Process 1–0 the mass of air is released to the atmosphere in a constant pressure process

GOVERNING EQUATIONS:-

Compression ratio (r) is the ratio of total cylinder volume to the clearance volume.

r = V1V2=r=V3V4">V1V2=r=V3V4V1V2=r=V3V4

We are going to find efficiency of system by using T_3,T_4.

Process 1-2:

This process is an isentropic (reversible adiabatic) process. For this process, the relation between T and V is as follows:

T2T1=(V1V2)γ-1=rγ-1…………..1">T2T1=(V1V2)γ−1=rγ−1…………..1T2T1=(V1V2)γ-1=rγ-1…………..1

Process 3-4:

This is also an isentropic process. The relation between T and V in this process is similar to the relation between T and V in process 1-2:

Here,

T3T4=(V3V4)γ-1=rγ-1……………2">T3T4=(V3V4)γ−1=rγ−1……………2T3T4=(V3V4)γ-1=rγ-1……………2

Air-standard efficiency of Otto cycle:

It is defined as the ratio between work done during Otto cycle to the heat supplied during Otto cycle.

So, nth=1-1rγ-1">nth=1−1rγ−1nth=1-1rγ-1 ,,,,,,by Eq 1,2

The following equation shows the relationship between crank angle and volume inside the crank chamber: -

VVc=1+0.5⋅(rc-1)⋅[R+1-cos⁡[θ-(R2-(sin⁡θ)2)0.5]">VVc=1+0.5⋅(rc−1)⋅[R+1−cos[θ−(R2−(sinθ)2)0.5]VVc=1+0.5⋅(rc-1)⋅[R+1-cos⁡[θ-(R2-(sin⁡θ)2)0.5]

coding:-

Coding: -First we have to define a function named as engine kinematics:


function [V] = engine_kinematics_otto(bore, stroke, c_rod, cr, start_crank, end_crank)
%defining the parameters of engine
a = stroke/2;
R = c_rod/a;
 
v_s = (pi/4)*bore^2*stroke;
v_c = v_s/(cr-1);
 
theta = linspace(start_crank,end_crank,100);
% here x,y,z are the three parts of equations 4(report)
x = 0.5*(cr-1);
y = R + 1 -cosd(theta);
z = (R^2 - sind(theta).^2).^0.5;
 
V = (1 + x*(y - z)).*v_c;
 
end

main code for finding thermal efficiency
clear all
close all
clc
 
% Inputs
gamma = 1.4;
t3 = 2300;
 	
% State Variables
p1 = 101325;
t1 = 500;
 
% Engine gepmetric parameters 
bore = 0.1;
stroke = 0.1;
c_rod = 0.15;
cr = 12;
 
% Calculating the swept volume and the clearance volume
v_swept = (pi/4)*bore^2*stroke;
v_clearance = v_swept/(cr-1);
 
v1 = v_swept + v_clearance;
v2 = v_clearance;
 
% State varaible 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_otto(bore, stroke, c_rod, cr, 180, 0);
P_compression = constant_c./V_compression.^gamma;
 
 
% State varaible 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_otto(bore, stroke, c_rod, cr, 0, 180);
P_expansion = constant_c./V_expansion.^gamma;
 
% State varaible at state point 4
v4 = v1;
 
% p3v3^gamma = p4v4^gamma | p4 = p3(v3/v4)^gamma
p4 = p3*(v3/v4)^gamma;
 
% p4v4/t4 = p3v3/t3 | t4 = (p4*v4*t3)/(p3*v3)
t4 = (p4*v4*t3)/(p3*v3);
 
figure(1)
hold on
plot(v1,p1,'*','color','r')
plot(V_compression,P_compression,'color','g')
plot(v2,p2,'*','color','r')
plot(v3,p3,'*','color','r')
plot(V_expansion, P_expansion,'color','g')
plot(v4,p4,'*','color','r')
plot([v2 v3],[p2 p3],'color','g')
plot([v1 v4],[p1 p4],'color','g')
xlabel('Volume')
ylabel('Pressure')
title('OTTO CYCLE')
 
% Thermal efficiency
term_1 = t4-t1;
term_2 = t3-t2;
Thermal_efficiency = (1-(term_1/term_2))*100

output:-

Conclusion:-

CONCLUSION :- In the process 1–2 the piston does work on the gas and in process 3–4 the gas does work on the piston during those isentropic compression and expansion processes, respectively. Processes 2–3 and 4–1 are isochoric processes; heat is transferred into the system from 2—3 and out of the system from 4—1 but no work is done on the system or extracted from the system during those processes. No work is done during an isochoric (constant volume) process because addition or removal of work from a system requires the movement of the boundaries of the system; hence, as the cylinder volume does not change, no shaft work is added to or removed from the system.