2021-01-23热解系数校核记录 DEMPIC

newCoeff3:根据殷维杰的方法设置热解

其方法总结如下

A: 得到分子式的方法不变

按照干燥无灰基碳氢氧质量比

总分子量452

C19.32 H27.57 O12.03
$$C_{19.32}{H}_{27.57}{O}_{12.03}$$

B: 固定碳系数的方法也不变

干燥无灰基挥发分比例：79.67%
则固定碳比例为：20.33%

热解方程碳单质系数
452*0.2033/12=7.66

C:得到气体的系数（重要）

原方程：

chem_eq = “BiomassCore --> 6.53CO + 2.95CO2 + 2.79CH4 + 10.245H2 +” &
" 6.19*Carbon

根据Yin的方法

CO2取2.95， 根据O元素守恒确定CO，根据C元素守恒确定CH4，根据H元素守恒确定H2

实践如下：

C.1 根据O元素守恒确定CO

$C_{19.32}{H}_{27.57}{O}_{12.03}->x\ast CO+2.95\ast CO2+y\ast CH4+z\ast H2+7.66\ast Carbon$

除去CO2剩余O：12.03-2.95*2=6.13

即CO系数为6.13

C.2 根据C元素守恒确定CH4

$C_{19.32}{H}_{27.57}{O}_{12.03}->6.13\ast CO+2.95\ast CO2+y\ast CH4+z\ast H2+7.66\ast Carbon$

除去CO和CO2和Carbon剩余C: 19.32-6.13-2.95-7.66=2.58

即CH4系数为2.58

C.3 根据C元素守恒确定CH4

$C_{19.32}{H}_{27.57}{O}_{12.03}->6.13\ast CO+2.95\ast CO2+2.58\ast CH4+z\ast H2+7.66\ast Carbon$

除去CH4剩余H: 27.57-2.58*4=17.25

即H2系数为(27.57-2.58*4)/2=8.625

最终确定热解方程为

$C_{19.32}{H}_{27.57}{O}_{12.03}-->6.13CO+2.95CO2+2.58CH4+8.625H2+7.66Carbon$

chem_eq = "BiomassCore --> 6.13CO  + 2.95*CO2 + 2.58*CH4 + 8.625*H2 +" &
        " 7.66*Carbon

D：进行模拟校准

调整系数需要注意几点：

1. 调整四种气体的比例，但是必须保证原子守恒
2. carbon的系数不能变
3. 生物质分子式不能变

同时模拟要注意：

1. 采用切片法后处理
2. 去除BC中raw biomass的碳单质，也就是保证BiomassCore=0.9094
3. 先采用PIC 20s模拟

根据newCoeff1校核的结果，保证四种气体的摩尔分数，也就是系数之比，如下：

名称	系数之比（%）
H2	37.0
CO	9.6
CO2	14.2
CH4	39.3

该比值作为结果可以借鉴，但是不能作为判断标准，因为尚未使用切片法校核

D.1 第一次校准结果（基准校核）

PIC newCoeff3_element 耗时13.56h模拟46.766s

运行46s结果

名称	实验值	TOP_SUM	AVG_ALL	PV全域	PV切片(直接molar)	PV切片(mass转molar,脚本转)	PV切片(mass转molar,filter转)	AVG_ALL(origin处理)	TOP_SUM(origin处理)
H2	43.6	48.24	46.08	49.00	51.07	48.42	48.42	46.08	48.24
CH4	11.5	10.39	11.18	10.08	9.38	10.38	10.38	11.18	10.39
CO	33.2	21.14	24.84	25.19	21.15	21.06	21.06	24.84	21.14
CO2	11.7	20.22	17.89	15.73	18.41	20.14	20.14	17.89	20.22

结论：

1. PV切片法，先算出质量分数再转换成摩尔分数，脚本转和filter转完全一致（证明filter转molar可行）
2. PV切片法 直接molar 与 PV切片法 mass 转molar，两者差异较大，最大差距2.65
3. PV切片法 与 PV全域法 对比， 两者差异较大，最大差距4.04
4. PV全域与monitor ALL_AVG 对比， 两者差距较大，最大差距2.92
5. monitor TOP_SUM与 monitor ALL 对比，两者差距较大，最大差距3.7（CO）
6. monitor TOP_SUM与 切片法（mass转molar）对比，两者差距较小，最大差距0.08
7. 无论是TOP SUM还是AVG_ALL 用脚本积分和用origin积分，结果完全一致（证明自写的python脚本完全正确）

排除PV切片法直接molar

最后选用monitor TOP_SUM 脚本积分法 确定后处理

D.2 第一次校准更改系数

以TOP_SUM 脚本处理结果为依据

