Skill-Lync Launch Pad – Your Gateway to a Core Engineering Job by 2025! Only 1̶0̶0̶ +50 Seats Available.

00D 12H 59M 42S

Menu

Executive Programs

Workshops

Projects

Blogs

Careers

Placements

Student Reviews

For Business

More

Academic Training

Informative Articles

Find Jobs

We are Hiring!

All Courses

Choose a category

Mechanical

Electrical

Civil

Computer Science

Electronics

Offline Program

All Courses

All Courses

logo

CHOOSE A CATEGORY

Mechanical

Electrical

Civil

Computer Science

Electronics

Offline Program

Top Job Leading Courses

Automotive

CFD

FEA

Design

MBD

Med Tech

Courses by Software

Design

Solver

Automation

Vehicle Dynamics

CFD Solver

Preprocessor

Courses by Semester

First Year

Second Year

Third Year

Fourth Year

Courses by Domain

Automotive

CFD

Design

FEA

Tool-focused Courses

Design

Solver

Automation

Preprocessor

CFD Solver

Vehicle Dynamics

Machine learning

Machine Learning and AI

POPULAR COURSES

coursePost Graduate Program in Hybrid Electric Vehicle Design and Analysis
coursePost Graduate Program in Computational Fluid Dynamics
coursePost Graduate Program in CAD
coursePost Graduate Program in CAE
coursePost Graduate Program in Manufacturing Design
coursePost Graduate Program in Computational Design and Pre-processing
coursePost Graduate Program in Complete Passenger Car Design & Product Development
Executive Programs
Workshops
For Business

Success Stories

Placements

Student Reviews

More

Projects

Blogs

Academic Training

Find Jobs

Informative Articles

We're Hiring!

phone+91 9342691281Log in
  1. Home/
  2. Yogessvaran T/
  3. Week 12 - Validation studies of Symmetry BC vs Wedge BC in OpenFOAM vs Analytical H.P equation

Week 12 - Validation studies of Symmetry BC vs Wedge BC in OpenFOAM vs Analytical H.P equation

Hagen-Poiseuille's law: In nonideal fluid dynamics, the Hagen–Poiseuille equation, also known as the Hagen–Poiseuille law, is a physical law that gives the pressure drop in an incompressible and Newtonian fluid in laminar flow flowing through a long cylindrical pipe of constant cross-section. The…

  • CFD
  • HTML
  • MATLAB

  • Yogessvaran T

    updated on 11 Sep 2022

Hagen-Poiseuille's law:

In nonideal fluid dynamics, the Hagen–Poiseuille equation, also known as the Hagen–Poiseuille law, is a physical law that

gives the pressure drop in an incompressible and Newtonian fluid in laminar flow flowing through a long cylindrical pipe of

constant cross-section. The pressure drop is given by:

ΔP=32μLVD2

where,

ΔP

 is the pressure drop.

 μ

 is the dynamic viscosity.

  L is the length of the pipe.

 V is the velocity.

 

 

At the entrance, the flow is said to be in the hydrodynamic entrance region where the velocity profile gradually takes the

shape of a parabola from a straight line. As soon as the parabolic velocity profile becomes steady the flow is said to be in its

fully developed state. Normally Hagen Poiseuille's Equation is applicable for a fully developed flow which helps us to calculate

the pressure loss.

Assumptions made:

Geometry: Axisymmetric.

Reynolds number: 2100.

Working fluid: Water.

Diameter of the pipe: 0.02 m.

Dynamic viscosity(μ)=0.89x10^-3 Pas

Kinematic viscosity(ν)=0.8927x10^-6 m^2/s

theta=10, 25, 45

 

EQUATIONS USED.

Average velocity

Vavg=ReμρD

 m/s 

Maximum velocity

Vmax=2⋅Vavg

 m/s

Pressure drop

ΔP=32μLVD2

 pa

kinematic pressure drop=ΔPρ

 

 

PROCEDURE.

Will start with copying the cavity file from the tutorial to the FOAM_RUN folder.

In the run folder only we are going to simulate and whatever we want to do the changes in the blockMeshDict file and

