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

Isotope physics of heat and particle transport with tritium in JET-ILW type-I ELMy H-mode plasmas

As part the DTE2 campaign in the JET tokamak, we conducted a parameter scan in T and D-T complementing existing pulses in H and D. For the different main ion masses, type-I ELMy H-modes at fixed plasma current and magnetic field can have the pedestal pressure varying by a factor of 4 and the total pressure changing from <?CDATA $\beta_{\mathrm{N}} = 1.0$?> β N = 1.0 to 3.0. We investigated the pedestal and core isotope mass dependencies using this extensive data set. The pedestal shows a strong mass dependence on the density, which influences the core due to the strong coupling between both plasma regions. To better understand the causes for the observed isotope mass dependence in the pedestal, we analysed the interplay between heat and particle transport and the edge localised mode (ELM) stability. For this purpose, we developed a dynamic ELM cycle model with basic transport assumptions and a realistic neutral penetration. The temporal evolution and resulting ELM frequency introduce an additional experimental constraint that conventional quasi-stationary transport analysis cannot provide. Our model shows that a mass dependence in the ELM stability or in the transport alone cannot explain the observations. One requires a mass dependence in the ELM stability as well as one in the particle sources. The core confinement time increases with pedestal pressure for all isotope masses due to profile stiffness and electromagnetic turbulence stabilisation. Interestingly, T and D-T plasmas show an improved core confinement time compared to H and D plasmas even for matched pedestal pressures. For T, this improvement is largely due to the unique pedestal composition of higher densities and lower temperatures than H and D. With a reduced gyroBohm factor at lower temperatures, more turbulent drive in the form of steeper gradients is required to transport the same amount of heat. This picture is supported by quasilinear flux-driven modelling using TGLF -SAT2 within Astra . With the experimental boundary condition TGLF -SAT2 predicts the core profiles well for gyroBohm heat fluxes <?CDATA ${\gt}15$?> > 15 , however, overestimates the heat and particle transport closer to the turbulent threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On Practical Aspects of Variational Consistency in Contact Dynamics

Usage of contact mechanics methodologies is a pervasive modeling requirement in dynamic simulations. While for some trivial problems, solutions taken from analytical geometry are available, use of a finite element framework is common to achieve formulation generality. This work explores two dynamic contact formulations: one based on the traditional node-to-segment (NTS) approach, and a variationally consistent segment-to-segment (STS) mortar formulation. The NTS formulation employed here enforces the constraints kinematically (i.e., the interpenetration is enforced to the solver tolerance), whereas the mortar approach uses Lagrange multipliers to enforce the contact constraints. Both approaches are implemented in the open-source finite element framework Multiphysics Object-Oriented Simulation Environment (MOOSE). The results highlight two relevant contact-interface-related dynamic phenomena in finite element simulations. First, stabilization of contact constraints is discussed, taking into account the evolution of the total energy in a benchmark problem. Second, the influence of finite element discretization on both of the aforementioned contact formulations is analyzed by exercising a large-deformation example with continuous relative sliding. Variationally consistent contact approaches such as the mortar formulation lead to improved energy preservation and avoid spurious excitation of the system's frequencies. This is especially relevant in settings where inertia and vibrations are of importance.

42 ENGINEERING↗
DOE OSTI2026

Pathways to High-Performance Salt Hydrate Thermochemical Energy Storage Materials and Systems

Thermochemical materials (TCMs) based on salt hydrates are promising for thermal energy storage as they combine high energy densities with low reaction temperatures. However, their adoption is hindered by poor structural integrity and degradation under hygrothermal cycling. Storage performance is governed not only by the chemical reaction, but also by the coupled thermo-chemo-mechanical behavior that evolves with cycling. Understanding and controlling this coupling across length scales (material-to-reactor) is necessary to improve TCM stability and lifetime. In this perspective, we discuss the shortcomings of current characterization approaches and emphasize the need for measuring transport properties and structural transformations using in situ techniques that capture the dynamic evolution of these materials. We also outline opportunities for multiscale modeling frameworks that link thermodynamics and mechanics, enabling predictive evaluation of composite architectures designed for cycling stability. We conclude by identifying research questions that must be addressed to transform TCMs into viable energy storage technologies.

Composites↗
DOE OSTI2021

Exomoons in Systems with a Strong Perturber: Applications to α Cen AB

