Week 8 - Simulation of a backward facing step in OpenFOAM

The entrance length is the distance that the fluid need to travel through before the full development of the flow. see the following image the entry length for laminar flow is L = 0.06*Re*D the study variable values that I chose Reyonlds number, re = 2100 Diameter, d = 0.004 radius = d/2 = 0.002 theta = 4 Density, rho…

the governing equation for conservative and non conservative forms The initial profile for a, rho, t and v as it is defined in CFD book by John Anderson.  x = linspace(0,3,n); % calculate intial profilesrho = 1-0.3146*x;t = 1-0.2314*x; % t = tempraturev = (0.1+1.09*x).*t.^0.5; % areaa = 1+2.2*(x-1.5).^2 Non conservative…

The structure dimensions are as it is described in the task. First, I divided the study area into five blocks which make it easy to draw the structure using blockMesh. The blocks corners will be numbed as the picture below. the following code explains how this structure could be build in blockMesh vertices( // Front Face…

First defining the boundary conditions, the space nodes size and important study parameters Then define “if else” to switch the method regarding iterative solver identifier, after that a “while” loop to get the error to the required accuracy. The space will be covered by two “for” loop…

I start by deriving the skewed right-sided difference, then skewed left-sided difference and then the central difference (see below the matrix A_right, A_left and A_center). My skewed schemes are fourth-order accurate : Order of approximation = 5 - 2 + 1 = 4 (for central difference) Order of approximation = 6 - 2 = 4 (for…