Skip to content

GitLab

Open
Created by Martin Heinrich@martinHeinrich

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

Possible fixes