UDF編寫問題
本人剛剛接觸fluent自定義函數(shù),對函數(shù)的編寫理解還很淺顯,寫了以下一段程序作為速度入口邊界條件,原風(fēng)速廓線方程為u=Uref*(z/Zref)^0.22,其中Uref=5m/s,Zref=0.5m,但是導(dǎo)入fluent以后,第一步迭代就出現(xiàn)divergence detected in AMG solver:x-momentum的錯誤提示,之前不用自定義函數(shù)時,計算是可以收斂的,請問這是什么原因?程序編寫有錯誤嗎?小弟先謝過了.
#include "udf.h"
#define a 0.22
#define Uref 5
#define Zref 0.5
DEFINE_PROFILE(x_velocity,thread,position)
{
real z,u,x[ND_ND];
face_t f;
begin_f_loop(f,thread)
{
F_CENTROID(x,f,thread);
z=x[1];
u=Uref*pow(z/Zref,a);
F_PROFILE(f,thread,position)=u;
}
end_f_loop(f,thread)
}