CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Post-Processing

SampleSurface problem

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 12, 2008, 13:26
Default I am trying to run the sampleS
  #1
New Member
 
Scott Ford
Join Date: Mar 2009
Location: Pennsylvania, USA
Posts: 7
Rep Power: 8
saf5029 is on a distinguished road
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?
saf5029 is offline   Reply With Quote

Old   June 12, 2008, 17:02
Default Hi Scott I think I have see
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,619
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
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
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   June 13, 2008, 14:42
Default Thanks for the input Niels.
  #3
New Member
 
Scott Ford
Join Date: Mar 2009
Location: Pennsylvania, USA
Posts: 7
Rep Power: 8
saf5029 is on a distinguished road
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
saf5029 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SampleSurface with polyhedral mesh braennstroem OpenFOAM Post-Processing 1 February 29, 2008 05:30
Gmsh and samplesurface touf Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 December 10, 2007 03:27
SurfaceWriters other than raw are not selectable in sampleSurface hannes OpenFOAM Bugs 4 August 28, 2007 09:04
Question about sampleSurface osimonsimon OpenFOAM Paraview & paraFoam 0 October 16, 2006 22:32
Zero values using sampleSurface steja OpenFOAM Paraview & paraFoam 2 June 15, 2006 04:23


All times are GMT -4. The time now is 17:51.