[saf5029]

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:stream&) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt in "/lib64/tls/libc.so.6"
#3 Foam::triSurfaceTools::surfaceNormal(Foam::triSurf ace const&, int, Foam::Vector<double> const&) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libtriSurface.so"
#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/libmeshTools.so"
#5 Foam::orientedSurface::orient(Foam::triSurface&, Foam::Vector<double> const&, bool) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libmeshTools.so"
#6 Foam::orientedSurface::orientedSurface(Foam::triSu rface const&, Foam::Vector<double> const&, bool) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libmeshTools.so"
#7 Foam::meshCutSurface::meshCutSurface(Foam::cellDec ompCuts const&) in "/home/saf5029/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#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/libc.so.6"
#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 this...my sampleSurfaceDict looks like the following:

surfaceFormat raw;
interpolationScheme cellPointFace;

surfaces
{
interpolatedIsoSurface
{
name freeSurface;
field gamma;
value 0.5;
}
);

fields
(
gamma
);

Can anybody give me any help?

[ngj]

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,

Niels

[saf5029]

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?

Thanks,
Scott