用于预测气流床煤气化的EDC模型参数研究

doi:10.11949/0438-1157.20221131

摘要/Abstract

摘要：

气流床煤气化炉是多相湍流反应器，伴随燃烧、变换等均相反应过程，可由eddy dissipation concept（EDC）模型模拟湍流-化学反应相互作用，但需修正模型参数。以1500 t/d多喷嘴对置水煤浆气化炉为对象，分析了Realizable κ-ε条件下气化炉EDC模型的参数修正范围，考察了EDC模型参数对合成气温度和组成的影响。结果表明，当C_D1=0.213时，0.272≤C_D2≤2.268。在此范围修正EDC模型参数，气化炉出口温度、CO和H₂体积分数相对误差比修正前降低，最小可分别达到-2.54%、-3.93%和0.74%，模型修正范围可靠。默认参数的EDC模拟温度和组成结果始终处于边界参数模型范围内，且该范围越靠近气化炉出口越窄，各模型模拟结果趋于一致。

关键词: 气流床, 多喷嘴, 浆料, 气化, 数值模拟

Abstract:

The entrained flow coal gasifier is a multi-phase turbulent reactor accompanied by homogeneous reaction processes such as combustion and transformation, which can be simulated by eddy dissipation concept (EDC) model for turbulence-chemical reaction interaction, but the model parameters need to be corrected. The effects of EDC model parameters on syngas temperature and composition were investigated by analyzing the parameter modification range of EDC model for gasifier under Realizable κ-ε condition with a 1500 t/d opposed multi-burner (OMB) coal-water slurry (CWS) gasifier. Four groups of EDC models with different parameters and (finite-rate/eddy-dissipation) ED model was numerically investigated and verified with the outlet industrial data. The results showed that when C_D1=0.213,0.272≤C_D2≤2.268. After modification of the model constants in this range, the relative errors of temperature, CO and H₂ volume fraction at the outlet were lower than original value, and the minimum were -2.54%, -3.93%, and 0.74%, respectively, while those of the ED model were 7.95%, 5.52% and -2.74%, respectively. The temperature and composition results of the EDC simulation with default parameters are always within the range of the boundary parameter model, and the range becomes narrower closer to the gasifier outlet, and the simulation results of each model tend to be consistent.

Key words: entrained flow reactor, opposed multi-burner, slurry, gasification, numerical simulation

中图分类号:

TQ 546

仇鹏, 韩洋, 许建良, 王辅臣, 代正华. 用于预测气流床煤气化的EDC模型参数研究[J]. 化工学报, 2023, 74(1): 428-437.

Peng QIU, Yang HAN, Jianliang XU, Fuchen WANG, Zhenghua DAI. Study of EDC parameters for predicting entrained flow coal gasification[J]. CIESC Journal, 2023, 74(1): 428-437.

图/表 16

参考文献 33

1	王辅臣, 于广锁, 龚欣, 等. 射流携带床气化炉内宏观混合过程研究(Ⅰ): 冷态浓度分布[J]. 化工学报, 1997, 48(2): 193-199.
	Wang F C, Yu G S, Gong X, et al. Investigation of macro-mixing process for jet-entrained gasifier (‍Ⅰ): Cold model concentration distribution[J]. Journal of Chemical Industry and Engineering(China), 1997, 48(2): 193-199.
2	王辅臣, 龚欣, 吴韬, 等. 射流携带床气化炉内宏观混合过程研究(Ⅱ): 停留时间分布[J]. 化工学报, 1997, 48(2): 200-207.
	Wang F C, Gong X, Wu T, et al. Investigation of macro-mixing process for jet-entrained gasifier (Ⅱ): Cold model residence time distribution[J]. Journal of Chemical Industry and Engineering(China), 1997, 48(2): 200-207.
