[saratunni07]

Kindly provide me suggestions related to a UDF that has been wrote in Visual Studio Code. I need suggestions on the code whether its right or need restructuring. It has to be hooked up in DPM user scalars for obtaining weber number so that i can understand the wall film interaction.

#include "udf.h"
#define DPM_SURFTEN 0.065

DEFINE_INIT(weber_setup,domain)
{
/* if memory for the particle variable titles has not been
* allocated yet, do it now */
if (NULLP(user_particle_vars)) Init_User_Particle_Vars();
/* now set the name and label */
strcpy(user_particle_vars[0].name,"weber-number");
strcpy(user_particle_vars[0].label,"Weber Number");
strcpy(user_particle_vars[1].name,"weber-number-0");
strcpy(user_particle_vars[1].label,"Weber Number 0");
}

/* update the user scalar variables */
DEFINE_DPM_SCALAR_UPDATE(weber_number,cell,thread, initialize,p)
{
if (initialize)
{
/* this is the initialization call, set:
* P_USER_REAL(p,0) contains the weber number, initialize to 0 */
P_USER_REAL(p,0) = 0.;
}
else
{
real p_vel=sqrt(P_VEL(p)[0]^2+P_VEL(p)[1]^2+P_VEL(p)[2]^2);
P_USER_REAL(p,0)=P_RHO(p)*(p_vel^2)*P_DIAM(p)/DPM_SURFTEN(p);
}
}

Kindly help.

Thanks Sarat