名称	实验值	TOP_SUM
H2	43.6	48.24
CH4	11.5	10.39
CO	33.2	21.14
CO2	11.7	20.22

根据模拟结果增减系数（尚未校核元素守恒）

进行元素守恒校核，保证原料分子式不变的情况下，原子守恒

校核结果为

最终得到系数

产物	系数
CO	9.82
CO2	1.105
CH4	2.9175
H2	7.95
Carbon	5.4775

最终确定热解方程为

$C_{19.32}{H}_{27.57}{O}_{12.03}-->9.82CO+1.105CO2+2.9175CH4+7.95H2+5.4775Carbon$

chem_eq = "BiomassCore --> 9.82*CO  + 1.105*CO2 + 2.9175*CH4 + 7.95*H2 +" &
        " 5.4775*Carbon

D.3 第二次校准结果

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-PGGRflx7-1611388948237)(E:\Pictures\markdown图片\21-1-20更改热解系数newCoeff3\校核2结果和更改系数)]

	实验值	TOP_SUM
H2	43.60	45.02
CH4	11.5	10.65
CO	33.20	29.14
CO2	11.70	15.19

元素守恒校核结果

最终确定热解方程为

$C_{19.32}{H}_{27.57}{O}_{12.03}-->10.45CO+0.79CO2+3.1CH4+7.585H2+4.98Carbon$

chem_eq = "BiomassCore --> 10.45*CO  + 0.79*CO2 + 3.1*CH4 + 7.585*H2 +" &
        " 4.98*Carbon

D.4 第三次校核结果

