Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability theory”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hybrid simulations of sub-cyclotron compressional and global Alfven Eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron frequency compressional (CAE) and global (GAE) Alfven Eigenmodes in low aspect ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam power in NSTX experiments, prompting further examination of the conditions for their excitation. Linear simulations are performed with the hybrid MHD-kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Three distinct types of modes are found in simulations -- co-CAEs, cntr-GAEs, and co-GAEs -- with differing spectral and stability properties. The simulations reveal that unstable GAEs are more ubiquitous than unstable CAEs, consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from simulations and estimated analytically for relevant sources not present in the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

Compressional Alfven Eigenmodes↗
DOE OSTI2023

Neutralizing topological obstructions to bubbles of nothing

Theories with compact extra dimensions can exhibit a vacuum instability known as a bubble of nothing. These decay modes can be obstructed if the internal manifold is stabilized by fluxes, or if it carries Wilson lines for background gauge fields, or if the instanton is incompatible with the spin structure. In each of these cases the decay can proceed by adding dynamical charged membranes or gauge fields. We give a general, bottom-up procedure for constructing approximate bubble of nothing solutions in models with internal spheres stabilized by flux and study the influence of the brane tension on the tunneling exponent, finding two branches of solutions that merge at a minimal superextremal value of the tension. In the case of Wilson operators and incompatible fermions, the relevant bubble is shown to be the Euclidean Reissner-Nordstrom black hole, and the ordinary decay exponent is modified by 1/g 2 effects. We examine the Dirac operator on this background and comment on the relevance for models of supergravity with gauged R-symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Feasibility of Prelithiation in LiFePO 4

Lithium iron phosphate (LiFePO 4 ) is widely applied as the cathode material for the energy storage Li-ion batteries due to its low cost and high cycling stability. However, the low theoretical specific capacity of LiFePO 4 makes its initial capacity loss more concerning. Therefore, lithium compensation by way of prelithiation and applications of sacrificial Li-rich additives in LiFePO 4 is imminent in elevating the energy density and/or prolonging the lifetime of the LiFePO 4 -based Li-ion batteries (LIBs). Prelithiation in LiFePO 4 is herein carried out by electrochemical and chemical methods and its feasibility is proved on the basis of the electrochemical evaluations such as the initial charge capacity and the cycling stability. In addition, the site of the pre-intercalated Li-ions is found via comprehensive physical characterizations and the density functional theory (DFT) calculations. Furthermore, these findings open a new avenue for elevating the energy density and/or prolonging the lifetime of the high-energy-density batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of Cr and Co on 99Tc retention in magnetite: A combined study of ab initio molecular dynamics and experiments

This work explores the effect of co-mingled dopants, Co(II) and Cr(III), on Tc(IV) incorporation and retention in magnetite when heat treated to 625 or 700 °C. Key trends in Tc retention in the high temperature regime were identified using a combination of density-functional-theory based ab initio molecular dynamics (AIMD) simulations, and batch experiments including solid phase characterization techniques, e.g. X-ray absorption spectroscopy. A stabilizing effect on Tc(IV) was observed when the number of Tc and Cr atoms are equal or when the magnetite surface is oversaturated with Tc and Cr inclusions. Here, oversaturation is hypothesized to force Cr from the magnetite surface to form a Cr2O3 phase, which may act as a protective layer that prevents Tc release. With the addition of Co, Tc(IV) is stabilized via redox processes. The presence of Cr in low concentrations interferes with this redox stabilization and Cr is preferentially stabilized as opposed to Tc. As a result, using Co as a stabilizing dopant for Tc may be compromised in the presence of Cr. When the relative concentration of Tc, Cr and Co is the same, or more Co atoms are added to high Cr incorporated systems, the formation of the Cr2O3 phase may be suppressed. Although waste streams with co-mingled Tc and Cr potentially may benefit from a Co dopant, since the formation of a Cr2O3 passivation layer may protect incorporated Tc from being released, the relative concentration of the three elements will be a critical parameter for maximizing effectiveness of this strategy.

Lee, Mal Soon↗
DOE OSTI2024

Pure-Water-Fed Forward-Bias Bipolar Membrane CO 2 Electrolyzer

