Novel simplification approach for large-scale structural models of coal

Three-dimensional molecules to two-dimensional lattices. Part 3: Reactive lattice simulations

Yesica E. Alvarez, Brian M. Moreno, Michael T. Klein, Justin Kyle Watson, Fidel Castro-Marcano, Jonathan P. Mathews

Research output: Contribution to journalArticle

4 Citations (Scopus)

Abstract

Various studies have used lattice structures combined with simplistic kinetics to explore the devolatilization of coal with the goal of predicting yields of gas, tar, and char. In our previous work, we have demonstrated the ability to reduce a three-dimensional (3D) large-scale (>50 000 atoms) atomistic representation of Illinois no. 6 bituminous coal to a coal-specific two-dimensional (2D) lattice with cross-links being depicted as linkage lines and aromatic clusters being depicted as nodes. Because there is a direct link between the full complexity of the atomistic representation and the 2D lattice, the cross-links within the coal topography can be identified. Using this structural information, a statistical simulation of the kinetically controlled liquefaction reactions of an Illinois no. 6 coal lattice was performed. The Illinois no. 6 coal lattice was composed of a distribution of singular nodes, dimers, and multi-linkage molecules that ranged from 3 nodes (or clusters) to a 16-node entity. Probabilities for the cleavage of coal interunit links were calculated on the basis of the mixture of 10 model compounds comprised of two benzene rings connected by a cross-link (e.g., diphenyl sulfide, benzyl phenyl ether, diphenyl ether, etc.) equivalent to those found in the Illinois no. 6 model. The mixture of model compounds was based on the concentrations of cross-links representative of the model (and, thus, the coal analysis). Cross-link-specific kinetic parameters (Arrhenius pre-exponential factor and activation energy) for each pyrolysis reaction were obtained from literature sources or tuned from composition data using the Kinetic Modeling Editor, a software tool used to generate and solve the balance equations in a kinetic model. At a heating rate of 5 C/min, thermolysis calculations from 360 to 490 C showed that the 2D lattice structure broke down extensively, generating mainly monomer cluster molecules. At 400 C with residence times of 600 and 1600 s, the number of cross-links between aromatic units was reduced by 40 and 60%, respectively, from the original population. At higher temperatures, the breakdown of the weaker and more kinetically favored cross-links (-CH 2S- and -OCH2-) occurred rapidly at the earlier stages of heating. The absolute probability of cleavage (Pi) of certain cross-links (-CH2CH2-, -CH2OCH2-, etc.) increased at longer residence times. Cross-links such as phenyl and doubly aromatic bound (-S- and -O-) remained essentially unreactive at the conditions explored.

Original languageEnglish (US)
Pages (from-to)2915-2922
Number of pages8
JournalEnergy and Fuels
Volume27
Issue number6
DOIs
StatePublished - Jun 20 2013

Fingerprint

Coal
Molecules
Kinetics
Ethers
Thermolysis
Bituminous coal
Tars
Tar
Liquefaction
Heating rate
Kinetic parameters
Dimers
Topography
Benzene
Pyrolysis
Activation energy
Monomers
Gases
Heating
Atoms

All Science Journal Classification (ASJC) codes

  • Chemical Engineering(all)
  • Energy Engineering and Power Technology
  • Fuel Technology

Cite this

