Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “favorability”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Decoupling size and surface effects of intermetallic CuPd nanocrystals for electrocatalytic nitrate reduction to ammonia

Nitrate pollution poses a major environmental challenge, but its electrochemical conversion to ammonia offers a sustainable waste-to-value solution. Here, in this study, we synthesized monodisperse, size-tunable B2-phase CuPd intermetallic nanocrystals (6–46 nm) and studied their performance in the electrochemical nitrate reduction reaction (eNO 3 RR). By using bromide ions to modulate Pd reduction and applying mild annealing, we achieved phase-pure B2 structures across all sizes. Catalytic testing revealed a volcano-like trend in ammonia yield, peaking at 33 nm nanocubes with a rate of 6.97 mol h −1 g −1 at −0.6 V vs. reversible hydrogen electrode (RHE). This optimum reflects a balance between the increased surface area of smaller particles and the enhanced exposure of active (100) facets in larger ones. Theoretical calculations indicated that the B2-CuPd (100) facet is favorable for nitrate adsorption, thereby supporting the high activity of nanocubes. Our results highlight the critical role of tuning both nanoparticle size and surface structure to maximize eNO 3 RR efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shifting modes of iron sulfidization at the onset of OAE-2 drive regional shifts in pyrite δ 34 S records

We report total reduced inorganic S isotope ratios (δ 34 S CRS ) shift toward more negative values across much of the southern North Atlantic just before the onset of the Cenomanian-Turonian Ocean Anoxic Event (OAE-2). At the same time, there is no parallel isotopic change in the significantly larger pool of kerogen (organic) S, which indicates that the distribution and S-isotope composition of sulfide in the environment likely did not drive the change in δ 34 S CRS . Here, we investigate possible explanations for the negative shift in δ 34 S CRS values and their divergence from organic S by isolating iron sulfides for morphological identification and grain-specific isotopic analysis using secondary ion mass spectrometry (SIMS). In pre- and syn-OAE-2 sedimentary rocks from Demerara Rise, we find four distinct morphologies of iron sulfides: pyrite framboids (1–20 μm diameter), irregular pyrite aggregates (1–38 μm diameter), large cemented pyrite aggregates (~60 μm diameter), and irregular and cemented aggregates of the pyrite polymorph marcasite (1–45 μm diameter). These different textural groups have distinct S-isotopic compositions that are largely consistent through the onset of OAE-2. As such, the secular change in bulk δ 34 S CRS values likely reflects the changing proportions of these phases stratigraphically across OAE-2. All textural groups feature resolvable intra-grain δ 34 S variability, suggesting that the environments in which they formed were characterized by dynamic sulfide δ 34 S values and/or by partial closed-system distillation. We use grain-specific δ 34 S distributions to rule out shoaling of the chemocline within the sediments as a mechanism for the observed decrease in δ 34 S CRS . Instead, we propose that changes in the reactivity of the iron species delivered to Demerara Rise over the ~200 kyr leading up to the onset of OAE-2 impacted the relative contributions of pyrite with S-isotope signatures reflecting the water column, shallow sediments, and deeper sediments to the bulk sedimentary δ 34 S CRS value. Specifically, the change in iron reactivity at the onset of OAE-2 favored the production of 34 S-depleted large, cemented aggregates and framboids at the expense of more 34 S-enriched irregular aggregates. Our results underscore that bulk δ 34 S CRS measurements integrate multiple reduced phases that form via distinct reaction mechanisms and potentially in different parts of the depositional environment. Grain-specific SIMS analyses dramatically enrich our ability to interpret pyrite isotopic patterns in the geologic record.

58 GEOSCIENCES↗
DOE OSTI2026

Mobility of heavy rare earth elements in magmatic-hydrothermal systems: Experimental determination of DyPO 4 solubility in supercritical fluids

