CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [Salome] Spherical slice in Salome for OpenFOAM (https://www.cfd-online.com/Forums/openfoam-meshing/127954-spherical-slice-salome-openfoam.html)

Jakob1 December 27, 2013 10:58

Spherical slice in Salome for OpenFOAM
 
I am trying to generate and mesh a slice of a spherical structure in Salome for an OpenFOAM simulation. I kinda managed to set it up, but I must have made some mistakes somewhere, since the converted unv produces some funky results.

What I did:
- Create symmetry face and revolute ( and explode)
- create NETGEN 1D2D submesh on symmetry face
- use mesh revolution to create prisms from submesh that match the geometry
- split pyramid into tetrahedra (pyramids cannot be exported through unv)
- create groups through filters (that got kinda messy)
- export unv and convert with ideasUnvToFoam

checkMesh gave me a whole bunch of errors:

Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.2.x                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.2.x-3addb4f901b3
Exec  : checkMesh
Date  : Dec 27 2013
Time  : 15:25:48
Host  : "jakob-SVS1312C5E"
PID    : 4078
Case  : /home/jakob/OpenFOAM/jakob-2.2.x/run/slice
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:          8210
    internal points:  1
    faces:            28044
    internal faces:  11854
    cells:            7980
    faces per cell:  4.99975
    boundary patches: 5
    point zones:      0
    face zones:      0
    cell zones:      0

Overall number of cells of each type:
    hexahedra:    0
    prisms:        7978
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    2
    polyhedra:    0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
 ***Unused points found in the mesh, number unused by faces: 1 number unused by cells: 1
  <<Writing 1 unused points to set unusedPoints
    Upper triangular ordering OK.
    Face vertices OK.
  *Number of regions: 2
    The mesh has multiple regions which are not connected by any face.
  <<Writing region information to "0/cellToRegion"

Checking patch topology for multiply connected surfaces...
    Patch              Faces    Points  Surface topology                 
    back                7979    4105    ok (non-closed singly connected) 
    top                34      69      ok (non-closed singly connected) 
    bottom              113      228      ok (non-closed singly connected) 
    front              7979    4105    ok (non-closed singly connected) 
    defaultFaces        85      167      multiply connected (shared edge) 
  <<Writing 2 conflicting points to set nonManifoldPoints

Checking geometry...
    Overall domain bounding box (-4.40873e-15 -24 -8.88178e-16) (23.9994 16.9706 4.16746)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
 ***Boundary openness (2.33305e-05 2.36133e-05 2.04112e-06) possible hole in boundary description.
 ***Open cells found, max cell openness: 0.573014, number of open cells 1
  <<Writing 1 non closed cells to set nonClosedCells
  <<Writing 3 cells with high aspect ratio to set highAspectRatioCells
    Minimum face area = 2.34156e-18. Maximum face area = 2.09338.  Face area magnitudes OK.
    Min volume = 1.33333e-300. Max volume = 0.453491.  Total volume = 1507.27.  Cell volumes OK.
    Mesh non-orthogonality Max: 147.38 average: 5.11099
 ***Number of non-orthogonality errors: 2.
  <<Writing 2 non-orthogonal faces to set nonOrthoFaces
 ***Error in face pyramids: 9 faces are incorrectly oriented.
  <<Writing 8 faces with incorrect orientation to set wrongOrientedFaces
    Max skewness = 0.842957 OK.
    Coupled point location match (average 0) OK.

Failed 5 mesh checks.

End

And in the OF case some defaultFaces appear in the boundary file.

Does anyone know what is going on and how to fix it? Like, what do I do with the defaultFaces (set to empty and ignore?) and where does that crazy local non-orthogonality originate from? What are wrongly oriented faces?

Would be very grateful for any help, this is how far my limited experience takes me :)

Heres the study: http://www.file-upload.net/download-...lice2.hdf.html


Oh and happy holidays !

Jakob1 January 10, 2014 04:51

I have still to solve my problem =/ Has noone encountered this before?

Can anyone think of a work-around? Like are there other, more appropriate meshing mechanisms?

Could really use some help!

EDIT: Solved. If anyone faces a similiar issue: There is a very useful Python script for Salome to OpenFOAM conversion: http://www.salome-platform.org/forum/forum_12/23863165


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