CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

problem when use foamToVTK

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 4, 2009, 11:11
Default problem when use foamToVTK
  #1
New Member
 
ouafa
Join Date: Jul 2009
Posts: 15
Rep Power: 16
ouafa is on a distinguished road
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
ouafa is offline   Reply With Quote

Old   November 4, 2009, 11:45
Default
  #2
ICL
New Member
 
Ali Sh
Join Date: Jun 2009
Location: London
Posts: 28
Rep Power: 16
ICL is on a distinguished road
Have you tried checkMesh before running your model?

Ali
ICL is offline   Reply With Quote

Old   November 4, 2009, 11:50
Default
  #3
New Member
 
ouafa
Join Date: Jul 2009
Posts: 15
Rep Power: 16
ouafa is on a distinguished road
thank you Ali,
yes i've tested checkMesh and it seems ok.
Ouafa
ouafa is offline   Reply With Quote

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

Ali
ICL is offline   Reply With Quote

Old   November 4, 2009, 12:35
Default
  #5
New Member
 
ouafa
Join Date: Jul 2009
Posts: 15
Rep Power: 16
ouafa is on a distinguished road
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.
ouafa is offline   Reply With Quote

Old   November 4, 2009, 12:48
Default
  #6
ICL
New Member
 
Ali Sh
Join Date: Jun 2009
Location: London
Posts: 28
Rep Power: 16
ICL is on a distinguished road
The code start the calculation, right?
If it is, please do foamToVTK -time (the first time here) and see the result.
ICL is offline   Reply With Quote

Old   November 4, 2009, 12:57
Default
  #7
New Member
 
ouafa
Join Date: Jul 2009
Posts: 15
Rep Power: 16
ouafa is on a distinguished road
yes the calculations are started but even with the option ( -time 1 ) i still get the same error!!!
ouafa is offline   Reply With Quote

Reply

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
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


All times are GMT -4. The time now is 08:33.