Search NASA⌕ Search

SEARCH · Search NASA

Results for “Decomposition kinetics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Elucidating the differences in oxidation of high-performance α- and β- diisobutylene biofuels via Synchrotron photoionization mass spectrometry

Abstract Biofuels are a promising ecologically viable and renewable alternative to petroleum fuels, with the potential to reduce net greenhouse gas emissions. However, biomass sourced fuels are often produced as blends of hydrocarbons and their oxygenates. Such blending complicates the implementation of these fuels in combustion applications. Variations in a biofuel’s composition will dictate combustion properties such as auto ignition temperature, reaction delay time, and reaction pathways. A handful of novel drop-in replacement biofuels for conventional transportation fuels have recently been down selected from a list of over 10,000 potential candidates as part of the U.S. Department of Energy’s (DOE) Co-Optimization of Fuels and Engines (Co-Optima) initiative. Diisobutylene (DIB) is one such high-performing hydrocarbon which can readily be produced from the dehydration and dimerization of isobutanol, produced from the fermentation of biomass-derived sugars. The two most common isomers realized, from this process, are 2,4,4-trimethyl-1-pentene (α-DIB) and 2,4,4-trimethyl-2-pentene (β-DIB). Due to a difference in olefinic bond location, the α- and β- isomer exhibit dramatically different ignition temperatures at constant pressure and equivalence ratio. This may be attributed to different fragmentation pathways enabled by allylic versus vinylic carbons. For optimal implementation of these biofuel candidates, explicit identification of the intermediates formed during the combustion of each of the isomers is needed. To investigate the combustion pathways of these molecules, tunable vacuum ultraviolet (VUV) light (in the range 8.1–11.0 eV) available at the Lawrence Berkeley National Laboratory’s Advanced Light Source (ALS) has been used in conjunction with a jet stirred reactor (JSR) and time-of-flight mass spectrometry to probe intermediates formed. Relative intensity curves for intermediate mass fragments produced during this process were obtained. Several important unique intermediates were identified at the lowest observable combustion temperature with static pressure of 93,325 Pa and for 1.5 s residence time. As this relatively short residence time is just after ignition, this study is targeted at the fuels’ ignition events. Ignition characteristics for both isomers were found to be strongly dependent on the kinetics of C 4 and C 7 fragment production and decomposition, with the tert-butyl radical as a key intermediate species. However, the ignition of α-DIB exhibited larger concentrations of C 4 compounds over C 7 , while the reverse was true for β-DIB. These identified species will allow for enhanced engineering modeling of fuel blending and engine design.

09 BIOMASS FUELS↗
DOE OSTI2020

CHAPTER 9. Bijels the Easy Way

Spinodal decomposition is not the only way to make a bijel. Indeed, while spinodal decomposition produces structures with a potentially useful morphology, it can be challenging to make bijels using this method and the resulting systems can be hard to process and manipulate. Furthermore, exploiting the functional properties of the assembled particle monolayer is extremely challenging. Here, we show how the assembly of nanoparticle surfactants at the liquid–liquid interface can be used to kinetically trap liquids into a wealth of complex structures without using spinodal decomposition. We apply liquid three-dimensional printing and moulding methods, along with patterned substrates with controllable wetting properties, to build all-liquid devices with applications in chemical synthesis, separation, and purification. The functional properties of the assembled nanomaterials can be exploited to produce interfacially structured liquids that are plasmonically and magnetically responsive. Finally, we conclude by arguing that, while the field shows great promise, efforts need to be made to translate liquid bicontinuous systems out of the laboratory and into meaningful, real-world applications, as well applications in more ‘exotic’ disciplines, such as synthetic biology.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polysulfone shock compressed above the decomposition threshold: Velocimetry and modeling of two-wave structures