Dysprosium (Dy) is part of the heavy rare earth elements (REE) and a critical component in the production of high-performance super magnets used in wind turbines, electric vehicles, and other green technologies. The HREE phosphate, xenotime, is commonly associated with hydrothermal systems, where the mobility, fractionation, and deposition of REE are controlled by the stability of aqueous complexes and solubility of REE minerals which depend on temperature, pressure, and fluid composition (salinity, pH, and ligand availability). Here, we conducted a series of synthetic DyPO 4 solubility experiments in NaCl-HCl-bearing aqueous solutions (pH of 2 at 25 °C) from 500 to 700 °C and 1.5 kbar using externally heated pressure vessels. DyPO 4 displays a strong prograde solubility with temperature across all studied fluid salinities (0–0.5 m NaCl). The DyPO₄ solubility isotherms at 500 °C, 600 °C, and 700 °C show a complex relationship with salinity where DyPO₄ solubility is highest (39.2–1552 ppm at 500–700 °C) in 0 m NaCl solutions, decreases strongly (21–331 ppm at 500–700 °C) in the 0.06 m NaCl solutions and remains relatively constant in the 0.25 m (42.8–248 ppm) and 0.5 m NaCl (40.2–124 ppm) solutions. Thermodynamic modeling suggests that pH increases as a function of temperature and salinity, and controls solubility due to the stability of chloride complexes (i.e., DyCl +2 and DyCl 2 + ) at 500–600 °C, transitioning to hydroxyl complexes (i.e., Dy(OH) 3 0 ) dominating at 700 °C and relatively alkaline pH. Comparison with NdPO 4 solubility data indicates that the light REE in monazite is significantly more soluble than the heavy REE in xenotime in saline fluids at high temperature (500–700 °C), implying efficient light/heavy REE fractionation in supercritical fluids. The lower solubility of DyPO 4 relative to NdPO 4 suggests that at high temperatures, acidic (pH 2–4), low-salinity (∼5 wt% NaCl) fluids will favor the precipitation of xenotime over monazite as it will cool from 700 °C to below 500 °C. These observations align with paragenetic mineral assemblages in NYF pegmatites (e.g., Baveno, Italy), where xenotime commonly immobilizes heavy REE during the high-temperature hydrothermal stages.

58 GEOSCIENCES↗
DOE OSTI2026

Ab initio thermodynamics of Ni and Co incorporation in Mg hydroxide, carbonate, and hydroxycarbonate minerals

Ni and Co are critical elements needed for modern technologies, and a better understanding of the ability of Mg-based minerals to incorporate these elements would benefit strategy development for Ni and Co recovery from mafic and ultramafic deposits. Here, in this work, we performed density functional theory (DFT) calculations of Ni and Co incorporation in six potential products of the carbonation of mafic and ultramafic silicates: brucite (Mg(OH) 2 ), magnesite (MgCO 3 ), nesquehonite (MgCO 3 ⸱3H 2 O), lansfordite (MgCO 3 ⸱5H 2 O), artinite (Mg 2 CO 3 (OH) 2 ⸱3H 2 O), and hydromagnesite (Mg 5 (CO 3 ) 4 (OH) 2 ⸱4H 2 O). The DFT results were used in an ab initio thermodynamics framework to explore the pH 2 O–pCO 2 conditions at which the Mg-based minerals were predicted to be thermodynamically stable and to quantify the Gibbs free energy of Ni and Co substitution at Mg sites. Among the six Mg-based minerals, brucite and magnesite were predicted to have the lowest Ni and Co substitution free energy. An analysis of the effect of temperature indicated that, at low temperature (<100 K), brucite more readily accommodated Ni and Co, while, at higher temperature (>335 K), magnesite more favorably incorporated Ni and Co. Between 100 K and 335 K, Ni was predicted to preferentially substitute for Mg in brucite and Co for Mg in magnesite, thus leading to a driving force for separating Ni and Co in conditions where brucite and magnesite both form. Insights gained in this work could therefore help select experimental conditions that either promote or inhibit incorporation of these critical elements into Mg-based mineral phases.

