Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic strength”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗
DOE OSTI2021

Comparison of Conventional and Constrained Variational Methods for Computing Large-Scale Budgets and Forcing Fields

Analyses of atmospheric heat and moisture budgets serve as an effective tool to study convective characteristics over a region and to provide large-scale forcing fields for various modeling applications. This paper examines two popular methods for computing large-scale atmospheric budgets: the conventional budget method (CBM) using objectively gridded analyses based primarily on radiosonde data and the constrained variational analysis (CVA) approach which supplements vertical profiles of atmospheric fields with measurements at the top of the atmosphere and at the surface to conserve mass, water, energy, and momentum. Successful budget computations are dependent on accurate sampling and analyses of the thermodynamic state of the atmosphere and the divergence field associated with convection and the large-scale circulation that influences it. Utilizing analyses generated from data taken during Dynamics of the Madden-Julian Oscillation (DYNAMO) field campaign conducted over the central Indian Ocean from October to December 2011, we evaluate the merits of these budget approaches and examine their limitations. While many of the shortcomings of the CBM, in particular effects of sampling errors in sounding data, are effectively minimized with CVA, accurate large-scale diagnostics in CVA are dependent on reliable background fields and rainfall constraints. Finally, for the DYNAMO analyses examined, the operational model fields used as the CVA background state provided wind fields that accurately resolved the vertical structure of convection in the vicinity of Gan Island. However, biases in the model thermodynamic fields were somewhat amplified in CVA resulting in a convective environment much weaker than observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Deciphering the non-linear impact of Al on chemical durability of silicate glass

The role of Al in aluminosilicate glasses remains somewhat a mystery: at low concentrations, it increases the resistance to hydrolysis of the glass, whereas at high concentrations an opposite effect is observed. To understand the origin of the phenomenon on a fundamental atomistic scale, we performed 577 MD simulations and applied potential mean force (PMF) calculations to estimate the activation barriers for hydrolysis and to statistically correlate them with local structural features of the glass. Models of pure silicate and aluminosilicate glasses are constructed and investigated. PMF simulation results are further validated by the experimental measurements and revealed that Al is very easy to dissociate, but it also increases the glass chemical durability through significantly increasing both the strength of Si and network connectivity of the glass. In contrast, at high Al concentration, preferential dissolution of Al weakens the silicate network, which it supposes to strengthen, and so the glass resistance becomes poor. Through PMF calculations, we evaluated the activation barriers for dissociating bonds around Al as 0.49 eV, which is less than a half of the energy to dissociate bonds around Si in pure silicate (1.22 eV) and around Si in aluminosilicate glass (1.34 eV), all these energy differences being statistically significant. Molecular structural level investigation revealed that Si with Al as a second neighbor in the glass network has a significantly higher activation energy for dissociation than Si in pure silicate glass. Furthermore, the proposed approach opens the way to the development of quantitative predictive models of glass durability.

36 MATERIALS SCIENCE↗
DOE OSTI2021

System Configuration Evaluation for Process Settling of Hanford Waste Solid Particles

