udf有问题,你试试这个
#include "udf.h"
DEFINE_PROFILE(inlet_velocity, t, i)
{
real x[ND_ND],y,vel;
face_t f;
begin_f_loop(f,t)
{
F_CENTROID(x, f, t);
y = x[1];
if(y>0.0393)
vel = 29.5;
else
vel = 29.5 * pow((y / 0.0393), 1/7);
F_PROFILE(f, t, i) = vel;
}
end_f_loop(f,t)
}