Week 3 - Solving second order ODEs

import numpy as np import matplotlib.pyplot as plt from scipy.optimize import curve_fit # curve fit function def func(t, a, b, c): return a*pow(t,2)+b*t+c # reading thermodynamics file def read_file(): temperature=[] cp=[] for line in open('data','r'): values=line.split(',') temperature.append(float(values[0])) cp.append(float(values[1]))…

# ODE solver Pendulum Problem import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt import math def model(theta,t,b,g,l,m): # function definition theta1=theta[0] theta2=theta[1] dtheta1_dt=theta2 # Second order diffierential equation simplified dtheta2_dt=-(b/m)*theta2-(g/l)*math.sin(theta1)…

""" Otto Cycle Simulator """ import math import matplotlib.pyplot as plt def engine_kinematics (bore,stroke,conn_rod,cr,start_crank,end_crank): """ Engine Kinematics """ # Geometric paramters a=stroke/2 R=conn_rod/a # Volume parameters V_s = math.pi*(1/4)*pow(bore,2)*stroke V_c = V_s/(cr-1) sc=math.radians(start_crank)…

Please chech the attached

Please find the attached

Please find the attached excel file.

Cyclone separator Introduction A division of multiphase flows is the discrete phase model. The difference between the Lagrangian and Eulerian views is that the Eulerian perspective assumes a finite volume element in the fluid flow route while the Lagrangian view bases its analysis of fluid behavior on particle tracking…

Parametric study on gate valveGate valve is a device used to control the fluid flow.The flow coefficient - Cv - or the flow factor - Kv - are commonly used to specify capacities of control valves.Flow Coefficient - CvThe flow coefficient - Cv - is based on the imperial units system and is defined as the flow of water through…

https://screenrec.com/share/gSN9p2LeUh

https://screenrec.com/share/ldjs3xfvLz   Introduction Flow around the Ahmed body is a well-known CFD benchmark for external aerodynamics. The Ahmed body is a generic car model first proposed by [Ahmed et al., (1984)]. All ground vehicles can be termed as bluff bodies moving close to the road surface. The analyzed…

https://screenrec.com/share/RPx6tnygAj https://screenrec.com/share/AXaF6BgWbc

https://screenrec.com/share/sKPVaxz7RN https://screenrec.com/share/zd4qtEB7yT

https://screenrec.com/share/LB14XzOKRV                                                                            …

https://screenrec.com/share/iGASwjCMeQ clear close clc f=fopen('THERMO.dat','r'); % file open for i= 1:5 % skipping the header fgetl(f); end for j=1:53 % using the the command for loop to pick each substance data % extraction of local temperatures a=fgetl(f); % picking the Oxygen temperatures line b=strsplit(a,' '); %…

https://screenrec.com/share/KCjfEtreIg clear close clc fprintf('RANKINE CYCLE SIMULATOR') fprintf('1-2: Isentropic Expansion in the Turbine') fprintf('2-3: Constant Pressure Heat Rejection in the Condenser') fprintf('3-4: Isentropic Compression in the Pump') fprintf('4-1: Constant Pressure Heat Addition in the Boiler')…

https://youtu.be/iS-55hSRgls https://screenrec.com/share/YdmOqrhL5B   clear close clc % inputs m=1 b=0.05 l=1 g=9.81 % initial condition theta_0=[0,3] % time points t_span=linspace(0,20,300) % solve ode [t,results]=ode45(@(t,theta)ode_fun(t,theta,b,g,l,m),t_span,theta_0) figure(1) plot(t,results(:,1)) hold on plot(t,results(:,2))…

https://screenrec.com/share/tdMVfiCzwp https://youtu.be/MZ6dWjCifMY clear all close all clc % inputs L1= 1; L2= 0.5; Theta1= linspace(0,90,10); Theta2= linspace(-180,180,50); ct=1; for i=1:length(Theta1); THETA1=Theta1(i); for j=1:length(Theta2); THETA2=Theta2(j); x0=0 y0=0 x1=L1*cosd(THETA1); y1=L1*sind(THETA1); x2=x1+L2*cosd(THETA2);…