Tip Surface elevation in time

1 Introduction

Working with free surface flow, one might want to track a point at the interface in time, possibly as some kind of buoy.

A common method to achieve this (e.g. in interFoam), is to sample the value of alpha1 along a line (x,y,z(t)) and averaging z(t) based on the sampled value alpha1(z). Approaches following this procedure can be found at:

• http://www.cfd-online.com/Forums/openfoam-post-processing/64462-how-monitor-free-surface-elevation-vs-time.html
• http://www.cfd-online.com/Forums/openfoam/64888-help-modified-probe-file-wave-elevation-vs-time.html

This wiki-page will describe a different method I used, based on first sampling the iso-surface of alpha1=0.5 (which I need anyhow) and subsequentlly sampling points from this iso-surface using Python.

2 Iso-surface

Sampling the iso-surface of alpha1 is straight-forward to perform during runtime. It can be obtained by putting the following sampling function at the end of the controlDict of case.

functions
(
   freeSurface
   {   
       type            surfaces;
       functionObjectLibs
       (   
           "libsampling.so" 
       );  
       outputControl   outputTime;
       outputInterval  1;  
       surfaceFormat  vtk;
       fields
       (   
           alpha1
       );  
       surfaces
       (   
           freeSurface
           {   
               type        isoSurfaceCell;
               isoField    alpha1;
               isoValue    0.5;
               interpolate false;
               regularise  false;
           }   
       );  
       interpolationScheme cell;
   }  
);