3	刘海峰, 王辅臣, 吴韬, 等. 撞击流反应器内微观混合过程的研究[J]. 华东理工大学学报(自然科学版), 1999(3): 228-232.
	Liu H F, Wang F C, Wu T, et al. Micromixing process in counterflowing gasifier[J]. Journal of East China University of Science and Technology (Natural Science Edition), 1999(3): 228-232.
4	Magnussen B F. On the structure of turbulence and a generalized eddy dissipation concept for chemical reaction in turbulent flow[C]//19th Aerospace Sciences Meeting. Virginia: AIAA, 1981: 42.
5	Magnussen B F. Modelling of NO _x and soot formation by the eddy dissipation concept[C]//1st Topic Oriented Technical Meeting. Amsterdam, Holland: 1989.
6	Rehm M, Seifert P, Meyer B. Theoretical and numerical investigation on the EDC-model for turbulence-chemistry interaction at gasification conditions[J]. Computers & Chemical Engineering, 2009, 33(2): 402-407.
7	He D, Yu Y S, Kuang Y C, et al. Analysis of EDC constants for predictions of methane MILD combustion[J]. Fuel, 2022, 324: 124542.
8	Ertesvåg I S. Analysis of some recently proposed modifications to the eddy dissipation concept (EDC)[J]. Combustion Science and Technology, 2020, 192(6): 1108-1136.
9	Ertesvåg I S. Scrutinizing proposed extensions to the eddy dissipation concept (EDC) at low turbulence Reynolds numbers and low Damköhler numbers[J]. Fuel, 2022, 309: 122032.
10	Li X Y, Dai Z H, Wang F C. Characteristic chemical time scale analysis of a partial oxidation flame in hot syngas coflow[J]. Energy & Fuels, 2017, 31(4): 4382-4390.
11	Xu Y T, Dai Z H, Li C, et al. Numerical simulation of natural gas non-catalytic partial oxidation reformer[J]. International Journal of Hydrogen Energy, 2014, 39(17): 9149-9157.
12	Caudal J, Fiorina B, Labégorre B, et al. Modeling interactions between chemistry and turbulence for simulations of partial oxidation processes[J]. Fuel Processing Technology, 2015, 134: 231-242.
13	Förster T, Voloshchuk Y, Richter A, et al. 3D numerical study of the performance of different burner concepts for the high-pressure non-catalytic natural gas reforming based on the Freiberg semi-industrial test facility HP POX[J]. Fuel, 2017, 203: 954-963.
14	Voloshchuk Y, Vascellari M, Hasse C, et al. Numerical study of natural gas reforming by non-catalytic partial oxidation based on the Virtuhcon Benchmark[J]. Chemical Engineering Journal, 2017, 327: 307-319.
15	Liu X, Zheng H T, Liu Q. Research of methane reforming and combustion characteristics in chemically recuperated gas turbine[J]. Industrial & Engineering Chemistry Research, 2014, 53(5): 1940-1946.
16	Li X Y, Dai Z H, Xu Y T, et al. Inverse diffusion flame of CH₄-O₂ in hot syngas coflow[J]. International Journal of Hydrogen Energy, 2015, 40(46): 16104-16114.
17	Richter A, Vascellari M, Nikrityuk P A, et al. Detailed analysis of reacting particles in an entrained-flow gasifier[J]. Fuel Processing Technology, 2016, 144: 95-108.
18	Chen T W, Zhang Q, Wang J F, et al. Simulation of industrial-scale gas quenching process for partial oxidation of nature gas to acetylene[J]. Chemical Engineering Journal, 2017, 329: 238-249.
19	Xiang Y L, Lin Q, Cai L, et al. Study of the effect mechanism of municipal solid waste gasification conditions on the production of H₂ and CO using modelling technique[J]. Journal of Environmental Management, 2019, 230: 301-310.
20	Cai L, He T Z, Xiang Y L, et al. Study on the reaction pathways of steam methane reforming for H₂ production[J]. Energy, 2020, 207: 118296.
21	Wang B, Qiu J Y, Guo Q H, et al. Numerical study on the effects of homogeneous reactions on the composition distributions of syngas in radiant syngas cooler[J]. Applied Thermal Engineering, 2022, 210: 118307.
22	Jones W P, Launder B E. The prediction of laminarization with a two-equation model of turbulence[J]. International Journal of Heat and Mass Transfer, 1972, 15(2): 301-314.
23	Shih T H, Zhu J, Lumley J L. A new Reynolds stress algebraic equation model[J]. Computer Methods in Applied Mechanics and Engineering, 1995, 125(1/2/3/4): 287-302.
24	Ertesvåg I S, Magnussen B F. The eddy dissipation turbulence energy cascade model[J]. Combustion Science and Technology, 2000, 159(1): 213-235.
25	Gran I R, Magnussen B F. A numerical study of a bluff-body stabilized diffusion flame (part 2): Influence of combustion modeling and finite-rate chemistry[J]. Combustion Science and Technology, 1996, 119(1/2/3/4/5/6): 191-217.
26	Sun Z H, Dai Z H, Zhou Z J, et al. Numerical simulation of industrial opposed multiburner coal-water slurry entrained flow gasifier[J]. Industrial & Engineering Chemistry Research, 2012, 51(6): 2560-2569.
27	代正华, 刘海峰, 于广锁, 等. 四喷嘴对置式撞击流的数值模拟[J]. 华东理工大学学报, 2004, 30(1): 65-68.
	Dai Z H, Liu H F, Yu G S, et al. Numeric simulation of the flow field in four opposed impinging jets[J]. Journal of East China University of Science and Technology, 2004, 30(1): 65-68.
