Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “formation 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Concentric core-shell tracks and spectroscopic properties of SrTiO 3 under intense electronic excitation

A deeper understanding of the intrinsic link between the irradiation-induced microstructures and the corresponding spectroscopic properties is becoming increasingly attractive for the research fields of materials science, physics, and information technology. In this work, the structural damage response, spectroscopic features, and associated physical mechanisms of SrTiO 3 single crystals under swift heavy ion irradiation were comparatively analyzed using a combination of experimental and theoretical approaches. Corresponding to 0.11–5.00 MeV/u ion irradiation with electronic energy loss ranging from 4.0 to 29.3 keV/nm, the inelastic thermal spike calculations combined with molecular dynamics simulations are compared with the experimental observations, revealing the track fine structures (individual spherical defects with a disordered region, and discontinuous and continuous tracks consisting of an amorphous core and a disordered outer shell) and demonstrating the dominant effects of the deposition energy and lattice temperature on track damage formation and evolution; thus, two essential thresholds for defect formation of ~0.60 eV/atom and amorphous region formation of ~1.81 eV/atom were identified to better describe the concentric core-shell track structure. The enhanced nanohillock formation in the surface region is attributed to the combined action of kinetic and potential energy depositions, and the fluence dependence in regulating the hillock dimensions is also presented. The measured refractive index profiles along different crystal axes further indicate the anisotropy of irradiation-induced lattice expansion. With increasing ion fluence and damage level, the intrinsic bandgap gradually decreased, and the concentration of additional radiative recombination centers (1.81 eV, 2.31 eV, and 2.42 eV) accordingly increased in the SrTiO 3 crystal, providing the possibility of regulating related defects to achieve tunable/selective photon emission for more critical technological applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Engineering bimetallic interfaces and revealing the mechanism for carbon dioxide electroreduction to C 3+ liquid chemicals

Reduction reaction of CO 2 (CO 2 RR) to liquid C 3+ chemicals is a net-zero-carbon process and can increase local resiliency to power outages and fuel consumption. The mechanism and the catalyst design rules of CO 2 RR-to-C 3+ are unknown. Engineering bimetallic interface (e.g., Pd/Au) to tune the intermediate adsorption is promising for promoting C 3+ formation. Our density functional theory (DFT) calculations find that *CH 2 could be the key intermediate and C 1 -CH 2 coupling could be the rate-limiting step to generate C 3+ . High CO surface coverages can promote the bimetallic interfacial sites, lower the energetics of the C 1 -CH 2 coupling step, and enhance the C 3+ formation. We further construct a volcano plot of C 1 -CH 2 kinetics as a function of the binding strength of key intermediate *CH 2 via engineering the d-band center of the interfacial site. Our findings could guide the rational design of bimetallic interfaces and their near-surface microenvironment for enhancing CO 2 RR-to-C 3+ .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Competitive Adsorption of NH 3 and H 2 O in Metal–Organic Framework Materials: MOF-74

Elucidating the interaction between coadsorbed H 2 O and NH 3 in metal–organic frameworks (MOFs) is of paramount importance to uncover mechanistic details of their competitive coadsorption behavior as well as to guide the design of new materials for enhanced NH 3 adsorption in humid environments. Nevertheless, molecular competition between NH 3 and H 2 O within the confined nanopores of MOFs was rarely explored and is poorly understood due to challenges in characterization. Here, we combine in situ infrared spectroscopy with ab initio calculations to unveil the competition of NH 3 and H 2 O for occupying active adsorption sites in the representative MOF-74 material by analyzing the kinetics and energetics of the molecular exchange process. We find that at a high NH 3 /H 2 O ratio, the incoming NH 3 is capable of displacing metal-bound H 2 O and moving it to secondary adsorption sites due to the stronger binding of NH 3 compared with H 2 O. Interestingly, the reverse process of H 2 O displacing metal-bound NH 3 is also possible upon increasing water concentration. Our calculations show that H 2 O exchanging the preabsorbed NH 3 at the metal site is driven not only by a reduced kinetic barrier but also by a favorable energetical state resulting from the formation of water clusters at metal sites and intermolecular H-bonding between the metal-coordinated H 2 O and displaced NH 3 . Our finding emphasizes that the description of molecular occupation in MOFs at equilibrium cannot simply be established by comparing molecules’ binding energies at their strongest binding sites derived by single-component measurements; rather, intermolecular interactions can greatly affect molecular distribution at equilibrium. Furthermore, we show that vibrational modes of adsorbed NH 3 are markedly perturbed upon contact with water molecules, accompanied by a large frequency shift (>30 cm –1 ) and considerable intensity decrease, which arises from the freezing of NH 3 vibrations by coadsorbed H 2 O. As a result, the mechanistic insight obtained through our study sheds light on molecular coadsorption processes in MOFs and helps to assess NH 3 removal efficiency of MOFs containing open-metal sites under realistic conditions, particularly in the presence of humidity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Potential Energy Hotspot: Effects of Impact Velocity, Defect Geometry, and Crystallographic Orientation

