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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent模拟堆积床内焦炭颗粒燃烧，出现错误提示“pressure correction”
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 3300  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

nwj0502

至尊木虫 (知名作家)

[求助] fluent模拟堆积床内焦炭颗粒燃烧，出现错误提示“pressure correction”

大家好，我想用fluent模拟堆积床内焦炭颗粒（直径0.005m）燃烧过程传热传质，采用气固双能量方程，组分方程。
遇到的问题：不加气固能量源项和气体组分方程，给定固体初始温度1000K，气体初始温度300K，模拟气固换热是没有问题的；但是加上气固能量源项后，点开始计算，给定气体温度1473K模拟点火，气体想固体传热，当固体温度达到着火点923K时（这时开始发生焦炭燃烧反应），fluent就会提示错误“pressure correction”，是不是能量源项太大？需要对源项进行线性化处理，具体怎么处理呢？还请fluent高手给予指点，小弟不胜感激。


模拟焦炭燃烧过程传热传质

焦炭燃烧反应：C+O2→CO2
焦炭燃点：923K
反应速率：Rc=3.14*（dc）^2*kc *CO2*Nc
                 式中，dc是焦炭颗粒直径，kc是综合反应速率常数，CO2是氧气摩尔浓度，Nc是单位体积焦炭颗粒数；
初始条件和边界条件：
Tg，0=1473K，Ts=340K，Co2,0= Co2，Co2,0= Cco2，Cn2,0= Cn2
点火温度1473K(Tg=1473K)，时间60s；保温温度1043K(Tg=1043K)，时间90s；接下来抽风温度300K(Tg=300K)（压力出口，-10KPa）

以下是我的部分udf代码：
real kfc;                               //传质系数,m/s；
    real Sh,Re,Sc;                           //舍伍德数，颗粒雷诺数，施密特准数；
    real mu_lam=C_MU_L(c,t);                //气体黏度；
    real den_gas=C_R(c,t);                    //气体密度，kg/m3；
    real Nc=9.92*pow(10.0,5.0);               //单位体积焦炭颗粒数，1/m3;
    real mole_total;                         //气体总的摩尔数
    real P_total;                            //气体总压,pa
    real P_O2;                             //氧气分压,pa
    real DO2;                              //氧气扩散系数,m2/s;
    real Co;                                //氧气摩尔浓度，mol/m3;
    real O2_mole_frac;                       //氧气摩尔分数；
    real kc1;                               //综合反应速率常数，m/s;
    real k1;                                //反应速率常数,m/s;
    real R1;                                //焦炭完全燃烧的反应速率，mol/m3s;

    P_total=ABS_P(C_P(c,t),op_pres);           //气体总压?
    mole_total=C_YI(c,t,0)/32.0+C_YI(c,t,1)/44.0+C_YI(c,t,2)/28.0; //
    O2_mole_frac=C_YI(c,t,0)/(32.0*mole_total);
    P_O2=P_total*O2_mole_frac;              //氧气分压
    Co=P_O2/R*C_T(c,t);

    DO2=2.00*pow(10.0,-4.0)*pow(0.001,1.78)/P_total;
    Sc=mu_lam/(den_gas*DO2);              //施密特准数舍伍德准数
    Re=ND_MAG(C_U(c,t),C_V(c,t),C_W(c,t))*DIAMETER*den_gas/mu_lam;   //雷诺数
    Sh=2.0/POROUSITY+0.75*pow(Re,0.5)*pow(Sc,0.333)/POROUSITY;      //舍伍德数
    kfc=Sh*DO2/DIAMETER;                                     //传质系数
    k1=6.53*pow(10.0,7.0)*pow(C_UDSI(c,t,0),0.5)*exp(-44000.0/(R*C_UDSI(c,t,0))); //反应速率常数
    kc1=(kfc+k1)/kfc*k1;

    if (C_UDSI(c,t,0)<923.0)

        R1=0;

    else

        R1=3.14*pow(DIAMETER,2.0)*kc1*Co*Nc;

return R1;       //
DEFINE_SOURCE(source_energy_gas,c,t,dS,eqn)
{
    real source;
    real delta_H=-2.79*pow(10.0,5.0);
    real a=0.5;  //反应热分配系数
    source=htc(c,t)*(C_UDSI(c,t,0)-C_T(c,t))-(1-a)*(delta_H*RC(c,t));
    dS[eqn]=0;
    return source;
}

DEFINE_SOURCE(source_energy_solid,c,t,dS,eqn)
{
    real source;
    real delta_H=-2.79*pow(10.0,5.0);
    real a=0.5;  //反应热分配系数
    source=htc(c,t)*(C_T(c,t)-C_UDSI(c,t,0))-a*(delta_H*RC(c,t));
    dS[eqn]=0;
    return source;
}

