Week 5.1 - Mid term project - Solving the steady and unsteady 2D heat conduction problem

DescriptionConsider a given convergent-divergent nozzle, the flow at the inlet is fed from a large gas reservoir where the gas is essentially stationary. The reservoir pressure and temperature are Po and To, respectively. The subsonic Mach number at the inlet is very small. Thus, the pressure and temperature at the…

Project2-EDA on Vehicle Insurance Customer Data PFA  

Project 1 - English Dictionary App & Library Book Management System PFA File  

Project 4 In this challenge, we will be running a tutorial from $FOAM_TUTORIALS/compressible/sonicFoam/shockTube and replacing the controlDict, fvSchemes and fvSolutions with that from $FOAM_TUTORIALS/basic/scalarTransportFoam/pitzDaily Instead of keeping temperatures on one side as high and the other side as low, we will…

Creating a custom IOObject in Visual Studio Code In this challenge, we will create a custom IOobject in visual studio code to compute the total pressure flowfield for an incompressible case, by copying pitzDaily from ~incompressible/simpleFoam and pasting it in the run directory. The total pressure ptotal is…

Adding a temperature dictionary to the incompressible icoFoam case to solve the energy equation A pdf file is attached that contains the code in LATEX. In this project, we are adding a temperature dictionary to the incompressible icoFoam solver to solve for the energy equation. The following changes are done: First, copy…

Creating your own OpenFOAM solver In this challenge, the objective is to create a custom OpenFOAM software by modifying some details in the codes and in the files. At first, the solver is copied from ~applications/solvers/incompressible/icoFoam and the folder is renamed as scalarFoam. Once done, rename the code as scalarFoam.C…

Creating a test case in OpenFOAM In this challenge, we will find the value assigned to kinematic viscosity stored in a dictionary file transportProperties and add a variable to it and print the values of the variables. Procedure   #include "fvCFD.H" //using namespace Foam; int main(int argc, char *argv[]) { //Foam::argList::noBanner();…

FVM Fundamentals 101 Derivation In this challenge, we shall look at the derivation of second order accurate first order scheme, followed by applying FVM formulation to the equation and describing about TVD schemes in detail. The taylor series expansion of `f_(i+1)` about the point i is given b` `f_(i+1)=f_i + f'_iDeltax…

                                                      Week 4 Challenge : Compiling C++ libraries In this challenge, we will copy a specified folder from the library into our run folder…

Senstivity Analysis Assignment Aim In this experiment, we will be working on a source code provided in the video by doing minor changes like: -> writing a code that takes all the reactions from the GRI mechanism and calculating 10 most sensitive reactions -> The sensitivity parameter created should be with respect…

                                                                Auto ignition using Cantera Aim In this challenge, we will work on auto ignition using cantera by…

