[HagawaArashi]

hello everyone! I encounter a problem about post-processing.
I want to use "volTensorField grad(U)" to obtain the viscous dissipation field.
first, I calculate the grad(U), and then use the tensor calculus in OpenFOAM to obtain the volScalarField viscousDissipation.
But the code only can perform in "0/", neglecting the result file at latest time.
As a result, the "viscousDissipation" field can only be found in "0/", others are empty.
So how can I perform the code also in other time named files?
Sincerely wish some advice!Thank you!

And my code is follows, it can be complied successfully without any error:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include "fvCFD.H"

int main(int argc, char *argv[])
{
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"

volVectorField U
(
IOobject
(
"U",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
),
mesh
);

volScalarField viscousDissipation
(
IOobject
(
"viscousDissipation",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("viscousDissipation", dimensionSet(1,-1,-3,0,0,0,0),0.0)
);

volTensorField Symmetry_ =fvc::grad(U);
forAll(mesh.C(), cellI)
{
Symmetry_[cellI] = symm(Symmetry_[cellI]);
}

forAll(mesh.C(), cellI)
{
viscousDissipation[cellI] = 2 * operator||(Symmetry_[cellI], Symmetry_[cellI]);
}

viscousDissipation.write();
return 0;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

[olesen]

Take a look at various utilities that use the timeSelector and thus have the "-latestTime" (and other time selection methods). You need to place a loop over selected times.



See, for example, applications/utilities/postProcessing/miscellaneous/postChannel/postChannel.C as one of the simplest examples that I can think of.