Critical elements↗
DOE OSTI2020

Born effective charges and electric polarization in bulk ε-Fe 2 O 3 : An ab-initio approach

Limitation in the synthesis of single-crystal has restricted the complete understanding of the magnetoelectric properties of bulk ε-Fe 2 O 3 . In order to understand the electric polarization of bulk ε-Fe 2 O 3 , we have employed a first-principles Berry-phase method and the density functional perturbation theory to calculate the spontaneous electric polarization and Born effective charges (BECs) of ε-Fe 2 O 3 . The polarization quanta of the ferroelectric bulk ε-Fe 2 O 3 and its polarization lattice are also reported. The polarization switching through a paraelectric transition phase of a low barrier is found to be favorable. The polarization and BECs calculated by both methods are in good agreement.

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

Analysys of the performance a PEM-type electrolyzer in variable energy supply conditions

Here in this article, the main variables related to the performance of a proton exchange membrane-type electrolyzer are analyzed along with the risk of explosion for different operating conditions, notably when the electrolyzer works with a variable solar energy power supply and at high pressure. A novel mathematical model was developed that considers all mass transport phenomena through the membrane electrode assembly (MEA) for water, hydrogen, and oxygen, as well as the water evaporation in the anode and cathode channels. The main findings of this work are that for a solar energy supply operating at high pressure, the electrolyzer cannot function safely throughout the day, since in the hours when solar irradiation is lower, low current densities favor high hydrogen concentrations at the anode. Finally, it was found that greater thicknesses of the membrane increased the explosion risk due to the increase in the MEA temperature and therefore the hydrogen solubility.

08 HYDROGEN↗
DOE OSTI2026

Temperature-dependent changes in gas chromatographic separation metrics for trihexyl(tetradecyl)phosphonium-based ionic liquid stationary phases and comparison to conventional polysiloxane stationary phases

Here, to assess whether the trihexyl(tetradecyl)phosphonium chloride ([P 66614 + ][Cl - ]) ionic liquid (IL), employed as a gas chromatographic stationary phase, exhibits temperature-dependent phase transitions or structural heterogeneity, retention thermodynamics for probe molecules of varied chemical structure were investigated by gas chromatography (GC) over the temperature range of 25 to 64 °C using van’t Hoff analysis. All van’t Hoff plots for n-alkanes (C8–C12) and other probe analytes were linear, indicating the absence of detectable phase-transition behavior. The [P 66614 + ][Cl - ] IL was further compared with the [P 66614 + ] bis[(trifluoromethyl)sulfonyl]imide ([P 66614 + ][NTf 2 - ]) IL and commercial polysiloxane stationary phases in terms of temperature-dependent resolution behavior and chromatographic selectivity over the temperature range of 28–52 °C. The resolution factors of analytes in mixtures obtained using gas chromatography–mass spectrometry (GC–MS) and the summation of resolution values decreased smoothly with temperature on both IL stationary phases and on the dimethylphenylcyano-substituted polymeric stationary phase (OV-1701), whereas the poly(50% diphenyl/50% dimethyl siloxane) stationary phase (SPB-50) exhibited irregularities in its temperature-dependent summation of resolution profile, showing a distinct maximum in the temperature range of 34–40 °C. Selectivity values for the IL stationary phases were largely determined by the anion, with the [P 66614 + ][Cl - ] IL providing enhanced separation between polar and aromatic analytes, while the [P 66614 + ][NTf 2 - ] IL favored dispersive and π-π interactions that resulted in a reversal of elution order for some aliphatic–aromatic and alkyne-containing analyte pairs. Capillary columns prepared using different surface activation methods and coated with identical [P 66614 + ][Cl - ] IL stationary phase showed indistinguishable temperature-dependent trends in summation of resolution value and selectivity. An evaluation of analyte-specific selectivity provides insight into the characteristic separation behavior of the stationary phases and enables identification of analyte classes that are more efficiently separated on the [P 66614 + ]-based IL stationary phases compared to conventional stationary phase materials.

