[vgchem2017]

Hello all,

I am working on liquid-solid counter current system with liquid entering from bottom and solid entering at the top. Solid enters only after 10s of flow time. I wrote a udf for solid inlet velocity as follows:

#include "udf.h"
DEFINE_PROFILE(solid_y_velocity,mixture_thread,i)
{
int phase_domain_index;
real flow_time = CURRENT_TIME;
face_t f;
Thread *thread_s;
phase_domain_index = 1;
thread_s = THREAD_SUB_THREAD(mixture_thread, phase_domain_index) ;
{
if (flow_time <= 10.)
{
begin_f_loop(f,mixture_thread)
{
F_PROFILE(f,mixture_thread,i) = 0.;
}
end_f_loop(f,mixture_thread)
}
else
{
begin_f_loop(f,mixture_thread)
{
F_PROFILE(f,mixture_thread,i) = 0.001245;
}
end_f_loop(f,mixture_thread)
}
}
}

I don't face any issue in compiling the udf but after 10 seconds of simulation still velocity of solids is zero.

Kindly help me out with this.

Thank You

[CeesH]

Have you tried supplying the solids thread instead of mixture?

[vgchem2017]

Thank for your reply

Is this the way you are asking me to do?

#include "udf.h"
DEFINE_PROFILE(solid_y_velocity,mixture_thread,i)
{
int phase_domain_index;
real flow_time = CURRENT_TIME;
face_t f;
Thread *thread_s;
phase_domain_index = 1;
thread_s = THREAD_SUB_THREAD(mixture_thread,phase_domain_inde x);
{
if (flow_time <= 10.)
{
begin_f_loop(f,thread_s)
{
F_PROFILE(f,thread_s,i)= 0.;
}
end_f_loop(f,thread_s)
}
else
{
begin_f_loop(f,thread_s)
{
F_PROFILE(f,thread_s,i) = 0.001245;
}
end_f_loop(f,thread_s)
}
}