Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Stationary Crossflow Breakdown due to Mixed Mode Spectra of Secondary Instabilities

Numerical simulations are used to study laminar breakdown characteristics associated with stationary crossflow instability in the boundary-layer flow over a subsonic swept-wing configuration. Previous work involving the linear and nonlinear development of individual, fundamental modes of secondary instability waves is extended by considering the role of more complex, yet controlled, spectra of the secondary instability modes. Direct numerical simulations target a mixed mode transition scenario involving the simultaneous presence of Y and Z modes of secondary instability. For the initial amplitudes investigated in this paper, the Y modes are found to play an insignificant role during the onset of transition, in spite of achieving rather large, O(5%), amplitudes of RMS velocity fluctuation prior to transition. Analysis of the numerical simulations shows that this rather surprising finding can be attributed to the fact that the Y modes are concentrated near the top of the crossflow vortex and exert relatively small influence on the Z modes that reside closer to the surface and can lead to transition via nonlinear spreading that does not involve interactions with the Y mode. Finally, secondary instability calculations reveal that subharmonic modes of secondary instability have substantially lower growth rates than those of the fundamental modes, and hence, are less likely to play an important role during the breakdown process involving complex initial spectra.

Li, Fei↗
NASA NTRS2018

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗
NASA NTRS2021

Wake Instability Behind Isolated Trip Near the Leading Edge of the BOLT-II Configuration

The BOLT-II (Holden Mission) configuration is an extended version of the BOLT flight article and depicts hypersonic boundary-layer transition in the presence of multiple and potentially interacting instability mechanisms. Several numerical studies of the boundary-layer instabilities over these configurations have been reported in the recent literature, including our previous studies of the modal instability characteristics of boundary-layer streaks adjacent to the minor-axis symmetry plane of the BOLT configuration and the wake instabilities behind a diamond planform (“pizza-box”) trip along the symmetry plane on the secondary side of the BOLT-II configuration. The present work extends the latter study to a scaled version of the same trip that is located in the region of nonzero crossflow in the vicinity of the leading edge at X/L = 0.5. Collectively, the trips along the symmetry plane and near the leading edge are the focus of the NASA roughness experiment on the secondary side of the BOLT-II configuration. The laminar basic state computation at the nominal flight design condition of Re ∞ = 5.44 x 10 6 /m and Re ∞ = 2.5 x 10 6 /m shows that the leading-edge trip with k/δ ≈ 0.70 and planform-side-length-to-height ratio of b/k = 3.0 induces multiple asymmetric, longitudinal streaks within the trip wake. The most prominent streak among these resembles a finite amplitude crossflow vortex, and it supports the amplification of multiple families of unstable modes. The application of multi-dimensional instability analysis to the wake flow reveals that the most amplified unstable mode can achieve a peak N-factor of up to 20 by X/L = 0.80, indicating that the onset of transition is more than likely to occur within eighty percent of the model length. To the best of our knowledge, the present study represents the first analysis including nonparallel and curvature effects on hypersonic tripwake instabilities in the presence of boundary-layer crossflow over a three-dimensional configuration.

Boundary layer transition↗
DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Droplet breakup by multimodal nonlinear Rayleigh Taylor instability

A droplet impacted by a shock wave will undergo a process of fragmentation due to the development of interfacial hydrodynamic instabilities. The interface experiences variable acceleration and shear that result in the development of both inertial (Rayleigh–Taylor) and shear (Kelvin–Helmholtz) instabilities. These perturbations grow in time and drive the fragmentation and breakup of the deformed droplet. Experiments are performed on nominally 0.86 mm water droplet subjected to a Mach 7.6 detonation wave, resulting in a high Weber number (~36,000) breakup event. Perturbation growth is measured from a series of high-speed (> MHz) shadowgraph images. It is proposed that, given the size of the large-scale perturbations observed in experiments, these instabilities are growing in the non-linear regime and can be described by bubble-merger models for nonlinear mixing. Calculations are performed for the growth rates and size of these instabilities using deformation and external flow models to establish the time-dependent boundary conditions. The concurrence of the measured perturbation widths and the predictions of the simple model lend credence to the theory. This novel approach serves to open a new avenue in the characterization of droplet breakup via hydrodynamic instabilities.

Droplet breakup↗
DOE OSTI2025

The gyrokinetic field invariant and electromagnetic temperature-gradient instabilities in ‘good-curvature’ plasmas

Curvature-driven instabilities are ubiquitous in magnetised fusion plasmas. By analysing the conservation laws of the gyrokinetic system of equations, we demonstrate that the well-known spatial localisation of these instabilities to regions of ‘bad magnetic curvature’ can be explained using the conservation law for a sign-indefinite quadratic quantity that we call the gyrokinetic field invariant. Its evolution equation allows us to define the local effective magnetic curvature whose sign demarcates the regions of ‘good’ and ‘bad’ curvature, which, under some additional simplifying assumptions, can be shown to correspond to the inboard (high-field) and outboard (low-field) sides of a tokamak plasma, respectively. We find that, given some reasonable assumptions, electrostatic curvature-driven modes are always localised to the regions of bad magnetic curvature, regardless of the specific character of the instability. More importantly, we also deduce that any mode that is unstable in the region of good magnetic curvature must be electromagnetic in nature. As a concrete example, we present the magnetic-drift mode, a novel good-curvature electromagnetic instability, and compare its properties with the well-known electron-temperature-gradient instability. Finally, we discuss the relevance of the magnetic drift mode for high-β fusion plasmas, and in particular its relationship with microtearing modes.

