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

foamToVTK ERROR, help me!

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 1, 2012, 11:53
Default foamToVTK ERROR, help me!
  #1
New Member
 
Francesco
Join Date: Jun 2012
Location: Rome, Italy
Posts: 23
Rep Power: 5
Bombolati is on a distinguished road
Hi at all!
I've done a simulation with reactingFoam. Now I've a folder for every writing point that is fixed in the controlDict file. The next step is postprocess the case with paraFoam but I'm working in a parallel environment and there is only paraview installed. So I thought to use the foamToVTK utility to visualize the case with paraview. When I run this utility this error encountered:

--> FOAM FATAL ERROR:
More than one patch accessing the same transform but not of the same sign.
patch:CameraFL1_shadow transform:0 sign:-1 current transforms ( -1 0 0)

From function Foam::label Foam::globalIndexAndTransform::addToTransformIndex
(
const label,
const label,
const bool
) const

in file lnInclude/globalIndexAndTransformI.H at line 240.

FOAM exiting


I've done the mesh in gambit ad I've exported it with fluent3DMeshToFoam utility. The mesh has some surfaces defined as cyclic, CameraFL1_shadow is one of this.
Someone know how this error encountered?
Thanks a lot.
Bombolati is offline   Reply With Quote

Old   October 1, 2012, 17:21
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Bombolati,

I'm not sure at which level you're doing parallel processing, but if in the cluster you can load the whole case in a single machine for post-processing, then you could use ParaView's internal ".foam" reader. Run this command inside your case folder:
Code:
touch case.foam
Then open the "case.foam" file in ParaView. You then have a choice to load either the reconstructed or parallel case.

As for the problem with foamToVTK: did you check the mesh before running the case?
Code:
checkMesh -constant
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   October 2, 2012, 04:38
Default
  #3
New Member
 
Francesco
Join Date: Jun 2012
Location: Rome, Italy
Posts: 23
Rep Power: 5
Bombolati is on a distinguished road
Hi Bruno and thanks a lot for you answer. I used the first command to create the file <<case.foam>> and all is ok. Next I tried to use the utility <<check mesh>> and some problems come out. I post here below the result:

Create time

Create polyMesh for time = constant

Time = constant

Mesh stats
points: 1527113
faces: 17663239
internal faces: 17406505
cells: 8766526
boundary patches: 20
point zones: 0
face zones: 4
cell zones: 1

Overall number of cells of each type:
hexahedra: 0
prisms: 0
wedges: 0
pyramids: 3640
tet wedges: 0
tetrahedra: 8762886
polyhedra: 0

Checking topology...
Boundary definition OK.
Cell to face addressing 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
CameraFL2_shadow 9308 4792 ok (non-closed singly connected)
CameraFL2 9308 4792 ok (non-closed singly connected)
CameraFL1_shadow 19082 9725 ok (non-closed singly connected)
CameraFL1 19082 9725 ok (non-closed singly connected)
MixerFL_shadow 1820 1908 ok (non-closed singly connected)
MixerFL 1820 1908 ok (non-closed singly connected)
IniettoreAriaFL_shadow15463 7907 ok (non-closed singly connected)
IniettoreAriaFL 15463 7907 ok (non-closed singly connected)
CameraFSup 40844 20654 ok (non-closed singly connected)
IngAria 12544 6427 ok (non-closed singly connected)
CombIng2 1684 891 ok (non-closed singly connected)
CombIng1 1682 890 ok (non-closed singly connected)
CameraFL 10446 5408 ok (non-closed singly connected)
MixerFint 7254 3817 ok (non-closed singly connected)
MixerFext 11222 5788 ok (non-closed singly connected)
IniettoreAriaFint 12196 6261 ok (non-closed singly connected)
TestaPiano 1959 1040 ok (non-closed singly connected)
IniettoreAriaFext 36654 18584 ok (non-closed singly connected)
Bordo 12707 6474 ok (non-closed singly connected)
CameraFinf 16196 8269 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-1.29063e-15 -2.17066e-14 -7.10543e-15) (98 98 57)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-3.55779e-15 -4.11417e-15 2.20632e-15) OK.
Max cell openness = 3.30458e-16 OK.
Max aspect ratio = 8.87394 OK.
Minumum face area = 0.00211186. Maximum face area = 1.39301. Face area magnitudes OK.
Min volume = 5.43598e-05. Max volume = 0.441604. Total volume = 231604. Cell volumes OK.
Mesh non-orthogonality Max: 67.2498 average: 19.7284
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.954274 OK.
**Error in coupled point location: 45673 faces have their 0th vertex not opposite their coupled equivalent. Average mismatch 38.8115.
<<Writing 45673 faces with incorrectly matched 0th vertex to set coupledFaces

Failed 1 mesh checks.

End

What do you think?
Bombolati is offline   Reply With Quote

Old   October 2, 2012, 15:58
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Bombolati,

Quote:
Originally Posted by Bombolati View Post
**Error in coupled point location: 45673 faces have their 0th vertex not opposite their coupled equivalent. Average mismatch 38.8115.
<<Writing 45673 faces with incorrectly matched 0th vertex to set coupledFaces
That is enough proof for me to conclude that this is the reason why foamToVTK is complaining!

You'll find a usage example for fluent3DMeshToFoam: http://openfoamwiki.net/index.php/Fl...y_step_example
As you can see, in that example it also indicates that you should check the mesh after conversion

I'm not familiar with fluent3DMeshToFoam, so my advice is to look in the forum for the keywords:
Code:
fluent3DMeshToFoam cyclic
I'm guessing there is a protocol that must be followed for cyclics to be properly converted.

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Reply

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
foamToVTK Astarta OpenFOAM Bugs 2 September 2, 2011 13:44
Command foamToVTK gruber OpenFOAM 0 July 26, 2010 09:44
foamToVTK sameer_kumar OpenFOAM Post-Processing 1 May 6, 2010 21:17
FoamTOVTK yapalparvi OpenFOAM Post-Processing 3 August 12, 2009 08:19
FoamToVTK error with OF 13 melanie OpenFOAM Paraview & paraFoam 1 May 22, 2006 04:40


All times are GMT -4. The time now is 18:30.