24小时热门版块排行榜    

查看: 1541  |  回复: 1

xgp198749

新虫 (初入文坛)

[求助] 求助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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangxueqi

新虫 (初入文坛)

【答案】应助回帖

求问楼主怎么解决的?感谢!
2楼2017-11-27 11:09:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xgp198749 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 BMTdnp4OJO9f 2026-11-06 4/200 2026-11-07 19:03 by Xc2FxHsHkBNc
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 BMTdnp4OJO9f 2026-11-06 5/250 2026-11-07 17:56 by Xc2FxHsHkBNc
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 BMTdnp4OJO9f 2026-11-06 5/250 2026-11-07 17:45 by Xc2FxHsHkBNc
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 5sJFDKp2kwuI 2026-11-06 4/200 2026-11-07 16:50 by Xc2FxHsHkBNc
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:17 by Xc2FxHsHkBNc
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:06 by Xc2FxHsHkBNc
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见