postProcess: crash with surface sampling
Summary
Running the postProcess utility to extract a vtk surface (e.g. of a phase fraction alpha) results in a crash after the first time step has been processed. I noticed this bug in OpenFOAM v2006, not sure about earlier versions. However, this problem does not appear in OpenFOAM 7 and OpenFOAM 8.
Steps to reproduce
Run a simple interFoam multiphase tutorial case (for example damBreakWithObstacle) and use the following surface file to extract the free surface in VTK format with:
postProcess -func surface
The utility postProcess will create a VTK surface for the first time step but will crash at the second. When using the tool only for the second time step via
postProcess -func surface -time 0.02
it works. Letting it run for several time steps starting with the second one like
postProcess -func surface -time 0.02:
works for the second time step, but it eventually crashes on the third time step.
What is the current bug behaviour?
Sample surfaces in VTK for all time steps.
What is the expected correct behavior?
Only the first time step gets sampled in a series of time steps.
Relevant logs and/or images
Create time
Create mesh for time = 0
Sampled surface:
cuttingPlane -> vtk
Time = 0
Reading fields:
Executing functionObjects
Time = 0.02
Reading fields:
Executing functionObjects
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3 void Foam::volPointInterpolation::interpolateInternalField<double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>&) const at ??:?
#4 void Foam::volPointInterpolation::interpolate<double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>&) const at ??:?
#5 Foam::sampledIsoSurface::getIsoFields() const at ??:?
#6 Foam::sampledIsoSurface::updateGeometry() const at ??:?
#7 Foam::sampledSurfaces::performAction(unsigned int) at ??:?
#8 Foam::functionObjects::timeControl::write() at ??:?
#9 Foam::functionObjectList::execute() at ??:?
#10 ? at ??:?
#11 ? at ??:?
#12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#13 ? at ??:?
Speicherzugriffsfehler
Environment information
- OpenFOAM version : OpenFOAM v2006
- Operating system : Debian 10
- Hardware info :
- Compiler : GCC 9.1.0