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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » udf代码求助
北京石油化工学院2026年研究生招生接收调剂公告
52/1返回列表
查看: 1540  |  回复: 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

新虫 (小有名气)
引用回帖:
12楼: Originally posted by CFD_light at 2015-05-26 11:03:36
case发散？udf通过了？...

出现Divergence detected in AMG solver: species-0的错误，udf应该是通过了，不过有这样的warning：temp definition shadows previous definition
temp definition shadows previous definition
temp definition shadows previous definition
temp definition shadows previous definition
temp definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
mw definition shadows previous definition
yi definition shadows previous definition
肯定是哪里有问题了，不过我真是小白。
赞一下 回复此楼
13楼2015-05-26 14:02:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

huangkai1

新虫 (小有名气)
真诚期待得到大家的帮助指导，快不能毕业了。
赞一下 回复此楼
2楼2015-05-25 16:58:03
已阅   回复此楼   关注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的回帖
查看全部 21 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856，材料与化工321分求调剂 +4 大馋小子 2026-03-27 4/200 2026-03-27 13:07 by yhmsz
[考研] 266求调剂 +7 阳阳哇塞 2026-03-27 8/400 2026-03-27 13:04 by zzll406
[考研] 275求调剂 +5 jjjjjjjjjjl 2026-03-27 5/250 2026-03-27 12:52 by lature00
[考研] 329求调剂 +4 星野? 2026-03-26 4/200 2026-03-27 12:00 by 不吃魚的貓
[考研] 一志愿211，335分，0856，求调剂院校和导师 +4 倾____萧 2026-03-27 5/250 2026-03-27 11:52 by zhshch
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[考研] 359求调剂 +4 王了个楠 2026-03-25 4/200 2026-03-27 08:43 by 不吃魚的貓
[考研] 349求调剂 +4 李木子啊哈哈 2026-03-25 4/200 2026-03-26 22:49 by fmesaito
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +5 星沉uu 2026-03-26 6/300 2026-03-26 19:02 by macy2011
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 0703化学调剂，求导师收 +7 天天好运来上岸� 2026-03-24 7/350 2026-03-24 20:26 by peike
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 277分求调剂，跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
信息提示
关闭
请填处理意见
关闭