[Ariyan]

Dear Friends

I have developed a simple Drag force model code. my computational domain includes 4 different zones. the drag formulation is the same for all zones, except the constants will change. To specify which formulation applies to each zone, I plan to utilize conditional statements with the "IF" condition based on the zone_ID. could you please advise on how and where I should consider "IF" statement?


#include "udf.h"

DEFINE_EXCHANGE_PROPERTY(DRAG,c,mix_thread,s_col,f _col)
{
Thread *gas, *liq;
real r_vel_g,z_vel_g,r_vel_l,z_vel_l,abs_u,slip_r,slip_ z,rho_g,void_g,mu_g,drag_f;

/* find the threads for the gas(primary) */
/* find the threads for the liquid(secondary) */
gas=THREAD_SUB_THREAD(mix_thread,s_col);
liq=THREAD_SUB_THREAD(mix_thread,f_col);

/* find phase velocity and properties */
r_vel_g=C_U(c,gas);
z_vel_g=C_V(c,gas);
r_vel_l=C_U(c,liq);
z_vel_l=C_V(c,liq);
slip_r=r_vel_g-r_vel_l;
slip_z=z_vel_g-z_vel_l;

/* absolute slip */
abs_u=sqrt(pow(slip_r,2)+pow(slip_z,2));

/* properties */
rho_g=C_R(c,gas);
mu_g=C_MU_L(c,gas);
void_g=C_VOF(c,gas);

/* drag calculation */

/* Zone_ID=1 */
drag_f=121564*mu_g/void_g+11*rho_g*abs_u;

/* Zone_ID=2 */
drag_f=141266*mu_g/void_g+21*rho_g*abs_u;

/* Zone_ID=3 */
drag_f=181786*mu_g/void_g+11*rho_g*abs_u;

/* Zone_ID=4 */
drag_f=201454*mu_g/void_g+51*rho_g*abs_u;

return drag_f;
}

[AlexanderZ]

you may check in Fluent GUI ID of each zone
each ID will be mix_thread value
so if your zones have ID (1, 5, 10, 15)
if statement could look

/* Zone_ID=1 */
if (mix_thread == 1)
drag_f=121564*mu_g/void_g+11*rho_g*abs_u;
else if (mix_thread == 5)
/* Zone_ID=5 */
drag_f=141266*mu_g/void_g+21*rho_g*abs_u;
else if (mix_thread == 10)
/* Zone_ID=10 */
drag_f=181786*mu_g/void_g+11*rho_g*abs_u;
else if (mix_thread == 15)
/* Zone_ID=15 */
drag_f=201454*mu_g/void_g+51*rho_g*abs_u;