28	李贤斌, 郭庆华, 张杰, 等. 多喷嘴对置式气化炉内气固两相流动与炉壁的颗粒捕捉特性[J]. 高校化学工程学报, 2014, 28(5): 957-964.
	Li X B, Guo Q H, Zhang J, et al. The gas-solid flow and characteristics of particles capture by the wall in an opposed multi-burner gasifier[J]. Journal of Chemical Engineering of Chinese Universities, 2014, 28(5): 957-964.
29	Li C, Dai Z H, Xu J L, et al. Numerical study of the particle residence time and flow characters in an opposed multi-burner gasifier[J]. Powder Technology, 2015, 286: 64-72.
30	Choi Y C, Li X Y, Park T J, et al. Numerical study on the coal gasification characteristics in an entrained flow coal gasifier[J]. Fuel, 2001, 80(15): 2193-2201.
31	Li C, Dai Z H, Li W F, et al. 3D numerical study of particle flow behavior in the impinging zone of an opposed multi-burner gasifier[J]. Powder Technology, 2012, 225: 118-123.
32	倪建军, 梁钦锋, 代正华, 等. 撞击流气化炉内气粒两相流动的数值模拟[J]. 化工学报, 2009, 60(4): 864-871.
	Ni J J, Liang Q F, Dai Z H, et al. Numerical simulation of gas-particle two-phase flow in impinging streams gasifier[J]. CIESC Journal, 2009, 60(4): 864-871.
33	于遵宏, 沈才大, 王辅臣, 等. 水煤浆气化炉气化过程的三区模型[J]. 燃料化学学报, 1993, 21(1): 90-95.
	Yu Z H, Shen C D, Wang F C, et al. Three-region model for gasification process in coal-slurry-gasifier[J]. Journal of Fuel Chemistry and Technology, 1993, 21(1): 90-95.