fusion plasma↗
DOE OSTI2025

Saturation and decay of an induced scattering instability driven by static field asymmetries in a non-neutral plasma

Experimental data on the saturation and decay of a previously described induced scattering instability are presented. The instability is driven by static background field asymmetries and initially produces a single growing mode (ω 1 , k 1 ). As the instability develops, several harmonically related (ω n =nω 1 , k n =nk 1 ) eigenmodes of the plasma column grow up, resulting in a spatially localized moving potential perturbation with δϕ/T ≈ 20%. After a period of rapid growth, the mode amplitudes saturate and vary slowly for roughly 10 5 ω 1 −1 . This quasi-equilibrium ends with an unusual change in the eigenmode phases, wherein δϕ changes from positive to negative. Plasma heating during this period increases the damping of all modes and rapidly quenches the instability. This instability is of particular interest because the coupling with the field asymmetry produces an enhanced level of radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1966

Nonlinear Aspects of Combustion Instability in Liquid Propellant Rocket Motors Sixth Yearly Progress Report, 1 Jun. 1965 - 31 May 1966

This is the sixth yearly report on the research at Princeton University on the problem of combustion instability in liquid-propellant rocket motors. These investigations have been sponsored by the National Aeronautics and Space Administration under NASA Contract NASr-217 and, in the past, under NASA Grant NsG-99-60. The nonlinear aspects of the problem have been emphasized in both the theoretical and experimental portions of the program. Presently, the theoretical research primarily involves studies of possible mechanisms or forcing functions of the instabilities, the nonlinear gasdynamic oscillations associated with instability in rocket motors, the effect of the exhaust nozzle on the stability, and the effect of damping devices on the stability. In conjunction with the theoretical program, an experimental program of an equally broad range is being conducted. This program includes basic experiments both with rocket motors and "cold" chambers which have been designed to study various combustion processes in order to provide understanding which would be useful in the construction of a theoretical model. It also includes studies of the effects of exhaust nozzles and damping devices upon the instability as well as controlled testing of liquid rocket motors with parametric variations intended to determine the factors which strongly influence the incidence of nonlinear combustion instability.

Combustion instability↗
NASA NTRS2021

Boundary-layer Instability Measurements on a Cone at Freestream Mach 3.5

An experimental study was conducted in the NASA (National Aeronautics and Space Administration) Langley Supersonic Low-Disturbance Tunnel to investigate naturally-occurring instabilities in a supersonic boundary layer on a 7◦ half-angle cone at nominal freestream conditions: Mach 3.5, total temperature of 299.8K, and unit Reynolds numbers (millions per m) of 9.89, 13.85, 21.77, and 25.73. Instability measurements were acquired under noisy-flow and quiet-flow conditions. Pitot-pressure and calibrated hot-wire measurements were obtained using a model-integrated traverse system to document the model flow field. In noisy-flow conditions, growth rates and mode shapes achieved good agreement between the measured results and linear stability theory (LST). The corresponding N factor at transition from LST is N ≈ 3.9. Under quiet-flow conditions, the most unstable first-mode instabilities as predicted by LST were measured, but this mode was not the dominant instability measured in the boundary layer. Instead, the dominant instabilities were less-amplified, low-frequency disturbances predicted by LST, and grew according to linear theory. These low-frequency unstable disturbances were initiated by freestream acoustic disturbances through a receptivity process believed to occur near the branch I location of the cone. Under quiet-flow conditions, the boundary layer remained laminar up to the last measurement station for the largest unit Reynolds number, implying a transition N factor of N > 8.5.

Boundary Layer↗
DOE OSTI2024

Visualizing plasmons and ultrafast kinetic instabilities in laser-driven solids using X-ray scattering

Ultra-intense lasers that ionize atoms and accelerate electrons in solids to near the speed of light can lead to kinetic instabilities that alter the laser absorption and subsequent electron transport, isochoric heating, and ion acceleration. These instabilities can be difficult to characterize, but X-ray scattering at keV photon energies allows for their visualization with femtosecond temporal resolution on the few nanometer mesoscale. Here, we perform such experiment on laser-driven flat silicon membranes that shows the development of structure with a dominant scale of 60 nm in the plane of the laser axis and laser polarization, and 95 nm in the vertical direction with a growth rate faster than 0.1 fs -1 . Combining the XFEL experiments with simulations provides a complete picture of the structural evolution of ultra-fast laser-induced plasma density development, indicating the excitation of plasmons and a filamentation instability. Particle-in-cell simulations confirm that these signals are due to an oblique two-stream filamentation instability. These findings provide new insight into ultra-fast instability and heating processes in solids under extreme conditions at the nanometer level with possible implications for laser particle acceleration, inertial confinement fusion, and laboratory astrophysics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

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

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

