Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamical evolution and 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

In situ probing of structure and deagglomeration of SnO 2 colloids via small-angle X-ray scattering

Transforming dry nanopowders into stable colloidal dispersions remains challenging due to the cohesive forces between the nanoparticles (NPs) which promote agglomeration. Effective dispersion and deagglomeration of these agglomerates is a critical process in the formulation and preparation of nanoparticle-based functional materials via the colloidal route. Understanding the deagglomeration dynamics provides information to improve microstructural quality in various applications by enabling engineering of agglomerate size, structure and morphology. However, the deagglomeration process dynamics with respect to the evolution of the fractal agglomerate structures, particularly for very small NPs, is still poorly understood and requires further investigation. This study employs in situ small-angle X-ray scattering (SAXS) to investigate the sonication-induced deagglomeration of SnO 2 NPs. Electrostatically stabilized SnO 2 colloids with varying primary particle size (6–21 nm) are investigated in a specifically designed in situ cell using synchrotron-based SAXS to study the influence of sonication time and intensity on the nanoscaled agglomerates. Complete structural analysis via SAXS reveals a direct correlation between changes in agglomerate size and structure, size-dependent deagglomeration behavior and a dependence on the overall energy introduced during sonication into the dispersion regardless of the actual power as well as ultrasonic process parameters in case of SnO 2 NPs. The results suggest that control of the dispersion process during ultrasonic deagglomeration results in tailoring agglomerates with respect to size and structure.

Deagglomeration↗
DOE OSTI2025

In-situ irradiation of uranium carbide

Uranium carbide (UC) is a leading candidate fuel for Generation IV reactors due to its high uranium density and thermal conductivity. However, its irradiation performance—particularly gas bubble swelling and defect dynamics—remains poorly characterized. Using in-situ transmission electron microscopy (TEM), we irradiated UC with 300 keV Xe + and 1 MeV Kr 2+ ions at temperatures up to 900 °C to quantify swelling behavior and dislocation loop evolution. The swelling remained below 0.6 % across all temperatures, suggesting the dimensional stability of UC under irradiation at these temperatures. Dislocation loops grew faster in UC than in UO 2 or UN, correlating with its lower homologous temperature. Notably, nanograin structures emerged in thin regions of the lamellar, mirroring phenomena previously observed in UO 2 and ZrC. These results address critical knowledge gaps in the radiation tolerance of UC and provide insight into its suitability for advanced reactor systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Time-dependent experimental identification of inter-ELM microtearing modes in the tokamak edge on DIII-D

In a series of discharges on the DIII-D tokamak, fast vertical plasma jogs are used to induce current perturbations in the steep gradient region of the H-mode edge. These current perturbations directly influence the edge q profile, decoupling the resonant location and instability drive of pedestal-localized microtearing modes (MTMs.) By exploiting this effect, we develop and apply a new experimental technique to track the dynamical frequency evolution of MTMs in the pedestal region, providing a compelling validation of the MTM model. The frequency of potential MTMs is calculated as the Doppler-shifted electron diamagnetic frequency at rational q=m/n surfaces, showing remarkable agreement with chirped frequency behavior of n=3, 4 and 5 modes detected with fast magnetics. Data is collected throughout multiple ELM cycles in order to build robust statistics describing the time-dependent frequency evolution of MTMs, which can be explained by examining the recovery of pedestal gradients after an ELM event. MTMs have a dominant transport contribution in the electron thermal channel, so the presented results indicate that reduced models of pedestal transport must be electromagnetic in nature and constructed with accurate calculations of MTM stability; inclusion of this physics is essential for accurate predictions of the electron temperature pedestal profile. Supporting measurements of mode saturation, propagation direction and transport fingerprints are made to support the dynamic frequency determination, unambiguously and experimentally identifying MTMs in the pedestal region of DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dislocation nano-hydrides in nickel: Nucleation, evolution and effects on dislocation behaviors

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width d eq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: d eq decreases for θ > 30°, d eq remains unchanged for θ = 30°, and d eq increases for θ < 30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. In conclusion, these findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Dislocations↗
DOE OSTI2024

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗
DOE OSTI2020