Direct Feed High-Level Waste (DFHLW) is a potential flowsheet operations approach to initiating high-level waste (HLW) vitrification prior to completion of the Hanford Waste Treatment and Immobilization Plant (WTP) Pretreatment Facility. A settle/decant process has been proposed to concentrate solids prior to delivery to the WTP HLW Facility during DFHLW operations, wherein the solids in a settled layer would be remixed with the supernatant liquid remaining after decanting operations to provide the feed at required solids concentrations. Settling would be used in lieu of purpose-built filtration or other solids separation equipment. Pacific Northwest National Laboratory (PNNL) is providing baseline technical support to the Washington River Protection Solutions (WRPS) Flowsheet Integration group. To support planning for DFHLW, WRPS previously requested that PNNL evaluate the current data set available to predict the time needed for HLW solids to settle and the solids concentration and strength of that settled layer, to identify gaps in the understanding and predictive capability of HLW solids waste settling times, and to provide scoping estimates of the potential settling times. Eight technical gaps were identified for predicting settling times and characteristics of the formed sediment layers. In addition to the data gaps, an overarching observation was made that there is significant variation in behavior of settling rate and settled layer data. The settling time required to concentrate solids via a settle/decant process was determined from the limited data to have a difference of potentially more than a factor of 5,000 in the estimated settling times, varying from 0.2 to 1,060 days for example depending on process vessel depth and final sediment solids concentration. In contrast, successful processes of liquid-forward output streams resulting from in-tank settling and decanting forward liquid have been reported for operations conducted at the Hanford Site. The purpose of this current report is to further support DFHLW planning by evaluating double-shell tank (DST) and alternate vessel equipment and operational configurations to enable optimization of the settle/decant process to concentrate solids. Hanford waste processing behavior specific to liquid feed availability following a slurry transfer in a DST is summarized, including process stream characteristics and process equipment configurations. The performance of DST process equipment configurations is evaluated for possible improvements using computational fluid dynamics (CFD) and simple analytical models. Potential new vessel design(s) specific to enabling effective settle/decant processes, and cursory summary of other separate and inline solids separations processes, are also provided. The CFD results indicated that improvement in outflow solids concentration was promoted by a reduction in the slurry flow rate, angling the distributor nozzles downward, and lifting the transfer pump. The solid-liquid analysis evaluating particle trajectory confirmed that the potential for particle ingestion (in the transfer pump) was decreased with increased radial separation between the inlet and outlet (transfer pump inlet), decreased inlet flow, and decreased liquid density and viscosity for a neutrally buoyant inlet flow. An assessment was also made of the potential for inflow configuration changes to result in the discrete mounding or piling of solids within the tank. Based on the characterization of the settled waste to date, HLW sediments will be unlikely to sustain a substantial angle of repose to facilitate significant variations in the elevation of the settled solids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Dislocation cross-slip in precipitation hardened Mg–Nd alloys

An integrated experimental characterization and molecular dynamics (MD) simulation approach was used to explore the dislocation-precipitate interactions in a dilute Magnesium–Neodymium (Mg–Nd) precipitation hardened alloy. In situ indentation in a transmission electron microscope (TEM) and postmortem TEM characterization of deformed samples and MD simulations revealed that basal < a > type dislocations interacted differently with β 1 (Mg 3 Nd) and β"(Mg 3-7 Nd) precipitates. For β" precipitates, the basal dislocations directly shear the precipitates. For β 1 precipitates, such shearing becomes much more difficult because it requires the creation of antiphase boundaries in the ordered lattice of β 1 precipitates. Screw dislocations were observed to cross-slip from basal to the prismatic plane, which could be parallel to the broad facet of β 1 precipitates. It is postulated that double cross-slip (basal to prismatic to basal) via the Hirsch mechanism may enable screw dislocations to overcome the β 1 precipitates. MD simulations also revealed that an edge dislocation is unable to bypass the precipitate via the generation of screw dislocation segments plus the double cross-slip mechanism. The edge dislocation can bow around the precipitate, and the segment that adopts a screw character can cross-slip to the prismatic plane but with increasing applied stress, cross-slips back to the original basal plane to continue glide via precipitate shearing. Finally, the implication of glide dislocation - β 1 -precipitate interaction mechanisms on the strength and ductility of β 1 precipitate dominant microstructures is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Biomimetic Light Harvesting Complexes Based on Self-Assembled Dye-DNA Nanostructures (Final Technical Report)

A fundamental limitation in the development of artificial molecular light harvesting structures is that generally if one places pigment molecules close enough to achieve very high energy transfer efficiency and absorbance cross section, the dyes interact in such a way that self-quenching occurs. This limits both the optical cross section and the overall efficiency of light capture in such systems. Nature overcomes this problem and achieves efficient energy transfer without self-quenching by controlling the environment, vibration characteristics, orientation and dynamics of the pigments involved in light harvesting. Natural light harvesting systems have evolved highly intricate chromophore architectures of protein-bound organic molecules that absorb and funnel solar energy with high efficiency and fast transfer rates using elegant excitonic circuitries. Nature achieves this precise control over the properties of pigment complexes by using specific proteins as templates to guide the formation of very well-defined pigment aggregates, groups of pigment molecules that are in close enough contact to undergo excitonic interactions resulting in new electronic properties that were not present in the monomer. In synthetic systems, pigment aggregates are well known, and J-aggregates and H-aggregates have been long studied for their unique properties. In Nature, the chlorosomes of green bacteria form what appear to be classical J-aggregates and the light harvesting complexes of purple nonsulfur bacteria form ring structures that have many of the properties of J-aggregates such as a red-shifted absorbance spectrum and excited states that are delocalized over multiple pigment molecules. These photosynthetic systems leverage the spectral and enhanced transport properties of J-like aggregates to achieve fast exciton transport within the antenna and thus high light-harvesting efficiency. The close interaction of the pigments results in substantial electronic coupling that lifts the degeneracy of excitonic states and therefore modifies the absorption and emission spectra of the monomer. In the case of J-aggregation, lower energy excitonic states acquire significant oscillator strength, resulting in a characteristic bathochromic shift of the absorption and fluorescence and often superradiant delocalized exciton states.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Unified Universal Control and Coordination of Inverter-Based Resources, and Validation for a PV + Battery Hybrid Plant

