CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Paraview & paraFoam (http://www.cfd-online.com/Forums/openfoam-paraview/)
-   -   Visualizing checkMesh results in Paraview (http://www.cfd-online.com/Forums/openfoam-paraview/119813-visualizing-checkmesh-results-paraview.html)

Claudio June 25, 2013 12:58

Visualizing checkMesh results in Paraview
 
Hi there,

I'm working on a case where the grid is built using enGrid. It all seems fine, but when I run CheckMesh (OpenFOAM) it fails 4 checks, and writes the culprits to 4 different cellSets.
How do I display these in Paraview?

I tried copy the entire constant folder over, with the sets subfolder in it, but no luck.

Thanks for any help.

Claudio

Code:

Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:          900745
    internal points:  900745
    faces:            9630792
    internal faces:  9630792
    cells:            4748020
    boundary patches: 0
    point zones:      0
    face zones:      0
    cell zones:      0

Overall number of cells of each type:
    hexahedra:    0
    prisms:        269504
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    4478516
    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                 

Checking geometry...
    Overall domain bounding box (-15 -10 -10) (25 10 10)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (0 0 0) OK.
 ***Open cells found, max cell openness: 1, number of open cells 21134
  <<Writing 21134 non closed cells to set nonClosedCells
  <<Writing 245 cells with high aspect ratio to set highAspectRatioCells
    Minumum face area = 1.27259977e-07. Maximum face area = 10.29648759.  Face area magnitudes OK.
    Min volume = 1.666666667e-300. Max volume = 9.232429517.  Total volume = 31465.12518.  Cell volumes OK.
    Mesh non-orthogonality Max: 179.7735569 average: 15.80337767
  *Number of severely non-orthogonal faces: 8959.
 ***Number of non-orthogonality errors: 20686.
  <<Writing 29645 non-orthogonal faces to set nonOrthoFaces
 ***Error in face pyramids: 21214 faces are incorrectly oriented.
  <<Writing 21198 faces with incorrect orientation to set wrongOrientedFaces
 ***Max skewness = 216.0434447, 549 highly skew faces detected which may impair the quality of the results
  <<Writing 549 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Failed 4 mesh checks.

End


tomf June 26, 2013 04:43

Hi Claudio,

You can use the foamToVTK utility to make VTK files of any set:

Code:

foamToVTK -faceSet nonOrthoFaces -time 0
You can than open your case in ParaView and open the vtk files for each set. Just check the names of the set from the checkMesh log and run foamToVTK for each faceSet. Do not forget the -time 0 or otherwise it may overwrite your old VTK files.

Regards,
Tom


All times are GMT -4. The time now is 10:36.