Case	原料	湍流方程	次级参数	文献
1~12	CH₄	Realizable κ-ε	C_τ =0.4083~8，C_γ =2.1377~13	[6]
13~15	CH₄	Realizable κ-ε	C_τ =0.4083，C_γ =2.14~4.4	[11]
16	CH₄	DNS/Standard κ-ε	C_τ =0.41，C_γ =4.6	[12]
17	CH₄	SST κ-ω	C_τ =0.4083，C_γ =10.7	[13]
18~23	CH₄	Realizable κ-ε	C_τ =0.4083，C_γ =4.3~21.45	[14]
—	CH₄	Standard κ-ε	not modified	[15]
—	CH₄	Standard κ-ε	not modified	[16]
—	coal	Realizable κ-ε	not modified	[17]
—	CH₄	Standard κ-ε	not modified	[10]
—	CH₄	Standard κ-ε	not modified	[18]
—	solid waste	Standard κ-ε	not modified	[19]
—	CH₄	Standard κ-ε	not modified	[20]
—	coal	Realizable κ-ε	not modified	[21]

Case	原料	湍流方程	次级参数	文献
1~12	CH₄	Realizable κ-ε	C_τ =0.4083~8，C_γ =2.1377~13	[6]
13~15	CH₄	Realizable κ-ε	C_τ =0.4083，C_γ =2.14~4.4	[11]
16	CH₄	DNS/Standard κ-ε	C_τ =0.41，C_γ =4.6	[12]
17	CH₄	SST κ-ω	C_τ =0.4083，C_γ =10.7	[13]
18~23	CH₄	Realizable κ-ε	C_τ =0.4083，C_γ =4.3~21.45	[14]
—	CH₄	Standard κ-ε	not modified	[15]
—	CH₄	Standard κ-ε	not modified	[16]
—	coal	Realizable κ-ε	not modified	[17]
—	CH₄	Standard κ-ε	not modified	[10]
—	CH₄	Standard κ-ε	not modified	[18]
—	solid waste	Standard κ-ε	not modified	[19]
—	CH₄	Standard κ-ε	not modified	[20]
—	coal	Realizable κ-ε	not modified	[21]

