Hi

I am trying to implement of plasma actuator body force to ANSYS Fluent. My work is based on Singh and Roy's body force approximation. Equation is attached as a picture.

(Reference: Singh, Kunwar Pal, and Subrata Roy. "Force approximation for a plasma actuator operating in atmospheric air." Journal of Applied Physics 103.1 (2008): 013305. )

I wrote a UDF for this equation and I used that for a flat plate. Problem is, I do not get any reasonable results. When I use 1000 V, I got velocity magnitude in order of 10^5, in other hand when I use 1 KV (the same 1000 V but with changing the unit), I got velocity in order of 10^-2!!!

I do not understand, what is mu work is wrong?IS it about Voltage unit or something else?

I appreciate that if you let me your advise.

Best Regards.

# include "udf.h"

#define Fx0 2.6

#define Fy0 2.0

#define Phi0 1

#define Phi4 pow(Phi0,4)

#define Betax 800000

#define Betay 10000000

#define X0 0.015

#define Y0 0.001

DEFINE_SOURCE(xmom, c, t, dS, eqn)

{

real Z[ND_ND];

real a1,a2,a3,X,Y,source;

C_CENTROID(Z, c, t);

X=Z[0];

Y=Z[1];

a1=((X-X0)-(Y-Y0))/Y;

a2=pow(a1,2.0);

a3=Betax*(Y-Y0)*(Y-Y0);

source=Fx0*Phi4*exp(-a2-a3);

return source;

}

DEFINE_SOURCE(ymom, c, t, dS, eqn)

{

real Z[ND_ND];

real a1,a2,a3,X,Y,source;

C_CENTROID(Z, c, t);

X=Z[0];

Y=Z[1];

a1=(X-X0)/Y;

a2=pow(a1,2.0);

a3=Betay*(Y-Y0)*(Y-Y0);

source=Fy0*Phi4*exp(-a2-a3);

return source;

}