Search NASA⌕ Search

SEARCH · Search NASA

Results for “Instability waves”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High-resolution particle-in-cell simulations of two-dimensional Bernstein–Greene–Kruskal modes

We present two-dimensional (2D) particle-in-cell (PIC) simulations of 2D Bernstein–Greene–Kruskal modes, which are exact nonlinear steady-state solutions of the Vlasov–Poisson equations, on a 2D plane perpendicular to a background magnetic field, with a cylindrically symmetric electric potential localized on the plane. PIC simulations are initialized using analytic electron distributions and electric potentials from the theory. We confirm the validity of such solutions using high-resolutions up to a 20482 grid. We show that the solutions are dynamically stable for a stronger background magnetic field, while keeping other parameters of the model fixed, but become unstable when the field strength is weaker than a certain value. When a mode becomes unstable, we observe that the instability begins with the excitation of azimuthal electrostatic waves that ends with a spiral pattern.

Physics↗
DOE OSTI2024

An extended variational method for the resistive wall mode in toroidal plasma confinement devices

The external-kink stability of a toroidal plasma surrounded by a rigid resistive wall is investigated. The well-known analysis of Haney and Freidberg is rigorously extended to allow for a wall that is sufficiently thick that the thin-shell approximation does not necessarily hold. A generalized Haney–Freidberg formula for the growth-rate of the resistive wall mode is obtained. Thick-wall effects do not change the marginal stability point of the mode but introduce an interesting asymmetry between growing and decaying modes. Growing modes have growth-rates that exceed those predicted by the original Haney–Freidberg formula. On the other hand, decaying modes have decay-rates that are less than those predicted by the original formula. The well-known Hu–Betti formula for the rotational stabilization of the resistive wall mode is also generalized to take thick-wall effects into account. Increasing wall thickness facilitates the rotational stabilization of the mode, because it decreases the critical toroidal electromagnetic torque that the wall must exert on the plasma. On the other hand, the real frequency of the mode at the marginal stability point increases with increasing wall thickness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The radial phase variation of reversed-shear and toroidicity-induced Alfvén eigenmodes in DIII-D

The eigenfunction of an instability contains information about energy flow in the wave. Here, the amplitude and phase of electron cyclotron emission radiometer data from hundreds of DIII-D reversed shear Alfvén eigenmodes (RSAE) and toroidicity-induced Alfvén eigenmodes (TAE) are analyzed along the outboard horizontal midplane. The radial phase profile can be flat, linearly rising or falling, convex or concave; in other words, a wide variety of shapes is observed. For a particular mode, often the radial phase profile remains approximately constant as the mode evolves in time but sometimes it changes rapidly. Many TAEs and some RSAEs have phase profiles that are rather flat where the mode amplitude is largest but rise steadily by ~2π at large major radius. Rapid phase changes are observed when the frequencies of an RSAE and TAE overlap and the modes couple. The phase profile depends weakly on the fast-ion gradient that would appear in the absence of wave-induced transport. Linear and quadratic fits to the phase profiles, together with many plasma parameters, are assembled into RSAE and TAE databases. In both cases, large variability is observed. For RSAEs, the strongest phase dependencies are on electron temperature T e , RSAE mode frequency, and the density of carbon impurities. For TAEs, the strongest dependencies are on beam power and major radius of the mode. In general, the average RSAE radial phase profile is essentially flat but the TAE profile has non-zero slope and curvature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Formation and Diagnostics of Collisionless Shocks in Laser-Produced Plasmas (Final Report)

Collisions of supersonic flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large; therefore, most shock waves in astrophysics are “collisionless” and form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at DOE high energy density (HED) facilities can achieve the conditions necessary for the formation of collisionless shocks and will provide a unique avenue for studying the nonlinear physics of shock waves. This project investigates the parameters and experimental configurations needed for the creation of collisionless shocks in the laboratory experiments. This study is being performed with state of the art 3D kinetic plasma simulations. Recent advances in ab-initio particle-in-cell (PIC) simulations of astrophysical shocks enable us to simulate shock formation on timescales and in geometries relevant to the laboratory experiments. The main goal of the proposal is to ascertain which signatures of shock physics can be extracted from the realistic experiments, accounting for limitations imposed by finite spatial and time scales in the lab. To aid the interpretation of current and future experimental data, PIC simulation data have been analyzed using the same diagnostic techniques as used in the laboratory, such as proton radiography and neutronics, to unambiguously predict the experimental signatures of relevant collisionless shock phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Collaborative Research: Particle Simulation of Collisionless Magnetic Reconnection Under Finite Guide Field

