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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multiple incommensurate magnetic states in the kagome antiferromagnet Na 2 Mn 3 Cl 8

The kagome lattice can host exotic magnetic phases arising from frustrated and competing magnetic interactions. However, relatively few insulating kagome materials exhibit incommensurate magnetic ordering. Here, we present a study of the magnetic structures and interactions of antiferromagnetic Na 2 Mn 3 Cl 8 with an undistorted Mn 2+ kagome network. Using neutron-diffraction and bulk magnetic measurements, we show that Na 2 Mn 3 Cl 8 hosts two different incommensurate magnetic states, which develop at T N1 = 1.6 K and T N2 = 0.6 K. Magnetic Rietveld refinements indicate magnetic propagation vectors of the form q=(q x , q y , $\frac{3}{2}$), and our neutron-diffraction data can be well described by cycloidal magnetic structures. By optimizing exchange parameters against magnetic diffuse-scattering data, we show that the spin Hamiltonian contains ferromagnetic nearest-neighbor and antiferromagnetic third-neighbor Heisenberg interactions, with a significant contribution from long-ranged dipolar coupling. This experimentally determined interaction model is compared with density-functional-theory simulations. Using classical Monte Carlo simulations, we show that these competing interactions explain the experimental observation of multiple incommensurate magnetic phases and may stabilize multi-q states. Here, our results expand the known range of magnetic behavior on the kagome lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Complex Density Wave Orders and Quantum Phase Transitions in a Model of Square-Lattice Rydberg Atom Arrays

In this work, we describe the zero-temperature phase diagram of a model of a two-dimensional square-lattice array of neutral atoms, excited into Rydberg states and interacting via strong van der Waals interactions. Using the density-matrix renormalization group algorithm, we map out the phase diagram and obtain a rich variety of phases featuring complex density wave orderings, upon varying lattice spacing and laser detuning. While some of these phases result from the classical optimization of the van der Waals energy, we also find intrinsically quantum-ordered phases stabilized by quantum fluctuations. These phases are surrounded by novel quantum phase transitions, which we analyze by finite-size scaling numerics and Landau theories. Our work highlights Rydberg quantum simulators in higher dimensions as promising platforms to realize exotic many-body phenomena.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

An active, stable cubic molybdenum carbide catalyst for the high-temperature reverse water-gas shift reaction

Although technologically promising, the reduction of carbon dioxide (CO 2 ) to produce carbon monoxide (CO) remains economically challenging owing to the lack of an inexpensive, active, highly selective, and stable catalyst. We show that nanocrystalline cubic molybdenum carbide (α-Mo 2 C), prepared through a facile and scalable route, offers 100% selectivity for CO 2 reduction to CO while maintaining its initial equilibrium conversion at high space velocity after more than 500 hours of exposure to harsh reaction conditions at 600°C. The combination of operando and postreaction characterization of the catalyst revealed that its high activity, selectivity, and stability are attributable to crystallographic phase purity, weak CO-Mo 2 C interactions, and interstitial oxygen atoms, respectively. Mechanistic studies and density functional theory (DFT) calculations provided evidence that the reaction proceeds through an H 2 -aided redox mechanism.

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

NSTX-U National Research Program: White Paper in Response to Call from FESAC Sub-Committee

Both scientific and technical innovation is needed for the realization of an attractive engineering solution for a timely and cost-effective Pilot Plant, the design and construction of which is the overarching recommendation of the FESAC Long Range Plan, and the 2021 NASEM Pilot Plant reports, which underpin the Bold Decadal Vision. The two most significant plasma physics gaps to close for a Compact Pilot Plant (CPP) are core confinement improvement and heat flux mitigation, neither of which have been closed in an integrated fashion for any planned fusion power production device. High core confinement and stability are essential for producing majority self-driven plasmas in CPPs with reduced size and auxiliary heating power requirements, with an improvement in confinement being the major driver for cost reduction of a CPP. The National Spherical Tokamak Experiment - Upgrade (NSTX-U) is a unique low aspect ratio research facility that will address the fundamental challenge of developing the science and technology basis for a CPP design that integrates high core and edge confinement with the ability to mitigate very high incident heat fluxes. NSTX-U capabilities will enable the high performance, already achieved on NSTX, to extend into physics regimes much closer to those anticipated in Spherical Tokamak (ST)-based CPPs. These confinement and stability properties will be assessed by a full complement of diagnostics and analysis tools, which will also aid in the development of the underlying theory and predictive models needed for further optimization. Both conventional and transformative heat flux mitigation methods, such as liquid lithium plasma-facing components, will be developed and tested in-situ in NSTX-U at incident heat fluxes of ~100 MW/m 2 , and will inform plans and reduce risk for a subsequent major upgrade to the device to fully heated, high-Z wall and full liquid lithium divertor capability, a technology that potentially could then be implemented on any magnetic confinement device at any aspect ratio. NSTX-U research is fully complementary to programs performed on other STs, nationally and internationally. Furthermore, NSTX-U research has a direct connection to the private sector by informing design choices for future power production facilities being developed by these companies. The NSTX-U program will operate as a national User Facility, with collaborating researchers, engineers, and graduate students from 19 outside institutions, and open to participation and experiments led by researchers from both public and private entities. The research program will advance workforce development through training of young scientists, engineers, and technicians, and it will also serve for further diagnostic innovation, especially for high heat flux and high-Z wall environments, and implementation of advanced artificial intelligence (AI) for plasma and heat flux control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Randomized Algorithms for Linear Solvers