Coupling renewable electricity to reduce carbon dioxide (CO 2 ) electrochemically into carbon feedstocks offers a promising pathway to produce chemical fuels sustainably. While there has been success in developing materials and theory for CO 2 reduction, the widespread deployment of CO 2 electrolyzers has been hindered by challenges in the reactor design and operational stability due to CO 2 crossover and (bi)carbonate salt precipitation. Herein, we design asymmetrical bipolar membranes assembled into a zero-gap CO 2 electrolyzer fed with pure water, solving both challenges. By investigating and optimizing the anion-exchangelayer thickness, cathode differential pressure, and cell temperature, the forward-bias bipolar membrane CO 2 electrolyzer achieves a CO faradic efficiency over 80% with a partial current density over 200 mA cm –2 at less than 3.0 V with negligible CO 2 crossover. In addition, this electrolyzer achieves 0.61 and 2.1 mV h –1 decay rates at 150 and 300 mA cm –2 for 200 and 100 h, respectively. Postmortem analysis indicates that the deterioration of catalyst/polymer– electrolyte interfaces resulted from catalyst structural change, and ionomer degradation at reductive potential shows the decay mechanism. Furthermore, all these results point to the future research direction and show a promising pathway to deploy CO 2 electrolyzers at scale for industrial applications.

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

Majorana zero modes in Y -shape interacting Kitaev wires

Motivated by the recent experimental realization of minimal Kitaev chains using quantum dots, we investigate the Majorana zero modes (MZM) in Y-shape Kitaev wires. We solve the associated Kitaev models analytically at the sweet spot (t h = Δ) and derive the exact form of MZM wave-functions in this geometry. We observed exotic multi-site MZMs located near the junction center, on the nearby edge sites of each leg. This result is important for potential braiding of Majoranas and the performance of Y-junctions made from arrays of quantum dots. Furthermore, we study the stability of local (single-site) and multi-site MZMs modes in the presence of Coulomb repulsion, using density matrix renormalization group theory. Our local density-of-states calculation shows that these multi-site MZMs are as equally topologically protected as the single-site MZMs when in the presence of Coulomb repulsion or when away from the sweet-spot.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structures and localized vibrational states of defects in graphite by tight-binding calculations

The structural and vibrational properties of pristine graphite and point defects in graphite are studied by tight-binding (TB) calculations using a three-center TB potential model. We showed that the three-center TB potential without ‘‘ad hoc’’ van der Waals interaction corrections can accurately describe the inter-layer distance of graphite and the lowest-energy structures and stabilities of typical point defects in graphite. The results from our TB calculations are in good agreement with those from density-functional theory calculations with van der Waals interaction corrections. We also investigated the vibrational properties to gain better understanding on the localization of vibrational states induced by the point defects. Our calculation results show that although localized or quasi-localized vibrational modes can be found in all defected graphite, the localization induced by Frenkel pair, dual-vacancy, and dual-interstitial defects is much stronger. Atomic displacements associated with the localized vibrational modes induced by these three point defects are also analyzed

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

Structures of glass-forming liquids by x-ray scattering: Glycerol, xylitol, and D-sorbitol

In this work, synchrotron X-ray scattering has been used to investigate three liquid polyalcohols of different sizes (glycerol, xylitol and D-sorbitol) from above the glass transition temperatures $T_g$ to below. We focus on two structural orders: the association of the polar OH groups by hydrogen bonds (HBs) and the packing of the non-polar hydrocarbon groups. We find that the two structural orders evolve very differently reflecting the different natures of bonding. Upon cooling from 400 K, the O$\cdots$O correlation at 2.8 Å increases significantly in all three systems, indicating more HBs, until kinetic arrest at $T_g$; the increase is well described by an equilibrium between bonded and non-bonded OH with ΔH = 9.1 kJ/mol and ΔS = 13.4 J/mol/K. When heated above $T_g$, glycerol loses the fewest HBs per OH, followed by xylitol and by D-sorbitol, in the same order the number of OH groups per molecule increases (3, 5, and 6). The pair correlation functions of all three liquids show exponentially damped density modulations of wavelength 4.5 Å, which are associated with the main scattering peak and with the intermolecular C$\cdots$C correlation. In this respect, glycerol is the most ordered with the most persistent density ripples, followed by D-sorbitol and by xylitol. Heating above $T_g$ causes faster damping of the density ripples with the rate of change being the slowest in xylitol, followed by glycerol and by D-sorbitol. Given the different dynamic fragility of the three liquids (glycerol being the strongest and D-sorbitol the most fragile), we relate our results to the current theories of the structural origin for the difference. We find that the fragility difference is better understood on the basis of the thermal stability of HB clusters than that of the structure associated with the main scattering peak.

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

