Learning Forum

[YasserSelima]

try this ..
#include "udf.h"




#define Q 7.0/* particle electric charge*/

#define E 3.0/* electric field magnitude */

#define TSTART 0/* field applied at t = tstart*/
DEFINE_DPM_BODY_FORCE(electric_force, p, i)
{
real bforce = 0;
if (P_TIME(p) >= TSTART)
{
if (i == 0) {
bforce = Q * E*P_VEL(p)[0];
} else if (i == 1) {
bforce = -Q * E*P_VEL(p)[1];
} else {
bforce = 0.0;
}
}
/* an acceleration should be returned */
return (bforce / P_MASS(p));
}

Let me know if you are still getting the same error