[sailor79]

Hi folks,

I want to access the faces of a certain patch in a certain region. I do some first tests with the code from the Thread area of a patch

// Find the patchID of the patch by name
label patchID = mesh.boundaryMesh().findPatchID("movingWall");

// Create a polyPatch for looping
const polyPatch& myPatch = mesh.boundaryMesh()[patchID];

// Initialize patchArea
scalar patchArea = 0.0;

// Loop trhough all faces on the polyPatch, adding their magnitude surface
// area vectors
forAll(myPatch, faceI)
{
    patchArea += mesh.magSf().boundaryField()[patchID][faceI];
}

But as I am working with chtMultiRegionFoam and the desired patch is an interface patch created by splitMeshRegions (region1_to_region2) I am not able to access the faces. The patchID becomes "-1" as the patch is defined in the boundary file within the regions-directory.

How do I have to modify the mentioned code?

Thanks for your support

Michael

[sailor79]

I figured out how to access the mesh of a region. Maybe the following information will help other foamers:

You just have to add these lines to the main program:

    timeSelector::addOptions();
#   include "addRegionOption.H"

and replace

#   include "createMesh.H"

by

#   include "createNamedMesh.H"

Michael