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

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Which of the steps to run
castellatedMesh true;
snap            true;
addLayers       true;

geometry
{
    surface.stl
    {
        type triSurfaceMesh;
        name surface;
    }
    PcRegion
    {
        type   cylinder;
        point1 ( 0  0  0 );
        point2 ( 0  0  0.02731 );
        radius 0.009;
    }
    nozzleRegion
    {
        type   cylinder;
        point1 ( 0  0  0.02731 );
        point2 ( 0  0  0.03731 );
        radius 0.008;
    }
}


// Settings for the castellatedMesh generation.
castellatedMeshControls
{

    maxLocalCells 100000;

    maxGlobalCells 5000000;

    minRefinementCells 10;

    nCellsBetweenLevels 3;//require 3

    features
    (
        {
            file "surface.eMesh";
            level 3;
       }
    );

    refinementSurfaces
    {
        surface
        {
            level (2 3);
        }
    }

    // Resolve sharp angles
    resolveFeatureAngle 30;//require 30

    refinementRegions
    {
        
        PcRegion
        {
            mode inside;
            levels ((0.0 3));
        }
        nozzleRegion
        {
            mode inside;
            levels ((0.0 3));
        }
    }

    locationInMesh (0 0 0.01);

    allowFreeStandingZoneFaces false;
}



// Settings for the snapping.
snapControls
{

    nSmoothPatch 5;

    tolerance 2;

    nSolveIter 30;

    nRelaxIter 5;

    nFeatureSnapIter 20;

    implicitFeatureSnap false;

    explicitFeatureSnap true;

    multiRegionFeatureSnap false;
}



// Settings for the layer addition.
addLayersControls
{

    relativeSizes true;

    layers
    {
    }

    expansionRatio 1.0;

    finalLayerThickness 0.3;

    minThickness 0.1;

    nGrow 0;

    featureAngle 30;

    nRelaxIter 3;

    nSmoothSurfaceNormals 1;

    nSmoothNormals 3;

    nSmoothThickness 10;

    maxFaceThicknessRatio 0.5;

    maxThicknessToMedialRatio 0.3;

    minMedialAxisAngle 90;

    nBufferCellsNoExtrude 0;

    nLayerIter 50;

    nRelaxedIter 20;
}

meshQualityControls
{

    maxNonOrtho 65;


    maxBoundarySkewness 20;
    maxInternalSkewness 4;

    maxConcave 80;

    minVol 1e-20;

    minTetQuality 1e-30;

    minArea -1;

    minTwist 0.02;

    minDeterminant 0.001;

    minFaceWeight 0.02;

    minVolRatio 0.01;

    minTriangleTwist -1;

    nSmoothScale 4;

    errorReduction 0.75;
}

mergeTolerance 1e-6;


// ************************************************************************* //