Case No.	C_τ	C_γ	C_D1	C_D2	C_R	α_L	Re_T(1)	Re_T(2)	L^*/η	τ^/τ_η*	Re^*	Re_T atε₂/ε=0.1	ε₂/ε atRe_T=50	ε₂/ε atRe_T=1	Re_T atγ_λ=0.75
1	0.4083	2.1377	0.134	0.50	11.2	27.84	20.9	0.2	1.4	0.4	2.49	225.5	0.2848	0.8276	66.0
2	0.4083	3.2066	0.060	0.50	25.2	140.97	105.7	0.2	2.1	0.4	5.60	1141.8	0.5559	0.9192	334.1
3	0.4083	5.56	0.020	0.50	75.7	1274.20	955.7	0.2	3.7	0.4	16.83	10321.0	0.8206	0.9724	3020.3
4	0.4083	8.45	0.009	0.50	174.9	6797.76	5098.3	0.2	5.6	0.4	38.87	55061.8	0.9178	0.9879	16113.2
5	0.4083	13	0.004	0.50	413.9	38081.33	28561.0	0.2	8.7	0.4	92.00	308458.8	0.9644	0.9949	90266.9
6	8.0000	13	0.071	192.00	21.1	38081.33	28561.0	64.0	169.8	8.0	1802.67	308458.8	0.9644	0.9949	90266.9
7	8.0000	4	0.750	192.00	2.0	341.33	256.0	64.0	52.3	8.0	170.67	2764.8	0.6836	0.9473	809.1
8	8.0000	5	0.480	192.00	3.1	833.33	625.0	64.0	65.3	8.0	266.67	6750.0	0.7832	0.9660	1975.3
9	8.0000	6	0.333	192.00	4.5	1728.00	1296.0	64.0	78.4	8.0	384.00	13996.8	0.8438	0.9762	4096.0
10	8.0000	7	0.245	192.00	6.1	3201.33	2401.0	64.0	91.4	8.0	522.67	25930.8	0.8826	0.9825	7588.3
11	8.0000	10	0.120	192.00	12.5	13333.33	10000.0	64.0	130.6	8.0	1066.67	108000.0	0.9406	0.9914	31604.9
12	8.0000	13	0.071	192.00	21.1	38081.33	28561.0	64.0	169.8	8.0	1802.67	308458.8	0.9644	0.9949	90266.9
13	0.4083	2.14	0.134	0.50	11.2	27.96	21.0	0.2	1.4	0.4	2.49	226.5	0.2855	0.8279	66.3
14	0.4083	3.2	0.060	0.50	25.1	139.81	104.9	0.2	2.1	0.4	5.57	1132.5	0.5546	0.9189	331.4
15	0.4083	4.4	0.032	0.50	47.4	499.75	374.8	0.2	2.9	0.4	10.54	4047.9	0.7298	0.9563	1184.6
16	0.4100	4.6	0.029	0.50	51.6	596.99	447.7	0.2	3.1	0.4	11.57	4835.7	0.7495	0.9599	1415.1
17	0.4083	10.7	0.005	0.50	280.4	17477.28	13108.0	0.2	7.1	0.4	62.33	141566.0	0.9479	0.9925	41427.6
18	0.4083	4.3	0.033	0.50	45.3	455.84	341.9	0.2	2.9	0.4	10.07	3692.3	0.7191	0.9542	1080.5
19	0.4083	6.4	0.015	0.50	100.3	2236.96	1677.7	0.2	4.3	0.4	22.30	18119.4	0.8613	0.9791	5302.4
20	0.4083	8.55	0.008	0.50	179.0	7125.30	5344.0	0.2	5.7	0.4	39.80	57714.9	0.9197	0.9882	16889.6
21	0.4083	12.95	0.004	0.50	410.7	37498.84	28124.1	0.2	8.6	0.4	91.30	303740.6	0.9641	0.9948	88886.1
22	0.4083	15	0.003	0.50	551.1	67500.00	50625.0	0.2	10.0	0.4	122.49	546750.0	0.9732	0.9962	160000.0
23	0.4083	21.45	0.001	0.50	1126.9	282259.08	211694.3	0.2	14.3	0.4	250.48	2286298.6	0.9868	0.9981	669058.6
24	0.8695	2.4745	0.213	2.268	7.0	49.99	37.5	0.8	3.5	0.9	7.10	404.9	0.3819	0.8682	118.5
25	0.3011	1.4562	0.213	0.272	7.0	6.00	4.5	0.1	0.7	0.3	0.85	48.6	0.0976	0.6666	14.2

