Search NASA⌕ Search

SEARCH · Search NASA

Results for “mechanical instabilities”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stable low diffusion flux splitting schemes on unstructured meshes

Shock instabilities are shown to manifest in modern low-diffusion flux-vector splitting (FVS) schemes when used on unstructured meshes, or situations where shocks do not align with the mesh lines. These instabilities occur irrespective of the Mach number of the shock. Three types of dissipative mechanisms that suppress these instabilities are presented. These mechanisms are carefully designed in order to affect only problematic regions of the flux-splittings. The AUSM + and LDFSS schemes are stabilized using the proposed modifications. It is shown that the added dissipation improves the shock behavior of AUSM and LDFSS on unstructured meshes. It is also shown that the AUSM + -up scheme is prone to the “carbuncle” instability, a specific type of shock instability, when used on unstructured meshes. The modifications proposed in this work do not lead to carbuncle instabilities for the problems considered here. Furthermore, the modified schemes are shown to satisfy certain properties that are crucial for accurate shear layer computations, such as stationary contact preservation. Using benchmark problems, it is demonstrated that despite the diffusion added for stabilization, these schemes are not overly diffusive. Furthermore, due to these advantages, the modified FVS schemes presented here are promising candidates for high-speed compressible flow computations on unstructured meshes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Effect of plasma initialization on 3D PIC simulation of Hall thruster azimuthal instability

Abstract The lack of understanding of the azimuthal instability and the resulting electron anomalous transport limits further improvement of Hall thrusters. Compared to theoretical and experimental approaches, the numerical particle-in-cell (PIC) simulation is a suitable and powerful tool, which has been widely applied to investigate the azimuthal instability, and great progress has been made in the past decades. However, PIC simulations are intrinsically computationally expensive, and it is realized that the Hall thruster azimuthal instability has a three dimensional nature. Therefore, massive 3D PIC simulation must be carried out to completely reveal the mechanism of the instability. In this paper, the effect of plasma initialization on 3D PIC simulation of Hall thruster azimuthal instability is studied as a starting point. It is found that by initializing with ion density and velocity fitting functions to the steady-state simulation results, a faster convergence can be obtained and the computational time can be reduced by about 1.5 times. Typical fitting functions of ion density, drifting velocity, and temperature are given, and the influence of different initialization profiles is presented.

Physics↗
DOE OSTI2021

Quantification of Variation in Combustion Instability Amplitude in a Multi-Nozzle Can Combustor

In this work, we quantify the level of variation in unstable combustion oscillation amplitudes and identify the source of this variation in a multi-nozzle can combustor. At conditions where pollutant emissions are reduced, lean-premixed combustors can undergo thermoacoustic instability when pressure and heat release rate oscillations couple. A commonly used method for suppressing instability is fuel staging, which is a method where fuel is unevenly distributed between nozzles in a multi-nozzle combustor. Our work follows others who have characterized the effect of fuel staging on combustion instability and the mechanisms by which it works during both steady-state and transient operation. One of the outcomes from our previous work was that certain instability operating points display a high level of pressure oscillation amplitude variation. Instead of oscillating at a constant limit-cycle amplitude, the pressure oscillation amplitude varies significantly in time. In this work, we use the concept of permutation entropy to quantify the level of variation in the pressure fluctuation amplitude. We correlate the level of variation with a number of state variables over a range of operating conditions, 291 test cases in all. These state variables include mixture equivalence ratio; transient timescale, amplitude, and direction; hardware temperatures; gas temperatures; and thermoacoustic damping and growth rates. Significant pressure oscillation amplitude variation occurs when the thermoacoustic damping rate is low. The damping and growth rates can be low for a number of reasons, but they are highly correlated with the metal temperature of the centerbody, where the flames are anchored; lower temperatures result in lower damping rates during stable operation and lower growth rates during unstable operation. These results show the importance of the thermal boundary condition on the time-dependent behavior of the thermoacoustic instability.

Westfall, Seth↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2024

Investigation of Tail-dominated Instability in the Fermilab Recycler Ring

In our recent operational run, a single bunch, tail-dominated instability was observed in the Fermilab Recycler ring. This instability exclusively occurs in the vertical plane when the chromaticity is close to zero. In this study, we conduct a detailed analysis of this instability under different operational parameters. We investigate the impact of space charge on the head-tail motion and propose potential interpretations of the underlying mechanism of the instability. Moreover, we explore methods to mitigate this instability in the future.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An optimization-based approach to tailor the mechanical response of soft metamaterials undergoing rate-dependent instabilities

