| 查看: 1512 | 回复: 1 | ||
[求助]
求助UDF 已有1人参与
|
我用fluent中的UDF编的程序在点编译（compile）的时候为什么显示下面的错误啊，“line 50: non-integer subscript expression: float.”第50行( m_0=C_UDSI(c,t,M0) 是一个宏调用啊，提示的意思是下标不是整数，而是浮点型，可是我定义的就不是整数啊，定义什么类型怎么还会有限制呢，太不理解了。多谢各位帮帮忙，我就这么点金币，都给了。 #include "udf.h" #include "mem.h" #include "sg_udms.h" #include "math.h" real pi=3.141593; real T=300.0; real KB=1.38e-23; /*Boltzmann constant*/ real rou_air=1.225; /*air density*/ real rou_particle=4.0e3; /*particle density*/ real m_0,m_1,m_2; /*moment variables*/ enum { M0, M1, M2, N_REQUIRED_UDS }; /*functions for calculating Mk*/ double m_k(real k) { return pow(m_0,(1./2.)*pow(k,2.)-(3./2.)*k+1.)*pow(m_1,-pow(k,2.)+2.*k)*pow(m_2,(1./2.)*pow(k,2.)-(1./2.)*k); } /*Make sure there are enough user defined-scalars*/ DEFINE_ADJUST(adjust,domain) { if(n_uds<N_REQUIRED_UDS) Internal_Error("not enough user-defined scalars allocated”）; } DEFINE_DIFFUSIVITY(Moment_diffusivity,c,t,i) { real D_mol; /*molecular diffusivity*/ real D_tur; /*turbulent diffusivity*/ real D_eff; /*effective diffusivity*/ real v_g; /*mean particle volume*/ real d_p; /*mean particle diameter*/ real f; /*friction coefficient*/ real c; /*slip corretion factor*/ real l=6.91e-8; /*mean free path of air*/ real mu=1.86e-05; /*air viscosity*/ real A1=1.257, A2=0.4, A3=0.55; real Sct=0.7; /*turbulent schmidt number*/ m_0=C_UDSI(c,t,M0); m_1=C_UDSI(c,t,M1); m_2=C_UDSI(c,t,M2); v_g=pow(m_1,2.)/(pow(m_0,3./2.)*pow(m_2,1./2.)); d_p=pow((6./pi)*v_g,1./3.); c=1.+2.*l*(A1+A2*exp(-A3*d_p/l))/d_p; f=3.*pi*mu*d_p/c; D_mol=KB*T/f; D_tur=C_MU_T(c,t)/rou_air*Sct; D_eff=D_mol+D_tur; return D_eff; } [ Last edited by xgp198749 on 2013-11-12 at 13:36 ] |
回复此楼» 猜你喜欢
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有11人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
面上函评意见出来了，像什么等级？
已经有14人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
-
面上意见出来了
已经有23人回复
2楼2017-11-27 11:09:35
投票:
10