student666 |
July 15, 2014 08:33 |
defining two cellZone in SHM
Hi all,
I'm trying to define two cellzone in SHMdict: one is rotating, the other is porous zone.
here the SHM dict.
Code:
/*--------------------------------*- 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
|