Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Stability”

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

Molecular Dynamics Simulations of Calcite Fracture in Water

Calcite (CaCO 3 ) is one of the most common minerals in geologic and engineered systems. It is often in contact with aqueous solutions, causing chemically assisted fracture that is critical to understanding the stability of subsurface systems and manmade structures. Calcite fracture was evaluated with reactive molecular dynamics simulations, including the impacts of crack tip geometry (notch), the presence of water, and surface hydroxyl groups. Chemo-mechanical weakening was assessed by comparing the loads where fracture began to propagate. Our analyses show that in the presence of a notch, the load at which crack growth begins is lower, compared to the effect of water or surface hydroxyls. Additionally, the breaking of two adjacent Ca–O bonds is the kinetic limitation for crack initiation, since transiently broken bonds can reform, not resulting in crack growth. In aqueous environments, fresh (not hydroxylated) calcite surfaces exhibited water strengthening. Manual addition of H + and/or OH – species on the (104) calcite surface resulted in chemo-mechanical weakening of calcite by 9%. Achieving full hydroxylation of the calcite surface was thermodynamically and kinetically limited, with only 0.17–0.01 OH/nm 2 surface hydroxylation observed on the (104) surface at the end of the simulations. In conclusion, the limited reactivity of pure water with the calcite surface restricts the chemo-mechanical effects and suggests that reactions between physiosorbed water and localized structural defects may be dominating the chemo-mechanical process in the studies where water weakening has been reported.

Calcite↗
DOE OSTI2023

Stability and Fidelity of Power-Hardware-in-the-Loop Using NREL's CGI

The presentation focuses on dynamic analysis of power hardware-in-the-loop (PHIL) experiments using National Renewable Energy Laboratory Flatirons Campus' grid simulator - 7MVA Controllable Grid Interface (CGI). A control diagram of typical PHIL experiment is described, and transfer functions of critical components are explained. An impedance-based definition of PHIL experiment bandwidth is defined using open loop transfer functions of the system. Further an example resistive divider circuit is introduced that is used to explain implication of PHIL on accuracy of simulation by comparing measured results against ideal expected signals. At the end the real-life PHIL simulation using Maui power system model is described showing various oscillation modes observed and discussing the meaning of these in context of dynamic bandwidth of simulation.

grid simulator↗
DOE OSTI2025

Optimizers for stabilizing likelihood-free inference

A growing number of applications in particle physics and beyond use neural networks as unbinned likelihood ratio estimators applied to real or simulated data. Precision requirements on the inference tasks demand a high-level of stability from these networks, which are affected by the stochastic nature of training. We show how physics concepts can be used to stabilize network training through a physics-inspired optimizer. In particular, the energy conserving descent (ECD) optimization framework uses classical Hamiltonian dynamics on the space of network parameters to reduce the dependence on the initial conditions while also stabilizing the result near the minimum of the loss function. We develop a version of this optimizer known as , which has few free hyperparameters with limited ranges guided by physical reasoning. We apply to representative likelihood-ratio estimation tasks in particle physics and find on average that it out-performs the widely used Adam optimizer. We expect that ECD will be a useful tool for wide array of data-limited problems, where it is computationally expensive to exhaustively optimize hyperparameters and mitigate fluctuations with ensembling.

Monte Carlo methods↗
DOE OSTI2021

Determination of Thermal Expansion, Defect Formation Energy, and Defect-Induced Strain of α-U Via ab Initio Molecular Dynamics

