[google9002]

Hi,everyone
I wanna use a udf to define the wall shear on the boundary, this udf could be loaded and initialized right, but when I begin to calculate, Fluent tells me that Divergence detected in AMG solver: k -> Increasing relaxation sweeps!
When this udf unloaded, it works well, the solution is not so good though.
By the way, i use VOF model here.
And my code is below here, is there any problem with it?
I need your help...


#include <stdio.h>
#include "udf.h"
#include "mem.h"
#include <math.h>
DEFINE_PROFILE(wall_shear_x, thread, position)
{
real u, v, w, segm, dens,text,velo,bata;
face_t f;
cell_t c;
Thread *t0;
Domain *d;
d = Get_Domain(3);
int Zone_ID = 4;
thread = Lookup_Thread(d, Zone_ID);
begin_f_loop(f, thread)
{
u = F_U(f, thread);
v = F_V(f, thread);
w = F_W(f, thread);
segm = F_P(f, thread);
if (NNULLP(THREAD_STORAGE(thread, SV_DENSITY)))
dens = F_R(f, thread);
else
{
c = F_C0(f, thread);
t0 = F_C0_THREAD(f, thread);
dens = C_R(c, t0);
}
velo = sqrt(u*u + v*v + w*w);
bata = -u / velo;
text = segm*0.2 + 0.037731222*dens*velo*velo;
F_PROFILE(f, thread, position) = text*bata;
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(wall_shear_y, thread, position)
{
real u, v, w, segm, dens, text, velo, bata;
face_t f;
cell_t c;
Thread *t0;
Domain *d;
d = Get_Domain(3);
int Zone_ID = 4;
thread = Lookup_Thread(d, Zone_ID);
begin_f_loop(f, thread)
{
u = F_U(f, thread);
v = F_V(f, thread);
w = F_W(f, thread);
segm = F_P(f, thread);
if (NNULLP(THREAD_STORAGE(thread, SV_DENSITY)))
dens = F_R(f, thread);
else
{
c = F_C0(f, thread);
t0 = F_C0_THREAD(f, thread);
dens = C_R(c, t0);
}
velo = sqrt(u*u + v*v + w*w);
bata = -v / velo;
text = segm*0.2 + 0.037731222*dens*velo*velo;
F_PROFILE(f, thread, position) = text*bata;
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(wall_shear_z, thread, position)
{
real u, v, w, segm, dens, text, velo, bata;
face_t f;
cell_t c;
Thread *t0;
Domain *d;
d = Get_Domain(3);
int Zone_ID = 4;
thread = Lookup_Thread(d, Zone_ID);
begin_f_loop(f, thread)
{
u = F_U(f, thread);
v = F_V(f, thread);
w = F_W(f, thread);
segm = F_P(f, thread);
if (NNULLP(THREAD_STORAGE(thread, SV_DENSITY)))
dens = F_R(f, thread);
else
{
c = F_C0(f, thread);
t0 = F_C0_THREAD(f, thread);
dens = C_R(c, t0);
}
velo = sqrt(u*u + v*v + w*w);
bata = -w / velo;
text = segm*0.2 + 0.037731222*dens*velo*velo;
F_PROFILE(f, thread, position) = text*bata;
}
end_f_loop(f, thread)
}