fluent模拟堆积床内焦炭颗粒燃烧，出现错误提示“pressure correction”
压力修正.jpg


fluent模拟堆积床内焦炭颗粒燃烧，出现错误提示“pressure correction”-1
出现错误.jpg
回复此楼

» 收录本帖的淘帖专辑推荐

模拟相关

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-16 00:00:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nwj0502

至尊木虫 (知名作家)
这么晚大家都睡了吧，自己顶一下吧，希望好心人来帮助，在这个问题上纠结两个礼拜了
赞一下 回复此楼
2楼2014-06-16 00:08:55
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nwj0502

至尊木虫 (知名作家)
怎么没人回复呢，热心的虫子帮帮呗

[ 发自小木虫客户端 ]
赞一下 回复此楼
3楼2014-06-16 08:53:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gearboy78

禁虫 (小有名气)
本帖内容被屏蔽
4楼2014-06-16 09:20:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nwj0502

至尊木虫 (知名作家)
引用回帖:
4楼: Originally posted by gearboy78 at 2014-06-16 09:20:43
具体调试是很费时间和精力的事情，只能自己去做，不要期望别人能帮你。就好像写一篇论文一样，你觉得会有人帮你写吗？

我不是要别人帮我调程序，因为刚学UDF不久，对UDF理解不深，我是想让大家看看程序，源项部分处理有什么不对的地方
赞一下 回复此楼
5楼2014-06-16 09:47:07
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gearboy78

禁虫 (小有名气)
本帖内容被屏蔽
6楼2014-06-16 12:35:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nwj0502

至尊木虫 (知名作家)
引用回帖:
6楼: Originally posted by gearboy78 at 2014-06-16 12:35:32
不调试直接看源码的话只能看语法是否错误，你编译通过了语法应该是没问题的，其它就是系统问题了，不去加载调试很难看出问题。

嗯，谢谢你，最近自己也一直在调试，找问题

[ 发自小木虫客户端 ]
赞一下 回复此楼
7楼2014-06-16 12:50:06
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gearboy78

禁虫 (小有名气)
本帖内容被屏蔽
8楼2014-06-16 13:23:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

longma135

金虫 (小有名气)
点火温度1473K(Tg=1473K)，时间60s；保温温度1043K(Tg=1043K)，时间90s；接下来抽风温度300K(Tg=300K)（压力出口，-10KPa）

我想问下上面这部分怎么做到的啊？
赞一下 回复此楼
9楼2014-06-18 11:01:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nwj0502

至尊木虫 (知名作家)
引用回帖:
9楼: Originally posted by longma135 at 2014-06-18 11:01:14
点火温度1473K(Tg=1473K)，时间60s；保温温度1043K(Tg=1043K)，时间90s；接下来抽风温度300K(Tg=300K)（压力出口，-10KPa）

我想问下上面这部分怎么做到的啊？

点火，保温，抽风的边界条件是通过UDF来完成的，具体帮助里有
赞一下 回复此楼
10楼2014-06-18 15:27:08
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nwj0502 的主题更新
101/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703调剂，一志愿天津大学319分 +16 haaaabcd 2026-04-05 16/800 2026-04-06 18:42 by lijunpoly
[考研] 277工科求调剂 +8 1915668 2026-04-04 9/450 2026-04-06 15:52 by lin-da
[考研] 工科 22408 267求推荐 +3 wanwan00 2026-04-05 4/200 2026-04-06 15:36 by lin-da
[考研] 290求调剂085701 +11 1314捧花 2026-04-02 11/550 2026-04-06 13:31 by zllcz
[考研] 302分求调剂 一志愿安徽大学085601 +9 zyx上岸！ 2026-04-04 9/450 2026-04-06 13:05 by 猪会飞
[考研] 0703化学调剂325分 +12 15771691647 2026-04-04 13/650 2026-04-06 12:00 by lijunpoly
[考研] 295求调剂 +8 FZAC123 2026-04-03 8/400 2026-04-05 17:46 by 蓝云思雨
[考研] 306分材料与化工求调剂 +7 黎吧啦啦你很有� 2026-04-03 7/350 2026-04-05 17:18 by Hdyxbekcb
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 11408,335分，本科211，求调剂，可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 085601，一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 一志愿南农090401，268，求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
[考研] 278求调剂 +6 Yy7400 2026-04-03 6/300 2026-04-04 09:53 by zhangdingwa
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[硕博家园] 考研调剂 +5 骆驼男人 2026-04-01 5/250 2026-04-01 14:28 by syjjj0321
[考研] 085601英二数二求调剂 总分325 +4 余航航 2026-03-31 4/200 2026-03-31 17:38 by 唐沐儿
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
信息提示
关闭
请填处理意见
关闭