Uranium (U) is often alloyed with molybdenum (Mo) or zirconium (Zr) in order to stabilize its high-temperature body-centered cubic phase for use in nuclear reactors. However, in all metallic fuel forms, the α phase of U remains in some fraction. This phase decomposition due to temperature or compositional variance can play an outsized role on fuel performance and microstructural evolution. Relatively little is known about fundamental point defect properties in α -U at non-zero temperatures, from either computational or experimental studies. This work performs the first thorough evaluation of the α phase of U via ab initio molecular dynamics (AIMD). A number of thermophysical properties are calculated as a function of temperature, including equilibrium lattice parameters, thermal expansion, and heat capacity. These results indicate a two-region behavior, with the transition at 400 K. The thermal expansion/contraction in the a / b direction occurs rapidly from 100 up to 400 K, after which a more linear and gradual change in the lattice constant takes place. The volumetric expansion matches experiments quantitatively, but the individual lattice constant expansion only matches experiments qualitatively. Point defect formation energies and induced lattice strains are also determined as a function of temperature, providing insight on defect populations and the fundamentals of irradiation growth in α -U. Interstitials induce significantly more strain than vacancies, and the nature of that strain is highly dependent on the individual lattice directions. The direction of point defect-induced lattice strain is contrary to the irradiation growth behavior of α -U. This work shows that isolated point defects cannot be the primary driving force responsible for the significant irradiation-induced growth of α -U observed experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Atomic-scale 3D structural dynamics and functional degradation of Pt alloy nanocatalysts during the oxygen reduction reaction

Pt-based electrocatalysts are the primary choice for fuel cells due to their superior oxygen reduction reaction (ORR) activity. To enhance ORR performance and durability, extensive studies have investigated transition metal alloying, doping, and shape control to optimize the three key governing factors for ORR: geometry, local chemistry, and strain of their surface and subsurface. However, systematic optimization remains incomplete, as it requires an atomic-scale understanding of these factors and their dynamics over potential cycling, as well as their relationship to ORR activity. Here, we implement neural network-assisted atomic electron tomography to measure the 3D atomic structural dynamics and their effects on the functional degradation of PtNi alloy catalysts. Our results reveal that PtNi catalysts undergo shape changes, surface alloying, and strain relaxation during cycling, which can be effectively mitigated by Ga doping. By combining geometry, local chemistry, and strain analysis, we calculated the changes in ORR activity over thousands of cycles and observed that Ga doping leads to higher initial activity and greater stability. These findings offer a pathway to understanding 3D atomic structural dynamics and their relation to ORR activity during cycling, paving the way for the systematic design of durable, high-efficiency nanocatalysts.

Jeong, Chaehwa↗
DOE OSTI2026

Influence of the drive design on ablation front hydrodynamic instability growth in a capsule implosion at the National Ignition Facility

A critical aspect of inertial confinement fusion experiments lies in the control of ablation front instabilities during the implosion process. The growth of Rayleigh–Taylor and Richtmyer–Meshkov instabilities, seeded by target defects, can significantly degrade the performance of the implosion. Here, this study explores the influence of drive design on ablation front stability using the hydro-growth radiography platform at the National Ignition Facility. This platform allows the measurement of the ablation front hydrodynamic instability growth. Because these measurements are restricted to early convergence, their primary role is to constrain the initial instability growth and benchmark the simulations used to infer shell integrity and performance at peak velocity and ignition. Three ignition designs, Hybrid-E (HyE), SQ-n (“S” for scaling and “Q” for quality), and High temperature, High thickness (HiT) were analyzed. The results demonstrated a high dependence of the growth factor on the choice of design, capsule scale as well as the hohlraum conditions, mainly the radiative temperature and the gold M-band emission. HyE and SQ-n display a similar growth factor range, but their evolution is different due to the design differences in hohlraum conditions. HiT has the lowest growth factor of the studied designs, which is a result of a higher radiative temperature. These results highlight the importance of design choices in controlling instability dynamics. These insights inform future drive design strategies to enhance stability and efficiency in fusion ignition experiments.

Physics↗
DOE OSTI2020

Improved Single-Ion Conductivity of Polymer Electrolyte via Accelerated Segmental Dynamics

