[floquation]

While regenerating a sampledSurface during the first dynamic mesh refinement (dynamicRefineFvMesh), I receive the following fatal error.
The solver is interDyMFoam, using a 3D uniform blockMesh. The OpenFoam version is 4.0.
I should mention that the functionObject below does work for version 2.4.0.

Fatal Error:

PIMPLE: iteration 1
Selected 400 cells for refinement out of 65536.
Refined from 65536 to 68336 cells.
Forcing regeneration of surface awInterface


--> FOAM FATAL ERROR:
Incorrect size.
mesh: nCells:68336 points:74235
cVals:65536 <-- pre-refinement #cells
cCoords:68336 <-- post-refinement #cells
snappedCc:68336
pVals:74235
pCoords:74235
snappedPoint:74235


    From function void Foam::isoSurface::generateTriPoints(const volScalarField&, const scalarField&, const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&, const Foam::Field<Type>&, const Foam::DynamicList<T>&, const labelList&, const labelList&, Foam::DynamicList<T>&, Foam::DynamicList<int>&) const [with Type = Foam::Vector<double>; Foam::volScalarField = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>; Foam::scalarField = Foam::Field<double>; Foam::labelList = Foam::List<int>]
    in file sampledSurface/isoSurface/isoSurfaceTemplates.C at line 546.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  void Foam::isoSurface::generateTriPoints<Foam::Vector<double> >(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::Field<double> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::Field<Foam::Vector<double> > const&, Foam::DynamicList<Foam::Vector<double>, 0u, 2u, 1u> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::DynamicList<Foam::Vector<double>, 0u, 2u, 1u>&, Foam::DynamicList<int, 0u, 2u, 1u>&) const at ??:?
#3  Foam::isoSurface::isoSurface(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::Field<double> const&, double, bool, double) at ??:?
#4  Foam::sampledIsoSurface::updateGeometry() const at ??:?
#5  Foam::SurfacesRepository::updateSurface(Foam::word const&, Foam::fvMesh const&) at ??:?
#6  Foam::createSampledSurface::redoSurface() at ??:?
#7  Foam::functionObjectList::updateMesh(Foam::mapPolyMesh const&) at ??:?
#8  Foam::polyMesh::updateMesh(Foam::mapPolyMesh const&) at ??:?
#9  Foam::fvMesh::updateMesh(Foam::mapPolyMesh const&) at ??:?
#10  Foam::dynamicRefineFvMesh::refine(Foam::List<int> const&) at ??:?
#11  Foam::dynamicRefineFvMesh::update() at ??:?
#12  ? at ??:?
#13  __libc_start_main in "/lib64/libc.so.6"
#14  ? at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
Aborted

swak4Foam functionObject:
A swak4Foam functionObject was used for generating the sampledSurface, but I don't think that is relevant for the error?
I presume that the calculation of the sampledSurface is part of OpenFoam, and swak4Foam just calls the appropriate function?

createInterface {
  type createSampledSurface;
  outputInterval 1;
  outputControl timeStep;
  surfaceName awInterface;
  surface {
    type isoSurface;
    isoField alpha.water;
    isoValue 0.5;
    interpolate true;
  }
}

No interpolation:
A possible solution is to use "interpolate false;", but that doesn't work for my case:

--> FOAM FATAL IO ERROR:
Non-interpolated iso surface not supported since triangles span across cells.

Known Bug?
According to the swak4Foam ReadMe, this might be a known bug.
However, the date on that document is 2012 and the functionObject does work in v2.4.0, so I'm not sure if this is the same thing.
Besides, I have a dynamicRefineMesh, which is not really a moving mesh, is it?