Magnetic reconnection is believed to be one of the most fundamental physical processes in laboratory and space plasmas. Fundamental physics of collisionless reconnection under general conditions of finite guide magnetic field, realistic electron-to-ion mass ratio, and in 3D geometries is still not adequately understood. The major goal of this project is to study the physics of collisionless magnetic reconnection in a current sheet under a finite guide field and with a realistic ion-to-electron mass ratio mi/me. Under the support of this DoE grant, we have made significant progress in (1) the investigation of linear and nonlinear current sheet instabilities under a wide range of guide field using the 3-D gyro kinetic electron and fully kinetic ion (GeFi) particle simulation model and (2) the investigation of fully nonlinear reconnection process with a strong guide field with the GeFi simulations. Furthermore, (3) we have extended the study to the realistic geometry of the magnetopause using 3D global hybrid simulations. This study, from local to global scales and for 3-D physics at plasma boundaries, has provides advanced theoretical understanding of the fundamental physics of magnetic reconnection in laboratory and space plasma regimes on multi-scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Limit Cycle Oscillating Detonation Wave Behavior Analysis Within a Rotating Detonation Engine

Limit cycle oscillation (LCO) detonation wave behaviors are presented and analyzed for test times exceeding 20 s in a water-cooled rotating detonation engine (RDE). LCO detonation waves exhibit cyclic acceleration and deceleration, resulting in oscillating wave spacing at unique process conditions. In previous RDE studies, similar behaviors have been studied as microsecond-scale instabilities leading to ascending or descending modal transitions. In the current work, however, LCO waves are considered a persistent wave mode, occupying unique portions of the operational envelope adjacent to those of their equally spaced counterparts. These occurrences of LCO waves are repeatable, enduring behaviors. A method to generate shifted contour surfaces specifically intended to extract and analyze wave spacing variation through time, termed limit cycle oscillation visualization (LCOV) surfaces, is presented. LCOV surfaces transform data into the reference frame of a primary traveling wave and are used to analyze quasi-steady, short-timescale, and transitional LCO modes. Results are leveraged to understand the relationship between fill height, wave strength, local wave acceleration, and subsequent LCO wave spacing for individual wave sets. In conclusion, quasi-steady LCO waves display wave spacing oscillations between equal spacing values associated with ±1 wave across runs exceeding 18 s.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Experimental Investigation of Interactions Between Two Closely Spaced Azimuthal Modes in a Multinozzle Can Combustor

Thermoacoustic instabilities in annular or circular combustors are often coupled with azimuthal modes. These modes can be characterized by a spin ratio (SR), which quantifies the dominant mode between two counter-rotating waves, and a phase difference between them, which is directly related to the orientation of the antinodal line. This study investigates the instability amplitudes, SR, and phase difference of two closely spaced (3% of their mean frequency), yet distinct azimuthal modes; one is the first azimuthal (1A) mode, and the other is a combination of first azimuthal and first longitudinal (1A1L) mode. Each mode itself consists of two peaks that are spaced even more closely in frequency (0.8%). Furthermore, distinct harmonics at 2× and 3× of these frequencies, presumably associated with nonlinearities, are also evident in the spectra. Each mode is bandpass filtered in spectrum to analyze them separately. For the 1A mode, the SR and phase difference exhibit a variety of behaviors—including quasi-periodic standing waves, spinning waves, and intermittency—depending on operating conditions such as thermal power and azimuthal fuel staging. Similar trends are observed for the 1A1 L mode. Moreover, there is clear coupling between the 1A and 1A1 L modes, as their SRs are almost synchronized during the quasi-periodic standing wave. This synchronization is observed in phase differences as well, but not in the instability amplitude. For spinning dominant wave conditions, the SRs of each mode have similar average values, but they fluctuate in a seemingly random fashion. For the phase difference, both average and fluctuation are not correlated. In contrast, the instability amplitudes are strongly correlated, with modulation of the 1A mode leading to that of the 1A1 L mode. Furthermore, these results clearly indicate that complex coupling occurs across closely spaced frequencies under instability conditions, coupling that must be understood in order to capture limit cycle dynamics.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗
DOE OSTI2021

Wave-kinetic approach to zonal-flow dynamics: Recent advances

The basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of the wave-kinetic theory. Recently, this framework has been reexamined from first principles, which has led to more accurate yet still tractable “improved” wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow “curvature” (the second radial derivative of the flow velocity) on the dynamics and stability of drift waves and zonal flows. Here, we overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Role of ion acoustic instability in magnetic reconnection

We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Atomic structure and formation mechanism of a newly discovered charge density wave in the m = 2 monophosphate tungsten bronze

The m = 2 member of the monophosphate tungsten bronze family had been thought to be the only one without an electronic instability at low temperature. In this paper, we report the discovery of a charge density wave phase in this compound, with a transition temperature of 290 K and an incommensurate modulation vector q = 0.245b* + ξc*, which reaches a lock-in with a commensurate vector at 130 K. The presence of this new phase is confirmed by diffraction and resistivity measurements. Pre-transitional dynamics are investigated using diffuse and inelastic X-ray scattering, revealing a clear Kohn anomaly. We analyze both structural and electronic contributions to the phase transition, providing a comprehensive picture of the mechanism driving this newly identified instability.

Minelli, Arianna [ORNL] (ORCID:0000000283103712)↗
DOE OSTI2021

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗
DOE OSTI2022

