[matiashess]

Hi all, I'm making use of DPM to model a recirculating batch reactor with liquid-solid heterogeneous flow. I need help with a UDF that works on boundary conditions by injecting the same number of particles that left the domain to maintain constant mass.
My boundary conditions were on output "pressure outlet" and "velocity inlet" at the inlet.
Really I would greatly appreciate if you help me with the UDF.


Greetings my dear!

[matiashess]

pleasee someone who can help me? any contribution will help me to advance my project

[macfly]

Hi,

The following example is not about particles, but it shows how you can do some calculations on a boundary at each time step and apply the results somewhere else.

- Area averaged temperature is calculated at the outlet.
- The calculated average temperature is augmented by 1 K and will be applied as the inlet temperature at next time step.

#include "udf.h"
real T_mean = 300.0;  /* defined outside because will be used in multiple DEFINE macros */

DEFINE_EXECUTE_AT_END(execute_at_end)
{
Domain *domain;
Thread *thread;
face_t face;
real area[ND_ND];
real total_area = 0.0;
real total_area_ave_temp = 0.0;
int ID = 18; /* outlet ID displayed in Fluent boundary conditions panel */
domain = Get_Domain(1);
thread = Lookup_Thread(domain, ID); /*change "12" to the actual id of your b.c. */
	begin_f_loop(face, thread)
		F_AREA(area, face, thread);
		total_area += NV_MAG(area);
		total_area_ave_temp += NV_MAG(area)*F_T(face, thread);
	end_f_loop(face, thread)
T_mean = total_area_ave_temp/total_area;
printf("Area averaged T on boundary %d = %f K\n", ID, T_mean);
}

DEFINE_PROFILE(T_profile, thread, position)
{
face_t face;
real T_prof = T_mean + 1;  
	begin_f_loop(face, thread)
		F_PROFILE(face, thread, position) = T_prof;
	end_f_loop(face, thread)
}

I guess it's not a big deal to modify the code for your case, but I don't have time for that. Hope it helps.

[matiashess]

Thanks for your reply François, i'm working on that udf.