/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.0                                   |
|   \\  /    A nd           | Web:      http://www.openfoam.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version         2.0;
    format          ascii;
    class           dictionary;
    object          snappyHexMeshDict;
}

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

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

geometry {
    clarkypt_le.stl {
        type triSurfaceMesh;
        name leadingEdge;
    }

    clarkypt_body.stl {
        type triSurfaceMesh;
        name foil;
    }

    clarkypt_te.stl {
        type triSurfaceMesh;
        name trailingEdge;
    }

    innerBox {
    	type searchableBox;
    	min (-0.12 -0.12 -0.07);
    	max (2 0.12 0.19);
    }
};



// Settings for the castellatedMesh generation.
castellatedMeshControls {
	maxLocalCells 1000000;
	maxGlobalCells 2000000;
	minRefinementCells 0;
	nCellsBetweenLevels 3;
	refinementSurfaces {
		foil {
			level (6 6);
		}

		leadingEdge {
			level (7 7);
		}

		trailingEdge {
			level (7 7);
		}
	}
	refinementRegions {
		innerBox {
			mode inside;
			levels ((5 5));
		}
	}
	resolveFeatureAngle 60;
	locationInMesh (-1 0 0);
	features ();
	allowFreeStandingZoneFaces false;
}



// Settings for the snapping.
snapControls {
	nSmoothPatch 6;
	tolerance 10;
	nSolveIter 60;
	nRelaxIter 10;
}



// Settings for the layer addition.
addLayersControls {
	relativeSizes false;
	layers{
		leadingEdge_OBJECT {
			nSurfaceLayers 5;
		}

		foil_OBJECT {
			nSurfaceLayers 5;
		}

		trailingEdge_OBJECT {
			nSurfaceLayers 5;
		}
	}
	expansionRatio 1.2;
	finalLayerThickness 0.0001;
	minThickness 0.0000311;
	nGrow 0;
	featureAngle 180;
	nRelaxIter 5;
	nSmoothSurfaceNormals 1;
	nSmoothNormals 3;
	nSmoothThickness 10;
	maxFaceThicknessRatio 0.5;
	maxThicknessToMedialRatio 0.3;
	minMedianAxisAngle 130;
	nBufferCellsNoExtrude 0;
	nLayerIter 50;
}



// Generic mesh quality settings. At any undoable phase these determine
// where to undo.
meshQualityControls {
	maxNonOrtho 60;
	maxBoundarySkewness 20;
	maxInternalSkewness 4;
	maxConcave 80;
	minFlatness 0.5;
	minVol -1e13;
	minArea -1;
	minTwist 0.05;
	minDeterminant 0.001;
	minFaceWeight 0.05;
	minVolRatio 0.01;
	minTriangleTwist -1;
	nSmoothScale 4;
	errorReduction 0.75;
	minTetQuality 0.01;
}


// Advanced

// Flags for optional output
// 0 : only write final meshes
// 1 : write intermediate meshes
// 2 : write volScalarField with cellLevel for postprocessing
// 4 : write current intersections as .obj files
debug 0;


// Merge tolerance. Is fraction of overall bounding box of initial mesh.
// Note: the write tolerance needs to be higher than this.
mergeTolerance 1E-6;


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