Revealing the Drag Instability in One-fluid Nonideal Magnetohydrodynamic Simulations of a 1D Isothermal C-shock

C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in nonideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the post-shock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Effects of the Wavelength of the Plasma Waves on Cross-Field Electron Transport in Partially Magnetized Plasmas

In this article, the effects of increasing plasma density on the cross-field electron transport driven by kinetic instabilities are investigated. A 2-D collisionless particle-in-cell (PIC) simulation is performed on cases with only singly charged ions and with a mixture of singly and doubly charged ions, at different levels of plasma density. It is observed using the PIC simulations that the level of plasma density affects the plasma waves that are driven by the microinstabilities, including the electron cyclotron drift instability (ECDI) and ion–ion two-stream instability (IITSI). Here, the simulation results suggest that the enhancement of cross-field electron transport is correlated with the excitation of smaller wavelength modes and the increased amplitude of the plasma waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effects of density turbulence on helicon wave propagation in the core plasmas

Radio frequency wave propagation can be significantly affected by density irregularities, such as filaments in the scrape-off layer or instabilities in the core plasma. In this study, we examine the impact of edge turbulence on helicon wave propagation using the Petra-M simulation tool. To analyze the effect of edge turbulence, we utilize a realistic background plasma derived from XGC simulations, which includes spatial density fluctuations at the edge. This focus is particularly relevant for the DIII-D configuration, characterized by edge density turbulence of the core plasmas in a wide pedestal QH-mode. We focus on helicon wave propagation in the core plasma since the slow mode cannot propagate into the core plasma when the density is higher than the lower hybrid resonance, where a wide pedestal QH-mode occurs. The simulation results indicate that edge density fluctuations have a substantial impact on helicon wave coupling due to scattering in the core. We specifically demonstrate that the toroidal mode number and the level of density fluctuations are significant factors influencing scattering. While we have minimized the excitation of slow modes from the antenna, we still observe mode-converted slow modes resulting from the incoming helicon waves. The insights derived from these simulations will inform upcoming tokamak experiments regarding helicon antenna coupling in long pulse scenarios.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗
DOE OSTI2021

Acoustic velocities, elasticity, and pressure-induced elastic softening in compressed neodymium

For the first time, the acoustic velocities and elasticity of double hexagonal close-packed neodymium (Nd) were measured at high pressure using ultrasonic interferometry in conjunction with synchrotron X-ray techniques. Here, both the compressional and shear wave velocities increased with an increase in pressure of up to ~5.4 GPa before they exhibited almost constant values upon further compression. This pressure-induced elastic wave velocity softening can be attributed to the structural instability in dhcp-Nd. From the current sound velocities and density data, we obtained new bulk and shear moduli as well as their pressure dependence for polycrystalline dhcp-Nd, yielding $B_0$=33.2(9) GPa, $G_0$=14.1(5) GPa, $\partial$$B/$$\partial$$P$ = 3.6 (4) and $\partial$$=G/$$\partial$$P$ = 1.3 (2) while other elasticity-related/mechanical properties, including Young's modulus, Debye temperature, Poisson's ratio, and Grüneisen parameter, which are important parameters for the application of Nd in extreme environments, were also derived.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Altermagnetism from coincident Van Hove singularities: application to κ -Cl

Realizing two-dimensional (2D) altermagnets is important for spintronics applications. Here we propose a microscopic template for stabilizing 2D altermagnetism through Van Hove singularities that are coincident in both energy and momentum. These coincident Van Hove singularities are a generic consequence of non-symmorphic symmetries in nine 2D space groups. Due to nontrivial symmetry properties of the Hamiltonian, these coincident Van Hove singularities allow new hopping interactions between the Van Hove singularities that do not appear in analogous Van Hove singularity based patch models for cuprates and graphene. We show these new interactions can give rise to various weak coupling, and BCS-based instabilities, including altermagnetism, nematicity, inter-band d-wave superconductivity, and orbital altermagnetic order. We apply our results to quasi-2D organic κ -Cl in which altermagnetism is known to appear.

electronic properties and materials↗
DOE OSTI2025

Synoptic and Mesoscale Modulation of Dynamic and Thermodynamic Impacts on Central Arctic Sea Ice During MOSAiC

The primary objectives of this proposed project are to determine 1) to what extent dynamic and thermodynamic air‐ice interactions in the Central Arctic are driven by atmospheric synoptic and mesoscale phenomena, and 2) the key processes involved in these interactions. These phenomena include large‐scale features such as fronts, high‐ and low‐pressure systems, low‐level jets, and cloud bands, but may also include smaller‐scale phenomena such as mesoscale instabilities (e.g., symmetric instability or slantwise convection), transient meso‐lows and gravity waves. Dynamic air‐ice interactions include direct atmospheric stress and atmospheric stress divergence/deformation; thermodynamic air‐ice interactions involve the net surface energy budget, especially in association with non‐stratocumulus clouds generally occurring with stronger dynamical forcing.

58 GEOSCIENCES↗