[student666]

Hi all,

I'm trying to define two cellzone in SHMdict: one is rotating, the other is porous zone.

here the SHM dict.

   /*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      snappyHexMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Which of the steps to run
castellatedMesh true;    // make basic mesh ?
snap            true;    // decide to snap back to surface ?
addLayers       true;   // decide to add viscous layers ?


geometry // Load in STL files here
{
    inlet.stl {type triSurfaceMesh; name inlet;}
    outlet.stl {type triSurfaceMesh; name outlet;}
    wallComp.stl {type triSurfaceMesh; name wallComp;}
    wallAsp.stl {type triSurfaceMesh; name wallAsp;}
    boccaglio.stl {type triSurfaceMesh; name boccaglio;}
    blade.stl {type triSurfaceMesh; name blade;}
    wallMan.stl {type triSurfaceMesh; name wallMan;}
    wallExt.stl {type triSurfaceMesh; name wallExt;}
    tubi.stl {type triSurfaceMesh; name tubi;}
    dominio.stl {type triSurfaceMesh; name volume;}

    porous {type searchableBox; min (-0.235 -0.3922 -0.087); max ( -0.165 -0.2072 0.0875);name porous}
    rotating {type searchableCylinder; point1 (0. -0.0072 0.); point2 ( 0. -0.0279 0.); radius 0.078; name rotating;}
};

castellatedMeshControls
{
    maxLocalCells 1500000;  //max cells per CPU core
    maxGlobalCells 3000000; //max cells to use before mesh deletion step
    minRefinementCells 10;  //was 0 - zero means no bad cells are allowed during refinement stages
    maxLoadUnbalance 0.10;
    nCellsBetweenLevels 1;  // expansion factor between each high & low refinement zone

    // Explicit feature edge refinement
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

    features // taken from STL from each .eMesh file created by "SurfaceFeatureExtract" command
    (
        {file "inlet.eMesh"; level 0;} 
        {file "outlet.eMesh"; level 0;}
        {file "wallComp.eMesh"; level 3;}
        {file "wallAsp.eMesh"; level 2;} 
        {file "boccaglio.eMesh"; level 5;}
        {file "wallMan.eMesh"; level 2;} 
        {file "blade.eMesh"; level 5;}
    {file "tubi.eMesh"; level 5;}
    {file "wallExt.eMesh"; level 0;}
                
    );

    // Surface based refinement
    // ~~~~~~~~~~~~~~~~~~~~~~~~

    refinementSurfaces // Surface-wise min and max refinement level
    {
    inlet {level (0 0);}
        outlet {level (0 0);}
    wallComp {level (3 3);}
    wallAsp {level (3 3);}
    boccaglio {level (5 5);}
    wallMan {level (3 3);}
    blade {level (5 5);}
    tubi {level (5 5);}
    wallExt {level (0 0);}
    
        
    rotating {level (5 5); regions{}; faceZone MRF; cellZone MRF; zoneInside true;}
    porous {level (5 5); regions{}; faceZone porous; cellZone pourous; zoneInside true;}
    }  


    resolveFeatureAngle 60;  //was80 Resolve sharp angles // Default 30
    refinementRegions        // In descending levels of fine-ness
    {volume {mode distance; levels ((0.0008 3) (0.002 2) (0.01 1));}}//((0.0008 4) (0.002 3) (0.01 2));}} // was ((0.001 4) (0.003 3) (0.01 2))
    locationInMesh (0. -0.05 0.);  //to decide which side of mesh to keep **
    allowFreeStandingZoneFaces true;
}

ECC ...

I successfully obtain a good mesh, but I only have one cellZone at all, the porous one, why?
I have a look deeper into forum, but I didn't find anything similar to my question.

Can someone help me please?

Regards

Michele