CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Native Meshers: blockMesh

Axisymmetric blockMesh error. Are Arcs needed?

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

LinkBack Thread Tools Display Modes
Old   November 2, 2010, 12: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, 28 views)
File Type: txt (5.4 KB, 15 views)
andygdl is offline   Reply With Quote


Thread Tools
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 14:53
BlockMesh not making arcs Disco_Caine OpenFOAM Meshing Format & General Technical 3 July 6, 2010 10:16
Blockmesh cavity error message tonitoney OpenFOAM Installation 2 March 17, 2008 12:59
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 04:10.