As renewable energy deployment grows, hybrid power plants (HPPs) combining photovoltaic (PV) and battery systems must evolve to offer both energy and grid stability services. These systems typically include a mix of grid-following (GFL) and grid-forming (GFM) inverters, presenting unique coordination and control challenges. This Department of Energy–funded project developed and validated a Unified Universal Control and Coordination (UUCC) framework for such PV + battery hybrid plants, enabling seamless and stable operation, including ultrafast black start, autonomous synchronization, and robust frequency and voltage regulation, under different grid conditions. The project significantly advanced the understanding of inverter-based resource (IBR) control by developing and validating three complementary system-level approaches for hybrid GFL/GFM operation: 1. A combined Virtual Resistance (VR)-based GFL and Virtual Oscillator Control (VOC)-based GFM method, where each inverter type is governed by a specialized control strategy. Together, these achieve stable, fast-response coordination, eliminating inrush current and enabling smooth black start and grid synchronization across a wide range of grid strengths. 2. A Deadbeat-based UUCC strategy, which uses discrete-time, switching-cycle-level control for both GFL and GFM inverters. This approach replaces traditional PI/PLL control with a control parameter-free, high-bandwidth framework that supports stable LVRT and instantaneous synchronization under all conditions. 3. A benchmark comparison with Siemens’ commercial GFM microgrid controller, which provided a fast baseline platform. The commercial approach decoupled v & f control was implemented on a commercial microgrid controller.The baseline commercial benchmark helped highlight superior transient response and black start performance offered by the deadbeat and VOC approaches. These technical contributions offer substantial improvements over conventional inverter control schemes, which often rely on slow phase-locked loop (PLL)-based synchronization, require careful control parameters tuning, and prone to unstable in weak grids with GFL inverters and in stiff grid with GFM inverters therefore challenging for hybrid GFL+GFM under all grid conditions. The deadbeat-based UUCC framework enables simpler, faster, and more robust operation of hybrid IBR systems using wide-bandgap (WBG) devices such as SiC power semiconductors. The rapid expansion of hybrid distributed energy resources (DERs), including residential and commercial PV-BESS installations such as Tesla Powerwall, PV with vehicle-to-grid (V2G) capability, and other integrated configurations, presents complex operational challenges for medium-voltage radial distribution feeders. These networks are subject to frequent disturbances such as faults, switching operations, rapid reclosing sequences, and feeder reconfigurations, all of which introduce dynamic stress on IBRs. In addition, planned feeder segmentation and deliberate islanding for resilience will require DERs that can autonomously perform blackstart, establish voltage and frequency references, and resynchronize with the main grid. The advanced deadbeat-based UUCC control and blackstart functionalities developed in this project directly address these requirements, enabling decentralized and autonomous operation of inverter-dominated DERs in distribution systems under a wide range of fault and reconfiguration scenarios. From a public benefit perspective, these innovations enable more reliable and cost-effective integration of renewable energy into distribution networks. The ability to autonomously black start and stabilize grids under varying grid conditions support accelerates recovery from outages and support decentralized resilient energy systems. By reducing system complexity and improving performance, this project lays critical groundwork for future inverter-dominated power grids that are clean, reliable, and accessible to all.

14 SOLAR ENERGY↗
DOE OSTI2022

The role of information structures in game-theoretic multi-agent learning

Multi-agent learning (MAL) studies how agents learn to behave optimally and adaptively from their experience when interacting with other agents in dynamic environments. The outcome of a MAL process is jointly determined by all agents’ decision-making. Hence, each agent needs to think strategically about others’ sequential moves, when planning future actions. The strategic interactions among agents makes MAL go beyond the direct extension of single-agent learning to multiple agents. With the strategic thinking, each agent aims to build a subjective model of others decision-making using its observations. Such modeling is directly influenced by agents’ perception during the learning process, which is called the information structure of the agent’s learning. As it determines the input to MAL processes, information structures play a significant role in the learning mechanisms of the agents. This review creates a taxonomy of MAL and establishes a unified and systematic way to understand MAL from the perspective of information structures. Here we define three fundamental components of MAL: the information structure (i.e., what the agent can observe), the belief generation (i.e., how the agent forms a belief about others based on the observations), as well as the policy generation (i.e., how the agent generates its policy based on its belief). In addition, this taxonomy enables the classification of a wide range of state-of-the-art algorithms into four categories based on the belief-generation mechanisms of the opponents, including stationary, conjectured, calibrated, and sophisticated opponents. We introduce Value of Information (VoI) as a metric to quantify the impact of different information structures on MAL. Finally, we discuss the strengths and limitations of algorithms from different categories and point to promising avenues of future research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