Radiation Damage in Materials—Helium Effects

Despite its scarcity in terrestrial life, helium effects on microstructure evolution and thermo-mechanical properties can have a significant impact on the operation and lifetime of applications, including: advanced structural steels in fast fission reactors, plasma facing and structural materials in fusion devices, spallation neutron target designs, energetic alpha emissions in actinides, helium precipitation in tritium-containing materials, and nuclear waste materials. The small size of a helium atom combined with its near insolubility in almost every solid makes the helium–solid interaction extremely complex over multiple length and time scales. This Special Issue, “Radiation Damage in Materials—Helium Effects”, contains review articles and full-length papers on new irradiation material research activities and novel material ideas using experimental and/or modeling approaches. These studies elucidate the interactions of helium with various extreme environments and tailored nanostructures, as well as their impact on microstructural evolution and material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling and Analysis of Wind Turbine Wake Vortex Evolution Due to Time-Constant Spatial Variations in Atmospheric Flow

Modern utility-scale wind turbines are evolving toward larger, lighter, and more flexible designs to meet the growing demand for renewable energy while minimizing logistical costs. However, these advancements in lightweight design result in heightened aeroelastic sensitivity, leading to complex interactions which affect the rotor’s capacity to withstand aerodynamic loading and the cascading effects that manifest in the wake’s vortex-structure evolution under variable atmospheric conditions. In this paper, we analyze the influence of stream-wise fluctuating atmospheric flow conditions on wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex–wake interactions for a range of conditions with spatially variant atmospheric flow. CODEF solutions for turbine operation in wind conditions featuring only one parameter fluctuation are compared to wind conditions with several wind parameter variations in combination. By isolating individual inflow variations and comparing them to multi-parameter scenarios, we determine the contributions of each atmospheric factor to rotor dynamics, wake evolution, and downstream wind farm interactions. The purpose of this paper is to analyze the effects of spatial variations in atmospheric flow on the topological evolution of wind turbine vortex wakes, which constitutes a gap in the current understanding of wind turbine wake dynamics. The insights gained from this study are particularly valuable for the development of wind farm control strategies aimed at mitigating the adverse effects of wake interactions, enhancing energy capture, and improving the overall stability of wind farm operations. With these insights, we aim to contribute to the development of modeling and simulation tools to optimize utility-scale wind power plants operating in diverse atmospheric environments.

Farrell, Alayna (ORCID:000000023555720X)↗
DOE OSTI2022

Thermal dynamics of P2-Na 0.67 Ni 0.33 Mn 0.67 O 2 cathode materials for sodium ion batteries studied by in situ analysis

Cobalt-free layered Na 0.67 Ni 0.33 Mn 0.67 O 2 is an attractive cathode material for sodium ion batteries. The thermal stability of cathode materials is crucial to their practical applications. In the current work, we investigate structural and morphological evolution in P2-type Na 0.67 Ni 0.33 Mn 0.67 O 2 cathode materials during annealing via in situ synchrotron X-ray diffraction and transmission electron microscopy. Insights are obtained from the two complementary in situ characterizations (at different length scales) in term of the thermal stability of P2-Na 0.67 Ni 0.33 Mn 0.67 O 2 cathode materials. Here, the results indicate that the hexagonal P2 phase remains unchanged during the heat treatment process, and thermally driven expansion/contraction of lattice parameters exhibit an anisotropic character in the a and c directions. In addition, interfaces/grain boundaries play an important role in the structural stability, which leads to the distinct morphological evolution between the polycrystalline and single-crystal particles.

25 ENERGY STORAGE↗
DOE OSTI2023

Scenario optimization for the tokamak ramp-down phase in RAPTOR: Part B. safe termination of DEMO plasmas

