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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » udf代码求助
北京石油化工学院2026年研究生招生接收调剂公告
211/3返回列表123下一页
查看: 1542  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

huangkai1

新虫 (小有名气)

[求助] 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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-05-25 16:57:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangkai1

新虫 (小有名气)
真诚期待得到大家的帮助指导，快不能毕业了。
赞一下 回复此楼
2楼2015-05-25 16:58:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CFD_light

铁虫 (正式写手)
我们可以帮助
回复此楼
模拟的最高境界：同形不同理
3楼2015-05-25 17:44:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nickma

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +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

新虫 (小有名气)
引用回帖:
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的回帖

gyctju

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
mw数组赋值了吗？
赞一下 回复此楼
6楼2015-05-25 20:44:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nickma

至尊木虫 (著名写手)
调试了下好像可以，方法见站短
赞一下 回复此楼
7楼2015-05-25 21:06:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangkai1

新虫 (小有名气)
引用回帖:
6楼: Originally posted by gyctju at 2015-05-25 20:44:20
mw数组赋值了吗？

没有，估计可能问题在这里，我以为是直接传递求解器的数值，谢谢你，祝您事业一帆风顺！

[ 发自小木虫客户端 ]
赞一下 回复此楼
8楼2015-05-25 23:18:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CFD_light

铁虫 (正式写手)
解决了？
回复此楼
模拟的最高境界：同形不同理
9楼2015-05-26 10:36:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangkai1

新虫 (小有名气)
引用回帖:
9楼: Originally posted by CFD_light at 2015-05-26 10:36:18
解决了？

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

[ 发自小木虫客户端 ]
赞一下 回复此楼
10楼2015-05-26 10:55:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huangkai1 的主题更新
211/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
[考研] 南昌大学材料专硕311分求调剂 +10 77chaselx 2026-03-20 10/500 2026-03-27 11:28 by wangjy2002
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
[考研] 一志愿华理，数一英一285求A区调剂 +8 AZMK 2026-03-25 10/500 2026-03-26 22:37 by 学员8dgXkO
[考研] 321求调剂 +5 材料cailiao 2026-03-21 5/250 2026-03-26 20:41 by fmesaito
[考研] 289求调剂 +17 硕星赴 2026-03-23 17/850 2026-03-26 16:18 by 不吃魚的貓
[考研] 279 分 求调剂 +3 睡个好觉_16 2026-03-24 3/150 2026-03-26 15:43 by zzll406
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 考研调剂 +6 来好运来来来 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[考研] 299求调剂 +7 某某某某位 2026-03-21 8/400 2026-03-25 20:34 by 热情沙漠
[考研] 0854AI CV方向招收调剂 +4 章小鱼567 2026-03-23 4/200 2026-03-25 17:04 by CoderLoser
[考研] B区考研调剂 +4 yqdszhdap－ 2026-03-22 5/250 2026-03-25 08:51 by baoball
[考研] 工科0856求调剂 +5 沐析汀汀 2026-03-21 5/250 2026-03-23 17:56 by 海瑟薇-
[考研] 一志愿东华大学化学070300，求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
信息提示
关闭
请填处理意见
关闭