Numerical Study of the Pyrolysis of Ellipsoidal Low ... - ACS Publications

A discrete element method (DEM) model is used to describe the packing density distribution of ellipsoidal briquettes. The model is validated against t...
0 downloads 0 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Numerical Study of the Pyrolysis of Ellipsoidal Low-Rank Coal Briquettes Yuting Zhuo, Tianyu Wang, Changxing Li, and Yansong Shen Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b03224 • Publication Date (Web): 01 Jan 2018 Downloaded from http://pubs.acs.org on January 17, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

1

Numerical Study of the Pyrolysis of Ellipsoidal

2

Low-Rank Coal Briquettes

3

Yuting Zhuo, Tianyu Wang, Changxing Li, Yansong Shen*

4

School of Chemical Engineering, University of New South Wales, Sydney, NSW 2052,

5

Australia

6

KEYWORDS: Low-rank coal, pyrolysis, briquette, ellipsoid, CFD, DEM

7

ABSTRACT

8

Low-rank coal (LRC) upgrading is essential to convert LRC to a more thermal-efficient and

9

environmentally friendly fuel before utilization in coal-based industries. Briquetting and

10

pyrolysis is the dominant route where briquettes are usually in shape of ellipsoid through a

11

two-roll briquetting process. In this work, an integrated numerical model is developed to

12

predict the pyrolysis process of ellipsoidal LRC briquettes in a packed bed pyrolyzer. A

13

Computational Fluid Dynamics (CFD) model is developed to describe the flow and

14

thermochemical behaviours related to the pyrolysis of ellipsoidal LRC briquettes including

15

dewatering, pyrolysis, and other homogeneous and heterogeneous chemical reactions. A

16

Discrete Element Method (DEM) model is used to describe the packing density distribution

17

of ellipsoidal briquettes. The model is validated against the measurements in a pilot-scale test

18

rig in terms of temperature history and gas species yields. The typical in-furnace phenomena

19

are illustrated including flow field, temperature field, and products evolution. Then the

20

effects of key variables including briquettes properties and heating condition on pyrolysis

1 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

21

behaviour are investigated quantitatively. The effects of key parameters including briquettes

22

moisture, final pyrolysis temperature, and briquettes aspect ratio on pyrolysis are studied and

23

the optimal values are identified under the given conditions. For example, more H2 and less

24

CH4 are generated when the initial briquettes’ moisture content is increased in the range of 2%

25

to 15%, but H2 is then decreased and CH4 is increased in the range of 15% to 20%. The

26

maximum packing density is obtained when briquettes’ aspect ratio is 2.0. However, the

27

highest temperature increase rate is observed when the aspect ratio is 1.7. An appropriate

28

final pyrolysis temperature, 973K in this study, is suggested to balance pyrolysis rate and

29

energy consumption. This model provides a cost-effective tool for optimizing the design and

30

operation of pyrolyzers for ellipsoidal LRC briquettes.

2 ACS Paragon Plus Environment

Page 2 of 45

Page 3 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

31

1

32

The low-rank coals (LRCs), typically sub-bituminous and lignite, are geographically

33

dispersed, abundant, and accounts for almost half of the world's coal reserves

34

LRCs have some undesirable characteristics, such as high moisture content, low heat value,

35

and high mineral content etc., leading to limited use in industries. The upgraded LRCs have

36

been proven to improve plant efficiency, enhance safety, and reduce greenhouse gas

37

emissions3,4. On the other hand, as non-renewable premium coking coal reserve decreases

38

and the greater awareness of potential environmental risks of coal utilization, there is also a

39

need to convert LRCs to a thermal-efficient and environmentally friendly fuel for wider and

40

cleaner applications with higher economic value5. Thus, these necessitate LRCs upgrading

41

before effective applications in industries in power plant and metallurgical industries.

42

To date, the major LRC upgrading techniques include blending, briquetting, drying, cleaning

43

(removal of minerals) and chemical upgrading6,7. Out of these techniques, briquetting-

44

pyrolysis is a cost-effective and dominant route. LRC briquetting is a physical process where

45

excess moisture is removed by extruding coal fines with pressure and likely temperature,

46

producing briquettes in certain shapes, typically in ellipsoid shape using two-roll briquetting

47

machines. It has been widely investigated in terms of products yield and physical properties

48

worldwide8,9. On the other hand, pyrolysis is a thermochemical decomposition process where

49

the coal particles or briquettes are heated up without oxygen present, and produce tar, gases,

50

and semi-coke. Compared to raw coal combustion, pyrolysis is able to convert LRCs to more

51

diverse, specifically designed and profitable fuels. Moreover, in conventional industrial fixed

52

bed pyrolyzer, lump coal is usually the major material, allowing gas heat carrier and pyrolysis

53

gases pass through the bed smoothly and reduce the risk of explosion when semi-coke is

54

discharged after pyrolysis. Hence, fine coal cannot be directly utilized in such pyrolysis

55

process in the past, limiting LRC’s application. Combining briquetting and pyrolysis

INTRODUCTION

3 ACS Paragon Plus Environment

1,2

. However,

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

56

techniques, more fine coals can be used efficiently and the friability and emission problems

57

can be mitigated at a competitive cost against the traditional fuels10. Therefore, it is important

58

to understand and optimize the pyrolysis of LRC briquettes, typically in the shape of ellipsoid,

59

including the packing of ellipsoidal briquettes and their pyrolysis behaviour in the packed bed

60

pyrolyzer.

61

The pyrolysis of LRCs has been studied by different methods, experimentally and

62

numerically11. Extensive experimental studies have been conducted to investigate coal

63

pyrolysis process. For example, Li et al.12 reviewed the recent experimental studies of

64

pyrolysis and gasification behaviour of Victoria brown coal13,14. Due to the harsh thermal and

65

chemical conditions inside the pyrolysis reactor, it is difficult to measure the detail internal

66

information in a continuous operating pyrolyzer, such as temperature and gas evolution, but

67

such information plays a significant role in optimizing design and operation.

68

As the development of computer technology, numerical modelling has emerged as an

69

effective tool for understanding and optimizing complex reaction systems including the

70

pyrolysis of LRCs15–18. Some 2D CFD models were developed for describing the

71

carbonization process in coke ovens19–21. In these models, fine coal bed was treated as porous

72

medium and the flow pattern and evolution of coke oven gas was investigated. Moreover,

73

some multi-fluid models are developed to study the influence of operation conditions on coal

74

pyrolysis and gasification processes22. For example, Zhang et al.23,24 developed a CFD model

75

for municipal solid waste gasification process in an industrial scale fixed bed, and discussed

76

the effects of operation conditions on energy consumption and syngas composition. In

77

addition, some CFD-DEM models

78

pyrolysis and gasification process at particle sale, where particles were in the shape of sphere.

79

On the other hand, in the packing processes, packing density is widely used as a parameter in

80

