Curve fitting using Python

Centrifugal Pump Design and Analysis : Skill-Lync (skill-lync.com)

Aim: To study the flow through a centrifugal pump, observe the fluid flow pattern and arrive at a pressure ratio   Background: Centrifugal pumps are used to transport fluids by the conversion of rotational kinetic energy to the hydrodynamic energy of the fluid flow. The rotational energy typically comes from…

Aim: To model the fluid flow through a flowbench with a valve and analysing the results obtained.   Background: A flowbench is a device used for testing internal aerodynamic qualities of an engine component. It is mainly used for testing intake and exhaust ports. Here, we are using a generalised flow bench geometry…

Aim: The goal of the project is to simulate the flow over a NACA symmetric airfoil for various angles of attack and interpret the results.   Background: An airfoil or aerofoil is the cross-sectional shape of a wing, blade, or sail. An airfoil-shaped body moving through a fluid produces an aerodynamic force. The component…

Aim: The goal of the project is to simulate the external flow on a cylinder for varying Reynolds Numbers and observe velocity and pressure profiles.   Background: Reynold's Number: The Reynolds number is the ratio of inertial forces to viscous forces within a fluid that is subjected to relative internal…

Aim: The goal of the project is to simulate the flow in a cylindrical pipe for multiple Reynolds Number values and observe the flow development and velocity profiles for the above cases.   Background: Reynold's Number: The Reynolds number is the ratio of inertial forces to viscous forces within a fluid…

Aim: The goal of the project is to understand the construction of an IC Engine valve train, model the same, carry out a multi-body simulation in Solidworks.   Background: An engine valve train is a mechanical system that operates the valves in an engine that allows for the fresh charge and exhaust gas flow in and…

Aim: The goal of the study is to simulate the piston movement and study the linear displacement of the piston for variations in the piston-conrod assembly.   Background: A piston is a component of reciprocating engines, reciprocating pumps, gas compressors, hydraulic cylinders, and pneumatic cylinders, among other…

Aim: The goal of the project is to conduct a multi-body simulation on a Planetary geartrain in solidworks.   Background: A planetary gear train, also known as an Epicyclic geartrain, consists of two gears mounted so that the center of one gear revolves around the center of the other. A carrier connects the centers…

Aim: To conduct a motion study for an internal gearing Geneva Mechanism.   Background:  The Geneva drive or Maltese cross is a gear mechanism that translates a continuous rotation movement into intermittent rotary motion. The rotating drive wheel is usually equipped with a pin that reaches into a slot located…

Aim: The goal of the project is to understand frequency analysis on a rotating body, the types, and the results that would be obtained from the said analysis. Introduction:  Natural frequencies are the frequencies with which a body oscillates when no repeated or continuous load is being acted upon it. Based on the…

Aim: The goal of the analysis is to understand the stress, strain, and displacement when a plat with a set number of through-holes is subjected to load.   Introduction: We are provided with 2 plates with different hole geometries to analyze for stress, strain, and displacement. The plates will be subjected to equal…

Aim: To perform buckling analysis on a cyclonic separate and decide the optimal height at which a support structure must be added to stiffen the structure.   Background: The goal of the project is to analyze the buckling that a structure undergoes upon application of a load. We will first be isolating the stand and…

Aim: To analyse beams of different geometries and suggest the best beam geometry based on stress, displacement, and strain.   Introduction: 3 beam geometries are analysed here by applying a load of 1500 N on a 400mm segment of the beam. The geometries used are: Solid rectangle section I-beam U-profile Following parameters…

Aim: To find the roots of the equation using Newton-Rhapson method in Python.   Introduction: Newton-Rhapson is an iterative solving technique used to arrive at the solution by continuous solution of the equation: We start with an assumption (here referred to as guess value) as x0 and proceed to find x1 using the…

Aim: Goal of the project is to develop a data visualiser tool in python for data analysis. Background: The input provided is an engine data file extracted from CONVERGE. The file consists of multiple data columns like pressure, volume, crank angle, viscosity, etc.. It is necessary to extract data from the file to analyse…

The project aims at finding the best fit curve for a given data and finding the error between the original and said curve. The data being used is that of Specific heat (Cp) vs Temperature. The plot will be made to show the variation of Cp with Temp first.  Code: import numpy as nm from numpy import array from scipy.integrate…

Project aims at finding the best fit curve for a given data and finding the error between the original and said curve. The data being used is that of Specific heat (Cp) vs Temeperature. Plot will be made to show variation of Cp with Temp first.  Code: % Data cpdata = load('data'); temp = cpdata(:,1); Cp = cpdata(:,2);…

This project explores the simulation of the motion of a pendulum for the equation of damped vibration of a pendulum. Below is the equation (taken from "Challenges" page): The above equation is solved using ODE function, and then odeint is used to get an array of results vs time (position vs time). These positions…

Project involves plotting an Otto Cycle graph in MATLAB The goal is to plot all strokes for a petrol operating by the otto cycle from intake, compression, power to exhaust, including both adiabatic curves. Description: The Otto Cycle is the idealised cycle that describes the functioning of a typical spark ignition piston…

Objective: The project involves simulating a 2-joint robotic arm.   Code and Description: import math import matplotlib.pyplot as plt import numpy as nm # Lengths l1 = 1 l2 = 0.5 # Angles ang1 = nm.radians(nm.linspace(0, 90, 19)) ang2 = nm.radians(nm.linspace(0, 90, 19)) # Origin x0 = 0 y0 = 0 count = 1 for a1 in…

Objective: To calculate the drag force for various drag coefficient and velocity values for a certain geometry.   Introduction: Drag force is the force exerted by the air as the vehicle moves through it. The force exerted is directly proportional to the square of velocity, so as the vehicle accelerates and velocity…

Aim: To parse the NASA Thermodynamc Properties data file, calculate thermodynamic properties and plot the characteristics. Introduction: NASA's thermodynamic properties file consists of multiple species of gases, their respective operating temperature ranges, and coefficients to calculate thermodynamic properties.…

The project aims at optimization of a stalagmite function and finding the global maxima. For this, we will be using the concept of genetic algorithm. Following are the steps used to optimise the function: Code the function given in a separate function code. This is the function to be optimized. Define the search space…

This project explores the simulation of the motion of a pendulum for the equation of damped vibration of a pendulum. Below is the equation (taken from "Challenges" page): An attempt is made to solve the above equation, and then call it using ode45 to get an array of results vs time (position vs time). These positions are…

Project involves plotting an Otto Cycle graph in MATLAB The goal is to plot all strokes for a petrol operating by the otto cycle from intake, compression, power to exhaust, including both adiabatic curves. Description: The Otto Cycle is the idealised cycle that describes the functioning of a typical spark ignition piston…

Code and Description: The project involves simulating a 2-joint robotic arm. % Lengths l1 = 1; l2 = 0.5; % Angles ang1 = linspace (0, 90, 20); ang2 = linspace (0, 90, 20); count = 1; % loop for i = 1:length(ang1) ANG1 = ang1(i); for j = 1:length(ang2) ANG2 = ang2(j); x0 = 0; y0 = 0; x1 = l1*cosd(ANG1); y1 = l1*sind(ANG1);…