[guohf]

Hello everyone,

I write UDF to loop all injection particles, and the particles will be trap by the wall when the distance between the particle and the wall equal or less than the radius of the particle.
The udf can be 'compiled' and compute. However, it cannot trap particles. Please help me and thanks.

This is my UDF:


#include "stdio.h"
#include "udf.h"
#include "dpm.h"
#include "mem.h"
#include "sg.h"
#include "math.h"
#include "surf.h" /* for macros: RP_Cell() & RP_Thread() */
#include "sg_mphase.h"



#define wall_id 9 /* wall face ID*/
#define dim 2 /* dimension number 2D or 3D */

DEFINE_ADJUST(trapparticle, domain)

{
cell_t pc;
face_t f;
Thread *ft, *tcell;
Injection *Ilist;
Injection *I;
Particle *p;


real fp_dis, x, y, radius,particlenum=0, total_volume=0.0;
real xc[ND_ND], xf[ND_ND];

int i; /*working counter for script*/

ft = Lookup_Thread(domain, wall_id); /* wall face*/

Ilist = Get_dpm_injections();
loop(I, Ilist)
{
loop(p, I->p) /* Standard ANSYS FLUENT Looping Macro to get particle streams in an Injection */
{
pc = P_CELL(p); /* Get the cell and thread that the particle is currently in */
tcell = P_CELL_THREAD(p); /* Particle Structure using new macros*/

radius=P_DIAM(p); /* Particle radius*/
x = P_POS(p)[0];
y = P_POS(p)[1];
C_UDMI(pc,tcell,0) = 0;

begin_f_loop(f,ft)
{
F_CENTROID(xf,f,ft); /*face centerial */
fp_dis=fabs(sqrt(SQR(x-xf[0]) + SQR(y-xf[1]))); /* the distance between particle centre and wall face*/

/* if the distance less than particle radius,particle velocity equal 0 and save particle cell 1 and particle position */
if(fp_dis <= radius)
{
for (i=0; i<dim; i++)
{
P_VEL(p)[i] = 0.0;
}
C_UDMI(pc,tcell,0) = 1;
C_UDMI(pc,tcell,1) = P_POS(p)[0];
C_UDMI(pc,tcell,2) = P_POS(p)[1];
particlenum +=1; /* if particle is trapped, increase particle-number */
MARK_PARTICLE(p, P_FL_REMOVED); /* Remove particle after accretion */

}
else
{
C_UDMI(pc,tcell,0) = 0;
C_UDMI(pc,tcell,1) = 0.0;
C_UDMI(pc,tcell,2) = 0.0;
particlenum +=0;

}
}
end_f_loop(f,ft)

C_UDMI(pc,tcell,3) = particlenum;


}
}
}

[touyet]

Did you get Answer of this problem? please share your UDF of this Method (Collection on Outer Boundary of Particle)

[ctw987]

hi touyet! You seem to like the filter too much.

[ctw987]

hi touyet!! Do you want anonymous?

[bloodflow]

Originally Posted by guohf

Hello everyone,

I write UDF to loop all injection particles, and the particles will be trap by the wall when the distance between the particle and the wall equal or less than the radius of the particle.
The udf can be 'compiled' and compute. However, it cannot trap particles. Please help me and thanks.

This is my UDF:


#include "stdio.h"
#include "udf.h"
#include "dpm.h"
#include "mem.h"
#include "sg.h"
#include "math.h"
#include "surf.h" /* for macros: RP_Cell() & RP_Thread() */
#include "sg_mphase.h"



#define wall_id 9 /* wall face ID*/
#define dim 2 /* dimension number 2D or 3D */

DEFINE_ADJUST(trapparticle, domain)

{
cell_t pc;
face_t f;
Thread *ft, *tcell;
Injection *Ilist;
Injection *I;
Particle *p;


real fp_dis, x, y, radius,particlenum=0, total_volume=0.0;
real xc[ND_ND], xf[ND_ND];

int i; /*working counter for script*/

ft = Lookup_Thread(domain, wall_id); /* wall face*/

Ilist = Get_dpm_injections();
loop(I, Ilist)
{
loop(p, I->p) /* Standard ANSYS FLUENT Looping Macro to get particle streams in an Injection */
{
pc = P_CELL(p); /* Get the cell and thread that the particle is currently in */
tcell = P_CELL_THREAD(p); /* Particle Structure using new macros*/

radius=P_DIAM(p); /* Particle radius*/
x = P_POS(p)[0];
y = P_POS(p)[1];
C_UDMI(pc,tcell,0) = 0;

begin_f_loop(f,ft)
{
F_CENTROID(xf,f,ft); /*face centerial */
fp_dis=fabs(sqrt(SQR(x-xf[0]) + SQR(y-xf[1]))); /* the distance between particle centre and wall face*/

/* if the distance less than particle radius,particle velocity equal 0 and save particle cell 1 and particle position */
if(fp_dis <= radius)
{
for (i=0; i<dim; i++)
{
P_VEL(p)[i] = 0.0;
}
C_UDMI(pc,tcell,0) = 1;
C_UDMI(pc,tcell,1) = P_POS(p)[0];
C_UDMI(pc,tcell,2) = P_POS(p)[1];
particlenum +=1; /* if particle is trapped, increase particle-number */
MARK_PARTICLE(p, P_FL_REMOVED); /* Remove particle after accretion */

}
else
{
C_UDMI(pc,tcell,0) = 0;
C_UDMI(pc,tcell,1) = 0.0;
C_UDMI(pc,tcell,2) = 0.0;
particlenum +=0;

}
}
end_f_loop(f,ft)

C_UDMI(pc,tcell,3) = particlenum;


}
}
}

Did anyone have any progress to this? I am attempting a similar thing