The presence of a stellar companion can place constraints on occurrence and orbital evolution of satellites orbiting exoplanets, i.e., exomoons. In this work we revise earlier orbital stability limits for retrograde orbits in the case of a three-body system consisting of a star, planet, and satellite. The revised limit reads a{sub sat}{sup crit}≈0.668(1-1.236e{sub p}) for e {sub p} ≤ 0.8 in units of the Hill Radius and represents the lower critical orbit as a function of the planetary eccentricity e {sub p}. A similar formula is determined for exomoons hosted by planets in binary star systems, where e {sub p} is replaced with the components of free and forced eccentricity from secular orbit evolution theory. By exploring the dynamics of putative exomoons in α Centauri AB we find that the outer stability limit can be much less than half the Hill Radius due to oscillations in the planetary orbital eccentricity caused by the gravitational interaction with the binary star. We show, furthermore, how the resulting truncation of the outer stability limit can affect the outward tidal migration and potential observability of exomoons through transit-timing variations (TTVs). Typical TTV (rms) amplitudes induced by exomoons in binary systems are ≲10 minutes and appear more likely for planets orbiting the less massive stellar component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Electrolyte Role in SEI Evolution at Si in the Pre-lithiation Stage vs the Post-lithiation Stage

The formation and evolution of the dynamic solid electrolyte interphase (SEI) at the Si anode/electrolyte interface are yet to be completely understood to solve irreversible capacity loss and increase battery cycle life. Herein, the evolution of SEI and its dynamic properties at the Si anode/electrolyte interface are investigated in two electrolyte systems, a 1.2 M LiPF6 in EC: EMC 3:7 (wt%) electrolyte (referred to as Gen2) and a 1.2 M LiTFSI in EC: EMC 3:7 (wt%) electrolyte (referred to as LiTFSI). Two lithiation stages are studied: the pre-lithiation (pre-Li) SEI stage and the post-lithiation (post-Li) stage. Findings reveal at the pre-Li, SEI formation starts at an early potential and contributes to the greater mass gain in the Si/Gen2, and it is dominated by the formation of a non-uniform F- and P-rich layer in Si/Gen2, in contrast to a homogeneous F- and C-containing layer at the Si/LiTFSI interphase. The initially formed SEI in LiTFSI further benefits the charge transfer kinetics. At the post-Li stage, a more substantial SEI evolution is observed at Si/LiTFSI. This paper offers a foundational understanding of the SEI dynamic evolution with electrolyte dependence. Findings from this report offer important insights into solving the complex SEI stability issues on Si.

25 ENERGY STORAGE↗
DOE OSTI2021

Origin of increased helium density inside bubbles in Ni (1–$x$) Fe $x$ alloys

Due to virtually no solubility, He atoms implanted or created inside materials tend to form bubbles, which are known to damage material properties through embrittlement. Higher He density in nano-sized bub- bles was observed both experimentally and computationally in Ni (1–$x$) Fe $x$ -alloy samples compared to Ni. The bubbles in the Ni (1–$x$) Fe $x$ -alloys were observed to be faceted, whereas in elemental Ni they were more spherical. Molecular dynamics simulations showed that stacking fault structures formed around bubbles at maximum He density. Higher Fe concentrations stabilize stacking fault structures, suppress evolution of dislocation network around bubbles and suppress complete dislocation emission, leading to higher He density.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2026

Chemical and Structural Insights into Solid Electrolyte Interphase Evolution for Sodium Metal Electrodes

The solid electrolyte interphase (SEI) critically governs the reversibility of sodium metal batteries, through dynamically mediating ion transport and interfacial reactions. However, its kinetic evolution under operating conditions, and how it influences interfacial stability, remains poorly understood. Here, in this study, we reveal that the SEI undergoes coupled chemical and mechanical changes during sodium plating and stripping, leading to spatial and temporal heterogeneity that drives interfacial degradation. Synchrotron operando grazing-incidence wide-angle X-ray scattering and soft X-ray absorption spectroscopy capture the sequential formation and dissolution of inorganic SEI phases (NaF, NaH, NaOH, Na 2 PO 3 F), accompanied by depth-dependent alterations in organic SEI components. Mesoscale modeling connects this evolving SEI heterogeneity to localized current density fluctuations and stress accumulation at the Na interface, identifying pathways to electrically isolated sodium formation. These findings show that SEI instability fundamentally limits reversibility in sodium metal batteries, and that controlling SEI chemistry–mechanics coupling is essential to achieving its durability.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Unraveling the Crystallization Kinetics of 2D Perovskites with Sandwich-Type Structure for High-Performance Photovoltaics

