close all clear all clc   %inputs l1=1; l2=0.5; theta1=linspace(0,90,10); theta2=linspace(0,90,10);   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); y2= y1+ l2*sind(THETA2);   %plotting…

clear all close all clc   %inputs b=0.05; g=9.81; l=1; m=1;   %initial condition theta_0=[0;3];   %time points t_span=linspace(0,20,500);   %solve ODE [t,results]= ode45(@(t,theta)ode_func(t,theta,b,g,l,m),t_span,theta_0);   plot(t,results(:,1)) hold on plot(t,results(:,2)) xlabel('time') ylabel('plot')…

clear all close all clc   %defining our search space x=linspace(0,0.6,150); y=linspace(0,0.6,150);   %creating a 2 dimensional mesh [xx yy]=meshgrid(x,y);   %evaluating stalagmite function for i=1:length(xx) for j=1:length(yy) input_vector(1)=xx(i,j); input_vector(2)=yy(i,j); f(i,j)=stalagmite(input_vector);…

close all clear clc   %Inputs R= 8.314;   %Reading the data file f1 = fopen('THERMO.dat','r');   %Reading the first line of the file L1 = fgetl(f1);   %Getting the global minimum, maximum and mid temperatures L2 = fgetl(f1); t = strsplit(L2); global_min = str2double(t{2}); global_mid = str2double(t{3});…

close all clear all clc % displaying rankine cycle simulator and its processess disp(" Rankine cycle simulator") fprintf('\n 1-2 is the isentropic expansion process') fprintf('\n 2-3 is the constant pressure heat rejection process') fprintf('\n 3-4 is the isentropic compression process') fprintf('\n 4-1 is the constant…

close all clear all clc % displaying rankine cycle simulator and its processess disp(" Rankine cycle simulator") fprintf('\n 1-2 is the isentropic expansion process') fprintf('\n 2-3 is the constant pressure heat rejection process') fprintf('\n 3-4 is the isentropic compression process') fprintf('\n 4-1 is the constant…

close all clear clc   %Inputs R= 8.314;   %Reading the data file f1 = fopen('THERMO.dat','r');   %Reading the first line of the file L1 = fgetl(f1);   %Getting the global minimum, maximum and mid temperatures L2 = fgetl(f1); t = strsplit(L2); global_min = str2double(t{2}); global_mid = str2double(t{3});…

AIM: To set up steady-state simulations to compare the mixing effectiveness when hot inlet temperature is 360C & the Cold inlet is at 190C using k-epsilon and k-omega SST model for the following model. Case 1  Short mixing tee with a hot inlet velocity of 3m/s. Momentum ratio of 2, 4.  Case 2  Long…

1) Aim - - Simulate the flow over the cylinder and explain the phenomenon of Karman vortex street.- First,simulate the flow with steady and unsteady case and calculate the Strouhal Number.- Second,calculate the coefficient of drag and lift over a cylinder by setting the Reynolds number to 10,100,1000,10000…

Aim: External flow simulation over an Ahmed body. Objective: The objective of this project is to determine the aerodynamic forces on the Ahmed body such as drag and lift coefficient and to perform the grid independence test. The expected results will include 1. Velocity and pressure contours.  2. The drag coefficient…

Aim: Steady-state CHT analysis on Exhaust port at an inlet velocity of 5m/sec Objective: The objectives will mainly focus on Give a brief description of why and where a CHT analysis is used. Maintain the y+ value according to the turbulence model and justify the results.  Calculate the wall/surface heat transfer coefficient…

PARAMETRIC STUDY ON GATE VALVE INTRODUCTION A gate valve, also called sluice valve, often installed in a pipe fitting, opens by lifting a barrier called gate out of the path of the fluid. Gate valves require very little space along the pipe axis and hardly restrict the flow of fluid when the gate is fully opened. The gate…

Simulating combustion of natural gas  Aim: Perform a combustion simulation on the combustor model using methane as fuel. Objectives: Part I Perform a combustion simulation on the combustor model and plot the variation of the mass fraction of the different species’ in the simulation using line probes at different…

Aim: To perform analysis on cyclone separator and calculate the separation efficiency and pressure drop.   Objective: To write a few words about any four empirical models used to calculate the cyclone separator efficiency.  To perform an analysis on a given cyclone separator model by varying the particle…

RAYLEIGH TAYLOR INSTABILITY Objectives: What are some practical CFD models that have been based on the mathematical analysis of Rayleigh Taylor waves? In your own words, explain how these mathematical models have been adapted for CFD calculations. Perform the Rayleigh Taylor instability simulation for 2 different mesh…

Aim: Steady-state CHT analysis on Exhaust port at an inlet velocity of 5m/sec Objective: The objectives will mainly focus on Give a brief description of why and where a CHT analysis is used. Maintain the y+ value according to the turbulence model and justify the results.  Calculate the wall/surface heat…

CHT ANALYSIS ON A GRAPHICS CARD OBJECTIVE To perform steady-state conjugate heat transfer analysis on Graphics Card To find the effect of different velocities on the temperature  INTRODUCTION Conjugate Heat Transfer (CHT) Conjugate heat transfer is defined as the heat transfer between two domains by exchange of thermal…

Surface meshing on a Pressure valve Given Data To clean up the geometry and mesh with the below criteria Target length = 1mm, 3 mm and 5 mm  Element type = Tria To Build So slit the geometry into 3 cases to show the differentiation of mesh output Case1 Target length-5 mm  Quality criteria min length- 3mm max…

  CFD meshing on turbocharger  Aim: For the given model of a turbocharger, check for the geometrical errors to make appropriate volumes. Create and assign PIDs accordingly. Create surface mesh and use that to create a volumetric mesh.  Objective: 1. Perform surface mesh with the given target lengths…