[decibelle]

Hello everybody,

I would like to understand how surfaceToCell works in topoSet for a cellSet source.
surfaceToCell is composed of of this arguments:
- surface
- outsidepoint
- cut
- inside
- outside
- near
- curvature

I don't understand the meaning of these arguments except "surface".

Someone can help me?

Thanks in advance

[alexeym]

Hi,

$ src
$ pwd
$HOME/OpenFOAM/OpenFOAM-4.1/src
$ find . -name 'surfaceToCell.*'
./meshTools/sets/cellSources/surfaceToCell/surfaceToCell.C
./meshTools/sets/cellSources/surfaceToCell/surfaceToCell.H
./meshTools/lnInclude/surfaceToCell.C
./meshTools/lnInclude/surfaceToCell.H
$ less ./meshTools/sets/cellSources/surfaceToCell/surfaceToCell.H

If you scroll a little bit, you find these comments:

        //- Name of surface file
        const fileName surfName_;

        //- Points which are outside
        const pointField outsidePoints_;

        //- Include cut cells
        const bool includeCut_;

        //- Include inside cells
        const bool includeInside_;

        //- Include outside cells
        const bool includeOutside_;

        //- Determine inside/outside purely using geometric test
        //  (does not allow includeCut)
        const bool useSurfaceOrientation_;

        //- If > 0 : include cells with distance from cellCentre to surface
        //  less than nearDist.
        const scalar nearDist_;

        //- If > -1 : include cells with normals at nearest surface points
        //  varying more than curvature_.
        const scalar curvature_;

Correspondence between variable names and dictionary arguments can be found in ./meshTools/sets/cellSources/surfaceToCell/surfaceToCell.C (though they are quite obvious).

[decibelle]

Originally Posted by alexeym

Hi,

$ src
$ pwd
$HOME/OpenFOAM/OpenFOAM-4.1/src
$ find . -name 'surfaceToCell.*'
./meshTools/sets/cellSources/surfaceToCell/surfaceToCell.C
./meshTools/sets/cellSources/surfaceToCell/surfaceToCell.H
./meshTools/lnInclude/surfaceToCell.C
./meshTools/lnInclude/surfaceToCell.H
$ less ./meshTools/sets/cellSources/surfaceToCell/surfaceToCell.H

If you scroll a little bit, you find these comments:

        //- Name of surface file
        const fileName surfName_;

        //- Points which are outside
        const pointField outsidePoints_;

        //- Include cut cells
        const bool includeCut_;

        //- Include inside cells
        const bool includeInside_;

        //- Include outside cells
        const bool includeOutside_;

        //- Determine inside/outside purely using geometric test
        //  (does not allow includeCut)
        const bool useSurfaceOrientation_;

        //- If > 0 : include cells with distance from cellCentre to surface
        //  less than nearDist.
        const scalar nearDist_;

        //- If > -1 : include cells with normals at nearest surface points
        //  varying more than curvature_.
        const scalar curvature_;

Correspondence between variable names and dictionary arguments can be found in ./meshTools/sets/cellSources/surfaceToCell/surfaceToCell.C (though they are quite obvious).

Thank you for your help.

but it is not yet clear for me.

If we take for example the motorbike tutorials, the outsidepoints are outside the moto but inside the domain, outside the moto and outside the domain or inside the moto?
About the cut cells, how it is possible to choose the cells cut by motorbike surface?
Inside, I think it's inside the motorbike but, in this case there are not cells inside.
Outside I think it's outside the motorbike but in the domain.
near I understand
and I don't understand the function of curvature.

[alexeym]

Hi,

Now it is not clear to me, what you are asking. You started from cellSet source in topoSet, i.e. you were selecting CELLS. Now you talk about points (though maybe it was my mistake to include this property in example).

You have surface. You can select cells of a computation domain, that are inside this surface, outside this surface, and cut by this surface. To define meaning of "inside" and "outside", you can use either reference point, or surface orientation.

Curvature works only for nearDist. And it means the following:

            Info<< "    Selecting cells with cellCentre closer than "
                << nearDist_ << " to surface and curvature factor"
                << " less than " << curvature_ << endl;

You can go further into source, to see exact mathematical expression.

[decibelle]

in my previous message I wanted to apply surfaceToCell in the motorbike case to understand. But I'm agree with you it wasn't clear. Do you know this tutorial?

[alexeym]

It is MUCH easier to understand if you really try to apply this topoSet to your problem and then post your cognitive dissonances.

I know the tutorial.

[HPE]

For future reference, the following tutorials in v1912 contain examples of `surfaceToCell`:

mesh/foamyHexMesh/flange/system/topoSetDict-background
multiphase/interFoam/laminar/waves/mangroveInteraction/system/setFieldsDict

    surfaceToCell
    {
        file            "./constant/triSurface/seaweed.stl";
        outsidePoints   ((2 0.25 0.15));
        includeCut      true;
        includeInside   true;
        includeOutside  false;
        nearDistance    -1;

        curvature       -100;

        fieldValues
        (
            volScalarFieldValue MangrovesIndex 1
        );
    }

[ofWorld]

Hi Alexey,

I extruded a 0.5 mm thick solid shell from the external surface of a fluid mesh region (see attached shell image).

The external surface is named as "fluidToSolid", and in the extruded solid, both the inner and outer walls are named by "fluidToSolid".

I'd like to separate the two surfaces (inner wall vs. outer wall) using topoSet and then rename them using createPatch.

No matter how I played with the settings in topoSetDict, I cannot have expected results. Below are my codes. Do you have any idea?

Thanks a lot for your help in advance.
Guanqi

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

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

actions
(
    {
        name    inner-cells;
        type    cellSet;
        action  new;
        source  surfaceToCell;
        sourceInfo
        {
	    file "constant/geometry/fluidToSolid.stl";
	    outsidePoints (
(0.151308 0.00186969 0.208)
);
            includeCut false;
            includeInside true;
            includeOutside false;
            nearDistance -1; // cells with centre near surf
                             // (set to -1 if not used)
            curvature -100; // cells within nearDistance
                            // and near surf curvature
                            // (set to -100 if not used)
        }
    }
    {
        name    inner-wall;
        type    faceSet;
        action  new;
        source  cellToFace;

        sourceInfo
        {
            set inner-cells;
            option  all;
        }
    }
    {
        name    inner-wall;
        type    faceSet;
        action  subset;
        source  patchToFace;
        sourceInfo
        {
            name "fluidToSolid";
        }
    }
);

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

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       polyBoundaryMesh;
    location    "1e-05/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

5
(
    fluidGraftInlet
    {
        type            wall;
        inGroups        List<word> 1(wall);
        nFaces          32;
        startFace       17454;
    }
    fluidVeinInlet
    {
        type            wall;
        inGroups        List<word> 1(wall);
        nFaces          56;
        startFace       17486;
    }
    fluidVeinOutlet
    {
        type            wall;
        inGroups        List<word> 1(wall);
        nFaces          44;
        startFace       17542;
    }
    fluidToSolid
    {
        type            wall;
        inGroups        List<word> 1(wall);
        nFaces          17480;
        startFace       17586;
    }
    solid_to_fluid
    {
        type            patch;
        nFaces          0;
        startFace       35066;
    }
)

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