surfaceZonesInfo.C

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Copyright (C) 2013-2025 OpenFOAM Foundation
     \\/     M anipulation  |
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.
 
    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.
 
    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
\*---------------------------------------------------------------------------*/
 
#include "surfaceZonesInfo.H"
#include "searchableSurface.H"
#include "searchableSurfaceList.H"
#include "polyMesh.H"
#include "dictionary.H"
 
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
const Foam::NamedEnum<Foam::surfaceZonesInfo::areaSelectionAlgo, 4>
Foam::surfaceZonesInfo::areaSelectionAlgoNames
{
    "inside",
    "outside",
    "insidePoint",
    "none"
};
 
const Foam::NamedEnum<Foam::surfaceZonesInfo::faceZoneType, 3>
Foam::surfaceZonesInfo::faceZoneTypeNames
{
    "internal",
    "baffle",
    "boundary"
};
 
 
// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
Foam::surfaceZonesInfo::surfaceZonesInfo
(
    const searchableSurface& surface,
    const dictionary& surfacesDict
)
:
    faceZoneName_(),
    cellZoneName_(),
    zoneInside_(NONE),
    zoneInsidePoint_(point::min),
    faceType_(INTERNAL)
{
    // Global zone names per surface
    if (surfacesDict.readIfPresent("faceZone", faceZoneName_))
    {
        // Read optional entry to determine inside of faceZone
 
        word method;
        bool hasSide =
            surfacesDict.readIfPresent("mode", method)
         || surfacesDict.readIfPresent("cellZoneInside", method);
 
        if (hasSide)
        {
            zoneInside_ = areaSelectionAlgoNames[method];
            if (zoneInside_ == INSIDEPOINT)
            {
                zoneInsidePoint_ =
                    surfacesDict.lookup<point>("insidePoint", dimLength);
            }
        }
 
        // Read optional cellZone name
 
        if (surfacesDict.readIfPresent("cellZone", cellZoneName_))
        {
            if
            (
                (
                    zoneInside_ == INSIDE
                 || zoneInside_ == OUTSIDE
                )
            && !surface.hasVolumeType()
            )
            {
                IOWarningInFunction
                (
                    surfacesDict
                )   << "Illegal entry zoneInside "
                    << areaSelectionAlgoNames[zoneInside_]
                    << " for faceZone "
                    << faceZoneName_
                    << " since surface is not closed." << endl;
            }
        }
        else if (hasSide)
        {
            IOWarningInFunction
            (
                surfacesDict
            )   << "Unused entry zoneInside for faceZone "
                << faceZoneName_
                << " since no cellZone specified."
                << endl;
        }
 
        // How to handle faces on faceZone
        word faceTypeMethod;
        if (surfacesDict.readIfPresent("faceType", faceTypeMethod))
        {
            faceType_ = faceZoneTypeNames[faceTypeMethod];
        }
    }
}
 
 
Foam::surfaceZonesInfo::surfaceZonesInfo
(
    const word& faceZoneName,
    const word& cellZoneName,
    const areaSelectionAlgo& zoneInside,
    const point& zoneInsidePoint,
    const faceZoneType& faceType
)
:
    faceZoneName_(faceZoneName),
    cellZoneName_(cellZoneName),
    zoneInside_(zoneInside),
    zoneInsidePoint_(zoneInsidePoint),
    faceType_(faceType)
{}
 
 
Foam::surfaceZonesInfo::surfaceZonesInfo(const surfaceZonesInfo& surfZone)
:
    faceZoneName_(surfZone.faceZoneName()),
    cellZoneName_(surfZone.cellZoneName()),
    zoneInside_(surfZone.zoneInside()),
    zoneInsidePoint_(surfZone.zoneInsidePoint()),
    faceType_(surfZone.faceType())
{}
 
 
Foam::labelList Foam::surfaceZonesInfo::getUnnamedSurfaces
(
    const PtrList<surfaceZonesInfo>& surfList
)
{
    labelList anonymousSurfaces(surfList.size());
 
    label i = 0;
    forAll(surfList, surfi)
    {
        if (surfList[surfi].faceZoneName().empty())
        {
            anonymousSurfaces[i++] = surfi;
        }
    }
    anonymousSurfaces.setSize(i);
 
    return anonymousSurfaces;
}
 
 
Foam::labelList Foam::surfaceZonesInfo::getNamedSurfaces
(
    const PtrList<surfaceZonesInfo>& surfList
)
{
   labelList namedSurfaces(surfList.size());
 
    label namedI = 0;
    forAll(surfList, surfi)
    {
        if
        (
            surfList.set(surfi)
         && surfList[surfi].faceZoneName().size()
        )
        {
            namedSurfaces[namedI++] = surfi;
        }
    }
    namedSurfaces.setSize(namedI);
 
    return namedSurfaces;
}
 
 
Foam::labelList Foam::surfaceZonesInfo::getClosedNamedSurfaces
(
    const PtrList<surfaceZonesInfo>& surfList,
    const searchableSurfaceList& allGeometry,
    const labelList& surfaces
)
{
    labelList closed(surfList.size());
 
    label closedI = 0;
    forAll(surfList, surfi)
    {
        if
        (
            surfList.set(surfi)
         && surfList[surfi].cellZoneName().size()
         && (
                surfList[surfi].zoneInside() == surfaceZonesInfo::INSIDE
             || surfList[surfi].zoneInside() == surfaceZonesInfo::OUTSIDE
            )
         && allGeometry[surfaces[surfi]].hasVolumeType()
        )
        {
            closed[closedI++] = surfi;
        }
    }
    closed.setSize(closedI);
 
    return closed;
}
 
 
Foam::labelList Foam::surfaceZonesInfo::getUnclosedNamedSurfaces
(
    const PtrList<surfaceZonesInfo>& surfList,
    const searchableSurfaceList& allGeometry,
    const labelList& surfaces
)
{
    labelList unclosed(surfList.size());
 
    label unclosedI = 0;
    forAll(surfList, surfi)
    {
        if
        (
            surfList.set(surfi)
         && !allGeometry[surfaces[surfi]].hasVolumeType()
        )
        {
            unclosed[unclosedI++] = surfi;
        }
    }
    unclosed.setSize(unclosedI);
 
    return unclosed;
}
 
 
Foam::labelList Foam::surfaceZonesInfo::getAllClosedNamedSurfaces
(
    const PtrList<surfaceZonesInfo>& surfList,
    const searchableSurfaceList& allGeometry,
    const labelList& surfaces
)
{
    labelList closed(surfList.size());
 
    label closedI = 0;
    forAll(surfList, surfi)
    {
        if
        (
            surfList.set(surfi)
         && surfList[surfi].cellZoneName().size()
         && allGeometry[surfaces[surfi]].hasVolumeType()
        )
        {
            closed[closedI++] = surfi;
        }
    }
    closed.setSize(closedI);
 
    return closed;
}
 
 
Foam::labelList Foam::surfaceZonesInfo::getInsidePointNamedSurfaces
(
    const PtrList<surfaceZonesInfo>& surfList
)
{
    labelList closed(surfList.size());
 
    label closedI = 0;
    forAll(surfList, surfi)
    {
        if
        (
            surfList.set(surfi)
         && surfList[surfi].cellZoneName().size()
         && surfList[surfi].zoneInside() == surfaceZonesInfo::INSIDEPOINT
        )
        {
            closed[closedI++] = surfi;
        }
    }
    closed.setSize(closedI);
 
    return closed;
}
 
 
Foam::labelList Foam::surfaceZonesInfo::addCellZonesToMesh
(
    const PtrList<surfaceZonesInfo>& surfList,
    const labelList& namedSurfaces,
    polyMesh& mesh
)
{
    labelList surfaceToCellZone(surfList.size(), -1);
 
    cellZoneList& cellZones = mesh.cellZones();
 
    forAll(namedSurfaces, i)
    {
        label surfi = namedSurfaces[i];
 
        const word& cellZoneName = surfList[surfi].cellZoneName();
 
        if (cellZoneName != word::null)
        {
            label zonei = cellZones.findIndex(cellZoneName);
 
            if (zonei == -1)
            {
                zonei = cellZones.size();
                cellZones.setSize(zonei+1);
                cellZones.set
                (
                    zonei,
                    new cellZone
                    (
                        cellZoneName,   // name
                        labelList(0),   // addressing
                        cellZones       // cellZones
                    )
                );
            }
 
            surfaceToCellZone[surfi] = zonei;
        }
    }
 
    // Check they are synced
    List<wordList> allCellZones(Pstream::nProcs());
    allCellZones[Pstream::myProcNo()] = cellZones.toc();
    Pstream::gatherList(allCellZones);
    Pstream::scatterList(allCellZones);
 
    for (label proci = 1; proci < allCellZones.size(); proci++)
    {
        if (allCellZones[proci] != allCellZones[0])
        {
            FatalErrorInFunction
                << "Zones not synchronised among processors." << nl
                << " Processor0 has cellZones:" << allCellZones[0]
                << " , processor" << proci
                << " has cellZones:" << allCellZones[proci]
                << exit(FatalError);
        }
    }
 
    return surfaceToCellZone;
}
 
 
Foam::labelList Foam::surfaceZonesInfo::addFaceZonesToMesh
(
    const PtrList<surfaceZonesInfo>& surfList,
    const labelList& namedSurfaces,
    polyMesh& mesh
)
{
    labelList surfaceToFaceZone(surfList.size(), -1);
 
    faceZoneList& faceZones = mesh.faceZones();
 
    forAll(namedSurfaces, i)
    {
        label surfi = namedSurfaces[i];
 
        const word& faceZoneName = surfList[surfi].faceZoneName();
 
        label zonei = faceZones.findIndex(faceZoneName);
 
        if (zonei == -1)
        {
            zonei = faceZones.size();
            faceZones.setSize(zonei+1);
            faceZones.set
            (
                zonei,
                new faceZone
                (
                    faceZoneName,   // name
                    labelList(0),   // addressing
                    boolList(0),    // flipmap
                    faceZones       // faceZones
                )
            );
        }
 
        surfaceToFaceZone[surfi] = zonei;
    }
 
    // Check they are synced
    List<wordList> allFaceZones(Pstream::nProcs());
    allFaceZones[Pstream::myProcNo()] = faceZones.toc();
    Pstream::gatherList(allFaceZones);
    Pstream::scatterList(allFaceZones);
 
    for (label proci = 1; proci < allFaceZones.size(); proci++)
    {
        if (allFaceZones[proci] != allFaceZones[0])
        {
            FatalErrorInFunction
                << "Zones not synchronised among processors." << nl
                << " Processor0 has faceZones:" << allFaceZones[0]
                << " , processor" << proci
                << " has faceZones:" << allFaceZones[proci]
                << exit(FatalError);
        }
    }
 
    return surfaceToFaceZone;
}
 
 
// ************************************************************************* //