Chromatographic selectivity↗
DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗
DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An optimization-based approach to tailor the mechanical response of soft metamaterials undergoing rate-dependent instabilities

An optimization-based design framework is proposed to tune the response of soft metamaterials involving both geometric instabilities and nonlinear viscoelastic material behavior. Designing the response of soft metamaterials to harness instabilities and undergo large, tailored configuration changes will enable advancements in soft robotics, shock and vibration mitigation, and flexible electronics. In line with the metamaterial concept, the response of these materials is governed to a large extent by the geometric and topological makeup of their small-scale features. However, the link between structure and response is less intuitive for soft metamaterials due to their reliance upon highly nonlinear responses triggered by geometric instabilities. This is further complicated by the effects of viscoelastic relaxation, which recent studies have shown to alter the emergence of instabilities in non-intuitive ways. Here, these effects are accounted for in our framework to achieve various design objectives, including tailored force–displacement response and maximized energy absorption from both geometric and material effects. To fully automate this process, it is essential to have a completely robust equation solver for forward problems involving instabilities and viscoelastic relaxation. We achieve this by casting the search for stable mechanical equilibrium — i.e. the forward problem — as a minimization problem and utilize a trust region algorithm to robustly handle instabilities and follow energetically-favorable equilibrium paths through critical points.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Techno-economic analyses of several redox flow batteries using levelized cost of energy storage

Development of inexpensive long-duration energy storage supports widespread deployment of variable renewable energy resources onto the electricity grid. Flow batteries are a promising class of devices for long-duration energy storage. Techno-economic modeling is needed to evaluate commercial feasibility of existing technologies and to help guide research and development of new technologies. Levelized cost of storage is a useful metric that accounts for capital and operating costs and energy throughput over the life of a project. This metric is used to compare the economic prospects of lithium ion to eight aqueous and two hypothetical nonaqueous flow batteries in four use cases. Finally, flow batteries with inexpensive active materials are favored in long-duration applications.

25 ENERGY STORAGE↗
DOE OSTI2022

Hydration of alumina (Al 2 O 3 ) toward advancing aluminum particles for energy generation applications

Transforming metal particle combustion may require alteration of the metal oxide passivation shell surrounding the metal core. One approach for aluminum (Al) relies on the hydrated form of the metal oxide to incite surface reactions. Here, this study explores the conditions required for hydrating alumina (Al 2 O 3 ) particles and then extends those conditions toward hydration of the Al 2 O 3 passivation layer surrounding an Al core particle. By raising the pH of the water slurry to 11 and controlling temperature and time in slurry (i.e., aging), aluminum hydroxide Al(OH) 3 formation from Al 2 O 3 particles was optimized. The procedure was then extended to Al nanoparticles (nAl). Even though heating and extended aging time in slurry proved advantageous for Al 2 O 3 particles, these conditions were discarded for nAl particles because they favored formation of AlOOH, a less desirable hydrate. Through microscopy, spectroscopy, X-ray diffraction, and thermal analyses, results indicate that for a pH range of 11.26–11.56, the original Al 2 O 3 shell on Al particles transformed into Al(OH) 3 . Results indicate a feasible path forward towards producing new shell chemistries that may result in more directed energy from metal particle combustion.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Conversion of acenaphthalene to phenalene via methylation: A theoretical study