Aim To solve the following set of equations as mentioned below by using multivariate newton raphson method and assuming appropriate initial conditions. Theory The Newton-Raphson method (also known as Newton's method) is a way to quickly find a good approximation for the root of a real-valued function `f(x)=0`. It uses…

Literature review: ODE Stability Aim In this challenge, we will be working on a literature survey to answer the following questions alongside writing a code that will help us in answering these as follows: When does a ODE become unstable? How does stability affect engineering simulations? Can stability condition be derived…

Compact Notation Derivation for a simple Mechanism Aim In this challenge, we will derive the reaction rate ODEs & production of each species for the simple reaction mechanism given below and also use the compact notation that was discussed in the Week 5 lecture videos (Part 1&2) and write down the desired matrices. …

Combustion Efficiency Calculation after Preheating AIM In this assignment, we shall calculate the effect of preheating on the adiabatic flame temperature by changing the pre-heating temperature(PHT) from 298K to 600K. Also, we will cross-verify whether preheating helps with fuel saving by plotting the effect of pre-heat…

Adiabatic flame temperature calculation using python In this assignment, we are going to write a code to calculate AFT(Adiabatic Flame Temperature) for a given reaction mechanism. Also, we will try to observe the variation of the AFT with varying equivalent ratios, with varying heat losses , with number if carbon atoms,etc.…

Emission Characterization on a CAT3410 engine In this project, we will simulate two geometries of engine(open-v and omega) types and obtain useful performance parameters and animations showing different plots. Also, we shall characterise the emissions of NOx,soot and UHC by suitable models like Hiroy's soot model,etc.…

FINAL TEST PFI QUESTIONS Q1. What is the Compression ratio for the engine? The compression ratio (CR) is defined as the ratio of the volume of the cylinder and its head space (including the pre-combustion chamber, if present) when the piston is at the bottom of its stroke to the volume of the head space when…

Simulation of full hydro setup of a Port Fuel Injection(PFI) engine In this project, we shall continue the previous project by setting up a full-hydro case setup in converge and post process the results file and visualise the meshing process. PFI-a look PFI engines use a three-way catalytic converter, exhaust sensors and…

Surface preparation and Boundary Flagging of a Port Fuel Injection engine(PFI) In this project, we will prepare the surface of a PFI engine by cleaning the given CAD file and marking the boundaries followed by which we shall simulate a no hydro simulation in converge and postprocessing in paraview. PFI-a look PFI engines…

Literature review - RANS derivation and analysis In this challenge, we shall derive reynolds averaged navier stokes equation and come up with the expression for reynolds stress. Also we shall see what is reynolds stress, and what is turbulent viscosity and how it is different from molecular viscosity. Theory All flows…

Shock tube simulation project In this project, we will be working on simulating flow through a shock tube by constructing a geometry and setting up the case in Converge CFD software and  Plot the pressure and temperature history in the entire domain and explain the result Plot the cell count as a function of time…

Conjugate Heat Transfer Simulation In this simulation, we will setup a Conjugate heat transfer(CHT) analysis of flow through a solid pipe with a thickness and Setup a flow simulation through a pipe Inlet Reynolds number should be 7,000 Add Y+ to your post.in file Start with an initial base grid Run grid dependence test…

Prandtl Meyer Shock problem In this simulation, we will understand what a shock is (expansion fan) and understand its significance in the flow properties alterations. Also, we will Do a literature search on what BC's are typically used for shock flow problems  We will understand what is a shock wave. In your own words,…

Transient simulation of flow in a throttle body The aim of this challenge is to perform transient simulation of flow over a throttle body and to  Set up and run the transient simulation for flow over throttle body. Post process the results and show pressure and velocity contours. Show the mesh (i.e surface with edges)…

                                      Steady state simulation of flow in a throttle body The aim of this challenge is to perform steady state simulation of flow over a throttle body and to  Set up and run the steady state…

                                                             Simulating Combustion of Natural Gas Table of contents Aim Theory Procedure Results and discussions Aim…

                                                                   Parametric study on Gate valve Aim To perform a parametric study on the given model of gate valve by varying the opening from 10% to 80% and performing the following calculations: To obtain the mass flow rates at the outlet for each design point. To calculate…

                                          Simulating Cyclone separator with Discrete Phase Modelling Table Of Contents 1. Aim 2. Theory 3. Procedure 4. Results 5. Animation videos 6. Conclusions Aim: 1. To perform analysis…

                                               Conjugate Heat transfer analysis over a graphics card Table of contents ->Aim ->CHT(Conjugate Heat Transsfer) ->Procedure ->Fluent Settings ->Results…

Challenge: In this challenge, we are going to conduct a CHT(aka conjugate heat transfer) analysis in an exhaust port in Ansys Fluent. Here, we will try to explain what is Y+ and what are its use, calculate the surface heat transfer coefficient on the internal solid surface and show the velocity. CHT(Conjugate Heat Transfer):…

Challenge In this challenge, Rayleigh taylor instability simulation is performed for 2 different mesh sizes and an extra simulation with a user defined material incorporated in fluent. The following tasks are carried out in this challenge: 1. To compare the results by showing the animations of all the cases. 2. To observe…

Challenge In this challenge, an ahmed body is simulated in FLUENT at a velocity of 25m/s with the default air properties. As the body is perfectly symmetric, we can run the simulation by considering the only half body. This is the best practice where we can save on the number of cells and get the results faster as well.  The…

FLOW OVER A CYLINDER   Challenge In this challenge, flow over a cylinder is simulated in fluent and certain characteristics like vortex sheet, boundary layer flow separation,Reynolds number,etc are calculated. Also we will simulate the flow with the steady and unsteady case and calculate the Strouhal Number for Re= 100.…

MIXING TEE JOINT FLOW SIMULATION IN ANSYS FLUENT In this challenge, a mixing tee is simulated with two different cases-a longer pipe and a shorter pipe. A CFD analysis is carried out over the mixing tee subjected to two cases which are enlisted below: Case 1: Short mixing tee with a hot inlet velocity of 3 m/s and with…

AIM: To continue the previous challenge by generalising the automation code for generating the blockmesh file for various angles(`theta=10^o ,25^o,45^o`) under various symmetry and wedge conditions subjected to initial velocity and pressure conditions CONTINUATON: A symmetric condition exists when a section of a pipe extends…

Challenge In this challenge, an axi-symmetric flow is simulated by applying the wedge boundary condition. The flow is laminar and flows through a constant cross-section pipe at an inlet Reynold's number of 2100. The results are then validated with the Hagen- Poiseuille's equations for the fully developed flow. Generally,…

AIM To simulate an incompressible laminar viscous flow through the backward facing step geometry and to perform a transient simulation. The test case should also be performed with 2 cases of grading factor and the results are enlisted. THEORY AND PROCEDURE: A laminar flow is one in which the reynolds number for that external…

AIM To derive the conservative and non-conservative forms of the governing equations and to simulate the isentropic flow through a convergent-divergent nozzle by writing a code and plotting sufficient graphs. THEORY SETUP There are three steps needed to be followed here  -> The governing equations will be couched…

AIM To understand the steady and unsteady heat conduction over a domain by considering iterative solvers like Jacobian, Gauss-Seidel and Successive_Over_Relaxation and to write a MATLAB program to do the same. THEORY & BOUNDARY CONDITIONS The general,transient 2D equation of conduction is `(delT)/(delt)=alpha((del^2T)/(delx^2)+(del^2T)/(dely^2))`…

AIM To derive the fourth-order approximations of the second-order derivative by following central difference, skewed right side difference, and skewed left side difference schemes and to write a program to determine the error in each case.   DERIVATION We know that the Taylor's series expansion of a function `f(x)`…