This the program to plot P-V diagram of otto cycle with it's efficiency. In otto cycle heat addition and heat rejection takes place at constant volume and compression and expansion takes place by reversible adiabatic process. p1,p2,p3,p4 are the pressure\'s at sec 1,2,3 and 4, v1,v2,v3,v4 are the volume\'s at at sec…

#Pipe Dimension\'s Pipe length(L) = 500mm = 0.5m Pipe Diameter(D) = 50mm = 0.05m   #Fluid properties Working fluid = water Density(rho) = 1000(kg/m^3) Dynamic viscosity(meu) = 8.90*10^-4 Pa.s #Boundary condition\'s Inlet velocity Outlet : static pressure Wall type : real wall   #Computation Renoldy\'s number(Re)=rho*V*D/meu…

This the Program for the crank and connecting rod anniamtion for Diesel and Petrol Engines using Matlab As we know that the rotation of the Crankshaft for Petrol engine is in the Clockwise Direction and anticlockwise for Diesel engine with respect to that I have written this program. You can check the animation Directly…

  %Drag force vs the velocity close all clear all clc   %constants   c_d=0.3   %density kg/m^3   rho=1000   %Velocity m/s   Vel=[20 30 40 50]   %frontal area m^2   A=0.9   %Dragforce newton   F_D=0.5*(Vel.*Vel)*A*rho*c_d   %ploting the valuesclear all   plot(Vel,F_D,\'linewidth\',3,\'color\',\'b\')…

b=damping coefficient., g = gravity in m/s2, L = length of the pendulum in m, m = mass of the ball in kg,   clear all close all clc %constants b=0.5 g=9.81 l=1 m=1 %intial condition theta_0=[0;5] %time t_span=linspace(0,20,100); %solve ode [t,x] = ode45(@(t1,theta) exam(t1,theta,b,g,l,m),t_span,theta_0); %origin O=[0,0];…

This the program to find out the root of any function using Newton-Raphson method. You need to just specify a fuction w.r.t 'x' ,intial condition(x0) and the no.of itertations max no.of iterations is 3 and min is 1. the final answer will be upto 3 decimal places. clear all close all clc %diff w.r.t to x syms x; %function…

clear all close all clc   %state variables p1 = 10000 t1=313 gama=1.4 t3=2150 %engine parameters bore=0.12; stroke=0.2; con_rod=0.7; cr=6;   %calculating clerence volume and stroke volume swept_vol=(3.14/4)*bore*bore*stroke clearence_vol=(swept_vol/(cr-1))   v1=clearence_vol+swept_vol   v2=clearence_vol…

%to find row vector A = [1 2 3 4 ] B= [3 4 5 6] %req row vector A.*B=[3 8 15 24] %to find coloum vector C=[1       2       3      4]   D=[3       4      5      6] %req coloum vector C.*D=[3          …

https://youtu.be/m-trexRaji4   clear all close all clc   %constants b=0.5 g=9.81 l=1 m=1     %intial condition   theta_0=[0;5] %time t_span=linspace(0,20,100);   %solve ode   [t,x] = ode45(@(t1,theta) exam(t1,theta,b,g,l,m),t_span,theta_0);   %origin O=[0,0]; axis(gca,\'equal\');…

This is the program to prase the NASA thermodynamic file and calculate the required thermodynamic properties. In main function,at first I am reading the file using fopen and then reading the first string and similarly reading the temperatures using fgetl command and converting it to numbers using str2num command,Since…