[hconel]

Hello,

I'm running a case with dynamic mesh and topoSet. The solver I'm using is pimpleDyMFoam. The mesh gets refined as I expect during the run, but the problem is, when I try to continue the job from the latestTime, I run into problems. I get the following error in log.pimpleDyMFoam:

[4] 
[4] 
[4] --> FOAM FATAL ERROR: 
[4] Cannot find directory "polyMesh/sets" in times 173.612294099971336 down to 173.612294099971336
[4] 
[4]     From function Foam::word Foam::Time::findInstance(const Foam::fileName&, const Foam::word&, Foam::IOobject::readOption, const Foam::word&) const
[4]     in file db/Time/findInstance.C at line 175.
[4]

This error is raised for all processes which are supposed to include the set of cells defined by topoSet.

To overcome this, I have run topoSet in parallel before running pimpleDyMFoam to create the sets directory in time directories, but this time, I get the following error:

    Reading/initialising field UMean
[62] 
[62] 
[62] --> FOAM FATAL IO ERROR: 
[62] size 541 is not equal to the given value of 538
[62] 
[62] file: /cfd/honel/OpenFOAM/honel-5.x/c.adapt/processor62/173.612294099971336/UMean.boundaryField.procBoundary62to61 from line 91 to line 92.
[62] 
[62]     From function Foam::Field<Type>::Field(const Foam::word&, const Foam::dictionary&, Foam::label) [with Type = Foam::Vector<double>; Foam::label = long int]
[62]     in file /cfd/honel/OpenFOAM/OpenFOAM-5.x/src/OpenFOAM/lnInclude/Field.C at line 302.
[62] 
FOAM parallel run exiting
[62]

Now, I can not reconstruct the mesh and continue from there, because dynamicMesh tool uses the base mesh to determine refinement levels, and my current mesh is already refined. There should be a way to overcome this problem, but I could not find.

My constant/dynamicMeshDict is as follows:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh   dynamicRefineFvMesh;

// How often to refine
refineInterval  4;

// Field to be refinement on
field           nut;//refineHere;

// Refine field inbetween lower..upper
lowerRefineLevel 0.200;
upperRefineLevel 0.999;

// If value < unrefineLevel unrefine
unrefineLevel   0.05;

// Have slower than 2:1 refinement
nBufferLayers   1;

// Refine cells only up to maxRefinement levels
maxRefinement   2;

// Stop refinement if maxCells reached
maxCells        2000000;

// Flux field and corresponding velocity field. Fluxes on changed
// faces get recalculated by interpolating the velocity. Use 'none'
// on surfaceScalarFields that do not need to be reinterpolated.
correctFluxes
(
    //(phi none)
    //(nHatf none)
    //(rhoPhi none)
    //(alphaPhi none)
    //(ghf none)
);

// Write the refinement level as a volScalarField
dumpLevel       true;