# Creating curves in blockMesh (An Example)

In this post, I’ll give a simple example of how to create curves in blockMesh. For this example, we’ll look at the following basic setup:

As you can see, we’ll be simulating the flow over a bump defined by the curve:

$y=H*\sin\left(\pi x \right)$

First, let’s look at the basic blockMeshDict for this blocking layout WITHOUT any curves defined:

/*--------------------------------*- C++ -*----------------------------------*\
=========                 |
\\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
\\    /   O peration     | Website:  https://openfoam.org
\\  /    A nd           | Version:  6
\\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       dictionary;
object      blockMeshDict;
}

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

convertToMeters 1;

vertices
(
(-1 0 0)    // 0
(0 0 0)     // 1
(1 0 0)     // 2
(2 0 0)     // 3
(-1 2 0)    // 4
(0 2 0)     // 5
(1 2 0)     // 6
(2 2 0)     // 7

(-1 0 1)    // 8
(0 0 1)     // 9
(1 0 1)     // 10
(2 0 1)     // 11
(-1 2 1)    // 12
(0 2 1)     // 13
(1 2 1)     // 14
(2 2 1)     // 15
);

blocks
(
hex (0 1 5 4 8 9 13 12) (20 100 1) simpleGrading (0.1 10 1)
hex (1 2 6 5 9 10 14 13) (80 100 1) simpleGrading (1 10 1)
hex (2 3 7 6 10 11 15 14) (20 100 1) simpleGrading (10 10 1)
);

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(0 8 12 4)
);
}
outlet
{
type patch;
faces
(
(3 7 15 11)
);
}
lowerWall
{
type wall;
faces
(
(0 1 9 8)
(1 2 10 9)
(2 3 11 10)
);
}
upperWall
{
type patch;
faces
(
(4 12 13 5)
(5 13 14 6)
(6 14 15 7)
);
}
frontAndBack
{
type empty;
faces
(
(8 9 13 12)
(9 10 14 13)
(10 11 15 14)
(1 0 4 5)
(2 1 5 6)
(3 2 6 7)
);
}
);

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

This blockMeshDict produces the following grid:

It is best practice in my opinion to first make your blockMesh without any edges. This lets you see if there are any major errors resulting from the block topology itself. From the results above, we can see we’re ready to move on!

So now we need to define the curve. In blockMesh, curves are added using the edges sub-dictionary. This is a simple sub dictionary that is just a list of interpolation points:

edges
(
polyLine 1 2
(
(0	0       0)
(0.1	0.0309016994    0)
(0.2	0.0587785252    0)
(0.3	0.0809016994    0)
(0.4	0.0951056516    0)
(0.5	0.1     0)
(0.6	0.0951056516    0)
(0.7	0.0809016994    0)
(0.8	0.0587785252    0)
(0.9	0.0309016994    0)
(1	0       0)
)

polyLine 9 10
(
(0	0       1)
(0.1	0.0309016994    1)
(0.2	0.0587785252    1)
(0.3	0.0809016994    1)
(0.4	0.0951056516    1)
(0.5	0.1     1)
(0.6	0.0951056516    1)
(0.7	0.0809016994    1)
(0.8	0.0587785252    1)
(0.9	0.0309016994    1)
(1	0       1)
)
);

The sub-dictionary above is just a list of points on the curve $y=H\sin(\pi x)$. The interpolation method is polyLine (straight lines between interpolation points). An alternative interpolation method could be spline.

The following mesh is produced:

Hopefully this simple example will help some people looking to incorporate curved edges into their blockMeshing!

Cheers.

This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software via http://www.openfoam.com, and owner of theOPENFOAM®  andOpenCFD®  trademarks.

Categories

Tags