CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   2D axisymmetric Mesh problem in OpenFoam (https://www.cfd-online.com/Forums/openfoam/74964-2d-axisymmetric-mesh-problem-openfoam.html)

Javed April 13, 2010 06:22

2D axisymmetric Mesh problem in OpenFoam
 
I am creating 2D axisymmetric geometry using follwoing blockMesh file in OpenFoam...But not working:confused:
It will be great if u plz let help me for the following pbm..

blocks
(
hex (0 1 2 3 0 1 6 7) (20 20 1) simpleGrading (1 1 1)
);

edges
(
);

patches
(
wall fixed
(
(3 7 6 2)
)
patch inlet
(
(0 0 7 3)
)
patch outlet
(
(2 6 1 1)
)
empty center
(
(1 1 0 0)
)
wedge front
(
(0 3 2 1)
)
wedge back
(
(0 1 6 7)
)
);

but getting follwoing error..
Creating block mesh from
"/home/javed/OpenFOAM/javed-1.6/run/wedge-new/constant/polyMesh/blockMeshDict"


Creating blockCorners

Creating curved edges

Creating blocks

Creating patches

Creating block mesh topology

Default patch type set to empty
#0 Foam::error::printStack(Foam::Ostream&) in "/home/javed/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/javed/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::wedgePolyPatch::initTransforms() in "/home/javed/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::wedgePolyPatch::wedgePolyPatch(Foam::word const&, int, int, int, Foam::polyBoundaryMesh const&) in "/home/javed/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::polyPatch::addwordConstructorToTable<Foam::w edgePolyPatch>::New(Foam::word const&, int, int, int, Foam::polyBoundaryMesh const&) in "/home/javed/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::polyPatch::New(Foam::word const&, Foam::word const&, int, int, int, Foam::polyBoundaryMesh const&) in "/home/javed/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::polyMesh::polyMesh(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 "/home/javed/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#8 Foam::blockMesh::createTopology(Foam::IOdictionary &) in "/home/javed/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/blockMesh"
#9 Foam::blockMesh::blockMesh(Foam::IOdictionary&) in "/home/javed/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/blockMesh"
#10 main in "/home/javed/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/blockMesh"
#11 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#12 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122
Floating point exception
..

Thanks for the help...

hamcer June 9, 2010 04:17

hi Javed,

Perhaps it is late , but I think your pb is in the block direction (change direction see cavity case) or mistakes in coordinate points.

hope that help you.

drrbradford August 25, 2011 11:53

It looks like you've tried to create a wedge yourself:

Code:

hex (0 1 2 3 0 1 6 7)
Look at the makeAxialMesh function which can create a wedge mesh from a 2D plane mesh relatively quickly and easily.

http://www.cfd-online.com/Forums/ope...symmetric.html

Hisham September 22, 2011 02:23

It is most probably an OpenFOAM version problem. The pasted text is, I believe, from a Wiki page for an older version. Currently, the blockMesh syntax is different. Check with a tutorial to see the changes you need to apply to your boundaries sub dictionary. It should work smoothly thereafter!

Hisham


All times are GMT -4. The time now is 22:41.