High-Pressure Synthesis of Bulk Cobalt Cementite, Co 3 C

Transition metal carbides find use in a wide range of advanced high-resilience applications including high-strength steels, heat shields, and deep-earth drills. However, carbides of the mid-to-late transition metals remain difficult to isolate and characterize on account of their metastability, which precludes the preparation of high-quality bulk single crystal samples using traditional solid-state methods. Herein, we report a combined computational and experimental survey of the cobalt-carbon binary system under high pressures and demonstrate that pressure offers a route toward the bulk synthesis of the metastable cementite-type cobalt carbide, Co 3 C, which under ambient conditions can only be prepared in low-dimensional thin film or nanoparticle forms. First-principles calculations reveal two competitive low-energy stoichiometric phases under ambient pressures - Pnnm-Co 2 C (Fe 2 C-type) and Pnma-Co 3 C (Fe 3 Ctype) - consistent with the known low-dimensional phases that have been studied for their promising magnetic properties. However, the calculated formation enthalpy of Pnma-Co 3 C decreases steadily with the applied pressure, while that of Pnnm-Co 2 C increases. Here, we pursue these results using high-pressure laser-heated synthesis methods coupled with in situ X-ray diffraction and observe the formation of Pnma-Co 3 C above 4.8 GPa. We determine the experimental bulk modulus of Co 3 C to be K o = 237 GPa (K p = 4.0). First-principles calculations of the phonon modes in Co 3 C reveal dynamical instabilities at ambient pressure that are absent under compression. These results offer a promising new route for the synthesis of rare-earth-free magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. Here, we apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.

14 SOLAR ENERGY↗
DOE OSTI2023

Measurement of ion temperature and toroidal flow during magnetic reconnection with a large guide field

Here, ion temperature and toroidal flow along the guide field direction are measured using a new ion tomographic diagnostic on the Magnetic Reconnection eXperiment (MRX) during magnetic reconnection with a guide field strength of about 1.4 and 2.1 times the strength of the reconnecting component. Strong toroidal flows, beyond what has been measured in anti-parallel and lower guide field conditions on MRX, are observed. Sustained ion heating with no discernible structure within the measurement region is also observed. Probe measurements including Langmuir and Mach probe measurements are made to support the tomographic inversion of line-integrated measurements, as well as to provide local measurements of plasma parameters. Measurements of toroidal velocity and ion temperature are supported with time series data. Energy flow into and out of the X-line region is estimated using a guiding center framework and presented in the Appendix of this manuscript, suggesting an outsized role played by parallel electric field in energizing ions. The guiding center approximation is not well satisfied in the region of interest; however, the estimates provide a springboard for future, further experimentation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Critical role of solvent-modulated hydrogen-binding strength in the catalytic hydrogenation of benzaldehyde on palladium

Protic and non-protic solvents influence the rate of benzaldehyde hydrogenation on Pd by one order of magnitude in the following sequence of increasing reactivity, dioxane < tetrahydrofuran < water < methanol. Despite these differences, the reaction pathway does not change; the majority of turnovers occurs via a Langmuir-Hinshelwood mechanism by stepwise addition of sorbed H to sorbed benzaldehyde, first to the carbonyl O and then to the C atom of the formyl group, forming benzyl alcohol. An analysis of the solvation energies shows that both ground and transition states are destabilized by the solvents compared to the same reaction at the gas-solid interface. The destabilization extent of the reacting organic substrates in both states are similar and, therefore, compensate each other, making the net kinetic effects inconsequential. In consequence, the marked reactivity differences arise solely from the differences in the surface stabilization of adsorbed hydrogen adatoms moderated by the four solvents-more weakly bound hydrogen adatoms, because of their solvation by solvent, exhibit higher rates of their hydrogenation reaction in kinetically relevant attack on the C atom.

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

Theory of the phase transition in random unitary circuits with measurements