算例位置: /home/ubuntu1/disk_DATA/Works/MFiX_cases/save_prj_1121/DEMPIC/PRJ DEMPIC_rerun_90s/newCoeff3_element/PIC_base
处理的csv文件: GAS_YIELD_TOP_SUM.csv
积分时间区间(s): 0.0 - 30.000009000000002
列名： [’# “Time”’ ‘y_g(1)’ ‘y_g(3)’ ‘y_g(4)’ 'y_g(6) ']
43.66368072428311
11.239591253274488
30.790435934368798
14.30629208807359

H2
CH4
CO
CO2

最终确定热解方程为

$C_{19.32}{H}_{27.57}{O}_{12.03}-->10.99CO+0.522CO2+3.14CH4+7.505H2+4.67Carbon$

chem_eq = "BiomassCore --> 10.99*CO  + 0.522*CO2 + 3.14*CH4 + 7.505*H2 +" &
        " 4.67*Carbon

D.5 第四次校核

算例位置: /home/ubuntu1/disk_DATA/Works/MFiX_cases/save_prj_1121/DEMPIC/PRJ DEMPIC_rerun_90s/newCoeff3_element/PIC_base
处理的csv文件: GAS_YIELD_TOP_SUM.csv
积分时间区间(s): 0.0 - 30.000807000000002
列名： [’# “Time”’ ‘y_g(1)’ ‘y_g(3)’ ‘y_g(4)’ 'y_g(6) ']
43.517253624890785
11.156975540847842
31.428828516021262
13.896942318240118

D.6 第五次校核

算例位置: /home/ubuntu1/disk_DATA/Works/MFiX_cases/save_prj_1121/DEMPIC/PRJ DEMPIC_rerun_90s/newCoeff3_element/PIC_base
处理的csv文件: GAS_YIELD_TOP_SUM.csv
积分时间区间(s): 0.0 - 30.000698999999997
列名： [’# “Time”’ ‘y_g(1)’ ‘y_g(3)’ ‘y_g(4)’ 'y_g(6) ']
43.0254275302088
11.322001714000981
32.24497618010058
13.407594575689629

D.7 校核结束 热解系数最终结果

最终结果

热解方程

校核结束最终确定热解方程为

$C_{19.32}{H}_{27.57}{O}_{12.03}-->11.37CO+0.33CO2+3.213CH4+7.36H2+4.41Carbon$

Pyrolysis {
    chem_eq = "BiomassCore --> 11.37*CO + 0.33*CO2 + 3.213*CH4 + 7.36*H2 +" &
        " 4.41*Carbon"
    dh = 10000000.0
    fracdh(0) = 0.5
    fracdh(1) = 0.5
}

附录：热解校核excel文件和结果excel文件

DEMPIC最终全结果汇总

链接：https://pan.baidu.com/s/11pvtPIepkKNzOMGVLI72wg
提取码：nv0z
复制这段内容后打开百度网盘手机App，操作更方便哦–来自百度网盘超级会员V4的分享

DEMPIC最终热结系数校核方法

链接：https://pan.baidu.com/s/1cE1Nv_OdyfXyJRpOkJOK4w
提取码：x9i5
复制这段内容后打开百度网盘手机App，操作更方便哦–来自百度网盘超级会员V4的分享

附录：已发表论文DOI

Li, C., Eri, Q. Comparison between two Eulerian-Lagrangian methods: CFD-DEM and MPPIC on the biomass gasification in a fluidized bed. Biomass Conv. Bioref. (2021). https://doi.org/10.1007/s13399-021-01384-2

https://doi.org/10.1007/s13399-021-01384-2
或
https://link.springer.com/article/10.1007/s13399-021-01384-2

附录：已废弃的校核记录

2020-1-11

2 completed cases:

rerun1028 PIC

/home/ubuntu1/disk_DATA/Works/MFiX_cases/save_prj_1121/DEMPIC/PRJ DEMPIC_rerun_90s/rerun1028_90s_PIC/
serial
1day 2h 2m=90s

2DonePar_Gasi

/home/ubuntu1/disk_720/Lu2020Complex/Lu2020_2DPatel/2DOnePar_Gasi/
serial
3day 10h 56m=60s
测试是否都堆在x_p(31)
发现有些颗粒是x_p(31)=1，有些颗粒则x_p(31)=0.2左右，而char在0.7

2020-1-12

###改论文后处理：遇到问题，PIC传热之和2.7，DEM传热之和1.9.
PIC无Q_pp Q_wall Q_rad， 只有Q_conv和Q_rxns
干脆DEM也取消

###测试提高模型精度
按照DEM和PIC平均值 认为是模拟值
相对误差分母取模拟值
CO过低，为25.4，exp 33.2，+30.7%
CO2过高，为15.6，exp 11.7，-25%
CH4过低，为9.8，exp 11.5，+17.3%
H2过高，为49.3，exp 43.6，-11.5%

chem_eq = “BiomassCore --> 6.53CO + 2.95CO2 + 2.79CH4 + 10.245H2 +” &
" 6.19*Carbon

452.49421 BiomassCore分子量
保证分子量不变，剩下的给carbon

改为
chem_eq = “BiomassCore --> 8.53CO + 2.21CO2 + 3.27CH4 + 9.07H2 +” &
" 3.79*Carbon

算例：newCoeff1

改变点：

1. 改变热解系数
2. 取消DEM的辐射，壁面导热和颗粒导热

失败结果

名称	newCoeff1	实验	原PIC模拟
H2	44.3	43.6	49.0
CO	12.0	33.2	25.0
CO2	29.8	11.7	15.8
CH4	14.2	11.5	10.2

重新根据碳氢氧确定分子式

碳氢氧比例：
实验值给定干燥无灰基为：51.3 6.1 42.6
总分子量452

C 45251.3%/12=19.32
H 4526.1%=27.57
O 452*42.6%/16=12.03

C19.32 H27.57 O12.03

碳系数

干燥无灰基挥发分比例：79.67%
则固定碳比例为：20.33%

热解方程碳单质系数
452*0.2033/12=7.66

算例：newCoeff2

newCoeff1错误总结：调整一个系数的时候，其他气体系数应当按比例缩小，而不是全部缩小到char上

而保持char的系数不变，给定为7.66

剩余四种气体分子量为452-7.66*12=360.08

各个系数原始摩尔比例为：
CO:CO2:CH4:H2=6.53:2.95:2.79:10.245=29% :13.1% :12.4% :45.5%
转化为质量分数
48.4:34.4:11.8:5.4
分配分子量
174.28
123.72
42.55
19.53
得到新系数（保持摩尔比例不变）
CO 6.22
CO2 2.81
CH4 2.66
H2 9.77

先暂且不改变系数比例进行newCoeff2

原子不守恒
原料C19.32 H27.57 O12.03
产物含char C19.35 H30.18 O11.84

微调保证元素守恒

先调整CO保证O CO 6.41
再调整CH4保证C CH4 2.44
最后调整H2保证H H2 8.905

newCoeff2最终系数为

1. 方程系数
   chem_eq = “BiomassCore --> 6.41CO + 2.81CO2 + 2.44CH4 + 8.905H2 +” &
   " 7.66*Carbon

2. 微调BiomassCore分子量为452.31

进行PIC验证

发现之前是后处理错误了，反而是第一次改变coeff的结果成功了

但是仍然需要一些微调
H2让其下降1%(绝对值)
CO2让其下降1.5%（绝对值）
CO让其上升2%（绝对值）
CH4不变

同时，需要满足元素守恒
同时，需要去掉raw中的carbon

于是

目前要做的事情：

1. 保证热解方程元素守恒而且精度足够
2. 去掉BC中的carbon
3. 改变统计方法为monitor切片法