2D perovskite solar cells with high stability and high efficiency have attracted significant attention. A systematical static and dynamic structure investigation is carried out to show the details of 2D morphology evolution. A dual additive approach is used, where the synergy between an alkali metal cation and a polar solvent leads to high-quality 2D perovskite films with sandwich-type structures and vertical phase segregation. Such novel structure can induce high-quality 2D slab growth and reduce internal and surface defects, resulting in a high device efficiency of 16.48% with enhanced continuous illumination stability and improved moisture (55-60%) and thermal (85°C) tolerances. Transient absorption spectra reveal the carrier migration from low n to high n species with different kinetics. An [PbI 6 ] 4- octagon coalescence transformation mechanism coupled with metal and organic cations wrapped is proposed. By solvent vapor annealing, a recrystallization and reorientation of the 2D perovskite slabs occurs to form an ideal structure with improved device performance and stability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

SEI Formation and Lithium-Ion Electrodeposition Dynamics in Lithium Metal Batteries via First-Principles Kinetic Monte Carlo Modeling

The stabilization and enhanced performance of lithium metal batteries (LMBs) depend on the formation and evolution of the Solid Electrolyte Interphase (SEI) layer as a critical component for regulating the Li metal electrodeposition processes. This study employs a first-principles kinetic Monte Carlo (kMC) model to simulate the SEI formation and Li + electrodeposition processes on a lithium metal anode, integrating both the electrochemical electrolyte reduction reactions and the diffusion events giving place to the SEI aggregation processes during battery charge and discharge processes. The model replicates the competitive interactions between organic and inorganic SEI components, emphasizing the influence of the cycling regime. Results indicate that grain boundaries within the SEI facilitate faster lithium-ion transport compared to crystalline regions, crucial for improving the performance and stability of LMBs. The findings underscore the importance of dynamic SEI modeling for further development of next-generation high-energy-density batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Distinctive features of structural evolution and thermodynamic response in wide-bandgap semiconductors driven by intense electronic excitation

Radiation-tolerant material selection requires balancing lattice rigidity, defect dynamics, and electronic stability, as shown by covalent SiC outperforming ionic Ga 2 O 3 and GaN under extreme environments. Responding to intense electronic excitation, irradiation-driven phase segregation (β → δ/κ in Ga 2 O 3 ) and core–shell track (disordered structure in GaN), accompanied by elemental redistribution, contrastingly, exceptional radiation tolerance manifested by comparatively minimal lattice distortion (0.17 % strain variation) was demonstrated in SiC. These differential responses are primarily attributed to two fundamental mechanisms: (i) thermodynamic driving forces governing defect migration and phase separation, and (ii) the synergistic effects of robust covalent bonding composition coupled with efficient defect recombination processes. Here, the stronger electron–phonon (e-ph) coupling in Ga 2 O 3 (4.34 × 1018 W m −3 K −1 ) and GaN (3.55 × 10 18 W m −3 K −1 ) enhances lattice energy deposition, triggering thermal spikes (ΔT ≫ T m ) and structural transition behaviors, whereas weaker e-ph coupling in SiC (3.69 × 10 18 W m −3 K −1 ), relatively high thermodynamic parameters and efficient energy dissipation suppress thermal spikes to maintaining lattice integrity. The photoresponse degradation driven by enhanced radiative recombination is dominant in N-doped SiC, while V-doped systems achieve defect-mediated photoconduction optimization characterized by abrupt current transitions, matching fluorescence yield evolutions, and directly connecting defect engineering to optoelectronic performance.

Intense electronic excitation↗
DOE OSTI2025

Cu Evolution over Bimetallic Cu‐Y/Beta Zeolite Under H 2 and Ethanol Atmospheres: Unveiling the Role of Diatomic Metal–Metal Interactions

Understanding the dynamic evolution of Cu species under varying environmental conditions is critical for addressing challenges related to the activity and the stability of copper‐based catalysts in thermo‐, photo‐, and electrocatalysis. However, metal–metal interactions between dual single atoms and their effects on Cu evolution after exposure to different environmental molecules remain underexplored. Herein, we synthesized bimetallic Cu‐Y/Beta catalysts with dual single‐atom Cu and Y sites and monometallic Cu‐Beta catalysts with isolated Cu sites in dealuminated Beta zeolites. By varying Cu and Y compositions, diatomic interactions were studied under H 2 and ethanol atmospheres. With 6 wt% Y loading, approximately 0.4 wt% of Cu species in Cu‐Y/Beta remained partially oxidized as Cu(I) after reduction in pure H 2 at 350 °C, in contrast to the full transition to metallic Cu observed in Cu‐Beta. Combining X‐ray absorption spectroscopy with kinetic studies revealed that metallic Cu became the predominant species after reduction with H 2 as Cu loading increased from 0.4 to 1.7 wt%, quadrupling the initial ethanol dehydrogenation rate and demonstrating the dominant role of Cu(0) sites. In conclusion, scanning transmission electron microscopy and density functional theory simulations indicated spatial proximity between dual single‐atom Cu and Y sites and elucidated Cu speciation controlled by diatomic interactions.