pressure and velocity files and controlDict files we will be doing here.

creating the block in the OpenFoam here we are creating the wedge after making the required changes in the blockMeshDict

file to check the geometry run the blockMesh command.

After creating the wedge assign the boundary conditions, and also sitting the simulation end time and time-step in such a

way that the courant number should be less than 0.5 and also in the transportProperties enter the kinematic viscosity.

Run the simulation using the icoFoam solver after calculating the time needed for the fluid to travel from inlet to outlet

After the simulation, use paraFoam to view the flow inside the wedge and use plot over line command to plot the velocity and

shear stress changes at a point after the hydrodynamical entry region so that we analyse the fully developed region of the

fluid. 

Use Matlab to parse the blockMeshDict so that the same file can be used for any wedge values.

 

MATLAB CODE TO GENERATE blockMeshDict FILE for wedge boundary condition.

clear all
close all
clc

%  inputs
% assumptions made the fluid is assumed as water
Re=2100; % given reynolds number
rho=997; % density of the water 
D=0.02; % diiameter of the pipe
r=D/2; % radius of the pipe
mu=0.89*10^-3; % dynamic viscosity of the water at 25 degree celcius
v=0.88927*10^-6; % kinematic viscosity of the water at 25 degree celcius
theta=input('Enter the value of the wedge angle= ')
L_h=0.06*Re*D; % hydrodynamic lenght 

L=L_h+0.5;
% Analytical solutions 
V_avg=(Re*mu)/(rho*D);  % average velocity
V_max=2*V_avg;   % maximum velocity
P_drop=(32*mu*L*V_avg)/D^2;  % pressure drop
P_kinematic=(P_drop)/rho;  % kinematic pressure drop

% vertices 
v0=[0 0 0]
v1=[0 r*cosd(theta/2) r*sind(theta/2)]
v2=[0 r*cosd(theta/2) -r*sind(theta/2)]
v3=[L 0 0]
v4=[L r*cosd(theta/2) r*sind(theta/2)]
v5=[L r*cosd(theta/2) -r*sind(theta/2)]


% PRINTING THE blockMeshDict FILE

H1='/*--------------------------------*- C++ -*----------------------------------*\';
H2='  =========                 |';
H3='  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox';
H4='  \\    /   O peration     | Website:  https://openfoam.org';
H5='    \\  /    A nd           | Version:  8';
H6='     \\/     M anipulation  |';
H7='\*---------------------------------------------------------------------------*/';
H8='FoamFile';
H9='{'
H10='    version     2.0;';
H11='    format      ascii;';
H12='    class       dictionary;';
H13='    object      blockMeshDict;';
H14='}'
H15='// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //';

% opening the file
f1=fopen('blockMeshDict_parsedNEW.txt','wt');

fprintf(f1,'%s \n',H1);
fprintf(f1,'%s \n',H2);
fprintf(f1,'%s \n',H3);
fprintf(f1,'%s \n',H4);
fprintf(f1,'%s \n',H5);
fprintf(f1,'%s \n',H6);
fprintf(f1,'%s \n',H7);
fprintf(f1,'%s \n',H8);
fprintf(f1,'%s \n',H9);
fprintf(f1,'%s \n',H10);
fprintf(f1,'%s \n',H11);
fprintf(f1,'%s \n',H12);
fprintf(f1,'%s \n',H13);
fprintf(f1,'%s \n',H14);
fprintf(f1,'%s \n',H15);
fprintf(f1,'\n');

fprintf(f1,'%s \n','convertToMeters 1;');

