DPM UDF - Segmentation Violation & Stack Backtrace

Prashanth

Hi All

I am working on my Masters Thesis which involves FLUENT modeling of particle deposition. I am a greenhorn when it comes to UDF's and the following is the first UDF I have ever created. The UDF is compiling without any errors and I am also able to choose it as a Boundary Condition in my setup. When I run the particle tracking though, I get the following error:

Stack backtrace generated for node id 3 on signal 11 : Please include this information with any bug report you file on this issue

/nfs/14/osu5206/.bashrc: line 1: module: command not found Primitive Error at Node 3: received Segmentation Violation

I have the code below if that would help in debugging the problem. The BC is basically calculating the particle normal velocity and critical velocity, comparing the both and determining if particle sticks based on this condition. Any help would be greatly appreciated as I have a 5 day deadline to get this working.

#include "udf.h" #include "dpm.h" #define nu_s 0.27 #define nu_p 0.27 #define YMIN 0.0 #define YMAX 0.4 #define UMEAN 1.0 #define B 1./7. #define DELOVRH 0.5 #define VISC 1.7894e-05 #define RGAS (UNIVERSAL_GAS_CONSTANT/MW) #define Tdatum 288.15 #define NUM_UDM 2

DEFINE_DPM_BC(deposition_bc, p, thread, f, f_normal, dim) { Domain *d; cell_t c; Thread *t; d=Get_Domain(1);

real Tavg,Ep,vcr; /*real norm_coeff=1.; real tang_coeff=1.;*/ real vn=0.; real R=287.; real normal[3]; int i,idim=dim; real x[ND_ND];

for (i=0.; i<idim; i++)

normal[i] = f_normal[i];

if(p->type==DPM_TYPE_INERT)

{

/*computing normal velocity*/

for(i=0.;i<idim;i++)

vn += p->state.V[i]*normal[i];

/*computing critical velocity*/

Tavg = (F_T(f,t)+P_T(p))/2.;

Ep = (3.*(pow(10.,20.)))*exp(-0.02365*Tavg);

vcr = pow(((2.*Ep)/P_DIAM(p)),(10./7.));

/*particle continues its path - no sticking*/

if (vn > vcr)

{

return PATH_ACTIVE;

}

/*particle sticks*/

else

{

thread_loop_c(t,d)

{

begin_c_loop(c,t)

{

C_UDMI(c,t,0) += P_MASS(p); /*mass of particles deposited*/

}

end_c_loop(c,t)

}

return PATH_ABORT;

} } }