clear all close all clc x=pi/4; dx=linspace(pi/40,pi/60); % f(x)=exp(x)*cos(x); % exact_derivative_second_order = -2*exp(x)*sin(x); exact_derivative = -2*exp(x)*sin(x); % calculating the coefficients of the equation in forward differencing % scheme % obtaining 6 equations for 6 unknowns from the taylor table, we use matrix…

clear all close all clc b = 0.05; g = 9.81; L = 1; m = 1; % time span t_span = linspace(0,20,500); % initial conditions theta_0 = [0;3]; % [ displacement ; velocity ] [t, results] = ode45(@(t,theta) ode_func(t,theta,b,g,L,m),t_span,theta_0); displacement = results(:,1); velocity = results(:,2); figure('windowstate','maximized','Name','Second…

clear all close all clc %% Determining minima values for the stalactites % Creating an array of 1D xi = linspace(0,1,200); yi = linspace(0,1,150); % Creating a mesh of x and y [X,Y] = meshgrid(xi,yi); % Creating a dummy matrix of length(yi)*length(xi) f =nan(size([length(yi) length(xi)])); % Specifying co-ordinates within…

  clear all close all clc % open THERM data in the read mode t1 = fopen('THERMO.dat','r'); % Reading header, won't be assigning to anything, will clear off at program end fgetl(t1); % Reading next temperature line gl_temp = fgetl(t1); % Splitting the string global_temp_split = strsplit(gl_temp); % Converting to number,segregating…

clear all close all clc diary Rankine_Cycle_Results.txt fprintf('\t\t RANKINE CYCLE SIMULATOR\n\n') fprintf('1-2 : Isentropic Expansion in Turbine\n') fprintf('2-3 : Constant pressure Heat Rejection by the Condensor\n') fprintf('3-4 : Isentropic Compression in Compressor\n') fprintf('4-1 : Constant Pressure Heat Addition…

clear all close all clc %% Define dimensions of the domain L = 1; % Nodes on each axis, only n is enough as uniform mesh is considered n = 50; % Define x and y array x = linspace(0,1,n); y = linspace(0,1,n); % Size of mesh grid dx = L/(n-1); dy = dx; % Define K as mesh grid in constant and does'nt change K = ((2*(dx^2+dy^2))…

Main Script calling non conservative and conservative functions. clear all close all clc %% Defining the common parameters L = 3; n = [61,31]; gamma = 1.4; C = 0.5; nt = 1:1400; %% The subsonic-supersonic isentropic nozzle solution: Non-Conservation Form [drho_dt,x,A,rho,V,T_throat,p,M,rho_throat,mfr,p_throat,M_throat]…

Finite Volume Method: A Literature Review The finite volume method (FVM) is a discretization technique for partial differential equations, especially those that arise from physical conservation laws. FVM uses a volume integral formulation of the problem with a finite partitioning set of volumes to discretize the equations.…

In this challenge we are asked to simulate a flow of fluid through a Backward facing step, in which the fluid is Incompressible, Laminar and viscous in nature. The simulation shall be Transient in nature; thus, we shall have to discretize an unsteady Reynolds scalar transport equation for a fluid flow. Based on the above…

Week 11 - Simulation of Flow through a pipe in OpenFoam In this challenge, we shall simulate an axi-symmetric flow by applying a wedge boundary condition. Problem Statement: The flow is axi-symmetric. Wedge Boundary condition is used. Flow is laminar and flows through a constant cross section pipe. Reynolds number is 2100.…

Week 12 - Validation studies of Symmetry BC vs Wedge BC in OpenFOAM vs Analytical H.P equation In this project, we shall not be performing analysis on the complete circular cross section of a pipe, rather we shall be considering different sizes of the section of circular cross section of pipe for analysis and simulation.…

Aim     To set up steady state simulations in order to compare mixing of hot and cold fluid through outlet of a geometric tee model. Given 2 cases, were in two versions of a mixing tee have been modelled. Need to use k-epsilon and k-omega SST model for case 1 and based on the judgement, use only one method accordingly…

Aim: To study vortex shedding over a cylinder Objective: Simulation of steady and unsteady state Calculate lift and drag coefficient Calculate Strouhal number for Reynolds number 100 in both steady and unsteady case Perform steady state simulation for various cases of Reynolds number. Introduction: Flow over cylinder in…

The aim, Introduction, Theory The Ahmed body, even though it does not replicate the exact shape of a vehicle body used in the current industry, the shape of the rear of the vehicle is imperative in capturing the forces, pressure and velocity in the wake region of any generic vehicle. The airflow around the body captures…

Introduction, aim and objective Conjugate heat transfer (CHT) analysis is used were there is a temperature variation during heating or cooling of a material as a result of interaction between solid and fluid phases. It is a transfer of heat in solids to fluids and vice versa Some of the examples include HVAC, Heat Exchangers,…

OBJECTIVE: The objective of this challenge is to simulate a two-dimensional Rayleigh-Taylor instability problem. INTRODUCTION: In Multiphase flows, the dynamic variables like velocity, pressure, viscosity, density etc are the driving members for movement of fluids under gravitational field. One such example is of Rayleigh-Taylor…

In this simulation, we shall be performing a study of conjugate heat transfer using steady state conditions on a CAD model of a graphics card. The various steps involved shall include performing necessary steps in CFD SpaceClaim, to make the geometry complacent for further analysis, meshing individual components according…

SIMULATING CYCLONE SEPARATOR Aim: To perform an analysis on cyclone separator and calculate the separation efficiency and pressure drop. Objective: To write a few words about any four empirical models used to calculate cyclone separator efficiency. To perform an analysis on a given cyclone separator model by varying particle…

Week 9 – Parametric Study on Gate Valve   AIM: To perform a parametric study on the gate valve simulation by setting the opening from 10% to 80%.   Objective: Obtain the mass flow rates at the outlet for each design point. Calculate the flow coefficient and flow factor for each opening and plot the graph.…

In this simulation, we have performed a combustion simulation on the combustor model and plotted the variation of the mass fraction of the different species in the simulation using line probes at different locations of the combustor. Plots for CO2, H2O, CH4, N2, O2, NOx emissions & Soot formation have been attached.…