Week 10 - Simulating Combustion of Natural Gas.

Objective: To calculate Adiabatic Flame Temperature for combustion with variation in equivalence ratio, type of fuel, with heat loss using the Newton-Raphson method and Cantera. Theory: Adiabatic flame temperature(AFT):          The final temperature of the product, When the combustion takes place…

Objective:              To simulate the combustion in Ansys fluent with the parametric study. Theory: Combustion is a chemical process during which heat energy is generated.  During combustion, the Hydro-carbon and Oxygen combined when a spark is introduced chemical reaction takes…

Objective:                 To simulate the Gate valve with different opening positions using parametric study in Ansys fluent. Theory: Flow coefficient: `C_v=Q√((SG)/(∇P))` where,        Q is the rate of flow (expressed in US gallons per minute),…

 Objective:                    To simulate cyclone separator with Discrete Phase Modelling in Ansys fluent. Theory: Empirical models used to calculate the cyclone separator efficiency: 1.Iozia and Leith Model Iozia and Leith logistic model is a modified version of…

Objective:            To simulate the steady-state conjugate heat transfer analysis on the model of a  graphics card. Theory: Conjugate heat transfer: The term Conjugate heat transfer is used to describe that when there is a heat transfer between solid to fluid and vice versa. The conjugate…

Objective:                To simulate the Multiphase flow of rayleigh Taylor instability in Ansys fluent. Theory: Multiphase flows take place in nuclear geothermal power plants, food processing industries, and many others. One of the basic examples is Rayleigh Taylor Instability.…

Objective:   To simulate the Exhaust manifold in Ansys fluent for analyzing the heat transfer and the flow of fluid inside. Conjugate heat transfer: The term Conjugate heat transfer is used to describe that when there is a heat transfer between solid to fluid and vice versa. The conjugate heat transfer (CHT)…

Objective:           External flow simulation over an Ahmed body in Ansys fluent. Theory: Ahmed's body and its importance The Ahmed body is a simplified vehicle model or a bluff model body with basic aerodynamic properties of a vehicle and it was developed for investigating the influence of…

Objective: Simulate the flow over a cylinder and explain the phenomenon of Karman vortex street.    Theory: Karman vortex street: It is a repeated pattern of swirl, caused due to vortex shedding. Vortex shedding is an oscillating flow that takes place when a fluid such as air or water flows past a bluff…

Case 1  Short mixing tee with a hot inlet velocity of 3m/s. Momentum ratio of 2, 4.  Length = 212.85mm Case 2  Long mixing tee with a hot inlet velocity of 3m/s. Momentum ratio of 2, 4. Length = 216mm Hot air temperature is 36 degree Celsius Cold air temperature is 19 degree Celsius V1 = 3m/s,  V2 =…

  dt based on Courant number: Courant number C = U*(dt/dx) Where, U - velocity = 1 m/s assumed dt - time steps to be find dx = x/no. of cells along x axis dx = 0.1/200 = 0.0005 courant number should be less than one (C < 1) thus, using the above equation dt < 5e-5 .   Grading factor (f): f = size of the…

                          Symmetry vs Wedge vs HP equation Objective: To write a Matlab program that takes an angle as input and generates a blockMesh file for the given angle of 10,25,45. To compare the results obtained from solving the Hagen-Poiseuille…

Objective: To write a program in MATLAB that can generate the computational mesh automatically for any wedge angle and grading scheme. To determine, whether the velocity profile matches between the computational method and analytical method. To show that the velocity profile is fully developed after the entry length. To…

    clear all close all clc % 2D HEAT CONDUCTION % (1,2,3) steady state % 4,5,6,7 for unsteady / transient state % (1,2,3) or (5,6,7) for (jacobi,gs,SOR) % 4 explicit % 5,6,7 implicit Z = 7; %input l = 1; %length b = 1; %breadth %number of nodes along x & y axis nx = 10; ny = 10; % Discretization x = linspace(0,l,nx);…

clear all close all clc disp(\'RANKINE CYCLE\'); disp(\'1-2 : isentropic expansion\') disp(\'2-3 : constant pressure heat rejection\') disp(\'3-4 : isentropic compression\') disp(\'4-1 : constant pressure heat addition\') %input parameters t1 = 400; p1 = 30; p2 = 0.05; %---pt 1--- h1 = XSteam(\'h_pt\',p1,t1); s1 = XSteam(\'s_pt\',p1,t1);…

*NOT ABLE TO UPLOAD IN HTML LAST 10 LINES ARE NOT PASTING clear all close all clc file = fopen(\'THERMO.dat\',\'r\'); line1 = fgetl(file); line2 = fgetl(file); split_temp = findstr(line2, \' \'); t_gmin = str2num(line2(split_temp(1):split_temp(4))); t_gmid = str2num(line2(split_temp(5):split_temp(6))); t_gmax = str2num(line2(split_temp(7):end));…

error: without the hold off in line 33 we can\'t get the animation clearly https://www.youtube.com/watch?v=JbXtB-WIv38 clear all close all clc %input Parameters b = 0.2; g = 9.81; l = 1; m = 1; x0 = 0; y0 = 0; % initial condition theta_0 = [1 0]; %time points t_span = linspace(0,20,900); % solve ODE [t, results] = ode45(@(t,theta)…