Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Free electron model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Bespoke Liquid/Liquid Interfaces (Final Technical Report)

The goal of DE-SC0001815 was to advance the basic science of liquid:liquid interface formation, to develop a deeper understanding of the mechanisms of phase separation and the essential relationships between solution composition, organization and dynamics that underlie the kinetic regime of solvent extraction. This included learning how interfacial organization and dynamics alters the properties of the primary coordination sphere of ions and the free energy of transport of ions complexes across a phase boundary. We relied primarily upon classical molecular dynamics studies to determine the equilibrium ensembles of these complex systems, but also utilized ab-initio MD and cluster-based density functional theory (DFT) calculations when more detailed investigation of the electronic structure was needed. We continued development of graph-theory based analyses to elucidate hierarchical correlations and expanded into geometric topology methods to quantify the collectively organized structures that can organize at a liquid/liquid interface during solute transport. One of the main conclusions was from the observation of two distinct mechanisms for solute transport - those that derive from amplifications of interfacial heterogeneity and surface roughness, and those wherein surface roughness has been dampened and instead collectively organized macrostructures work to bring solutes into the organic phase. It was our aim to create a concrete chemical model of the underlying driving forces behind interfacial primary and secondary structure formation and to map out the energetic features of solute transport so that tailored liquid/liquid can be developed that have characteristic kinetic features associated with mass transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Single-shot in-line x-ray phase-contrast imaging of void-shockwave interactions in fusion energy materials

Recent breakthroughs in nuclear fusion, specifically the report of reactions exceeding scientific breakeven at the National Ignition Facility (NIF), highlight the potential of inertial fusion energy (IFE) as a sustainable and virtually limitless energy source. However, further progress in IFE requires characterization of defects in ablator materials and how they affect fuel capsule compression. Voids within the ablator can degrade energy yield, but their impact on the density distribution has primarily been studied through simulations, with limited high-resolution experimental validation. To address this, we used the x-ray free-electron laser (XFEL) at the matter in extreme conditions (MECs) instrument at the Linac coherent light source (LCLS) to capture 2D x-ray phase-contrast (XPC) images of a void-bearing sample with a composition similar to inertial confinement fusion (ICF) ablators. By driving a compressive shockwave through the sample using MEC's long-pulse laser system, we analyzed how voids influence shockwave propagation and density distribution during compression. To quantify this impact, we extracted phase information using two phase retrieval algorithms. First, we applied the contrast transfer function (CTF) method, paired with Tikhonov regularization and a fast optimization approach to generate an initial phase estimate. We then refined the result using a projected gradient descent (PGD) method that works directly with the sample's refractive index. Comparing these results with radiation adaptive grid Eulerian (xRAGE) radiation hydrodynamic simulations enables identification of model validation needs or improvements. By calculating phase maps in situ, it becomes possible to reconstruct areal density maps, improving understanding of laser-capsule interactions and advancing IFE research.