Polysulfone was shock loaded to pressures of 14.6–26.2 GPa in a series of gas gun-driven plate-impact experiments measuring material response with embedded electro-magnetic particle velocity gauges and optical velocimetry. The embedded electro-magnetic particle velocity gauges did not show a distinct two-wave structure but did show rounding that suggested a reaction but not a distinct separation of the reactants and product waves. In contrast, the transmission experiments fielded with optical velocimetry, with product pressures ranging from 21.0 to 26.2 GPa, showed well-defined two-wave structures due to shock-driven chemical decomposition of the polymer to products at a higher density. Distinct two-wave structures have not previously been published in literature during polymer compression; here, we observed these two-wave structures at both the polymer/lithium fluoride and polymer/polymethylmethacrylate interfaces. Hydrodynamic simulations of the experiments were performed using a newly constructed SESAME equation of state (EOS) for the polymer, a thermochemical EOS for the decomposition products, and an Arrhenius reaction rate model for the kinetics of the reaction. Finally, simulation results also demonstrated two-wave structures but were unable to quantitatively reproduce either the embedded gauge or optical velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Energizing Robust Sulfur/Lithium Electrochemistry via Nanoscale-Asymmetric-Size Synergism

Sluggish redox kinetics and dendrite growth perplex the fulfillment of efficient electrochemistry in lithium–sulfur (Li–S) batteries. The complicated sulfur phase transformation and sulfur/lithium diversity kinetics necessitate an all-inclusive approach in catalyst design. Herein, a compatible mediator with nanoscale-asymmetric-size configuration by integrating Co single atoms and defective CoTe 2–x (Co SA -CoTe 2–x @NHCF) is elaborately developed for regulating sulfur/lithium electrochemistry synchronously. Substantial electrochemistry and theoretical analyses reveal that CoTe 2–x exhibits higher catalytic activity in long-chain polysulfide transformation and Li 2 S decomposition, while monodispersed Co sites are more effective in boosting sulfur reduction kinetics to regulate Li 2 S deposition. Such cascade catalysis endows Co SA -CoTe 2–x @NHCF with the all-around service of “trapping-conversion-recuperation” for sulfur species during the whole redox reaction. Furthermore, it is demonstrated by in situ transmission electron microscopy that initially formed electronic-conductive Co and ionic-conductive Li 2 Te provide sufficient lithiophilic sites to regulate homogeneous Li plating and stripping with markedly suppressed dendrite growth. Consequently, by coupling the Co SA -CoTe 2–x @NHCF interlayer and Li@Co SA -CoTe 2–x @NHCF anode, the constructed Li–S full batteries deliver superior cycling stability and rate performance, and the flexible pouch cell exhibits stable cycling performance at 0.3 C. In conclusion, the gained insights into the synergistic effect of asymmetric-size structures pave the way for the integrated catalyst design in advanced Li–S systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The role of ternary alloying elements in eutectoid transformation of U-10Mo alloy part II. In and Ex-situ Neutron diffraction-based assessment of eutectoid phase transformation kinetics in U-9.8Mo-0.2X alloy (X = Cr, Ni or Co)

HExploring the effects of minor ternary alloying additions, typically impurity elements, on the phase stability of U-10Mo is important for preventing undesirable phase decomposition during processing or during service. This work examines the influence small ternary additions of Cr, Ni, and Co. Both in-situ and ex-situ neutron diffraction measurements made during and after high temperature (450 – 525°C) exposures were used to better define the influence of these elements on the time-temperature-transformation (TTT) behavior of U-10Mo, providing information which is complementary to electron microscopy investigations of the same alloy systems performed in the first part of this work. Minor additions of Ni and Co decrease the ?-phase stability at all temperatures investigated. Signatures of U6X (X = Ni or Co) compounds were shown to be present in amounts of up to 6 wt% in the heat treated alloys, suggesting that the initial precipitation of this phase may catalyze further ?-phase decomposition. On the other hand, the Cr containing alloys were observed to have nearly the same, and in some cases slower, phase transformation kinetics when compared to the binary U-10Mo control samples. The results of the study have enabled preliminary estimates of the TTT curves for the ternary alloys.

phase transformation, kinetics, discontinuous prec↗
DOE OSTI2023

Modeling–Experiment–Theory Analysis of Reactions Initiated from Cl + Methyl Formate