Single-ion conducting polymer electrolytes (SICPEs) have many advantages for solid-state battery applications, while low conductivities remain the bottleneck for their practical applications. In this paper, a strategy that can significantly improve the ionic conductivity of SICPEs based on the concept of accelerated segmental dynamics was demonstrated by the evaluation of the thermal property, rheological behavior, morphology analysis, molecular dynamics simulation, and conductivity performance. With both “soft” poly(dimethylsiloxane) (PDMS) backbone and poly(ethylene glycol) (PEG) side chains, the obtained SICPE possesses faster segmental relaxation and higher lithium-ion conductivity. With lithium “transference” number close to unity, the obtained SICPE also shows excellent electrochemical stability against lithium metal electrodes. The clear relationship established between the segmental dynamics of polymer electrolyte and its ionic conductivity should contribute to achieve a solid electrolyte with improved ionic conductivity toward the next-generation solid-state battery.

Segmental Dynamics↗
DOE OSTI2025

Mixed lipid bilayers enable enhanced stability and activity retention of Lipase A in low-pH environments

Lipase A (LipA) from Bacillus subtilis is a versatile and industrially relevant enzyme, but its activity is compromised under acidic conditions due to aggregation and deactivation. In this study, we investigated the use of mixed lipid bilayers to stabilize and improve activity retention of LipA at low pH. Attachment to lipid bilayers, particularly those containing cationic lipids, greatly enhanced the long-term stability of LipA in acidic conditions while also leading to improvements in activity. Tethering to cationic bilayers not only shifted the apparent pH activity profile toward more acidic conditions, but also significantly enhanced activity retention upon incubation at pH 6. Notably, this protective effect persisted even without direct tethering, indicating that reversible, non-covalent interactions with the bilayer surface are sufficient for long-term stability. Circular dichroism further revealed that the secondary structure of LipA was retained, while dynamic light scattering suggested that activity loss in solution was primarily due to aggregation of the native state. Together, these findings support a model in which lipid bilayers mitigate aggregation and stabilize LipA through transient interactions and electrostatic modulation of the local surface environment. Furthermore, this approach exemplifies a simple, tether-free strategy for enhancing enzyme performance in acidic or destabilizing conditions, with implications for biocatalysis, biosensing, and therapeutic delivery.

Biocatalysis↗
DOE OSTI2021

Borehole breakout modeling in arkose and granite rocks

The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.

58 GEOSCIENCES↗
DOE OSTI2025

Conformationally Adaptable Extractant Flexes Strong Lanthanide Reverse-Size Selectivity

Chemical selectivity is traditionally understood in the context of rigid molecular scaffolds with precisely defined local coordination and chemical environments that ultimately facilitate a given transformation of interest. By contrast, nature leverages dynamic structures and strong coupling to enable specific interactions with target species in otherwise complex media. Taking inspiration from nature, we demonstrate unconventional selectivity in the solvent extraction of light over heavy lanthanides using a conformationally flexible ligand called octadecyl acyclopa (ODA). This novel ligand forms pseudocyclic molecular complexes with lanthanide ions at organic/aqueous interfaces, revealed by vibrational sum frequency generation spectroscopy. These complexes are extracted into the organic phase, where femtosecond structural dynamics are probed by two-dimensional infrared spectroscopy and ab initio molecular dynamics simulations to mechanistically frame the macroscopic selectivity trends. We find larger-than-expected structural fluctuations and bond lengths for heavy Ln–ODA complexes that arise from an inability of ODA to contort around the smaller ions to satisfy all would-be bonding interactions, despite forming some individually strong bonds. This finding contrasts with the binding of ODA with lighter lanthanides where, despite individually weaker bonds, collective interactions manifest that minimize structural fluctuations and give rise to enhanced thermodynamic stability. Furthermore, these results point to a new paradigm where conformational dynamics and cumulative bonding interactions can be used to facilitate unconventional chemical transformations.

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

Exploring the High-Pressure Phases of Carbon through X-ray Diffraction of Dynamic Compression Experiments on Sandia’s Z Pulsed Power Facility