Ab initio calculations of the C 13 H 10 and C 13 H 9 potential energy surfaces related to the reaction of 1-acenaphthyl and methyl radicals and secondary isomerization of C 13 H 9 primary radical products have been performed at the chemically accurate G3(MP2,CC)//B3LYP/6-311G** level of theory to unravel the mechanism of conversion of acenaphthalene to phenalene or phenalenyl radical + H. The computed energetics and molecular parameters were utilized in Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) calculations of reaction rate constants and relative product yields. The 1-acenaphthyl + CH 3 reaction is predicted to proceed by a fast radical-radical recombination mechanism and to predominantly produce collisionally stabilized 1-methylacenaphthalene B1 or a C 13 H 9 benzylic radical A1 + H with exothermicities of 106.8 and 25.2 kcal/mol, respectively. The A1 + H channel is preferable at higher temperatures, whereas the stabilization of 1-methylacenaphthalene is favored at higher pressures. The radical A1 can nearly irreversibly interconvert to phenalenyl radical AP via a 22.7 kcal/mol exothermic isomerization process involving formal insertion of the CH 2 group into a C-C bond of the five-member ring leading to the expansion of this ring to a six-member ring. The rate constants for the H addition reaction to phenalenyl radical to form phenalene and the reverse H loss from phenalene were also computed and the strength of the weakest C-H bond in the CH 2 group of phenalene is evaluated as 62.2 kcal/mol. The analysis of the reaction kinetics allowed us to deduce a mechanism for the conversion of 1-acenaphthyl radical to phenalene or phenalenyl radical + H via methylation involving the formation of A1 via a well-skipping channel or stabilization and dissociation of (or H abstraction from) B1 followed by isomerization of A1 to phenalenyl AP, which can add an H atom producing phenalene. As a result, rate constants for the significant elementary reactions are fitted to modified Arrhenius expressions and are proposed for kinetic modeling of the expansion of a five-member ring on a zigzag edge of PAH to a six-member ring by methylation.

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

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗
DOE OSTI2020

Understanding how chemical structure affects ignition-delay-time $\phi$-sensitivity

$\phi$-sensitivity is the change in ignition delay time (IDT) with respect to the fuel-to-air equivalence ratio ($\phi$). High $\phi$-sensitivity is a desirable fuel property for applications in advanced compression ignition and multi-mode engine designs. Understanding how $\phi$-sensitivity depends on chemical structure is essential for selecting promising biofuels from the ever-growing list of proposed candidates. Here, we investigate the effect of chemical structure on $\phi$-sensitivity with experiment, simulation, and theory. Experimental Advanced Fuel Ignition Delay Analyzer (AFIDA) measurements for 2,4-dimethylpentane and diisopropyl ether provide evidence that branching and functional groups strongly impact $\phi$-sensitivity. Further insights into this dependence are obtained with 0-D kinetic simulations with existing mechanisms for n-pentane, diethyl ether, 3-pentanone, n-heptane, 2-methylhexane, 2,4-dimethylpentane, and 2,2,3-trimethylbutane. Quantum mechanical (QM) G4 calculations of low-temperature reactions help explain the observed experimental and simulation trends. Specifically, these QM calculations provide theoretical estimates of the ketohydroperoxide (KHP) dissociation rates, the HO 2 formation rates from peroxy radical (ROO), and the “cross-over” temperatures, i.e., the temperature at which ROO dissociation is favored compared to hydroperoxyl radical (QOOH) formation. Each of these reaction rates is compared to the n-alkane reference point to determine the impact of branching and different functional groups. Although kinetic mechanisms typically assume that KHP dissociation rates are invariant of chemical environment, our QM results suggest that this rate can span a range of roughly two orders of magnitude. We also discuss the importance of including the peroxy-hydroperoxy (OO-OOH) hydrogen transfer reaction for branched ethers. Finally, the insights gained assist in proposing a highly $\phi$-sensitive compound, namely, isopropyl propyl ether.

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

An experimental and kinetic modeling study of the ignition of 2-methyl decane