An optimization-based design framework is proposed to tune the response of soft metamaterials involving both geometric instabilities and nonlinear viscoelastic material behavior. Designing the response of soft metamaterials to harness instabilities and undergo large, tailored configuration changes will enable advancements in soft robotics, shock and vibration mitigation, and flexible electronics. In line with the metamaterial concept, the response of these materials is governed to a large extent by the geometric and topological makeup of their small-scale features. However, the link between structure and response is less intuitive for soft metamaterials due to their reliance upon highly nonlinear responses triggered by geometric instabilities. This is further complicated by the effects of viscoelastic relaxation, which recent studies have shown to alter the emergence of instabilities in non-intuitive ways. Here, these effects are accounted for in our framework to achieve various design objectives, including tailored force–displacement response and maximized energy absorption from both geometric and material effects. To fully automate this process, it is essential to have a completely robust equation solver for forward problems involving instabilities and viscoelastic relaxation. We achieve this by casting the search for stable mechanical equilibrium — i.e. the forward problem — as a minimization problem and utilize a trust region algorithm to robustly handle instabilities and follow energetically-favorable equilibrium paths through critical points.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Further developments on observations of the Taylor instability in a dusty plasma

We previously reported the observation of a hydrodynamic-type interchange instability occurring sporadically at the curved interface between two dust clouds of different densities, in a three-dimensional dusty plasma [Pacha et al., Phys. Plasmas 19, 014501 (2012)]. The mechanism for the instability was not conclusively identified. Further analysis of the video images has been performed, showing that the instability is associated with the appearance of large amplitude dust acoustic waves or shocks. Evidence of a dust acoustic shock at the boundary, which could drive a Richtmyer–Meshkov-type instability, will be presented. We also discuss the possibility, suggested by Avinash and Sen [Phys. Plasmas 22, 083707 (2015)], that the observed instability is due to a Rayleigh–Taylor instability driven by a component of the gravitational acceleration, which is perpendicular to the curved boundary.

Pacha, K. A.↗
DOE OSTI2022

Generation and Control of Self-Organized Nonlinear Kinetic Structures in High Energy Density Plasmas in the Presence of Intense Magnetic Fields and Ultrashort Laser Pulses

In the project we investigated mechanisms to sculpt phase space by controlling plasma instabilities. We found that linear theory, which applies in the initial growth phase of an instability, allows for initial conditions that can suppress instabilities. This provides a possible mechanism to release the energy available to an instability on demand. We developed various computational techniques that can reduce computational cost of simulations while maintaining or improving physical fidelity. We considered the Buneman instability in the presence of a ponderomotively driven plasma wave to examine the interplay between an instability and steady-state plasma response.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multimodal In Situ Characterization Uncovers Unexpected Stability of a Cobalt Electrocatalyst for Acidic Sustainable Energy Technologies

An accelerated development of durable and affordable sustainable energy technologies is often hindered by a limited understanding of how non-precious materials within these systems degrade. In acidic proton exchange membrane fuel cells and water electrolyzers, metallic cobalt (Co) is considered an unstable component that is often combined with precious metals or other stabilizers. To understand the mechanisms behind Co instability, we employ an experimental platform that quantifies dissolution with on-line inductively coupled plasma mass spectrometry and product formation with electrochemical mass spectrometry during electrochemical testing, along with ex- situ characterization. Under varied conditions (electrocatalysis, time, gas-type saturation, and ion concentration), windows of Co stability are observed that are different than predicted with classical chemical thermodynamics, suggesting new stabilization and degradation mechanisms than previously understood. Notably, Co is active for the hydrogen evolution reaction (HER), with prolonged stability that is ~300 mV different than thermodynamically projected. Additionally, in an oxygenated environment, Co concurrently performs the HER and oxygen reduction reaction (ORR) yet undergoes different morphology changes and dissolution mechanisms. Interestingly, in the absence of electrocatalysis, there is a 22x decrease in dissolution in an oxygen-free environment, proposing a route to decrease Co losses during device shutdown protocols. Lastly, under more extreme operating conditions, Co becomes stable after a substantial amount of dissolution, suggesting that high concentrations of Co 2+ ions in the microenvironment induce the formation of a stable CoHO 2 surface. Altogether, these results can be leveraged to improve the design and development of more robust and cost-effective sustainable energy technologies, as well as promote strategic strategies for prolonged material utilization.

Dissolution↗
DOE OSTI2022

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Adaptive Control of Distributed Energy Resources for Distribution Grid Voltage Stability

Volt-VAR and Volt-Watt functionality in photovoltaic (PV) smart inverters provide mechanisms to ensure system voltage magnitudes and power factors remain within acceptable limits. However, these control functions can become unstable, introducing oscillations in system voltages when not appropriately configured or maliciously altered during a cyberattack. In the event that Volt-VAR and Volt-Watt control functions in a portion of PV smart inverters in a distribution grid are unstable, the proposed adaptation scheme utilizes the remaining and stably-behaving PV smart inverters and other Distributed Energy Resources to mitigate the effect of the instability. The adaptation mechanism is entirely decentralized, model-free, communication-free, and requires virtually no external configuration. Here we provide a derivation of the adaptive control approach and validate the algorithm in experiments on the IEEE 37 and 8500 node test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Decoding the formation of hammerhead ion populations observed by Parker Solar Probe

