[pezarello86]

hello everybody,
i am trying to simulate erosion in a pipe with more elbows and i found erosion values very high; the mesh is such that the yplus is in compliance with the turbulence model.
the erosion model used is that of Finnie (implemented in Fluent) but my idea is to use that of tabakoff by mean a UDF.
could you someone help me? i have a udf but does not work.
#include "udf.h"
DEFINE_DPM_EROSION(Tabakoff_Grant, tp, t, f, normal, alpha, Vmag, mdot)
{
#define pi 3.14159
Thread *t0;
cell_t c0;

real K1=5.42e-8; /*constant*/
real K12=0.293328; /*constant*/
real K3=6.44e-14; /*constant*/
real K4=0.00525;
real K2;
real alpha0Deg, alpha0;
real V1,V2,V3;
real E,Rt,f_Vpn,Vpn,f_alpha, erosion,teta;
real A[ND_ND], area;

V1=1/(pow(K1,0.5));
V2=1/(pow(K3,1/4));
V3=1/K4;
alpha0Deg=30;

begin_f_loop(f, t)
{
F_AREA(A,f,t);
area = NV_MAG(A);

alpha0=alpha0Deg*pi/180;
Rt=1-(Vmag/V3)*sin(alpha);
Vpn=(Vmag/V2)*sin(alpha);
f_Vpn=K3*pow(Vpn,4);
if(alpha<2*alpha0)
K2=1;
else
K2=0;
teta=(alpha/alpha0)*pi/2;
f_alpha=pow((1+K2*K12*sin(teta)),2);
erosion=K1*f_alpha*pow(Vmag,2)*pow(cos(alpha),2)*( 1-pow(Rt,2))+f_Vpn;
F_STORAGE_R(f,t,SV_DPMS_EROSION)=erosion;
}
end_f_loop(f, t)
}