24小时热门版块排行榜

返回列表

steelbody

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 107.5
帖子: 20
在线: 7.2小时
虫号: 4641502
注册: 2016-04-29
专业: 流体力学

[求助] 载入UDF之后无法初始化，求助高手

载入速度，湍动能以及湍流耗散率的udf后无法初始化，求教，下面附上UDF，谢谢
#include "udf.h"
#define YMIN 0.0 /* constants */
#define YMAX 0.17
#define UMEAN 1.0
#define B 1./7.
#define DELOVRH 0.5
#define VISC 1.7894e-05
#define CMU 0.09
#define VKC 0.41
/* profile for x-velocity */
DEFINE_PROFILE(x_velocity, t, i)
{
real y, del, h, x[ND_ND], ufree; /* variable declarations */
face_t f;
h = YMAX - YMIN;
del = DELOVRH*h;
ufree = 2.147;
begin_f_loop(f, t)
{
F_CENTROID(x,f,t);
y = x[1];
if (y <= del)
F_PROFILE(f,t,i) = ufree*pow(y/del,B);
else
F_PROFILE(f,t,i) = ufree*pow((h-y)/del,B);
}
end_f_loop(f, t)
}
/* profile for kinetic energy */
DEFINE_PROFILE(k_profile, t, i)
{
real y, del, h, ufree, x[ND_ND];
real ff, utau, knw, kinf;
face_t f;
h = YMAX - YMIN;
del = DELOVRH*h;
ufree = 2.147;
ff = 0.045/pow(ufree*del/VISC,0.25);
utau=sqrt(ff*pow(ufree,2.)/2.0);
knw=pow(utau,2.)/sqrt(CMU);
kinf=0.002*pow(ufree,2.);
begin_f_loop(f, t)
{
F_CENTROID(x,f,t);
y=x[1];
if (y <= del)
F_PROFILE(f,t,i)=knw+y/del*(kinf-knw);
else
F_PROFILE(f,t,i)=knw+(h-y)/del*(kinf-knw);
}
end_f_loop(f, t)
}
/* profile for dissipation rate */
DEFINE_PROFILE(dissip_profile, t, i)
{
real y, x[ND_ND], del, h, ufree;
real ff, utau, knw, kinf;
real mix, kay;
face_t f;
h = YMAX - YMIN;
del = DELOVRH*h;
ufree = 2.147;
ff = 0.045/pow(ufree*del/VISC,0.25);
utau=sqrt(ff*pow(ufree,2.)/2.0);
knw=pow(utau,2.)/sqrt(CMU);
kinf=0.002*pow(ufree,2.);
begin_f_loop(f, t)
{
F_CENTROID(x,f,t);
y=x[1];
if (y <= del)
kay=knw+y/del*(kinf-knw);
else
kay=knw+(h-y)/del*(kinf-knw);
if (VKC*y < 0.085*del)
mix = VKC*y;
else
mix = 0.085*del;
F_PROFILE(f,t,i)=pow(CMU,0.75)*pow(kay,1.5)/mix;
}
end_f_loop(f,t)
}
出现问题如下

MBZ0BDJ)6NIO09USKTKQ[$9.png