Genetic Algorithm to solve Stalagmite function using MATLAB

Grid  Dependency  Test:  For Valve Lift = 5 mm and Inlet Total Pressure = 303975Pa, three different Mesh SetUp were solved. Mesh SetUp:           (1)  Number of cells in X Axis  =  16                  …

1. Opening the \'therm.dat\' file in reading mode. 2. Ignoring the title line, and storing the global max, mid and min in parameters. 3. Defining temperture range in between global max and min values. 4. Ignoring the comments (that starts with \'!\'). 5. Taking the species namein \'inp\' variable. 6. Searching the species…

import matplotlib.pyplot as mol import numpy as np import sys import math import time value=[] dx=[] k=0 W=[] Work=0 t=0 i=1 res=\'y\' try: f=open(\'engine_data.out\',\'r\') except FileNotFoundError: print(\'File not recognized. Please provide a valid CONVERGE output file\') sys.exit(1) for line in f: if \'#\' not in line:…

Flow Simulation over a NACA airfoil is carried out at a velocity of 600m/s in x-direction, for different angle of attack, viz. 0,2,4,6,8 and 10. In the following model, only 0 degree angle of attack is shown(but for results it was changed over and over again).  https://drive.google.com/file/d/11rI_5kAaZtmLemt56M0zAcReNxOLuJ8w/view?usp=sharing…

INPUTS: Fluid = Air   Fluid Properties: Density,ρ= 1.225 kg/m3 Dynamic Viscosity, μ=1.81×10−5 kg/(m. s) OR 18.1 μPa   Deminsions of the Pipe: Diameter, D = 5mm Length, L = 25mm   Meshing: Number of Grids in X direction (Length) = 50 Number of Grids in Y direction = 10 Number…

INPUTS: Fluid = Air   Fluid Properties: Density, ρ= 1.225 kg/m3 Dynamic Viscosity, μ=1.81×10−5 kg/(m. s) OR 18.1 μPa   Deminsions of the Pipe: Diameter, D = 5mm Length, L = 25mm   Meshing: Number of Grids in X direction (Length) = 50 Number of Grids in Y direction = 10 Number…

1. What does popt and pcov mean? A. On using curve_fit command, it returns the optimum coefficients in the function defination which is stored in popt; and pcov stores the estimated covariance of popt in a matrix form (the diagonals provides the variance). 2. What does np.array(temperature) do? A. This command will call…

#otto cylce #PV plot and thermal efficiency import matplotlib.pyplot as mol import math #inputs p1=101325 t1=500 t3=2300 gamma=1.4 #geometric paprameters bore=0.1 stroke=0.1 con_rod=0.15 cr=12 a=stroke/2 R=con_rod/a #volume vs=(math.pi/4)*pow(bore,2)*stroke vc=vs/(cr-1) v1=vs+vc def engine_kinematics(bore,stroke,con_rod,cr,start,end):…

import matplotlib.pyplot as asd import math l1=1 l2=1.5 n=10 start=0 end=math.pi theta1=[] theta2=[] for i in range(0,n): temp=start+i*(end-start)/(n-1) theta1.append(temp) theta2.append(temp) x0=0 y0=0 ct=0 for i in theta1: for j in theta2: x1=l1*math.cos(i) y1=l1*math.sin(i) x2=x1+l2*math.cos(j) y2=y1+l2*math.sin(j)…

#program to calculate drag force of a cyclist import matplotlib.pyplot as asd #inputs #drag coefficient array cd=[0.8,1,1.2,1.4,1.6,1.8,2] #frontal area(m^2) A=0.2 #density of air(kg/m^2) den=1.25 #velocity array v=[1,2,3,4,5,6,7,8,9,10] dragforce_1=[] dragforce_2=[] for a in v: dragforce_1.append(0.5*A*den*cd[0]*a*a)…

Genetic Algorithm is an optimisation technique which works on natural selection of randomly generated values. It first selects a particular set of values from the total values, then performs some function to it, to produce optimised value, and then repeat the same procedure to some extend.  In MATLAB, ga is inbuilt…

Flow Simulation over a NACA airfoil is carried out at a velocity of 600m/s in x-direction, for different angle of attack, viz. 0,2,4,6,8 and 10. In the following model, only 0 degree angle of attack is shown(but for results it was changed over and over again).  https://drive.google.com/file/d/11rI_5kAaZtmLemt56M0zAcReNxOLuJ8w/view?usp=sharing…

