CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Native Meshers: blockMesh (
-   -   BlockMesh error with growing mesh size (

kian September 21, 2007 00:08

Hello, I just upgraded to 1

I just upgraded to 1.4.1 by compiling the source.
I have noticed that blockMesh does not like more than 962 cells in the x-direction in the sonicFoam/shockTube example (1 block), and it gives this error if you specify more than 962 cells:

--> FOAM FATAL ERROR : face 0 in patch 0 does not have neighbour cell face: 4(0 2 3 1)#0

I don't know if this is related to the problem or not, since I had not started to really use version 1.4.1, but I had made changes to the source of sonicFoam and executed Allwmake in the applications directory. I saw some errors in vtkFoam application (different gcc versions), but the modified sonicFoam application works fine.

Thank you for your help.

Kian Mehravaran

msrinath80 September 21, 2007 00:11

The error message seems to sug
The error message seems to suggest that there is a problem in your blockMeshDict. Paste it here. Someone might be able to identify the source of the problem.

kian September 21, 2007 10:01

Dear Srinath, Thank you for y
Dear Srinath,
Thank you for your reply.
It is the one in sonicFoam/shockTube. I haven't changed the blockMeshDict. As I said, the error goes away if the mesh size is 962, instead of the original 1000! I feel that something is wrong with my blockMesh. The full error is below the letter.
Thanks, Kian Mehravaran
--> FOAM FATAL ERROR : face 0 in patch 0 does not have neighbour cell face: 4(0 1 3 2)#0 Foam::error::printStack (Foam: in "/home/kian/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccSPOpt/"
#1 Foam::error::abort() in "/home/kian/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccSPOpt/"
#2 Foam::polyMesh::facePatchFaceCells(Foam::List<foam ::face> const&, Foam::List<foam::list<int> > const&, Foam:: List<foam::list<foam::face> > const&, int) const in "/home/kian/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccSPOpt/libOpen"
#3 Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Field<foam::vector<float> > const&, Foam::List<foam::ce> const&, Foam::List<foam::list<foam::face> > const&, Foam::List<foam::word> const&, Foam::List<foam::word> const&, Foam::word const&, Foam::List<foam::word> const&, bool) in "/home/kian/OpenFOAM/OpenFOAM-1.4.1/lib/linu xGccSPOpt/"
#4 main in "/home/kian/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccSPOpt/blockMesh"
#5 __libc_start_main in "/lib/tls/"
#6 __gxx_personality_v0 at ../sysdeps/i386/elf/start.S:122

msrinath80 September 21, 2007 12:09

Wierd. It works perfectly fine
Wierd. It works perfectly fine for me here. Compare your blockMeshDict with this one: blockMeshDict

kian September 24, 2007 16:00

Thank you for posting your blo
Thank you for posting your blockMeshDict, mine is exactly the same. Weird indeed.

Another weird problem that I had, reminded me of the difference our systems might have:
I compiled mine with single precision.
I noticed that when I use a deltaT lower than 1.0E-05, it is basically regarded as zero. In other words, time is always zero!
I checked the deltaT.value() in the code (sonicFoam), and it contains the correct value.

All times are GMT -4. The time now is 11:21.