particle packing characterization. The effect of particle shape on packing density has been

25–28

were developed for investigating coal or biomass

4 ACS Paragon Plus Environment

Page 4 of 45

Page 5 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

81

studied for years, for example, Zhou et al.29,30 and Li et al.31 Recently, the investigation

82

regarding the packing of non-spherical particles by DEM approach is increasingly gaining

83

popularity32,33. They showed that ellipsoid packing is quite different from the packing of

84

spheres. To date, few studies have been found in the literature on the modelling of the

85

pyrolysis or combustion of ellipsoidal briquettes34. In the present work, an integrated

86

numerical model is developed for describing the pyrolysis behaviour of ellipsoidal LRC

87

briquettes in a packed bed reactor. In this integrated model, a CFD model is developed for

88

describing the flow and thermochemical behaviour related to pyrolysis including dewatering,

89

pyrolysis, and other homogeneous and heterogeneous reactions, and a DEM model is used for

90

describing the packing density distribution of ellipsoidal briquettes. The model is validated

91

against the measurements in a pilot-scale test rig in terms of temperature history and gas

92

yields. The typical in-furnace phenomena are illustrated including spatial distribution of flow

93

field and temperature profile, and products evolutions. The effects of key variables on

94

pyrolysis are then investigated quantitatively including briquettes properties and external

95

heating conditions.

96

2

97

2.1 Experiment setup

98

A pilot-scale fixed bed pyrolyzer is built at Sinosteel, China and used to investigate pyrolysis

99

performance of LRC briquettes. The schematic of pilot-scale pyrolysis test system is shown

100

in Figure 1. The geometry information is shown in Figure 2. Briquettes are introduced into

101

the reactor at room temperature. The remaining air is purged by N2 with a rate of 2 L/min for

102

10 minutes. The reactor is heated up under the control system using a predefined scheme.

103

Four thermocouples are placed inside the wall for monitoring the heating wall’s practical

104

temperature and feedback to the temperature control system to adjust the temperature. One

105

more thermocouple is placed at the centre axis of the bed to measure the bed centre

EXPERIMENT DETAILS

5 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 45

106

temperature during the pyrolysis process. The gas products including vapour and tar releasing

107

from the reactor will outflow from the upper exit, and then are separated by the tar trapper

108

and dryer for being collected and weighted. The composition of the remained gases is

109

collected and analyzed by gas chromatography. 3.5kg (dry base) briquettes are used in each

110

test. The geometry and operational conditions will be used as boundary condition for the

111

numerical model, and the measured data will be used for model validation.

112

2.2

113

In this study, briquettes made from Shenmu LRCs through the two-roll briquetting process

114

are used in the experiment. The average particle size of an ellipsoidal briquette is shown in

115

Figure 3, where briquette’s semimajor axis length, a, is 25mm and semiminor axis length, b,

116

is 14.7mm. The average weight of each briquette is 12g with bulk density being 1280 kg/m3.

117

The proximate and ultimate analyses of the briquettes are shown in Table 1.

118

3

119

3.1

120

In the present work, a CFD model is developed for describing the pyrolysis process in the

121

fixed bed pyrolyzer. In this model, Euler-Euler multiphase approach is employed in the CFD

122

model. The gas flow is treated as a laminar flow of ideal gas

123

composition analyses in the experiment, 15 kinds of gas species are defined in this model for

124

the involved homogeneous and heterogeneous chemical reactions. The solid phase is treated

125

as a continuous granular flow. The flow of granular materials in the fixed bed pyrolyzer is

126

treated as plastic flow35–37. The radiation heat transfer is described by the discrete ordinate

127

model. On the other hand, a DEM model is used for describing the packing of the ellipsoidal

128

briquettes. The model is detailed elsewhere 31. The results of the DEM model will be used as

Properties of the briquettes

MODEL DETAILS Gas-solid flow and heat transfer

6 ACS Paragon Plus Environment

21

. According to the gas

Page 7 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

129

the initial condition of the CFD model. The governing equations of continuity, momentum,

130

and energy for solid and gas phases are listed in Table 2.

131

The Gidaspow drag model38 is used to calculate the drag force. This paper studies ellipsoidal

132

LRC briquettes rather than spherical particles. It is necessary to introduce the effect of

133

particle shape on drag force and heat transfer. Nanda et al.

134

describing the non-spherical particles in terms of drag force and heat transfer by modifying

135

drag force coefficient and Nusselt number, respectively. That is, the drag force coefficient Cd

136

for the ellipsoid particles proposed by Nanda et al.

137

interphase momentum exchange is as follows: r r  α sα g ρ g v s − v g 3  K sg = C D α g −2.65  ds 4 r r  ρ gα s v s − v g  α s (1 − α g ) µ g + 1.75  K sg = 150 ds α g ds2  CD =

39

39

proposed a modified model for

is used. The coefficient K sg of the

If α g ≥ 0.8

(7) If α g < 0.8

24 e 0 .4 9 (1.0 5 + 0.1 52 R e 0 .6 87 e 0 .6 7 1 ) Re

(8)

138

The interphase heat exchange intensity between solid phase and gas phase is driven by the

139

temperature difference, and calculated as follows: Q sg = − Q gs = k sg (T s − T g ) k sg =

(9)

6 k gα sα gNu s dS 2

(10)

Nus = 2e0.3 + Pr 0.4 (0.4 Re0.5 e0.83 + 0.06 Re2/3 e0.1 )

(11)

140

Here, the modified Nusselt number modified by Nanda et al.39 is introduced for considering

141

the influence of ellipsoid particle’s aspect ratio on interphase heat transfer.

7 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 45

142

3.2

143

Moisture retained in briquettes could be divided into three categories, free water, bound water

144

and water vapour, which will be driven out when particles are heated up

145

process involves free water evaporation, bound water desorption and evaporation as well as

146

the chemical bound water separation24. In the present work, a simplified drying model is used

147

to predict the moisture evolution in pyrolysis process24,40, as below.

Drying

25

. This complex

r = A v k m ( C M oi − C H 2 O )

when Ts < 373.15K

(12)

r = Q a / H evap

when Ts ≥ 373.15K

(13)

