An official website of the United States government
A series of reactive molecular dynamics simulations is used to study the internal structure of incipient soot particles obtained from acetylene pyrolysis. The simulations were performed using ReaxFF potential at four different temperatures. The resulting soot particles are cataloged and analyzed to obtain statistics of their mass, volume, density, C/H ratio, number of cyclic structures, and other features. A total of 3324 incipient soot particles were analyzed in this study. Based on their structural characteristics, the incipient soot particles are classified into two classes, referred to as type 1 and type 2 incipient soot particles in this work. The radial distribution of density, cyclic (5-, 6-, or 7-member rings) structures, and C/H ratio inside the particles revealed a clear difference in the internal structure between type 1 and type 2 particles. These classes were further found to be well represented by the size of the particles, with smaller particles in type 1 and larger particles in type 2. The radial distributions of ring structures, density, and C/H ratio indicated the presence of a dense core region in type 2 particles. In contrast, no clear evidence of the presence of a core was found in type 1 particles. In type 2 incipient soot particles, the boundary between the core and shell was found to be around 50%–60% of the particle radius of gyration.
more »
« less
This content will become publicly available on October 1, 2025
Physical, chemical and morphological evolution of incipient soot obtained from molecular dynamics simulation of acetylene pyrolysis
Incipient soot particles obtained from a series of reactive molecular dynamics simulations were studied to understand the evolution of physical, chemical, and morphological properties of incipient soot. Reactive molecular dynamics simulations of acetylene pyrolysis were performed using ReaxFF potential at 1350, 1500, 1650, and 1800 K. A total of 3324 incipient soot particles were extracted from the simulations at various stages of development. Features such as the number of carbon and hydrogen atoms, number of ring structures, mass, C/H ratio, radius of gyration, surface area, volume, atomic fractal dimension, and density were calculated for each particle. The calculated values of density and C/H ratio matched well with experimental values reported in the literature. Based on the calculated features, the particles were classified in two types: type 1 and type 2 particles. It was found that type 1 particles show significant morphological evolution while type 2 particles undergo chemical restructuring without any significant morphological change. The particle volume was found to be well-correlated with the number of carbon atoms in both type 1 and type 2 particle, whereas surface area was found to be correlated with the number of carbon atoms only for type 1 particles. A correlation matrix comparing the level of correlation between any two features for both type 1 and type 2 particle was created. Finally, based on the calculated statistics, a set of correlations among various physical and morphological parameters of incipient soot was proposed.
more »
« less
- Award ID(s):
- 2144290
- PAR ID:
- 10584347
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Fuel
- Volume:
- 373
- Issue:
- C
- ISSN:
- 0016-2361
- Page Range / eLocation ID:
- 132197
- Subject(s) / Keyword(s):
- Soot Molecular Dynamics Radius of gyration C/H Ratio Surface Area
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Soot or black carbons are combustion-generated carbonaceous nanoparticles formed during the incomplete combustion of hydrocarbon fuels. The complexity of hydrocarbon systems often makes it difficult to investigate the fundamentals of soot formation experimentally. To address this, this study uses reactive molecular dynamics simulations with reactive force field (ReaxFF) potentials. The current work focuses on the formation and evolution of soot during acetylene pyrolysis. The analysis provides insights into the physicochemical aspects of soot formation and the maturation of incipient soot particles. In this work, we focus on the evolution and interdependence of features such as the number of carbon atoms, number of aromatic rings, mass, C/H ratio, the radius of gyration, atomic fractal dimension, surface area, volume, and density. Based on the physicochemical features, two distinct classes of nascent soot can be observed. These are termed type-1 and type-2 particles. The type-1 particles show significant morphological evolution, while the type-2 particles show chemical restructuring without significantly changing the morphology. Qualitative correlations of various degrees are also observed between some of these morphological features.more » « less
-
Abstract. Soot inception by acetylene pyrolysis at 1350–1800 K is investigated using reactive molecular dynamics. The composition and chemical structure of soot precursor molecules formed during inception are elucidated. During soot inception, increasing the process temperature leads to faster depletion of C2H2 molecules and faster formation of C2H3, C2H4, C2H6, CH4, and C2 with the concurrent appearance of H2 molecules. Small molecules consisting of 1–5 C atoms (C1–C5) are formed due to reactive collisions and grow further to larger hydrocarbon compounds consisting of 6–10 C atoms. At initial stages of inception, prior to the formation of incipient soot, three-member rings are formed, which are associated with the formation of compounds with fewer than 10 C atoms. Once incipient soot is formed, the number of C1–C10 compounds and the number of three-member rings drop, while the number of five- and six-member rings increases, indicating that the formation of larger rings is associated with the growth of soot clusters. The chemical structure of soot precursor molecules obtained by bond order analysis reveals that molecules with up to 10 C atoms are either linear or branched aliphatic compounds or may contain three-member rings fused with aliphatic components. Molecules with more than 10 C atoms often exhibit structures composed of five- or six-member C rings, decorated by aliphatic components. The identification of molecular precursors contributing to soot inception provides crucial insights into soot formation mechanisms, pinpointing potential pathways of soot formation during combustion.more » « less
-
Black carbon data including number concentrations of refractory black carbon particles, and number concentration of scattering particles from the University of California - Riverside Single Particle Soot Photometer (SP2) that was deployed on the NSF NCAR C-130 aircraft for the CAESAR campaign.more » « less
-
The rheology of suspensions of rings (tori) rotating in an unbounded low Reynolds number simple shear flow is calculated using numerical simulations at dilute particle number densities ( n ≪ 1 ). Suspensions of non-Brownian rings are studied by computing pair interactions that include hydrodynamic interactions modeled using slender body theory and particle collisions modeled using a short-range repulsive force. Particle contact and hydrodynamic interactions were found to have comparable influences on the steady-state Jeffery orbit distribution. The average tilt of the ring away from the flow-vorticity plane increased during pairwise interactions compared to the tilt associated with Jeffery rotation and the steady-state orbit distribution. Particle stresses associated with the increased tilt during the interaction were found to be comparable to the stresses induced directly by particle contact forces and the hydrodynamic velocity disturbances of other particles. The hydrodynamic diffusivity coefficients in the gradient and vorticity directions were also obtained and were found to be two orders of magnitude larger than the corresponding values in fiber suspensions at the same particle concentrations. Rotary Brownian dynamics simulations of isolated Brownian rings were used to understand the shear rate dependence of suspension rheology. The orbit distribution observed in the regime of weak Brownian motion, P e ≫ ϕ T − 3, was surprisingly similar to that obtained from pairwise interaction calculations of non-Brownian rings. Here, the Peclet number P e is the ratio of the shear rate and the rotary diffusivity of the particle and ϕ T is the effective inverse-aspect ratio of the particle (approximately equal to 2 π times the inverse of its non-dimensional Jeffery time period). Thus, the rheology results obtained from pairwise interactions should retain accuracy even for weakly Brownian rings ( n ≪ 1 and ϕ T − 3 ≪ P e ).more » « less