Workhorse minimally empirical dispersion-corrected density functional with tests for weakly bound systems: r 2 SCAN + rVV 10

SCAN+rVV10 has been demonstrated to be a versatile van der Waals (vdW) density functional that delivers good predictions of both energetic and structural properties for many types of bonding. Recently, the r 2 SCAN functional was devised as a revised form of SCAN with improved numerical stability. In this work, we refit the rVV10 functional to optimize the r 2 SCAN+rVV10 vdW density functional and test its performance for molecular interactions and layered materials. Our molecular tests demonstrate that r 2 SCAN+rVV10 outperforms its predecessor SCAN+rVV10 in both efficiency (numerical stability) and accuracy. This good performance is also found in lattice-constant predictions. In comparison with benchmark results from higher-level theories or experiments, r 2 SCAN+rVV10 yields excellent interlayer binding energies and phonon dispersions for layered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fluctuations in Hill’s equation parameters and application to cosmic reheating

Cosmic inflation provides a compelling framework for explaining several observed features of our Universe, but its viability depends on an efficient reheating phase that converts the inflaton’s energy into Standard Model particles. This conversion often proceeds through nonperturbative mechanisms such as parametric resonance, which is described by Hill’s equation. In this work, we investigate how stochastic fluctuations in the parameters of Hill’s equation can influence particle production during reheating. We show that such fluctuations can arise from couplings to light scalar fields and can significantly alter the stability bands in the resonance structure, thereby enhancing the growth of fluctuations and broadening the region of efficient energy transfer. Using random matrix theory and stochastic differential equations, we decompose the particle growth rate into deterministic and noise-induced components and demonstrate analytically and numerically that even modest noise leads to substantial particle production in otherwise stable regimes. Furthermore, these results suggest that stochastic effects can robustly enhance the efficacy of reheating across a wide swath of parameter space, with implications for early Universe cosmology, UV completions involving multiple scalar fields, and the resolution of the cosmological moduli problem.

Cosmology↗
DOE OSTI2020

Robust Control of Wave Energy Converters Using Unstructured Uncertainty

In the design of ocean wave energy converters, proper control design is essential to the maximization of the power generation performance for the device. However, in realistic applications, this control design must be undertaken in the presence of model uncertainty. This paper considers the use of robust control theory to optimize the nominal performance for a wave energy converter in stochastic waves, subject to the constraint that the controller be stability-robust to unstructured uncertainties. We formulate the problem as a multi-objective optimal control problem, in which the primary objective is the maximization of power generation for the nominal system, and the competing objective is the norm of the uncertainty input/output channel. This optimal control problem is nonconvex, and we therefore propose an iterative algorithm that can be used to arrive at a local optimal solution. This iterative approach employs the concept of Iterative Convex Overbounding, in the context of the classical Method of Centers. Here, the methodology is demonstrated on a model of a single, buoy- type wave energy converter.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

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

Stability and metallization of solid oxygen at high pressure

