24小时热门版块排行榜    

查看: 1950  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

duoraeman

银虫 (小有名气)

[求助] UDF载入fluent之后出错,求指点!

本人不怎么会编程,需要整一个点源反应速率的UDF,但是导入fluent之后总是显示错误,下面是我写的程序的代码,求高手指点!出错的两行已经标出。


#include "udf.h"
#include

/*----------------------------arrhenius速率常数--------------------------------*/
#define A1 1.9e+7         
#define E1 23815
#define B  0.0

#define MCO 28            
#define MO2 36
         
#define O2   0
#define CO2  1
#define N2   2
#define NUM_SPECS 3

real arrhenius_rate(real temp)
{
return
A1*pow(temp,B)*exp(-E1/temp);
}

/*-----------------------------O2反应速率---------------------------------------*/
real YANGQI_rate(cell_t c, Thread *cthread,real mw[],real yi[])
{

  real concenO2 = C_R(c, cthread)*mw[0]/yi[0];

return arrhenius_rate(C_T(c,cthread))*concenO2;
}

/*------------------------------O2质量源项---------------------------------------*/
DEFINE_SOURCE(O2_source,c,t,dS,eqn)
{
    real source_O2;
    source_O2=0.0-MO2*YANGQI_rate/1000.0;   
    dS[eqn]=0.0;
    return source_O2;
}

/*------------------------------CO质量源项--------------------------------------*/
DEFINE_SOURCE(CO_source,c,t,dS,eqn)
{
    real source_CO;
    source_CO=2.0*MCO*YANGQI_rate/1000.0;   
    dS[eqn]=0.0;
    return source_CO;
}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoqiang3819

金虫 (小有名气)

【答案】应助回帖

duoraeman(金币+50): 钱财乃身外之物,拿去吧~ 2011-12-07 10:51:37
xiegangmai: 非应助请勿选择“应助回帖”。 2011-12-08 23:04:23
私底下私聊  金币给我
6楼2011-12-07 10:50:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

duoraeman

银虫 (小有名气)

是不是YANGQI_rate不能直接引用?应该怎么写啊?求指教~~~~
2楼2011-12-06 16:09:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elonlucy

银虫 (小有名气)

【答案】应助回帖


xiegangmai(金币+1): 谢谢参与! 2011-12-06 22:54:04
你在最开始定义的时候怎么定义氧气和二氧化碳的呢?
是不是需要在开始的时候先定义下呢?

» 本帖已获得的红花(最新10朵)

3楼2011-12-06 16:37:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duoraeman

银虫 (小有名气)

引用回帖:
3楼: Originally posted by elonlucy at 2011-12-06 16:37:47:
你在最开始定义的时候怎么定义氧气和二氧化碳的呢?
是不是需要在开始的时候先定义下呢?

你是说YANGQI_rate吗?这个在函数体内定义了 real YANGQI_rate也还是不行
4楼2011-12-06 17:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见