SEARCH · Search NASA
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.
Important issues: Mass gauging; Stability dynamics of disconnected capillary surfaces; Slow capillary driven flow (i.e. wicking); Long-term material property evolution in micro-g; Dumping problem with freezing of dump lines.
We study the long term orbital evolution of a terrestrial planet under the gravitational perturbations of a giant planet. In particular, we are interested in situations where the two planets are in the same plane and are relatively close. We examine both possible configurations: the giant planet orbit being either outside or inside the orbit of the smaller planet. The perturbing potential is expanded to high orders and an analytical solution of the terrestrial planetary orbit is derived. The analytical estimates are then compared against results from the numerical integration of the full equations of motion and we find that the analytical solution works reasonably well. An interesting finding is that the new analytical estimates improve greatly the predictions for the timescales of the orbital evolution of the terrestrial planet compared to an octupole order expansion. Finally, we briefly discuss possible applications of the analytical estimates in astrophysical problems.
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.
A low-density gas jet injected into a high-density ambient gas is known to exhibit self-excited global oscillations accompanied by large vortical structures interacting with the flow field. The primary objective of the proposed research is to study buoyancy effects on the origin and nature of the flow instability and structure in the near-field of low-density gas jets. Quantitative rainbow schlieren deflectometry, Computational fluid dynamics (CFD) and Linear stability analysis were the techniques employed to scale the buoyancy effects. The formation and evolution of vortices and scalar structure of the flow field are investigated in buoyant helium jets discharged from a vertical tube into quiescent air. Oscillations at identical frequency were observed throughout the flow field. The evolving flow structure is described by helium mole percentage contours during an oscillation cycle. Instantaneous, mean, and RMS concentration profiles are presented to describe interactions of the vortex with the jet flow. Oscillations in a narrow wake region near the jet exit are shown to spread through the jet core near the downstream location of the vortex formation. The effects of jet Richardson number on characteristics of vortex and flow field are investigated and discussed. The laminar, axisymmetric, unsteady jet flow of helium injected into air was simulated using CFD. Global oscillations were observed in the flow field. The computed oscillation frequency agreed qualitatively with the experimentally measured frequency. Contours of helium concentration, vorticity and velocity provided information about the evolution and propagation of vortices in the oscillating flow field. Buoyancy effects on the instability mode were evaluated by rainbow schlieren flow visualization and concentration measurements in the near-field of self-excited helium jets undergoing gravitational change in the microgravity environment of 2.2s drop tower at NASA John H. Glenn Research Center. The jet Reynolds number was varied from 200 to 1500 and jet Richardson number was varied from 0.72 to 0.002. Power spectra plots generated from Fast Fourier Transform (FFT) analysis of angular deflection data acquired at a temporal resolution of 1000Hz reveal substantial damping of the oscillation amplitude in microgravity at low Richardson numbers (~0.002). Quantitative concentration data in the form of spatial and temporal evolutions of the instability data in Earth gravity and microgravity reveal significant variations in the jet flow structure upon removal of buoyancy forces. Radial variation of the frequency spectra and time traces of helium concentration revealed the importance of gravitational effects in the jet shear layer region. Linear temporal and spatio-temporal stability analyses of a low-density round gas jet injected into a high-density ambient gas were performed by assuming hyper-tan mean velocity and density profiles. The flow was assumed to be non parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results were delineated. A decrease in the density ratio (ratio of the density of the jet to the density of the ambient gas) resulted in an increase in the temporal amplification rate of the disturbances. The temporal growth rate of the disturbances increased as the Froude number was reduced. The spatio-temporal analysis performed to determine the absolute instability characteristics of the jet yield positive absolute temporal growth rates at all Fr and different axial locations. As buoyancy was removed (Fr . 8), the previously existing absolute instability disappeared at all locations establhing buoyancy as the primary instability mechanism in self-excited low-density jets.
Magnetic reconnection is a fundamental process operating in laboratory, space, and astrophysical plasmas. We understand today that it is the mechanism behind the often explosive release of the stored magnetic energy, and that it directly or indirectly powers as diverse a set of phenomena as solar eruptions, stellar flares, geomagnetic storms, the aurora, and even disruption events in fusion machines. It is also the mechanism powering the majority of the deleterious space effects collectively referred to as “space weather.” Owing to its importance, reconnection has been a research focus for many years, culminating in the research based on NASA’s Magnetospheric Multiscale mission. A major component of foundational research involved studies of the stability and disruption of current sheets by collisionless tearing, which still plays a major role in the dynamical evolution of the interfaces between magnetic fields of different directions. This presentation will illuminate these beginnings, with special focus on the seminal work of Academician A. A. Galeev, and also provide a summary of the present state of knowledge. It will conclude with an outlook toward future reconnection research targets.
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.
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.
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.
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.
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.
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.
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.
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.
The Cassini mission has provided the best opportunity to date to extensively study the seasonal variation in Titan’s atmosphere, with observations spanning almost half a Titan year (𝐿𝑠 = 293 − 93°). An important feature in the Titan middle-atmosphere is the formation of a polar vortex. Observations have shown that an initially well-developed northern vortex enriched with trace gas species gradually breaks down after spring equinox as a new vortex emerges in southern winter. Here we use Cassini CIRS observations to derive the temperature and composition of the middle-atmosphere. We use the gradient wind equation to first estimate the mean zonal winds, and then the Potential Vorticity (PV) throughout Titan’s atmosphere over the timespan of the Cassini mission. PV is a useful diagnostic quantity for studying the dynamics of polar vortices because it is materially conserved for adiabatic and frictionless flows, and can be inverted to find all other dynamical fields. Our results show the formation of a strong zonal jet in the winter hemisphere, with wind velocities reaching 220 m/s, which is consistent with previous studies. An annular PV structure is also observed over the winter poles, whereby a ring of PV encircles a local minima over the pole. Such distributions are often found to be unstable without a restoring force, yet they are seen here in numerous observations in both the northern and southern hemispheres. A comparison with the annular Martian vortices shows that latent heat release from condensation or subsidence-induced adiabatic heating may explain the origin and stability of the annulus. Finally, we investigate the evolution of the size of the vortices and the role of strong PV gradients as a dynamical mixing barrier for trace gas species across the vortex edge. We find that longer lived gases are less confined to the vortex than those with shorter photochemical lifetimes.
Further results of research, reported in several previous NASA Tech Briefs articles, were obtained on a mathematical formalism for postinstability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence). To recapitulate: Fictitious control forces are introduced to couple the dynamical equations with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in ordinary perceived three-dimensional space is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. Con sequently, the postinstability motion is represented by a set of functions describing the evolution of such statistical quantities as expectations and higher moments, and this representation is stable. The previously reported findings are analyzed from the perspective of the authors Stabilization Principle, according to which (1) stability is recognized as an attribute of mathematical formalism rather than of underlying physics and (2) a dynamical system that appears unstable when modeled by differentiable functions only can be rendered stable by modifying the dynamical equations to incorporate intrinsic stochasticity.
This algorithm provides a new way to improve the accuracy and asymptotic behavior of a low-dimensional system based on the proper orthogonal decomposition (POD). Given a data set representing the evolution of a system of partial differential equations (PDEs), such as the Navier-Stokes equations for incompressible flow, one may obtain a low-dimensional model in the form of ordinary differential equations (ODEs) that should model the dynamics of the flow. Temporal sampling of the direct numerical simulation of the PDEs produces a spatial time series. The POD extracts the temporal and spatial eigenfunctions of this data set. Truncated to retain only the most energetic modes followed by Galerkin projection of these modes onto the PDEs obtains a dynamical system of ordinary differential equations for the time-dependent behavior of the flow. In practice, the steps leading to this system of ODEs entail numerically computing first-order derivatives of the mean data field and the eigenfunctions, and the computation of many inner products. This is far from a perfect process, and often results in the lack of long-term stability of the system and incorrect asymptotic behavior of the model. This algorithm describes a new stabilization method that utilizes the temporal eigenfunctions to derive correction terms for the coefficients of the dynamical system to significantly reduce these errors.
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.
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.