Week 3: Flow over a backward facing step

Aim: To set up and simulate the Transient simulation of flow over the throttle body. Objective:  Setting up Transient state simulation for flow over a throttle body. To generate the mesh that gives more accurate results. Post-processing the obtained results using ParaView and show pressure and velocity contour. …

Aim: The main objective of this project is to set up the case to perform a simulation of steady-state flow over a Throttle valve when it is fully open and post-processing the result using ParaView.   Introduction: The report describes the steady-state flow simulation over a Throttle valve using converge CFD software,…

Aim:  To set up the backward-facing step case in converge environment and post-process the result in the Paraview Introduction: The flow is through the backward-facing step. In this project converge studio is used to set up the profile, the solver settings, and the initial conditions. The flow is driven by the change…

Aim:  To set up the channel flow case in converge environment and post-process the result in the Paraview Introduction: Channel flow is the internal flow over the walls which changes the hydrodynamic structure of the flow. In this project converge studio is used to set up the profile, the solver settings, and the…

Aim: To write a code in MATLAB to extract thermodynamic data from the file and calculate the cp, H and S for the species provided in the file and specified species.  Plot the thermodynamic properties with respect to the temperature in the respective folder which should be generated automatically.   Objective:…

  STALAGMITE FUNCTION: `f(x,y) = -f_(1,x)*f_(1,y)*f_(2,x)*f_(2,y)` where, `f_(1,x) = [sin(5.1*pi*x) + 0.5]^6` `f_(1,y) = [sin(5.1*pi*y) + 0.5]^6` `f_(2,x) = exp[-4*ln(2)*((x-0.0667)^2/0.64)]` `f_(2,y) = exp[-4*ln(2)*((y-0.0667)^2/0.64)]` Negative sign is introduced to find the maximum value   GENETIC ALGORITHM:…

FUNCTION: function [dy] = myODE_simple_pendulum(t,y,b,g,m,L) dy = zeros(2,1); % creating the matrix %Converting second order ODE into first order ODE dy(1) = y(2); %y' = y2 dy(2) = ((-b/m)*y(2)) - ((g/L)*(sin(y(1)))); % y" = y2' end CODE: clear all close all clc t = [0,20]; %Simulation time y0 =[0, 3]; %Intial conditions…

Code: clear all close all clc %Inputs l1 = 1; %Length of one arm l2 = 0.5; %Length of second arm theta1 = linspace(0,90,10); %Creating points for theta1 theta2 = linspace(0,90,10); %Creating points for theta2 ct = 1; %Variable for frames for i = 1:length(theta1) %For loop theta1_1 = theta1(i); for j = 1:length(theta2)…

Aim: Simulate an incompressible laminar viscous flow through the backwards facing step geometry Perform transient simulation Explain the simulation setup of the problem in OpenFoam Perform two case studies one without using grading factor and the other with the grading factor of 0.2 Objectives: Copy a similar OpenFoam…

Aim: Differentiate FVM and FDM and describe the need for interpolation schemes and flux limiters in FVM. FVM: Finite Volume Method is a discretization method to evaluate the partial differential equation in the form of algebraic equations. It approximates through the control volumes and calculates the flux of the parameters…

Objective: To simulate the isentropic flow through a quasi 1D subsonic-supersonic nozzle. To derive in both the conservative and non-conservative forms of the governing equations by solving using MacCormark's technique. Aim: To plot the steady-state distribution of primitive variables inside the nozzle To plot the time-wise…

Objective: To solve the 2D heat conduction equation in a steady and unsteady state with different solving techniques i.e., Jacobi, Gauss-seidel and Successive over relaxation with implicit and explicit equations. Boundary Conditions: Top boundary condition = 600 K Bottom boundary condition = 900 K Right boundary condition…

function [y,am] = expr(x) y = exp(x)*cos(x); am = exp(x)* (-2*sin(x)); end function error_c = central(x,dx,P,Q) %N = number of Nodes %P = order of accuracy %Q = order of derivative N = P+Q-1 ; %Number of nodes k = (N-1)/2; %constant A_c = zeros(N); %Defining matrix B_c = zeros(N,1); %Defining matrix B_c(Q+1,1) = 1; %Equating…