Case No.	C_τ	C_γ	C_D1	C_D2	C_R	α_L	Re_T(1)	Re_T(2)	L^*/η	τ^/τ_η*	Re^*	Re_T atε₂/ε=0.1	ε₂/ε atRe_T=50	ε₂/ε atRe_T=1	Re_T atγ_λ=0.75
1	0.4083	2.1377	0.134	0.50	11.2	27.84	20.9	0.2	1.4	0.4	2.49	225.5	0.2848	0.8276	66.0
2	0.4083	3.2066	0.060	0.50	25.2	140.97	105.7	0.2	2.1	0.4	5.60	1141.8	0.5559	0.9192	334.1
3	0.4083	5.56	0.020	0.50	75.7	1274.20	955.7	0.2	3.7	0.4	16.83	10321.0	0.8206	0.9724	3020.3
4	0.4083	8.45	0.009	0.50	174.9	6797.76	5098.3	0.2	5.6	0.4	38.87	55061.8	0.9178	0.9879	16113.2
5	0.4083	13	0.004	0.50	413.9	38081.33	28561.0	0.2	8.7	0.4	92.00	308458.8	0.9644	0.9949	90266.9
6	8.0000	13	0.071	192.00	21.1	38081.33	28561.0	64.0	169.8	8.0	1802.67	308458.8	0.9644	0.9949	90266.9
7	8.0000	4	0.750	192.00	2.0	341.33	256.0	64.0	52.3	8.0	170.67	2764.8	0.6836	0.9473	809.1
8	8.0000	5	0.480	192.00	3.1	833.33	625.0	64.0	65.3	8.0	266.67	6750.0	0.7832	0.9660	1975.3
9	8.0000	6	0.333	192.00	4.5	1728.00	1296.0	64.0	78.4	8.0	384.00	13996.8	0.8438	0.9762	4096.0
10	8.0000	7	0.245	192.00	6.1	3201.33	2401.0	64.0	91.4	8.0	522.67	25930.8	0.8826	0.9825	7588.3
11	8.0000	10	0.120	192.00	12.5	13333.33	10000.0	64.0	130.6	8.0	1066.67	108000.0	0.9406	0.9914	31604.9
12	8.0000	13	0.071	192.00	21.1	38081.33	28561.0	64.0	169.8	8.0	1802.67	308458.8	0.9644	0.9949	90266.9
13	0.4083	2.14	0.134	0.50	11.2	27.96	21.0	0.2	1.4	0.4	2.49	226.5	0.2855	0.8279	66.3
14	0.4083	3.2	0.060	0.50	25.1	139.81	104.9	0.2	2.1	0.4	5.57	1132.5	0.5546	0.9189	331.4
15	0.4083	4.4	0.032	0.50	47.4	499.75	374.8	0.2	2.9	0.4	10.54	4047.9	0.7298	0.9563	1184.6
16	0.4100	4.6	0.029	0.50	51.6	596.99	447.7	0.2	3.1	0.4	11.57	4835.7	0.7495	0.9599	1415.1
17	0.4083	10.7	0.005	0.50	280.4	17477.28	13108.0	0.2	7.1	0.4	62.33	141566.0	0.9479	0.9925	41427.6
18	0.4083	4.3	0.033	0.50	45.3	455.84	341.9	0.2	2.9	0.4	10.07	3692.3	0.7191	0.9542	1080.5
19	0.4083	6.4	0.015	0.50	100.3	2236.96	1677.7	0.2	4.3	0.4	22.30	18119.4	0.8613	0.9791	5302.4
20	0.4083	8.55	0.008	0.50	179.0	7125.30	5344.0	0.2	5.7	0.4	39.80	57714.9	0.9197	0.9882	16889.6
21	0.4083	12.95	0.004	0.50	410.7	37498.84	28124.1	0.2	8.6	0.4	91.30	303740.6	0.9641	0.9948	88886.1
22	0.4083	15	0.003	0.50	551.1	67500.00	50625.0	0.2	10.0	0.4	122.49	546750.0	0.9732	0.9962	160000.0
23	0.4083	21.45	0.001	0.50	1126.9	282259.08	211694.3	0.2	14.3	0.4	250.48	2286298.6	0.9868	0.9981	669058.6
24	0.8695	2.4745	0.213	2.268	7.0	49.99	37.5	0.8	3.5	0.9	7.10	404.9	0.3819	0.8682	118.5
25	0.3011	1.4562	0.213	0.272	7.0	6.00	4.5	0.1	0.7	0.3	0.85	48.6	0.0976	0.6666	14.2

项目	出口温度/℃	相对误差/%	CO体积分数 (干基)/%	相对误差/%	H₂体积分数 (干基)/%	相对误差/%	碳转化率/%
工业数据	1220	—	46.77	—	35.10	—	98.8
ED	1317	7.95	49.35	5.52	34.14	-2.74	98.1
EDC-case 1	1138	-6.72	42.15	-9.88	37.71	7.44	98.1
EDC-case 2	1144	-6.23	42.55	-9.02	37.56	7.01	98.1
EDC-case 24	1149	-5.82	43.59	-6.80	35.74	1.82	98.1
EDC-case 25	1189	-2.54	44.93	-3.93	35.36	0.74	98.1