Week 3.5 - Deriving 4th order approximation of a 2nd order derivative using Taylor Table method

Aim: Write a function that extracts the 14 co-efficient and calculates the enthalpy, entropy and specific heats for all the species in the data file. Calculate the molecular weight of each species and display it in the command window.Plot the Cp, Enthalpy and Entropy for the local temperature range (low…

Aim: Create a Rankine Cycle Simulator using OCTAVE/MATLABCode should calculate the state points of the Rankine Cycle (any type of your choice) based on user inputs. Plot the corresponding T-s and h-s plots for the given set of inputs. Title and axes labels are a must, while legends could be used if necessary.Attach…

Aim:  To write a program that solves the following ODE. This ODE represents the equation of motion of a simple pendulum with damping     In the above equation, g = gravity in m/s2, L = length of the pendulum in m, m = mass of the ball in kg, b=damping coefficient. Write a program in Matlab\Octave that will…

Aim: Write a program in Matlab to simulate the forward kinematics of a 2R Robotic Arm, as shown in the video.      Errors Encountered when executing the code:   1. Error Code:   Output: ·        Explanation:        Typo error, while…

  Aim: Describe Ahmed's body and its importance. Setup the domain and setup assuming body is travelling at 110 kmph Run the setup till convergence and plot for static pressure and coefficient of drag and lift to check for convergence Turn on Broadband- modelling and detail the acoustic field study over Ahmed body.…

Aim: To simulate the motion of blades in a generic turbomachinery component using Steady & Unsteady approach.     Simulation setup: Initially we have three zone mesh files, so first we will read one file and then append the other two file in fluent. File >> read the mesh file >> mesh model will be displayed. Domain >>…

What are the approaches by which you would model moving zones? Two types of motoin involved in fluid flow are as follows Rotational Motion Rectilinear Motion Rotational motion in CFD is modelled by the following approaches. Moving Reference Frame (MRF) Moving Mesh Method (MMM)     Navier stokes equation is solved separately…

Turbulence Model: Reynolds No.: It is the ratio of the inertial force to the viscous flow of the fluid `Re=(rho*V*l)/mu`   The fluid flow is classified into different categories based on the Reynolds no.                                …

Aim: A pitot static tube      a.You are required to do a literature review of pitot-static tube. Also establish a relationship between Total (stagnation) pressure and static pressure     b.You are required to create surfaces and fill it separately for Total pressure port and static pressure…

  Aim: Simulate the flow through converging-diverging nozzle for the following boundary conditions,  At Inlet, Pressure = 1atm and Temperature = 300K. The outlet pressure= 0.027 atm. In the operating conditions, consider the reference pressure to be zero. Plot: The pressure, velocity and temperature variations…

Aim: For the given model, check for the geometrical errors and delete surfaces which are unwanted for surface wrap. After clearing geometry, merge all 3 models and perform surface wrap. Target Length for wrap is 4mm.   Surface wrap: It is a type of surface meshing techniques which is used in the complex geometry having…

Aim: For the given Tesla Cyber Truck model, check the geometrical errors on half portion. Create and assign PID’s, surface mesh the half model and do symmetry to the other side. Target length for the different parts of a model are as follows Body: 3mm Rims and Tyres: 2mm Window: 2mm     Procedure: Open…

Aim: For the given BMW car model, check the geometrical errors on half portion. Create and assign PID’s, surface mesh the half model and do symmetry to the other side. Target length for the different parts of a model are as follows Body: 3mm Lights: 4mm Grills: 1mm Rims and Tyres: 2mm Optical Rear View Mirror: 2mm…

Aim: For the given Turbocharger model, check the geometrical errors to make appropriate volumes. Create and assign PID’s, surface mesh with the given target length as per PID’s. Blade stage-1: 1mm Blade Stage-2: 1mm Impeller: 2mm Shaft rotor: 1mm Turbo casing: 5mm Compressor casing: 5mm Inlet casing cover:…

Aim: To check for the geometrical errors and perform topology clean up accordingly. Set three different target lengths as three different cases (5mm, 3mm, 1mm) and mesh the model.     Procedure: Open ANSA, change to CFD mode then click on the “?” symbol another window is opened in which scroll down…

Aim: To perform the 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 and also plot CO2, H2O, CH4, N2, O2, NOx emissions & soot formation. Addition of water from 5% to 30% by mole…

Aim: To study the flow through the gate valve for different opening position positions and also studying flow coefficient and flow factors for various mass flow rate at the outlet.   Theory: For measuring the valve capcity of control valves, the coefficient of Cv or its Metric equivalent Kv is adopted. The Valve flow…

Aim: To simulate the flow through the cyclone separator and find its efficiency for varied velocity and particle sizes.   Theory: Cyclones are devices that employ a centrifugal force generated by a spinning gas stream to separate particles from the carrier gas. Their simple design, low capital cost, and nearly maintenance-free…

Aim: To Simulate the flow over Graphic card and study Conjugate heat transfer method. Create Heat generation source terms for the processor and to study the temperature around the surrounding components and its effect.   Theory: In order to study air flow around the graphic card, we first need to select the material…

Aim: To simulate the Multiphase simulation and study the Rayleigh Taylor Instability. To also study the relationship between flow simulation and mesh size and Atwood Number   Theory: Rayleigh Taylor Instability (RTI) it arises when the two fluid of different density are subjected to a acceleration due to density gradient…

Aim: To simulate flow through the exhaust manifold wherein 4 inlet port is connected to 1 outlet port. Velocity, Pressure, heat transfer coefficient and Nusselt number for the flow are simulated for the manifold.   Conjugate Heat Transfer: In real world the heat transfer takes from solid and fluid, in solids conduction…

Aim : To simulate 2D flow over Ahmed body to study Coefficient of Drag, Lift, Velocity , Pressure, Streamlines and wake generation around it. Using Ansys Fluent as the simulation tool we will look into the Boundary layer formation around the Ahmed body.   Ahmed body and its Importance: Ahmed body is a bluff body which…

Aim : To simulate 2D flow over cylinder for varying Reynold number and observing the Coefficient of Drag, Lift, Strouhal no., Avg. vortex velocity for steady & transient flow. Further to study the von kormann vortex shedding phenomenon using Ansys Fluent. Theory: When fluid flow over a cylinder wake generation starts…

Aim:  To study T section which gives the desired outlet temperature, whether short or long. which Turbulence model captures the result efficiently for given momentum ratio. Ansys Workbench: Geometry: The geometry of the project is loaded in here and control volume is extracted using space claim. It is used as a inbuilt…

In this project we try to compare the results of Wedge and symmetry boundary condition in openFoam. Assumption of Pipe Flow: We assume the fluid flow to be laminar in nature i.e. Re = 2100 (Reynolds Number). At 30*C the water properties are density=995.6 Kg/m^3, Kinematic viscosity =0.801 X 10^(-6), Dynamic viscosity =…

This project deals with the flow of water in the cylindrical pipe. The flow in the pipe is possible only due to the pressure difference. By using the hagen Poiseuille equation the pressure difference can be found                                …

Finite Difference Method (FDM):                                                The computational domain is divided into node and then PDE ‘s are discretized and converted into algebraic solution…

Aim: Perform a transient simulation over backward facing step geometry to study the effects of different grading factor and plot there result using plot over line command. Code setup: Inside the terminal window Type 'tut' : it is shortcut for directly going to the tutorial directory of open Foam. Type ' cd incompressible/icoFoam/cavity'…

Description:     Here, we try to study quasi 1D supersonic nozzle. Inorder to do that we consider a convergent-divergent(CD) nozzle, which changes its crossection area along the x direction only. Further we assumw that air flows from a very large reservior so velocity is small, pressure and density of air…

Explanation: 2D conduction equation is given by, `(delT)/(delt)=alpha*((del^2T)/(delx^2)+(del^2T)/(dely^2))`   so, by using the given values in the question the main program code is created. Main Program: clear all close all clc %input nx=30; ny=nx; x=linspace(0,1,nx); dx=x(2)-x(1); y=linspace(0,1,ny); dy=y(2)-y(1);…

Explanation: In order to solve the 2nd order differential equation using Taylor table method we need to find no. of node necessary for symmetric and asymmetric stencil. For Symmetric stencil: ( central scheme )         no. of nodes = (order of accuracy) + (order of derivative) – 1    …