Solving the steady and unsteady 2D heat conduction problem using explicit and impicit methods in MATLAB.

AIM – Determine the 10 most sensitive reactions with respect to temperature in stoichiometric methane air combustion using Cantera THEORY - Sensitivity analysis investigates the effect of parameter change on the solution of mathematical models. In chemical kinetics, models are usually based on differential equations…

AIM – To simulate auto ignition in a methane air mixture using Cantera in Python interface. Study the effect of initial temperature and pressure on the ignition delay in the same mixture. Study the effect of temperature on molar concentration of H2O, O2 and OH in the same system at 5 atm. THEORY - Auto ignition occurs…

AIM – To solve a system of coupled Ordinary Differential Equations using Multivariate Newton Raphson Method THEORY – Newton Raphson method for a single variable is used for root finding of non linear equations which are given as – `f(x)=0` It is an iterative method which approximates the value of x by…

AIM – Study the different numerical methods for solving ODEs to understand the stability issues in the solvers. Ordinary Differential Equations (ODEs) – Differential equations are among the most important mathematical tools used in producing models in the physical sciences, biological sciences, and engineering.…

AIM - Develop a compact notation for the given reaction mechanism so that it can be used in a computer program. The reactions in the mechanism are - `CO + O_2 => CO_2 + O ....................R1` `O + H_2O => OH + OH ....................R2` `CO + OH => CO_2 + H ....................R3` `H + O_2 => OH + O ....................R4`…

Aim – Study the effect of air preheating temperature on the Adiabatic Flame temperature in Methane – Air Combustion. Analyzing the effect of preheating temperature on the combustion efficiency assuming constant exit temperature in the furnace (Combustion products – CO2, H2O and N2)   Theory –…

Answer to the questions -  1. Moles method and its calculation approach The moles function defines the number of moles of the object A (comprised of the species O2 and N2). By basic chemistry, we arrive at the number of moles of A - `x_(O_2) = n_(O_2)/(n_(O_2)+n_(N_2)) = n_(O_2)/n_(air)` At STP, the oxidizer air contains…

Aim – Determination of Adiabatic Flame temperature of Methane in a Constant Volume process at different equivalence ratio using Python and comparing it with the results obtained from Cantera. Determination of Flame temperature of different Hydrocarbons (Ethane, Ethene and Acetylene) with heat loss in a constant pressure…

Aim – Volumetric combustion simulation in a 2D combustor model using Ansys Fluent. Studying the impact of water addition in the fuel inlet on pollutants (NO and Soot). Introduction – In this work, a combustion simulation has been performed on a 2D combustor model using Ansys Fluent. Steady state temperature,…

Aim – Three dimensional parametric study on a gate valve using Ansys Fluent. Introduction – Gate valves work by inserting a rectangular gate or wedge into the path of a flowing fluid. They are operated by a threaded stem which connects the actuator (generally a hand wheel or motor) to the stem of the gate.…

Aim – Three dimensional transient CFD study for prediction of collection efficiency of a cyclone separator. Introduction – Cyclone separators or simply cyclones are separation devices (dry scrubbers) that use the principle of inertia to remove particulate matter from flue gases. A high speed rotating (air)…

Aim – Three dimensional CFD analysis of air cooling of a graphics card setup.Introduction – A Conjugate heat transfer (CHT) study has been performed for flow of air over a graphics card setup (consisting of a processor, fins and base plate) in Ansys Fluent.The impact of inlet velocity on the average temperatures…

Aim – Two dimensional CFD analysis of Rayleigh Taylor Instability. Introduction – The transition to turbulence begins usually with the instability of laminar state. This has been an area of interest from the last two centuries. Mathematicians and engineers have many a times offered remarkable insights from…

Aim – Three dimensional CFD analysis of flow through an engine exhaust port. Introduction – Understanding of the gas exchange process in an internal combustion engine (ICE) is essential in order to construct efficient engines. In a Diesel engine, of the total available energy, approximately 30-40% of the energy…

Aim – Three dimensional CFD analysis of flow over an Ahmed body using Ansys 2020 R2. Introduction – The wake flow behind a car is the major contributor to the drag and also poses severe problems during numerical and experimental investigations. The point of separation of the flow determines the size of the…

Aim – Two dimensional CFD analysis of flow over a cylinder using Ansys 2020 R2. Introduction – Flow over a cylinder is one of the most extensively fluid dynamics problems studied, both experimentally and numerically. Intensive investigations are still performed due to the complex dynamics of the wake flow.…

Aim – Three dimensional CFD analysis of turbulent mixing of hot and cold fluid streams in T - junction using Ansys 2020 R2. Introduction – The turbulent mixing problem is one of the most studied CFD problems, both computationally and experimentally. Its applications are widespread in the chemical, petrochemical,…

OBJECTIVES - Write a MATLAB code to take angle as an input and generate blockMeshDict file for simulation through pipe. Use symmetry BC and analyze the results for wedge angles of 10, 25 and 45 degrees. RESULTS -   %MATLAB code to generate blockMeshDict file for simulation of Hagen %Poiseulli flow in openFoam clear…

Objectives -  Simulate the laminar incompressible flow through a pipe using openFoam. Develop a MATLAB code for generating the blockMeshDict file to be used for this simulation. Compare the results of the simulation with the theoretical results using the Hagen Poiseuill's equations Methodology -  The icofoam…

In the FVM, interpolation schemes are used to determine the values of properties at cell faces after discretization. Flux limiter functions are used with a class of discretization schemes called as Total Variation Diminishing schemes to achieve oscillation free solutions in CFD. These schemes were developed to obtain higher…

OBJECTIVE - The aim of this project is to design and analyze incompressible,transient, laminar flow in a 2D backstep channel using openFOAM. METHODOLOGY - 1. Geometry creation and Meshing A similar tutorial to study a steady flow in a backstep channel is already available in tutorials folder of openFoam.  However,…

OBJECTIVES - Developing MATLAB based solvers to simulate quasi 1D supersonic nozzle flows using both the non conservative and conservative forms of the governing equations.  Determine the steady state distributions of the flow field variables (density, temperature, velocity etc.) for both the forms. Determine the…

The general 2D transient heat conduction equation is given as below -  `(delT)/(delt) = alpha((del^2T)/(delx^2) + (del^2T)/(dely^2))` where `alpha`is the thermal diffusivity of the material (Assumed to be copper) and T is the temperature in Kelvin. The above equation needs to be solved numerically using the Jacobi,…

Logic -  The objective has been achieved by reducing the second order DE to two first order DE. In this way, we can utilize the MATLAB solver ode45 to then solve both these equations. The result is then a 2x1 column vector whose 1st column is the Angular displacement while the 2nd column is the angular velocity. An…

The entire report and the MATLAB Codes have been attached using google drive.  Please Refer it. The MATLAB Code and the results are -  %Comparison of the different 4th order schemes for the second derviative of % f(x) = exp(x)*cos(x) % f\'(x) = exp(x)*cos(x) - exp(x)*sin(x) % f\"(x) = -2*exp(x)*sin(x) clear all…