An optimized plasma current ramp-down strategy is critical for safe and fast termination of plasma discharges in a tokamak demonstration fusion reactor (DEMO), both in planned and emergency scenarios, avoiding plasma disruptions and excessive heat loads to the first wall. Plasma stability limits and machine-specific technical requirements constrain the stable envelope through which the plasma must be navigated. Large amounts of auxiliary heating are required throughout the ramp-down phase, to avoid a radiative collapse in the presence of intrinsic tungsten and seeded xenon impurities, as quantitatively estimated in this work. As the plasma current is reduced, the current density becomes increasingly peaked, reflected by a growing value of the internal inductance $\ell$ $i3$ , resulting in reduced controllability of the vertical position of the plasma. The feasibility of different plasma current ramp-down rates is tested by applying an automated optimization framework embedding the RAPTOR core transport solver. Optimal time traces for plasma current $I$ $p$ $(t)$ and plasma elongation $κ(t)$ are proposed, to satisfy an $I$ $p$ -dependent upper limit on the plasma internal inductance, as obtained from vertical stability studies using the CREATE-NL code, as well as a constraint on the time evolution of $q$ 95 , to avoid an ideal MHD mode. A negative current density near the plasma edge is observed in our simulations, even for the most conservative $I$ $p$ ramp-down rate, indicating significant transient dynamics due to a large resistive time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Stability of Neptune’s Distant Resonances in the Presence of Planet Nine

Trans-Neptunian Objects (TNOs) in the scattered disk with 50 ≲ a ≲ 100 au are thought to cluster near Neptune’s n:1 resonances (e.g., 3:1, 4:1, and so on). While these objects spend lengthy periods of time at large heliocentric distances, if their perihelia remain less than around 40 au, their dynamical evolution is still largely coupled to Neptune’s. Conversely, around a dozen extreme TNOs with a ≳ 250 au and detached perihelia seem to exist in a regime where they are too distant to be affected by the giant planets and too close for their dynamics to be governed by external forces. Recent work suggests that the apparent alignment of these orbits in physical space is a signature of gravitational shepherding by a distant massive planet. In this paper, we investigate the evolution of TNOs in each of Neptune’s n:1 resonances between the 3:1 and 14:1. We conclude that both resonant and nonresonant objects beyond the 12:1 near ∼157 au are removed rather efficiently via perturbations from the hypothetical Planet Nine. Additionally, we uncover a population of simulated TNOs with a ≲ 100 au, 40 ≲ q ≲ 45 au, and low inclinations that experience episodes of resonant interactions with both Neptune and Planet Nine. Finally, we simulate the evolution of observed objects with a > 100 au and identify several TNOs that are potentially locked in n:1 resonances with Neptune, including the most distant known resonant candidates, 2014 JW{sub 80} and 2014 OS{sub 394},which appear to be in the 10:1 and 11:1 resonances, respectively. Our results suggest that the detection of similar remote objects might provide a useful constraint on hypotheses invoking the existence of additional distant planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Transport barrier and spinning blob dynamics in the tokamak edge

Abstract In this work, we investigate the dynamics of plasma blobs in the edge of magnetic confinement devices using a full-f gyrokinetic particle-in-cell code with X-point geometry. In simulations, the evolution of a seeded blob is followed as it approaches a naturally-forming zonal shear layer near the separatrix, where the blob is stabilized by a large spin induced by the self-consistent adiabatic electron response, and blob bifurcation and trapping are observed during the cross-field propagation of blobs. A new theoretical explanation in both the zonal free and zonal shear layer is constructed, where the dominant E × B spin motion is included. A theoretical condition for a transport barrier induced by the interaction between spinning blobs and the zonal shear layer is obtained, and its scaling is verified with simulations. The new theoretical framework, especially the transport barrier, is applicable to explain and predict various experimental phenomena. In particular, the transport barrier condition calculated with experimental parameters demonstrates that the blob radial transport for H mode is smaller than L mode in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Demonstration of a quantum-classical coprocessing protocol for simulating nuclear reactions

