[EnigmaFluids]

Hello,

I am using OpenFOAM 2.3.0 probeLocations function to gather property data for plotting over time. Using this produces a 0/ directory and a 0.001/ directory (the first saved time step). Inside the 0/ directory we have property files like temperature T, pressure p, etc. Inside the 0.001/ directory we have a single rho file. This is because rho is not calculated at time 0. However I am building a script to compile all the sample data for plotting in excel and it would be convenient if I could calculate the rho property for time 0 so that all the property files end up in the 0/ directory. Does anyone know of a convenient way to do this. Is there a subfunction in foamCalc that I've missed?

Thank you for your time.

More background: The reason I am using such an old version of OpenFOAM is that my research lab has created a few custom solvers that have not been updated to work with the new library structure of OpenFOAM 4+

[ano]

Hi EnigmaFluids,


I didn't try myself, but perhaps you get it, if you modify your solver in a manner that it outputs rho from the beginning. After sourcing your Openfoam version, you would have to type in your terminal:

sol

which brings you to the source code for the solvers. Go to the directory which is called like your solver and open 'createFields.H'


Below the definition of rho, which will look similar to:

    volScalarField rho
    (
        IOobject
        (
            "rho",
            runTime.timeName(),
            mesh
        ),
        thermo.rho()
    );

add a line

rho.write();

Afterwards close the file and recompile by typing

wmake

Running this modified solver should give you the rho at the initial time step. In case you don't want to rerun the whole case, you can change your controlDict to rerun only from time 0 to 0.

[EnigmaFluids]

Ano,

This was very helpful. Thanks!