problem when use foamToVTK
 

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::printStack(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::polyBoundaryMesh::whichPatch(int) const in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#3 Foam::pointPatchInterpolation::makePatchPatchWeigh ts() in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::pointPatchInterpolation::pointPatchInterpola tion(Foam::fvMesh const&, Foam::pointMesh const&) in "/home/ouafa.bouloumou/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::volPointInterpolation::volPointInterpolation (Foam::fvMesh const&, Foam::pointMesh 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

Have you tried checkMesh before running your model?

Ali

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

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

Ali

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.

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

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