/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    object      extrudeToRegionMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Name of the region to create
region          liquidFilm;

// Specification of faces to extrude. Either faceZones (either exclusively
// internal faces or boundary faces) or faceSets (boundary faces only). These
// faces will be connected to the bottom patch of the extruded region. If
// shadow faces are also specified then these will be connected to the top
// patch of the extruded region.

// The faceZones to extrude
faceZones       (f0 f1);
// faceZonesShadow (f0Shadow f1Shadow);

// The faceSets to extrude
// faceSets        (f0 f1);
// faceSetsShadow  (f0Shadow f1Shadow);

// Set intrude to yes/true to extrude into/overlapping the primary region mesh
// i.e. intrude rather than extrude
// Used to create film regions which overlap the primary region mesh
// intrude         yes;

// Adapt the original mesh to have mapped patches that connect to the extruded
// region? If so, then extruded internal faces become baffles with mapped
// patches, and extruded boundary faces are repatched onto mapped patches.
// Otherwise the original mesh is not modified. The extruded mesh will still
// have mapped patches created which will need to be changed.
adaptMesh       true;

//- Extrusion settings
extrudeModel    linearNormal;
// extrudeModel   linearDirection;
// extrudeModel   wedge;
// extrudeModel   linearRadial;
// extrudeModel   radial;
// extrudeModel   sigmaRadial;

nLayers         10;

expansionRatio  0.9;

linearNormalCoeffs
{
    thickness       0.05;
}

linearDirectionCoeffs
{
    direction       (0 1 0);
    thickness       0.05;
}

wedgeCoeffs
{
    axisPt      (0 0.1 -0.05);
    axis        (-1 0 0);
    angle       360;  // For nLayers=1 assume symmetry so angle/2 on each side
}

linearRadialCoeffs
{
    R              0.1;
    // Optional inner radius
    Rsurface        0.01;
}

radialCoeffs
{
    // Radii specified through interpolation table
    R               table ((0 0.01)(3 0.03)(10 0.1));
}

sigmaRadialCoeffs
{
    RTbyg           1;
    pRef            1;
    pStrat          1;
}


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