我用的ansys fluent 14.0,自定義了一個(gè)換熱系數(shù),計(jì)算公式是下面那個(gè)公式,我編的udf是這樣的,上面那部分是換熱系數(shù)的,下面那部分是流體粘度的:
#include "udf.h"
#include "mem.h"
#define dp 0.02
#define cp 1561.70
#define y1 0.33
#define y2 0.60
DEFINE_PROFILE(heat_coe,t,i)
{
cell_t c;
real Pr;
real Re;
real x1,x2;
begin_c_loop(c,t)
{
Pr=C_MU_L(c,t)*cp/C_K_L(c,t);
Re=C_R(c,t)*C_V(c,t)*dp/C_MU_L(c,t);
x1=Pr;
x2=Re;
F_PROFILE(c,t,i)=
C_K_L(c,t)*(2.0+1.1*pow(x1,y1)*pow(x2,y2))/dp;
}
end_c_loop(c,t)
}
#include "udf.h"
DEFINE_PROPERTY(cell_viscosity,c,t)
{
real mu;
real tempr=C_T(c,t);
mu=exp((0-4.343)-2.0143*(log(tempr-273.15)-5.011));
return mu;
}
不知道有沒有錯(cuò)誤呢?因?yàn)閷?dǎo)進(jìn)fluent算的時(shí)候老是算若干步后就出現(xiàn)錯(cuò)誤,就是什么divergence detected in AMG solver:temperature。初始化的時(shí)候也有錯(cuò)誤。但是如果我把這個(gè)換熱系數(shù)設(shè)置為常數(shù)而不是我這個(gè)udf就不會有報(bào)錯(cuò),所以我想請各位朋友幫我看看我這個(gè)udf哪里出了錯(cuò)??
我是fluent新手,對udf更是生疏,雖然很簡單的一個(gè)udf,對于我來說確實(shí)不知道錯(cuò)誤在哪了啊~~跪求各位幫幫我吧~~~小女子感激不盡啊 ~~~
![]()
未命名4.jpg |