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

A Cooperative Game Theory-based Approach to Under-frequency Load Shedding Control

This paper proposes a cooperative game theory-based under-frequency load shedding (UFLS) approach for frequency stability and control in power systems. UFLS is a crucial factor for frequency stability and control especially in power grids with high penetration of renewable energy sources and restructured power systems. Conventional UFLS methods, most of which are off-line, usually shed fixed amounts of predetermined loads based on a predetermined schedule which can lead to over or under curtailment of load. This paper presents a co-operative game theory-based two-stage strategy to effectively and precisely determine locations and amounts of loads to be shed for UFLS control. In the first stage, the total amount of loads to be shed, also referred to as deficit in generation or the disturbance power, is computed using the initial rate of change of frequency (ROCOF) referred to the equivalent inertial center. In the second stage, the Shapley value, one of the solution concepts of cooperative game theory, is used to determine load shedding amounts and locations. The proposed method is implemented on the reduced 9-bus 3-machine Western Electricity Coordinating Council (WECC) system and simulated on Real-time Digital Simulators (RTDS). The results show that the proposed UFLS approach can effectively return the system to normal state after disturbances.

42 ENGINEERING↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase Stability, Band Gap Tuning, and Rashba Splitting in Selenium-Alloyed Bournonite: CuPbSb(S 1– x Se x ) 3

Recently, bournonite (CuPbSbS 3 ) has been identified as a potential ferroelectric photovoltaic (PV) material with a 1.3 eV band gap, which falls in an appropriate range for single junction PV devices. Progress in applying bournonite has yielded several studies reporting successful thin-film processing, the most recent of which demonstrated a 2.65% power conversion efficiency for a bournonite-based PV device. In an effort to explore bournonite band gap engineering for PV and other prospective applications (e.g., thermoelectric, spintronic), we here consider the solid-state synthesis of selenium-alloyed bournonite, CuPbSb(S 1-x Se x ) 3 , across the full range of x (0.0 ≤ x ≤ 1.0), and report phase purity for 0.0 ≤ x ≤ 0.5. We characterize the crystal structure and band gap of the samples using X-ray diffraction and diffuse reflectance spectroscopy, determining a band gap decrease for single-phase samples from 1.25 eV at x = 0.0 to 1.06 eV at x = 0.5. Formation energies determined using dispersion corrected hybrid density functional theory (DFT) support experimental findings and are consistent with the stability (instability) of the x = 0.5 (x = 1.0) structure relative to starting materials. The computed band structures show a decreasing trend in the band gap with increasing x, again consistent with experiment, with the states at the conduction (valence) band minima (maxima) being principally derived from Pb (S/Se) states. Spin texture analysis and detailed comparison of spin splitting parameters show that Se-alloying modifies the band gap while maintaining Rashba spin splitting character within the electronic structure. Rashba splitting is most pronounced for the conduction band minimum along the Γ-X path in the Brillouin zone. Energy separation between spin states is maintained at ΔE~0.2 eV for the various Se contents, with σz as the spin direction. Furthermore, a smaller spin splitting occurs along Γ-Z (decreasing value with increasing x), with σ x as the spin direction.

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

Effect of Impurities on the Redox Properties of Goethite

Iron oxide minerals regulate the flux of electrons in the environment and are important hosts for trace and minor, yet critical, elements. Here, we present the first evidence of a direct link between the local coordination environments of Ni and Zn and the redox properties of their host phase goethite (α-FeOOH), the most abundant Fe(III) (oxyhydr)oxide at Earth’s surface. Here, we used aqueous redox measurements to show that the redox potential E H 0 , and hence the mineral’s stability, follows the order: pure goethite ≥ Zn-goethite > Ni-goethite. Parallel X-ray absorption and scattering measurements demonstrate, using quantum-informed analysis, that the local coordination environment of the smaller impurity, Ni, causes more bulk strain energy than Zn, which nearly accounts for the difference in E H 0 between Ni- and Zn-goethite. Our theory-informed, experimental study reveals how two common impurities affect the stability of goethite with implications for the biogeochemical reactivity of Fe(III) (oxyhydr)oxide in mediating elemental and electron fluxes in the environment.

42 ENGINEERING↗
DOE OSTI2025

Thermodynamic Stability and Site‐Specific Distribution of Graphitic and Pyridinic Nitrogen in Graphene Moiré on Ru(0001)

