[bloodflow]

Hi All, I am trying to deactivate particles within a certain region using the following UDF:

/******* UDF to deactivate particles in specified region *********/

#include "udf.h"
#include "dpm.h"

DEFINE_DPM_SCALAR_UPDATE(deactivation_position,cel l,thread,initialize,p)
{
real x[ND_ND];
real y;
real z;
real position;
x=P_POS(p)[i];
y=x[1];
z=x[2];
position=(y*y+z*z)
if (position>0.00000169)
{
MARK_PARTICLE(p, P_FL_REMOVED);
}

}


However when interpreting the UDF I get the error: "line 12: P_POS: undeclared variable"

I'm not sure how this simple UDF has gone wrong....

Any help is much appreciated

[pakk]

I don't know if this solves your problem, but I corrected some mistakes below.

/******* UDF to deactivate particles in specified region *********/
 
#include "udf.h"
#include "dpm.h"

DEFINE_DPM_SCALAR_UPDATE(deactivation_position,cell,thread,initialize,p)
{
real x[ND_ND];
real y;
real z;
real position;
x=P_POS(p);
y=x[0];
z=x[1];
position=(y*y+z*z);
        if (position>0.00000169)
        {
        MARK_PARTICLE(p, P_FL_REMOVED);
        }

 }

[bloodflow]

Originally Posted by pakk

I don't know if this solves your problem, but I corrected some mistakes below.

/******* UDF to deactivate particles in specified region *********/
 
#include "udf.h"
#include "dpm.h"

DEFINE_DPM_SCALAR_UPDATE(deactivation_position,cell,thread,initialize,p)
{
real x[ND_ND];
real y;
real z;
real position;
x=P_POS(p);
y=x[0];
z=x[1];
position=(y*y+z*z);
        if (position>0.00000169)
        {
        MARK_PARTICLE(p, P_FL_REMOVED);
        }

 }

Still Has the same problem.

[AlexanderZ]

compile UDF

best regards

[pakk]

/******* UDF to deactivate particles in specified region *********/

#include "udf.h"
#include "dpm.h"

DEFINE_DPM_SCALAR_UPDATE(deactivation_position,cell,thread,initialize,p)
{
real y;
real z;
real position;
y=TP_POS(p)[1];
z=TP_POS(p)[2];
position=(y*y+z*z);
if (position>0.00000169)
{
MARK_PARTICLE(p, P_FL_REMOVED);
}

}

[bloodflow]

Originally Posted by pakk

/******* UDF to deactivate particles in specified region *********/

#include "udf.h"
#include "dpm.h"

DEFINE_DPM_SCALAR_UPDATE(deactivation_position,cell,thread,initialize,p)
{
real y;
real z;
real position;
y=TP_POS(p)[1];
z=TP_POS(p)[2];
position=(y*y+z*z);
if (position>0.00000169)
{
MARK_PARTICLE(p, P_FL_REMOVED);
}

}

That is fantastic!

Now I am going to try and include a velocity component to this, and look at best way to analyse the results. If I am exporting particle tracking data as a .xml from fluent, is there a way to view these ones that I am removing? Or do I have to look at this with a UDF?