ArticlePyrolysis of vulcanized styrene-butadiene rubber via ReaxFF molecular dynamics simulation
Graphical abstract
Introduction
Styrene-Butadiene Rubber (SBR) is a general-purpose synthetic rubber and a main component in passenger car tires. The production of scrap tires increases with the development of automotive industry. However, only a small portion of scrap tires was recycled or reused, and the large part was treated unprofessionally and caused “black pollution”. Pyrolysis of scrap tires is an efficient and eco-friendly process that allows producing fuels and chemicals from scrap tires [1], [2], and also provides heat and electrical energy due to the high calorific value [1], [2], [3], [4], [5]. The products from thermal decomposition of scrap tires include gas and liquid products which can be used as energy sources, and char that can be transformed into carbon black after treatment. Additionally, sulfur is widely used in tires manufacture as an additive [1], [6]. The addition of sulfur forms a network structure by bringing in poly-, di-, and mono-sulfide bonds [7] to enhance the mechanical properties of rubber. In the pyrolysis process, the sulfur which is originally in organic forms releases into pyrolysis products such as oil, tar and gases. The removal of sulfur during rapid pyrolysis is of high importance for industrialization of scrap tires.
The sulfur usually remains in the oil after pyrolysis process resulting in low quality products [8], [9]. Thus, the desulfurization during pyrolysis is necessary to recycle vulcanized rubber effectively [10]. Previous studies indicate that sulfur is involved in a serial of reactions in the pyrolysis process [9], and sulfur removal is accomplished by H2S release [11]. The pyrolysis product is strongly dependent on the temperature [12], heating rate, residence time and mass transfer, which are related with the reactor configurations [9], [13]. Among these factors, the temperature is a common factor discussed in their works. However, it is still a challenge to examine the effect of a factor (independent variable) in the pyrolysis experiment. During the last decade, rapid development of computer power has made molecular simulation indispensable in the study on many chemical processes such as the thermal decomposition of organic materials.
Reactive force field [14], [15] molecular dynamics (ReaxFF MD) simulation as an efficient computation tool has been used in simulating the oxidation of hydrocarbons. Importantly, single factor experiment is able to be performed with the help of this theoretical method. The bond order is employed in this force field to determine the bond formation and breakage. Unlike other simulation methods, ReaxFF MD can simulate the formation of radicals, the bonding and breaking dynamically. The accuracy of ReaxFF force field was insured by quantum mechanics (QM) data, and ReaxFF MD can simulate large-scale systems with lower computational expense than QM methods [16]. The application of ReaxFF has been expanded constantly, which has been used in simulating the pyrolysis of complex compound such as coal [17], [18], [19] and biomass [20]. Lately, the natural rubber and its pyrolysis mechanism was studied through ReaxFF MD combined with thermogravimetry-infrared spectroscopy experiments [21]. The simulation results are in good agreement of experimental pyrolysis data. Therefore, ReaxFF MD simulation is practicable to explain the pyrolysis mechanism of vulcanized SBR at various reaction conditions.
Our previous work [22] used density function theory (DFT) and ReaxFF MD simulation to study bond dissociation energies of SBR polymer chain and the pyrolysis process of pure SBR. In this work, we build the SBR models with cross-linking structures firstly. And the ReaxFF MD simulation is used to investigate the pyrolysis mechanism of vulcanized SBR at different temperatures and densities. Especially, we focus on the transfer of sulfur during the pyrolysis which is key to the devulcanization [7], and the sulfur products distribution are well concluded. We intend to provide a detailed theoretic understanding and yield guidelines on the devulcanization during pyrolysis of vulcanized SBR.
Section snippets
Computational Method
The bond dissociation energies of C-S and S-S bonds are carried out using the Gaussian 09 series of programs [23]. All simulations were performed by using the MD software Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) software. Inspired by the cross-linking of polymers by united-atom model in MD simulation [24], [25], we built cross-linked models based on all-atom SBR chains. The hydrogen atoms must be taken into calculation to meet the need of pyrolysis simulations. And
The analysis of vulcanized SBR structures
The vulcanized SBR models were accomplished by a cross-linking process (Section 2.1) in MD simulations with UFF. The initial density of all samples is set as 0.1 g⋅cm−3. Note that the rubber is not fully filled in the pyrolysis reactor in the industrial production process, and thus, we prepare the samples with different densities to simulate different solid contents at the beginning. The SBR chains contain different number of sulfur atoms ranging from 1 to 8 for different vulcanized SBR
Conclusions
In this study, we built the cross-linked SBR model and performed the ReaxFF molecular dynamics simulation to study the pyrolysis process of the vulcanized SBR. The effect of temperature and pressure on the pyrolysis product distribution is greatly emphasized, especially the formation and transfer of sulfur containing compounds. Based on simulation results and comparison with corresponding experiments, the main conclusions are as follows:
- 1.
Sulfur as a cross-linking agent decreases the mobility of
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgements
The authors would like to express appreciation for the support of National Key Research and Development Program of China (Grant No. 2018YFC1902601).
References (37)
- et al.
Waste tyre pyrolysis – A review
Renew. Sustain. Energy Rev.
(2013) - et al.
Carbon black recovery from waste tire pyrolysis by demineralization: Production and application in rubber compounding
Waste Manage.
(2019) - et al.
Composition of oils derived from the batch pyrolysis of tyres
J. Anal. Appl. Pyrolysis.
(1998) - et al.
The production of fuel oil and combustible gas by catalytic pyrolysis of waste tire using waste heat of blast-furnace slag
Energy Convers. Manage.
(2017) - et al.
Fuel production from pyrolysis of natural and synthetic rubbers
Fuel
(2017) - et al.
Recycling of rubber wastes by devulcanization
Resour., Conserv. Recycl.
(2018) Kinetics of sulphur compounds in waste tyres pyrolysis
J. Anal. Appl. Pyrol.
(2000)- et al.
Clean pyrolysis oil from a continuous two-stage pyrolysis of scrap tires using in-situ and ex-situ desulfurization
Energy
(2017) - et al.
TG-FTIR and Py-GC/MS study on pyrolysis mechanism and products distribution of waste bicycle tire
Energy Convers. Manage.
(2018) - et al.
ReaxFF and DFT study on the sulfur transformation mechanism during the oxidation process of lignite
Fuel
(2016)
Combustion of an Illinois No. 6 coal char simulated using an atomistic char representation and the ReaxFF reactive force field
Combust. Flame
Studying the mechanisms of natural rubber pyrolysis gas generation using RMD simulations and TG-FTIR experiments
Energy Convers. Manage.
Crystallization of cross-linked polyethylene by molecular dynamics simulation
Polymer
Study of coal hydropyrolysis and desulfurization by ReaxFF molecular dynamics simulation
Fuel
Pyrolysis simulations of Fugu coal by large-scale ReaxFF molecular dynamics
Fuel Process. Technol.
A review of dipentene (dl-limonene) production from waste tire pyrolysis
J. Anal. Appl. Pyrolysis.
Formation of styrene monomer, dimer and trimer in the primary volatiles produced from polystyrene pyrolysis in a wire-mesh reactor
Fuel
Studies on the thermal degradation of cis-1,4-polyisoprene
Fuel
Cited by (13)
The mechanism of the influence of radiation heat flux on the combustion behavior of raw rubbers
2024, Polymer Degradation and StabilityEnd-of-life tyre conversion to energy: A review on pyrolysis and activated carbon production processes and their challenges
2023, Science of the Total EnvironmentModelling the dynamic physical properties of vulcanised polymer models by molecular dynamics simulations and machine learning
2023, Computational Materials ScienceThe reaction mechanism and sulfur evolution during vulcanized nature rubber pyrolysis in the atmosphere of H<inf>2</inf>O: A ReaxFF molecular dynamics study
2022, Polymer Degradation and StabilityCitation Excerpt :The bond order and non-bonded cutoffs were selected as 0.3 and 4.5 Å, respectively. The ReaxFF force field of the C/H/O/N/Si/S system was chosen for the systems under study, which has been successfully applied to the pyrolysis of vulcanized rubber [8,12]. To further verify the accuracy of the force field for this work, we calculated the activation energy of NR pyrolysis based on experimental and simulated data, which showed little difference, 164.75 kJ/mol and 156.72 kJ/mol, respectively.