Abstract Graphene‐like materials are of interest for large‐scale hydrogen storage applications due to their lightweight, durable, and scalable properties. Nitrogen‐doping minimizes kinetic limitations in diffusion and recombination on surfaces, however, the role of graphitic nitrogen (GN) and pyridinic nitrogen (PN) is not well understood. Nitrogen‐doped graphene is synthesized on Ru(0001) using chemical vapor deposition (CVD) of pyridine and ion irradiation. Scanning tunneling microscopy (STM), x‐ray photoelectron spectroscopy (XPS), and density functional theory (DFT) are used to identify the structure, location, and thermodynamic stability of nitrogen species within the graphene moiré. CVD of pyridine results in a low nitrogen concentration (<0.1at%), while the post‐growth nitrogen ion irradiation allows us to increase the concentration further. The concentration of GN and PN is controlled by varying the ion dose and annealing temperature. Comparison of measured and simulated STM images of GN and PN yield an excellent agreement, allowing us to confidently establish that GN is preferentially located near the center of the Atop region, while PN is located in the valley region of the graphene moiré. This report explicitly confirms the site assignments and provides a foundation for the site synthesis and analysis of structural and electronic properties that drive the reactivity of N‐doped graphene.

Gedara, Buddhika S. A. [Physical and Computational↗
DOE OSTI2023

Small bimetallic clusters Agn-1M (M = Au, Co, Cu, Ni, Pd, Pt; n = 3, 9, 15): Density functional theory and genetic algorithm

We investigated the effect of size and composition on the properties of bimetallic nanoclusters. The geometric structures, stabilities, and electronic properties of size-selected Ag n-1 M (M = Au, Co, Cu, Ni, Pd, Pt; n = 3, 9, 15) bimetallic nanoclusters are systematically analyzed using spin-polarized density functional theory (DFT) within the generalized gradient approximation (GGA). We determine the most stable geometries for these clusters using a genetic algorithm (GA) in combination with DFT. Our results show that doping pure silver clusters with an M atom (transition metal), referred to as a “guest atom”, increases the stability as compared to pure Ag n (n = 3, 9, 15) clusters. The results for various properties including formation energy per atom, electronic structure, magnetic moments, and vibrational density of states (VDOS) are evaluated as a function of both size and composition of the system. The adsorption of selected bimetallic clusters on hydroxylated alumina substrate shows weak binding and minor changes in geometric properties except for Ag 8 Pt.

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

Photodissociation and Theory to Investigate Uranium Oxide Cluster Cations

Uranium oxide cluster cations of the form U n O m + are produced by laser vaporization of a depleted uranium rod in a pulsed supersonic expansion. Ions are mass-analyzed and mass-selected with a time-of-flight spectrometer and studied with UV laser multiphoton dissociation. Cations of the stoichiometry UO 2 (UO 3 ) n + were observed as photofragments from all photodissociated cluster cations. (UO 3 ) n + clusters were also observed to result from dissociation of larger (UO 3 ) n + clusters, with UO 3 neutral as a common leaving group. Electronic structure calculations were used to investigate the stability of the prominent uranium oxide cluster cations using density functional theory (DFT) with the hybrid B3LYP exchange-correlation functional and at the CCSD(T) level and cc-pVnZ-PP basis sets (n = D,T), including diffuse orbitals as computational expense and availability permitted. Clustering energies, relative energies and dissociation energies of the cations are reported. The lowest energy neutral (UO 3 ) n clusters up to n = 3 are rings, n = 4 and 5 are chains with very low energy rings, and n = 6 is 3D. Finally, the lowest energy structures for UO 2 (UO 3 ) n + are composed of uranyl-like UO 2 + units bound by bridging oxygens to other UO 2 2+ units for n = 2 and 3, and for n = 4 a more complex 3D structure is predicted.

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

The role of excess vacancies in stabilizing solute clusters in low-alloy steels

The formation of solute clusters in irradiated low-alloy steels, such as reactor pressure vessel steels, is a critical cause of radiation hardening and embrittlement. However, the chemical interactions among solute elements and excess vacancies remain to be fully understood. This study employs density functional theory, cluster expansion, and lattice-based Monte Carlo simulations to examine the stability and morphology of nano-size coherent solute clusters with excess vacancies. The findings reveal that excess vacancies are crucial in stabilizing and promoting the growth of Mn-Ni-Si-vacancy clusters. A minimum of seven vacancies is required for stable nucleation and growth of these clusters. These Mn-Ni-Si clusters act as defect sinks, effectively trapping and absorbing mobile vacancies generated under irradiation. Furthermore, phosphorus (P) preferentially dissolves in Mn-Ni-Si clusters due to chemical coupling with vacancies. In conclusion, this study offers new insights into solute-vacancy interactions, enhancing our understanding of solute-defect cluster formation and embrittlement in low-alloy steels.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Trends in electrocatalytic activity and stability of transition-metal nitrides

Transition metal nitrides (TMNs) are a class of electrocatalyst support materials similar to transition metal carbides (TMCs) with the advantage of avoiding the issues arising from graphitic carbon surface deposits during synthesis. Inspired by previous studies suggesting that TMCs could be used to reduce Pt loading for the hydrogen evolution reaction (HER), this work explored the feasibility of TMN-supported Pt and Au as HER electrocatalysts. Here this study established a volcano-like trend between electrochemical HER activity and hydrogen-binding energy (HBE) calculated from density functional theory. The Pt/TiN and Au/TiN materials were used to extend knowledge from well-characterized thin films to powder catalysts. In situ X-ray absorption spectroscopy (XAS) measurements provided additional characterization of the Pt/TiN and Au/TiN catalysts under HER conditions. Trends in the electrochemical stability of TMNs were also investigated over a wide range of potentials and pH values, which can be used to guide future studies for TMN-supported electrocatalysts.

58 GEOSCIENCES↗
DOE OSTI2021

First-principles calculations for understanding microstructures and mechanical properties of co-sputtered Al alloys

Recent experimental studies show that co-sputtering solutes with Al, together, can refine columnar grain size around few tens of nanometers and promote the formation and enhance the stability of planar defects such as stacking faults (SFs) and grain boundaries (GBs) in Al alloys. These crystal defects and fine columnar grains result in high strength, enhanced strain hardening and thermal stability of Al alloys. Using first-principles density-functional theory (DFT) calculations, we studied the role of eleven solutes in tailoring kinetics and energetics of adatoms and clusters on Al {111} surface, stable and unstable stacking fault energies, and kinetic energy barriers for the migration of defects. The calculations show that most solutes can effectively refine columnar grain size by decreasing the diffusivity of adatoms and surface clusters. These solutes do not necessarily decrease the stacking fault energy of Al alloys, but reduce the formation energy of faulted surface clusters and increase the energy barriers for the recovery of faulted surface clusters. Correspondingly, the formation of SFs is kinetically promoted during sputtering. Furthermore, solutes are segregated into the core of Shockley partial dislocations and play a pinning effect on SFs, SF arrays and twin boundaries, enhancing the thermal stability of these crystal defects. These findings provide insights into the design of high-strength Al alloys for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗
DOE OSTI2026

Synthesis of epitaxial TaO$_2$ thin films on Al$_2$O$_3$ by suboxide molecular-beam epitaxy and thermal laser epitaxy

Tantalum dioxide (TaO2) is a metastable tantalum compound. Here, we report the epitaxial stabilization of TaO2 on Al2O3 (1-102) (r-plane sapphire) substrates using suboxide molecular-beam epitaxy (MBE) and thermal laser epitaxy (TLE), demonstrating single-oriented, monodomain growth of anisotropically strained thin films. Microstructural investigation is performed using synchrotron X-ray diffraction and scanning transmission electron microscopy. The tetravalent oxidation state of tantalum is confirmed using X-ray absorption and photoemission spectroscopy as well as electron energy-loss spectroscopy. Optical properties are investigated via spectroscopic ellipsometry and reveal a 0.3 eV Mott gap of the tantalum 5d electrons. Density-functional theory and group theoretical arguments are used to evaluate the limited stability of the rutile phase and reveal the potential to unlock a hidden metal-insulator transition concomitant with a structural phase transition to a distorted rutile phase, akin to NbO2. Our work expands the understanding of tantalum oxides and paves the way for their integration into next-generation electronic and photonic devices.

FOS: Physical sciences↗
DOE OSTI2025

Computer-aided design of stability enhanced nicotinamide cofactor biomimetics for cell-free biocatalysis

Cell-free biocatalysis (CFB) is an efficient and environmentally friendly method to synthesize molecules such as pharmaceuticals, biochemicals, and biofuels through the in vitro use of enzyme cascades. These enzymes often require redox cofactors to drive chemical reactions. Natural redox cofactors (NAD(P)H) are expensive to isolate, motivating synthetic nicotinamide cofactor biomimetics (NCBs) as a cost-effective solution. A select handful of NCBs have been identified as potential NAD(P)H alternatives with comparable or improved redox capabilities, however, they display a tendency to degrade in common buffers. In this study, a library of 132 NCB candidates is systematically generated, over 85% of which have not been characterized in the literature, to expand the diversity of currently explored NCBs. The decomposition mechanism of NCBs in phosphate is evaluated using density functional theory (DFT), revealing protonation at the nicotinamide C5 position as a reporter of cofactor stability. Based on this result, we trained a linear regression model on DFT calculated descriptors to predict NCB stability in phosphate buffer, achieving mean absolute error (MAE) and root mean squared error (RMSE) values within computational accuracy. Analysis of key atomic descriptors and qualitative trends in our dataset informed the design of novel NCB candidates we propose with optimized stability. This work enables researchers to predict the relative stability of NCBs before synthesis, thereby streamlining the process to make CFB more affordable and viable at industry scales.

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

Metal Identity and Coordination Environment Modulate Single-Atom Catalyst Stability During Electrocatalysis

A major hurdle to the implementation of single-atom catalysts (SACs) in real-world systems is a poor understanding of their stability under operating conditions, which is particularly relevant due to the high surface free energy of SACs. Here, we evaluated the aggregation behavior of a suite of SACs varied by metal identity (Fe, Co, Ni, and Cu) during electrocatalytic nitrate reduction using in situ X-ray absorption spectroscopy. The metal center had significant influence on reconstruction, where under identical applied reductive potentials, SACs underwent varying levels of reconstruction, ranging from no discernible change to complete reduction into metallic nanoparticles. Such in situ experiments revealed Cu SACs to be the most susceptible to aggregation, prompting a deeper investigation into how coordination environment (O-, B-, and N-graphene) affected Cu SAC aggregation. We further conducted density functional theory calculations to elucidate the relationship between Cu SAC structure and stability. This work deconvolutes the relationship between SAC architecture and stability, which is essential to evaluate and explain for the realization of SACs for electrocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulated annealing of reduced magnetohydrodynamic systems

Theory of simulated annealing (SA), a method for equilibrium and stability analyses for Hamiltonian systems, is reviewed. The SA explained in this review is based on a double bracket formulation that derives from Hamiltonian structure. In addition to general theoretical aspects, the explicit formulation as well as numerical applications are presented. Both finite and infinite degree-of-freedom systems are treated, in particular, the heavy top, a toy model mimicking low-beta reduced magnetohydrodynamics (MHD) and low- and high-beta reduced MHD. Furthermore, the numerical results successfully demonstrate the usefulness of SA for equilibrium and stability analyses. At the same time, the results raise some future issues that are discussed in the paper.

Poisson Bracket↗
DOE OSTI2025

On the distribution of stacking faults at dissociated medium-angle grain boundaries: Crystallographic geometry and metastability

Grain boundaries in FCC metals with low stacking-fault energy can form in dissociated configurations of stacking faults. Perhaps the most studied example is 9R stacking at boundaries near Σ3{112}, where the distribution of stacking faults is related to the emission of Shockley partial dislocations. Here, we combine atomic-resolution electron microscopy, atomistic simulations, and dislocation theory to demonstrate that boundaries vicinal to Σ33a support the stabilization of dissociated 9R stacking within a narrow range of inclinations. This boundary is interesting since its misorientation (20.05°) lies in the medium-angle regime, just past the upper misorientation limit for low-angle boundaries, motivating questions for how best to describe it in terms of dislocations. Our HAADF-STEM observations of thin film bicrystals, supported by atomistic modeling, reveal that this inclination dependence arises from specific geometric constraints on the arrangement of Shockley partial dislocations at the interface. Quantification of stacking-fault distributions across multiple boundaries indicates that the density and spacing of faults closely follow the ideal 9R motif, with subtle variations reflecting the complex energy landscape of these boundaries. Through energetic analysis, we establish the presence of competing metastable states enabled by variations in stacking sequences, emphasizing the significant role of crystallographic geometry. We generalize our analysis as a function of misorientation, showing how 9R at the Σ33a boundary is related to previous observations and calculations of HCP at a 29.7° boundary. This study provides a crystallographically grounded framework connecting dislocation structures, stacking-fault distributions, and metastability at grain boundaries in FCC metals.

Atomistic modeling↗
DOE OSTI2026

High-throughput validation of phase formability and simulation accuracy of Cantor alloys

High-throughput methods enable accelerated discovery of novel materials in complex systems such as high-entropy alloys, which exhibit intricate phase stability across vast compositional spaces. Computational approaches, including Density Functional Theory (DFT) and calculation of phase diagrams (CALPHAD), facilitate screening of phase formability as a function of composition and temperature. However, the integration of computational predictions with experimental validation remains challenging in high-throughput studies. In this work, we introduce a quantitative confidence metric to assess the agreement between predictions and experimental observations, providing a quantitative measure of the confidence of machine learning models trained on either DFT or CALPHAD input in accounting for experimental evidence. The experimental dataset was generated via high-throughput in-situ synchrotron X-ray diffraction on compositionally varied FeNiMnCr alloy libraries, heated from room temperature to ~1000 °C. Agreement between the observed and predicted phases was evaluated using either temperature-independent phase classification or a model that incorporates a temperature-dependent probability of phase formation. This integrated approach demonstrates where strong overall agreement between computation and experiment exists, while also identifying key discrepancies, particularly in FCC/BCC predictions at Mn-rich regions to inform future model refinement.

36 - MATERIALS SCIENCE↗