The carbon phase diagram is rich with polymorphs which possess very different physical and optical properties ideal for different scientific and engineering applications. An understanding of the dynamically driven phase transitions in carbon is particularly important for applications in inertial confinement fusion, as well as planetary and meteorite impact histories. Experiments on the Z Pulsed Power Facility at Sandia National Laboratories generate dynamically compressed high-pressure states of matter with exceptional uniformity, duration, and size that are ideal for investigations of fundamental material properties. X-ray diffraction (XRD) is an important material physics measurement because it enables direct observation of the strain and compression of the crystal lattice, and it enables the detection and identification of phase transitions. Several unique challenges of dynamic compression experiments on Z prevent using XRD systems typically utilized at other dynamic compression facilities, so novel XRD diagnostics have been designed and implemented. We performed experiments on Z to shock compress carbon (pyrolytic graphite) samples to pressures of 150–320 GPa. The Z-Beamlet Laser generated Mn-Heα (6.2 keV) X-rays to probe the shock-compressed carbon sample, and the new XRD diagnostics measured changes in the diffraction pattern as the carbon transformed into its high-pressure phases. Quantitative analysis of the dynamic XRD patterns in combination with continuum velocimetry information constrained the stability fields and melting of high-pressure carbon polymorphs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spontaneous formation of robust two-dimensional perovskite phases

The two-dimensional on three-dimensional (2D/3D) perovskite bilayer heterostructure can improve the stability and performance of perovskite solar cells. We show that the 2D/3D perovskite stack in a device evolves dynamically during its end-of-life decomposition. Initially phase-pure 2D interlayers can evolve differently, resulting in different device stabilities. We show that a robust 2D interlayer can be formed using mixed solvents to regulate its crystallinity and phase purity. The resulting 2D/3D devices achieved 25.9% efficiency and had good durability, retaining 91% of their initial performance after 1074 hours at 85°C using maximum power point tracking.

14 SOLAR ENERGY↗
DOE OSTI2024

Emergent symmetry in TbTe 3 revealed by ultrafast reflectivity under anisotropic strain

Here, we report ultrafast reflectivity measurements of the dynamics of the order parameter of the charge density wave (CDW) in TbTe 3 under anisotropic strain. We observe an increase in the frequency of the amplitude mode with increasing tensile strain along the a-axis (which drives the lattice into a > c, with a and c the lattice constants), and similar behavior for tensile strain along c (c > a). This suggests that both strains stabilize the corresponding CDW order and further support the near equivalence of the CDW phases oriented in a- and c-axis, in spite of the orthorhombic space group. The results were analyzed within the time-dependent Ginzburg–Landau framework, which agrees well with the reflectivity dynamics. Our study presents an ultrafast approach to assess the stability of phases and order parameter dynamics in strained systems.

Kim, Soyeun↗
DOE OSTI2025

Unravelling the dynamics of the maturation protein in MS2 bacteriophage via molecular simulations

The MS2 bacteriophage capsid serves as a model system for studying viral structure and function. Mature MS2 virus consists of 178 capsid proteins and a single maturation protein (MP), which is essential for host receptor binding and infection initiation. Despite its critical role, the dynamic behavior of the capsid with the MP remains poorly understood. To address this, we conducted 0.5 µs all-atom molecular dynamics (MD) simulations of the MS2 capsid with and without the MP, revealing key insights into its structural dynamics. Our simulations showed that MP exhibits high flexibility, particularly in the “tip” and “side-loop” regions, which undergo significant motions that likely enhance its ability to engage with the F-pilus receptor. Detailed analysis of MP conformational states revealed that loop rearrangements around H357 enable transient switching between “semi-closed” and “open” conformations, suggesting a conformational selection mechanism for pilus binding. Additionally, ion interaction analyses revealed distinct sodium and chloride binding patterns, where sodium ions were mostly found at the outer capsid shell, while chloride ions interacted with the basic residues on the RNA-facing side. We also found that the presence of the MP enhances salt-bridge interactions, contributing to increased capsid stability, yet it does not significantly alter the pore sizes of pentameric and hexameric units. Together, these findings provide new insights into the functional role of the MP, highlighting its contribution to capsid stability and host receptor engagement. This study offers a foundation for understanding capsid dynamics relevant to viral infectivity and may guide future rational strategies aimed at disrupting host-virus interactions.

