[ec91]

Let’s suppose that I have a volScalarField that must be zero everywhere except for cells adjacent to wall patches, where it is expected to be 1 (only near walls, not other types of patches).

Is there a way to loop over cells and identify those being next to a wall?

Thank you.

[clapointe]

Yes there is -- I think in the past I've done something similar by looping over patches, checking if the patches are walls, and if so looping over patch face cells. It would look something like :

    //Patch list
    const fvPatchList& patches = mesh.boundary();

    //Patch loop
    forAll(patches, patchi)
    {
        //Current patch
        const fvPatch& currPatch = patches[patchi];

        //Is wall?
        if (mesh.boundary()[patchi].type() == "wall")
        {
            //Patch face cells
            forAll(currPatch,facei)
            {
                    //Current face cell
                    label faceCelli = currPatch.faceCells()[facei];
                    //Do something here with faceCelli
            }
        }
    }

Hopefully this will be enough to get you started.

Caelan

[olesen]

Originally Posted by clapointe

Yes there is -- I think in the past I've done something similar by looping over patches, checking if the patches are walls, and if so looping over patch face cells. It would look something like :

    //Patch list
    const fvPatchList& patches = mesh.boundary();

    //Patch loop
    forAll(patches, patchi)
    {
        //Current patch
        const fvPatch& currPatch = patches[patchi];

        //Is wall?
        if (mesh.boundary()[patchi].type() == "wall")
        {
            //Patch face cells
            forAll(currPatch,facei)
            {
                    //Current face cell
                    label faceCelli = currPatch.faceCells()[facei];
                    //Do something here with faceCelli
            }
        }
    }

Hopefully this will be enough to get you started.

Caelan

Don't use literal strings to check the boundary condition type!! You want isA<wallPolyPatch> which checks based on the class inheritance.

[clapointe]

Absolutely -- good catch. This was old code and should be updated if used.

Caelan