반갑습니다.
허나 구지 프로그래밍 코드를 이용하여 계산식을 짤 필요 없이 OpenFOAM에서 제공하는 forces, forceCoeffs 를 이용하는 것이 좋을 듯합니다. OKUC 혹은 Google에서도 정보를 찾아볼 수 있을 듯 하니 한번 검색해서 자세히 알아보시는 것도 좋을 듯 합니다.
참고하시라고 forces, forceCoeffs의 실례를 위해 airfoil force 계산에 관한 controlDict를 아래에 첨부하였습니다. 이는 system/controlDict의 functions에서 수정하여 실행하실 수 있습니다. 주석을 통해 각각이 무슨 기능인지 직접 학습하시기 바라며 이를 이용하면 쉽게 해결하실 수 있을거라 생각합니다.

FoamFile
{
    version         2.0;
    format          ascii;
    class           dictionary;
    object          controlDict;
}

application       icoFoam;
startFrom         startTime;
startTime         0.0;
stopAt            endTime;
endTime           30;
deltaT            0.0025;
writeControl      timeStep;
writeInterval     20;
purgeWrite        0;
writeFormat       ascii;
writePrecision    6;
writeCompression  uncompressed;
timeFormat        general;
timePrecision     6;
runTimeModifiable yes;

functions
(
    forces
    {
        type forces;              // use the forces class
        functionObjectLibs ("libforces.so"); // "plugin" library to load
        patches (naca0012); // Name of patches to integrate forces
        outputControl outputTime; // When to write the output
        rhoName rhoInf;           // Name of the reference density
        rhoInf 1.0;               // Reference density for fluid
        CofR (0.33333 0 0);       // Origin for moment calculations
    }
    forceCoeffs
    {
        type forceCoeffs;     // compute force-coefficients
        functionObjectLibs ("libforces.so");
        patches (naca0012);
        outputControl outputTime;
        rhoName rhoInf;
        rhoInf 1.0;
        CofR (0.33333 0 0);
        liftDir (0 1 0);      // lift direction (parallel to U_inf)
        dragDir (1 0 0);      // drag direction (normal to U_inf)
        pitchAxis (0 0 1);    // rotational moment axis
        magUInf 1;            // magnitude of U_inf (for normalization)
        lRef 1.0;             // reference length
        Aref 1.0;             // reference area
    }
);

좋은 결과 있길 바라겠습니다.