24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3287  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
回复此楼

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

模拟相关

» 猜你喜欢

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

已阅   关注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的回帖
查看全部 10 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 求调剂 有机化学考研356分 +11 Nadiums 2026-03-25 12/600 2026-03-31 18:24 by oooqiao
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
[考研] 生物学296求调剂 +8 汤圆包 2026-03-29 12/600 2026-03-31 17:05 by 18828373951
[考研] 271求调剂 +5 勒布朗@ 2026-03-31 7/350 2026-03-31 16:26 by hypershenger
[考研] 材料调剂 5+5 想要一壶桃花水 2026-03-25 11/550 2026-03-31 13:34 by Huaxue_Wang
[考研] 274求调剂 +6 xiao爱同学 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 335求调剂 +3 321* 2026-03-31 3/150 2026-03-31 09:00 by 功夫疯狂
[考研] 材料与化工调剂一志愿大连海事085600,349 +6 吃的不少 2026-03-30 6/300 2026-03-31 04:17 by fmesaito
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +9 @taotao 2026-03-29 9/450 2026-03-30 22:29 by 我是小康
[考研] 11408总分309,一志愿东南大学求调剂,不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 322求调剂 +10 宋明欣 2026-03-27 10/500 2026-03-30 18:47 by 544594351
[考研] 284求调剂 +14 junqihahaha 2026-03-26 15/750 2026-03-30 14:12 by 探123
[考研] 296求调剂 +10 彼岸t 2026-03-29 10/500 2026-03-30 10:50 by 探123
[硕博家园] 招收生物学/细胞生物学调剂 +4 IceGuo 2026-03-26 5/250 2026-03-29 01:25 by griffith2014
[考研] 304求调剂 +6 曼殊2266 2026-03-27 6/300 2026-03-28 14:10 by 唐沐儿
[考研] 085701环境工程,267求调剂 +16 minht 2026-03-26 16/800 2026-03-28 12:16 by zllcz
[考研] 一志愿郑大085600,310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 一志愿吉大071010,316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
信息提示
请填处理意见