Research articleReaxFF-based molecular dynamics simulation of the initial pyrolysis mechanism of lignite
Introduction
Coal, which is the principal traditional source of energy, occupies a dominant position in world energy consumption, especially in China [1]. With the amount consumption of high rank coal, lignite has attracted considerable interest due to its vast reserves, low selling price, high chemical reactivity, and low pollution-forming impurities [2,3]. However, self-defection (e.g., low calorific value, high moisture content, and a high tendency of spontaneous combustion) makes it unfit for direct utilization [4]. Upgrading lignite to higher quality products by means of pyrolysis is considered to be a promising way to achieve economically efficient utilization of lignite [5,6]. Solid residues (coke or char) with high calorific value can be utilized for direct combustion or preparation of slurry fuels [7]. Pyrolysis gases contain large amounts of combustible gases, such as H2, CO, and hydrocarbon gas [8]. Tars produced from lignite pyrolysis are suitable for liquid fuels and valuable chemicals [9].
Numerous efforts have been devoted to study the pyrolysis behavior of lignite. Liu et al. [10] investigated the interrelation between lignite structural characteristics and product distributions. They concluded that the formation of pyrolysis products is related to the breakage of chemical bonds and the decomposition of functional groups. Xu et al. [11] reported pyrolysate distributions and kinetics characteristics of Zhaotong lignite. They found the optimal reaction conditions for generation of tar and pyrolysis gas. Meng et al. [12] studied the morphology and chemical structure of low-temperature pyrolysis chars. They found that temperature has an apparent influence on product yields, surface morphology, and the evolution of chemical functional groups of chars. He et al. [13] analyzed the thermal degradation of Shengli lignite using TG–GC–MS. The results indicated that benzene series are the major gaseous products, and small molecular structures are the main components of aliphatic hydrocarbons. According to experimental studies, some progress has been made in understanding the characteristics of pyrolysis products and the chemical reaction mechanism of lignite pyrolysis. However, it is quite difficult to understand the complicated thermal decomposition reactions in depth through experimental methods alone because of the heterogeneity of lignite and the complexity of the pyrolysis process. Furthermore, vast free radicals are generated in an exceedingly short time, and these are difficult to detect in laboratories [14]. Undoubtedly, computational approaches would provide a promising platform for further researching the lignite pyrolysis mechanism.
Quantum chemistry (QC) modeling is applicable for studying chemical reactions with high accuracy. However, QC-based methods are computationally expensive and intensive, thus they have rarely been applied in models with more than 100 atoms [15]. Classical molecular dynamics (MD) has the ability to model molecular systems that contain more than 10,000 atoms, but it is not suitable for exploring the cleavage and production of bonds [16]. Fortunately, the Reactive Force Field (ReaxFF), introduced by van Duin et al., can address the chemical interactions of atoms and molecules [17]. ReaxFF-MD combines the advantages of MD and ReaxFF, which can be used to study the chemical reactions of large and complicated molecular systems. Additionally, ReaxFF-MD has been proven to maintain basically the calculation precision of QC but with much reduced computational costs [18,19].
In the last decade, ReaxFF-MD has been successfully used to simulate the combustion and pyrolysis characteristics of complex compounds, such as coal [20], biomass [21], oil shale kerogen [22] and char [23], with the goal of trying to understand the reaction process in depth. Zhan et al. [24] investigated the reaction mechanism of subbituminous coal via ReaxFF-MD. The calculation results indicated that the pyrolysis process is initiated by breaking unstable CC and CO bonds, which is followed by intramolecular hydrogen transformation. Castro-Marcano et al. [25] adopted ReaxFF-MD to study the combustion processes of coal char. They found that aromatic structures are more liable to oxidize and combust under fuel lean and higher temperature conditions. Zheng et al. [26] employed ReaxFF-MD to analyze the overall reaction stages for Liulin bituminous coal. The simulation results showed that the pyrolysis processes are roughly divided into four stages according to the breakage behavior of chemical bonds. These four stages are the activation stage, primary pyrolysis stage, secondary pyrolysis stage, and recombination dominated stage. Chen et al. [27] studied the pyrolysis and combustion behavior of biomass using ReaxFF-MD. The simulation results indicated that environment and temperature play a significant influence on the dissociation of chemical bonds and product distributions. Extensive studies of the chemical reaction mechanism using ReaxFF-MD simulations have been reported. Nevertheless, coal type is a critical factor that greatly affects the pyrolysis process. Compared with the structural features of bituminous coal and subbituminous coal, the molecular structure of lignite contains more oxygen-containing functional groups, aliphatic side chains, and single aromatic ring structures [28]. Published studies on the pyrolysis behavior of lignite are rare [29]. Further research is needed to better understand the complex pyrolysis process, especially the relationship between pyrolysis products and lignite structure.
The main purpose of this work is to explore the pyrolysis behavior of lignite using ReaxFF-MD simulations. The paper is organized as follows: First, the characteristics of pyrolysis products, especially tar and gas, are analyzed in detail. Second, the transformation behavior of carbon, hydrogen, and oxygen during the pyrolysis process are revealed. Finally, two typical temperatures (2000 K and 3000 K) are selected for an in-depth investigation of the pyrolysis reaction mechanism. The elucidation of these problems may provide more constructive information of high-efficiency clean technology for converting lignite.
Section snippets
Computational details
ReaxFF [17] is an empirical bond-order-based reactive force field and can explicitly describe chemical reactions within complex systems. The relationship between bond order and bond energy plays an important role in ReaxFF. Bond order is got from interatomic distances and is updated continually at every iteration, thus ensuring connectivity changes. The bonded interactions (i.e., bonds, angles, and torsions) are bond-order dependent so that energies related to these terms disappear upon bond
Analysis of lignite pyrolysis products
Undoubtedly, temperature is a critical factor that affects the generation and evolution of pyrolysis products. To evaluate the influence of temperature on product distributions, a set of ReaxFF-MD constant temperature simulations were conducted at temperatures 1600–3000 K for 250 ps. To be consistent with previous work [39,40], pyrolysis products from ReaxFF-MD simulations are classed as five types: inorganic gas, organic gas, heavy tar, light tar, and char. Inorganic small molecules (e.g., H2,
Conclusions
In this work, the initial pyrolysis mechanism of lignite was investigated using ReaxFF-MD constant temperature simulations at 1600–3000 K for 250 ps. The results show that lignite pyrolysis process begins with the dissociation of the weaker bridge bonds such as Cal-O and Cal-Cal. Moreover, increasing temperature is favorable for the decomposition of macromolecular structure and the breakage of bridge bonds at the lower temperatures (1600–2200 K). At relatively high temperatures (2200–3000 K),
Declaration of Competing Interest
The authors declare no competing financial interest.
Acknowledgments
The authors gratefully acknowledge the financial support received from the National Natural Science Foundation of China (Grant No. 51676032) and Program for Changjiang Scholars and Innovative Research Team in University (No. IRT13052).
References (52)
- et al.
Correlations of chemical properties of high-alkali solid fuels: a comparative study between Zhundong coal and biomass
Fuel
(2018) - et al.
Improving gas/particle flow deflection and asymmetric combustion of a 600 MWe, supercritical down-fired boiler by increasing its upper furnace height
Energy
(2017) - et al.
Characteristics comparison of tar from lignite pyrolysis with inherent and simulated moisture for adopting a two-stage gasification process
Fuel
(2019) - et al.
Suppression of crosslinking combination of carboxyl functional groups with NaOH on combustion performance of Shengli lignite
Fuel
(2018) - et al.
Effect of pyrolysis temperature on lignite char properties and slurrying ability
Fuel Process. Technol.
(2015) - et al.
Gas evolution characteristics of lignite during low-temperature pyrolysis
J. Anal. Appl. Pyrolysis
(2013) - et al.
The structure and pyrolysis product distribution of lignite from different sedimentary environment
Appl. Energy
(2016) - et al.
Low temperature pyrolysates distribution and kinetics of Zhaotong lignite
Energy Convers. Manag.
(2016) - et al.
TG–GC–MS study of volatile products from Shengli lignite pyrolysis
Fuel
(2015) - et al.
DFT study on the sulfur migration during benzenethiol pyrolysis in coal
J. Mol. Struct. (THEOCHEM)
(2010)
Ignition and combustion of nano-sized aluminum particles: a reactive molecular dynamics study
Combust. Flame
Initial reaction mechanisms of cellulose pyrolysis revealed by ReaxFF molecular dynamics
Fuel
Primary understanding of non-isothermal pyrolysis behavior for oil shale kerogen using reactive molecular dynamics simulation
Int. J. Hydrog. Energy
Preliminary understanding of initial reaction process for subbituminous coal pyrolysis with molecular dynamics simulation
Fuel
Combustion of an Illinois No. 6 coal char simulated using an atomistic char representation and the ReaxFF reactive force field
Combust. Flame
The molecular representations of coal–a review
Fuel
Early maturation processes in coal. Part 2: reactive dynamics simulations using the ReaxFF reactive force field on Morwell Brown coal structures
Org. Geochem.
Experimental and dynamics simulation studies of the molecular modeling and reactivity of the Yaojie oil shale kerogen
Fuel
Sulfur removal from petroleum coke during high-temperature pyrolysis, analysis from TG-MS data and ReaxFF simulations
J. Anal. Appl. Pyrolysis
Investigation of N behavior during coal pyrolysis and oxidation using ReaxFF molecular dynamics
Fuel
Reactive molecular dynamics simulation and chemical kinetic modeling of pyrolysis and combustion of n-dodecane
Combust. Flame
ReaxFF simulations of hydrothermal treatment of lignite and its impact on chemical structures
Fuel
Theoretical study of the CO formation mechanism in the CO2 gasification of lignite
Fuel
Molecular dynamics simulations of Zhundong coal pyrolysis using reactive force field
Fuel
Lignite upgrading by multi-stage fluidized bed pyrolysis
Fuel Process. Technol.
Analysis of tars formed during co-pyrolysis of coal and biomass at high temperature in carbon dioxide atmosphere
J. Anal. Appl. Pyrolysis
Cited by (69)
Insights into the pyrolytic coking process of RP-3 fuel from ReaxFF molecular dynamics
2024, Chemical Engineering ScienceExperimental and reactive molecular dynamics simulation of municipal sludge extract pyrolysis
2024, Journal of the Energy InstituteTransformation simulation of N-containing functional groups in coal pyrolysis and combustion processes by using ReaxFF
2024, Chemical Engineering ScienceMolecular structure model construction and pyrolysis mechanism study on low-rank coal by experiments and ReaxFF simulations
2024, Journal of Analytical and Applied Pyrolysis