Capsid stability↗
DOE OSTI2025

Compatibilization Strategy and Mechanism for Co-stabilizing Commingled Plastics and Pyrolyzed Rubber in Asphalt

Hot mix asphalt mixture is considered the ideal approach to reuse waste plastics in high-value applications because of its very high amount of usage in highway construction. However, the differences in polarity and density between polymers and asphalt lead to polymer coalescence and therefore the poor storage stability of modified asphalt. These challenges are exalted when recycling commingled plastics. This study introduced an innovative compatibilization strategy and mechanism for co-stabilizing commingled plastics and pyrolyzed rubber in asphalt. Commingled plastics were first grafted with maleic anhydride for surface activation, followed by reactive kneading with pyrolyzed rubber and crosslinking agent to form an integrated thermoplastic elastomer (ITPE) for asphalt modification. The mechanical, thermal, and interfacial behaviors of the ITPE were evaluated through tensile testing, thermogravimetric analysis, and scanning electron microscopy. The storage stability and rheological properties of the modified binder blends were evaluated through the cigar tube test and dynamic shear rheometer testing. Results demonstrated a successful formation of imide bonds in the ITPE, which can improve the strength, ductility, and thermal stability of rubber–plastic composites. Appropriate utilization of crosslinking agents can improve both rutting and fatigue resistance of ITPE-modified asphalt with good storage stability because of the co-existence of rigid plastic and soft rubbery regimes and the formation of a crosslink network. Furthermore, excessive content of crosslinker led to severe phase separation and reduced storage stability of modified binder blends. Extra crosslinker tended to float in asphalt because of its low density and caused an excessive formation of the crosslink network in the top section of the asphalt.

asphalt binder modifiers↗
DOE OSTI2024

Impedance Scan of Inverter-Based Resources and Diesel Generator for Stability Analysis

Impedance-based methods are widely used for power system stability analysis with inverter-based resources (IBRs), e.g., assessing dynamic interactions between the power grid and an IBR, control interactions between multiple IBRs, and the sub-synchronous oscillation and damping phenomenon. Since it is difficult to get a numerical model 100% matching with the hardware IBR, using the hardware inverter directly to obtain its output impedance has become a prominent approach nowadays. Therefore, this article presents the impedance scan using hardware IBRs, and also a hardware diesel generator as it still stays with the grid before the grid completely goes to renewable. The devices under test (DuTs) for the impedance scan includes two 3-..phi.., 480 V, 60 Hz commercial grid-forming IBRs (one of 250 kVA and another of 125 kVA rating) in series with ..delta..-Y transformers, one 3-..phi.., 480 V, 60 Hz commercial grid-following IBR (of 125 kVA rating), and a 3-..phi.., 480 V, 60 Hz commercial diesel generator (of 187.5 kVA rating). Using voltage signals perturbed with sub-, inter-, and higher harmonic components, and measuring the current response, the positive-sequence impedances are computed via an offline-based post-analysis. Moreover, best-fit transfer functions are estimated that closely resemble the measured data points of the positive-sequence impedances. Based on the observations from various outcomes of the hardware experiments, this article also provides some fundamental insights on the equivalent positive-sequence impedance of a combination of multiple hardware components by comparing the estimated and the empirically computed impedances. A comparative insight on the damping capability of the DuTs using the positive-sequence impedances of the hardware is also discussed.

current measurement↗
DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