Curve fitting

An epicyclic gear train consists of two gears mounted so that the centre of one gear revolves around the centre of the other. A carrier connects the centres of the two gears and rotates to carry one gear, called the planet gear or planet pinion, around the other, called the sun gear or sun wheel. Following inputs are considered…

In this project-forward kinematics of a 2R robotic arm  we determine the workspace of a robotic arm with two links . The two links can have different constraints and lengths. In the code first we give the input parameters like length of arms [L1,L2],and the angle with which the arms can be rotated[theta1,theta2].…

IC Engine Valve Train A valvetrain or valve train is a mechanical system that controls operation of the valves in an internal combustion engine, whereby a sequence of components transmits motion throughout the assembly.  Conventional reciprocating internal combustion engine uses valves to…

In this project the analysis of the reciprocating motion of the piston is studied for various positions of the wrist pin axis with respect to the axis of rotation of the crank. The offset between these axis is made positive or negative and Linear displacement of the Piston Head is studied. Modeling of all the parts.  1.…

The Geneva drive 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 in the otherwheel (driven wheel) that advances it by one step at a time. …

In determining the minimum cushion pressure needed to break a given thickness of ice using an air cushion vehicle, Muller (“Ice Breaking with an Air Cushion Vehicle”) in Mathematical Modeling: Classroom Notes in Applied Mathematics, SIAM 1987) they derived the equation:- Equation is used to find the minimum…

  Compatability check The code checks whether the input file name given exsists with extension '.out' or not, if a non-compatible file is provided as an input. It shows something like "File not recognized. Please provide a valid CONVERGE output file" Data visualizer  Script takes column numbers as the input and…

    For Cubic ploynomial #Curve fitting #importing libraries import numpy as np import matplotlib.pyplot as plt from scipy.optimize import curve_fit #polynomial used in the program def function(t,a,b,c,d): return a*pow(t,3)+b*pow(t,2)+c*t+d #reads the file and returns the values of the temperature and specific…

In determining the minimum cushion pressure needed to break a given thickness of ice using an air cushion vehicle, Muller (“Ice Breaking with an Air Cushion Vehicle”) in Mathematical Modeling: Classroom Notes in Applied Mathematics, SIAM 1987) they derived the equation:- Equation is used to find the minimum…

An Otto cycle is an idealized thermodynamic cycle that describes the functioning of a spark ignition engine. The processes of Otto cycle are Process 1–2 is an adiabatic (isentropic) compression. Process 2–3 is a constant-volume heat transfer. Process 3–4 is an adiabatic (isentropic)…

A program in python to simulate the forward kinematics of a 2R Robotic Arm  In this program we determine the workspace of a robotic arm with two links . The two links can have different constraints and lengths. In the code First we give the input parameters like length of arms [l1,l2] The angle with which the arms…

Program in Matlab that will simulate the pendulum motion using second order Ordinary Differential Equation. solving second orders ODEs parameters used and values L=1 metre, m=1 kg, b=0.5. g=9.81 m/s2. Simulating the motion between 0-20 sec, for angular displacement=0,angular velocity=3 rad/sec2 at time t=0. clear all close…

Python program to calculate drag force against a cyclist. The factors involved are  Velocity-is an array of values Drag coefficient-cd=0.8 Frontal Area-area=0.2 Density of air-rho=1.2 Program outline Library function to create plots is first imported. The input parameters and the factors are defined. Preallocation…

Minimizing a simple non-linear function subjected to a constraint using Lagrange Multipliers. Given the non linear function as F=5-(x-2)^2 -2(y-1)^2 which after simplification becomes F=-X^2-2Y^2+4X+4Y-1 subjected to the following constraint G= x+4y-3 Thus the function and the constraint can be writtten using the lagranges…

In this challenge, the 3D model for a centrifugal pump is created and then flow simulation is run and various parameter variations are studied. Centrifugal pump design dimensions in mm Diameter of the Volute casing is 90 mm. Thickness of the Volute casing is 20 mm. Rotar blade diameter is from 40mm to 80mm. Diameter of…

In this project, we will create a 3D model of a flowbench and run flow analysis. And obtain a plot of lift vs mass flow rate, also find the pressure drop in the system. Flow bench is a device used to design the intake port. Usually in the flowbench we take the port which is required to be tested. Here in this project we…

    %Parameters drag_coeff=linspace(.5,1,25); %contact area A=.2 %density rho=1.2 %velocity v=linspace(1,2,25); dragforce=0.5*rho*A*v.*v.*drag_coeff; figure(1) plot(v,dragforce,'color','r') xlabel('velocity','color','r') ylabel('dragforce','color','r') figure(2) plot(drag_coeff,dragforce,'color','g') xlabel('drag_coeff','color','g')…

Simulating flow over a cylinder for different reynolds numbers To calculate the baseline case reynolds number. we have the following formula Re = rho*U_inlet*D/mu D = cylinder diameter mu = dynamic viscosity rho = density U_inlet=Inlet velocity Thus for Re = 20, 40 and 100 we are runing the transient flow simulation also…

Modeling of  flow over a NACA0017 air foil and to compare the lift and drag forces on the foil for angle of attacks 0,2,4,6,8,10. The NACA airfoil series is controlled by 4 digits NACA  MPXX which designate the camber, position of the maximum camber and thickness. If an airfoil number is e.g.NACA 1128…

OBJECTIVES To run a pipe flow simulation with an inlet Reynolds number of 100,1000 and 10,000. And for each of these cases 1. Place line probes at 95%, 90% and 85% of the pipe length. 2. Compare the normalized velocity profile at each of these locations. 3. Normalize the velocity profile by the inlet velocity. Reynolds…

RANKINE CYCLE SIMULATOR The Rankine cycle is a model used to predict the performance of steam turbine systems. It was also used to study the performance of reciprocating steam engines. The Rankine cycle is an idealized thermodynamic cycle of a heat engine that converts heat into…

PARSING NASA THERMODYNAMIC DATA NASA came up with polynomials that can be used to evaluate thermodynamic properties such as Cp, H and S using polynomials. They have also documented the co-efficients that are required to evaluate these polynomials for 1000+ species.The code below parses the therm.dat file to make use of…

Genetic algorithm is a method used to solve constrained or unconstrained optimization problems based on natural selection,the process is inspired by Charls darwins theory of evolution.It selects individuals at random from the current population as parents and use them to produce children for next generation.Elite children…

Code to fit a linear and cubic polynomial curve for the given Cp data. In order to obtain the curve fit i have used polyfit(x,y,r) and polyval(a,b) commands ,in the above commands \'r\'- degree of the polynomial.  In the code below i have used polyfit command three times with r=1,2&3       …

%Air standard Cycle [Otto cycle] %Create a PV diagram of otto cycle %Thermal efficiency of the engine for various compression ratio vales close all clear all clc %inputs gamma=1.4 t3=2000 %physical parameters bore_dia=0.1 stroke_leng=0.1 con_rod=0.1 cr=12 %calculate swept and clearance volume v_swept=(pi/4)*bore_dia^2*stroke_leng…