Week 5 - Curve fitting

"""Task to write a script basic performance of Engine with given engine_data Basic performance calculation Calculate the area under the P-V diagram. P is in column 2 and V is in column 8.Calculate the power output of this engine. Assume that RPM is 1500This engine consumed 20 micro grams of fuel (per 1 stroke cycle). Calculate…

import numpy as np import matplotlib.pyplot as pltfrom scipy.optimize import curve_fit #Curve fit funtion for linear polynomialdef func(t,a,b): return a*t +b #reading the thermodynamic data fliedef read_file(): temperature =[] cp =[] for line in open('data','r'): values =line.split(',') temperature.append(float(values[0]))…

https://youtu.be/x98n8BcsuaU   import numpy as npfrom scipy.integrate import odeintimport matplotlib.pyplot as pltimport math # function that returns dz/dtdef model(theta,t,b,g,l,m): theta1 = theta[0] #displacement theta2 = theta[1] #velocity dtheta1_dt = theta2 dtheta2_dt = -(b/m)*theta2-(g/l)*math.sin(theta1) dtheta_dt…

"""Otto_cycle"""import mathimport matplotlib.pyplot as plt def engine_kinematics(bore, stroke,con_rod, cr, start_crank, end_crank): #geometric parameters a = stroke/2 R = con_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) ec= math.radians(end_crank) num_values…