CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Post-Processing (
-   -   SampleSurface problem (

saf5029 June 12, 2008 13:26

I am trying to run the sampleS
I am trying to run the sampleSurface post processing utility on a case with two phase flow around a cylinder (using lesInterFoam), to capture the free surface elevation in the form of a data file (like x,y,z values). I wanted to capture the raw data of the gamma isosurface for a gamma value of 0.5 but i get the following error when executing the sampleSurface command...

#0 Foam::error::printStack(Foam: in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#2 __restore_rt in "/lib64/tls/"
#3 Foam::triSurfaceTools::surfaceNormal(Foam::triSurf ace const&, int, Foam::Vector<double> const&) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#4 Foam::orientedSurface::propagateOrientation(Foam:: triSurface const&, Foam::Vector<double> const&, bool, int, Foam::Vector<double> const&, Foam::List<int>&) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#5 Foam::orientedSurface::orient(Foam::triSurface&, Foam::Vector<double> const&, bool) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#6 Foam::orientedSurface::orientedSurface(Foam::triSu rface const&, Foam::Vector<double> const&, bool) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#7 Foam::meshCutSurface::meshCutSurface(Foam::cellDec ompCuts const&) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#8 Foam::interpolatedIsoSurface::correct(bool, Foam::volPointInterpolation const&, Foam::dictionary const&, Foam::fieldsCache<double> const&) in "/home/saf5029/OpenFOAM/saf5029-1.4.1/applications/bin/linux64GccDPOpt/sampleSur face"
#9 main in "/home/saf5029/OpenFOAM/saf5029-1.4.1/applications/bin/linux64GccDPOpt/sampleSur face"
#10 __libc_start_main in "/lib64/tls/"
#11 __gxx_personality_v0 in "/home/saf5029/OpenFOAM/saf5029-1.4.1/applications/bin/linux64GccDPOpt/sampleSur face"
Floating point exception

I'm not really sure how to decipher sampleSurfaceDict looks like the following:

surfaceFormat raw;
interpolationScheme cellPointFace;

name freeSurface;
field gamma;
value 0.5;


Can anybody give me any help?

ngj June 12, 2008 17:02

Hi Scott I think I have see
Hi Scott

I think I have seen a similar message. Problem was that for the initial time step the surface was perfectly aligned with a grid line, thus some problems occurs, but temporarily removing it from the time directories made the execution successful.

Good luck,


saf5029 June 13, 2008 14:42

Thanks for the input Niels.
Thanks for the input Niels.

I actually figured out that because of my static mesh, the sampleSurface utility was having problems. The sampleSurface.C file essentially contains a bug that won't run static meshes. So I commented out some things and got it to run.

However, my newest problem is that after I execute the utility, all my files are empty like this...

# "gamma" POINT_DATA 0
# x y z "gamma"


Anybody have some advice?


All times are GMT -4. The time now is 07:39.