8.1 Moving meshes and sampledSurfaces
It seems that with moving meshes sampledSurfaces don’t get updated.
This seems to be a problem with OpenFOAM itself (the regular surfaces-functionObject doesn’t get updated. This is currently investigated

Question:
Has anyone encountered this issue before?
Is there any sort of solution known?
Is there any reason in particular that it does work in v2.4.0, but not in v4.0?

Thanks in advance.

[cyw]

Originally Posted by floquation

While regenerating a sampledSurface during the first dynamic mesh refinement (dynamicRefineFvMesh), I receive the following fatal error.
The solver is interDyMFoam, using a 3D uniform blockMesh. The OpenFoam version is 4.0.
I should mention that the functionObject below does work for version 2.4.0.

Fatal Error:

PIMPLE: iteration 1
Selected 400 cells for refinement out of 65536.
Refined from 65536 to 68336 cells.
Forcing regeneration of surface awInterface


--> FOAM FATAL ERROR:
Incorrect size.
mesh: nCells:68336 points:74235
cVals:65536 <-- pre-refinement #cells
cCoords:68336 <-- post-refinement #cells
snappedCc:68336
pVals:74235
pCoords:74235
snappedPoint:74235


    From function void Foam::isoSurface::generateTriPoints(const volScalarField&, const scalarField&, const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&, const Foam::Field<Type>&, const Foam::DynamicList<T>&, const labelList&, const labelList&, Foam::DynamicList<T>&, Foam::DynamicList<int>&) const [with Type = Foam::Vector<double>; Foam::volScalarField = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>; Foam::scalarField = Foam::Field<double>; Foam::labelList = Foam::List<int>]
    in file sampledSurface/isoSurface/isoSurfaceTemplates.C at line 546.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  void Foam::isoSurface::generateTriPoints<Foam::Vector<double> >(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::Field<double> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::Field<Foam::Vector<double> > const&, Foam::DynamicList<Foam::Vector<double>, 0u, 2u, 1u> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::DynamicList<Foam::Vector<double>, 0u, 2u, 1u>&, Foam::DynamicList<int, 0u, 2u, 1u>&) const at ??:?
#3  Foam::isoSurface::isoSurface(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::Field<double> const&, double, bool, double) at ??:?
#4  Foam::sampledIsoSurface::updateGeometry() const at ??:?
#5  Foam::SurfacesRepository::updateSurface(Foam::word const&, Foam::fvMesh const&) at ??:?
#6  Foam::createSampledSurface::redoSurface() at ??:?
#7  Foam::functionObjectList::updateMesh(Foam::mapPolyMesh const&) at ??:?
#8  Foam::polyMesh::updateMesh(Foam::mapPolyMesh const&) at ??:?
#9  Foam::fvMesh::updateMesh(Foam::mapPolyMesh const&) at ??:?
#10  Foam::dynamicRefineFvMesh::refine(Foam::List<int> const&) at ??:?
#11  Foam::dynamicRefineFvMesh::update() at ??:?
#12  ? at ??:?
#13  __libc_start_main in "/lib64/libc.so.6"
#14  ? at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
Aborted

swak4Foam functionObject:
A swak4Foam functionObject was used for generating the sampledSurface, but I don't think that is relevant for the error?
I presume that the calculation of the sampledSurface is part of OpenFoam, and swak4Foam just calls the appropriate function?

createInterface {
  type createSampledSurface;
  outputInterval 1;
  outputControl timeStep;
  surfaceName awInterface;
  surface {
    type isoSurface;
    isoField alpha.water;
    isoValue 0.5;
    interpolate true;
  }
}

No interpolation:
A possible solution is to use "interpolate false;", but that doesn't work for my case:

--> FOAM FATAL IO ERROR:
Non-interpolated iso surface not supported since triangles span across cells.

Known Bug?
According to the swak4Foam ReadMe, this might be a known bug.
However, the date on that document is 2012 and the functionObject does work in v2.4.0, so I'm not sure if this is the same thing.
Besides, I have a dynamicRefineMesh, which is not really a moving mesh, is it?

Question:
Has anyone encountered this issue before?
Is there any sort of solution known?
Is there any reason in particular that it does work in v2.4.0, but not in v4.0?

Thanks in advance.

Have You solved it? Ia also met this problem.