Here, we present a theory of the entanglement transition tuned by measurement strength in qudit chains evolved by random unitary circuits and subject to either weak or random projective measurements. The transition can be understood as a nonanalytic change in the amount of information extracted by the measurements about the initial state of the system, quantified by the Fisher information. To compute the von Neumann entanglement entropy $\textit{S}$ and the Fisher information $\mathcal{F}$, we apply a replica method based on a sequence of quantities $\tilde{S}^{(n)}$ and $\mathcal{F}^{(n)}$ that depend on the $\textit{n}$th moments of density matrices and reduce to $\textit{S}$ and $\mathcal{F}$ in the limit $\textit{n}$ → 1. These quantities with $\textit{n}$ ≥ 2 are mapped to free energies of a classical spin model with $\textit{n}$! internal states in two dimensions with specific boundary conditions. In particular, $\tilde{S}^{(n)}$ is the excess free energy of a domain wall terminating on the top boundary, and $\mathcal{F}^{(n)}$ is related to the magnetization on the bottom boundary. Phase transitions occur as the spin models undergo ordering transitions in the bulk. Taking the limit of large local Hilbert space dimension $\textit{q}$ followed by the replica limit $\textit{n}$ → 1, we obtain the critical measurement probability $p_c$ = 1/2 and identify the transition as a bond percolation in the 2D square lattice in this limit. Finally, we show there is no phase transition if the measurements are allowed in an arbitrary nonlocal basis, thereby highlighting the relation between the phase transition and information scrambling. We establish an explicit connection between the entanglement phase transition and the purification dynamics of a mixed state evolution and discuss implications of our results to experimental observations of the transition and simulability of quantum dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bimetallic-Derived Catalytic Structures for CO 2 -Assisted Ethane Activation

In recent years, the simultaneous upgrading of CO 2 and ethane has emerged as a promising approach for generating valuable gaseous (CO, H 2 , and ethylene) and liquid (aromatics and C 3 oxygenates) chemicals from the greenhouse gas CO 2 and large-reserved shale gas. The key challenges for controlling product selectivity lie in the selective C-H and C-C bond cleavage of ethane with the assistance of CO 2 . Bimetallic-derived catalysts likely undergo alloying or oxygen-induced segregation under reaction conditions, thus providing diverse types of interfacial sites, e.g., metal/support (M/M'O x ) interface and metal oxide/metal (M'O x /M) inverse interface, that are beneficial for selective CO 2 -assisted ethane upgrading. The alloying extent can be initially predicted by cohesive energy and atomic radius (or Wigner-Seitz radius), while the preference for segregation to form the on-top suboxide can be approximated using the work function, electronegativity, and binding strength of adsorbed oxygen. Furthermore, bimetallic-derived catalysts are typically supported on high surface area oxides. Modifying the reducibility and acidity/basicity of the oxide supports and introducing surface defects facilitate CO 2 activation and oxygen supplies for ethane activation. Using in situ synchrotron characterization and density functional theory (DFT) calculations, we found that the electronic properties of oxygen species influence the selective cleavage of C-H/C-C bonds in ethane, with electron-deficient oxygen over the metal (or alloy) surface promoting nonselective bond scission to produce syngas and electron-enriched oxygen over the metal oxide/metal interface enhancing selective C-H scission to yield ethylene. Here, we further demonstrate that the preferred structures of the catalyst surfaces, either alloy surfaces or metal oxide/metal inverse interfaces, can be controlled through the appropriate choice of metal combinations and their atomic ratios. Through a comprehensive comparison of experimental results and DFT calculations, the selectivity of C-C/C-H bond scission is correlated with the thermodynamically favorable bimetallic-derived structures (i.e., alloy surfaces or metal oxide/metal inverse interfaces) under reaction conditions over a wide range of bimetallic catalysts. These findings not only offer structural and mechanistic insights into bimetallic-derived catalysts but also provide design principles for selective catalysts for CO 2 -assisted activation of ethane and other light alkanes. This Account concludes by discussing challenges and opportunities in designing advanced bimetallic-derived catalysts, incorporating new reaction chemistries for other products, employing precise synthesis strategies for well-defined structures with optimized site densities, and leveraging time/spatial/energy-resolved in situ spectroscopy/scattering/microscopy techniques for comprehensive structural analysis. The research methodologies established here are helpful for the investigation of dynamic alloy and interfacial structures and should inspire more efforts toward the simultaneous upgrading of CO 2 and shale gas.

