Week 7: Shock tube simulation project

Aim: Simulation of Steady and unsteady flow over cylinder body. Objective: Calculate the drag and lift coefficient over a cylinder by setting the Reynolds number to 10,100,1000,10000 & 100000 using a steady solver. Discuss the effect of Reynolds number on the coefficient of drag. Calculate the strouhal number for the…

Aim: Simulation of flow-through mixing tee using Ansys fluent. Objective: To determine the changes in flow properties for different lengths of geometry, and momentum ratio. Perform mesh independent study for any geometry. Use different turbulence models, and determine which is suitable for accurate results. Introduction:…

Aim: Simulating closed-cycle analysis for diesel engines with different bowl geometries.   Objective: Simulate combustion in an isolated chamber for two different piston bowl geometries. Compare the results from both the pistons and find out which among them is more suitable.   Introduction: CI engine: In compression-ignition…

Questions based on PFI Engine:   Q1. What is the Compression ratio for the engine? It is the ratio of total volume to the clearance volume of the cylinder. Rc = Total volume/Clearance volume = (Vs + Vc) / Vc     = 5.73*10-4 / 5.704*10-5     = 10.04:1     Q2. Modify the compression ratio…

AIM: FULL HYDRO case set up of port fuel injection for IC engine. OBJECTIVE: Simulate the SI engine using Converge CFD software and discuss the results. Calculate engine performance parameters such that we can make design changes. INTRODUCTION: Spray modelling: Converge can model both gaseous and liquid sprays. Converge…

Aim: Surface preparation and Boundary Flagging of IC engine geometry and running no hydro simulation. Objective: Performing a no hydro simulation for IC engine. Introduction: No hydro simulation is one in which we will not be solving any governing equation but will generate a mesh for the whole body.  This simulation…

Aim: Deriving a Reynolds averaged Navier stokes equation. Objective: Derive equation using Reynolds decomposition. Explain the end results of derivation. Introduction: Governing equation for the boundary layer: We can solve these three equation by direct numerical substitution this method is called DNS. But there is a…

Aim: Shock tube simulation using Converge CFD package and post processing. Objective: Setup a transient shock tube simulation in Converge. Plot the pressure and temperature history in the entire domain. Plotting the cell count as function of time. Introduction: Shock tube: A simple shock tube is a tube, rectangular…

Aim: Conjugate heat transfer simulation using Converge Cfd and post processing the results. Objectives: Setup a flow simulation through a pipe with inlet Reynolds number should be 7,000 for air. Add supercycling parameter asnd also Y+ in our post.in file. Simulate the results for different base grids for grid dependence…

Aim: To simulate prandlt meyer shock problem using converge CFD package, and post processing. Objective: Set the case for 2 different drichlet boundary conditions(specified value) for velocity at inlet. 100 m/s at 286.1k  at which speed of sound is approximately 360m/s. 680 m/s at 286.1k at which speed of sound is…

Aim : Transient state simulation of from over throttle body using Converge Cfd. Objectives: Set a case for transient floe in converge Cfd. Throttle rotation rate: At 0 : 0 degree 0 to 2ms : 25 degree 2ms to 4ms : 25 degree 4ms to 8ms : 0 degree Post process in paraview and provide detailed explanation for case with contours,…

Aim : Simulation of flow over throttle body using converge Cfd and post processing using paraview. Objectives : Set up and run the steady state simulation for flow over throttle body. Post process the results and show pressure and velocity contours. Show the mesh (i.e surface with edges) Show the plots for pressure, velocity,…

Aim: Simulating flow over backward facing step. Objective: Create a backward facing step geometry in converge. Set a case and simulate for three different base mesh sizes: dx = dy = dz =2e-3m dx = dy = dz =1.5e-3m dx = dy = dz =1.0e-3m  Introduction: Many industrial application involves flow seperation and reattachment,…

Aim: To set a channel flow in Converge and post-processing the results, also plot for velocity, cell count, and mass flow rate at inlet and outlet. Provide an animation video for velocity distribution using ParaView. Objective: Running the case for three different base mesh sizes 3 different base mesh sizes are case1: dx…

Aim : To find 10 most sensitive reactions from the set of reaction for the combustion of fuel. Objective: Write a code which calculate sensitivity of all the reaction and generates 10 most sensitive reactions. Plot the reactions with there sensitivity, change the equivalence ratio and discuss the effect of it on sensitivity…

Aim:To determine ignition delay for a given fuel at different temperature and pressure. Objective: Determine ignition delay for varying pressure from 1 to 5atm and temperature contant at 1250K. Determine ignition delay for varying temperature from 950 to 1450K and fixed pressure of 5atm. Determine the mole fraction of…

Aim:Solving the set of ordinary differential equation using multivariate newton raphson method. Objective: Solving the set of ordinary differential equation using Implicit Euler Method/forward Differencing.Run the simulation for 10 minutes Plot the number of iteration and values for variable y1,y2, an y3. Introduction:…