The phase diagram of oxygen is investigated for pressures from 50 to 130 GPa and temperatures up to 1200 K using first-principles theory. A metallic molecular structure with the P6 3 /mmc symmetry (η' phase) is determined to be thermodynamically stable in this pressure range at elevated temperatures above the ε(O 8 ) phase. Crucial for obtaining this result is the inclusion of anharmonic lattice dynamics effects and accurate calculations of exchange interactions in the presence of thermal disorder. Overall, we present analysis of electronic, structural, and thermodynamic properties of solid oxygen at 0 K and finite temperature with hybrid exchange functionals, including a comparison with available experimental data.

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

Understanding AMOC stability: the North Atlantic Hosing Model Intercomparison Project

Abstract. The Atlantic meridional overturning circulation (AMOC) is an important part of our climate system. The AMOC is predicted to weaken under climate change; however, theories suggest that it may have a tipping point beyond which recovery is difficult, hence showing quasi-irreversibility (hysteresis). Although hysteresis has been seen in simple models, it has been difficult to demonstrate in comprehensive global climate models. Here, we outline a set of experiments designed to explore AMOC hysteresis and sensitivity to additional freshwater input as part of the North Atlantic Hosing Model Intercomparison Project (NAHosMIP). These experiments include adding additional freshwater (hosing) for a fixed length of time to examine the rate and mechanisms of AMOC weakening and whether the AMOC subsequently recovers once hosing stops. Initial results are shown from eight climate models participating in the Sixth Coupled Model Intercomparison Project (CMIP6). The AMOC weakens in all models as a result of the freshening, but once the freshening ceases, the AMOC recovers in half of the models, and in the other half it stays in a weakened state. The difference in model behaviour cannot be explained by the ocean model resolution or type nor by details of subgrid-scale parameterisations. Likewise, it cannot be explained by previously proposed properties of the mean climate state such as the strength of the salinity advection feedback. Instead, the AMOC recovery is determined by the climate state reached when hosing stops, with those experiments where the AMOC is weakest not experiencing a recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Hydrogen generation via ammonia decomposition on highly efficient and stable Ru-free catalysts: approaching complete conversion at 450 °C

We report Hydrogen (H 2 ) is a prospective zero-carbon and high-energy-density fuel alternative to fossil fuels for generating power and clean energy. Ammonia (NH 3 ) is a promising H 2 (17.7%) carrier, which can easily overcome the challenges associated with H 2 storage and transportation. Thermocatalytic ammonia decomposition reaction (ADR) is an effective way to produce clean H 2 but it relies on the use of expensive and rare ruthenium (Ru)-based catalysts at elevated temperatures (>500 °C), hence is not sustainable and economically feasible. Herein, we report a synergistic strategy to design a heterostructured Ru-free catalyst, consisting of CoNi alloy nanoparticles well-dispersed on a MgO–CeO 2 –SrO mixed oxide support with potassium promotion. The resulting K–CoNi alloy –MgO–CeO 2 –SrO catalyst presents 97.7% and 87.5% NH 3 conversion efficiency at 450 °C at gas hourly space velocities (GHSVs) of 6000- and 12 000-mL h -1 g cat -1 , respectively. At 500 °C, the H 2 production rate (57.75 mmol g cat -1 min -1 ) becomes comparable to that of most of the reported Ru-based catalysts. The catalyst stability has been successfully demonstrated in both a fixed-bed reactor under high pressure (120 h at 5.0 bar) and a membrane reactor prototype (600 h at 1.5 bar) at 500 °C. High-temperature in situ XPS analysis, temperature-programmed desorption/reduction, and density functional theory calculations have been carried out to elucidate the possible active sites and performance enhancement mechanisms. This work highlights the importance of constructing optimal interfaces between active metal nanoparticles and oxide support for boosting the NH 3 to H 2 conversion efficiency and long-term stability.

25 ENERGY STORAGE↗
DOE OSTI2022

The solution structures and relative stability constants of lanthanide–EDTA complexes predicted from computation

Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La 3+ , Eu 3+ , Gd 3+ , and Lu 3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln–EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ~7 and ~11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln–EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln 3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ~11), as well as for a single Ln 3+ ion (La) with varying ligands (EDTA at pH ~7 and ~11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln–EDTA complex resulting in a reduction of the complex stability. As a result, increased Ln–ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature.

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

Spontaneous magnon decays from nonrelativistic time-reversal symmetry breaking in altermagnets

Quasiparticles are central to condensed matter physics, but their stability can be undermined by quantum many-body interactions. Magnons, i.e., quasiparticles in quantum magnets, are particularly intriguing because their properties are governed by both real and spin space. While crystal symmetries may be low, spin interactions often remain approximately isotropic, limiting spontaneous magnon decay. Textbook wisdom holds that collinear Heisenberg magnets follow a dichotomy: ferromagnets host stable magnons, while antiferromagnetic magnons may decay depending on dispersion curvature. Up to now, relativistic spin-orbit coupling and noncollinear order that connect spin space to real space were shown to introduce more complex magnon instability mechanisms. Here, we show that even in nonrelativistic isotropic collinear systems, this conventional dichotomy is disrupted in altermagnets. Altermagnets, a newly identified class of collinear magnets, exhibit compensated spin order with nonrelativistic time-reversal symmetry breaking and even-parity band splitting. Using kinematic analysis, nonlinear spin-wave theory, and quantum simulations, we reveal that even weak band splitting opens a decay phase space, driving quasiparticle breakdown. Additionally, 𝑑-wave altermagnets form a rare “island of stability” at the Brillouin-zone center. Furthermore, our findings establish a quasiparticle stability trichotomy in collinear Heisenberg magnets and position altermagnets as a promising platform for unconventional spin dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