In energetic materials, the localization of energy into “hotspots” is known to dictate the initiation of chemical reactions and detonation. Recent all-atom simulations have shown that more energy is localized as internal potential energy (PE) than can be inferred from the kinetic energy (KE) alone. The mechanisms associated with pore collapse and hotspot formation are known to depend on pore geometry and dynamic material response such as plasticity. Therefore, we use molecular dynamics (MD) simulations to characterize shock-induced pore collapse and the subsequent formation of hotspots in 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a highly anisotropic molecular crystal, for various defect shapes, shock strengths, and crystallographic orientations. We find that the localization of energy as PE is consistently larger than the KE in cases with significant plastic deformation. Additionally, an analysis of MD trajectories reveals the underlying molecular- and crystal-level processes that govern the effect of orientation and pore shape on PE localization. We find that the regions of highest PE relate to the areas of maximum plastic deformation, while KE is maximized at the point of impact. Comparisons against octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) reveal less energy localization in TATB, which could be a contributing factor to the latter’s insensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anion Effects on the Interfacial Structure and Bulk Physical Properties in Choline-Based Hydrogen-Bonded Electrolytes

Electrolytes consisting of choline (Ch) salts with acetate (Ac) and oxalate (Ox) anions were studied in mixtures with ethylene glycol (EG), in comparison to chloride (Cl)- and bis(trifluorosulfonyl)imide (TFSI)-containing systems to understand the impact of anion size and charge density on the bulk electrolyte and interfacial properties. These represent hydrogen-bonded concentrated electrolytes where the Ch salts are the H-bond acceptors (HBAs) and EG is the H-bond donor (HBD) in 1:2 and 1:4 HBA:HBD ratios. Increased H-bonding in Ch 2 Ox with EG compacts the liquid compared to ChAc, consistent with the increased dissociation energies of the HBA:HBD clusters with Ox –2 calculated by density functional theory (DFT), thus leading to denser and more viscous liquids. Electrochemical impedance spectroscopy (EIS) and surface-enhanced Raman spectroscopy (SERS) were used to investigate interfacial charge density and the electrode surface species. Differential capacitance curves demonstrate weak dependence on potential except for ChAc:EG, where Ac – is found to adsorb on the Au surface, in particular at +0.055 V vs Ag quasi-reference while no specific adsorption was observed on glassy carbon (GC). In situ SERS measurement confirms the formation of a surface layer in the ChAc:EG systems resulting from a Au–acetate interaction. Despite the structural similarities, Ox –2 did not show specific binding or any dependence of capacitance on potential as a result of strongly solvating EG molecules via H-bonds. Furthermore, this study shows that the charge density of the HBAs impacts both bulk properties and the electrical double-layer formation in concentrated H-bonded electrolytes that is relevant to electrochemical kinetics and morphology in electrodeposition processes and energy storage applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Doping- and Strain-Dependent Electrolyte-Gate-Induced Perovskite to Brownmillerite Transformation in Epitaxial La 1– x Sr x CoO 3–δ Films