fprintf(f1,'%s \n','vertices');
fprintf(f1,'%s \n','(');
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v0(1),v0(2),v0(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v1(1),v1(2),v1(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v2(1),v2(2),v2(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v3(1),v3(2),v3(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v4(1),v4(2),v4(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v5(1),v5(2),v5(3));
  
fprintf(f1,'%s \n',');');

fprintf(f1,'%s \n','blocks');
fprintf(f1,'%s \n','(');
fprintf(f1,'%s %s \n',blanks(5),'hex (0 2 1 0 3 5 4 3) (20 1 200) simpleGrading (1 1 1)');
fprintf(f1,'%s \n',');');

fprintf(f1,'%s \n','edges');
fprintf(f1,'%s \n','(');
fprintf(f1,'arc 1 2 (%f %f %f) \n',0,r,0);
fprintf(f1,'arc 4 5 (%f %f %f) \n',L,r,0);
fprintf(f1,'%s \n',');');

fprintf(f1,'%s \n','boundary');
fprintf(f1,'( \n');
fprintf(f1,'%s %s \n',blanks(5),'inlet');
fprintf(f1,'%s %s \n',blanks(5),'{');
fprintf(f1,'%s %s \n',blanks(10),'type patch;');
fprintf(f1,'           %s \n','faces');
fprintf(f1,'          %s \n','(');
fprintf(f1,'           (%d %d %d %d)\n',0, 1, 2, 0);
fprintf(f1,'          );\n');
fprintf(f1,'      }\n');
fprintf(f1,'      %s\n','outlet');
fprintf(f1,'      {\n');
fprintf(f1,'        %s\n','type patch;');
fprintf(f1,'         %s\n','faces');
fprintf(f1,'         (\n');
fprintf(f1,'            (%d %d %d %d)\n',3, 4, 5, 3);
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,'     %s\n','wedge_front');
fprintf(f1,'     {\n');
fprintf(f1,'         %s\n','type wedge;');
fprintf(f1,'         %s\n','faces');
fprintf(f1,'         (\n');
fprintf(f1,'             (%d %d %d %d)\n',2, 5, 3, 0);
fprintf(f1,'\n');            
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,'     %s\n','wedge_Back');
fprintf(f1,'     {\n');
fprintf(f1,'%s\n','         type wedge;');
fprintf(f1,'%s\n','         faces');
fprintf(f1,'         (\n');
fprintf(f1,'\n');            
fprintf(f1,'             (%d %d %d %d)\n',0, 3, 4, 1);
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,'     %s\n','top');
fprintf(f1,'     {\n');
fprintf(f1,'         %s\n','type wall;');
fprintf(f1,'         %s\n','faces');
fprintf(f1,'         (\n');
fprintf(f1,'             (%d %d %d %d)\n',1, 4, 5, 2);
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,'     %s\n','axis');
fprintf(f1,'     {\n');
fprintf(f1,'         %s\n','type empty;');
fprintf(f1,'         %s\n','faces');
fprintf(f1,'         (\n');
fprintf(f1,'             (%d %d %d %d)\n',0, 3, 3, 0);
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,' );\n');
fprintf(f1,'\n'); 
fprintf(f1,' %s\n','mergePatchPairs');
fprintf(f1,'(\n');
fprintf(f1,' );\n');
fprintf(f1,'\n'); 
fprintf(f1,' // ************************************************************************* //\n');


fclose(f1);

MATLAB CODE TO GENERATE blockMeshDict FILE for symmetry boundary condition.

clear all
close all
clc

%  inputs
% assumptions made the fluid is assumed as water
Re=2100; % given reynolds number
rho=997; % density of the water 
D=0.02; % diiameter of the pipe
r=D/2; % radius of the pipe
mu=0.89*10^-3; % dynamic viscosity of the water at 25 degree celcius
v=0.88927*10^-6; % kinematic viscosity of the water at 25 degree celcius
theta=input('Enter the value of the wedge angle= ')
L_h=0.06*Re*D; % hydrodynamic lenght 

L=L_h+0.5;
% Analytical solutions 
V_avg=(Re*mu)/(rho*D);  % average velocity
V_max=2*V_avg;   % maximum velocity
P_drop=(32*mu*L*V_avg)/D^2;  % pressure drop
P_kinematic=(P_drop)/rho;  % kinematic pressure drop

% vertices 
v0=[0 0 0]
v1=[0 r*cosd(theta/2) r*sind(theta/2)]
v2=[0 r*cosd(theta/2) -r*sind(theta/2)]
v3=[L 0 0]
v4=[L r*cosd(theta/2) r*sind(theta/2)]
v5=[L r*cosd(theta/2) -r*sind(theta/2)]