@article{cda86788a95b4495b5fd46175402ce25,
title = "Novel simplification approach for large-scale structural models of coal: Three-dimensional molecules to two-dimensional lattices. Part 3: Reactive lattice simulations",
abstract = "Various studies have used lattice structures combined with simplistic kinetics to explore the devolatilization of coal with the goal of predicting yields of gas, tar, and char. In our previous work, we have demonstrated the ability to reduce a three-dimensional (3D) large-scale (>50 000 atoms) atomistic representation of Illinois no. 6 bituminous coal to a coal-specific two-dimensional (2D) lattice with cross-links being depicted as linkage lines and aromatic clusters being depicted as nodes. Because there is a direct link between the full complexity of the atomistic representation and the 2D lattice, the cross-links within the coal topography can be identified. Using this structural information, a statistical simulation of the kinetically controlled liquefaction reactions of an Illinois no. 6 coal lattice was performed. The Illinois no. 6 coal lattice was composed of a distribution of singular nodes, dimers, and multi-linkage molecules that ranged from 3 nodes (or clusters) to a 16-node entity. Probabilities for the cleavage of coal interunit links were calculated on the basis of the mixture of 10 model compounds comprised of two benzene rings connected by a cross-link (e.g., diphenyl sulfide, benzyl phenyl ether, diphenyl ether, etc.) equivalent to those found in the Illinois no. 6 model. The mixture of model compounds was based on the concentrations of cross-links representative of the model (and, thus, the coal analysis). Cross-link-specific kinetic parameters (Arrhenius pre-exponential factor and activation energy) for each pyrolysis reaction were obtained from literature sources or tuned from composition data using the Kinetic Modeling Editor, a software tool used to generate and solve the balance equations in a kinetic model. At a heating rate of 5 C/min, thermolysis calculations from 360 to 490 C showed that the 2D lattice structure broke down extensively, generating mainly monomer cluster molecules. At 400 C with residence times of 600 and 1600 s, the number of cross-links between aromatic units was reduced by 40 and 60{\%}, respectively, from the original population. At higher temperatures, the breakdown of the weaker and more kinetically favored cross-links (-CH 2S- and -OCH2-) occurred rapidly at the earlier stages of heating. The absolute probability of cleavage (Pi) of certain cross-links (-CH2CH2-, -CH2OCH2-, etc.) increased at longer residence times. Cross-links such as phenyl and doubly aromatic bound (-S- and -O-) remained essentially unreactive at the conditions explored.",
author = "Alvarez, {Yesica E.} and Moreno, {Brian M.} and Klein, {Michael T.} and Watson, {Justin Kyle} and Fidel Castro-Marcano and Mathews, {Jonathan P.}",
year = "2013",
month = "6",
day = "20",
doi = "10.1021/ef4001105",
language = "English (US)",
volume = "27",
pages = "2915--2922",
journal = "Energy & Fuels",
issn = "0887-0624",
publisher = "American Chemical Society",
number = "6",

}

Novel simplification approach for large-scale structural models of coal : Three-dimensional molecules to two-dimensional lattices. Part 3: Reactive lattice simulations. / Alvarez, Yesica E.; Moreno, Brian M.; Klein, Michael T.; Watson, Justin Kyle; Castro-Marcano, Fidel; Mathews, Jonathan P.

In: Energy and Fuels, Vol. 27, No. 6, 20.06.2013, p. 2915-2922.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Novel simplification approach for large-scale structural models of coal

T2 - Three-dimensional molecules to two-dimensional lattices. Part 3: Reactive lattice simulations

AU - Alvarez, Yesica E.

AU - Moreno, Brian M.

AU - Klein, Michael T.

AU - Watson, Justin Kyle

AU - Castro-Marcano, Fidel

AU - Mathews, Jonathan P.

PY - 2013/6/20

Y1 - 2013/6/20

