[sap_92]

I am running a simple CFD simulation of flow in a backward-facing channel. I would like to monitor the average value of the square of velocity over the whole domain (sum of the square of velocity at all grid points/ number of grid points) as the solution progresses. I would like to write this quantity in a log file along with residual information. If it is not easy to write along with residual information, is there any other way that I can write the user-defined quantity at every time step in a separate log file. Basically, I would like to use the history of this quantity to alter some of the simulation parameters while the simulation is running.

Thank you for your help.

[crubio.abujas]

Hello Suraj!

One possibility is to create your own functionObject to do so. I don't know exactly what residuals are you interested on, but as a concept trial I've implemented a functionObject that returns the average square U of the entire domain along with the residual of U in each direction and write them into postProcessing/0/squareU. Right now the residuals are hardcoded, so if you need another field you must change it in the code.

I let you a link to the source code here
You just have to position yourself in the folder and execute wmake. Once compile just add the following code to controlDict functions dictionary.

    squareU
    {
        type           squareU;
        libs           ("libsquareUFunctionObject.so");
        writeControl   timeStep;
        writeInterval  1;
    }

Other option could be add a codedFixed function object to print the U² field and then process it with pyFoam. But that may require to manage multiple files, and so.

Concerning the changes in the controlDict, I think there shall be possible to code some rules to automatically change them, but probably you want to control it yourself by hand. Nevertheless if you think you get a strategy to work with you can check how to code it.

Let me know if it works!

[RGS]

Originally Posted by sap_92

I am running a simple CFD simulation of flow in a backward-facing channel. I would like to monitor the average value of the square of velocity over the whole domain (sum of the square of velocity at all grid points/ number of grid points) as the solution progresses. I would like to write this quantity in a log file along with residual information. If it is not easy to write along with residual information, is there any other way that I can write the user-defined quantity at every time step in a separate log file. Basically, I would like to use the history of this quantity to alter some of the simulation parameters while the simulation is running.

Thank you for your help.

Hello Suraj,


I think it is possible to get a log file of the residuals using swak4foam - at least the residuals that are reported by the solver and appear in the terminal. I hope this helps.

[sap_92]

Originally Posted by crubio.abujas

Hello Suraj!

One possibility is to create your own functionObject to do so. I don't know exactly what residuals are you interested on, but as a concept trial I've implemented a functionObject that returns the average square U of the entire domain along with the residual of U in each direction and write them into postProcessing/0/squareU. Right now the residuals are hardcoded, so if you need another field you must change it in the code.

I let you a link to the source code here
You just have to position yourself in the folder and execute wmake. Once compile just add the following code to controlDict functions dictionary.

    squareU
    {
        type           squareU;
        libs           ("libsquareUFunctionObject.so");
        writeControl   timeStep;
        writeInterval  1;
    }

Other option could be add a codedFixed function object to print the U² field and then process it with pyFoam. But that may require to manage multiple files, and so.

Concerning the changes in the controlDict, I think there shall be possible to code some rules to automatically change them, but probably you want to control it yourself by hand. Nevertheless if you think you get a strategy to work with you can check how to code it.

Let me know if it works!

When I executed the wmake command in the folder, I got an error that Residuals.H file is not in the directory. I think you forgot to attach the Residual.H file in the source code . Can you please attach it to the source code? Thank you very much for your help.

[crubio.abujas]

Residual.H should be defined within the finiteVolume dependencies (inside Make/options). At least that it the case for OF6. Which version of OF are you using? Maybe other versions have different management of this dependency.

[sap_92]

Originally Posted by crubio.abujas

Residual.H should be defined within the finiteVolume dependencies (inside Make/options). At least that it the case for OF6. Which version of OF are you using? Maybe other versions have different management of this dependency.

I am using OF5. I downloaded the file, extracted it, and then entered into the extracted folder. Then, I just typed wmake in the terminal. Am I a missing something? There is a residuals.H file in the OpenFOAM-5.x/src/functionObjects/utilities in my installation directory.
Thank you.

[crubio.abujas]

It seems that on OF7 they added a new Residual class (note the capitalization) to ease the manage of those, but OF5 is lacking for it. You've done it right, it just that the code was not adapted for your specific version.

I've tested this new version in OF5 and worked for me. The steps to compile it are the same. Let me know if you find further problems!

[sap_92]

Originally Posted by crubio.abujas

It seems that on OF7 they added a new Residual class (note the capitalization) to ease the manage of those, but OF5 is lacking for it. You've done it right, it just that the code was not adapted for your specific version.

I've tested this new version in OF5 and worked for me. The steps to compile it are the same. Let me know if you find further problems!

Thank you very much. It worked. This is exactly what I was looking for. Can you please point me the place where I can make changes so that residuals for all equations will be written in a squareU.dat file. For example, I would like to know the history of all residuals which appears in a log file when you run the simulation (simpleFoam > log).

Thank you very much for your help.

[crubio.abujas]

To change the code just modify the write() function inside squareU.C, You need to know the type of data (U-> vector, p->scalar).

bool Foam::functionObjects::squareU::write()
{
    logFiles::write();

    writeTime(file());
    writeSquareU();        
    writeResidual<vector>("U");


    // Added pressure residuals

    writeResidual<scalar>("p");

    file() << endl;
    return true;
}

In order to remember it is a good idea to change the header as well. Inside the same file change writeFileHeader()

void Foam::functionObjects::squareU::writeFileHeader(const label i)
{
    if (Pstream::master())
    {
        writeHeader(file(), "squareU");
        writeCommented(file(), "Time\tsquareU");
        writeFileHeader<vector>("U");
        writeFileHeader<scalar>("p");
    }

    file() << endl;
}

Once change, recompile and run!