24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1410  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856,材料与化工321分求调剂 +3 大馋小子 2026-03-27 3/150 2026-03-27 12:43 by wxiongid
[考研] 329求调剂 +4 星野? 2026-03-26 4/200 2026-03-27 12:00 by 不吃魚的貓
[考研] 求调剂 +3 刘柯@ 2026-03-24 4/200 2026-03-27 11:28 by shangxh
[论文投稿] Journal of Mechanical Science and Technology +3 Russ_ss 2026-03-25 5/250 2026-03-27 10:49 by 陆小果画大饼
[考研] 求调剂 +5 芦lty 2026-03-25 6/300 2026-03-26 20:49 by 不吃魚的貓
[考研] 303求调剂 +7 元夕元 2026-03-20 8/400 2026-03-26 20:38 by 不吃魚的貓
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 求调剂 +8 Auroracx 2026-03-22 8/400 2026-03-26 19:55 by 不吃魚的貓
[考研] 340求调剂 +3 Amber00 2026-03-26 3/150 2026-03-26 18:57 by 不吃魚的貓
[考研] 中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培 +5 YangTyu1 2026-03-26 5/250 2026-03-26 18:27 by 猫咪猫咪呀
[考研] 085602 289分求调剂 +8 WWW西西弗斯 2026-03-24 8/400 2026-03-26 16:33 by 不吃魚的貓
[考研] 考研一志愿苏州大学初始315(英一)求调剂 +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
信息提示
请填处理意见