[vkgupta1211]

Hi
I want to compute the total volume of cells where the secondary phase volume fraction is from 0.5 to 1. I have two fluid zones (ID:29 & 30).

here is the UDF. Compilation is successful, but on running, simulations gets crashed. Can somebody help me?

#include "udf.h"
/* UDF to calculate the total volume of cells with volume fraction between 0.5 and 1 */
DEFINE_ADJUST(sum_vol,d)
{
#if !RP_HOST
Thread *t;
cell_t c;
t=Lookup_Thread(d,30);
Thread **pt=THREAD_SUB_THREADS(t);
Thread *tp=THREAD_SUB_THREAD(t,1);
Domain *sd=DOMAIN_SUB_DOMAIN(d,1);
real total_vol=0.0;
mp_thread_loop_c(t,d,pt)
{
if (FLUID_THREAD_P(t))
{
begin_c_loop(c,t)
{
if ((C_VOF(c,tp)>=0.5) && (C_VOF(c,tp)<=1.0))
total_vol+=C_VOLUME(c,t);
C_UDMI(c,t,0)=total_vol;
}
end_c_loop(c,t)
}
}
Message0("\nNumber of particles trapped on the interface = %d\n\n",total_vol);
#endif
}