clear all close all clc l1 = 1; l2= 2; x0=0; y0 = 0; theta1=linspace(0,180,20); theta2=linspace(0,180,20); n=1; for i=1:length(theta1) for j=1:length(theta2) x1=l1*cosd(theta1(i)); y1=l1*sind(theta1(i)); x2=x1+l2*cosd(theta2(j)); y2=y1+l2*sind(theta2(j)); plot([x0 x1],[y0 y1],[x1 x2],[y1 y2]) axis([-3.5 3.5 -0.5 3.5])…

https://projects.skill-lync.com/projects/Pipe-Flow-Simulation-using-Solidworks-83728  

Air is flown through a pipe of diameter 5mm and length 25mm. The simulation is studied assuming density of air 1.225 kg/m^3 and dynamic viscosity 0.00181 Pa-s.  Three probes at 85%,90% and 95% of the pipe length are drawn. https://drive.google.com/file/d/1oNWRh2pgLm76yixQFUVjcc8wCW3lFsM9/view?usp=sharing For Renold's…

#otto cylce #PV plot and thermal efficiency import matplotlib.pyplot as mol import math #inputs p1=101325 t1=500 t3=2300 gamma=1.4 #geometric paprameters bore=0.1 stroke=0.1 con_rod=0.15 cr=12 a=stroke/2 R=con_rod/a #volume vs=(math.pi/4)*pow(bore,2)*stroke vc=vs/(cr-1) v1=vs+vc def engine_kinematics(bore,stroke,con_rod,cr,start,end):…

Genetic Algorithm is an optimisation technique which works on natural selection of randomly generated values. It first selects a particular set of values from the total values, then performs some function to it, to produce optimised value, and then repeat the same procedure to some extend.  In MATLAB, ga is inbuilt…

import matplotlib.pyplot as asd import math l1=1 l2=1.5 n=10 start=0 end=math.pi theta1=[] theta2=[] for i in range(0,n): temp=start+i*(end-start)/(n-1) theta1.append(temp) theta2.append(temp) x0=0 y0=0 ct=0 for i in theta1: for j in theta2: x1=l1*math.cos(i) y1=l1*math.sin(i) x2=x1+l2*math.cos(j) y2=y1+l2*math.sin(j)…

#program to calculate drag force of a cyclist import matplotlib.pyplot as asd #inputs #drag coefficient array cd=[0.8,1,1.2,1.4,1.6,1.8,2] #frontal area(m^2) A=0.2 #density of air(kg/m^2) den=1.25 #velocity array v=[1,2,3,4,5,6,7,8,9,10] dragforce_1=[] dragforce_2=[] for a in v: dragforce_1.append(0.5*A*den*cd[0]*a*a)…

https://drive.google.com/file/d/1JbfT2UQf9-EW467pybZmJrPUwDq-ot7S/view?usp=sharing  

1. https://drive.google.com/file/d/1Wwf88XvFT53zRFlyQtYJs9lQsmW3YWFb/view?usp=sharing 2. Higher the order of equation between x and y, the better is the curve.  As the precision of direction change for higher order equation is much high thereby resulting better curve fit. 

https://drive.google.com/file/d/1xbj6DWBueld20C9PI0G1JbOiz2jLij8v/view?usp=sharing  

https://drive.google.com/file/d/13J93G5_mtoDviuk_wxLFWbl4EChWvgNG/view?usp=sharing https://drive.google.com/file/d/1_OMO-mph098wNzNX8iSolIgHjXDKwuZS/view?usp=sharing https://drive.google.com/file/d/17HXPXW6LgTTSRCDSx2iTnOqc8f1Moywy/view?usp=sharing  

https://drive.google.com/file/d/1tS8QP7U9zaQ547MDUGFyZS-o6EICal1G/view?usp=sharing

https://drive.google.com/file/d/1-sGjTouT9-lgAuyxYo211Fs3H-EfLPpP/view?usp=sharing

https://drive.google.com/file/d/1dXkQuShZtcfbND1Y9c7Saaf5Egs-YVlY/view?usp=sharing

https://drive.google.com/drive/search?q=week%201

%program to calculate drag force on a vehicle clear all close all clc %density (kg/m^3) d=.75 %frontal area (m^2) a=.8 %air velocity (m/s^2) v=[1:0.5:5] %drag coeff (unitless) cd=[.04 .09 0.42 0.47 .5] %variation of drag force for cd=.04 (N) fd=.5*d*a*v.^2*cd(1) figure(1) plot(v,fd,'color','c') xlabel('velocity') ylabel('drag…

  %row vector a a=[1 2 44 54 58 8]; %column vector b b=[1;2;3;4;5];  

  %row vector a a=[1 2 44 54 58 8]; %column vector b b=[1;2;3;4;5];