Recently, randomized algorithms in numerical linear algebra, specifically those centered around random sketching, have gained traction in primarily theoretical research due to their potential to significantly reduce problem dimensionality at the cost of an O(1) multiplicative distortion factor. It has been assumed that this sketching can be done efficiently, but thorough investigation into how precisely to do it has been neglected. Moreover, the theory-based community has argued for sketching’s ability to reduce computational cost via complexity analysis, but has not researched how it affects the stability of the algorithms. At Sandia, efficient linear solvers that scale well on modern HPC architectures while maintaining stability are imperative for practical applications. In this LDRD, we developed a random sketching strategy that is substantially faster than existing ones, and demonstrate its superior performance in practice on a NVIDIA H100 GPU. Moreover, we show how this can be used to significantly outperform existing linear least squares solvers while improving the solver’s stability as well. Additionally, we demonstrate how this sketching strategy can be used to make a fast, stable QR factorization that can subsequently be used in s-step and block Krylov solvers. Finally, we incorporate a sketching-based block orthogonalization scheme into s-step GMRES, which is stable and faster than existing approaches on the Perlmutter supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Machine Learning-Based Model Predictive Control of Two-Time-Scale Systems

In this study, we present a general form of nonlinear two-time-scale systems, where singular perturbation analysis is used to separate the dynamics of the slow and fast subsystems. Machine learning techniques are utilized to approximate the dynamics of both subsystems. Specifically, a recurrent neural network (RNN) and a feedforward neural network (FNN) are used to predict the slow and fast state vectors, respectively. Moreover, we investigate the generalization error bounds for these machine learning models approximating the dynamics of two-time-scale systems. Next, under the assumption that the fast states are asymptotically stable, our focus shifts toward designing a Lyapunov-based model predictive control (LMPC) scheme that exclusively employs the RNN to predict the dynamics of the slow states. Additionally, we derive sufficient conditions to guarantee the closed-loop stability of the system under the sample-and-hold implementation of the controller. A nonlinear chemical process example is used to demonstrate the theory. In particular, two RNN models are constructed: one to model the full two-time-scale system and the other to predict solely the slow state vector. Both models are integrated within the LMPC scheme, and we compare their closed-loop performance while assessing the computational time required to execute the LMPC optimization problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Gas Hydrate Saturation Estimation Form Acoustic Log Data in the 2018 Alaska North Slope Hydrate-01 Stratigraphic Test

Completed in December 2018, the Alaska North Slope Hydrate 01 stratigraphic test well provides a wealth of logging-while-drilling (LWD) data for strata to below the base of gas hydrate stability (BGHS). This well is intended to be the first of three wells drilled for a comprehensive long-term gas hydrate production test conducted by the National Energy Technology Laboratory, the Japan Oil, Gas and Metals National Corporation, and the U.S. Geological Survey (USGS). The Hydrate 01 stratigraphic test well confirmed the presence of gas hydrate in two sand reservoirs within the hydrate stability zone, indicating the suitability of this location for a long-term gas hydrate production test.<p>The USGS, using an effective-medium-theory rock-physics approach, has estimated gas hydrate saturations from compressional (P) and shear (S) wave log data acquired in the Hydrate 01 well. We assume that gas hydrate occurs as pore-filling load-bearing material (i.e., part of the grain matrix). For Unit D, approximately 500 feet above the BGHS, both P-wave and S-wave acoustic logs indicate moderate gas hydrate saturations with S-wave results slightly lower than those for P-waves. For the Unit B, located just above the BGHS, we obtain moderate to high gas hydrate saturation estimates from both sonic logs. Our P-wave saturation estimates agree well with results from electrical-resistivity-based estimates, whereas estimates from nuclear magnetic resonance LWD data generally suggest 5 to 10 percent higher saturations; our S-wave results suggest lower saturations. These differences likely indicate complexities in the form of gas hydrate occurrence within the sediment pore space, potentially including differences between hydrate occurrence in Units D and B.</p>

Haines, Seth↗
DOE OSTI2024

Phase diagram of the three-dimensional subsystem toric code

