[luckysandhu1431]

Hello

I have a 3D periodic geometry of a concentric annulus with liquid-solid two phase flow. I want to introduce the Krieger & Dougherty Model as a UDF for granular viscosity.

I wrote the following code but it gives me an "received a fatal signal (Segmentation Fault)" error when I try to initialise the solution (Hybrid initialisation).

#include "udf.h"
#include "sg_mphase.h"

DEFINE_PROPERTY(gran_vis_kd,c,t)
{
Thread *part = THREAD_SUB_THREAD(t,3);
real mu_mix;
real B = 2.5;
real phi_m = 0.63;
real mu_liq = 10.0;
real phi = C_VOF(c,part);
mu_mix = mu_liq*pow(1-(phi/phi_m),-B*phi_m);
return mu_mix;
}

I have never used C or Fluent UDFs before so I am not sure what I am doing wrong.

I am really thankful to you all for the assistance.