% PRINTING THE blockMeshDict FILE

H1='/*--------------------------------*- C++ -*----------------------------------*\';
H2='  =========                 |';
H3='  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox';
H4='  \\    /   O peration     | Website:  https://openfoam.org';
H5='    \\  /    A nd           | Version:  8';
H6='     \\/     M anipulation  |';
H7='\*---------------------------------------------------------------------------*/';
H8='FoamFile';
H9='{'
H10='    version     2.0;';
H11='    format      ascii;';
H12='    class       dictionary;';
H13='    object      blockMeshDict;';
H14='}'
H15='// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //';

% opening the file
f1=fopen('blockMeshDict_parsedNEW.txt','wt');

fprintf(f1,'%s \n',H1);
fprintf(f1,'%s \n',H2);
fprintf(f1,'%s \n',H3);
fprintf(f1,'%s \n',H4);
fprintf(f1,'%s \n',H5);
fprintf(f1,'%s \n',H6);
fprintf(f1,'%s \n',H7);
fprintf(f1,'%s \n',H8);
fprintf(f1,'%s \n',H9);
fprintf(f1,'%s \n',H10);
fprintf(f1,'%s \n',H11);
fprintf(f1,'%s \n',H12);
fprintf(f1,'%s \n',H13);
fprintf(f1,'%s \n',H14);
fprintf(f1,'%s \n',H15);
fprintf(f1,'\n');

fprintf(f1,'%s \n','convertToMeters 1;');

fprintf(f1,'%s \n','vertices');
fprintf(f1,'%s \n','(');
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v0(1),v0(2),v0(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v1(1),v1(2),v1(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v2(1),v2(2),v2(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v3(1),v3(2),v3(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v4(1),v4(2),v4(3));
fprintf(f1,'%s (%f %f %f) \n',blanks(5),v5(1),v5(2),v5(3));
  
fprintf(f1,'%s \n',');');

fprintf(f1,'%s \n','blocks');
fprintf(f1,'%s \n','(');
fprintf(f1,'%s %s \n',blanks(5),'hex (0 2 1 0 3 5 4 3) (40 1 500) simpleGrading (1 1 1)');
fprintf(f1,'%s \n',');');

fprintf(f1,'%s \n','edges');
fprintf(f1,'%s \n','(');
fprintf(f1,'arc 1 2 (%f %f %f) \n',0,r,0);
fprintf(f1,'arc 4 5 (%f %f %f) \n',L,r,0);
fprintf(f1,'%s \n',');');

fprintf(f1,'%s \n','boundary');
fprintf(f1,'( \n');
fprintf(f1,'%s %s \n',blanks(5),'inlet');
fprintf(f1,'%s %s \n',blanks(5),'{');
fprintf(f1,'%s %s \n',blanks(10),'type patch;');
fprintf(f1,'           %s \n','faces');
fprintf(f1,'          %s \n','(');
fprintf(f1,'           (%d %d %d %d)\n',0, 1, 2, 0);
fprintf(f1,'          );\n');
fprintf(f1,'      }\n');
fprintf(f1,'      %s\n','outlet');
fprintf(f1,'      {\n');
fprintf(f1,'        %s\n','type patch;');
fprintf(f1,'         %s\n','faces');
fprintf(f1,'         (\n');
fprintf(f1,'            (%d %d %d %d)\n',3, 4, 5, 3);
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,'     %s\n','wedge_front');
fprintf(f1,'     {\n');
fprintf(f1,'         %s\n','type symmetry;');
fprintf(f1,'         %s\n','faces');
fprintf(f1,'         (\n');
fprintf(f1,'             (%d %d %d %d)\n',2, 5, 3, 0);
fprintf(f1,'\n');            
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,'     %s\n','wedge_back');
fprintf(f1,'     {\n');
fprintf(f1,'%s\n','         type symmetry;');
fprintf(f1,'%s\n','         faces');
fprintf(f1,'         (\n');
fprintf(f1,'\n');            
fprintf(f1,'             (%d %d %d %d)\n',0, 3, 4, 1);
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,'     %s\n','top');
fprintf(f1,'     {\n');
fprintf(f1,'         %s\n','type wall;');
fprintf(f1,'         %s\n','faces');
fprintf(f1,'         (\n');
fprintf(f1,'             (%d %d %d %d)\n',1, 4, 5, 2);
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,'     %s\n','axis');
fprintf(f1,'     {\n');
fprintf(f1,'         %s\n','type empty;');
fprintf(f1,'         %s\n','faces');
fprintf(f1,'         (\n');
fprintf(f1,'             (%d %d %d %d)\n',0, 3, 3, 0);
fprintf(f1,'          );\n');
fprintf(f1,'     }\n');
fprintf(f1,' );\n');
fprintf(f1,'\n'); 
fprintf(f1,' %s\n','mergePatchPairs');
fprintf(f1,'(\n');
fprintf(f1,' );\n');
fprintf(f1,'\n'); 
fprintf(f1,' // ************************************************************************* //\n');


