[PierceH]

Hello,


I'm trying to add an actuator disk into my simulation which involves creating an X-momentum source code. In order to create a Gaussian distribution over every cell (p) where the actuator disk lies I've created a for loop however upon implementing the adjust function to fluent the solution will not iterate once after 5 minutes (Fluent is not crashing).
Any suggestions?



See the code here


DEFINE_ADJUST(adjust1,d)
{
Thread *t;
cell_t c;
real xc[ND_ND];
real variable;
thread_loop_c (t,d)
{
begin_c_loop (c,t)

C_CENTROID(xc,c,t);

xp = 0.0125;
epsilon = 0.052632;
sum = 0;
delta = 1;
n=38;

for (i=1 ; i=n ; 1)
{
yp = -10-(i/38);

r = pow((xc[0]-xp),2) + pow((xc[1]-yp),2);
term1 = exp(-pow(r/epsilon,2));
eta = term1*(1/M_PI*epsilon*epsilon*delta);
sum = sum + eta;
}

C_UDSI(c,t,0) = sum*500000;
C_UDSI(c,t,1) = r;

end_c_loop(c,t)
}

[AlexanderZ]

check this condition

for (i=1 ; i=n ; 1)

use C language