[farzadmech]

Dear Oliveira
Thanks for your piece of code, but it integrate over points and not cells. Do you know how can I use this code but integrate over cells instead of nodes?

Thanks,
Farzad

Originally Posted by t.oliveira

This is an updated version that works with OpenFOAM 5.x:

functions
{
    flux
    {
        functionObjectLibs ("libfieldFunctionObjects.so");
        type             surfaceFieldValue;

        enabled          true;
        writeControl     outputTime;
        writeInterval    1;
        writeFields      false;

        name myName;

        // regionType: patch/faceZone/sampledSurface
        regionType sampledSurface;

        sampledSurfaceDict
        {
            type                  plane;
            planeType             pointAndNormal;
            pointAndNormalDict
            {
                point  (0.1 0 0);
                normal (1 0 0);
            };
            interpolate false;
        }

        // Operation: areaAverage/sum/weightedAverage/...
        operation    areaNormalIntegrate;
        fields
        (
            U
        );
    }
}

Use it executing:

postProcess -fields "(U)" -dict system/integrateOnSurfaceDict

Best wishes,

Thomas