problem when use foamToVTK

ouafa · November 4, 2009, 11:11

hi,

I'm trying to visualize my mesh witch i already verify with chekMesh and every things seems to be ok, using foamToVTK, since I use OpenFoam for simulation. the mesh was created under gmsh and converted by gmshToFoam command.

I want to do this before actually computing anything but it crach with the following error message:

Time 0
volScalarFields : p nuTilda k omega
volVectorFields : U
volTensorFields : R

Internal : "/home/ouafa.bouloumou/work/RANScanal/jonction1/VTK/jonction1_1.vtk"
Original cells:4430284 points:4558642 Additional cells:0 additional points:0

Cannot find face 13418476 in any of the patches
8
(
wall2
wall6
inlet2
wall4
wall5
inlet1
outlet
wall1
)

It seems your patches are not consistent with the mesh : internalFaces:13163224 total number of faces:13418480#0 Foam::error:

rintStack(Foam::Ostream&) in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam:

olyBoundaryMesh::whichPatch(int) const in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#3 Foam:

ointPatchInterpolation::makePatchPatchWeigh ts() in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam:

ointPatchInterpolation:

ointPatchInterpola tion(Foam::fvMesh const&, Foam:

ointMesh const&) in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::volPointInterpolation::volPointInterpolation (Foam::fvMesh const&, Foam:

ointMesh const&) in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 main in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/foamToVTK"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/foamToVTK"

From function label polyBoundaryMesh::whichPatch(const label faceIndex) const
in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 402.

FOAM aborting

Aborted

Any ideas? Could you help me plz?
Ouafa

ICL · November 4, 2009, 11:45

Have you tried checkMesh before running your model?

Ali

ouafa · November 4, 2009, 11:50

thank you Ali,
yes i've tested checkMesh and it seems ok.
Ouafa

ICL · November 4, 2009, 11:59

To me, there is definitely something inconsistent in your mesh. Check your mesh again and try to use a sketch for better understanding.

Ali

ouafa · November 4, 2009, 12:35

ok, i re-check the mesh and i get this:

Time = constant
Mesh stats
points: 1080154
faces: 3160732
internal faces: 3082172
cells: 1040484
boundary patches: 8
point zones: 0
face zones: 1
cell zones: 57
Number of cells of each type:
hexahedra: 1040484
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0
Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
wall2 23722 24113 ok (non-closed singly connected)
wall6 192 250 ok (non-closed singly connected)
inlet2 960 1009 ok (non-closed singly connected)
wall4 10290 10578 ok (non-closed singly connected)
wall5 10290 10578 ok (non-closed singly connected)
inlet1 4116 4257 ok (non-closed singly connected)
outlet 4116 4257 ok (non-closed singly connected)
wall1 24682 25026 ok (non-closed singly connected)
Checking geometry...
Overall domain bounding box (-0.1 -0.04 -0.05) (0.15 0.04 0.05)
Mesh (non-empty) directions (1 1 1)
Mesh (non-empty, non-wedge) dimensions 3
Boundary openness (6.0135e-17 1.02144e-14 3.04769e-15) OK.
Max cell openness = 2.75831e-16 OK.
Max aspect ratio = 33.6809 OK.
Minumum face area = 1.91799e-08. Maximum face area = 1.71079e-05. Face area magnitudes OK.
Min volume = 7.91111e-12. Max volume = 7.74162e-09. Total volume = 0.00101256. Cell volumes OK.
Mesh non-orthogonality Max: 42.2568 average: 4.83306
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.586609 OK.
Mesh OK.
End

maybe it is about the "non-closed singly connected" topology?
shoud i redefine all my patches?
ouafa.

ICL · November 4, 2009, 12:48

The code start the calculation, right?
If it is, please do foamToVTK -time (the first time here) and see the result.

ouafa · November 4, 2009, 12:57

yes the calculations are started but even with the option ( -time 1 ) i still get the same error!!!