Ignition delay times (IDTs) of 2-methyl decane (C 11 H 24 -2) are measured in a high-pressure shock tube and in a rapid compression machine at equivalence ratios in the range 0.5–2.0 at 90 % dilution, at temperatures in the range 600–1430 K and at pressures of 15 and 30 bar. To clarify the effect of the branched methyl group on fuel oxidation, IDTs of n-undecane (nC 11 H 24 ) are also measured at similar conditions to those measured for C 11 H 24 -2. A new chemical kinetic mechanism, using C3MechV4.0.1 as the core chemistry, is developed and validated against the new experimental data. The thermodynamic properties of the fuel (RH), alkyl (Ṙ), alkyl peroxy (RȮ 2 ), hydroperoxy-alkyl (Q̇OOH), and peroxy hydroperoxy alkyl (Ȯ 2 QOOH) radicals are updated in both the C 11 H 24 -2 and nC 11 H24 models using THERM25. A reaction path flux analysis for C 11 H 24 -2 at different temperatures was conducted. Compared to nC 11 H 24 , C 11 H 24 -2 shows slower reactivity. At low and intermediate temperatures, the chain propagation pathway Q̇OOH ↔ C 11 cyclic ether + ȮH is favored for C 11 H 24 -2, while the chain branching pathway Q̇OOH ↔ Ȯ 2 QOOH ↔ C 11 carbonyl hydroperoxide (KHP) + ȮH is suppressed, leading to lower reactivity compared to nC 11 H 24 . At high temperatures, the presence of the branched methyl group inhibits the direct decomposition of the fuel, resulting in reduced C 2 H 4 formation, which in turn suppresses the reactivity of the fuel.

2-methyl decane↗
DOE OSTI2022

Leaching Model of Radionuclides in Metal-Organic Framework Particles

Metal-organic frameworks (MOFs) have been used to sequester radionuclides and seal them inside of porous scaffolds using postsynthetic modification procedures. Experiments show that certain Zr-MOF with different capping linkers significantly affects the radionuclide release kinetics. In this work, we developed a leaching model of radionuclides in Zr-MOF particles. The model assumes that uranyl species occupy two energetically favored sites: the metal node and the MOF pores. For a given overall concentration of uranyl species, the partitions of uranyl species at the metal nodes and within the pores are determined by their chemical potentials. The model also considers the effect of particle surface and concentration on chemical potentials and diffusivity. The effects of spatial and structural dependent chemical potentials and diffusivity as well as particle sizes on leaching kinetics are investigated with the model. Predicted and measured uranyl leaching kinetics in Zr-MOF particles under batch experiments are compared. The results demonstrate the model’s capability for exploring the mechanisms of leaching and provide guidance for material design.

Metal-Organic Framework, Diffusion, uranyl, Leachi↗
DOE OSTI2021

2D Ni 0.25 Mn 0.75 O 2 : A high-performance cathode for multivalent ion batteries

Although Li-ion batteries have driven portable energy storage in recent decades, there is increasing concern about their safety, cost, and abundance of constituents. Multivalent ion batteries (MVIBs) have the potential to remedy these issues, but they are limited by the currently known MVIB cathodes, which fail to deliver unanimously favorable voltage, energy density, and diffusion kinetics. Here we used density functional theory (DFT) to model the performance of Li, Na, Mg, Ca, and Al ions when paired with 2D Ni 0.25 Mn 0.75 O 2 , a novel cathode that uses increased layer separation to improve on the kinetics of its 3D analog. Our calculations yielded maximum voltages of 3.38 V for Na and 2.7 V for Ca, outperforming 2D Na x MnO 2 and Na x NiO 2 . Diffusion barriers for Li, Na, and Ca are below 300 meV, comparable to existing battery technology and the endpoint 2D cathodes; meanwhile, Mg and Al have prohibitively high diffusion barriers, implying their incompatibility with this cathode. Lastly, density of states calculations and Bader charge analysis show that the cathode becomes conducting following ion adsorption, which is necessary for high-rate performance. 2D Ni 0.25 Mn 0.75 O 2 maintains performance seen with other 2D transition metal oxides while increasing cathode conductivity, indicating that it is a promising candidate for experimental investigation with Li, Na, and Ca ions.

25 ENERGY STORAGE↗