[psbharatkumar89]

Hi all,

I am a beginner in OpenFOAM and I am trying to create a cylinder using wedge extrusion in ExtrudeMesh. I created a slab using blockMesh and I want to extrude the face of the slab around an axis. But I couldnot do it succesfully. I am getting the same error no matter what changes I do. Pls help. My error is

"--> FOAM FATAL ERROR:
wedge inlet plane aligns with a coordinate plane.
The wedge plane should make a small angle (~2.5deg) with the coordinate plane
and the the pair of wedge planes should be symmetric about the coordinate plane.
Normal of wedge plane is (1.22465e-16 1 -8.34056e-34) , implied coordinate plane direction is (0 1 -0)

From function virtual void Foam::wedgePolyPatch::calcGeometry(Foam::PstreamBu ffers&)
in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 109"



My blockMeshDict and extrudeMeshDict are as follows for your reference. Thanks in advance.

blockMeshDict:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// General macros to create 2D/extruded-2D meshes


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

convertToMeters 0.001;

vertices
(
(20 0 0)
(65 0 0)
(64.9901001852 1.1344064184 0)
(19.9969539031 0.3490481287 0)

(20 0 40)
(65 0 40)
(64.9901001852 1.1344064184 40)
(19.9969539031 0.3490481287 40)
);

edges
(
line 0 1
line 3 2
line 7 6
line 4 5

arc 0 3 (19.9992384613 0.17453071 0)
arc 1 2 (64.9975249992 0.5672248074 0)
arc 5 6 (64.9975249992 0.5672248074 40)
arc 4 7 (19.9992384613 0.17453071 40)

line 0 4
line 1 5
line 2 6
line 3 7
);

blocks
(
hex (0 1 2 3 4 5 6 7)
(18 1 18)
simpleGrading
(
( //X dir
(0.2 0.3 4)
(0.6 0.4 1)
(0.2 0.3 0.25)
)
1 //Ydir
( //Z dir
(0.2 0.3 4)
(0.6 0.4 1)
(0.2 0.3 0.25)
)
)
);

boundary
(
inlet
{
type cyclic;
neighbourPatch outlet;
faces ((0 1 5 4));
}
outlet
{
type cyclic;
neighbourPatch inlet;
faces ((3 7 6 2));
}
bottom
{
type wall;
faces ((0 3 2 1));
}
top
{
type patch; //(or wall)
faces ((4 5 6 7));
}
ext_wall
{
type wall;
faces ((0 4 7 3));
}
inner_wall
{
type wall;
faces ((6 5 1 2));
}
);

ExtrudeMeshDict:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object extrudeMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// What to extrude:
// patch : from patch of another case ('sourceCase')
// mesh : as above but with original case included
// surface : from externally read surface

//constructFrom mesh;
constructFrom patch;
//constructFrom surface;

// If construct from patch/mesh:
sourceCase ;
sourcePatches (inlet);

// If construct from patch: patch to use for back (can be same as sourcePatch)
exposedPatchName outlet;

// If construct from surface:
//surface "surface_name.stl";

// Flip surface normals before usage. Valid only for extrude from surface or
// patch.
flipNormals false;

//- Linear extrusion in point-normal direction
//extrudeModel linearNormal;

//- Linear extrusion in specified direction
//extrudeModel linearDirection;

//- Wedge extrusion. If nLayers is 1 assumes symmetry around plane.
extrudeModel wedge;

//- Extrudes into sphere around (0 0 0)
//extrudeModel linearRadial;

//- Extrudes into sphere around (0 0 0) with specified radii
//extrudeModel radial;

//- Extrudes into sphere with grading according to pressure (atmospherics)
//extrudeModel sigmaRadial;

//nLayers 10;
nLayers 180;

//expansionRatio 1.0; //0.9;


sectorCoeffs
{
axisPt (0 0 0);
axis (0 0 1);
angle 360; // For nLayers=1 assume symmetry so angle/2 on each side
}


/*linearNormalCoeffs
{
thickness 0.1;
}*/

/*
linearDirectionCoeffs
{
direction (0 1 0);
thickness 0.05;
}
*/

/*
linearRadialCoeffs
{
R 0.1;
// Optional inner radius
Rsurface 0.01;
}
*/

/*
radialCoeffs
{
// Radii specified through interpolation table
R table ((0 0.01)(3 0.03)(10 0.1));
}
*/

/*
sigmaRadialCoeffs
{
RTbyg 1;
pRef 1;
pStrat 1;
}
*/

// Do front and back need to be merged? Usually only makes sense for 360
// degree wedges.
mergeFaces true; //true;

// Merge small edges. Fraction of bounding box.
//mergeTol 0;