CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[blockMesh] Axisymmetric blockMesh error. Are Arcs needed?

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   November 2, 2010, 11:25
Default Axisymmetric blockMesh error. Are Arcs needed?
New Member
Join Date: Sep 2010
Posts: 1
Rep Power: 0
andygdl is on a distinguished road

I'm having problems meshing a fairly complicated axisymmetric nozzle. I've had no problems in the past with axisymmetric meshes, but this one is causing problems. The error is as follows:

Create time

Creating block mesh from

Creating blockCorners

Creating curved edges
#0  Foam::error::printStack(Foam::Ostream&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/"
#1  Foam::sigFpe::sigFpeHandler(int) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/"
#2  Uninterpreted:
#3  Foam::coordinateRotation::calcTransform(Foam::Vector<double> const&, Foam::Vector<double> const&, Foam::coordinateRotation::axisOrder const&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/"
#4  Foam::coordinateRotation::coordinateRotation(Foam::Vector<double> const&, Foam::Vector<double> const&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/"
#5  Foam::coordinateSystem::coordinateSystem(Foam::word const&, Foam::Vector<double> const&, Foam::Vector<double> const&, Foam::Vector<double> const&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/"
#6  Foam::cylindricalCS::cylindricalCS(Foam::word const&, Foam::Vector<double> const&, Foam::Vector<double> const&, Foam::Vector<double> const&, bool) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/"
#7  Foam::arcEdge::calcAngle() in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/blockMesh"
#8  Foam::arcEdge::arcEdge(Foam::Field<Foam::Vector<double> > const&, Foam::Istream&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/blockMesh"
#9  Foam::curvedEdge::addIstreamConstructorToTable<Foam::arcEdge>::New(Foam::Field<Foam::Vector<double> > const&, Foam::Istream&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/blockMesh"
#10  Foam::curvedEdge::New(Foam::Field<Foam::Vector<double> > const&, Foam::Istream&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/blockMesh"
#11  Foam::blockMesh::createTopology(Foam::IOdictionary&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/blockMesh"
#12  Foam::blockMesh::blockMesh(Foam::IOdictionary&) in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/blockMesh"
#13  main in "/home/andy/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/blockMesh"
#14  __libc_start_main in "/lib/"
#15  _start at /usr/src/packages/BUILD/glibc-2.10.1/csu/../sysdeps/i386/elf/start.S:122
Floating point exception

The problem goes away if I replace the 'arcs' with lines - i.e. the arcs put in for rotation about the axis (not the arcs that define the shape of the geometry!).

I've checked, the interpolation points all seem in order. I don't think its any problem with 'crazy' surfaces bounded by 2 splines and 2 arcs, as I've used a very similar geometry before and not had any problems.

The obvious question is: 'do I actually need to make them arcs?'. It is only 5 degrees, would it make any difference to the solution if they were just lines?

Otherwise, does anyone have an idea what the problem is? Thanks very much for any help!

I've included the that doesn't work and the blockMeshDict.yay that works without arcs

ps. I realise I didn't need to define 'line a b'... I was young and crazy!
Attached Files
File Type: txt blockMeshDict.yay.txt (4.6 KB, 31 views)
File Type: txt (5.4 KB, 16 views)
andygdl is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[blockMesh] blockMesh gives "Inconsistent point locations" for axisymmetric problem irving OpenFOAM Meshing & Mesh Conversion 7 January 27, 2021 12:31
[blockMesh] Trouble using blockMesh for axisymmetric geometry twinturbotom OpenFOAM Meshing & Mesh Conversion 1 February 8, 2016 20:59
[swak4Foam] Installing swak4Foam to OpenFOAM in mac Kaquesang OpenFOAM Community Contributions 22 January 21, 2013 11:51
[Technical] BlockMesh not making arcs Disco_Caine OpenFOAM Meshing & Mesh Conversion 3 July 6, 2010 10:16
Axisymmetric cases Test with blockmesh and icoFoam frackowi OpenFOAM Running, Solving & CFD 2 July 11, 2007 19:02

All times are GMT -4. The time now is 13:44.