N2 - Various studies have used lattice structures combined with simplistic kinetics to explore the devolatilization of coal with the goal of predicting yields of gas, tar, and char. In our previous work, we have demonstrated the ability to reduce a three-dimensional (3D) large-scale (>50 000 atoms) atomistic representation of Illinois no. 6 bituminous coal to a coal-specific two-dimensional (2D) lattice with cross-links being depicted as linkage lines and aromatic clusters being depicted as nodes. Because there is a direct link between the full complexity of the atomistic representation and the 2D lattice, the cross-links within the coal topography can be identified. Using this structural information, a statistical simulation of the kinetically controlled liquefaction reactions of an Illinois no. 6 coal lattice was performed. The Illinois no. 6 coal lattice was composed of a distribution of singular nodes, dimers, and multi-linkage molecules that ranged from 3 nodes (or clusters) to a 16-node entity. Probabilities for the cleavage of coal interunit links were calculated on the basis of the mixture of 10 model compounds comprised of two benzene rings connected by a cross-link (e.g., diphenyl sulfide, benzyl phenyl ether, diphenyl ether, etc.) equivalent to those found in the Illinois no. 6 model. The mixture of model compounds was based on the concentrations of cross-links representative of the model (and, thus, the coal analysis). Cross-link-specific kinetic parameters (Arrhenius pre-exponential factor and activation energy) for each pyrolysis reaction were obtained from literature sources or tuned from composition data using the Kinetic Modeling Editor, a software tool used to generate and solve the balance equations in a kinetic model. At a heating rate of 5 C/min, thermolysis calculations from 360 to 490 C showed that the 2D lattice structure broke down extensively, generating mainly monomer cluster molecules. At 400 C with residence times of 600 and 1600 s, the number of cross-links between aromatic units was reduced by 40 and 60%, respectively, from the original population. At higher temperatures, the breakdown of the weaker and more kinetically favored cross-links (-CH 2S- and -OCH2-) occurred rapidly at the earlier stages of heating. The absolute probability of cleavage (Pi) of certain cross-links (-CH2CH2-, -CH2OCH2-, etc.) increased at longer residence times. Cross-links such as phenyl and doubly aromatic bound (-S- and -O-) remained essentially unreactive at the conditions explored.

AB - Various studies have used lattice structures combined with simplistic kinetics to explore the devolatilization of coal with the goal of predicting yields of gas, tar, and char. In our previous work, we have demonstrated the ability to reduce a three-dimensional (3D) large-scale (>50 000 atoms) atomistic representation of Illinois no. 6 bituminous coal to a coal-specific two-dimensional (2D) lattice with cross-links being depicted as linkage lines and aromatic clusters being depicted as nodes. Because there is a direct link between the full complexity of the atomistic representation and the 2D lattice, the cross-links within the coal topography can be identified. Using this structural information, a statistical simulation of the kinetically controlled liquefaction reactions of an Illinois no. 6 coal lattice was performed. The Illinois no. 6 coal lattice was composed of a distribution of singular nodes, dimers, and multi-linkage molecules that ranged from 3 nodes (or clusters) to a 16-node entity. Probabilities for the cleavage of coal interunit links were calculated on the basis of the mixture of 10 model compounds comprised of two benzene rings connected by a cross-link (e.g., diphenyl sulfide, benzyl phenyl ether, diphenyl ether, etc.) equivalent to those found in the Illinois no. 6 model. The mixture of model compounds was based on the concentrations of cross-links representative of the model (and, thus, the coal analysis). Cross-link-specific kinetic parameters (Arrhenius pre-exponential factor and activation energy) for each pyrolysis reaction were obtained from literature sources or tuned from composition data using the Kinetic Modeling Editor, a software tool used to generate and solve the balance equations in a kinetic model. At a heating rate of 5 C/min, thermolysis calculations from 360 to 490 C showed that the 2D lattice structure broke down extensively, generating mainly monomer cluster molecules. At 400 C with residence times of 600 and 1600 s, the number of cross-links between aromatic units was reduced by 40 and 60%, respectively, from the original population. At higher temperatures, the breakdown of the weaker and more kinetically favored cross-links (-CH 2S- and -OCH2-) occurred rapidly at the earlier stages of heating. The absolute probability of cleavage (Pi) of certain cross-links (-CH2CH2-, -CH2OCH2-, etc.) increased at longer residence times. Cross-links such as phenyl and doubly aromatic bound (-S- and -O-) remained essentially unreactive at the conditions explored.

UR - http://www.scopus.com/inward/record.url?scp=84879380896&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84879380896&partnerID=8YFLogxK

U2 - 10.1021/ef4001105

DO - 10.1021/ef4001105

M3 - Article

VL - 27

SP - 2915

EP - 2922

JO - Energy & Fuels

JF - Energy & Fuels

SN - 0887-0624

IS - 6

ER -