Q a = A s ( h s ' ( T s − T g ) + Q ra d

(14)

148

where Av is the specific area of particles, km is the interphase mass transfer coefficient

149

between solid and gas phase. Qa is a summation of heat absorbed by the solid phase.

150

3.3

151

The briquettes are assumed to be consist of four components, including water, volatile matter,

152

fixed carbon, and ash. The homogeneous reaction rate is calculated using the finite rate model,

153

and the unreacted shrinking core model

154

calculation. Reactions considered in this CFD model and their reaction kinetics are given in

155

Table 3. For example, many pyrolysis or devolatilization models have been proposed, such as

156

one step kinetic rate model42, two competing model43 and complex models whereby DAEM44,

157

FLASHCHAIN45, CPD46 and FG-DVC47. In the present model, the two steps model is used

158

to balance accuracy and computational load, where the coal will thermal decompose firstly to

159

primary tar, gas, and char at first and then primary tar will subsequently crack and convert to

160

the secondary tar and other hydrocarbon gases. T 1 is defined as the activation temperature at

161

which the first step of devolatilization is activated. It is estimated using the correlation

162

modified by Gregory and Littlejohn48. The CFD model is developed based on the framework

Reaction models

41

is applied for heterogeneous reactions rate

8 ACS Paragon Plus Environment

Page 9 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

163

of a commercial software package ANSYS-Fluent v17.2. The reaction rate, modified drag

164

force model, and modified heat transfer model introduced previously are written as UDF in

165

Fluent.

166 167

4

168

4.1

169

The properties of the coals used in the base case are the same as used in the experiment

170

(Table 1). The thermal properties of briquettes and operating conditions are shown in Table 4.

171

The simulation geometry is set as same as the experiment (Figure 2). It should be noted that

172

for evaluating the influence of final pyrolysis temperature on pyrolysis three heating schemes

173

for heating wall are defined. These heating schemes are used in the experimental studies to

174

make sure that the coal bed temperature can increase at a mild rate and therefore avoid crack

175

when the temperature increase too fast.

176

4.2

177

In order to control the briquettes bed height in the experiment, dense packing is conducted for

178

keeping the bed height as 0.355m in each test to guarantee the entire bed stay at the flat-

179

temperature zone of the furnace. The DEM model

180

process for obtaining the briquettes’ packing density distribution. The DEM code used in this

181

study has been validated in several cases29-31, where the simulation results were found

182

comparable with the experimental data. Based on the DEM simulation, the bed is divided into

183

8 sections and each section’s average packing density is calculated, as shown in Figure 4.

184

Based on this DEM simulation result, the correlation between bed height and packing density

185

is obtained as Equ 15. It will be used as the initial condition of the CFD-based pyrolysis

186

model.

SIMULATION CONDITIONS Materials properties, and boundary conditions for CFD model

Packing density

29–31

is applied to simulate the packing

9 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Pc = 0.543 + 0.874 y − 7.028 y 2 + 24.866 y3 − 32.408 y 4

Page 10 of 45

(15)

187

5

188

5.1 Model validation

189

The model is validated against the measurements in the experiments in term of temperature

190

history and gas yields. Figure 5 compares the predicted and measured temperature history at

191

the monitor point. The monitored point is set to the same position with the thermal probe

192

located in the experimental reactor (see Figure 2). Figure 6 illustrates the comparison of

193

generated gases composition between the simulation results and experimental data. It is

194

indicated in Figure 5 and Figure 6 that both the predicted temperature history and the gas

195

components show good agreements with the experimental data. This confirms the validity of

196

this model. This model slightly underestimates the mole fraction of CH4 and overestimates

197

the H2, this can be explained by the simplified devolatilization reaction assumption.

198

Generally, this deviation of predicted results is acceptable for understanding the

199

characteristics pyrolysis process.

200

5.2

201

In this section, the key in-furnace phenomena related to the pyrolysis of ellipsoidal LRC

202

briquettes are described including flow field, temperature field and gas distributions.

203

Figure 7 shows the temperature profiles of the briquettes and gas phases at the selected time

204

in the base case. Note that the wall temperature follows scheme 2 (see Table 4). It can be

205

observed that with the wall temperature increases, the briquettes bed temperature rises up

206

gradually from near wall section to the centre. With the reaction proceeds, the briquettes

207

temperature at upper part rises faster than that of the lower part, due to the intensive heat

208

transfer by convection between gas and solid phase. The lower temperature is found near the

209

bottom due to the adiabatic setting of the wall, consistent with the practical operation. At

RESULTS AND DISCUSSION

Typical phenomena of the pyrolysis of ellipsoidal LRC briquettes

10 ACS Paragon Plus Environment

Page 11 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

210

10000s, almost all briquettes are heated up to predefined final pyrolysis temperature (973K)

211

except for the bottom bed. It is indicated that the temperature history at the monitor point (see

212

Figure 5) has a plateau between around 1500s and 2100s, due to the endothermal dewatering,

213

then the temperature increase at a higher rate until the devolatilization reaction takes place.

214

The second plateau emerges after about 11000s where the briquettes bed’s temperature levels

215

off at 973K that reach the target heating temperature of the wall as the predefined scheme

216

(Table 5).

217

Coal pyrolysis is a complex process including gas-solid flow, heat transfer and a series of

218

chemical reactions, such as drying, devolatilization, water-gas shifting, etc. Devolatilization

219

is the key step in the entire pyrolysis process. Figure 8 upper part shows the remaining

220

volatile matter content distribution in the bed at the selected time. The initial VM content is

221

set to 32.7%. It is indicated that at the beginning stage, VM content decreases slowly. This is

222

because, the temperature of the centre part of briquettes bed is lower than that close to the

223

wall, as a result, the VM content will not change much until it reaches the activating

224

temperature. It is also observed that the remained VM content depends on the distance to the

225

wall. With increasing temperature of bed (Figure 7), the devolatilization reaction is favoured

226

by higher temperature. The heat exchange in the near-wall area and the roof of bed is more

227

intensive than other parts, resulting in the faster decrease of VM. Typically, the remaining

228

amount of volatile matter is used to characterize the extent of pyrolysis. At around 11000s all

229

the volatile matter is driven out. This means the pyrolysis reaction nearly finishes. Figure 8

230

middle and bottom part illustrate the evolutions of primary tar and secondary tar, respectively.

231

It can be observed that the primary tar resulting from the first step reaction in the

232

devolatilization model shows a lower mass fraction along the centreline initially and near the

233

bottom of reactor finally. This is because a gas flow stream exists along the centreline (Figure

234

9) driven by the temperature difference in briquettes bed. And the gas products in the second 11 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

235

step of devolatilization reaction tend to go through briquette bed from high temperature area

236

to low temperature area. Therefore, the primary tar mass fraction shows a lower value, while

237

the secondary tar shows a higher mass fraction in these areas.

238

Moreover, it is found in Figure 8 that as reaction proceeds, tar cracking becomes the

239

dominant reaction in the gas phase, and consequently, secondary tar concentration gradually

240

increases in the high-temperature area. It can be found that at the beginning stage of pyrolysis,

241

secondary tar emerges at the roof of reactor where the temperature is much higher than the

242

bed. As heating proceeds, more primary tar is converted to secondary tar in the briquettes bed.

243

This is induced by the increasingly higher gas temperature as well as the primary tar

244

concentration, where the latter accelerates the tar cracking rate and subsequently convert the

245

primary tar into other light hydrocarbon gases.

246

Typically, the in-furnace information of gas and moisture released from briquettes bed in an

247

industrial plant are difficult to measure. In order to produce results that can be more readily

248

and useful in practice, the mass flow rate of generated pyrolysis gases is monitored at the

249

outlet, as shown in Figure 10. Two peaks can be observed during the whole process. The first

250

peak is found ranging from 0 to 2500s, resulting from the heat transfer from the wall to the

251

briquettes bed and the intensive moisture evaporation. Before the secondary peak, the

252

pyrolysis gas release as a mild rate, because the bed average temperature has not reached the

253

point where the devolatilization rate becomes faster. The secondary peak is found around

254

8000s where the mass flow approach to the highest level and then decrease dramatically as

255

the remaining VM has been completely released. To observe the pyrolysis gas species

256

distributions at the middle stage of the process, profiles of some important gas species at

257

7500s (where mass flow rate reaches the peak) are shown in Figure 11. It is indicated that the

258

vapour is nearly disappeared at this moment, whereas CH4 and H2 show a relatively higher

259

concentration within the bed and CO2 shows a higher concentration near the roof. The latter 12 ACS Paragon Plus Environment

Page 12 of 45

Page 13 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

260

can be explained by the secondary tar cracking reaction takes place at that part and produce

261

more CO2.

262

5.3 Effects of key variables on ellipsoidal LRC pyrolysis

263

Pyrolysis performance may be directly affected by many conditions including charge

264

properties and operating conditions. In this section, the effects of briquettes’ moisture content,

265

ellipsoid’s aspect ratio, and final heating temperature on pyrolysis are investigated

266

quantitatively at the bed centre. Each operating condition is varied while keeping other

267

parameters constant. The varying parameters are listed in Table 5. Case 1 is the base case

268

where boundary conditions and charge properties are set as same as the experimental

269

conditions.

270

Effects of briquettes’ moisture content

271

Figure 12 shows the effect of initial moisture content of the briquettes on bed temperature

272

evolution (a), gas products composition (b), VM content (c) and mass flow rate at the outlet

273

(d), respectively, by means of comparing Cases 1-5. It can be found in Figure 12(a) that the

274

temperature increase rate gradually decreases in preheating stage and the dewatering time is

275

prolonged when briquettes bed has a higher initial briquettes moisture content. In particular,

276

the temperature difference is significant between Case 1 and Case 5 at the time when

277

dewatering is finished in Case 5, because more energy is needed for evaporating the

278

increased moisture in the briquettes. Even though in the later stage, the temperature

279

increasing rate in all cases becomes faster than preheating stage, yet Case 5 shows faster rate

280

than other cases, since less amount of VM remains and briquettes’ specific heat and thermal

281

conductivity are varied with the composition. Therefore, at the later stage of pyrolysis

282

process, the temperature difference between different cases is not as large as the initial stage.

283

Figure 12(b) indicates that CH4 mole fraction shows a slight decrease as the moisture increase

284

in Cases 1, 2, and 3, while H2 mole fraction shows an opposite trend due to more vapour 13 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

285

existing in the gas phase and favouring the reactions. However, Case 4 and Case 5 show

286

lower mole fraction of H2 and higher CH4 mole fraction compared to previous cases. That is,

287

more H2 and less CH4 are generated when the initial briquettes’ moisture content is increased

288

in the range of 2% to 15%, but then H2 decreased and CH4 increased in the range of 15% to

289

20%. This is because, the relative low bed temperature is obtained when briquette moisture

290

content is high, and thus the water-gas shift reaction rate is decreased.

291

Figure 12(c) compares the evolutions of volatile matter content remaining in the briquettes

292

for the briquettes beds with different initial moisture contents. In each case, VM mass

293

fraction shows a slight decrease in the beginning, due to the vapour condensation at the centre

294

bed. Then VM mass fraction increases dramatically as the moisture begins to release from the

295

briquettes bed until the dewatering process approach to complete. A stable stage occurs due

296

to the time gap between bed reaches evaporation temperature and devolatilization reaction

297

activated temperature. Therefore, before the start of devolatilization VM mass fraction keeps

298

constant for a while. On the other hand, these cases are compared. As higher moisture content

299

extends the dewatering time and therefore the beginning of devolatilization at the centre bed

300

is deferred when the briquettes with higher moisture is used. Moreover, the temperature is

301

increased after drying in all cases, where such increase in Case 1 is earlier than Case 5 by

302

over 2000s in the early stage, whereas the two cases reach the final temperature at the similar

303

time, within 1000s. Therefore, the temperature increases in Case 5, with a higher moisture

304

content, is faster and catches up in the later stage (see Figure 12(a)). This means that the

305

temperature is the critical factor to determine pyrolysis behaviour. Figure 12(d) indicates that

306

as the increase of moisture in the briquettes, at the beginning stage the outgoing mass flow

307

dramatically rises up due to the more vapour generated during the dewatering process. On the

308

other hand, the second peak emerges is deferred compared to other cases. This is due to the

309

relatively lower VM amount as well as the lower bed temperature. Even though higher 14 ACS Paragon Plus Environment

Page 14 of 45

Page 15 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

310

vapour generated in Case 5, the time difference between vapour and pyrolysis gases

311

generated also results in gas-water shift reaction doesn’t become the dominate homogeneous

312

reaction. Hence, the produced H2 does not show a significant increase when briquettes

313

moisture increases at the same time (see Figure 12(b)).

314

Effects of ellipsoidal briquette’s aspect ratio

315

In order to investigate the influence of ellipsoid shape on pyrolysis performance, in this

316

section, apart from the Case 1 (aspect ratio 1.7), four more aspect ratios are investigated and

317

compared. Four more packing density distribution results are obtained through the DEM

318

simulations based on these aspect ratios, as shown in Figure 13. Overall, because of the wall

319

effects, at the roof and bottom part, the packing density shows a little bit smaller that of other

320

parts. On the other hand, in comparison, it is indicated that the packing density is varied with

321

the aspect ratio. Specifically, the case with aspect ratio being 1.0 has the lowest packing

322

density distribution, while the highest value is obtained when the aspect ratio is 2.0. That is,

323

for the equivalent volume ellipsoidal briquettes, the packing density of charge is influenced

324

by the particle shape. The fitting curves of the packing density as a function of bed height are

325

obtained and then applied to the CFD model as initial conditions. In the CFD model, the

326

practical diameter of the briquettes is used to calculate the equivalent volume diameter. When

327

changing the aspect ratio, the equivalent diameter is fixed.

328

Figure 14 shows the effect of briquettes aspect ratio on bed temperature evolution, gas

329

products composition, the remaining VM content evolution and mass flow rate at the outlet,

330

respectively, by means of comparing between Cases 1 and 8-11. It is indicated in Figure 14(a)

331

that when other operation conditions are fixed, Case 1 with aspect ratio being 1.7 gives the

332

fastest heating rate after dewatering process than other cases, followed by Case 9, Case 11,

333

Case 8, Case 10. It can be found that the temperature difference between Case 1 and Case 10

334

is about 20K at a given time. It is indicated that the heating rate is not only determined by the 15 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

335

packing density but also determined by the heat exchange intensity that varies with the aspect

336

ratio. That is, the appropriate aspect ratio is beneficial to heat transfer as well as packing

337

process. In addition, the VM mass fraction evolution shows the same trend as temperature

338

varies (Figure 14(c)). However, the composition of generated gases does not show much

339

difference between these cases as aspect ratio changes. Moreover, Figure 14(d) indicates that

340

Case 1 has the smallest peak value of mass flow rate. This phenomenon indicates the

341

variation of aspect ratio has impacts on internal heat transfer and subsequent mass transfer.

342

Such difference will be more prominent if the gas heating is applied to the pyrolysis process.

343

They will be reported elsewhere.

344

Effects of different final pyrolysis temperature

345

In order to investigate impacts of final pyrolysis temperature on pyrolysis performance, three

346

cases with different final temperature are examined in Figure 15, namely Case 1 (973K),

347

Case 6 (873K), Case 7 (1073K), in terms of bed temperature history (a), gas products

348

composition (b), the remaining VM content evolution (c), mass flow rate at the outlet (d). It is

349

indicated that Case 6 with a lower final pyrolysis temperature shows a much slower VM

350

decrease compared to other cases. The VM decrease rate has a similar trend, which can be

351

found in Figure 15(c). In the comparison of gas products’ outflow rate, it can be found that

352

the Case 6 with 873K shows a significant difference with other cases, that is, the peak value

353

of mass flow rate of Case 7 is 2 times larger than Case 6. This means that higher final

354

pyrolysis temperature will enhance the devolatilization rate significantly. These results

355

indicate that the lower final pyrolysis temperature dramatically extends the time needed for

356

devolatilization reaction. It is important to note that the determination of final pyrolysis

357

temperature should consider the fact that the appropriate temperature can satisfy the demand

358

of pyrolysis completion with less time and energy consumption.

16 ACS Paragon Plus Environment

Page 16 of 45

Page 17 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

359

6

360

An integrated DEM-CFD numerical model is developed to investigate the packing and

361

pyrolysis process of ellipsoidal LRC briquettes. The model is validated against the

362

experimental measurements in a pilot scale fixed bed reactor in terms of briquettes

363

temperature history and gas products distribution. The typical in-furnace phenomena are

364

described, including packing density distribution, flow field, temperature field and gas

365

species distributions. Then the effects of some key variables on pyrolysis behaviour are

366

investigated in terms of LRC briquettes’ initial moisture content, heating condition, and

367

ellipsoid’ aspect ratio. Under the current conditions, some optimal values are found for this

368

pyrolyzer.

369

The main conclusions from the present study are summarized as follows:

370

CONCLUSIONS



The packing density distribution of the ellipsoidal briquettes varies with the briquettes

371

aspect ratio and then affects the heat and mass transfers in the pyrolysis process. As

372

briquette’s initial moisture content is higher, the bed temperature increase and volatile

373

release are slowed down and dewatering process is prolonged. At the latter stage of

374

pyrolysis, cases with higher initial moisture content give faster bed temperature

375

increase rates. But excessive moisture in briquettes will lead to the relative low bed

376

temperature and cause the water-gas shift reaction rate decrease, generating less H2

377

and more CH4.

378



Maximum packing density is obtained when briquettes’ aspect ratio is 2.0. However,

379

the highest temperature increase rate is observed when the aspect ratio is 1.70. The

380

heat transfer between bed and gas is significantly influenced by packing structure and

381

briquette aspect ratio.

17 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

382



Excessive low temperature (873K) leads to a much slower devolatilization reaction

383

rate so as to dramatically extend the reaction time. But the higher temperature (1073K)

384

does not favour the devolatilization too much. Therefore, an appropriate final

385

temperature, 973K in this study, is identified for lowering energy consumption in

386

pyrolysis process.

387

This model provides a cost-effective tool for understanding and optimizing the pyrolysis of

388

LRC ellipsoidal briquettes in a packed bed pyrolyzer.

389

NOMENCLATURE

Av

specific surface area (m-1)

As

Surface area of the briquettes (m2)

C

mass concentration (kg m-3)

CP

specific heat of briquette (J kg-1K-1)

CF

specific heat of fixed carbon (J kg-1K-1)

CV

specific heat of volatile matter (J kg-1K-1)

CA

specific heat of ash (J kg-1K-1)

CMoi

specific heat of moisture (J kg-1K-1)

d

equivalent particle diameter (m)

D

diffusion coefficient (m2 s-1)

R

universal gas constant (J mol-1 K-1)

ur g

gravitational acceleration (m s-2)

h

specific enthalpy (J kg-1)

λ

thermal conductivity (Wm-1K-1)

Heva

evaporation heat of the solid material (J kg-1)

18 ACS Paragon Plus Environment

Page 18 of 45

Page 19 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Kc

equilibrium constant

k

heat-transfer coefficient (W m-2 K-1)

km

mass-transfer coefficient (m s-1)

kri

kinetics rate of heterogeneous reaction I (m s-1)

K

interphase momentum exchange coefficient (kg m-3 s-1)



m

mass-transfer rate (kg m-3 s-1)

Nu

Nusselt number

Mi

molar weight (kg kmol-1)

WV

volatile mass fraction of briquette

p

pressure (Pa)

q

heat flux (W m-2)

Q

intensity of heat exchange (W m-3)

ri

reaction rate (kmol m-3 s-1)

rmi

heterogeneous reaction rate (kmol m-3 s-1)

rki

kinetic rate of heterogeneous reaction i (m s-1)

R

universal gas constant (J mol-1 K-1)

S

source term

Sh

Sherwood number

Pr

Prandtl number

e

aspect ratio

t

time (s)

T

temperature(K)

r v

velocity (m/s)

19 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

CD

drag force coefficient

vi

stoichiometric coefficient of reactant i

Yi

mass fraction of the ith species

Pc

packing density

y

bed height (m)

Greek Symbols α

volume fraction

ρ

density (kg m-3)

µ

dynamic viscosity (Pa s)

φ

angle of internal friction

λ

thermal conductivity (W m-1 K-1)

τ

stress tensor (Pa)

Subscripts

C

fixed carbon

CO

carbon monoxide

CH4

methane

CO2

carbon dioxide

volatile

volatile matter

primary tar

primary tar

g

gas phase

ra d

radiation

H2

hydrogen

H2O

vapour

20 ACS Paragon Plus Environment

Page 20 of 45

Page 21 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Moi

moisture

s

solid phase

Vi

mass fraction of solid species

390

ACKNOWLEDGEMENTS

391

The authors wish to thank Australian Research Council (LP150100112) and Coal Energy

392

Australia for the financial support of this project. Sino Steel is acknowledged for providing

393

the experimental data used in the model validation.

394

REFERENCES

395

(1)

396

Liang, P.; Zhang, Y.; Jiang, W.; Wei, A.; Liu, T.; Wu, J. Energy and Fuels 2015, 29 (11), 7727–7733.

397

(2)

Song, H.; Liu, G.; Wu, J. Energy Conversion and Management 2016, 126, 1037–1046.

398

(3)

Yu, J.; Tahmasebi, A.; Han, Y.; Yin, F.; Li, X. Fuel Processing Technology 2013, 106,

399 400

9–20. (4)

401

Yu, Y.; Xu, M.; Yao, H.; Yu, D.; Qiao, Y.; Sui, J.; Liu, X.; Cao, Q. Proceedings of the Combustion Institute 2007, 31 II, 1947–1954.

402

(5)

Liao, J.; Yu, A. B.; Shen, Y. Powder Technology 2017, 314, 550–556.

403

(6)

Katalambula, H.; Gupta, R. Energy & Fuels 2009, 23, 3392–3405.

404

(7)

Li, X.; Ashida, R.; Makino, M.; Nishida, A.; Yao, H.; Miura, K. Energy and Fuels

405 406 407

2014, 28 (9), 5690–5695. (8)

Osman, H.; Jangam, S. V.; Lease, J. D.; Mujumdar, A. S. Drying Technology 2011, 29 (15), 1763–1783.

21 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

408

(9)

Guy, P. J.; Perry, G. J. Fuel 1992, 71 (10), 1083–1086.

409

(10)

Stevenson, G. G.; Perlack, R. D. Energy Policy 1989, 17 (3), 215–227.

410

(11)

Zhu, X.; Sheng, C. Fuel Processing Technology 2010, 91 (8), 837–842.

411

(12)

Li, C.-Z. Fuel 2007, 86 (12), 1664–1683.

412

(13)

Hayashi, J. I.; Takahashi, H.; Iwatsuki, M.; Essaki, K.; Tsutsumi, A.; Chiba, T. Fuel

413

2000, 79 (3–4), 439–447.

414

(14)

Yip, K.; Wu, H.; Zhang, D.-K. Energy & Fuels 2007, 21, 419–425.

415

(15)

Schwarz, M. P. 10th International Conference on CFD in Oil & Gas, Metallurgical

416

and Process Industries 2014.

417

(16)

Huang, Q.; Zhang, Y.; Yao, Q.; Li, S. Powder Technology 2017, 317, 449–457.

418

(17)

Liu, C.; Liu, Z.; Zhang, T.; Huang, X.; Guo, J.; Zheng, C. Energy and Fuels 2017, 31

419

(3), 2596–2606.

420

(18) Lisandy, K. Y.; Kim, J. W.; Lim, H.; Kim, S. M.; Jeon, C. H. Fuel 2016, 185, 478–490.

421

(19)

Lin, W.; Feng, Y.; Zhang, X. Applied Thermal Engineering 2015, 81, 353–358.

422

(20)

Guo, Z.; Tang, H. China Particuology 2005, 3 (6), 373–378.

423

(21)

Slupik, L.; Fic, A.; Bulinski, Z.; Nowak, A. J.; Kosyrczyk, L.; Labojko, G. Fuel 2015,

424 425

150, 415–424. (22)

426 427

Zhong, W.; Yu, A.; Zhou, G.; Xie, J.; Zhang, H. Chemical Engineering Science 2016, 140, 16–43.

(23)

Zhang, Q.; Dor, L.; Fenigshtein, D.; Yang, W.; Blasiak, W. Applied Energy 2012, 90

22 ACS Paragon Plus Environment

Page 22 of 45

Page 23 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

428

(1), 106–112.

429

(24)

Zhang, Q.; Dor, L.; Yang, W.; Blasiak, W. Energy and Fuels 2011, 25 (9), 4129–4137.

430

(25)

Mahmoudi, A. H.; Hoffmann, F.; Peters, B. Journal of Analytical and Applied

431

Pyrolysis 2014, 106, 9–20.

432

(26)

433

(27) Wang, T.; Tang, T.; He, Y.; Yi, H. Chemical Engineering Journal 2016, 283, 127–140.

434

(28)

Q. F. Hou, Z. Y. Zhou, and A. B. Y. AIChE Journal 2012, 58 (5), 1422–1434.

435

(29)

Zhou, Z.; Zou, R.; Pinson, D.; Yu, A. AIP Conference Proceedings 2013, 1542, 357–

436 437

360. (30)

438 439

(31)

(32) Donev, A.; Stillinger, F. H.; Chaikin, P. M.; Torquato, S. Physical Review Letters 2004, 92(25 I) (255506–1). (33)

444 445

448

Zhou, Z. Y.; Pinson, D.; Zou, R. P.; Yu, A. B. Chemical Engineering Science 2011, 66 (23 (2011)), 6128–6145.

(34) Gera, D.; Mathur, M. P.; Freeman, M. C.; Robinson, A. Energy and Fuels 2002, 16 (6),

446 447

Li, C.; Zhou, Z.; Zou, R.; Pinson, D.; Yu, A. AIP Conference Proceedings 2013, 1542, 767–770.

442 443

Zhou, Z. Y.; Zou, R. P.; Pinson, D.; Yu, A. B. Industrial and Engineering Chemistry Research 2011, 50 (16), 9787–9798.

440 441

Ku, X.; Li, T.; Løvås, T. Chemical Engineering Science 2015, 122, 270–283.

1523–1532. (35)

J.T. Jenkins, S. C. C. Mechanics Applied to the Transport of Bulk Materials 1979, 79– 89. 23 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

449

(36)

Johnson, P. C.; Jackson, R. Journal of Fluid Mechanics 1987, 176, 67–93.

450

(37)

Zhang, Q.; Dor, L.; Biswas, A. K.; Yang, W.; Blasiak, W. Fuel Processing Technology

451 452

2013, 106, 546–554. (38)

453

Gidaspow, D.; Bezburuah, R.; Ding, J. Proceedings of the Seventh Engineering Foundatian Conference on Fluidization 1992, 75–82.

454

(39)

Kishore, N.; Gu, S. Chemical Engineering and Technology 2011, 34 (9), 1551–1558.

455

(40)

Yang, W.; Ponzio, A.; Lucas, C.; Blasiak, W. Fuel Processing Technology 2006, 87

456 457

(3), 235–245. (41)

458 459

Science 1993, 19 (6), 505–586. (42)

460 461

Hobbs, M. L.; Radulovic, P. T.; Smoot, L. D. Progress in Energy and Combustion

Badzioch, S.; Hawksley, P. Industrial & Engineering Chemistry Process Design and Development 1970, 9 (4), 521–530.

(43)

462

Kobayashi, H.; Howard, J. B.; Sarofim, A. F. Symposium (International) on Combustion 1977, 16 (1), 411–425.

463

(44)

Donald B. Anthony, Jack B. Howard, H. C. H. and H. P. M. Fuel 1976, 55, 121.

464

(45)

Niksa, S. Energy &Fuels 1991, 5, 673–683.

465

(46)

Grant, D. M.; Pugmire, R. J.; Fletcher, T. H.; Kerstein, A. R. Energy & Fuels 1989, 3

466 467

(2), 175–186. (47)

468 469

Solomon, P. R.; Hamblen, D. G.; Carangelo, R. M.; Serio, M. A.; Deshpande, G. V. ACS Division of Fuel Chemistry, Preprints 1987, 32 (3), 83–98.

(48)

Gregory, D. R.; Littlejohn, R. F. Br. Coal Util. Res. Assoc., Mon. Bull. 1965, 29:6. 24 ACS Paragon Plus Environment

Page 24 of 45

Page 25 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

470

(49)

Chan, W. C. R.; Kelbon, M.; Krieger, B. B. Fuel 1985, 64 (11), 1505–1513.

471

(50)

Boroson, M. L.; Howard, J. B.; Longwell, J. P.; Peters, W. A. AIChE Journal 1989, 35

472

(1), 120–128.

473

(51)

Grebenshchikova, G. Podzemnaya Gasifikatsiya Uglei 1957.

474

(52)

Lowry H.H. Chemistry of coal utilization; 1981; Vol. 2.

475

(53)

Jones, W. P.; Lindstedt, R. P. Combustion and Flame 1988, 73 (3), 233–249.

476

(54)

Yoon, H.; Wei, J.; Denn, M. AIChE Journal 1978, 24 (5), 885–903.

477

(55)

Hobbs, M. .; Radulovic, P. .; Smoot, L. D. AIChE Journal 1992, 38 (5), 681–702.

478 479 480 481 482 483 484

25 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

485 486

Figure 1. Schematic of the pilot-scale pyrolysis test system.

487

26 ACS Paragon Plus Environment

Page 26 of 45

Page 27 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

488 489

Figure 2. Dimensions of the pyrolyzer chamber (unit: mm) and position of the monitor point

490

in the middle of the chamber (red dot).

491 492

27 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

493 494

Figure 3. Particle size of an ellipsoidal briquette.

495 496 497 498 499

28 ACS Paragon Plus Environment

Page 28 of 45

Page 29 of 45

0.60 Packing density Polynomial Fit of "Packing Density"

0.59

Packing density

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

0.58 0.57 0.56 0.55 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40

500 501

Bed height (m)

Figure 4. Packing density distribution of the ellipsoids packed bed in the pyrolyzer.

502 503 504 505 506 507 508 509

29 ACS Paragon Plus Environment

Energy & Fuels

Experimental data Predict data

1000 900

Temperature (K)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 45

800 700 600 500 400 300 0

2000

4000

6000

8000

10000 12000

Time (s)

510 511

Figure 5. Comparison of temperature history between simulation results and experiment

512

measurements.

513 514 515 516 517 518 519

30 ACS Paragon Plus Environment

Page 31 of 45

35

2.1 % moisture Experimental data

35.2 32 29.6 27.6

30

15

25 20

10

18.5 16.6

15.915.1

15

5

10 5

2.9

1.1 1.5

4

0

0 CH4

CO2

H2

NH3

1.9

CO

C2H6

Remaining VM mass fraction %

20

40

2.1 % moisture

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

0

VM

Experiment

520 521

Figure 6. Comparison of gas yields between simulation results and experiment measurements.

522 523 524 525 526 527 528 529 530

31 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(a)

(b)

Page 32 of 45

(c)

(d)

531 532

Figure 7. Temperature profiles of the briquettes (upper row) and gas phase (bottom row) at

533

selected time in the base case: (a) t=2500s; (b) t=5000s; (c) t=7500s; (d) t=10000s.

534 535 536 537 538

32 ACS Paragon Plus Environment

Page 33 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

(a)

(b)

(c)

(d)

539 540

Figure 8. Evolution of briquettes’ remaining volatile matte, primary tar, and secondary tar in

541

the base case: (a) t=2500s; (b) t=5000s; (c) t=7500s; (d) t=10000s.

542

33 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 45

543 544

(a)

(b)

(c)

(d)

545

Figure 9. Evolution of gas flow field in the base case: (a) t=2500s; (b) t=5000s; (c) t=7500s;

546

(d) t=10000s.

547 548

34 ACS Paragon Plus Environment

Page 35 of 45

2.5x10

Mass flow rate (kg/s)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

-3

2.0x10-3 1.5x10-3 1.0x10-3 5.0x10-4

0.0 0

2500

5000

7500

10000

12500

15000

Time (s)

549 550

Figure 10. Mass flow rate versus time at the outlet in the base case.

551 552 553 554 555 556 557 558 559 560

35 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

561 562

Page 36 of 45

(a)

(b)

(c)

(d)

Figure 11. Profiles of some gas species at 7500s in the base case: (a), CH4; (b), CO2 ; (c), H2 and (d), H2O.

563 564 565

36 ACS Paragon Plus Environment

Page 37 of 45

1100

35

1000

30

900 800

Mole fraction (% )

Tem perature (k)

566

Case 1 Case 2 Case 3 Case 4 Case 5

700 600 500

CH4 CO2 H2 NH3 CO C2H6

25 20 15 10 5

400

0

300 0

2000 4000 6000 8000 10000 12000 14000

Case 1

Case 2

Case 3

Case 4

Case 5

Time (s)

(b)

(a)

3.0x10-3

0.30 0.25

Mass flow rate (kg/s)

Case 1 Case 2 Case 3 Case 4 Case 5

0.35

VM m ass fraction

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

0.20 0.15 0.10

Case 1 Case 2 Case 3 Case 4 Case 5

2.0x10-3

-3

1.0x10

0.05 0.00 0

2000 4000 6000 8000 10000 12000 14000

0.0

Time (s)

0

2000 4000 6000 8000 10000 12000 14000 Time (s)

(d)

(c) 567 568 569

Figure 12. Comparison of Cases 1-5 with different initial moisture contents of the briquettes, in terms of (a), temperature history, (b), gas products composition, (c), remaining volatile matter in briquettes, and (d), mass flow rate at the outlet.

570 571

37 ACS Paragon Plus Environment

Energy & Fuels

0.60 Case 1 Case 8 Case 9 Case 10 Case 11

0.59

Packing density

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 45

0.58 0.57 0.56 0.55 0.54 0.53 0.00

572 573

0.05

0.10

0.15

0.20

0.25

0.30

0.35

0.40

Bed height (m)

Figure 13. Packing density distribution for briquettes with different aspect ratio.

574 575 576 577

38 ACS Paragon Plus Environment

Page 39 of 45

1100

35

1000

30

900

Case Case Case Case Case

800 700 600

Mole fraction (%)

Temperature (K)

578

8 9 1 10 11

500

CH4 CO2 H2 NH3 CO C2H6

25 20 15 10 5

400 300

0 0

2000

4000

6000

Case 8 Case 9 Case 1 Case 10 Case 11 Cases with different aspect ratio

8000 10000 12000 14000

Time (s)

(a)

(b)

0.35

2.5x10

0.25 0.20

Mass flow rate (k g/s)

C ase 8 C ase 9 C ase 1 C ase 10 C ase 11

0.30

VM mass fraction

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

0.15 0.10

-3

Case Case Case Case Case

2.0x10-3

8 9 1 10 11

1.5x10-3 1.0x10-3 5.0x10-4

0.05 0.0

0.00 0

579 580 581

2000

4000

6000

0

8000 10000 12000 14000

2000 4000 6000 8000 10000 12000 14000

Time (s)

Tim e (s)

(c)

(d)

Figure 14. Comparison of Cases 1, 8-11 with different briquette aspect ratio, in terms of (a), temperature history, (b), gas products composition, (c), remaining volatile matter in briquettes, and (d), mass flow rate at the outlet.

582

39 ACS Paragon Plus Environment

Energy & Fuels

1100

35

1000

30

900

Mole fraction (%)

Temperature (k)

583

800 700

Case 6 Case 1 Case 7

600 500

CH4 CO2 H2 NH3 CO C2H6

25 20 15 10 5

400

0

300 0

2000

4000

6000

Case 6

8000 10000 12000 14000

Case 1

Case 7

Time (s)

(b)

(a)

0.35

3.0x10-3

Case 6 Case 1 Case 7

0.30

Mass flow rate (kg/s)

Mass fraction of VM

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 45

0.25 0.20 0.15 0.10

2.5x10

Case 6 Case 1 Case 7

-3

2.0x10-3

1.5x10

-3

1.0x10-3 5.0x10-4

0.05 0.00

0.0

0

2000 4000 6000 8000 10000 12000 14000

0

Time (s)

4000

6000

8000

10000

12000

14000

Time (s)

(c)

584 585 586

2000

(d)

Figure 15. Comparison of Cases 1, 6, 7 with different final pyrolysis temperature, in terms of (a), temperature history, (b), gas products composition, (c), remaining volatile matter in briquettes, and (d), mass flow rate at the outlet.

587 588 589

40 ACS Paragon Plus Environment

Page 41 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

590

Energy & Fuels

Table 1. Properties of the briquette investigated. Proximate Analysis (on a Dry Basis, except Moisture)

Moisture (%)

2.1

Fixed carbon (%)

61.2

Volatile (%)

32.7

Ash (%)

6.1 Ultimate Analysis (on a Dry basis)

Carbon (%)

81.1

Hydrogen (%)

4.7

Oxygen (%)

7.2

Nitrogen (%)

0.7

Sulphur (%)

0.2

591 592 593 594

41 ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

595

Page 42 of 45

Table 2. Governing equations of gas-solid flows in the CFD model. Gas phase: r ∂ (α g ρ g Yi ) + ∇ (α g ρ g Yi v g ) = m& i + Si ∂t

(1)

r r r ur r r r ∂ (α g ρ g v g ) + ∇ (α g ρ g v g v g ) = −α g ∇p + ∇τ g + α g ρ g g + K sg (v s − v g ) − m& sg v sg ∂t

(2)

r r r ∂ ∂p (α g ρ g hg ) + ∇ (α g ρ g v g hg ) = −α g + τ g : ∇ v g − ∇ q g + S g + Qsg + m& sg hsg ∂t ∂t

(3)

Solid phase: r ∂ (α s ρ sY j ) + ∇ (α s ρ sY j v s ) = m& j + S j ∂t

(4)

r r r ur r ∂ (α s ρ s v s ) + ∇ (α s ρ s v s v s ) = −α s ∇ p − ∇ps + ∇τ s + α s ρ s g − m& sg v sg ∂t

(5)

r r r ∂ ∂p (α s ρ s hs ) + ∇ (α s ρ s v s hs ) = −α s + τ s : ∇ v s − ∇ q s + S s + Qgs − m& sg hsg ∂t ∂t

(6)

596 597 598 599 600 601 602 603

42 ACS Paragon Plus Environment

Page 43 of 45 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

604

Energy & Fuels

Table 3. Reactions considered in this CFD model and their reaction kinetics.

Reaction

Reaction rate expression and kinetics

Coal → αPrimary tar +

r1 = 3.2 × 10 5 (1 − α g ) exp(

Reference 48,49

− 1.6 × 10 4 ) ρ volatile Ts

βGas + γChar T1=

2 6 7 .4 + 2 7 3 .1 5 W v 0 .25 5 5 50

Primary tar → Secondary tar + gas

r2 = 9.55 × 10 4 α g exp(

r3 = α g 0.03exp( r3 CO + H 2 O ← → CO 2 + H 2

4 C H 4 + H 2 O  r → CO + 3 H 2

1.12 × 10 4 ) ρ prim ary tar Tg

−6.03 × 10 4 CCO C H 2O ) − RTg M CO M H 2O

α 0.03exp( g

Kc

51,52

M CO2 M H 2

 1 91500  (11321 − 31.08Tg + 3Tg ln Tg 2.8 ×10−4 Tg 2 − ) Kc = exp  Tg   RTg

r4 = 3.0 × 108 × α gTg exp( −15083 / T g )CCH 4 C H 2O

km =

53

54

 1  Avρ i rmi =    v iM i  1 + 1 k m rki m1 C + H 2 O  r → CO + H 2

−6.03 × 10 4 ) CCO 2 C H 2 RTg

ShD ds

rk1 = 5.714T s exp( − 15600 / T s ) m2 C + CO 2  r→ 2CO

rk2 = 5.89 × 10 2 Ts exp( −26800 / Ts )

605 606 607

43 ACS Paragon Plus Environment

55

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

608

Page 44 of 45

Table 4. Boundary and operating conditions used in this model. Property

value

Thermal conductivity of

if Ts ≤ 673 λ = 0.23  −5 1.8 λ = 0.23 + 2.24 × 10 × (Ts − 673) if Ts > 673

briquettes

 C P = V1 C F + V 2 C V + V 3 C A + V 4 C M o i  2 −3 −6  C F = − 0.281 + 3.8 07 × 1 0 T s − 1.758 × 1 0 T s  −3  C V = 0.728 + 3.39 1 × 10 T s  −4  C A = 0.594 + 5.8 6 × 10 T s 3 C  M oi = 4.2 × 10

Specific heat of briquettes

Mass weighted

ds = 3.5cm

Equal volume briquette diameter

Wall temperature scheme 1

Wall temperature scheme 2

Wall temperature scheme 3

Twall = 873.15  Twall = 873.15 − 0.055 × ( t − 1800 )  Twall = 773.15 + 0.055 × (t − 3600) Twall = 873 

if t < 1800s

if 1800s ≤ t ≤ 3600s if 3600s