版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huangkai1

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 22.1
散金: 200
帖子: 173
在线: 166.1小时
虫号: 2694682
注册: 2013-10-01
专业: 天体力学和人造卫星动力学

[求助] udf代码求助已有3人参与

各位老师同学好，最近在在做毕设，遇到一个问题实在不知如何解决，特来求助，先谢谢各位的帮助。

我做的主要是固定床反应器的模拟，在不加载udf源项时，能够得到模拟结果，但当加载源项后，就出现Sample Text的错误，现附上自己的cas和代码，希望大家能够指导，献上所有金币。

#include "udf.h"
/*动力学参数常量*/
#define PRE_EXP1 865.609    /*指前因子单位mol/g.h.kpa2*/
#define PRE_EXP2 6275.871    /*指前因子单位mol/g.h.kpa2*/
#define PRE_EXP3 1.60375e+05  /*指前因子单位mol/g.h.kpa*/
#define PRE_EXP4 3277.5722 /*指前因子单位mol/g.h.kpa*/
#define PRE_EXP5 13481.3718 /*指前因子单位mol/g.h.kpa*/

#define ACTIVE1  6.04e+04 /*活化能 J/mol*/
#define ACTIVE2  6.24e+04 /*活化能 J/mol*/
#define ACTIVE3  6.66e+04 /*活化能 J/mol*/
#define ACTIVE4  2.64e+04 /*活化能 J/mol*/
#define ACTIVE5  5.42e+04 /*活化能 J/mol*/

#define BETA 0.0       /*温度指数*/

double arrhenius_rate1(real temp)
{
return PRE_EXP1*pow(temp,BETA)*exp(-ACTIVE1/(UNIVERSAL_GAS_CONSTANT*temp));
}
double arrhenius_rate2(real temp)
{
return PRE_EXP2*pow(temp,BETA)*exp(-ACTIVE2/(UNIVERSAL_GAS_CONSTANT*temp));
}
double arrhenius_rate3(real temp)
{
return PRE_EXP3*pow(temp,BETA)*exp(-ACTIVE3/(UNIVERSAL_GAS_CONSTANT*temp));
}
double arrhenius_rate4(real temp)
{
return PRE_EXP4*pow(temp,BETA)*exp(-ACTIVE4/(UNIVERSAL_GAS_CONSTANT*temp));
}
double arrhenius_rate5(real temp)
{
return PRE_EXP5*pow(temp,BETA)*exp(-ACTIVE5/(UNIVERSAL_GAS_CONSTANT*temp));
}

#define P_C8H10 0
#define O_C8H10 1
#define M_C8H10 2
#define C6H6    3
#define C2H4    4
#define H20       5
#define CH3OH  6
#define C7H8 7

double reaction_rate1(cell_t c, Thread *cthread,real mw[],real yi[])
{
real concenC7H8 = C_R(c,cthread)*yi[C7H8]/mw[C7H8];
real concenCH3OH = C_R(c,cthread)*yi[CH3OH]/mw[CH3OH];
return arrhenius_rate1(C_T(c,cthread))*pow(concenC7H8,1)*pow(concenCH3OH,1)*pow((UNIVERSAL_GAS_CONSTANT*C_T(c,cthread)),2);
}

double reaction_rate2(cell_t c, Thread *cthread,real mw[],real yi[])
{
real concenCH3OH = C_R(c,cthread)*yi[CH3OH]/mw[CH3OH];
return arrhenius_rate2(C_T(c,cthread))*pow(concenCH3OH,2)*pow((UNIVERSAL_GAS_CONSTANT*C_T(c,cthread)),2);
}

double reaction_rate3(cell_t c, Thread *cthread,real mw[],real yi[])
{
real concenC7H8 = C_R(c,cthread)*yi[C7H8]/mw[C7H8];
return arrhenius_rate3(C_T(c,cthread))*pow(concenC7H8,1)*pow((UNIVERSAL_GAS_CONSTANT*C_T(c,cthread)),1);
}

double reaction_rate4(cell_t c, Thread *cthread,real mw[],real yi[])
{
real concenP_C8H10 = C_R(c,cthread)*yi[P_C8H10]/mw[P_C8H10];
return arrhenius_rate4(C_T(c,cthread))*pow(concenP_C8H10,1)*pow((UNIVERSAL_GAS_CONSTANT*C_T(c,cthread)),1);
}

double reaction_rate5(cell_t c, Thread *cthread,real mw[],real yi[])
{
real concenP_C8H10 = C_R(c,cthread)*yi[P_C8H10]/mw[P_C8H10];
return arrhenius_rate5(C_T(c,cthread))*pow(concenP_C8H10,1)*pow((UNIVERSAL_GAS_CONSTANT*C_T(c,cthread)),1);
}

