Week 8 - BlockMesh Drill down challenge
Objective:
simulation of incompressible-laminar-viscous flow through the backward facing step geometry. by performing a transient simulation using icoFoam solver and velocity profile is ploted at a distance of 0.085 from starting node
Block diagram:
Steps to be follwed in open Foam:
after editing 1)BlockmeshDict 2) controlDict 3)P(pressure) 4) Velocity(V) type command para Foam to get the mesh and contours of domain
BlockMeshDict:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 0.01;
vertices
(
(0 0 0)//point 0
(8 0 0)//point 1
(8 0.5 0)//point 2
(0 0.5 0)//point 3
(8 1 0)//point 4
(0 1 0)//point 5
(20 1 0)//point 6
(20 0.5 0)//point 7
(20 0 0)//point 8
(20 -1 0)//point 9
(8 -1 0)//point 10
(0 0 0.1)//point 11
(8 0 0.1)//point 12
(8 0.5 0.1)//point 13
(0 0.5 0.1)//point 14
(8 1 0.1)//point 15
(0 1 0.1)//point 16
(20 1 0.1)//point 17
(20 0.5 0.1)//point 18
(20 0 0.1)//point 19
(20 -1 0.1)//point 20
(8 -1 0.1)//point 21
);
blocks
(
hex (0 1 2 3 11 12 13 14) (50 10 1) simpleGrading (1 1 1)
hex (3 2 4 5 14 13 15 16) (50 10 1) simpleGrading (1 1 1)
hex (2 7 6 4 13 18 17 15) (50 10 1) simpleGrading (1 1 1)
hex (1 8 7 2 12 19 18 13) (50 10 1) simpleGrading (1 1 1)
hex (10 9 8 1 21 20 19 12) (50 10 1) simpleGrading (1 1 1)
);
edges
(
);
boundary
(
inlet
{
type patch;
faces
(
(0 11 14 3)
(3 14 16 5)
);
}
outlet
{
type patch;
faces
(
(9 8 19 20)
(8 7 18 19)
(7 6 17 18)
);
}
frontAndBack
{
type empty;
faces
(
(0 3 2 1)
(3 5 4 2)
(10 1 8 9)
(1 2 7 8)
(2 4 6 7)
(11 12 13 14)
(14 13 15 16)
(21 20 19 12)
(12 19 18 13)
(13 18 17 15)
);
}
Noslipwalls
{
type wall;
faces
(
(5 16 15 4)
(4 15 17 6)
(0 1 12 11)
(10 21 12 1)
(10 9 20 21)
);
}
);
mergePatchPairs
(
);
// ************************************************************************* //2) ControlDict:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application icoFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 12;
deltaT 1e-2;
writeControl timeStep;
writeInterval 20;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable true;
// ************************************************************************* //Pressure:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
frontAndBack
{
type empty;
}
Noslipwalls
{
type zeroGradient;
}
}
// ************************************************************************* //Given problem nx=50,ny=10,nz=1 taken instead of nx=200,ny=20,nz=1
Simple meshgrid scheme (GF=1)
hex (0 1 2 3 11 12 13 14) (50 10 1) simpleGrading (1 1 1)
hex (3 2 4 5 14 13 15 16) (50 10 1) simpleGrading (1 1 1)
hex (2 7 6 4 13 18 17 15) (50 10 1) simpleGrading (1 1 1)
hex (1 8 7 2 12 19 18 13) (50 10 1) simpleGrading (1 1 1)
hex (10 9 8 1 21 20 19 12) (50 10 1) simpleGrading (1 1 1)
Mesh block diagram
pressure-variation diagram:
Velocity variation diagram:
Velocity profile:
comment:
velocity profile is at distance of 0.085 fnear step from inlet
Simple meshgrid scheme (GF=0.2)
hex (0 1 2 3 11 12 13 14) (50 10 1) simpleGrading (0.2 1.2 1)
hex (3 2 4 5 14 13 15 16) (50 10 1) simpleGrading (0.2 0.2 1)
hex (2 7 6 4 13 18 17 15) (50 10 1) simpleGrading (1.2 0.2 1)
hex (1 8 7 2 12 19 18 13) (50 10 1) simpleGrading (1.2 1.2 1)
hex (10 9 8 1 21 20 19 12) (50 10 1) simpleGrading (1.2 1.2 1)
Mesh block diagram
pressure-variation diagram:
Velocity variation diagram:
Velocity profile:
comment:
velocity profile is at distance of 0.085 fnear step from inlet
Simple meshgrid scheme (GF=0.5)
hex (0 1 2 3 11 12 13 14) (50 10 1) simpleGrading (0.5 1.5 1)
hex (3 2 4 5 14 13 15 16) (50 10 1) simpleGrading (0.5 0.5 1)
hex (2 7 6 4 13 18 17 15) (50 10 1) simpleGrading (1.5 0.5 1)
hex (1 8 7 2 12 19 18 13) (50 10 1) simpleGrading (1.5 1.5 1)
hex (10 9 8 1 21 20 19 12) (50 10 1) simpleGrading (1.5 1.5 1)
Simple meshgrid scheme (GF=0.5)
pressure-variation diagram:
Velocity variation diagram:
Velocity profile:
Simple meshgrid scheme (GF=0.8)
hex (0 1 2 3 11 12 13 14) (50 10 1) simpleGrading (0.8 1.8 1)
hex (3 2 4 5 14 13 15 16) (50 10 1) simpleGrading (0.8 0.8 1)
hex (2 7 6 4 13 18 17 15) (50 10 1) simpleGrading (1.8 0.8 1)
hex (1 8 7 2 12 19 18 13) (50 10 1) simpleGrading (1.8 1.8 1)
hex (10 9 8 1 21 20 19 12) (50 10 1) simpleGrading (1.8 1.8 1)
Simple meshgrid scheme (GF=0.8)
Pressure variation diagram:
Velocity variation diagram:
Velocity profile:
comment:
velocity profile is at distance of 0.085 fnear step from inlet
Gradingfactor doenot effect the profiles of velocity and pressure