Week 8 - Simulation of a backward facing step in OpenFOAM

Sensitivity Analysis: Aim: To write a code that takes all the reactions from the GRI-3.0 mechanism and calculates 10 most sensitive reactions with respect to the Temperature and show the reaction strings in the order of greatest sensitivity to the least sensitivity. Theory: Every reaction is represented by a global…

Combustion Efficiency Calculation after Preheating Aim: To calculate combustion efficiency of a furnace with preheating. Theory: A recuperator is a special purpose counterflow energy recovery heat exchanger positioned within the supply and exhaust air streams of an air handling system, or in the exhaust…

Handling Mixtures with Cantera: Aim: Use Quantity class in Cantera, to perform mixtures calculations. Theory: Balanced equation for stoichiometric combustion of methane is as follows: `CH_4 + 2(O_2 + 3.76N_2) -> CO_2 + 2H_2O + 7.52N_2` Composition of Air: Air is assumed to be contain only oxygen and Nitrogen. 21% oxygen…

ODE Stability: Aim: To understand when a Ordinary Differential Equation solved using Explicit method, becomes unstable. And understand its effect on engineering simulations. Description: In mathematics, an ordinary differential equation (ODE) is a differential equation containing one or more functions of one independent…

Auto ignition using Cantera PART I Aim: To simulate the Auto Ignition Time for Methane-air mixture under the following conditions, (simulation time = 10 secs). To plot the variation of Auto Ignition time of Methane with a constant temperature of 1250K and pressure varying from 1 to 5 atm.…

Multivariate Newton Rhapson Solver Aim: To solve set of coupled ordinary differential equations using Implicit Euler Method/Backward Differencing for discretizing, multivariate Newton Rhapson Method for solving the non-linear algebraic equations obtained from discretization and numerical differentiation for calculating…

Compact Notation Derivation for a simple Mechanism: Aim: We will derive the reaction rate ODEs & production of each species for the simple reaction mechanism given below by using the compact notation (indices for reactions and species). First we will find the net reaction rate (`dotq`) for all four reactions using…

1.Adiabatic flame temperature variation with Equivalence ratio (In constant volume-chamber): Assume that methane is contained in a constant volume chamber. Write a program in Python to plot the effect of equivalence ratio on the final adiabatic flame temperature.  Calculate the results with Cantera for varying equivalence…

Simulation of flow over a cylinder and explain the phenomenon of Karman vortex street Aim:    Part-1: To Simulate the flow with the steady and unsteady case and calculate the Strouhal Number for Re= 100.     Part-2: To calculate the coefficient of drag and lift over a cylinder by setting the…

Simulating Cyclone separator with Discrete Phase Modelling Aim: To perform analysis on cyclone separator and calculate the separation efficiency and pressure drop. 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…

GATE VALVE SIMULATION: Objective: In this challenge our objective is to perform parametric study on the gate valve by setting the opening from 10% to 80%. Thereby we will obtain mass flow rates at the outlet for each opening point. Calculate the flow coefficient and flow factor for each opening and plot the graph.…

 Simulating Combustion of Natural Gas Objective: In this challenge, we will be simulating combustion of non-premixed methane-air mixture. Introduction: Natural gas is a naturally occurring hydrocarbon gas mixture consisting of methane and commonly including varying amounts of other higher alkanes, and sometimes…

CHT Analysis on a Graphics card Objective: In this challenge our objective is to perform a steady state simulation of a conjugate heat transfer in a graphics card. We will be running the simulation for default mesh and refined mesh. The air flow over the surface is varied from 1m/s to 5m/s. From this simulation we want…

Rayleigh Taylor Instability Aim: In this our aim is to simulate a setup involving two phases using Volume of Fluid Method. To perform the Rayleigh Taylor instability simulation for 2 different mesh sizes with the base mesh being 0.5 mm. To run one more simulation with water and user-defined material(density = 400 kg/m3,…

CHT Analysis on Exhaust port Aim: To simulate Exhaust port by Conjugate Heat Transfer Analysis. Introduction:  The exhaust port is passage which transport the burned mixture during  produced during power stroke. These burnt gases are at high temperature and therefore be sent into the atmosphere without any issue.…

 External flow simulation over an Ahmed body Aim: Simulation of external flow over Ahmed Body. Introduction:   In this challenge, we will be simulating flow a special body, called Ahmed Body.  Ahmed body was developed by scientist named Ahmed, hence the name. The reason for using this particular…

Objective: To set up steady-state simulation to compare the mixing effectiveness when hot inlet temperature is 360C & the Cold inlet is at 190C, in the both versions of mixing tee (longer and shorter).  Simulations:  Case-1  a) Short mixing tee with a hot inlet velocity of 3m/s and momentum ratio of…

Finite Volume Method(FVM):        Finite Volume Method is a discretisation method, used for approximation of partial differential equations. It is a method for evaluating and representing partial differential equations in the form of algebraic equations. The key step of the control volume is the integration…

Objective:               In this challenge, we are going to simulate an incompressible-laminar-viscous flow through the backward facing step geometry.  This is a transient simulation. Procedure followed: Under this, we will be explaining each and every step followed while doing simulation.…

Simulation of a 1D Super-sonic nozzle flow using Macormack Method:                     In this exercise, we are going to simulate the isentropic flow through a quasi 1D subsonic-supersonic nozzle using Macormack Method. We will be using both the conservation and non-conservation…

Aim: Solving 2D heat conduction equation by Steady state analysis & Transient State Analysis by implicit and explicit method using iterative techniques like Jacobi, Gauss Seidel and SOR method. Inputs: Domain: square slab (1m x 1m) Top Boundary condition= 600 K Bottom Boundary condition = 900 K Left Boundary condition…

4th order approximations of the second-order derivative: We will be using Talor Table Method to formulate the numerical derivative. This method is a convenient and systematic way to deal with, while formulating higher order approximation. Taylor Table is derived from the Talor series and hence the name. It converts all…