[mediode]

Hello,

I would appreciate if anybody could give me a hint how to determine the area of an existing patch? (I would like to investigate the geometrical error produced by using snappyhexmesh.)

m*

[gwierink]

Hi mediode,

If you know the name of the patch (or the patchID straight away), then you can do the following. Here, the name of the patch is "movingWall", from the icoFoam cavity case.

// 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];
}

[mediode]

Thanks a lot for the answer. Sorry for the delay.
Best regards
M*