November 4, 2009, 11:11	problem when use foamToVTK	#1
ouafa New Member ouafa Join Date: Jul 2009 Posts: 15 Rep Power: 17	hi, I'm trying to visualize my mesh witch i already verify with chekMesh and every things seems to be ok, using foamToVTK, since I use OpenFoam for simulation. the mesh was created under gmsh and converted by gmshToFoam command. I want to do this before actually computing anything but it crach with the following error message: Time 0 volScalarFields : p nuTilda k omega volVectorFields : U volTensorFields : R Internal : "/home/ouafa.bouloumou/work/RANScanal/jonction1/VTK/jonction1_1.vtk" Original cells:4430284 points:4558642 Additional cells:0 additional points:0 Cannot find face 13418476 in any of the patches 8 ( wall2 wall6 inlet2 wall4 wall5 inlet1 outlet wall1 ) It seems your patches are not consistent with the mesh : internalFaces:13163224 total number of faces:13418480#0 Foam::error:rintStack(Foam::Ostream&) in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #2 Foam:olyBoundaryMesh::whichPatch(int) const in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #3 Foam:ointPatchInterpolation::makePatchPatchWeigh ts() in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #4 Foam:ointPatchInterpolation:ointPatchInterpola tion(Foam::fvMesh const&, Foam:ointMesh const&) in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #5 Foam::volPointInterpolation::volPointInterpolation (Foam::fvMesh const&, Foam:ointMesh const&) in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so" #6 main in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/foamToVTK" #7 __libc_start_main in "/lib64/libc.so.6" #8 Foam::regIOobject::readIfModified() in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/foamToVTK" From function label polyBoundaryMesh::whichPatch(const label faceIndex) const in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 402. FOAM aborting Aborted Any ideas? Could you help me plz? Ouafa

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
UDF compiling problem	Wouter	Fluent UDF and Scheme Programming	6	June 6, 2012 05:43
Incoherent problem table in hollow-fiber spinning	Gianni	FLUENT	0	April 5, 2008 11:33
natural convection problem for a CHT problem	Se-Hee	CFX	2	June 10, 2007 07:29
Adiabatic and Rotating wall (Convection problem)	ParodDav	CFX	5	April 29, 2007 20:13
Is this problem well posed?	Thomas P. Abraham	Main CFD Forum	5	September 8, 1999 15:52

November 4, 2009, 11:45		#2
ICL New Member Ali Sh Join Date: Jun 2009 Location: London Posts: 28 Rep Power: 17	Have you tried checkMesh before running your model? Ali

November 4, 2009, 11:50		#3
ouafa New Member ouafa Join Date: Jul 2009 Posts: 15 Rep Power: 17	thank you Ali, yes i've tested checkMesh and it seems ok. Ouafa

November 4, 2009, 11:59		#4
ICL New Member Ali Sh Join Date: Jun 2009 Location: London Posts: 28 Rep Power: 17	To me, there is definitely something inconsistent in your mesh. Check your mesh again and try to use a sketch for better understanding. Ali

November 4, 2009, 12:35		#5
ouafa New Member ouafa Join Date: Jul 2009 Posts: 15 Rep Power: 17	ok, i re-check the mesh and i get this: Time = constant Mesh stats points: 1080154 faces: 3160732 internal faces: 3082172 cells: 1040484 boundary patches: 8 point zones: 0 face zones: 1 cell zones: 57 Number of cells of each type: hexahedra: 1040484 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface topology wall2 23722 24113 ok (non-closed singly connected) wall6 192 250 ok (non-closed singly connected) inlet2 960 1009 ok (non-closed singly connected) wall4 10290 10578 ok (non-closed singly connected) wall5 10290 10578 ok (non-closed singly connected) inlet1 4116 4257 ok (non-closed singly connected) outlet 4116 4257 ok (non-closed singly connected) wall1 24682 25026 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (-0.1 -0.04 -0.05) (0.15 0.04 0.05) Mesh (non-empty) directions (1 1 1) Mesh (non-empty, non-wedge) dimensions 3 Boundary openness (6.0135e-17 1.02144e-14 3.04769e-15) OK. Max cell openness = 2.75831e-16 OK. Max aspect ratio = 33.6809 OK. Minumum face area = 1.91799e-08. Maximum face area = 1.71079e-05. Face area magnitudes OK. Min volume = 7.91111e-12. Max volume = 7.74162e-09. Total volume = 0.00101256. Cell volumes OK. Mesh non-orthogonality Max: 42.2568 average: 4.83306 Non-orthogonality check OK. Face pyramids OK. Max skewness = 0.586609 OK. Mesh OK. End maybe it is about the "non-closed singly connected" topology? shoud i redefine all my patches? ouafa.

November 4, 2009, 12:48		#6
ICL New Member Ali Sh Join Date: Jun 2009 Location: London Posts: 28 Rep Power: 17	The code start the calculation, right? If it is, please do foamToVTK -time (the first time here) and see the result.

November 4, 2009, 12:57		#7
ouafa New Member ouafa Join Date: Jul 2009 Posts: 15 Rep Power: 17	yes the calculations are started but even with the option ( -time 1 ) i still get the same error!!!