Physics↗
DOE OSTI2025

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2026

Study of fully coupled three-dimensional envelope instability using automatic differentiation

Automatic differentiation is a powerful tool for computing derivatives of simulation results with respect to given parameters. In this Letter, we have applied this tool to investigate the instability of a dynamical system governed by 21 ordinary differential equations. This second-order instability (named envelope instability) is driven by space-charge effects and has a significant impact on the operational regimes of particle accelerators. Our study delves into the three-dimensional envelope instability, incorporating both transverse and longitudinal coupling. Conventionally, analyzing this complex system would necessitate solving 441 ordinary differential equations, which is computationally intractable. However, by employing automatic differentiation, we were able to track only 21 equations. This approach allowed us to uncover an additional instability stopband, which arises from space-charge-induced coupling and has not been reported in previous studies. This research highlights the significant advantages of automatic differentiation in analyzing complicated dynamical systems involving a large number of ordinary differential equations.

Qiang, Ji [Lawrence Berkeley National Laboratory (↗
DOE OSTI2026

Neutrino flavor instabilities in neutron star mergers with moment transport: Slow, fast, and collisional modes

Determining where, when, and how neutrino flavor oscillations must be included in large-scale simulations of hot and dense astrophysical environments is an enduring challenge that must be tackled to obtain accurate predictions. Here, using an angular moment-based linear stability analysis framework, we examine the different kinds of flavor instabilities that can take place in the context of the postprocessing of a neutron star merger simulation, with a particular focus on the collisional flavor instability and a careful assessment of several commonly used approximations. First, neglecting anisotropies of the neutrino field, we investigate the extent to which commonly used monoenergetic growth rates reproduce the results obtained from a full multienergy treatment. Contrary to the large discrepancies found in core-collapse supernova environments, we propose a simple combination of energy-averaged estimates that reproduces the multienergy growth rates in our representative simulation snapshot. We then quantify the impact of additional physical effects, including nuclear many-body corrections, scattering opacities, and the inclusion of the vacuum term in the neutrino Hamiltonian. Finally, we include the neutrino distribution anisotropies, which allows us to explore, for the first time in a multienergy setting, the interplay between collisional, fast, and slow modes in a moment-based neutron star merger simulation. We find that, despite a dominance of the fast instability in most of the simulation volume, certain regions exhibit only a collisional instability, while others, especially at large distances, exhibit a slow instability that is largely underestimated if anisotropic effects are neglected.

neutrino oscillations↗
DOE OSTI2026

Machine Protection Studies of Fast Instability-Driven Beam Losses in the EIC Electron Storage Ring

Fast transverse beam instabilities represent one of the most challenging protection scenarios for the Electron-Ion Collider (EIC) electron storage ring (ESR). Coherent betatron oscillations driven by collective effects can grow rapidly and lead to substantial beam losses within only a few tens of turns, posing a threat to collimators, superconducting (SC) magnets, cryogenic systems, and detector components. A dedicated simulation framework based on distributed transverse excitations was developed to evaluate instability-driven beam losses in the ESR. A conservative instability rise time on the order of 10 turns was used to define a bounding protection scenario. Simulations show that the original two-collimator configuration is insufficient to localize losses under all operating conditions, whereas an optimized four-collimator configuration reduces loss leakage into protected regions by more than two orders of magnitude. Nevertheless, the most severe instability scenarios produce catastrophic beam impacts, with nearly half of the stored beam being intercepted by the collimators within a single turn, exceeding the preliminary estimated survivability limits of the baseline collimator design. The results demonstrate that collimation alone is insufficient to guarantee safe ESR operation and must be complemented by a fast machine-protection system. Dedicated turn-by-turn beam-position monitors and fast beam-loss monitors provide sufficient advance warning to detect instability growth and initiate a beam abort before equipment-damaging loss levels are reached. These findings establish quantitative requirements for the ESR machine-protection architecture and support the implementation of a multilayer protection strategy for both accelerator and detector systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗
NASA NTRS1971

Ion-wave current instabilities and anomalous resistivity.

A theory of ion-wave current instabilities which takes into account, in a self-consistent manner, the inhomogeneities generated by field-aligned currents in a collisionless plasma is presented. Diamagnetic current associated with the current-produced density gradient is included in the distribution. The theory predicts that for a given frequency, the current threshold for ion-wave current instabilities is, in general, much below the threshold of ion acoustic instability in a uniform plasma as given by Fried and Gould. The current threshold is essentially zero in the limit when the ion Landau damping effect is negligible, or equivalently, the ion-wave current instabilities are absolute if there are no limitations on the wavelength. This is true even in the absence of externally applied density gradients in contrast to Kadomtsev's drift-wave results. For dimensions of interest in laboratory plasmas, the predicted linear growth rate increases with increasing longitudinal wavelength and with decreasing wavelength parallel to the diamagnetic current. Under the conditions of an experiment on anomalous resistivity (the dimensions but not the geometry of the machine have been considered), there is good agreement between the predicted onset of ion-wave current instabilities and the experimental data on the onset of anomalous resistivity.

Kan, J. R.↗