/*定义能量方程源项*/
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
real delt_h1,delt_h2,delt_h3,delt_h4,delt_h5, source,mw[8],yi[8];
/*下面是五个反应的生成焓单位kJ/mol*/
delt_h1=-73.682;
delt_h2=-14.430;
delt_h3=-47.772;
delt_h4=59.252;
delt_h5=0.862;
source=(delt_h1*reaction_rate1(c,t,mw,yi)+delt_h2*reaction_rate2(c,t,mw,yi)+delt_h3*reaction_rate3(c,t,mw,yi)+delt_h4*reaction_rate4(c,t,mw,yi)+delt_h5*reaction_rate5(c,t,mw,yi))*1000*235.785;/*计算所有反应的反应热单位kj/h*/
dS=0;
return source;
}

/*定义C7H8组分方程源项*/
DEFINE_SOURCE(c7h8_source,c,t,dS,eqn)
{
      real source,mw[8],yi[8];
      source=(0-reaction_rate1(c,t,mw,yi)-reaction_rate3(c,t,mw,yi)+reaction_rate4(c,t,mw,yi))*mw[C7H8]*0.2385785;       /*计算甲苯的消耗速度单位g/h)*/
      dS=0;
return source;
}
/*定义CH3OH组分方程源项*/
DEFINE_SOURCE(ch3oh_source,c,t,dS,eqn)
{
      real source,mw[8],yi[8];
      source=(-reaction_rate1(c,t,mw,yi)-reaction_rate2(c,t,mw,yi))*mw[CH3OH]*0.235785;                               /*计算CH30H的消耗速度单位g/h*/
      dS=0;
return source;
}
/*定义C6H6组分方程源项*/
DEFINE_SOURCE(C6H6,c,t,dS,eqn)
{
      real source,mw[8],yi[8];
      source=0.5*(reaction_rate3(c,t,mw,yi))*mw[C6H6]*0.235785;                                                       /*计算CO2的生成速度单位g/h*/
      dS=0;
return source;
}

/*定义H20组分方程源项*/
DEFINE_SOURCE(h2o_source,c,t,dS,eqn)
{
      real source,mw[8],yi[8];
      source=(reaction_rate1(c,t,mw,yi)+reaction_rate2(c,t,mw,yi))*mw[H2O]*0.235785;/*计算H2O的生成速度单位g/h*/
      dS=0;
return source;
}

/*定义P-X组分方程源项*/
DEFINE_SOURCE(P_C810source,c,t,dS,eqn)
{
      real source,mw[8],yi[8];
      source=(reaction_rate1(c,t,mw,yi)-reaction_rate4(c,t,mw,yi)-reaction_rate5(c,t,mw,yi))*mw[P_C8H10]*0.235785; /*计算P-X的生成速度单位g/h*/
      dS=0;
return source;
}

/*定义O-X组分方程源项*/
DEFINE_SOURCE(O_C8H10source,c,t,dS,eqn)
{
      real source,mw[8],yi[8];
      source=0.5*reaction_rate5(c,t,mw,yi)*mw[O_C8H10]*0.235785; /*计算O-X的生成速度单位g/h*/
      dS=0;
return source;
}

/*定义M-X组分方程源项*/
DEFINE_SOURCE(M_C8H10source,c,t,dS,eqn)
{
      real source,mw[8],yi[8];
      source=0.5*reaction_rate5(c,t,mw,yi)*mw[M_C8H10]*0.235785; /*计算M-X的生成速度单位g/h*/
      dS=0;
return source;
}
/*定义GH组分方程源项*/
DEFINE_SOURCE(GHsource,c,t,dS,eqn)
{
      real source,mw[8],yi[8];
      source=0.5*(reaction_rate2(c,t,mw,yi)+reaction_rate4(c,t,mw,yi))*mw[C2H4]*0.235785; /*计算M-X的生成速度单位g/h*/
      dS=0;
return source;
}

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : fanyinqimoxing.cas

2015-05-25 16:55:43, 735.67 K

» 猜你喜欢

291求调剂已经有12人回复
317求调剂已经有5人回复
0703化学求调剂，各位老师看看我！！！已经有3人回复
339求调剂已经有5人回复
287求调剂已经有8人回复
一志愿北京化工大学 070300 学硕 336分求调剂已经有8人回复
085602 化工专硕 338分求调剂已经有5人回复
0856，材料与化工321分求调剂已经有4人回复
一志愿北京化工大学材料与化工 264分各科过A区国家线已经有10人回复
266求调剂已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

Fluent中编写UDF之后无法初始化。。。已经有5人回复
解释型UDF如何调用温度梯度已经有9人回复
udf程序改为并行已经有4人回复
动网格的UDF问题已经有21人回复
udf的编译出现的问题该怎么修改已经有12人回复
FLUENT中编写UDF时指针变量是否要初始化已经有4人回复
flunet UDF 求助已经有6人回复
udf添加能量源项已经有3人回复
fluent udf的应用已经有20人回复
Fluent中UDF的编译问题已经有25人回复
fluent UDF代码调用预定义程序的问题（如fortran程序）已经有15人回复
求助udf的输入输出已经有16人回复
【求助】FLUENT UDF 气液两相流汽蚀汽化潜热(问题详见内容) 已经有21人回复
边界条件UDF求助已经有4人回复
UDF源项的一个疑问已经有14人回复
求助：UDF模型的编译问题已经有15人回复
【求助】求UDF高手帮忙看看我的问题已经有18人回复
【讨论】请问大家用UDF的时候，有没有碰到如下代码，很重要的已经有7人回复
【求助】冷凝的UDF程序已经有15人回复

