FVM: The finite volume method (FVM) is a method for representing and evaluating partial differential equations in the form of algebraic equations. In the finite volume method, volume integrals in a partial differential equation that contain a divergence term are converted to surface integrals, using the divergence…

STEADY STATE: Discretization of the equation- For Successive Over-relaxation method, the following formula is used- Here, h^2,S(i,j) is the source term. But, as we are not considering the source term in our problem, this term is neglected.   %Steady State 2D conduction clear all close all clc for a=1:3 n=11; %number…

OBJECTIVE: This is a continuation of the previous project where we studied the Hagen-Poiseuille\'s equation, simulated laminar flow of an incompressible fluid through wedge-shaped pipe section and compared the analytical and observed values. In this part of the project, we will use the symmetry boundary conditions instead…

The objective of this challenge is to study the Macormack technique and set up a simulation case for the quasi 1-D supersonic flow through a nozzle.  Both, the conservative and non-conservative forms of the governing equations were set up. The characteristics of both these forms were studied. The case set-up and initial…

Derivation of fourth order approximations for second derivative using Taylor Table method                          Octave Plot for the CDS, Skewed RHS & LHS Explanation: As can be seen from the plot, the Central Difference(CD) scheme gives the least error.…

In this project, we simulate the laminar flow of an incompressible fluid (water) through a circular pipe. The objectives of this project are: Study the fluid flow through pipe based on the Hagen-Poiseuille equation. Create a wedge-shaped geometry in OpenFOAM and simulate the fluid flow. Write the blockMeshDict file using…

Objective: The objective of this challenge is to create a desired geometry and mesh in OpenFOAM and to run the simulation by applying boundary conditions. The grid dependency was also found out by using different grading factors of 0.2, 0.5 and 0.8. Procedure: The fluid used for this problem was incompressible. The icoFoam…

AIM: To simulate airflow around an Ahmed Body using STARCCM+ and compare the results with the given experimental data. PROCEDURE: The Ahmed Body was created using the in-built 3D CAD modeller in STARCCM+. The dimensions of the body were provided with two different slant angles of 25 and 35 degrees. The simulation was run…

To begin with the simulation set-up, we need to import the CAD and then close the geometry by filling up the holes at the windshield and two windows. This is done using the 'Fill Holes' option in the Surface Repair window.                            …

This mid-term project was a part of the course - Introduction to CFD using MATLAB and OpenFOAM   STEADY STATE: Discretization of the equation- For Successive Over-relaxation method, the following formula is used- Here, h^2,S(i,j) is the source term. But, as we are not considering the source term in our problem, this…

This mid-term project was done as a part of the course- Introduction to CFD using MATLAB and OpenFOAM   Steady State vs Unsteady State comparison: The plots obtained from the previous challenge are compared below. The time taken by MATLAB to complete the simulation was calculated using the tic-toc function. This time…

This mid-term project was done as a part of the course- Introduction to CFD using MATLAB and OpenFOAM   The Courant-Friedrichs-Lewy (CFL) Condition: The CFL condition is the condition for stability and to check the convergence criteria of the numerical solution. It is defined as-'The full numerical domain of dependence…

This project is a part of the course- Introduction to CFD using MATLAB and OpenFOAM   This is a continuation of the previous theoretical project where we studied the Hagen-Poiseuille's equation, simulated laminar flow of an incompressible fluid through wedge-shaped pipe section and compared the analytical and observed…

This project was done as a part of the course- External Aerodynamics using STARCCM+   A mesh is a representation of a larger geometric domain by smaller discrete cells. In CFD, meshes are used to compute solutions of partial differential equations. A mesh partitions space into elements over which the equations…

Aim: In this project, we aim to study the turbulence modelling by simulating a flow over a backward step with varying Reynolds numbers in STARCCM+. The concept of Y+ was also studied.   Calculations: The simulation was done for air at a temperature of 25 degree Celcius. Reynolds Number (Re) is given as: Re = (ρ…

This project was done as a part of the course - External Aerodynamics Simulations using STARCCM+   1. BCs used in External Aerodynamics: Boundary Conditions are the set of constraints to boundary value problems in computational fluid dynamics. These boundary conditions include inlet boundary conditions, outlet boundary…