Quantum computers hold great promise for exact simulations of nuclear dynamical processes (e.g., scattering and reactions), which are paramount to the study of nuclear matter at the limit of stability and in the formation of chemical elements in stars. However, quantum simulations of the unitary (real) time dynamics of fermionic many-body systems require a currently prohibitive number of reliable and long-lived qubits. Here we propose a co-processing algorithm for the simulation of real-time dynamics in which the time evolution of the spatial coordinates is carried out on a classical processor, while the evolution of the spin degrees of freedom is carried out on quantum hardware. We demonstrate this hybrid scheme with the simulation of two neutrons scattering at the Lawrence Berkeley National Laboratory's Advanced Quantum Testbed. After implementing error mitigation strategies to improve the accuracy of the algorithm in addition to a combination of circuit compression techniques and tomography as methods to elucidate the onset of decoherence, our results validate the principle of the proposed co-processing scheme. A generalization of this present scheme will open the way for (real-time) path integral simulations of nuclear scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗
DOE OSTI2025

Role of Mg:Al ratio and Pt promotion on mitigating the deactivation of Ni-based methane steam reforming catalysts

Ni-based catalysts for CH 4 steam reforming require fine-tuning to withstand deactivation under dynamic operation conditions relevant to emerging H 2 -driven technologies. This study investigates the impact of the Mg:Al ratio and Pt presence in catalyst composition on the activity and stability of industrially relevant Ni-based mono- and bimetallic catalysts during simulated daily start-up and shut-down cycles in various gas atmospheres. The evolution of the catalyst structure during extensive testing procedures was investigated in detail by complementary electron microscopy, in situ/operando XAS and XRD. The results obtained revealed that catalyst deactivation is promoted at high Mg:Al ratios and especially affects the monometallic catalysts. By converting CH 4 at lower temperatures, Pt regulates the extent of Ni oxidation and its incorporation into MgO lattice. Further prevention of catalyst deactivation was achieved by optimizing the reactor shut-down procedure to minimize the simultaneous exposure to high temperatures and H 2 O vapors. By flushing the reactor with N 2 only around the reaction extinction temperature, a fraction of Ni species is maintained in metallic state, which is beneficial for the long-term activity and reaction operation economy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Phase Transitions of Eutectic High Entropy Alloy AlCoCrFeNi 2.1 Under Shock Compression

High entropy alloys (HEAs) are a new class of metals that exhibit unique mechanical performance. Among HEAs, additively manufactured eutectic high entropy alloys (AM EHEAs) have recently emerged as candidate materials for use in extreme conditions due to their simultaneous high strength and ductility. However, the deformation and structural evolution of AM EHEAs under conditions of high pressure have not been well characterized, limiting their use in extreme applications. Dynamic compression experiments and molecular dynamics simulations are presented to study the structural evolution of AM EHEA AlCoCrFeNi 2.1 when compressed to pressures up to 400 GPa. In situ X-ray diffraction measurements capture the appearance of face-centered cubic and body-centered cubic phases at different pressure conditions, with pure- and mixed-phase regions. Furthermore, understanding the phase stability and structural evolution of the AM EHEA offers new insights to guide the development of high-performance complex materials for extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Non-autoregressive time-series methods for stable parametric reduced-order models

Advection-dominated dynamical systems, characterized by partial differential equations, are found in applications ranging from weather forecasting to engineering design where accuracy and robustness are crucial. There has been significant interest in the use of techniques borrowed from machine learning to reduce the computational expense and/or improve the accuracy of predictions for these systems. These rely on the identification of a basis that reduces the dimensionality of the problem and the subsequent use of time series and sequential learning methods to forecast the evolution of the reduced state. Often, however, machine-learned predictions after reduced-basis projection are plagued by issues of stability stemming from incomplete capture of multiscale processes as well as due to error growth for long forecast durations. To address these issues, we have developed a non-autoregressive time series approach for predicting linear reduced-basis time histories of forward models. In particular, we demonstrate that non-autoregressive counterparts of sequential learning methods such as long short-term memory (LSTM) considerably improve the stability of machine-learned reduced-order models. Further, we evaluate our approach on the inviscid shallow water equations and show that a non-autoregressive variant of the standard LSTM approach that is bidirectional in the principal component directions obtains the best accuracy for recreating the nonlinear dynamics of partial observations. Moreover-and critical for many applications of these surrogates-inference times are reduced by three orders of magnitude using our approach, compared with both the equation-based Galerkin projection method and the standard LSTM approach.

97 MATHEMATICS AND COMPUTING↗