| 查看: 1495 | 回复: 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 ] |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有48人回复
-
FileCode能看出啥？
已经有17人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有6人回复
-
filecode
已经有9人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有23人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有6人回复
-
为什么网上很多人说本周 12号出结果
已经有5人回复
-
帮忙看看fileCode
已经有11人回复
-
基金中了
已经有15人回复
2楼2017-11-27 11:09:35
100