Context. In situ observations by the Parker Solar Probe (PSP) have revealed new properties of the proton velocity distributions (VDs), including hammerhead features that suggest a non-isotropic broadening of the beams. Aims. The present work proposes a very plausible explanation for the formation of hammerhead proton populations through the action of a proton firehose-like instability triggered by the proton beam. Methods. We investigated a self-generated firehose-like instability driven by the relative drift of ion populations using a simplified moment-based quasi-linear (QL) theory. While simpler and faster than advanced numerical simulations, this toy model provided rapid insights and concisely highlighted the role of plasma micro-instabilities in relaxing the observed anisotropies of particle VDs in the solar wind and space plasmas. Results. The QL theory proposed here shows that the resulting transverse waves are right-hand polarized and have two consequences on the protons: (i) They reduce the relative drift between the beam and the core, but above all, (ii) they induce a strong perpendicular temperature anisotropy specific to the observed hammerhead ion beam. Moreover, the long-run QL results suggest that these hammerhead distributions are rather transitory states that are still subject to relaxation mechanisms, in which instabilities such as the one discussed here are very likely involved. This foundational work motivates future detailed studies using advanced methods.

Shaaban, Shaaban M. (ORCID:000000030465598X)↗
DOE OSTI2024

The radiation instability of thermally stable nanocrystalline platinum gold

Here, recent experimentally validated alloy design theories have demonstrated nanocrystalline binary alloys that are stable against thermally induced grain growth. An open question is whether such thermal stability also translates to stability under irradiation. In this study, we investigate the response to heavy ion irradiation of a nanocrystalline platinum gold alloy that is known to be thermally stable from previous studies. Heavy ion irradiation was conducted at both room temperature and elevated temperatures on films of nanocrystalline platinum and platinum gold. Using scanning/transmission electron microscopy equipped with energy-dispersive spectroscopy and automated crystallographic orientation mapping, we observe substantial grain growth in the irradiated area compared to the controlled area beyond the range of heavy ions, as well as compositional redistribution under these conditions, and discuss mechanisms underpinning this instability. These findings highlight that grain boundary stability against one external stimulus, such as heat, does not always translate into grain boundary stability under other stimuli, such as displacement damage.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interfaces in all solid state Li-metal batteries: A review on instabilities, stabilization strategies, and scalability

With technological advancements in electrochemical energy storage systems increasing at a spectacular rate, batteries equipped with a lithium anode hold the key towards unlocking high energy densities. While lithium-ion batteries with layered anodes (e.g. graphite) and liquid organic electrolytes have been ubiquitous in portable electronics, electric vehicles, and grid applications, all solid-state batteries that use the combination of a lithium anode and a solid-state electrolyte (SSE) will further advance the present technology. The underlying challenge that limits the successful development of all solid-state batteries (ASSBs) is dictated largely by the highly reactive interfaces at the anode/SSE and the cathode/SSE interface. In this comprehensive review, we present an overview of the following: (i) characterization of the electrode/SSE interface via multimodal characterization, which include X-ray-, electron-, neutron-, optical-, and computation-based methods: (ii) parasitic reactions that occur from chemical, mechanical or electrochemical instabilities and interfacial engineering strategies for improving the stability of these interfaces classified by the class of inorganic SSEs (sulfide-, NASICON-, and garnet-type SSEs); (iii) laboratory-to-industry scale processing perspectives of SSEs; (iv) scalability and manufacturing aspects of current interfacial strategies; and (v) the prospects of all ASSBs within the context of extreme fast charging capability. Here, this review seeks to highlight key efforts in the field of ASSBs, by focusing particularly on stabilizing the electrode/SSE interfaces, which will help to bridge fundamental studies to technological relevance.

25 ENERGY STORAGE↗
DOE OSTI2022

Low-Reynolds-number oscillating boundary layers on adiabatic slopes

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($Re\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

42 ENGINEERING↗
DOE OSTI2021

Identifying the microtearing modes in the pedestal of DIII-D H-modes using gyrokinetic simulations

Recent evidence points toward the microtearing mode (MTM) as an important fluctuation in the H-mode pedestal for anomalous electron heat transport. A study of the instabilities in the pedestal region carried out using gyrokinetic simulations to model an ELMy H-mode DIII-D discharge (USN configuration, 1.4 MA plasma current, and 3 MW heating power) is presented. The simulations produce MTMs, identified by predominantly electromagnetic heat flux, small particle flux, and a substantial degree of tearing parity. The magnetic spectrogram from Mirnov coils exhibits three distinct frequency bands—two narrow bands at lower frequency (~35–55 kHz and ~70–105 kHz) and a broader band at higher frequency (~300–500 kHz). Global linear GENE simulations produce MTMs that are centered at the peak of the ω* profile and correspond closely with the bands in the spectrogram. The three distinctive frequency bands can be understood from the basic physical mechanisms underlying the instabilities. For example (i) instability of certain toroidal mode numbers (n) is controlled by the alignment of their rational surfaces with the peak in the ω* profile, and (ii) MTM instabilities in the lower n bands are the conventional collisional slab MTM, whereas the higher n band depends on curvature drive. While many features of the modes can be captured with the local approximation, a global treatment is necessary to quantitatively reproduce the detailed band gaps of the low-n fluctuations. Notably, the transport signatures of the MTM are consistent with careful edge modeling by SOLPS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