CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Meshing Format & General Technical

How to get polyhedral mesh without additional cells when using foamToVTK

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 7, 2008, 11:54
Default I use salome to generate tetra
  #1
New Member
 
Martin/Run DU
Join Date: Mar 2009
Location: Hong Kong
Posts: 25
Rep Power: 8
chnrdu is on a distinguished road
Send a message via Skype™ to chnrdu
I use salome to generate tetrahedral mesh, then use ideasUnvToFoam to convert this mesh. I wanna use polyhedral mesh to compute some cases, so I use polyDualMesh to convert tetra to polyhedral mesh. That all is fine.

However, when I use foamToVTK to convert foam formated mesh to VTK formated mesh in order to show meshes in paraview 3.3.0, things are different, there are additional cells generated. I use this command:

foamToVTK . sonicLiquidPipeWithSymCylinDamperP100_mseq_poly -time 4.8828125e-06 -allPatches

It said:

...
Internal : "..../VTK/sonicLiquidPipeWithSymCylinDamperP100_mseq_poly_2. vtk"
Original cells:11079 points:60203 Additional cells:229990 additional points:11066

Combined patches : "..../VTK/allPatches/allPatches_2.vtk"
Combining patches:
patch 0 F1
patch 1 F2
patch 2 F3
patch 3 Walls
End

The boundaries (I load file <allpatches_2.vtk>) are shown correctly. It shows as


but the internal cells are not right, I think. It shows as


I'd like to shown it as real polyhedron without additional cells. How to get it?
__________________
rdu
------------------
Martin/Run Du
chnrdu is offline   Reply With Quote

Old   June 9, 2008, 03:18
Default As you correctly note, since V
  #2
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 18
olesen will become famous soon enough
As you correctly note, since VTK does not support polyhedral cells directly, foamToVTK tet decomposes polyhedral cells. If you really don't want to see the decomposed cells, you can try with converting to the EnSight format: the newer paraview versions (eg, 3.3 from CVS) manage to read the nfaced (polyhedral) cells and uses a convex point set to display them. This is generally slower than the tet decomposed equivalent and may fail when your cells have concavity. NOTE: the older paraview/VTK versions will segfault if you try to feed them an EnSight file with nfaced cells.
olesen is offline   Reply With Quote

Old   June 9, 2008, 22:29
Default Hi Mark, One question... is
  #3
Senior Member
 
Alexandre Pereira
Join Date: Mar 2009
Posts: 155
Rep Power: 8
alexandrepereira is on a distinguished road
Hi Mark,

One question... is there any way to import EnSight format into Paraview without the *%&"#)? "core dump"...?

Should i import it in "Little Endian" or "Big Endian" Format...?

As far i I can recall, Paraview 2.4.4. won't allow it,... Paraview 3.3.0 won't either...downloaded it today from kitware.com...

is there some kind of "trick" to be used in foamToEnsight? or it is just the way things are, not being able to use Paraview with EnSight Format...?

If I use OpenFOAM database format, or VTK export format, lots of streamlines get broken... but if i export it in EnSight Format it gets allright...

... the problem is that my EnSight demo license is almost expiring...

So I am looking for Paraview as one possible alternative...


Best Regards

Alex

Best regards

Alex
alexandrepereira is offline   Reply With Quote

Old   June 10, 2008, 03:43
Default Alex, I've currently broken
  #4
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 18
olesen will become famous soon enough
Alex,

I've currently broken my paraview 3.3 installation, but it did work with polyhedral without a core dump.
If you are using Linux, then the reader will initially read okay, but have the incorrect endian in the GUI. When you click Apply it will cause problems.

Another issue that you may have is that the current VTK EnSight reader incorrectly assumes that all the filenames/times are each on a single line (in the .case file). I submitted a patch, but haven't followed up on it.
olesen is offline   Reply With Quote

Old   June 10, 2008, 13:08
Default Hi Mark, Thanks for this in
  #5
Senior Member
 
Alexandre Pereira
Join Date: Mar 2009
Posts: 155
Rep Power: 8
alexandrepereira is on a distinguished road
Hi Mark,

Thanks for this info.

Alex
alexandrepereira 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
polyhedral cells bill Main CFD Forum 10 March 11, 2014 07:52
Polyhedral Mesh Diego Flores FLUENT 15 November 29, 2011 18:30
How to convert cells to polyhedral Sri FLUENT 1 September 12, 2007 07:43
Import error Gambit msh file with Cell Type 7 polyhedral cells philippose OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 June 1, 2007 09:14
polyhedral mesh guang ai CD-adapco 3 May 28, 2006 01:02


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