Much recent attention has focused on the voltage-driven reversible topotactic transformation between the ferromagnetic metallic perovskite (P) SrCoO 3-δ and oxygen-vacancy-ordered antiferromagnetic insulating brownmillerite (BM) SrCoO 2.5 . This is emerging as a paradigmatic example of the power of electrochemical gating (using, e.g., ionic liquids/gels), the wide modulation of electronic, magnetic, and optical properties generating clear application potential. SrCoO 3 films are challenging with respect to stability, however, and there has been little exploration of alternate compositions. Here, we present the first study of ion-gel-gating-induced P. BM transformations across almost the entire La 1-x SrxCoO 3 phase diagram (0 ≤ x ≤ 0.70), under both tensile and compressive epitaxial strain. Electronic transport, magnetometry, and operando synchrotron X-ray diffraction establish that voltage-induced P → BM transformations are possible at essentially all x, including x ≤ 0.50, where both P and BM phases are highly stable. Under small compressive strain, the transformation threshold voltage decreases from approximately +2.7 V at x = 0 to negligible at x = 0.70. Both larger compressive strain and tensile strain induce further threshold voltage lowering, particularly at low x. The P → BM threshold voltage is thus tunable, via both composition and strain. At x = 0.50, voltage-controlled ferromagnetism, transport, and optical transmittance are then demonstrated, achieving Curie temperature and resistivity modulations of similar to 220 K and at least 5 orders of magnitude, respectively, and enabling estimation of the voltage-dependent Co valence. The results are analyzed in the context of doping- and strain-dependent oxygen vacancy formation energies and diffusion coefficients, establishing that it is thermodynamic factors, not kinetics, that underpin the decrease in the threshold voltage with x, that is, with increasing formal Co valence. These findings substantially advance the practical and mechanistic understanding of this voltage-driven transformation, with fundamental and technological implications.

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

Methyl Termination of p-Type Silicon Enables Selective Photoelectrochemical CO 2 Reduction by a Molecular Ruthenium Catalyst

Methyl-terminated p-type silicon photoelectrodes selectively drive CO 2 reduction by a homogeneous [Ru(tpy)(Mebim-py)(NCCH 3 )] 2+ catalyst (tpy = 2,2′:6′,2″-terpyridine, Mebim-py = 1-methylbenzimidazol-2-ylidene-3-(2′-pyridine)). A 460 mV photovoltage is quantified for the photoelectrode. Under 1 sun illumination, this system achieves a Faradaic efficiency of 87% for CO at −1.7 V vs Fc +/0 , matching reports of the same catalyst at metallic electrodes operating at −2.1 V. When 5% water is introduced, the CH 3 -terminated Si photoelectrode remains stable, selectivity for CO is retained, and current density increases. Methyl termination suppresses the competitive hydrogen evolution observed for H-terminated Si photoelectrodes, which under the same conditions produce ca. 60% CO and 8% H 2 and have unstable performance. Furthermore, these results establish that a semiconductor photoelectrode can power a molecular CO 2 reduction catalyst without hydrogen evolution by the photoelectrode itself. Methyl termination of p-Si allows CO 2 reduction to kinetically outcompete proton reduction, revealing an important design principle for selective fuel formation.

Catalysts↗
DOE OSTI2025

Fast-Charging Li-Ion Battery Enabled by an Acetonitrile-Based Electrolyte

Fast charging remains a critical challenge for current-generation lithium-ion batteries (LIBs), particularly in electric vehicle applications. Here, in this study, we present a highly conductive electrolyte formulation based on a ternary solvent system consisting of acetonitrile (AN), fluoroethylene carbonate (FEC), and ethylene carbonate (EC), combined with a tailored additive, lithium difluoro­(oxalato)­borate (LiDFOB). This electrolyte demonstrates significantly enhanced ionic conductivity and a higher Li + transference number, enabling accelerated Li + transport kinetics. The synergistic effect of the solvents and the additive promotes the formation of a robust, low-resistance, inorganic-rich solid-electrolyte-interphase (SEI) that effectively passivates the graphite surface and suppresses AN decomposition. As a result, the electrolyte substantially reduces internal cell resistance and overpotential, both of which are critical for reliable fast charging. These findings highlight the essential role of rational electrolyte design in addressing the limitations of fast-charging LIBs.

