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

[snappyHexMesh] SnappyHexMesh and MultiRegion - get regions from Salome to OF

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 1, 2014, 11:07
Default SnappyHexMesh and MultiRegion - get regions from Salome to OF
  #1
dzi
Member
 
Join Date: Nov 2011
Location: Berlin
Posts: 31
Rep Power: 15
dzi is on a distinguished road
Hello Foamers

I tried many things and read a lot of snappyHexMesh threads but still I do not get it working to define regions in Salome and get them as patches in OF.
Mesh preparation in Salome:
Construct 2 bodies, here cylinder1 and cylinder2. The cylinders have a intersection.

I exploded the bodies into faces and meshed them with tetras.

I want to define regions at this state to find them in my case structure and set them as boundary.
Here cylinder2 is exploded to faces, the highlighted submesh is a face called cyl2botaussen.stl to be found as a patch.

After meshing, all faces are written out as cyl1*.stl
(total 3:
Code:
  cyl1bot.stl cyl1top.stl  cyl1wall.stl
) merged into one file
Code:
 cyl1.stl

and cyl2*.stl
(total 4:
Code:
cyl2botaussen.stl  cyl2botinner.stl cyl2top.stl  cyl2wall.stl
)
merged into
Code:
cyl2.stl
.

cyl2botinner and cyl1top are the interface between the cylinders, the mesh is matching/the same.

The salome-hdf and the resulting stl file one can be found here:
https://www.dropbox.com/sh/bkalhcqwv...pEOPQ2fz9KAMva

SHM:
For the next steps, I use this standard case structure. Meshing2 contains the mesh preparation folders. cyl1.stl and cyl2.stl are copied to /constant/triSurfaces.

The snappyHexmesh is invoked by the steps (=makeMesh.sh):

ideasUnvToFoam outbox.unv
surfaceFeatureExtract
snappyHexMesh -overwrite
splitMeshRegions -cellZones -overwrite
paraFoam -touchAll

After snappyhexmesh I get the expected 3 regions domain:
cylinder1, cylinder2, domain0
and the default patches inbetween.

My snappyHexmeshDict has the sections:

Code:
geometry
 {
     cyl1.stl
     {
         type triSurfaceMesh;
         name cylinder1;
     
     }
 

     cyl2.stl
     {
         type triSurfaceMesh;
     name cylinder2;
     regions
     {
         cyl2botaussen        { name cylinder2aussen; }
     }
     }
 };
 (...)
    features
     (
         {
             file "cyl1.eMesh";
             level 2;
         }
         {
             file "cyl2.eMesh";
             level 2;
         }      
     );
  (...)
 refinementSurfaces
     {
         cylinder1
         {
             // Surface-wise min and max refinement level
             level (2 3);
             faceZone cylinder1;
             cellZone cylinder1;
             cellZoneInside inside;
         }
          
         //cyl2.stl
     cylinder2
         {
           level (2 3);
           faceZone cylinder2;
           cellZone cylinder2;
           cellZoneInside inside;
 //        regions     -- must this be here?
 //          {
 //          cylinder2aussen { level (3 3);  }
 //          }
     }
     }
  (...)

I want to see „cylinder2aussen“ in paraview as a patch.

What is going wrong, why my patch „cylinder2aussen“ is lost after SHM and and how can I get it? (Do I need putting the desired patch in a separate stl-file? Different entries in the snappyHexMeshDict? In SurfaceFeatureExtractDict?) None of the tutorials and threads could really help me.

Maybe somebody can give me the missing link. Thanks for reading and support
dirk


All files and the complete case are here:
https://www.dropbox.com/s/awk8noaanm...ppyMRFEasy.tgz
(on Salome 7.2.0 and OF 2.2.2)





Attached Images
File Type: jpg salome1.jpg (16.9 KB, 152 views)
File Type: jpg salome2.jpg (50.5 KB, 148 views)
File Type: jpg salome3.jpg (41.9 KB, 169 views)
File Type: jpg pv_aftershm.jpg (72.1 KB, 175 views)
dzi is offline   Reply With Quote

Old   July 30, 2014, 06:39
Default
  #2
dzi
Member
 
Join Date: Nov 2011
Location: Berlin
Posts: 31
Rep Power: 15
dzi is on a distinguished road
... FYI there was an answer on another thread http://www.cfd-online.com/Forums/openfoam-meshing/135307-splitmeshregion.html#post503707
using
splitMeshRegions -cellZonesOnly -useFaceZones
but this issue still not solved for me yet.
dzi is offline   Reply With Quote

Old   September 4, 2014, 11:04
Default
  #3
dzi
Member
 
Join Date: Nov 2011
Location: Berlin
Posts: 31
Rep Power: 15
dzi is on a distinguished road
fyi, latest updates are in that thread:
http://www.cfd-online.com/Forums/openfoam-meshing/135307-splitmeshregion.html#post509069
dzi is offline   Reply With Quote

Reply

Tags
multiregion, salome, snappyhexmesh

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
[snappyHexMesh] Multiregion snappyhexmesh obiscolly50 OpenFOAM Meshing & Mesh Conversion 0 April 18, 2019 01:51
[snappyHexMesh] Unzoned cells by SnappyHexMesh for a Multiregion Mesh mujika OpenFOAM Meshing & Mesh Conversion 5 February 14, 2018 15:38
[snappyHexMesh] combine blockMesh features cyclicGgi and mixingPlane with multiregion snappyHexMesh AmRo OpenFOAM Meshing & Mesh Conversion 5 June 8, 2016 07:36
snappyHexMesh axisymmetric multiRegion Henning86 OpenFOAM Running, Solving & CFD 0 October 23, 2014 13:05
Simulation of SnappyHexMesh multiregion nithishgupta OpenFOAM Running, Solving & CFD 0 June 24, 2014 06:31


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