[Monsieur]

Hello I want to set a mass-flow outlet like the value from another boundary in my simulation. The source code worked well with a single phase flow, but if I use multiphase the mass flow is not like the value of the other face (reports say so). What can be the problem?


#include "udf.h"

int zone_ID = 4;
real my_mass_flow = 0.0; /* initial value 0 m/s */

DEFINE_ADJUST(get_my_mass_flow,domain)
{

#if !RP_HOST
Thread *t = Lookup_Thread(domain,zone_ID);
Domain *sub_domain;
face_t f;
real sum = 0.0;
my_mass_flow = 0.0;
int phase_domain;

sub_domain_loop(sub_domain, domain, phase_domain)
{
thread_loop_f(t, subdomain)
{
if (DOMAIN_ID(sub_domain) == 3)
{
/* Compute Mass Flow at Outlet1 */
begin_f_loop(f, t)
{
my_mass_flow += F_FLUX(f, t);
}
end_f_loop(f, t)
}
}
}
#endif

my_mass_flow = PRF_GRSUM1(my_mass_flow);
node_to_host_real_1(my_mass_flow);
}

DEFINE_PROFILE(apply_my_mass_flow,t,i)
{
face_t f;
begin_f_loop(f,t)
{
F_PROFILE(f,t,i) = my_mass_flow;
}
end_f_loop(f,t);
}

DEFINE_ON_DEMAND(print_my_mass_flow)
{
#if !RP_NODE
Message("My Mass Flow = %.5e\n",my_mass_flow);
#endif
}

DEFINE_EXECUTE_AT_END(my_exec_at_end)
{
print_my_mass_flow();
}