Aim: Write a program to check explicit and implicit schemes, and check the stability of ODE. Objective: To determine which scheme is stable. Plot for unstability condition. Define how can we derive stability for all problems. Introduction: When we are deriving finite difference approximation of a exact derivative we are…

Aim: Combustion efficiency calculation for recuperator at different preheating temperatures. Objectives: To determine the effect of preheating on the adiabatic flame temperature? Change pre-heating from 298 to 600K Assume that the furnace exit temperature is 1700K (For this case assume that the combustion products are…

Aim:To perform mixture calculations. Process: We will consider simple equation of fuel reacting with air producing water and carbon dioxide CH4 + 2(O2+3.76N2) = H2O+2CO2+7.52N2 To calculate mole fraction we have formula: X(o2) = n(O2) / n(total) where, n(total) = n(O2)+n(N2) n(total) = 2+7.52 n(total) = 9.52 X(O2) = 2/9.52…

Aim:To derive the reaction rate ODEs of a simple reaction mechanism. Reaction Mechanisms: CO + O2 <-> CO2 + O O + H2o <-> OH + OH CO + OH <-> CO2 + H H + O2 <-> OH + O There are four reactions given and the total number of species in the reactions are 7. To calculate a concentration of particular…

Aim: Adiabatic flame temperature calculation. Introduction: Adiabatic flame temperature: The final temperature that is obtained after the fuel is completely burned is called the adiabatic flame temperature. To calculate the adiabatic flame temperature of any process, the first step is to balance the combustion equation.…

Aim:Parsing the engine_data.out file and writing a code for automated plotting. Objective: Your script should take column numbers as the input and plot the respective columns as separate images. Each file should be saved by the name of the column. The plot labels should be extracted from the file. If, I request for a plot…

Aim: Breaking ice with air cushion vehicle. Objective: Use the newton raphson method to find out the value of pressure for h = 0.6 ft. Find the optimal relaxation factor for this problem with the help of a suitable plot. Tabulate the results of p for h = [0.6,1.2,1.8,2.4,3,3.6,4.2] assuming a suitable relaxation factor.…

Aim: Fitting a curve to given data and estimating goodness of fit. Objective: To write a program in python which uses a given data file and fit a curve. Trying different polynomials, split-wise method to calculate the parameters for perfect fit. Check goodness of fit using R2 method. Introduction: Curve fitting is…

Aim:Solving ordinary differential equation for simulation of simple pendulum. Objectives:To code for motion of simple pendulum and determine the variation of displacement and velocity with respect to time.To add animation for realistic data interpretation.To vary the inputs of the program and find out a change in…

Aim:  Write a code that can solve a otto cycle and plot of it Objectives: To write a basic code to solve Air standard cycle. To find out an efficiency of Air standard cycle and with trial and error method optimize engine efficiency. To plot a results obtained to visualize actual scenario. Introduction: An Otto…

  Finite Volume Method Finite volume method (FVM), like FEM, is based on an unstructured (e.g. triangular) mesh. Therefore, it is suitable for irregular and complex geometries. FVM has another advantage over FEM for fluid mechanic problems. So far, the numerical methods that we presented have been based on PDEs. By…

Aim: To simulate incompressible-laminar-viscous flow through the backward facing step geometry. Objective: Performing transient simulation of incompressible-laminar-viscous flow through the backward facing step geometry. Write a code in openfoam to perform above simulation. Explaining the effect of change in…

Aim:To simulate isentropic flow through quasi 1-dimensional subsonic_supersonic nozzle. Objective: Derive a equation for continuity, momentum, and energy for both conservative and non-conservative form. Write a matlab code for all the above derived equations. structure a matlab code according to MacCormack's technique.…

Aim: Deriving 2D heat conduction equation for steady, unsteady, implicit, and explicit conditions, Using different iterative techniques such as jacobi, gauss siedel, and successive over relaxation.   Introduction: To use navier stokes equation, we should rather divide a equation in several parts, and learn to solve…

Aim: Deriving 4th order approximations of the second-order derivative. using central difference, skewed right-sided difference, skewed left-sided difference. Second Order derivative of fourth order approximation,the method we will be using here is taylor’s table method Second order derivative is given as: Ә2Y/Әx2 …

Aim: Reading nasa thermodynamic data and calculating Cp,enthalpy,entropy and molecular weight of species Objective: To read the data from thermodynamic file and Select required data for calculation. To plot a calculated value for Cp,enthalpy,entropy and molecular weight of species. To create a folder in current directory…

Aim: Writing a program in matlab using genetic algorithm for global optimization of stalagmite function Objective: To generate a code in matlab for stalagmite function. To plot a stalagmite function and find maxima. To optimize a local maxima to global maxima using genetic algorithm in matlab. Plot a stalagmite function…

Aim:Solving ordinary differential equation for simple pendulum and adding animation for visualization of pendulum motion . Objectives: To gather required formulae's,and terminology to write a program on simple pendulum. To code for motion of simple pendulum and determine various plots. To add animation for realistic data…