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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent模拟流化床甲烷干重整反应udf使用浮点溢出
11/1返回列表
查看: 704  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者1173018735将赠送您 5 个金币

1173018735

新虫 (初入文坛)

[求助] fluent模拟流化床甲烷干重整反应udf使用浮点溢出

我在复现一篇用欧拉双流体方法复现流化床甲烷干重整反应的论文，我找不出来有什么问题，udf中公式的表达是没问题的，这个模型不使用UDF可以正常计算，但是一使用了UDF就会报错。而这个UDF在另外一个简单的模型里使用没有问题。有大佬帮我看看什么问题吗？
# include "udf.h"
# include "mem.h"

DEFINE_VR_RATE(rate, c, t, r, wk, yk, rate, rr_t)
{

#define UGC 8.314

        real T = C_T(c, t);
        real k1 = 1.29e6 * exp(-102065.0 / UGC / C_T(c, t));
        real k2 = 0.35e6 * exp(-81030.0 / UGC / C_T(c, t));
        real k3 = 6.95e3 * exp(-58893.0 / UGC / C_T(c, t));
        real k4 = 5.55e9 * exp(-166397.0 / UGC / C_T(c, t));
        real k5 = 1.34e15 * exp(-243835.0 / UGC / C_T(c, t));
        real Kco2_1 = 2.61e-2 * exp(37641.0 / UGC / C_T(c, t));
        real Kch4_1 = 2.60e-2 * exp(40684.0 / UGC / C_T(c, t));
        real Kco2_2 = 0.5771 * exp(9262.0 / UGC / C_T(c, t));
        real Kh2_2 = 1.494 * exp(6025.0 / UGC / C_T(c, t));
        real Kch4_3 = 0.21 * exp(-567.0 / UGC / C_T(c, t));
        real Kh2_3 = 5.18e7 * exp(-133210.0 / UGC / C_T(c, t));
        real Kh2_4 = 1.83e13 * exp(-216145.0 / UGC / C_T(c, t));
        real Kh2o_4 = 4.73e-6 * exp(97770.0 / UGC / C_T(c, t));
        real Kch4_4 = 3.49 * exp(-0 / UGC / C_T(c, t));
        real Kco_5 = 7.34e-6 * exp(100395.0 / UGC / C_T(c, t));
        real Kco2_5 = 2.81e7 * exp(-104085.0 / UGC / C_T(c, t));
        real Kp1 = 6.78e14 * exp(-259660.0 / UGC / C_T(c, t));
        real Kp2 = 56.4971 * exp(-36580.0 / UGC / C_T(c, t));
        real Kp3 = 2.95e5 * exp(-84400.0 / UGC / C_T(c, t));
        real Kp4 = 1.3827e7 * exp(-125916.0 / UGC / C_T(c, t));
        real Kp5 = 1.9393e9 * exp(-168527.0 / UGC / C_T(c, t));


        real mCO = C_YI(c, t, 0);
        real mH2 = C_YI(c, t, 1);
        real mH2O = C_YI(c, t, 2);
        real mCH4 = C_YI(c, t, 3);
        real mCO2 = C_YI(c, t, 4);

        real pOp = 101325;
        real sum = (mCH4 / 16 + mCO2 / 44 + mCO / 28 + mH2 / 2 + mH2O / 18);
        real PCH4 = mCH4 / 16 / sum * pOp;
        real PCO2 = mCO2 / 44 / sum * pOp;
        real PCO = mCO / 28 / sum * pOp;
        real PH2 = mH2 / 2 / sum * pOp;
        real PH2O = mH2O / 18 / sum * pOp;

if (!strcmp(r->name, "reaction-1")
    {
        *rate = k1 * Kco2_1 * Kch4_1 * (PCH4 / 101325) * (PCO2 / 101325) / pow((1 + Kco2_1 * (PCO2 / 101325) + Kch4_1 * (PCH4 / 101325)), 2) * (1 - pow(((PCO / 101325) * (PH2 / 101325)), 2) / (Kp1 * (PCH4 / 101325) * (PCO2 / 101325)));
        C_UDMI(c, t, 0) = *rate;
    }

    else if (!strcmp(r->name, "reaction-2")
    {
        *rate = k2 * Kco2_2 * Kh2_2 * (PCO2 / 101325) * (PH2 / 101325) / pow((1 + Kco2_2 * (PCO2 / 101325) + Kh2_2 * (PH2 / 101325)), 2) * (1 - (PCO / 101325) * (PH2O / 101325) / (Kp2 * (PCO2 / 101325) * (PH2 / 101325)));
        C_UDMI(c, t, 1) = *rate;
    }

    else if (!strcmp(r->name, "reaction-3")
    {
        *rate = k3 * Kch4_3 * ((PCH4 / 101325) - (pow((PH2 / 101325), 2) / Kp3)) / pow((1 + Kch4_3 * (PCH4 / 101325) + 1 / Kh2_3 * pow((PH2 / 101325), 1.5)), 2);
        C_UDMI(c, t, 2) = *rate;
    }

    else if (!strcmp(r->name, "reaction-4")
    {
        *rate = k4 / Kh2o_4 * ((PH2O / 101325) / (PH2 / 101325) - (PCO / 101325) / Kp4) / pow((1 + Kch4_4 * (PCH4 / 101325) + ((PH2O / 101325) / (PH2 / 101325)) / Kh2o_4 + pow((PH2 / 101325), 1.5) / Kh2_4), 2);
        C_UDMI(c, t, 3) = *rate;
    }

    else if (!strcmp(r->name, "reaction-5")
    {
        *rate = k5 / Kco_5 / Kco2_5 * ((PCO2 / 101325) / (PCO / 101325) - (PCO / 101325) / Kp5) / pow((1 + Kco_5 * (PCO / 101325) + ((PCO2 / 101325) / (PCO / 101325)) / (Kco_5 * Kco2_5)), 2);
        C_UDMI(c, t, 4) = *rate;
    }
    *rr_t = *rate;
}@wuming524
回复此楼

» 猜你喜欢
1楼 2024-05-01 15:48:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1173018735 的主题更新
11/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自科面上基金字体 +7 iwuli 2026-03-12 8/400 2026-03-19 08:55 by sjtu2012
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 11/550 2026-03-19 08:11 by 脱颖而出
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +10 yangfz 2026-03-17 10/500 2026-03-18 20:14 by walc
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考博] 环境领域全国重点实验室招收博士1-2名 +3 QGZDSYS 2026-03-13 5/250 2026-03-18 11:13 by QGZDSYS
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