[BasG]

Hello foamers,

I am trying to write the field average of U^2 - as a measure of field-averaged kinetic energy. So far I am able to write the field average of U by using:

PHP Code:

functions
{
    volFieldAverage_VA
    {
        // Mandatory entries (unmodifiable)
        type            volFieldValue;
        libs            ("libfieldFunctionObjects.so");

        // Mandatory entries (runtime modifiable)
        fields          (p U k);
        operation       volAverage;
        regionType      all;

        // Optional entries (runtime modifiable)
        weightField     alpha.water;

        // Optional (inherited) entries
        writeFields     no;
        writeToFile     true;
        writeControl    timeStep;
        writeInterval   1;
    }
} 


Additionally, I want to obtain the volAverage for the field U^2. Could you help me doing so?

Best,

Bas

[SHUBHAM9595]

You can use swak4Foam for that...The following snippet should give you an idea...you might also need to modify it if some syntax error occurs

kineticEnergyExpression
    {
        type expressionField;
        outputControl timeStep;
        outputInterval 1;
        fieldName kineticEnergy;
        variables "Ukin=sum(U*U*vol())/sum(vol());";
        expression "Ukin";
        autowrite true;
    }

expression "pAverage-p";
Finally, make sure to add the corresponding libraries in the controlDict file....something like

libs (
    "libOpenFOAM.so"
    "libsimpleFunctionObjects.so"
    "libsimpleSwakFunctionObjects.so"
    "libswakFunctionObjects.so"
    );

[BasG]

Thank you SHUBHAM9595 for your reply. Unfortunately I am running my models on a linux cluster owned by my institute and the swak4foam library is not yet installed. It will take quite a procedure to get it installed I am afraid. Are you aware of alternative ways to perform this task?

[SHUBHAM9595]

Alternatively, you can define a custom volScalarField in createFields.H (anywhere after the U ) as shown below

volScalarField KE
(
    IOobject
    (
        "KE",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mag(U)*mag(U)
);

AUTO_WRITE will ensure that it write the Kinetic energy field at each writing interval. Additionally, to update its value with each iteration you need to insert the following line at the end of UEqn.H

KE = mag(U)*mag(U);

Good luck !!!

[Yann]

Hi Bas,

If you are using a relatively new version of OpenFOAM, you might be able to do it with the "multiply" function object to create the U^2 field before your volAverage function.
Here is the syntax for OpenFOAM-9:

functions
 {


     Usquare
    {
        type            multiply;
        libs            ("libfieldFunctionObjects.so");
        fields          (U U);
        result          Usquare;
        executeControl  writeTime;
        writeControl    writeTime;
    }

    volFieldAverage_VA
    {
        // Mandatory entries (unmodifiable)
        type            volFieldValue;
        libs            ("libfieldFunctionObjects.so");

        // Mandatory entries (runtime modifiable)
        fields          (p U k Usquare);
        operation       volAverage;
        regionType      all;

        // Optional entries (runtime modifiable)
        weightField     alpha.water;

        // Optional (inherited) entries
        writeFields     no;
        writeToFile     true;
        writeControl    timeStep;
        writeInterval   1;
    }
}

Cheers,
Yann

[BasG]

Thanks SHUBHAM9595! This is useful. I am now able to write the field for U^2 for every map output time. Eventhough it is not for every timestep (as are the volume averages of the other quantities), it is still usefull as I can determine the domain average in Paraview.

@Yann, thanks for your reply. I tried this as well. Unfortunately the version I use (v1812) does not have this functionality yet.