Methyl formate (MF; CH 3 OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF—CH 3 OCO and CH 2 OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl 2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl 2 , HCl, CO 2 , CH 3 , CH 3 Cl, CH 2 O, and CH 2 ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400–750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH 3 OCO (R1a) and Cl + MF → HCl + CH 2 OCHO (R1b): k 1a,theory = 6.71 × 10 –15 ·T 1.14 ·exp(—606/T) cm 3 /molecule·s; k 1b,theory = 4.67 × 10 –18 ·T 2.21 ·exp(—245/T) cm 3 /molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH 3 OCO and CH 2 OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH 3 + CO 2 (R3) and CH 2 O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k 3,theory (P = ∞) = 2.94 × 10 9 ·T 1.21 ·exp(—6209/T) s –1 and k 5,theory (P = ∞) = 8.45 × 10 8 ·T 1.39 ·exp(—15132/T) s –1 over T = 300–1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH 2 OCHO + Cl, CH 2 OCHO + Cl 2 , and CH 3 + Cl 2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

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

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

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

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ignition delay time and speciation of dibutyl ether at high pressures

Dibutyl ether (DBE) is a promising biofuel due to its high cetane number (~ 100) and high volumetric energy density (31.6 MJ/L). It could either be used directly in compression ignition engines or blended with other conventional or renewable fuels. Oxidation and pyrolysis kinetics of DBE are not well known, particularly at high pressures. In this work, we have experimentally investigated the chemical kinetics of DBE in three domains: (a) ignition delay time measurements in a rapid compression machine over T = 550–650 K, P = 10, 20, 40 bar, φ = 0.5, 1; (b) ignition delay time measurements in a shock tube over T = 900–1300 K, P = 20, 40 bar, φ = 0.5, 1; (c) laser-based carbon monoxide speciation measurements in a shock tube during DBE pyrolysis and oxidation over T = 1100–1400 K, P = 20 bar. Pressure time-histories measured in RCM experiments exhibited unique 3-stage and 4-stage ignition behavior predominantly at fuel-lean conditions. Experimental data were compared with the predictions of two recent chemical kinetic models of DBE. Sensitivity analyses were carried out to identify key reactions which may have caused the discrepancy between experiments and simulations. We found that the rate of decomposition of DBE may need to be revisited to improve the oxidative and pyrolytic predictions of DBE kinetic model.

09 BIOMASS FUELS↗
DOE OSTI2021

Automated theoretical chemical kinetics: Predicting the kinetics for the initial stages of pyrolysis

Large scale implementation of high level computational theoretical chemical kinetics offers the prospect for dramatically improving the fidelity of combustion chemical modeling. To facilitate such efforts, we have developed a suite of codes, collectively referred to as AutoMech, that allow for the automatic prediction of the kinetics for large sets of reactions via ab initio transition-state-theory based master-equation calculations. The primary input is simply the mechanism, a dictionary relating chemically identifiable species descriptors (e.g., SMILES or InChIs) to species labels in the mechanism, and a specification of the electronic structure and transition state theory models to be implemented. Here we illustrate the current utility of AutoMech through a study of the initial stages of pyrolysis for 3 sets of fuels: sequences of alkanes, alcohols, and aldehydes. For simplicity, the analysis focuses on abstractions from the fuel by H, CH 3 , and OH, and the decomposition of the resulting radicals. Altogether, there are a total of 166 input channels in these sets (more than 363 forward reactions when expanded to the full set of elementary reactions). The code successfully produces high quality rate estimates (with apparent uncertainties less than a factor of two in limited comparisons with experiment) for > 95% of these. For the radical decomposition reactions, the analysis includes predictions for the pressure dependence of the kinetics. This wide-ranging exploration illustrates (i) the effect of different levels of prediction on the expected accuracy, (ii) the branching between abstractions at different sites for different abstractors, (iii) the dependence of the rates on the chemical structure, and (iv) the variation in radical stabilities across chemical families. These results, as well as the demonstrated feasibility of the methodology, should find further utility in the development of accurate rate expressions for arbitrary fuels. (c) 2020 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Elliott, Sarah N.↗
DOE OSTI2025

Multi-Scale Modeling for Plasma-Enhanced Ammonia Decomposition over Carbides and Nitrides

Ammonia is a carbon-free hydrogen carrier, but its decomposition typically requires high temperatures over costly Ru-based catalysts due to the large barrier for N≡N bond formation. We develop a multiscale framework combining density functional theory, zero-dimensional plasma kinetics, and microkinetic modeling to elucidate how non-thermal plasma (NTP) enables low-temperature NH 3 decomposition over Co-based carbides and nitrides, benchmarked against Ru and Co. Under thermal conditions, all catalysts are limited by N≡N bond formation, with Co 3 C(001) most active owing to its negatively charged surface, strong N* binding, and low activation barriers of N≡N bond formation. Plasma-induced vibrational excitation of NH 3 and its reactive radicals promotes a radical-driven •NH 2 –N* coupling pathway that dominates on Co 3 C(001) and Co 3 N(001), shifting the rate-limiting step to NH 3 (v1) dissociation, increasing turnover frequencies by up to 6 orders of magnitude, and reducing the temperature needed to reach a turnover frequency of 5 s –1 from >680 °C (Ru and Co under thermal condition) to 267 °C (Co 3 C) and 415 °C (Co 3 N). These results identify Co-based carbides and nitrides as promising plasma-active catalysts for energy-efficient hydrogen production from ammonia.

ammonia decomposition↗
DOE OSTI2021

Understanding the Stability of Etched or Platinized p-GaInP Photocathodes for Solar-Driven H 2 Evolution

The long-term stability in acidic or alkaline aqueous electrolytes of p-Ga 0.52 In 0.48 P photocathodes, with a band gap of ~1.8 eV, for the solar-driven hydrogen-evolution reaction (HER) has been evaluated from a thermodynamic, kinetic, and mechanistic perspective. At either pH 0 or pH 14, etched p-GaInP electrodes corroded cathodically under illumination and formed metallic In 0 on the photoelectrode surface. In contrast, under the same conditions, electrodeposition of Pt facilitated the HER kinetics and stabilized p-GaInP/Pt photoelectrodes against such cathodic decomposition. When held at 0 V versus the reversible hydrogen electrode, p-GaInP/Pt electrodes in either pH = 0 or pH = 14 exhibited stable current densities (J) of ~-9 mA cm -2 for hundreds of hours under simulated 1 sun illumination. During the stability tests, the current density-potential (J-E) characteristics of the p-GaInP/Pt photoelectrodes degraded due to pH-dependent changes in the surface chemistry of the photocathode. This work provides a fundamental understanding of the stability and corrosion mechanisms of p-GaInP photocathodes that constitute a promising top light absorber for tandem solar-fuel generators.

14 SOLAR ENERGY↗
DOE OSTI2022

Cascade Reaction of Ethanol to Butadiene over Multifunctional Silica-Supported Ag and ZrO 2 Catalysts

Although butadiene is currently a by-product of naphtha cracking, interest in producing butadiene from bio-based ethanol has increased because of the lower environmental impact of the ethanol to butadiene reaction. Furthermore, this work explores a multifunctional catalyst system composed of silica-supported Ag and ZrO 2 used for the cascade reaction of ethanol to butadiene at 573 K. The Ag and ZrO 2 components were synthesized on separate support particles enabling characterization of each component without interference from the other. High selectivity to butadiene (65%) at high ethanol conversion (75%) was achieved with an appropriate ratio of Ag and ZrO 2 in the reactor. Silver catalyzed the initial dehydrogenation of ethanol to acetaldehyde while ZrO 2 catalyzed the C-C coupling and subsequent dehydration reactions. The silica-supported ZrO 2 exhibited superior selectivity relative to bulk ZrO 2 in the Ag-promoted ethanol to butadiene reaction. Results from Zr K-edge X-ray absorption spectroscopy and UV-Vis spectroscopy showed that ZrO 2 was highly dispersed on the silica support over a range of loadings. Infrared spectroscopy of adsorbed pyridine, CO, and CO 2 , and kinetics of probe reactions 1-butene double bond isomerization, 2-propanol decomposition, and ethanol hydrogenation of acetone were used to compare the acid-base nature and chemical reactivity of silica-supported ZrO 2 to bulk ZrO 2 .

09 BIOMASS FUELS↗
DOE OSTI2021

Improving Photostability of Cesium-Doped Formamidinium Lead Triiodide Perovskite

Using time-resolved in situ X-ray photoelectron spectroscopy, we identify and suppress rapid degradation mechanisms for cesium-stabilized formamidinium lead iodide perovskite materials used in state-of-the-art photovoltaics. Accelerated degradation under high light intensity and heating reveals a photocatalytic reaction pathway responsible for rapid decomposition in iodide-rich compositions and a slower pathway for more stoichiometric samples. Using Avrami–Erofe’ev kinetic analysis, we find that the fast process is consistent with a 2D crystallization and growth mechanism fueled by excess halide salt at grain boundaries and surfaces. Moreover, the rate of decomposition varies dramatically with the wavelength of light used to illuminate the thin films. Furthermore, our results reveal the photodegradation mechanisms of PbI 2 and excess iodide and provide a path to increasing perovskite stability under photoexcitation.

14 SOLAR ENERGY↗
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Pyrolysis of bio-derived dioxolane fuels: A ReaxFF molecular dynamics study

Alkyl-substituted 1,3-dioxolanes, including 4,5-dimethyl-2-pentan-3-yl-1,3-dioxolane (Fuel 1), 4,5-dimethyl-2-pentyl-1,3-dioxolane (Fuel 2), and 2-(heptan-3-yl)-4,5-dimethyl-1,3-dioxolane (Fuel 3), have been recently suggested as potential biodiesels. In this paper, we investigate the initial pyrolysis of the alkyl-substituted 1,3-dioxolanes at high temperatures using ReaxFF molecular dynamics (MD) simulations. We analyze the decomposition rate, reaction mechanism, and product distribution in the pyrolysis of the three alkyl 1,3-dioxlanes. The three fuels primarily decompose to 4,5-dimethyl-1,3-dioxolane radical and hydrocarbons derived from the alkyl side-chains. The further decomposition of 4,5-dimethyl-1,3-dioxolane radical primarily leads to 2-C 4 H 8 and CO 2 within a few decomposition steps. The hydrocarbon product distribution is significantly affected by the molecular structure of the alkyl side-chain, which would have a strong influence on the sooting tendency of these fuels. The ReaxFF simulations predict that the order of sooting tendency would be Fuel 3 > Fuel 1 > Fuel 2, which agrees with the measured sooting tendency trend. Based on the pyrolysis mechanism identified by ReaxFF, we propose a new alkyl dioxolane, 4-hexyl-5-methyl-1,3-dioxolane (Fuel 4), which might produce even less soot, by modifying the molecular structure of Fuel 2. Our ReaxFF simulation shows that Fuel 4 produce much less C 4 H 8 , an effective non-aromatic soot precursor, than Fuel 2. Moreover, more carbon atoms are bonded to each oxygen atom in Fuel 4 than Fuel 2, which would help reduce soot yield by removing more carbon atoms from the soot-producing pool of species. The major decomposition pathways identified in this work can be used to develop chemical kinetic models for 1,3-dioxolane based compounds, as biodiesel components, applicable to combustion engine simulations. We also demonstrate that the chemical kinetic insight offered by ReaxFF simulations can be used to design new fuel molecules with more desired properties.

09 BIOMASS FUELS↗
DOE OSTI2020

In Situ Investigation of Chemomechanical Effects in Thiophosphate Solid Electrolytes

Solid-state batteries can suffer from catastrophic failure at high current densities due to solid electrolyte fracture, interface decomposition, or lithium filament growth. Failure is linked to chemomechanical material transformations that can manifest during electrochemical cycling. We systematically investigate how solid electrolyte microstructure and interfacial decomposition (e.g., interphase) affect failure mechanisms in lithium thiophosphates (Li3PS4, LPS) electrolytes. Kinetically metastable interphases are engineered with iodine doping, and microstructural control is achieved using milling and annealing processing techniques. In situ transmission electron microscopy reveals iodine diffusion to the interphase, and upon electrochemical cycling, pores are formed in the interphase region. In situ synchrotron tomography reveals that interphase pore formation drives edge fracture events, which are the origin of through-plane fracture failure. Fractures in thiophosphate electrolytes actively grow toward regions of higher porosity and are affected by heterogeneity in microstructure (e.g., porosity factor). This report provides fundamental design guidelines for high-performance solid-state batteries.

25 ENERGY STORAGE↗