Code to create PV diagram of an otto cycle and to output its thermal efficiency of the engine
Objective: 1) Code to create a PV diagram clear all;clc;%inputgamma=1.4;t3=2300;%state variablesp1=101325;t1=500;%engine geometric parametersbore=0.1;stroke=0.1;con_rod=0.15;cr=12;%calculating swept volumev_swept=(pi/4)*bore^2*stroke;v_clearance=v_swept/(cr-1);v1=v_swept+v_clearance;v2=v_clearance;%state variablesp2=p1*cr^gamma;t2=p2*v2*t1/(p1*v1);constant_c=p1*v1^gamma;V_compression=engine_kinematics(bore,stroke,con_rod,cr,180,0);P_compression=constant_c./V_compression.^gamma;v3=v2;p3=p2*t3/t2;constant_c=p3*v3^gamma;V_expansion=engine_kinematics(bore,stroke,con_rod,cr,0,180);P_expansion=constant_c./V_expansion.^gamma;v4=v1;p4=p3*(v3/v4)^gamma;figure(1)hold…