Copper↗
DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
DOE OSTI2024

Electrochemical Control of the Ultrafast Lattice Response of a Layered Semimetal

The unique layer-stacking in two-dimensional (2D) van der Waals materials facilitates the formation of nearly degenerate phases of matter and opens novel routes for the design of low-power, reconfigurable functional materials. Electrochemical ion intercalation between stacked layers offers a promising approach to stabilize bulk metastable phases and to explore the effects of extreme carrier doping and strain. However, in situ characterization methods to study the structural evolution and dynamical functional properties of these intercalated materials remains limited. Here a novel experimental platform is presented capable of simultaneously performing electrochemical lithium-ion intercalation and multimodal ultrafast characterization of the lattice using both electron diffraction and nonlinear optical techniques. Using the layered semimetal WTe 2 as a model system, the interlayer shear phonon mode that modulates stacking between 2Dlayers is probed, showing that small amounts of lithiation enhance the amplitude and lifetime of the phonon, contrary to expectations. This results from the dynamically fluctuating and anharmonic structure between nearly degenerate phases at room temperature, which can be stabilized by electronic carriers accompanying the inserted lithium ions. At high lithiation, the T d ’ structure emerges and quenches the phonon response. This work defines new approaches for using electrochemistry to engineer the dynamic structure of 2D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unraveling Morphology and Chemistry Dynamics in Fluoroethylene Carbonate Generated Silicon Anode Solid Electrolyte Interphase Across Delithiated and Lithiated States: Relative Cycling Stability Enabled by an Elastomeric Polymer Matrix

The silicon solid electrolyte interphase (SEI) faces cyclical cracking and reconstruction due to the ∼350% volume expansion. Understanding the SEI dynamic morphology and chemistry evolution from delithiated to lithiated states is thereby paramount to engineering a stable Si anode. Fluoroethylene carbonate (FEC) is a preferred additive with widely demonstrated enhancement of the Si cycling. Thus, insights into the dynamics of the FEC-SEI may provide hints toward engineering the Si interface. Herein, complementary ATR-FTIR, AFM, tip IR, and XPS probing reveal the presence of an elastomeric polycarbonate-like matrix in the FEC-generated SEI which is absent from the FEC-free SEI. Adding FEC to the baseline 1 M LiPF 6 in EC:EMC (1:1) electrolyte promotes formation of a thinner and more conformal SEI, and subdues morphology and chemistry changes between consecutive half-cycles. From AFM, morphological stabilization of the FEC-SEI occurs earlier. Furthermore, conventional SEI biproducts such as Li 2 CO 3 and LiEDC appear in reduced quantities in the FEC-SEI implying a reduced quantity of Li-consuming species. The thin polymeric FEC-SEI enables deeper (de)lithiation of silicon. In conclusion, the enhanced mechanical compliance, chemical invariance, and reduced Li inventory consumption of the FEC-SEI are logically the key features underlying the Si cycling enhancement by FEC.

25 ENERGY STORAGE↗
DOE OSTI2024

Thermally induced structural evolution and nanoscale interfacial dynamics in Bi-Sb-Te layered nanostructures

Layered chalcogenides, including Bi-Sb-Te ternary alloys and heterostructures, are renowned as thermoelectric and topological insulators and have recently been highlighted as plasmonic building blocks beyond noble metals. Here, we conduct joint in situ transmission electron microscopy and density functional theory calculations to investigate the temperature-dependent nanoscale dynamics and interfacial properties, identifying the role of native defects and edge configurations in the anisotropic sublimation of Bi 2 Te 3 -Sb 2 Te 3 heterostructures and Sb 2-x Bi x Te 3 alloys. We report structural dynamics, including edge evolution, layer-by-layer sublimation, and the formation and coalescence of thermally induced polygonal nanopores. These nanopores are initiated by preferential dissociation of tellurium, reducing thermal stability in heterostructures. Triangular and quasi-hexagonal configurations dominate nanopore structures in heterostructures. Our calculations reveal antisite defects (Te Sb and Te Bi ) as key players in defect-assisted sublimation. These findings enhance our understanding of nanoscale dynamics and assist in designing tunable low-dimensional chalcogenides.

Bi2Te3-Sb2Te3 heterostructure↗