[Eman.]

Hi,


I am doing a large LES simulation with OF 2.2.1, therefore the simulation exceeds the wall time and I need to continue it multiple times. However, the problem is that every time that I restart it, it doesn't read the old values of my field averages so it starts averaging those fields all over again. Any ideas how this problem can be solved?

Thanks in advance.

[piu58]

In system.controldict use

startFrom       latestTime;

[Eman.]

Originally Posted by piu58

In system.controldict use

startFrom       latestTime;

thanks for the reply. Of course, I already did that. The simulation starts from the latest time but averaging is done all over again.

[Bazinga]

I did this in v. 4.x without any problems. Can you post the file?

[Eman.]

Originally Posted by Bazinga

I did this in v. 4.x without any problems. Can you post the file?

Thanks a lot for the reply. I'm using 2.3 version and here is the ControlDict content:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
libs (
"libOpenFOAM.so"
"libsimpleSwakFunctionObjects.so"
"libswakFunctionObjects.so"
"libgroovyBC.so"
);


application pimpleFoam;

startFrom latestTime;

//startTime 0;

stopAt endTime;

endTime 10;

deltaT 2e-06;

writeControl timeStep;

writeInterval 500;

purgeWrite 50;

writeFormat ascii;

writePrecision 6;

writeCompression on;

timeFormat general;

timePrecision 6;

runTimeModifiable true;



functions
{
probes
{
type probes;
functionObjectLibs ("libsampling.so");
enabled true;
outputControl timeStep;
outputInterval 1;

fields
(
p
);

probeLocations
(
( -0.016522372 -0.015013 0.020)
( -0.016522372 0.002985 0.020)
( -0.016522372 -0.015013 0.021)
( -0.016522372 0.002985 0.021)
( -0.016522372 -0.015013 0.022)
( -0.016522372 0.002985 0.022)



( -0.016522372 -0.014 0.020)
( -0.016522372 0.0018 0.020)
( -0.016522372 -0.014 0.021)
( -0.016522372 0.0018 0.021)
( -0.016522372 -0.014 0.022)
( -0.016522372 0.0018 0.022)





);

}





fieldAverage1
{
type fieldAverage;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;

fields
(
U
{
mean on;
prime2Mean on;
base time;
}

p
{
mean on;
prime2Mean on;
base time;
}
);
}




surfaceSampling
{

type surfaces;

functionObjectLibs ("libsampling.so");
enabled true;
outputControl outputTime;

interpolationScheme cellPoint;

surfaceFormat vtk;

fields
(
U
p
UMean
pMean
);

surfaces
(
MainMidplane
{
//type plane;
//patches ( lowerWall );
//distance 1E-6;
//interpolate true;
//triangulate false;

type plane; // always triangulated
planeType pointAndNormal;
basePoint (-0.016522372 -0.01501396 0.0);
normalVector (1 0 0);
interpolate true;
}



SecondMidplane
{
//type plane;
//patches ( lowerWall );
//distance 1E-6;
//interpolate true;
//triangulate false;

type plane; // always triangulated
planeType pointAndNormal;
basePoint (-0.016522372 -0.00601399 0.0);
normalVector (0 1 0);
interpolate true;
}

);
}
}


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

[Bazinga]

There are flags like "resetOnRestart" and "resetOnOutput" that you could use which would lead to the behavior you mentioned but I didn't see them added to your file.

I checked my old files from 2.4.0 where I think it worked as intended and the only difference I found is the output line where I have

writeControl writeTime;
writeInterval 1;

instead of

outputControl timeStep;
outputInterval 1;

I do not exactly know the difference between the two.

Also, I have time start set (timeStart) but I do not think that this makes the difference.

If I have some minutes later, I will try to check, if it continues the averaging after a restart of the simulation.

[Eman.]

Originally Posted by Bazinga

There are flags like "resetOnRestart" and "resetOnOutput" that you could use which would lead to the behavior you mentioned but I didn't see them added to your file.

I checked my old files from 2.4.0 where I think it worked as intended and the only difference I found is the output line where I have

writeControl writeTime;
writeInterval 1;

instead of

outputControl timeStep;
outputInterval 1;

I do not exactly know the difference between the two.

Also, I have time start set (timeStart) but I do not think that this makes the difference.

If I have some minutes later, I will try to check, if it continues the averaging after a restart of the simulation.

Thanks a lot, I really appreciate it. I'll try using write control and see if it works.

[Bazinga]

So I just tested it on 2.4.0 with an old simulation of mine. It is continuing to average and is reading from the e.g. UMean file from the last time directory. I compared the UMean file of both time steps and they are basically identical. There is even an output in the terminal during the solving saying

fieldAverage fieldAverage:
Reading/initialising field UMean

Well, I guess you already checked it, but is there and output of the Mean fields in the time directories?

[Eman.]

Originally Posted by Bazinga

So I just tested it on 2.4.0 with an old simulation of mine. It is continuing to average and is reading from the e.g. UMean file from the last time directory. I compared the UMean file of both time steps and they are basically identical. There is even an output in the terminal during the solving saying

fieldAverage fieldAverage:
Reading/initialising field UMean

Well, I guess you already checked it, but is there and output of the Mean fields in the time directories?

I just checked it out. UMean is in all time directories but for some reasons, it doesn't read it when I continue the simulation instead, it recomputes it all over again.

[fumiya]

There was a bug in the fieldAverage function object:
Not continuous averaging when restarting a calculation
https://bugs.openfoam.org/view.php?id=898

I am afraid that you use the version with this bug.

Best regards,
Fumiya

[Eman.]

Originally Posted by fumiya

There was a bug in the fieldAverage function object:
Not continuous averaging when restarting a calculation
https://bugs.openfoam.org/view.php?id=898

I am afraid that you use the version with this bug.

Best regards,
Fumiya

Thanks a lot for the reply. I guess this is the problem, I need to compile another version on the cluster.