Week 7 - Simulation of a 1D Super-sonic nozzle flow simulation using Macormack Method

The Rayleigh–Taylor instability, or RT instability (after Lord Rayleigh and G. I. Taylor), is an instability of an interface between two fluids of different densities which occurs when the lighter fluid is pushing the heavier fluid. Examples…

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 locations of the combustor as shown…

Aim: To perform a parametric study on a gate valve at different valve opening settings.   Theory: Gate valves work by inserting a rectangular gate or wedge into the path of a flowing fluid. They are operated by a threaded stem which connects the actuator (generally a hand wheel or motor) to the stem of the gate.…

l. OBJECTIVE To perform steady-state conjugate heat transfer analysis on Graphics Card To find the effect of different velocities on the temperature   ll. INTRODUCTION 1. Conjugate Heat Transfer (CHT) Conjugate heat transfer is defined as the heat transfer between two domains by exchange of thermal energy. For a system…

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…

Importing the exhaust part geometry into spaceclaim Performing volume extract on the geometry using edge selection option Basic mesh generation This mesh has no conformance or it is a non conformal mesh. So shared topology has to be performed on the geometry. Enabling shared topology Regenerating the baseline mesh Creating…

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: To simulate the flow over a cylinder and explain the phenomenon of Karman vortex street and also to Simulate the flow with the steady and unsteady case and calculate the Strouhal Number for Re= 100.       2. Calculate the coefficient of drag and lift over a cylinder by setting the Reynolds number…

Aim: To check and solve all the geometrical errors, assign relevant PIDs, and perform surface mesh on the cyber truck geometry as well as the wind tunnel surfaces. Working Methodology : The working methodology that is followed to generate mesh on any geometry is : 1. Import the geometry. 2. Check for geometrical errors.…

Aim/Objective : To do surface wrap over assembled parts (engine,transmission and gear box) Target length for Wrap = 3mm Procedure : First the model is loaded on to ANSA Unwanted surface are deleted without damaging the outer shape of the geometries. Single cons are resolved mainly here as it would affect negetively on…

Aim: For the given model, check and solve all geometrical errors on half portion and Assign appropriate PIDs. Perform meshing with the given Target length and element Quality criteria. After meshing the half model, Do symmetry to the other side. Objective: - Target lengths for the different parts of a model are as…

Initial turbocharger model Using delete with PIDregion selection to remove the unwanted surfaces After surface deletion After deletion of extended surface Creating new surface using COONS option with loop selection After surface creation  Projecting the hot points to the cons using multi project Creating small surfaces…

Aim: To generate mesh for the given model and also check for the geometrical errors and perform Topology cleanup accordingly. To set three different target lengths as three different cases and mesh the model.  Target length = 1mm, 3 mm and 5 mm (i.e. 3 cases) Element type = Tria Performing Topo clean up Geometry check…

Introduction Mixing tees—in which two fluid streams with different physical and/or chemical properties mix—are widely used in the petrochemical industry. When there is a large temperature difference between two streams, large metal temperature fluctuations can occur in the region where the two streams meet…

Aim: To describe finite volume methods, differences between finite volume and finite difference methods along with the need for the interpolation schemes in the FVM and flux limiter. METHOD OF computational fluid dynamics: To solve the real-world problems,  three simples schemes for the numerical solution are employed…

OBJECTIVE. To simulate an axis-symmetric laminar flow through the constant cross-sectional of the pipe by applying the symmetry boundary condition. To simulating for the mentioned angles i.e for 10, 25, 45 degrees using both symmetry and wedge boundary conditions and validate them with HP equations.  Hagen-Poiseuille's…

Abstract The flow through a pipe is the case of internal flow. The surface is the duct in the case of air, it is a pipe in the case of water. The flow is driven because of the pressure difference created by the pumps, air handling units, fans, etc. They find their major application refrigeration and air conditioning industries.…

Backward-Facing Step (BFS) flow is one representative model for separation flows, which can be widely seen in aerodynamic flows (airfoil, spoiler, high attack angle process), engine flows, condensers, vehicles (cars, boat), heat transfer systems, and even the flow around buildings, etc. The flow separation after a simple…

Simulation of a 1D Super-sonic nozzle flow simulation using Macormack Method Nozzle A nozzle is a device designed to control the direction or characteristics of a fluid flow (specially to increase velocity) as it exits (or enters) an enclosed chamber or pipe. A nozzle is often a pipe or tube of…

Conductive Heat transfer The process by which heat is transferred from the hotter end to the colder end of object is known as conduction. Heat spontaneously flows from a hotter to a colder body. For example, heat is conducted from the hotplate of an electric stove to the bottom of a saucepan in contact with it. In the…

Taylor Series In mathematics, the Taylor series of a function is an infinite sum of terms that are expressed in terms of the function's derivatives at a single point. For most common functions, the function and the sum of its Taylor series are equal near this point. Taylor series…

RANKINE CYCLE The Rankine cycle is an idealized thermodynamic cycle describing the process by which certain heat engines, such as steam turbines or reciprocating steam engines, allow mechanical work to be extracted from a fluid as it moves between a heat source and heat sink. The…

Main script clear allclose allclc f=fopen('THERMO.dat','r'); %Temperature raange definitionT=linspace(200,3500,200);%Universal gas constantR = 8.314; for i =1:5 fgetl(f);end for k=1:53c = fgetl(f);d = strsplit(c,' ');species_name = d{1};z = length(d);t1 = str2double(d{z-3}); %lowtempt2 = str2double(d{z-2}); %hightempt3…

Main Script   clear all close all clc x=linspace(0,0.8,200); y=linspace(0,0.8,200); [xx yy] = meshgrid(x,y); num_cases=50; 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); end end %Surface Plot figure(1) surfc(xx,yy,-f) xlabel('x value')…

Main Scriptclear allclose allclc%To develop a matlab script to solve the ODE of a motion of a simple%pendulum with damping%ODE is given by%d/dt*(theta2) + b/m*(theta2) + g/l*(sin(theta1))%Defining the input parametersb=0.05;g=9.81;l=1;m=1;w = 3;%Angular displacement at inital position theta = 0;%initial conditiontheta_0=[0;3];%timepointst_span=linspace(0,20,500);%solve…

%Forward kinematics of a robotic arm with movie plot clear all close all clc %inputs %Length m %theta degrees l1=1; l2=0.5; theta1=linspace(0,90,4); theta2=linspace(0,180,4); ct=1; for i=1:length(theta1) THETA1= theta1(i); for j = 1:length(theta2) THETA2 = theta2(j); x= 0; y= 0; %Position of link 1 x1=l1*cosd(THETA1);…