Subsystem quantum error-correcting codes typically involve measuring a sequence of noncommuting parity check operators. They can sometimes exhibit greater fault tolerance than conventional codes, which use commuting checks. However, unlike subspace codes, it is unclear if subsystem codes—in particular their advantages—can be understood in terms of ground-state properties of a physical Hamiltonian. In this paper, we address this question for the three-dimensional subsystem toric code (3D STC), as recently constructed by Kubica and Vasmer [], which exhibits single-shot error correction. Motivated by a conjectured relation between single-shot properties and thermal stability, we study the zero- and finite-temperature phases of an associated noncommuting Hamiltonian. By mapping the Hamiltonian model to a pair of 3D Z 2 gauge theories coupled by a kinetic constraint, we find various phases at zero temperature, all separated by first-order transitions: There are 3D toric code-like phases with deconfined point-like excitations in the bulk, and there are phases with a confined bulk supporting a 2D toric code on the surface when appropriate boundary conditions are chosen. The latter is similar to the surface topological order present in 3D STC. However, the similarities between the single-shot correction in 3D STC and the confined phases are only partial: they share the same sets of degrees of freedom, but they are governed by different dynamical rules. Instead, we argue that the process of single-shot error correction can more suitably be associated with a path (rather than a point) in the zero-temperature phase diagram, a perspective, which inspires alternative measurement sequences enabling single-shot error correction. Moreover, since none of the above-mentioned phases survives at nonzero temperature, the single-shot error-correction property of the code does not imply thermal stability of the associated Hamiltonian phase. Published by the American Physical Society 2024

Li, Yaodong (ORCID:0000000337421944)↗
DOE OSTI2026

A physics informed model for the prediction of surface energies and the shape of metal particles

Surface energies of metal-based systems are important for determining the Wulff-constructed shapes of metal nanoparticles and understanding the stability. We have developed a physics informed method to predict the total energy of metal-based systems across a wide range of configurations. Our method has been tested against density functional theory (DFT) calculations for late transition metals. This method enables on-the-fly surface energy predictions based on the fundamental of materials property and allows for the Wulff construction of metal particles for a random number of elemental atoms and without the need for DFT calculations. By making a division between atoms in the different layers of the model system, we can improve the accuracy of the model, suggesting a dissimilarity between the electronic structure due to an alternating compression and expansion of atomic layers. We find that our model accurately and effectively provides valuable insights into the distribution and stability of nanoparticles.

Mandal, Shyama Charan (ORCID:0000000245882874)↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗
DOE OSTI2022

Stability Analysis of Substituted Cobaltocenium [Bis(cyclopentadienyl)cobalt(III)] Employing Chemistry-Informed Neural Networks

Cobaltocenium derivatives are promising components of the anion exchange membranes due to their excellent thermal and alkaline stability under the operating conditions of a fuel cell. Here we present an efficient modeling approach of assessing the chemical stability of substituted cobaltocenium CoCp 2 + based on the computed electronic structure enhanced by machine learning techniques. Within the aqueous environment, the positive charge of the metal cation is balanced by the hydroxide anion through formation of the CoCp 2 + OH¯ complexes, whose dissociation is studied within the implicit solvent employing density functional theory. The data set of about 118 species based on 42 substituent groups characterized by a range of electron- donating (ED) and electron-withdrawing (EW) properties is constructed and analyzed. Given 12 carefully chosen chemistry-informed descriptors of the complexes and relevant fragments, the stability of the complexes is found to strongly correlate with the energies of the highest occupied and lowest unoccupied molecular orbitals, modulated by a switching function of the Hirshfeld charge. The latter is used as a measure of the electron withdrawing-donating character of the substituents. Based on this observation from the conventional regression analysis, two fully connected, feed-forward neural network (FNN) models with different unit structures, called the chemistry-informed (CINN) and the quadratic (QNN) neural networks, are developed. Both models predict the bond dissociation energies of the cobaltocenium complexes with mean relative errors less than 5.40% and average absolute errors less than 0.94 kcal/mol. The results show the potential of QNN to efficiently capture more complex relationships. Here, the concept of incorporating the domain (chemical) knowledge/insight into the neural network structure paves the way to applications of machine learning techniques with small data sets, ultimately leading to better predictive models compared to the conventional regression analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Fluid and kinetic nonlinearities of near-acoustic plasma waves

