3D free surface and VOF

Orome

Hi, we are trying to simulate a 3D channel with a free surface. Our problem rise when we want to include the udf function outlet_pressure:

#include "udf.h"

#include "sg_mphase.h"

#define RHO_WATER 998.2

#define RHO_AIR 1.225

#define G 9.81

DEFINE_PROFILE(outlet_pressure,thread,nv)

{ face_t f;

real ymax=-1e4;

real x[ND_ND];

Thread *t_phase2=THREAD_SUB_THREAD(thread,1);

begin_f_loop(f,thread) {

F_CENTROID(x,f,thread);

if ((ymax < x[1]) && (F_VOF(f,t_phase2) > 0.5)) ymax=x[1]; }

end_f_loop(f,thread)

begin_f_loop(f,thread) {

F_CENTROID(x,f,thread);

F_PROFILE(f,thread,nv)=0;

if (x[1] < ymax) F_PROFILE(f,thread,nv)=G*(RHO_WATER-RHO_AIR)*(ymax-x[1]); }

end_f_loop(f,thread) }

Then Fluent gives the message

Error: fluent.6.0.20 received a fatal signal (SEGMENTATION VIOLATION). Error Object: #f

and it exit. We have no idea where is the error.

Thank you in advance for your help.