生物质一步热解（脱挥）方程的系数确定(2020-11-17 update)

最新推荐文章于 2022-12-01 12:01:39 发布

beidou111

最新推荐文章于 2022-12-01 12:01:39 发布

阅读量1k

点赞数 2

分类专栏： # 生物质

本文链接：https://blog.csdn.net/weixin_43940314/article/details/108609936

版权

生物质专栏收录该内容

8 篇文章 1 订阅

订阅专栏

文章目录

方法1：c3m方法
方法2：Basu 2011书中给出的方法
- Basu方法原理说明
方法3：根据热解结果倒推
附录1：未整理来源汇总
附录2：c3m中HPTR的matlab验证代码

方法1：c3m方法

1.1 c3m:HPTR实验法

1.1.1 HPTR法原理说明

2020-10-2
用matlab算了下HPTR的方法(见附录2)，比较简单。准确不准确就不知道了，因为里面只用到了灰分的质量分数和温度。
方法见下图
在这里插入图片描述

1.1.2 HPTR法验证

2020-10-7
今天确信这方法只看俩参数：
一个是生物质在混合物（煤炭与生物质）中的占比
一个是容器温度
其他一切参数都是没有用的。
这是随意定义的一个燃料的参数得到的结果
在这里插入图片描述
这是另一个燃料的结果，两者工业分析和元素分析完全不同
这个结果和我们用matlab代码计算的结果是非常相近的
证明了确实是这个算法计算的。

两种燃料定义的条件是一样的
在这里插入图片描述

所以得到了相同的结果。

但是如果对同一个燃料改变运行条件
例如temp改为900K
则得到了完全不同的结果
在这里插入图片描述

结论：HPTR计算热解系数，只和混合比例还有容器温度这俩参数有关，和燃料种类无关

所以这个模型非常不准确。实际上，燃料的种类已经被确定了，就是
在这里插入图片描述
完毕

1.2 C3M MGAS方法测试

MGAS方法只对元素分析敏感，与工业分析无关
以下是两种不同的燃料和结果
Fuel1
在这里插入图片描述

Fuel2

在这里插入图片描述

可见元素分析一样，但是工业分析不一样，得到的计算结果几乎一样。

1.3 C3M CPD方法

CPD也是只对元素分析敏感，对工业分析不敏感
并且，CPD不产生H2!
!
在这里插入图片描述

方法2：Basu 2011书中给出的方法

来自
[1] Basu P. 生物质气化与热解实用设计与理论
Biomass gasification and pyrolysis : practical design and theory. 2011.

第152页

Basu方法原理说明

在这里插入图片描述假定1 mol生物质被d mol 水蒸气 e mol O2， 3.76 mol N2气化
产物系数为n1~n7 共7个未知数。

a b c分别为生物质中H O N与C的原子数之比

根据原子守恒有下面四个式子
在这里插入图片描述

根据化学平衡常数有下面三个式子

在这里插入图片描述其中y是各组分摩尔分数，P是容器压力

R1 R2 R3分别代表三个碳转化反应

在这里插入图片描述而Ke是化学平衡常数，它可以表示为

在这里插入图片描述其中k_for是正向反应速率，k_back是逆向速率，C代表浓度，p q n 代表组分前系数

如图
在这里插入图片描述

所谓平衡常数就是正逆反应速率相等的时候，浓度乘积之比
在这里插入图片描述

已知各产物平衡时的浓度，就得到反应平衡常数，根据平衡常数，得到y之间关系，y之间关系就是摩尔分数关系，就是n之间关系。

因此7个方程7个未知数可以解。

方法3：根据热解结果倒推

来源：库晓珂老师的文章

[1] Ku X, Li T, Løvås T. Eulerian–Lagrangian Simulation of Biomass Gasification Behavior in a High-Temperature Entrained-Flow Reactor. Energy & Fuels 2014;28(8):5184-96.

在这里插入图片描述

Table 4就是挥发分脱除之后的四种气体的质量分数。只需要转换质量分数和摩尔分数，就得到了产物系数（即摩尔分数）。

转换摩尔分数和质量分数见：
https://blog.csdn.net/weixin_43940314/article/details/105819616

在这里插入图片描述

其中X是质量分数 Y是摩尔分数 MW是分子量

附录1：未整理来源汇总

fluent coal calculator
在这里插入图片描述

附录2：c3m中HPTR的matlab验证代码

matlab

%%
%%%%%%% 使用说明%%%%%%%
% 目标：根据X和T估算热解方程系数beta和alpha
% Yi是干燥基质量分数系数
% X是生物质在混合物(和煤)的质量分数
% T是固体最终温度(℃)
% 适用范围600-975℃
% 来源：c3m_theory 2.1.3.5, 文档见MFiX官网 c3m软件 https://mfix.netl.doe.gov/

%% 定义系数
clc;clear;
%给定X和T,需要手动改变
X=1;
T=800;

ash_dry=0.05;   %给定干燥基灰分


a0=[6.15E-01
5.65E-02
-1.21E-01
-2.05E-02
-2.49E-02
-2.46E-02
1.68E-01];

aX=[-8.01E-02
1.30E-01
-4.38E-01
-5.94E-02
2.49E-03
-2.32E-02
-4.89E-02];

aXX=[-3.99E-01
7.91E-02
-2.98E-03
-2.36E-02
3.84E-03
-4.51E-02
9.82E-02];

aT=[-1.35E-04
3.52E-05
2.73E-04
4.84E-05
4.29E-05
4.86E-05
-1.03E-04];

aXT=[-5.30E-04
-1.24E-04
7.85E-04
9.25E-05
-7.48E-06
3.75E-05
8.17E-05];

aXXT=[6.24E-04
-7.90E-05
5.78E-05
5.18E-05
-3.70E-06
8.52E-05
-1.12E-04];

%% 计算

Yi=a0+aX*X+aT*T+aXT*X*T+aXX*X*X+aXXT*X*X*T;

Ytar=1-sum(Yi);

Yi=[Yi;Ytar];

spcieName={
'Char'
'CO2'
'CO'
'CH4'
'H2'
'C2H4'
'H2O'
'Tar'};
spName = string(spcieName);

gamma_i=zeros(8,1); n_i=zeros(8,1);
gamma_char=(Yi(1)-ash_dry)/(1-ash_dry);

gamma_i(2:8)=Yi(2:end)/(1-ash_dry);

betad_i=gamma_i(2:8)/(1-gamma_char);

%假定挥发分分子量56 焦油144 
MW_VM=56;
MW_i=[12;44;28;16;2;28;18;144];
n_i(2:8)=MW_VM*betad_i./MW_i(2:end);

%% 输出结果:热解方程系数
%注：热解方程不包含焦炭，焦炭认为是剩余的未反应的固定碳。因此n_i第一个空出。
fid=fopen('result.txt','w');  
for ii=2:8
fprintf(fid,spName(ii));  
fprintf(fid,'\t%.4f \n',n_i(ii));
end
fclose(fid);