请问UDF给定速度型入口问题
-
自带速度型入口算例中,经过编译通过了,在速度入口边界条件中选择,计算发现仅在第一层网格内执行了,第二层网格已经是来流速度了,单位、网格均没有问题。
udf: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 = UMEAN*(B+1.); 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)