Week 4.1- Handling Mixtures with Cantera

Aim: TO Mesh the pressure vale with three different mesh sizes, by doing geometry cleanup and naming the surfaces accordingly .   Methodology: Firstly load the model and  check for single cons 1 and from faces-> new ->option  use the appropriate selection and fille the surfaces by selecting the loop…

ODE : Ordinary differential equation MVNR : multi variate newton raphson method  FDS: forward difference scheme Aim : To find the roots of the  first order 3 -ODE by using the multi variate newton raphson method Methodology: To apply the MVNR , for a system of 3 ODE's , we convert into an implicit form by…

Aim : In this assignment we need to setup the case for the single cylinder SI engine running this at 1800 rpm and finding the below calues from the report.  air flow rate BMEP BSFC In-cylinder pressure In the case setup the speed is set to 1800 rpm, the solution is ran using the default configurations and the report…

Aim : To prove the stability condition of ODE , by comparing the analytical solution with explicit and implicit schemes . Methodology : To prove this we have assumed an ODE     `y' =-100*y-e^t`    where the initial value of this ODE is    `y(0) =0`  that is at time t =0 , y =0 the…

Aim : In this challenge the production rates of the reactions given is solved by using chemical kinetics.   Procedure: To solve this first we identify the number of species and make matrices of reactants and products and formulate a reaction matrix , which will be eventually utilized to formulate the production rates.…

import cantera as ct import numpy as np import matplotlib.pyplot as plt gas_1 =ct.Solution("gri30.yaml") A =ct.Quantity(gas_1) T =np.linspace(298.15,600,40) Tad =[] B =ct.Quantity(gas_1) B.TPX = 298.15,101325,"CH4:1" for i in T: A.TPX = i,101325,"O2:1,N2:7.52" A.moles =2 M=A+B M.equilibrate('HP') M.T Tad.append(M.T) print(Tad)…

for question  4.(1) import cantera as ct gas = ct.Solution("gri30.yaml") A =ct.Quantity(gas) A.TPX =298.15,ct.one_atm,'O2:1,N2:3.76' print("cantera ",A.moles) A.moles =8.875 """ Analytical method , here we are assuming for 1 g of air , mass of air is 1 g and mass per one mole is given by molar mass of the mixture,…

the above program utilizes the cantera library in the python and calculates the Adiabatic flame Temperature by varying equivalence ratios, import cantera as ct import numpy as np import math import sys import csv import matplotlib.pyplot as plt """" RICH mixture CH4+(2/PHI)(O2+3.76N2)-> (4/PHI-3)CO2 +2H2O+ (7.52/PHI)…

The below image shows the various modules of GT Suite  A few have discussed breifly  Engine Performace and Acoustics : We can find the performance parameters of the engines such as Brake Thermal Efficiency ,Brake Power and Indicated Power , Volumetric efficiency ..etc  of the engine . Exhaust After treatment…

Aim: In this simulation we are using combustion model , For a methane fuel and air . The mass fractions of soot and other products  was observed at the outlet by varying amount of water composition in the fuel . Introduction:  Combustion is the chemical reaction which takes place between the fuel in the presence…

Aim : In this simulation  we control the flow through a pipe by a gate valve , and we observe the variation of flow though different opening levels of gate valve disc. The coefficient of flow is calculated and a parametric study is undergone by using the opening level as the input parameter,The simulation was conducted…

Aim:  The main aim of this project is to find the optimum performance of the cyclone separator by using two methods , firstly by changing the size of the anthrasite particles by fixing the velocity ,secondly by varying the velocity of the particle, inlet of the separator by fixing the particle size.The particle entry…

Aim : In this project our aim to find the possible way to reduce the heat which is developed in the electronic equipment, the heat generating component it the processor , since it  generates a high heat it has to be cooled such that we can safe guard  the soldering of the electronic circuits and hence preventing…

Rayleigh Taylor Instabilities Aim: The main aim of this project is to observe the Rayleigh instabilities between two fluids, lying one above the other with fluid with the higher density above. The comparison was made by changing the grid size and changing the heavier fluid with respect to time. Introduction:…

Conjugate Heat Transfer: In this simulation, there is a transfer of heat between solid and fluid mediums and there are two forms are heat transfers viz., Conduction and Convection. The typical applications include heat exchangers. etc Aim: To conduct a CHT analysis on an exhaust port, with  four inlets and one…

Aim: To simulate the velocity variation across the Ahmed body, using a Steady-state solver with K-epsilon turbulence model in fluent. The geometry was initially uploaded in space claim and two enclosures were created as shown below        Since the Ahmed body is symmetric, a symmetric plane is used…

geometry and mesh files     Unsteady vs Steady flow for a Reynolds nuber 100 It is observed that , the steady state required few number of iterations in comparison with unsteady flow, for a unsteady flow we have used the time step size of 0.1 and number of time steps are 1500 to get a converged solution, whereas…

Abstract:  In this simulation a mixing tee case is used and a comparison is made with varying the momentum ratios , varying the tee length , varying the turbulence models and also varying the grid size. And it is observed that the kOmegaSST  turbulence model  gives a better result which is closer to the…

The report is created in latex and the pdf file is uploaded below .

The latex code of the documentation and code used for total pressure , and the p_total plot is added in the latex below  \documentclass[12pt]{report} \usepackage{graphicx} \usepackage{caption} \author{Sai Krishnakanth} \title{Creating a custom IOObject in VsCode} \maketitle \begin{document} \section{Introduction}…

Scalar Transport Equation For solving the scalar transport equation  `(dels)/(delt)+nabla.(us)=0`  we are  using  the icoFoam solver  Copying the icoFoam solver from the solvers to applications directory     2. Moving the solvers to application directory which i created , and editing…

In this challange a comparison was made by using the 2 boundary conditions 1.Symmetry 2.Wedge Symmetry  In the Symmetry boundary condition the front and back side of the sector , is given the symmetry BC , and by varying the angle of the sector we compare wth the analytical and with each other The blockMeshDict file…

                                                       …

finite volume method :Finite volume method  is numerical method of integration , here we divide the domain in finite volumes and calculate the effect of the surfaces of the volume on a central node. We show the numerical descritization of the finite volume method using an ine dimensional  conduction equation …

For creating the geometry we are using 3 blocks , These blocks were created by  editing the system/blockMeshDict file /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ /…

1D Quasi Sub Sonic to Super Sonic Nozzle   Symbols and notation `rho` density `rho_o` reservoir density `rho^'=(rho)/(rho_o)`  dimensionless density V  velocity `V^'=V/(a_o)` Dimension less velocity `a_o` sound velocity T Temperature `T_o` Temperature of the reservoir `T^'=T/T_o` Dimensionless temperature…

Transient case: The 2D Conduction Equation  `(delT)/(delt)=alpha((del^2T)/(delx^2)+(del^2T)/(dely^2))` is a parabolic 2nd order partial differential equation , we can solve this equation by  using explicit and implicit method  First we descritise this equation by using finite difference method .We use forward…

For calculating the second order derivative of the given function at fourth order approximation , we need to  solve a set of finite difference eqautions, But firstly we need to find the number of nodes in  the stencil for each kind of scheme Central Difference Skewed right Skewed Left CENTRAL DIFFERENCE SCHEME…