Large amplitude, near-acoustic plasma waves have a temperature-dependent harmonic content and nonlinear frequency shifts. On cold plasmas, experiments show that fluid nonlinearities dominate, and the observed harmonic content is well predicted by a recent fluid analysis. In contrast, the measured nonlinear frequency shift is smaller than this theory prediction, suggesting the presence of other nonlinear effects neglected in this analysis. When the plasma is heated, kinetic effects become important, and we observe an increase in both the harmonic content and nonlinear frequency shift. Quantitative experiments on the stability of these large amplitude, near-acoustic plasma waves strongly support the idea that trapped particles reduce the amplitude threshold for instability below fluid models. Furthermore, at low temperatures, the broad characteristics of the parametric instability are determined by the frequency detuning of the pump and daughter wave and the wave-wave coupling strength, surprisingly consistent with simple cold fluid, three-wave theories. However, at higher temperatures, the pump wave becomes unstable at half the threshold pump wave amplitude, with similar exponential growth rates as for a cold plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unveiling the Role of Lithium Iodide in Stabilizing Solid Interfaces in All-Solid-State Li Metal Batteries

A critical challenge in all-solid-state lithium metal batteries (ASSLMBs) is achieving a stable interface between the lithium metal anode and the solid electrolyte. Leveraging its success in Li/I 2 batteries, lithium iodide has garnered significant attentions for its potential to enhance interfacial stability and overall cell performance in ASSLMBs. Here, we elucidate the role of lithium iodide in stabilizing the solid interface in all-solid-state Li metal batteries with a Li argyrodite electrolyte, particularly focusing on its influence on lithium deposition behavior and interfacial evolution. Through in situ optical imaging, we demonstrate more uniform lithium deposition on an iodide-contained argyrodite electrolyte compared to a chloride-based counterpart. Complementary density functional theory calculations attribute improved lithium plating behavior to the enhanced lithiophilicity and better ionic conductivity of lithium iodide at the solid interface, effectively reducing localized current density. In conclusion, these findings provide useful insights into the mechanisms through which lithium iodide enhances the interfacial stability in ASSLMBs.

Anions↗
DOE OSTI2024

Predicting the Slowing of Stellar Differential Rotation by Instability-driven Turbulence

Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.

Astronomy & Astrophysics↗
DOE OSTI2021

Inter-ELM pedestal evolution and the role of edge fluctuations in the C-Mod and DIII-D tokamaks

Edge localized modes (ELMs) lead to a cyclical behaviour of the H-mode pedestal density, temperature and pressure. Substantial evidence exists associating the ELM crash with violation of ideal MHD stability, via current driven kink/peeling modes and/or pressure gradient driven ballooning modes. Recovery from ELM crashes and buildup to the next ELM crash is an active and critical area of edge physics research for ITER projections. In one set of theories, the pedestal pressure profile in H-mode discharges is predicted to be limited by micro-instabilities. One such micro-instability is the kinetic ballooning mode (KBM), which is hypothesized to provide a ‘soft’ limit that regulates the edge transport, thereby restricting the local pressure gradient. In this picture, e.g. in the EPED model, the pedestal width expands slowly until the ideal MHD stability limit is reached. Experiments were recently performed on the C-Mod and DIII-D devices to search for instabilities correlated with the pedestal evolution between ELMs. The results show correlations between the onset of quasi-coherent fluctuations between type I ELMs and a critical temperature gradient. Linear gyrokinetic calculations have shown that these fluctuations have characteristics similar to those expected for KBMs. These results provide additional data toward validation of the EPED model, increasing confidence in ITER projections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Janus Nanocages of Platinum-Group Metals and Their Use as Effective Dual-Electrocatalysts

Janus nanocages with distinctive platinum-group metals on the outer and inner surfaces can naturally catalyze at least two different reactions. In this work, we report a general method based on successive deposition and then selective etching for the facile synthesis of such nanocages. Herein, we have fabricated 11 different types of Janus nanocages characterized by a uniform size and well-defined {100} facets, together with porous, ultrathin, asymmetric walls up to 1.6 nm thick. When tested as dual-electrocatalysts toward oxygen reduction and evolution reactions, the Janus nanocages based on Pt and Ir exhibited superior activities depending on the thickness and relative position of the metal layer. Density functional theory studies suggest that the alloy composition and surface structure of the nanocages both play important roles in enhancing the electrocatalytic activities by modulating the stability of key reaction intermediates.

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

Perspectives on Few-Body Cluster Structures in Exotic Nuclei

It is a fascinating phenomenon in nuclear physics that states with a pronounced few-body structure can emerge from the complex dynamics of many nucleons. Such halo or cluster states often appear near the boundaries of nuclear stability. As such, they are an important part of the experimental program beginning at the Facility for Rare Isotope Beams (FRIB). A concerted effort of theory and experiment is necessary both to analyze experiments involving effective few-body states, as well as to constrain and refine theories of the nuclear force in light of new data from these experiments. Here, as a contribution to exactly this effort, this paper compiles a collection of “perspectives” that emerged out of the Topical Program “Few-body cluster structures in exotic nuclei and their role in FRIB experiments” that was held at FRIB in August 2022 and brought together theorists and experimentalists working on this topic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