CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

UDF_Erosion

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 25, 2018, 08:33
Default UDF_Erosion
  #1
New Member
 
Giuseppe
Join Date: Jan 2013
Location: Italy
Posts: 15
Rep Power: 13
pezarello86 is on a distinguished road
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)
}
pezarello86 is offline   Reply With Quote

Reply

Tags
dpm fluent, erosion model, udf

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 17:23.