1楼 2015-05-25 16:57:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

huangkai1

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 22.1
散金: 200
帖子: 173
在线: 166.1小时
虫号: 2694682
注册: 2013-10-01
专业: 天体力学和人造卫星动力学

引用回帖:

9楼: Originally posted by CFD_light at 2015-05-26 10:36:18
解决了？

没有啊，还是有问题，真的要挂了，又调了一个上午还是说发散……………谢谢各位的帮助

[ 发自小木虫客户端 ]

赞一下

回复此楼

10楼2015-05-26 10:55:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 21 个回答

huangkai1

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 22.1
散金: 200
帖子: 173
在线: 166.1小时
虫号: 2694682
注册: 2013-10-01
专业: 天体力学和人造卫星动力学

真诚期待得到大家的帮助指导，快不能毕业了。

赞一下

回复此楼

2楼2015-05-25 16:58:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

nickma

至尊木虫 (著名写手)

应助: 17 (小学生)
金币: 10791.7
散金: 1987
红花: 3
帖子: 1436
在线: 2320.6小时
虫号: 336250
注册: 2007-04-01
专业: 体育史

【答案】应助回帖

感谢参与，应助指数 +1

DEFINE_SOURCE(C6H6,c,t,dS,eqn)
改为下面名称：
DEFINE_SOURCE(C6H6_source,c,t,dS,eqn)
C6H6已经定为常量了

赞一下

回复此楼

4楼2015-05-25 18:26:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

huangkai1

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 22.1
散金: 200
帖子: 173
在线: 166.1小时
虫号: 2694682
注册: 2013-10-01
专业: 天体力学和人造卫星动力学

引用回帖:

4楼: Originally posted by nickma at 2015-05-25 18:26:28
DEFINE_SOURCE(C6H6,c,t,dS,eqn)
改为下面名称：
DEFINE_SOURCE(C6H6_source,c,t,dS,eqn)
C6H6已经定为常量了

这个我程序改了，没有什么用，不过还是要谢谢你。

赞一下

回复此楼

5楼2015-05-25 19:30:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 21 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +7	丹青奶盖 2026-03-26	8/400	2026-03-27 11:24 by 丹青奶盖
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +3	辰沐5211314 2026-03-26	3/150	2026-03-27 10:24 by 尽舜尧1
[考研] 316求调剂 +5	江辞666 2026-03-26	5/250	2026-03-27 08:08 by hypershenger
[考研] 304材料求调剂 +4	钟llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 317求调剂 +7	蛋黄咸肉粽 2026-03-26	7/350	2026-03-27 02:29 by fmesaito
[考研] 294分080500材料科学与工程求调剂 +4	柳溪边 2026-03-26	4/200	2026-03-26 21:14 by XPU李庆
[考研] 321求调剂 +6	wasdssaa 2026-03-26	6/300	2026-03-26 20:57 by sanrepian
[考研] 329求调剂 +5	1() 2026-03-22	5/250	2026-03-26 20:40 by fmesaito
[考研] 085602化学工程求调剂。 +4	平乐乐乐 2026-03-26	4/200	2026-03-26 17:57 by fmesaito
[考研] 生物技术与工程 +3	1294608413 2026-03-25	4/200	2026-03-25 18:02 by 1294608413
[考研] 【2026考研调剂】制药工程 284分求相关专业调剂名额 +4	袁奂奂 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 282求调剂 +3	wcq131415 2026-03-24	3/150	2026-03-25 12:16 by userper
[考研] 调剂 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[考博] 26申博自荐 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 材料/农业专业，07/08开头均可，过线就行 +3	呵唔哦豁 2026-03-23	4/200	2026-03-23 22:30 by 汪！？！
[考研] 276求调剂 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在风里
[考研] 一志愿东华大学化学070300，求调剂 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc
[考研] 一志愿华中科技大学071000，求调剂 +4	沿岸有贝壳6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 一志愿深大，0703化学，总分302，求调剂 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 学员8dgXkO
[考研] 材料 271求调剂 +5	展信悦_ 2026-03-21	5/250	2026-03-21 17:29 by 学员8dgXkO

24小时热门版块排行榜

huangkai1

[求助] udf代码求助 已有3人参与

» 本帖附件资源列表

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

huangkai1

huangkai1

nickma

【答案】应助回帖

huangkai1

[求助] udf代码求助已有3人参与