[TimotheusFlavii]

Hi!

I have to simulate flow through a packed bed in which a reaction is catalyzed by catalysts on the spheres of the packed bed. The person who assigned this case wants that the inlet and the outlet are made periodic for the velocity profile, pressure and temperature, and he has defined a mass flow rate for me.
I know periodic boundary conditions and species transport are not compatible, and I therefore think a UDF is needed. I don't know however how to write this UDF while still defining a mass flow rate. Could someone help me?

[TimotheusFlavii]

I currently have written this UDF, but I don't know how I should specify the mass flow rate:

#include "udf.h"

DEFINE_PROFILE(x_velocity_inlet, thread, index)
{
real x[ND_ND];
face_t f;
real Uvel[ND_ND];
Domain *domain;
Thread *tf;

int Zone_ID=16;
domain=Get_Domain(1);
tf= Lookup_Thread(domain, Zone_ID);

begin_f_loop(f, tf)
{
Uvel[1]=F_U(f, tf);
}
end_f_loop(f, tf)

begin_f_loop(f, thread)
{
F_PROFILE(f, thread, index) = (Uvel[1]);
}
end_f_loop(f,thread)
}

DEFINE_PROFILE(y_velocity_inlet, thread, index)
{
real x[ND_ND];
face_t f;
real Vvel[ND_ND];
Domain *domain;
Thread *tf;

int Zone_ID=16;
domain=Get_Domain(1);
tf= Lookup_Thread(domain, Zone_ID);

begin_f_loop(f, tf)
{
Vvel[1]=F_V(f, tf);
}
end_f_loop(f, tf)

begin_f_loop(f, thread)
{
F_PROFILE(f, thread, index) = (Vvel[1]);
}
end_f_loop(f,thread)
}

DEFINE_PROFILE(z_velocity_inlet, thread, index)
{
real x[ND_ND];
face_t f;
real Wvel[ND_ND];
Domain *domain;
Thread *tf;

int Zone_ID=16;
domain=Get_Domain(1);
tf= Lookup_Thread(domain, Zone_ID);

begin_f_loop(f, tf)
{
Wvel[1]=F_W(f, tf);
}
end_f_loop(f, tf)

begin_f_loop(f, thread)
{
F_PROFILE(f, thread, index) = (Wvel[1]);
}
end_f_loop(f,thread)
}

DEFINE_PROFILE(pressure_inlet, thread, index)
{
real x[ND_ND];
face_t f;
real Pres[ND_ND];
Domain *domain;
Thread *tf;

int Zone_ID=16;
domain=Get_Domain(1);
tf= Lookup_Thread(domain, Zone_ID);

begin_f_loop(f, tf)
{
Pres[1]=F_P(f, tf);
}
end_f_loop(f, tf)

begin_f_loop(f, thread)
{
F_PROFILE(f, thread, index) = (Pres[1]);
}
end_f_loop(f,thread)
}

DEFINE_PROFILE(temperature_inlet, thread, index)
{
real x[ND_ND];
face_t f;
real Temp[ND_ND];
Domain *domain;
Thread *tf;

int Zone_ID=16;
domain=Get_Domain(1);
tf= Lookup_Thread(domain, Zone_ID);

begin_f_loop(f, tf)
{
Temp[1]=F_T(f, tf);
}
end_f_loop(f, tf)

begin_f_loop(f, thread)
{
F_PROFILE(f, thread, index) = (Temp[1]);
}
end_f_loop(f,thread)
}

[moun139]

I think there is a problem with V velocity !!