Hodge, D. S. [Colorado State Univ., Fort Collins, ↗
NASA NTRSDate not supplied

Wake Velocimetry of a Sphere-Cone Model in a Mach 10 Air Freestream using FLEET

Velocity measurements in the wake behind a sphere-cone-shaped vehicle wind tunnel model in a Mach 10 hypersonic air flow using femtosecond laser electronic excitation tagging (FLEET) are reported in this work. The FLEET measurements were performed at 1 kHz using a femtosecond laser centered at 800 nm and an ultrahigh-speed intensified camera system. The results reported here are for a single freestream condition corresponding to approximately Re ∞ /L = 3.6∙10 6 m -1 . The performance of the velocimetry instrument was evaluated in the quiescent test section at conditions relevant to the wake measurements. For velocity measurements in the wake, the FLEET measurement volume was translated to several radial locations from the model centerline in the vertical direction and at a streamwise location corresponding to approximately 12 mm from the payload. Reverse flow with mean velocities ranging from -12 m/s to -48 m/s were observed in the recirculating region of the wake, whereas hypersonic velocities in the range of 1079 m/s to 1183 m/s were observed in the free shear layer. A velocity uncertainty methodology is also outlined and applied for the reported measurements. The velocity data reported in this work is expected to play a significant role in the validation of computational codes modeling the hard-to-predict separated hypersonic wake.

FLEET↗
NASA NTRSDate not supplied

Wake Velocimetry of a Sphere-Cone Model in a Mach 10 Air Freestream using FLEET

Velocity measurements in the wake behind a sphere-cone-shaped vehicle wind tunnel model in a Mach 10 hypersonic air flow using femtosecond laser electronic excitation tagging (FLEET) are reported in this work. The FLEET measurements were performed at 1 kHz using a femtosecond laser centered at 800 nm and an ultrahigh-speed intensified camera system. The results reported here are for a single freestream condition corresponding to approximately Re ∞ /L = 3.6∙10 6 m -1 . The performance of the velocimetry instrument was evaluated in the quiescent test section at conditions relevant to the wake measurements. For velocity measurements in the wake, the FLEET measurement volume was translated to several radial locations from the model centerline in the vertical direction and at a streamwise location corresponding to approximately 12 mm from the payload. Reverse flow with mean velocities ranging from -12 m/s to -48 m/s were observed in the recirculating region of the wake, whereas hypersonic velocities in the range of 1079 m/s to 1183 m/s were observed in the free shear layer. A velocity uncertainty methodology is also outlined and applied for the reported measurements. The velocity data reported in this work is expected to play a significant role in the validation of computational codes modeling the hard-to-predict separated hypersonic wake.

FLEET↗
NASA NTRS2015

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗
NASA NTRS2014

Magnetic Field Amplification in the Thin X-Ray Rims of SN 1006

Several young supernova remnants, including SN 1006, emit synchrotron X-rays in narrow filaments, hereafter thin rims, along their periphery. The widths of these rims imply 50–100 μG fields in the region immediately behind the shock, far larger than expected for the interstellar medium compressed by unmodified shocks, assuming electron radiative losses limit rim widths. However, magnetic field damping could also produce thin rims. Here we review the literature on rim width calculations, summarizing the case for magnetic field amplification. We extend these calculations to include an arbitrary power-law dependence of the diffusion coefficient on energy, D ∝ Eμ. Loss-limited rim widths should shrink with increasing photon energy, while magnetic-damping models predict widths almost independent of photon energy. We use these results to analyze Chandra observations of SN 1006, in particular the southwest limb. We parameterize the FWHM in terms of energy as FWHM ∝ EmE γ. Filament widths in SN 1006 decrease with energy; mE ∼ −0.3 to −0.8, implying magnetic field amplification by factors of 10–50, above the factor of four expected in strong unmodified shocks. For SN 1006, the rapid shrinkage rules out magnetic damping models. It also favors short mean free paths (small diffusion coefficients) and strong dependence of D on energy (μ ≥ 1).

MAGNETIC FIELD AMPLIFICATION IN THE THIN X-RAY RIM↗
DOE OSTI2024

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Efficient prediction of attosecond two-colour pulses from an X-ray free-electron laser with machine learning

Abstract X-ray free-electron lasers are sources of coherent, high-intensity X-rays with numerous applications in ultra-fast measurements and dynamic structural imaging. Due to the stochastic nature of the self-amplified spontaneous emission process and the difficulty in controlling injection of electrons, output pulses exhibit significant noise and limited temporal coherence. Standard measurement techniques used for characterizing two-coloured X-ray pulses are challenging, as they are either invasive or diagnostically expensive. In this work, we employ machine learning methods such as neural networks and decision trees to predict the central photon energies of pairs of attosecond fundamental and second harmonic pulses using parameters that are easily recorded at the high-repetition rate of a single shot. Using real experimental data, we apply a detailed feature analysis on the input parameters while optimizing the training time of the machine learning methods. Our predictive models are able to make predictions of central photon energy for one of the pulses without measuring the other pulse, thereby leveraging the use of the spectrometer without having to extend its detection window. We anticipate applications in X-ray spectroscopy using XFELs, such as in time-resolved X-ray absorption and photoemission spectroscopy, where improved measurement of input spectra will lead to better experimental outcomes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

The X-ray surface brightness of Kepler's supernova remnant

The first X-ray images of Kepler's supernova remnant (SN Ophiuchi 1604) are presented, and consequences for SNR models are discussed. Observations made with the Einstein Observatory Imaging Proportional Counter and High Resolution Imager show the remnant to be circular, with a strong shell brighter in the north than in the south. A flux of 1.2 x 10 to the -10th ergs/sq cm per sec was measured in the 0.15-4.5 keV region, which corresponds to an X-ray luminosity of 1.0 x 10 to the 36th ergs/sec at a distance of 5 kpc and an interstellar medium density of 2.8 x 10 to the 21st/sq cm. The X-ray observations do not allow the determination of whether the SNR is in the adiabatic or free expansion phase, but in either case it is shown that the mean ISM density must be greater than about 0.1/cu cm. In addition, the density of the X-ray emitting gas must be high, and its electron temperature must be fairly low. The high ISM densities derived for Kepler's SNR and other SNRs thus suggest an atypical ISM, possibly influenced by mass lost from the pre-supernova star.

White, R. L.↗
NASA NTRS2009

Tin Whisker Electrical Short Circuit Characteristics

Existing risk simulations make the assumption that when a free tin whisker has bridged two adjacent exposed electrical conductors, the result is an electrical short circuit. This conservative assumption is made because shorting is a random event that has an unknown probability associated with it. Note however that due to contact resistance electrical shorts may not occur at lower voltage levels. In our first article we developed an empirical probability model for tin whisker shorting. In this paper, we develop a more comprehensive empirical model using a refined experiment with a larger sample size, in which we studied the effect of varying voltage on the breakdown of the contact resistance which leads to a short circuit. From the resulting data we estimated the probability distribution of an electrical short, as a function of voltage. In addition, the unexpected polycrystalline structure seen in the focused ion beam (FIB) cross section in the first experiment was confirmed in this experiment using transmission electron microscopy (TEM). The FIB was also used to cross section two card guides to facilitate the measurement of the grain size of each card guide's tin plating to determine its finish.

Courey, Karim J.↗
DOE OSTI2021

Clinical Natural Language Processing for Radiation Oncology: A Review and Practical Primer

Natural language processing (NLP), which aims to convert human language into expressions that can be analyzed by computers, is one of the most rapidly developing and widely used technologies in the field of artificial intelligence. Natural language processing algorithms convert unstructured free text data into structured data that can be extracted and analyzed at scale. In medicine, this unlocking of the rich, expressive data within clinical free text in electronic medical records will help untap the full potential of big data for research and clinical purposes. Recent major NLP algorithmic advances have significantly improved the performance of these algorithms, leading to a surge in academic and industry interest in developing tools to automate information extraction and phenotyping from clinical texts. Thus, these technologies are poised to transform medical research and alter clinical practices in the future. Radiation oncology stands to benefit from NLP algorithms if they are appropriately developed and deployed, as they may enable advances such as automated inclusion of radiation therapy details into cancer registries, discovery of novel insights about cancer care, and improved patient data curation and presentation at the point of care. However, challenges remain before the full value of NLP is realized, such as the plethora of jargon specific to radiation oncology, nonstandard nomenclature, a lack of publicly available labeled data for model development, and interoperability limitations between radiation oncology data silos. Successful development and implementation of high quality and high value NLP models for radiation oncology will require close collaboration between computer scientists and the radiation oncology community. Here, we present a primer on artificial intelligence algorithms in general and NLP algorithms in particular; provide guidance on how to assess the performance of such algorithms; review prior research on NLP algorithms for oncology; and describe future avenues for NLP in radiation oncology research and clinics.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An approach for the measurement of the bulk temperature of single crystal diamond using an X-ray Free Electron Laser

We present a method to determine the bulk temperature of a single crystal diamond sample at an X-Ray Free Electron Laser using inelastic X-ray scattering. The experiment was performed at the High Energy Density (HED) instrument at the European XFEL GmbH, Germany. The technique, based on inelastic X-ray scattering and the principle of detailed balance, was demonstrated to give accurate temperature measurements, within 8% for both room temperature diamond and heated diamond to 500 K. Here, the temperature was increased in a controlled way using a resistive heater to test theoretical predictions of the scaling of the signal with temperature. The method was tested by validating the energy of the phonon modes with previous measurements made at room temperature using inelastic X-ray scattering and neutron scattering techniques. This technique could be used to determine the bulk temperature in transient systems with a temporal resolution of 50 fs and for which accurate measurements of thermodynamic properties are vital to build accurate equation of state models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

An approach for the measurement of the bulk temperature of single crystal diamond using an X-ray free electron laser

We present a method to determine the bulk temperature of a single crystal diamond sample at an X-Ray free electron laser using inelastic X-ray scattering. The experiment was performed at the high energy density instrument at the European XFEL GmbH, Germany. The technique, based on inelastic X-ray scattering and the principle of detailed balance, was demonstrated to give accurate temperature measurements, within 8% for both room temperature diamond and heated diamond to 500 K. Here, the temperature was increased in a controlled way using a resistive heater to test theoretical predictions of the scaling of the signal with temperature. The method was tested by validating the energy of the phonon modes with previous measurements made at room temperature using inelastic X-ray scattering and neutron scattering techniques. This technique could be used to determine the bulk temperature in transient systems with a temporal resolution of 50 fs and for which accurate measurements of thermodynamic properties are vital to build accurate equation of state and transport models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Application and performance of the novel Elcomat 5000 autocollimator in slope-measuring profilometry of beamline optics

Autocollimators are used in profilometers for the precise form measurement of beam-shaping optics for synchrotron light sources and x-ray free-electron lasers. This application requires using an aperture stop, typically 2.5 mm in diameter, which limits the footprint of the autocollimator’s measuring beam on the optics to achieve a sufficiently high lateral resolution. Central to the development of profilometry has been the availability of commercial autocollimators, such as the Elcomat 3000 from Möller–Wedel Optical. Although autocollimators are generally not designed for use with small apertures, this device has proven capable and has become the de facto standard. Now that it has been replaced by the Elcomat 5000, we evaluate the performance of several autocollimators of this type at small apertures, including the characterization of instrument transfer functions with a chirped height profile. We demonstrate that the Elcomat 5000 is capable of repeatable angle measurements with a standard deviation of 0.014 arcsec at an aperture diameter of 1.6 mm, whereas the Elcomat 3000 achieves 0.03 arcsec at 2.5 mm, comparing models with modified reticle designs. Here, we also investigate the influence of optical aberrations of the autocollimator’s objective and their changes with the path length to the surface under test. We characterize the sensitivity of the angle measurement to changes in environmental parameters, in particular, barometric pressure. These efforts are aimed at approaching fundamental limits in the characterization of the shape of optical surfaces with autocollimator-based profilometers so that these are ready for characterizing the next generation of synchrotron and XFEL beamline optics.

angle metrology↗
DOE OSTI2024

Effects of mosaic crystal instrument functions on x-ray Thomson scattering diagnostics

Mosaic crystals, with their high integrated reflectivities, are widely employed in spectrometers used to diagnose high energy density systems. X-ray Thomson scattering (XRTS) has emerged as a powerful diagnostic tool of these systems, providing in principle direct access to important properties such as the temperature via detailed balance. However, the measured XRTS spectrum is broadened by the spectrometer instrument function (IF), and without careful consideration of the IF one risks misdiagnosing system conditions. Here, we consider in detail the IF of 40 and 100 μm mosaic Highly Annealed Pyrolytic Graphite crystals, and how the broadening varies across the spectrometer in an energy range of 6.7–8.6 keV. Notably, we find a strong asymmetry in the shape of the IF toward higher energies. As an example, we consider the effect of the asymmetry in the IF on the temperature inferred via XRTS for simulated 80 eV CH plasmas and find that the temperature can be overestimated if an approximate symmetric IF is used. We, therefore, expect a detailed consideration of the full IF will have an important impact on system properties inferred via XRTS in both forward modeling and model-free approaches.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Comparative determinant quantum Monte Carlo study of the acoustic and optical variants of the Su-Schrieffer-Heeger model

We compare the acoustic Su–Schrieffer–Heeger (SSH) model with two of its optical variants where the phonons are defined on either on the sites or bonds of the system. First, we discuss how to make fair comparisons between these models in any dimension by ensuring their dimensionless coupling λ and relevant phonon energies are the same. We then use determinant quantum Monte Carlo to perform nonperturbative and sign-problem-free simulations of all three models on one-dimensional chains at and away from half-filling. By comparing the results obtained from each model, we demonstrate that the optical and acoustic models produce near identical results within error bars for suitably chosen phonon energies and λ at half-filling. In contrast, the bond model has quantitatively different behavior due to its coupling to the q = 0 phonon mode. These differences also manifest in the total length of the chain, which shrinks for the bond model but not for the acoustic and optical models when λ ≠ 0. Furthermore, our results have important implications for quantum Monte Carlo modeling of SSH-like interactions, where these models are sometimes regarded as being interchangeable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Large-scale atomistic model construction of subbituminous and bituminous coals for solvent extraction simulations with reactive molecular dynamics

Large-scale atomistic models for complex polycyclic aromatic hydrocarbon systems help understand the chemical properties and behaviors of complex feedstocks such as coal or petroleum. However, the development and utilization of large-scale models remain limited due to the difficulty in achieving the varied structural characteristics necessary to capture stochastic nature of these feedstocks. Here we demonstrate a systematic workflow to construct stochastic molecular systems from a broad analytical suite: high-resolution transmission electron microscopy (HRTEM), carbon-13 nuclear magnetic resonance spectroscopy ( 13 C NMR), laser desorption ionization mass spectroscopy (LDI-MS), and elemental analysis. We present a model construction and analysis utility of a new Python-based module. We selected one subbituminous and three high-volatile bituminous coals to construct large-scale models (~40,000 atoms). The constructed models were utilized to examine the affinity for solvent extraction (naphthalene or tetralin) and the effect of structural properties (e.g., aromatic cluster size, functional groups, and cross-linking) in reactive molecular dynamics simulations. Complex chemical reactions were monitored with bond order transitions, intermediates formation, and mass distributions. Reactive molecular dynamics simulations suggest a plausible chemical extraction process and products for the complex fossil feedstocks. The results indicated that radical formations with bond breaking of bridging oxygens and carbons were required at high temperatures to facilitate hydrogeneration and extraction of gas molecules from radical-free molecules. We observed that aliphatic chains of tetralin were easily decomposed and combined with radicals to form small size of molecules with aryl bonding, mainly increasing molecules in the 500–1000 Da, while naphthalene had little impact on chemical extraction process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