fclose(f1);

 

For symmetry boundary conditions 

For angle 25 degrees the blockMeshDict file 

/*--------------------------------*- 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 1; 
vertices 
( 
      (0.000000 0.000000 0.000000) 
      (0.000000 0.009763 0.002164) 
      (0.000000 0.009763 -0.002164) 
      (3.020000 0.000000 0.000000) 
      (3.020000 0.009763 0.002164) 
      (3.020000 0.009763 -0.002164) 
); 
blocks 
( 
      hex (0 2 1 0 3 5 4 3) (40 1 500) simpleGrading (1 1 1) 
); 
edges 
( 
arc 1 2 (0.000000 0.010000 0.000000) 
arc 4 5 (3.020000 0.010000 0.000000) 
); 
boundary 
( 
      inlet 
      { 
           type patch; 
           faces 
          ( 
           (0 1 2 0)
          );
      }
      outlet
      {
        type patch;
         faces
         (
            (3 4 5 3)
          );
     }
     wedge_front
     {
         type symmetry;
         faces
         (
             (2 5 3 0)

          );
     }
     wedge_back
     {
         type symmetry;
         faces
         (

             (0 3 4 1)
          );
     }
     top
     {
         type wall;
         faces
         (
             (1 4 5 2)
          );
     }
     axis
     {
         type empty;
         faces
         (
             (0 3 3 0)
          );
     }
 );

 mergePatchPairs
(
 );

 // ************************************************************************* //

 controlDict file

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     icoFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         30;

deltaT          0.005;

writeControl    timeStep;

writeInterval   20;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;


// ************************************************************************* //

Pressure file

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     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;
    }

    wedge_front
    {
        type            symmetry;
    }

    wedge_back
    {
        type            symmetry;
    }

    top
    {
        type            zeroGradient;
    }

    axis
    {
        type            empty;
    }
}

// ************************************************************************* //

Velocity file

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform (0.09373 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    top
    {
        type            noSlip;
    }

    wedge_front
    {
        type            symmetry;
    }

    wedge_back
    {
        type            symmetry;
    }

    axis
    {
        type            empty;
    }
}

// ************************************************************************* //

Results:

Blockmesh

Velocity profile

Pressure profile

 

 

 

For symmetry boundary conditions 

For angle 10 degrees the blockMeshDict file 

/*--------------------------------*- 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 1; 
vertices 
( 
      (0.000000 0.000000 0.000000) 
      (0.000000 0.009962 0.000872) 
      (0.000000 0.009962 -0.000872) 
      (3.020000 0.000000 0.000000) 
      (3.020000 0.009962 0.000872) 
      (3.020000 0.009962 -0.000872) 
); 
blocks 
( 
      hex (0 2 1 0 3 5 4 3) (40 1 500) simpleGrading (1 1 1) 
); 
edges 
( 
arc 1 2 (0.000000 0.010000 0.000000) 
arc 4 5 (3.020000 0.010000 0.000000) 
); 
boundary 
( 
      inlet 
      { 
           type patch; 
           faces 
          ( 
           (0 1 2 0)
          );
      }
      outlet
      {
        type patch;
         faces
         (
            (3 4 5 3)
          );
     }
     wedge_front
     {
         type symmetry;
         faces
         (
             (2 5 3 0)

          );
     }
     wedge_back
     {
         type symmetry;
         faces
         (

             (0 3 4 1)
          );
     }
     top
     {
         type wall;
         faces
         (
             (1 4 5 2)
          );
     }
     axis
     {
         type empty;
         faces
         (
             (0 3 3 0)
          );
     }
 );

 mergePatchPairs
(
 );

 // ************************************************************************* //

Results:

Blockmesh

Velocity profile

Pressure profile

 

For symmetry boundary conditions 

For angle 45 degrees the blockMeshDict file 

/*--------------------------------*- 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 1; 
vertices 
( 
      (0.000000 0.000000 0.000000) 
      (0.000000 0.009239 0.003827) 
      (0.000000 0.009239 -0.003827) 
      (3.020000 0.000000 0.000000) 
      (3.020000 0.009239 0.003827) 
      (3.020000 0.009239 -0.003827) 
); 
blocks 
( 
      hex (0 2 1 0 3 5 4 3) (40 1 500) simpleGrading (1 1 1) 
); 
edges 
( 
arc 1 2 (0.000000 0.010000 0.000000) 
arc 4 5 (3.020000 0.010000 0.000000) 
); 
boundary 
( 
      inlet 
      { 
           type patch; 
           faces 
          ( 
           (0 1 2 0)
          );
      }
      outlet
      {
        type patch;
         faces
         (
            (3 4 5 3)
          );
     }
     wedge_front
     {
         type symmetry;
         faces
         (
             (2 5 3 0)

          );
     }
     wedge_back
     {
         type symmetry;
         faces
         (

             (0 3 4 1)
          );
     }
     top
     {
         type wall;
         faces
         (
             (1 4 5 2)
          );
     }
     axis
     {
         type empty;
         faces
         (
             (0 3 3 0)
          );
     }
 );

 mergePatchPairs
(
 );

 // ************************************************************************* //

Results:

Blockmesh

Velocity profile

Pressure profile

 

 

For wedge boundary conditions.

For the wedge angle 10 degrees the blockMeshDict file.

/*--------------------------------*- 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 1; 
vertices 
( 
      (0.000000 0.000000 0.000000) 
      (0.000000 0.009962 0.000872) 
      (0.000000 0.009962 -0.000872) 
      (3.020000 0.000000 0.000000) 
      (3.020000 0.009962 0.000872) 
      (3.020000 0.009962 -0.000872) 
); 
blocks 
( 
      hex (0 2 1 0 3 5 4 3) (40 1 500) simpleGrading (1 1 1) 
); 
edges 
( 
arc 1 2 (0.000000 0.010000 0.000000) 
arc 4 5 (3.020000 0.010000 0.000000) 
); 
boundary 
( 
      inlet 
      { 
           type patch; 
           faces 
          ( 
           (0 1 2 0)
          );
      }
      outlet
      {
        type patch;
         faces
         (
            (3 4 5 3)
          );
     }
     wedge_front
     {
         type wedge;
         faces
         (
             (2 5 3 0)

          );
     }
     wedge_back
     {
         type wedge;
         faces
         (

             (0 3 4 1)
          );
     }
     top
     {
         type wall;
         faces
         (
             (1 4 5 2)
          );
     }
     axis
     {
         type empty;
         faces
         (
             (0 3 3 0)
          );
     }
 );

 mergePatchPairs
(
 );

 // ************************************************************************* //

 

controlDict file

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     icoFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         10;

deltaT          0.005;

writeControl    timeStep;

writeInterval   20;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;


// ************************************************************************* //

 

Pressure file

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     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;
    }

    wedge_front
    {
        type            wedge;
    }

    wedge_back
    {
        type            wedge;
    }

    top
    {
        type            zeroGradient;
    }

    axis
    {
        type            empty;
    }
}

// ************************************************************************* //

 

Velocity file

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform (0.09373 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    top
    {
        type            noSlip;
    }

    wedge_front
    {
        type            wedge;
    }

    wedge_back
    {
        type            wedge;
    }

    axis
    {
        type            empty;
    }
}

// ************************************************************************* //

Results:

Blockmesh

Pressure profile

Velocity porfile

For wedge boundary conditions.

For the angle 25 degrees the blockMeshDict file.

/*--------------------------------*- 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 1; 
vertices 
( 
      (0.000000 0.000000 0.000000) 
      (0.000000 0.009763 0.002164) 
      (0.000000 0.009763 -0.002164) 
      (3.020000 0.000000 0.000000) 
      (3.020000 0.009763 0.002164) 
      (3.020000 0.009763 -0.002164) 
); 
blocks 
( 
      hex (0 2 1 0 3 5 4 3) (40 1 500) simpleGrading (1 1 1) 
); 
edges 
( 
arc 1 2 (0.000000 0.010000 0.000000) 
arc 4 5 (3.020000 0.010000 0.000000) 
); 
boundary 
( 
      inlet 
      { 
           type patch; 
           faces 
          ( 
           (0 1 2 0)
          );
      }
      outlet
      {
        type patch;
         faces
         (
            (3 4 5 3)
          );
     }
     wedge_front
     {
         type wedge;
         faces
         (
             (2 5 3 0)

          );
     }
     wedge_back
     {
         type wedge;
         faces
         (

             (0 3 4 1)
          );
     }
     top
     {
         type wall;
         faces
         (
             (1 4 5 2)
          );
     }
     axis
     {
         type empty;
         faces
         (
             (0 3 3 0)
          );
     }
 );

 mergePatchPairs
(
 );

 // ************************************************************************* //

Results:

Blockmesh

Pressure profile

Velocity profile

 

For wedge boundary conditions.

For the angle 45 degrees the blockMeshDict file.

/*--------------------------------*- 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 1; 
vertices 
( 
      (0.000000 0.000000 0.000000) 
      (0.000000 0.009239 0.003827) 
      (0.000000 0.009239 -0.003827) 
      (3.020000 0.000000 0.000000) 
      (3.020000 0.009239 0.003827) 
      (3.020000 0.009239 -0.003827) 
); 
blocks 
( 
      hex (0 2 1 0 3 5 4 3) (40 1 500) simpleGrading (1 1 1) 
); 
edges 
( 
arc 1 2 (0.000000 0.010000 0.000000) 
arc 4 5 (3.020000 0.010000 0.000000) 
); 
boundary 
( 
      inlet 
      { 
           type patch; 
           faces 
          ( 
           (0 1 2 0)
          );
      }
      outlet
      {
        type patch;
         faces
         (
            (3 4 5 3)
          );
     }
     wedge_front
     {
         type wedge;
         faces
         (
             (2 5 3 0)

          );
     }
     wedge_back
     {
         type wedge;
         faces
         (

             (0 3 4 1)
          );
     }
     top
     {
         type wall;
         faces
         (
             (1 4 5 2)
          );
     }
     axis
     {
         type empty;
         faces
         (
             (0 3 3 0)
          );
     }
 );

 mergePatchPairs
(
 );

 // ************************************************************************* //

 

 

Results:

Blockmesh

Pressure profile

Velocity profile

 

 

 

CONCLUSION.

We can observe there is no much change in the maximum velocities. The only difference observed is execution time a little

high than type-wedge. While the angle increases the error in the maximum velocity is decreased. Though there is not much

difference between analytical and numerical results and also with increasing the angle of the wedge, The Courant number is

slightly increasing. we can say that the symmetric boundary condition plays an important role in terms of accuracy in the

result

Leave a comment

Thanks for choosing to leave a comment. Please keep in mind that all the comments are moderated as per our comment policy, and your email will not be published for privacy reasons. Please leave a personal & meaningful conversation.

Please  login to add a comment

Other comments...

No comments yet!
Be the first to add a comment

Read more Projects by Yogessvaran T (50)

Week 14 challenge

Objective:

ASSEMBLY OF BUTTERFLY VALVE:- 1.All the parts that are saved in the required folder is opened in the Add components tool box. 2.Now using the Move option and Assembly Constraints option the different parts are joined together with the help of Align,touch,infer/Axis operations. 3. Finally,the assembly of butterfly valve…

calendar

18 Feb 2023 09:34 AM IST

    Read more

    Project - Position control of mass spring damper system

    Objective:

    To design a closed loop control scheme for a DC motor the following changes need to be done in the model in the previously created model. Speed is the controllable parameter, so we will set the reference speed in step block as 10,20, 40 whichever you want.  Subtract the actual speed from the reference speed to generate…

    calendar

    21 Jan 2023 10:29 AM IST

    • MATLAB
    Read more

    Project - Analysis of a practical automotive wiring circuit

    Objective:

    Identify each of the major elements in the above automotive wiring diagram.  Ans: Major Elements in the above automotive wiring diagram are - Genarator,                   Battery,              …

    calendar

    14 Dec 2022 03:37 AM IST

      Read more

      Week 6 - Data analysis

      Objective:

      -

      calendar

      04 Dec 2022 11:06 AM IST

        Read more

        Schedule a counselling session

        Please enter your name
        Please enter a valid email
        Please enter a valid number

        Related Courses

        coursecardcoursetype

        Post Graduate Program in CFD Solver Development

        4.8

        119 Hours of Content

        coursecard

        Introduction to OpenFOAM Development

        4.9

        18 Hours of Content

        coursecardcoursetype

        Post Graduate Program in Battery Technology for Mechanical Engineers

        4.8

        81 Hours of Content

        coursecard

        Simulation and Design of Power Converters for EV using MATLAB and Simulink

        4.9

        22 Hours of Content

        coursecard

        Introduction to Hybrid Electric Vehicle using MATLAB and Simulink

        4.8

        23 Hours of Content

        Schedule a counselling session

        Please enter your name
        Please enter a valid email
        Please enter a valid number

        logo

        Skill-Lync offers industry relevant advanced engineering courses for engineering students by partnering with industry experts.

        https://d27yxarlh48w6q.cloudfront.net/web/v1/images/facebook.svghttps://d27yxarlh48w6q.cloudfront.net/web/v1/images/insta.svghttps://d27yxarlh48w6q.cloudfront.net/web/v1/images/twitter.svghttps://d27yxarlh48w6q.cloudfront.net/web/v1/images/youtube.svghttps://d27yxarlh48w6q.cloudfront.net/web/v1/images/linkedin.svg

        Our Company

        News & EventsBlogCareersGrievance RedressalSkill-Lync ReviewsTermsPrivacy PolicyBecome an Affiliate
        map
        EpowerX Learning Technologies Pvt Ltd.
        4th Floor, BLOCK-B, Velachery - Tambaram Main Rd, Ram Nagar South, Madipakkam, Chennai, Tamil Nadu 600042.
        mail
        info@skill-lync.com
        mail
        ITgrievance@skill-lync.com

        Top Individual Courses

        Computational Combustion Using Python and CanteraIntroduction to Physical Modeling using SimscapeIntroduction to Structural Analysis using ANSYS WorkbenchIntroduction to Structural Analysis using ANSYS Workbench

        Top PG Programs

        Post Graduate Program in Hybrid Electric Vehicle Design and AnalysisPost Graduate Program in Computational Fluid DynamicsPost Graduate Program in CADPost Graduate Program in Electric Vehicle Design & Development

        Skill-Lync Plus

        Executive Program in Electric Vehicle Embedded SoftwareExecutive Program in Electric Vehicle DesignExecutive Program in Cybersecurity

        Trending Blogs

        Heat Transfer Principles in Energy-Efficient Refrigerators and Air Conditioners Advanced Modeling and Result Visualization in Simscape Exploring Simulink and Library Browser in Simscape Advanced Simulink Tools and Libraries in SimscapeExploring Simulink Basics in Simscape

        © 2025 Skill-Lync Inc. All Rights Reserved.

                    Do You Want To Showcase Your Technical Skills?
                    Sign-Up for our projects.