[parth_shah91]

We are trying to simulate 2D vortices in Ansys fluent using UDF.
we have written the code but it is not giving the desired output as shown in the given paper written by researchers from china.

I have written the following code in Define_init.
Please let me know where I am wrong

" #include "udf.h"

DEFINE_INIT(my_init_func, domain)
{
cell_t c;
Thread *t;
real xc[ND_ND];

/* loop over all cell threads in the domain */
thread_loop_c (t,domain)
{

/* loop over all cells */
begin_c_loop_all (c,t)
{
C_CENTROID(xc,c,t);
C_U(c,t) = (9.71*(sqrt(ND_SUM(pow(xc[0],2.),pow(xc[1],2.))) ))/ (ND_SUM(pow(xc[0] ,2.),pow(xc[1] ,2.))+0.36);

}
end_c_loop_all (c,t)
}
}
"

Please identify any errors in the code and suggest any changes. Any help will be greatly appreciated