[causuya]

Hi everyone,

I am using UDF to do a transient Dust Load Simulation for Air Cleaner(as SAE paper 2006-01-1316 describled), both DPM and UDF are new to me, so I hope I can get some help here.

I have particles with 4 different size injected in, and they will all pass through a internal surface boundary condition. I need compute their distribution on this surface.

So, first, I want to calculate the number of particles of each size, and then their mass.

Can anyone tell me how to calculate the number of DPM particles in UDF?

Thanks a million in advance.

[asal]

Hello causuya

Do you find any solution for calculate the number of DPM particles in UDF? I am stuck in a same situation.
thanks.

[gearboy]

Originally Posted by causuya

Hi everyone,

I am using UDF to do a transient Dust Load Simulation for Air Cleaner(as SAE paper 2006-01-1316 describled), both DPM and UDF are new to me, so I hope I can get some help here.

I have particles with 4 different size injected in, and they will all pass through a internal surface boundary condition. I need compute their distribution on this surface.

So, first, I want to calculate the number of particles of each size, and then their mass.

Can anyone tell me how to calculate the number of DPM particles in UDF?

Thanks a million in advance.

Here is an example. Remember to hook these two functions and check "interaction with continuous phase".


#include "udf.h"
#define FACE_ID 4 /*id of your internal face*/
#define TOTAL_INJECTIONS 4
char* Injection_names[TOTAL_INJECTIONS]={"injection-0","injection-1","injection-2","injection-3"}; //injection names in your setup
int total_particles[TOTAL_INJECTIONS]={0}; // particles from each injection
real total_mass[TOTAL_INJECTIONS]={0.0}; // mass from each injection
DEFINE_DPM_SCALAR_UPDATE(calc_size,c,t,initialize, p)
{
int i;
face_t f;
Thread*tf;
Domain*domain;
domain=Get_Domain(1);
tf=Lookup_Thread(domain,FACE_ID);
begin_f_loop(f,tf)
{
if(F_C0(f,tf)==c && THREAD_ID(t)==THREAD_ID(THREAD_T0(tf))) /* if c is neighbor to your internal face*/
{
Message0("current particle's inject name is %s, mass is %g\n",p->injection->name,P_MASS(p));
for(i=0;i<TOTAL_INJECTIONS;i++)
{
if(STREQ(Injection_names[i],p->injection->name))
{
total_particles[i]++;
(total_mass[i])+=P_MASS(p);
}
}
}
}
end_f_loop(f,tf);
}

DEFINE_ADJUST(report,domain)
{
int i;
for(i=0;i<TOTAL_INJECTIONS;i++)
{
Message0("-----------------------------------------------------\n");
Message0("total particles from %s is %d, total mass is %g\n",
Injection_names[i],total_particles[i],total_mass[i]);
Message0("-----------------------------------------------------\n");
total_particles[i]=0;
total_mass[i]=0.0;
}
}

[asal]

Dear Gearboy
Thanks a lot for your reply.
I try to compile the UDF which you offer and I got this error:

"The UDF library you are trying to load (libudf) is not compiled for 3ddp on the curent platform (win64).
The system cannot find the file specified. "

I try to follow the instruction here:

http://www.cfd-online.com/Wiki/Fluen...Fload_a_UDF.3F

but it doesn't work. any idea to fix it?
thanks.

[Erfan75]

Hello, Gearboy
could you please help me?
i need a udf that can calculate the number of trapped particle?
thanks alot

[Mraz]

Originally Posted by Erfan75

Hello, Gearboy
could you please help me?
i need a udf that can calculate the number of trapped particle?
thanks alot

Hi

Could you find out how to calculate the number of trapped particle?

I have a similar problem. That would be great if you help me with this.

Thanks

[DnyanMiri]

Originally Posted by Mraz

Hi

Could you find out how to calculate the number of trapped particle?

I have a similar problem. That would be great if you help me with this.

Thanks

You can see the number of trapped particles in the DPM summary report. You can go to particle tracking -> select injection -> summary -> display

[Mraz]

Originally Posted by DnyanMiri

You can see the number of trapped particles in the DPM summary report. You can go to particle tracking -> select injection -> summary -> display

Hi,

Thank you Dnyanesh for the reply. I am interested to have number of trapped particle over time. The summery just report the total number of trapped particle at the end of calculation. Do you have any suggestion for me?

Best,

[pakk]

Originally Posted by Mraz

Hi,

Thank you Dnyanesh for the reply. I am interested to have number of trapped particle over time. The summery just report the total number of trapped particle at the end of calculation. Do you have any suggestion for me?

Best,

Suggestion:
Write a UDF for a DPM boundary condition, that adds one line to a text file and then traps the particle. The line should contain at least the time of the event, perhaps more of you want more information. Then, use a different program (excel, Matlab, your choice) to do the post-processing.