03 NATURAL GAS↗
DOE OSTI2021

Basin scale distributions of dissolved manganese, nickel, zinc and cadmium in the Mediterranean Sea

Samples for dissolved trace metal concentrations were collected during GEOTRACES expedition GA04-N in summer/spring in the Mediterranean Sea, starting in the Atlantic Ocean and sampling both deep basins of the Mediterranean Sea. Outflow of Mediterranean Outflow Water leads to elevated concentrations of Mn, Ni and Zn in the Atlantic Ocean, but a concentration minimum in the Atlantic distribution of Cd. Nevertheless, when comparing the in- and outflow, the Mediterranean is a net source of Cd to the Atlantic Ocean. Surface concentrations of Mn, Ni, Zn and Cd are elevated in the Mediterranean relative to the Atlantic Ocean where Ni and Cd gradually increased along the eastward surface water flow path, Zn reached a homogenous concentration in the order of 1 to 1.5 nM, and Mn displayed a patchy surface distribution. The observed differences are attributable to the different dynamics of their biogeochemical cycling, notably the partitioning between the dissolved and particulate phases due to biological uptake, scavenging and possibly organic complexation. The elevated surface concentrations of Mn, Ni, Zn and Cd in the Mediterranean are derived from atmospheric deposition, where most likely Zn and Cd are mainly sourced from anthropogenic origin, Mn mostly from lithogenic origin and Ni from both anthropogenic and lithogenic origin. Dissolved Zn and Cd, as well as phosphate and nitrate, display striking inter-basin fractionations with elevated concentrations in the deep water of the western basin compared to the deep eastern basin, without a coinciding increase in the apparent oxygen utilization. Given that physical circulation or contribution from biogenic particulate metals cannot explain the elevated dissolved concentrations, an external non-biological source is required. This source is most likely a vertical flux of metal laden particles dissolving through the water column of the western Mediterranean where these particles, most likely from anthropogenic origin, are derived from either atmospheric deposition or particulate material deposited on the continental shelves that makes its way into the deep basin. Furthermore, to confirm or detect trends in dissolved metal concentrations in the deep basin, regular basin wide assessments of the trace metal distributions in the Mediterranean are needed. The distributions of Mn, Ni, Zn and Cd in the Sea of Marmara illustrate that all of these metals can be affected by anthropogenic surface sources and highlight the different susceptibilities of the dissolved metal distributions to supply from remineralization, and to removal through scavenging. This study provides a first baseline to assess future changes and underlines that the Mediterranean marine environment is susceptible to anthropogenic disturbances, with varying effects for different metals due to differing source strengths and biogeochemical cycles.

58 GEOSCIENCES↗
DOE OSTI2025

Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows

We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.

compressible flows↗
DOE OSTI2021

The many-body expansion for aqueous systems revisited: II. Alkali metal and halide ion – water interactions

We present a detailed study of the Many-Body Expansion (MBE) for alkali metal and halide ion-water interactions and quantify the effect of these ions on the strength of the surrounding aqueous hydrogen bonding environment. Building on our previous work on neutral water clusters [J. P. Heindel and S. S. Xantheas, J. Chem. Theor. Comput. 16 (11), 6843–6855 (2020)], we carry out the complete MBE for ion-water clusters, Z+/-(H2O)9, where Z = Li+, K+, Cs+, Cl-, Br-, I- and compare with the results for (H2O)10. The 2-B ion-water (I-W) interaction represents a larger percentage of the total cluster binding energy compared to a pure water cluster of the same size with the total 3-B term being smaller and of opposite sign (repulsive) whereas higher order terms are essentially negligible. The same oscillating behavior around zero for MBE terms higher than the 5-B with basis set that was reported for water clusters is also observed for the ion-water clusters considered with Basis Set Superposition Error (BSSE) corrections amending this like in the water cluster case. A remarkable, linear anti-correlation between the total 2-B ion-water and the total 2-B and 3-B water-water interactions is found, quantifying the effect of the different ions in disrupting and altering (weakening) the neighboring hydrogen bonded water network. Our results suggest a universal behavior of the two different families of ions (alkali metals and halides) for both the correlations of the various components of the total binding energies as well as the estimate of the 2-B BSSE correction which is reported to follow a common profile for ion-water and water-water interactions when cast in terms of reduced distances and energies. We expect our insights into the nature of BSSE will be important for future ab initio based, many-body molecular dynamics studies.

ion-water interactions, water-water interaction, m↗