Wu, Dezhen [Argonne National Laboratory (ANL), Arg↗
DOE OSTI2023

An investigation of shock formation vs shock mitigation of colliding plasma jets

Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reionization effective likelihood from Planck 2018 data

We release relike (reionization effective likelihood), a fast and accurate effective likelihood code based on the latest Planck 2018 data that allows one to constrain any model for reionization between 6 < z < 30 using five constraints from the CMB reionization principal components (PC). We tested the code on two example models which showed excellent agreement with sampling the exact Planck likelihoods using either a simple Gaussian PC likelihood or its full kernel density estimate. This code enables a fast and consistent means for combining Planck constraints with other reionization data sets, such as kinetic Sunyaev-Zeldovich effects, line-intensity mapping, luminosity function, star formation history, quasar spectra, etc., where the redshift dependence of the ionization history is important. Since the PC technique tests any reionization history in the given range, we also derive model-independent constraints for the total Thomson optical depth τ PC = $0.0619$$^{+0.0056}_{–0.0068}$ and its 15 ≤ z ≤ 30 high redshift component τ PC (15,30) < 0.020 (95% C.L.). Furthermore, the upper limits on the high-redshift optical depth is a factor of ~3 larger than those reported in the Planck 2018 cosmological parameter paper using the FlexKnot method and we validate our results with a direct analysis of a two-step model which permits this small high-z component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Self-Consistent Simulation of the Excitation of Compressional Alfvén Eigenmodes and Runaway Electron Diffusion in Tokamak Disruptions

Alfvénic modes in the current quench (CQ) stage of the tokamak disruption have been observed in experiments. In DIII-D the excitation of these modes is associated with the presence of high-energy runaway electrons (REs), and a strong mode excitation is often associated with the failure of RE plateau formation. In this work we present results of self-consistent kinetic-MHD simulations of RE-driven compressional Alfvén eigenmodes (CAEs) in DIII-D disruption scenarios, providing an explanation of the CQ modes. Simulation results reveal that high energy trapped REs can have resonance with the Alfvén mode through their toroidal precession motion, and the resonance frequency is proportional to the energy of REs. The mode frequencies and their relationship with the RE energy are consistent with experimental observations. Further, the perturbed magnetic fields from the modes can lead to spatial diffusion of REs including the nonresonant passing ones, thus providing the theoretical basis for a potential approach for RE mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Situating the Vector Density Approach Among Contemporary Continuum Theories of Dislocation Dynamics

Abstract For the past century, dislocations have been understood to be the carriers of plastic deformation in crystalline solids. However, their collective behavior is still poorly understood. Progress in understanding the collective behavior of dislocations has primarily come in one of two modes: the simulation of systems of interacting discrete dislocations and the treatment of density measures of varying complexity that are considered as continuum fields. A summary of contemporary models of continuum dislocation dynamics is presented. These include, in order of complexity, the two-dimensional statistical theory of dislocations, the field dislocation mechanics treating the total Kröner–Nye tensor, vector density approaches that treat geometrically necessary dislocations on each slip system of a crystal, and high-order theories that examine the effect of dislocation curvature and distribution over orientation. Each of theories contain common themes, including statistical closure of the kinetic dislocation transport equations and treatment of dislocation reactions such as junction formation. An emphasis is placed on how these common themes rely on closure relations obtained by analysis of discrete dislocation dynamics experiments. The outlook of these various continuum theories of dislocation motion is then discussed.

Engineering↗
DOE OSTI2020

Outflow from Outer-arm Starburst in a Grazing Collision between Galaxies

Gemini NIFS K-band spectra and Atacama Large Millimeter/submillimeter Array {sup 12}CO J=1→0, HCO{sup +}, and 100 GHz continuum observations are used to study a bright starburst clump on an outer arm of the interacting galaxy NGC 2207. This clump emits 23% of the total 24 μm flux of the galaxy pair and has an optically opaque dust cone extending out of its 170 pc core. The measured CO accounts for the dark cone extinction if almost all the gas and dust there are in front of the star clusters. An associated approaching CO outflow has v {sub z} ~ 16 km s{sup -1}, an estimated molecular mass 8 × 10{sup 6} M {sub ⊙}, and rises to heights ~0.9 kpc. A receding CO outflow on the far side with v {sub z} ~ 28 km s{sup -1} is less extensive. The observed star formation in the core over 10 Myr can supply the dark cone kinetic energy of roughly 2 × 10{sup 52} erg via supernovae and stellar winds. Other signs of intense activity are a variable radio continuum, suggesting an embedded supernova or other outburst; X-ray emission possibly from an X-ray binary or intermediate-mass black hole, depending on the extinction; and Brγ and He i lines with 82 km s{sup -1} line widths and fluxes consistent with excitation by embedded O-type stars. According to previous models, the retrograde encounter suffered by NGC 2207 caused the loss of angular momentum. This compressed its outer disk. We suggest that the resulting inward crashing stream hit a massive H i clump on the preexisting spiral arm and triggered the observed starburst.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Cool flame chemistry of diesel surrogate compounds: n -Decane, 2-methylnonane, 2,7-dimethyloctane, and n -butylcyclohexane

Elucidating the formation of combustion intermediates is crucial to validate reaction pathways, develop reaction mechanisms and examine kinetic modeling predictions. While high-temperature pyrolysis and oxidation intermediates of alkanes have been thoroughly studied, comprehensive analysis of cool flame intermediates from alkane autoxidation is lacking and challenging due to the complexity of intermediate species produced. Here, jet-stirred reactor autoxidation of four C 10 alkanes: n -decane, 2-methylnonane, 2,7-dimethyloctane, and n -butylcyclohexane, as model compounds of diesel fuel, was investigated from 500 to 630 K using synchrotron vacuum ultraviolet photoionization molecular beam mass spectrometry (SVUV-PIMS). Around 100 intermediates were detected for each fuel. The classes of molecular structures present during the autoxidation of the representative paraffinic functional groups in transport fuels, i.e., n -alkanes, branched alkanes, and cycloalkanes were established and were found to be similar from the oxidation of various alkanes. A theoretical approach was applied to estimate the photoionization cross sections of the intermediates with the same carbon skeleton as the reactants, e.g., alkene, alkenyl keto, cyclic ether, dione, keto-hydroperoxide, diketo-hydroperoxide, and keto-dihydroperoxide. These species are indicators of the first, second, and third O 2 addition reactions for the four C 10 hydrocarbons, as well as bimolecular reactions involving keto-hydroperoxides. Chemical kinetic models for the oxidation of these four fuels were examined by comparison against mole fraction of the reactants and final products obtained in additional experiments using gas chromatography analysis, as well as the detailed species pool and mole fractions of aforementioned seven types of intermediates measured by SVUV-PIMS. This works reveals that the models in the literature need to be improved, not only the prediction of the fuel reactivity and final products, but also the reaction network to predict the formation of many previous undetected intermediates.

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

A pseudo-kinetic model to simulate phase changes in gas hydrate bearing sediments

Herein, modeling of the phase transitions anticipated in gas hydrate bearing sediments (GHBS) is critical for a proper understanding of time-dependent changes in states and volumes (e.g. the production of methane from this type of soils). We propose a new pseudo-kinetic approach to simulate the typical phase changes anticipated in GHBS, using published experimental results involving gas hydrate dissociation that are the basis of a widely used kinetic model. The proposed pseudo-kinetic model is formulated in the pressure-temperature (P-T) plane and assumes a rate of gas hydrate dissociation (or formation) proportional to the distance between the current state and the phase boundary. The model consists of only one parameter and is simple to implement in numerical simulators. A similar concept is used to model ice formation/thawing phenomena, but based on the ice/liquid-water phase boundary. We implemented the pseudo-kinetic model in a fully coupled thermo-hydro-chemo-mechanical (THCM) finite element code and validated it against experimental results performed on the dissociation of synthetic gas hydrate. We also evaluated the pseudo-kinetic model using synthetic cases covering several scenarios associated with gas hydrate formation/dissociation and ice formation/thawing. The model successfully reproduced the gas production test from a natural GHBS core from Korea (scaled gas venting experiment over 14 h), and also the formation of gas hydrate and ice in permafrost in Alaska (over 2 × 10 6 years). -The analyses show the versatility of the proposed pseudo-kinetic approach by applying it to model the different types of phase transitions typically encounter in GHBS. The simple formulation, easy implementation in numerical simulator, and reduced number of parameters (only one per phase change) make this model an attractive option for simulating phase transformations in problems involving GHBS.

58 GEOSCIENCES↗
DOE OSTI2025

Bulk Anion Recognition Kinetically Holds Back Interfacial Adsorption

The competition between bulk and interfacial phenomena underlies many key processes in complex chemical phenomena and transport. While competitive processes are often framed in a thermodynamic context, opportunities to leverage transient species found away from equilibrium can provide a kinetic handle to achieve unconventional reaction outcomes. In this work, we outfit an iminoguanidinium headgroup capable of selective SO 4 2– complexation with alkyl tails of varying complexity to probe competitive bulk and interfacial reaction pathways and tune kinetic pathways for selective chemical separations. Using sum frequency generation (SFG) vibrational spectroscopy we unexpectedly find that adsorption of ligands to the air–aqueous interface was dramatically slowed down for species with increasingly hydrophobic tails. Underlying this phenomenon, we show that the formation of bulk colloidal species with differing propensities for SO 4 2– inhibited surface adsorption via a kinetic bottleneck in the exchange of molecular extractants with colloidal aggregates. This kinetic effect could open up avenues to access unconventional selectivity via complexation of strongly coordinating species in the bulk phase, allowing for more weakly coordinating species to transport via interfacial mechanisms. Furthermore, this work broadly probes nonequilibrium phenomena in chemical separations that arise through unexpected interfacial events that are neglected in traditional equilibrium descriptions.

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

Promotional Effect of ZnO and ZrO 2 in K-Doped Fe Catalysts for CO 2 Hydrogenation to Light Olefins

Promoters play a critical role in tuning the activity and selectivity of Fe catalysts in CO 2 hydrogenation to produce light olefins, which are key building blocks in the petrochemical industry. Herein, by a combined experimental and theoretical approach, we show that high and stable performance of Fe catalysts could be achieved by taking advantage of the promotional effect of both Zn and Zr. Structural characterization indicates that ZnO could improve the dispersion and reducibility of Fe oxides and facilitate the formation of active Fe carbide species, whereas ZrO 2 could stabilize the structure and catalytic performance, especially the selectivity of hydrocarbon products. In situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) experiments suggest that carbonate, bicarbonate, formate, and methoxy are essential intermediates in the CO 2 hydrogenation to hydrocarbon products, including paraffins and olefins. The conversion kinetics of each intermediate species are dependent on the type of promoters as well as the phase structure of the active Fe species. DFT calculations revealed a strong correlation between the formation energy of surface oxygen vacancies and that of Fe carbide species in promoted Fe oxides, in accordance with experimental results. Moreover, the calculated energy profiles of CO 2 hydrogenation over different catalysts indicate that Zn could promote the activation of CO 2 and its transformation to the oxygenate intermediates, while Zr could facilitate the conversion of oxygenates to hydrocarbon precursors. Finally, the discrepancies in the evolution trend of various intermediate species on promoted Fe catalysts in the transient DRIFTS experiments can be rationalized by the differences in energy barriers of elementary or rate limiting steps.

CO2 conversion↗
DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