[WolfgangS.]

hi guys!

I am doing my first steps to openfoam .. and i need yours help!

I like to simulate a laminar and turbulent flow over a circular cylinder (RANS k-w-Modell, LES). Later i want to simulate the magnus (flow over a rotating circlar cylinder) and inverse magnus effect (maby you know it from sports) for a studyproject.

First i made a big block (will be meshed with blockMesh) in OpenFoam and i want a smaler block, which will be meshed with sappyHexMesh. The cylinder will lie in smaler block.

So i started:

FoamFile
{ version 2.0; format ascii; class dictionary; object blockMeshDict;}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;
vertices
(
(110 0 -200) //0 front circle
(77.781746 77.781746 -200) //1 45 Degrees
(0 110 -200) //2
(-77.781746 77.781746 -200) //3
(-110 0 -200) //4
(-77.781746 -77.781746 -200) //5
(0 -110 -200) //6
(77.781746 -77.781746 -200) //7
(110 0 200) //8 back circle
(77.781746 77.781746 200) //9
(0 110 200) //10
(-77.781746 77.781746 200) //11
(-110 0 200) //12
(-77.781746 -77.781746 200) //13
(0 -110 200) //14
(77.781746 -77.781746 200) //15
(-1100 -1100 -200) //16 front big block
(3300 -1100 -200) //17
(3300 1100 -200) //18
(-1100 1100 -200) //19
(-1100 -1100 200) //20 back big block
(3300 -1100 200) //21
(3300 1100 200) //22
(-1100 1100 200) //23
);

blocks
(
hex (16 17 18 19 20 21 22 23) (1 1 1) simpleGrading (1 1 1)
);

edges
(
arc 0 1 (95.262794 55 -200) //front circle
arc 1 2 (55 95.262794 -200) //30 and 60 Degrees ( cos() and sin() )
arc 2 3 (-55 95.262794 -200)
arc 3 4 (-95.262794 55 -200)
arc 4 5 (-95.262794 -55 -200)
arc 5 6 (-55 -95.262794 -200)
arc 6 7 (55 -95.262794 -200)
arc 7 0 (95.262794 -55 -200)
arc 0 1 (95.262794 55 -200) //back circle
arc 1 2 (55 95.262794 -200)
arc 2 3 (-55 95.262794 -200)
arc 3 4 (-95.262794 55 -200)
arc 4 5 (-95.262794 -55 -200)
arc 5 6 (-55 -95.262794 -200)
arc 6 7 (55 -95.262794 -200)
arc 7 0 (95.262794 -55 -200)
);

patches
(
patch outlet
(
(18 22 21 17)
)
patch inlet
(
(16 20 23 19)
)
patch fixedWalls
(
(19 23 22 18)
(17 21 20 16)
)
patch frontAndBack
(
(16 19 18 17)
(20 21 22 23)
)
patch hole // I think it must be here, where is something wrong
(
(2 1 9 10)
(1 0 8 9)
(0 7 15 8)
(7 6 14 15)
(6 5 13 14)
(5 4 12 13)
(4 3 11 12)
(3 2 10 11)
)
);

mergePatchPairs
(
);

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


Unfortunately it doesn´t work . I always get this information:

Default patch type set to emptyface 0 in patch 4 does not have neighbour cell face: 4(2 1 9 10)#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"#1 Foam::error::abort() in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"#2 Foam::Ostream& Foam:perator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/blockMesh"#3 Foam:olyMesh::facePatchFaceCells(Foam::List<Foam ::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"#4 Foam:olyMesh:olyMesh(Foam::IOobject const&, Foam::Xfer<Foam::Field<Foam::Vector<double> > > const&, Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<Foam::word> const&, Foam::word const&, Foam::word const&, Foam::List<Foam::word> const&, bool) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"#5 Foam::blockMesh::createTopology(Foam::IOdictionary &) in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/blockMesh"#6 Foam::blockMesh::blockMesh(Foam::IOdictionary&) in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/blockMesh"#7 main in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/blockMesh"#8 __libc_start_main in "/lib64/libc.so.6"#9 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116 From function polyMesh::facePatchFaceCells(const faceList& patchFaces,const labelListList& pointCells,const faceListList& cellsFaceShapes,const label patchID) in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 125.FOAM aborting
What´s wrong ? I do not realy understand where the problem is? I think it must be in the "patch hole".

Please! Could you help me!

!!!!!!!!!! I have a further question: !!!!!!!!!!!!!

Maby this Mesh i have created will run, how could i rotate the cylinder? or will i need there a flow around the cylinder, which will simulate the rotation?
Would it be better to simulate a flow over a CAD created cylinder (.stl) ????

maby there is somebody with experience with cylinders!

Thank you in anticipation!!


Wolfgang