Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unstable 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

FICH: Feature instability cavitation history [Slides]

We seek to improve our ability to predict and understand the Richtmyer-Meshkov (RM) unstable dynamics of explosively shocked isolated features interfaced with vacuum or gases, and an HE products cavity. Explosively shocked features form as 2D or 3D RM phenomena. Liquid cavitation cause cavitation bubbles when shockwaves collide, and cavitation waves from wave interactions with isolated features form that source mass into the RM unstable sheets or spikes, and the cavitation waves may slice through bulk material into the HE products cavity.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Electron-scale Magnetic Holes Generation Driven by Whistler-to-Bernstein Mode Conversion in Fully Kinetic Plasma Turbulence

Magnetic holes (MHs) are coherent structures characterized by a strong and localized magnetic field amplitude dip, commonly observed in the heliosphere. These structures come in different sizes, from magnetohydrodynamic to kinetic scales. Subion-scale MHs are usually sustained by an electron current vortex and exhibit a strong electron temperature anisotropy, with higher temperatures perpendicular to the background magnetic field. Magnetospheric multiscale observations (MMSs) have revealed electron-scale MHs to be ubiquitous in the turbulent Earth’s magnetosheath and the solar wind, potentially playing an important role in the energy cascade and dissipation. Despite abundant observations, the origin of electron-scale MHs is still unclear and debated. In this work, we use fully kinetic simulations to investigate the role of plasma turbulence in generating electron-scale MHs. We find that the turbulence spontaneously produces electron-scale MHs via the following mechanism: first, large-scale turbulent velocity shears produce regions with high electron temperature anisotropy; these localized regions become unstable, generating oblique electron-scale whistler waves; as they propagate over the inhomogeneous turbulent background, whistler fluctuations develop an electrostatic component, turning into Bernstein-like modes; the strong electrostatic fluctuations produce current filaments that merge into an electron-scale current vortex; the resulting electron vortex locally reduces the magnetic field amplitude, finally evolving into an electron-scale MH. We show that MHs generated by this mechanism have properties consistent with MMSs and nontrivial kinetic features with a “mushroom”-shaped electron velocity distribution function. Our results have potential implications for understanding the formation and occurrence of electron-scale MHs in astrophysical turbulent and space environments, such as the Earth’s magnetosheath and the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Linear shaped-charge jet optimization using machine learning methods

Linear shaped charges are used to focus energy into rapidly creating a deep linear incision. The general design of a shaped charge involves detonating a confined mass of high explosive (HE) with a metal-lined concave cavity on one side to produce a high velocity jet for the purpose of striking and penetrating a given material target. This jetting effect occurs due to the interaction of the detonation wave with the cavity geometry, which produces an unstable fluid phenomenon known as the Richtmyer–Meshkov instability and results in the rapid growth of a long narrow jet. We apply machine learning and optimization methods to hydrodynamics simulations of linear shaped charges to improve the simulated jet characteristics. The designs that we propose and investigate in this work generally involve modifying the behavior of the detonation waves prior to interaction with the liner material. These designs include the placement of multiple detonators and the use of metal inclusions within the HE. In conclusion, we are able to produce a linear shaped-charge design with a higher penetration depth than the baseline case that we consider and accomplish this using the same amount of or less HE.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mitigation of resistive drift wave and ion temperature gradient instabilities by velocity shear

The effects of velocity shear on the resistive drift wave instability in the non-adiabatic limit and the toroidal ion temperature gradient instability are investigated for a plasma of inhomogeneous density/ion temperature, respectively. For the resistive drift wave, we find that the instability growth rate decreases monotonically with increasing magnitude of shear, but we find that complete stabilization is impossible. For the ion temperature gradient instability, we find that the standard WKB approximation is insufficient to describe the full behavior of the instability, and that an analysis of the localized eigenmode problem reveals two separate unstable solutions, which the WKB approximation does not predict. The impact of flow shear on these two new unstable solutions is discussed. In both resistive drift wave and ion temperature gradient instabilities, the sheared flow causes a shifting, tilting, and sharpening of the electrostatic potential eddies.

Physics↗
DOE OSTI2023

Stability of multi-population traffic flows

Traffic waves, known also as stop-and-go waves or phantom jams, appear naturally as traffic instabilities, also in confined environments as a ring-road. A multi-population traffic is studied on a ring-road, comprised of drivers with stable and unstable behavior. There exists a critical penetration rate of stable vehicles above which the system is stable, and under which the system is unstable. In the latter case, stop-and-go waves appear, provided enough cars are on the road. The critical penetration rate is explicitly computable, and, in reasonable situations, a small minority of aggressive drivers is enough to destabilize an otherwise very stable flow. This is a source of instability that a single population model would not be able to explain. Also, the multi-population system can be stable below the critical penetration rate if the number of cars is sufficiently small. Instability emerges as the number of cars increases, even if the traffic density remains the same (i.e., number of cars and road size increase similarly). This shows that small experiments could lead to deducing imprecise stability conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Interpreting and Stabilizing Machine-Learning Parametrizations of Convection

Neural networks are a promising technique for parameterizing subgrid-scale physics (e.g., moist atmospheric convection) in coarse-resolution climate models, but their lack of interpretability and reliability prevents widespread adoption. For instance, it is not fully understood why neural network parameterizations often cause dramatic instability when coupled to atmospheric fluid dynamics. This paper introduces tools for interpreting their behavior that are customized to the parameterization task. First, we assess the nonlinear sensitivity of a neural network to lower-tropospheric stability and the midtropospheric moisture, two widely studied controls of moist convection. Second, we couple the linearized response functions of these neural networks to simplified gravity wave dynamics, and analytically diagnose the corresponding phase speeds, growth rates, wavelengths, and spatial structures. To demonstrate their versatility, these techniques are tested on two sets of neural networks, one trained with a super-parameterized version of the Community Atmosphere Model (SPCAM) and the second with a near-global cloud-resolving model (GCRM). Additionally, even though the SPCAM simulation has a warmer climate than the cloud-resolving model, both neural networks predict stronger heating/drying in moist and unstable environments, which is consistent with observations. Moreover, the spectral analysis can predict that instability occurs when GCMs are coupled to networks that support gravity waves that are unstable and have phase speeds larger than 5 m s -1 . In contrast, standing unstable modes do not cause catastrophic instability. Using these tools, differences between the SPCAM-trained versus GCRM-trained neural networks are analyzed, and strategies to incrementally improve both of their coupled online performance unveiled.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Stability Criteria for Self-Propagating Reaction Waves in Co/Al Multilayers

The propagation of self-sustained formation reactions in sputter-deposited Co/Al multilayers is known to exhibit a design-dependent instability. Multilayers having thin bilayers (<55 nm period) exhibit stable propagating waves, whereas those with a larger period react unstably. The specific two-dimensional (2D) instability observed involves the transverse propagation of a band in front of a stalled front commonly referred to as a “spin band.” Previous finite-element studies have shown that these instabilities are thermodynamically driven by the forward conduction of heat away from the flame front. However, the magnitude of that loss is inherently tied to the bilayer design in traditional bimetallic multilayers, which couples any proposed stability criteria to a varying critical diffusion distance. Here, this work utilizes a recently developed class of materials known as “inert-mediated reactive multilayers” to decouple the thermodynamic and kinetic contributions to propagating wave stability by reducing the stored chemical energy density in normally stable bilayer designs. By depositing an inert product phase (B2-CoAl) within the mid-plane of Co and Al reactant layers, spin instabilities arise as a function of both diluted volume and critical diffusion distance. From there, a stability criterion is determined for Co/Al multilayers based on enthalpy loss from the reaction zone, and its physical significance is explored.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗
DOE OSTI2023

Quantum three-wave instability

For the three-wave interaction, the lowest-order nonlinear interaction in plasma dynamics, we describe how the quantum system, which is time-independent, nonchaotic, finite-dimensional, and Hermitian, can give rise to a linear instability corresponding to that in the classical system. Furthermore, we show that the instability is realized in the quantum regime as a cascade of the wave function in the space of occupation number states, and the unstable quantum system has a richer spectrum and a much longer recurrence time than the stable quantum system. The conditions for instability of the quantum three-wave interaction are described.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Synergetic Ferroelectricity and Superconductivity in Zero-Density Dirac Semimetals near Quantum Criticality

Here, we study superconductivity in a three-dimensional zero-density Dirac semimetal in proximity to a ferroelectric quantum critical point. We find that the interplay of criticality, inversion-symmetry breaking, and Dirac dispersion gives rise to a robust superconducting state at the charge-neutrality point, where no Fermi surface is present. Using Eliashberg theory, we show that the ferroelectric quantum critical point is unstable against the formation of a ferroelectric density wave (FDW), whose fluctuations, in turn, lead to a first-order superconducting transition. Surprisingly, long-range superconducting and FDW orders are found to cooperate with each other, in contrast to the more usual scenario of phase competition. Therefore, we suggest that driving charge neutral Dirac materials, e.g., Pb x Sn 1–x Te, through a ferroelectric quantum critical point may lead to superconductivity intertwined with FDW order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Consideration of Nonideal Detonation Regimes Influenced by Wave Modes in a Water-Cooled Rotating Detonation Engine Using OH* Chemiluminescence

Although inherently unstable, existing research in rotating detonation combustion supports its application in notionally steady processes resulting in greater availability compared to conventional, constant pressure combustion. Further improvements rely on a more in-depth understanding of system losses and identifying conditions which optimize device performance. Within this study, the presence and proportion of ideal and nonideal combustion regimes are compared across a variety of process conditions and wave modes. Here, large-scale data analysis seeks to summarize proportional heat release associated with commensal, parasitic, and detonative combustion averaged across individual traces of OH* chemiluminescent data acquired at the detonation plane. Means of regime partitioning based on the anatomy of the time-resolved OH* signal are proposed to ensure consistent analysis throughout the current and future studies concerning combustion regimes. Of particular interest is the possible influence of wave on the nonideal combustion in relative proportion to the desired detonation. Results showed improved percent detonation with increasing significance for the following trends: decreasing equivalence ratio, increasing wave count, decreasing wave velocity, and increasing detonation time. Increased wave number, brought on by decreased equivalence ratios and wave velocities, is thought to decrease fill region surface area, and therefore, decrease nonideal contact burning. Additional analysis is performed to consider possible trend variation due to the presence of stable galloping waves, which were found to have minimal influence on relative percent detonation values. The outcome of this study suggests operational states, which correspond to increased wave quantities for increased proportions of reactants consumed by the targeted detonative combustion regime.

42 ENGINEERING↗
DOE OSTI2026

Stability and Characteristics of Lower-hybrid Drift Waves: Dependence on Electron Beta and Cross-field Relative Drift

Lower-hybrid drift waves (LHDWs) are frequently observed microinstabilities in both space and laboratory plasmas. Despite decades of study, the relationship between electrostatic (ES-LHDW) and electromagnetic (EM-LHDW) variants and the plasma parameters controlling their stability remains unclear. Here, we systematically examine LHDW behavior by solving the local linear dispersion relation over a wide range of plasma and field conditions. Our results demonstrate that ES-LHDWs and EM-LHDWs are not distinct modes but rather two different regimes of the same drift wave whose character evolves smoothly with electron beta (β e ) and the cross-field electron drift velocity relative to ions, normalized to the ion sound speed (u 0x /C s ). The nature of the waves changes from electrostatic to electromagnetic when β e increases. Growth rates increase with u 0x /C s but decrease with β e , while the most unstable wavelength remains nearly universal, with kρ e ∼ 0.8 (k is the magnitude of the wave vector and ρ e is the electron gyroradius). We further present quasi-linear estimates of nonlinear saturation properties, including energy partition among electric fields, magnetic fields, and particle kinetic responses. We show that ES-LHDWs reach higher electric-field saturation amplitudes, whereas EM-LHDWs generate strong magnetic perturbations and parallel electric fields that may enable efficient particle heating. Comparisons with the classical model reveal that retaining electromagnetic effects is essential for accurate predictions of frequency, growth rate, and the propagation angle. These findings provide a unified framework for understanding LHDWs across diverse collisionless plasma environments, including current sheets of magnetic reconnection, shear layers, collisionless shocks, and boundary regions.

Solar coronal waves↗
DOE OSTI2024

A new method for solving the linearized 1D Vlasov–Poisson system yielding a new class of solutions

We describe a new method for solving the linearized 1D Vlasov–Poisson system by using properties of Cauchy-type integrals. Our method remedies critical flaws of the two standard methods, reveals a previously unrecognized Gaussian-in-time-like decay, and can also account for an externally applied electric field. The Landau approximation involves deforming the Bromwich contour around the poles closest to the real axis due to the analytically continued dielectric function, finding the long-time behavior for a stable system: Landau damping. Jackson's generalization encircles all poles while sending the contour to infinity, assuming its contribution vanishes, which is not true in general. This gives incorrect solutions for physically reasonable configurations and can exhibit pathological behavior, of which we show examples. The van Kampen method expresses the solution for a stable equilibrium as a continuous superposition of waves, resulting in an opaque integral. Case's generalization includes unstable systems and predicts a decaying discrete mode for each growing discrete mode, an apparent contradiction to both the Jackson solution and ours. We show, without imposing additional constraints, that the decaying modes are never present in the time evolution due to an exact cancellation with part of the continuum. Our solution is free of integral expressions, is obtained using algebra and Laurent series expansions, does not rely on analytic continuations, and results in a correct asymptotically convergent form in the case of infinite sums. The analysis used can be readily applied in higher-dimensional, electromagnetic systems and also provides a new technique for evaluating certain inverse Laplace transforms.

Physics↗
DOE OSTI2020

Experiments and Simulations on the Turbulent, Rarefaction Wave Driven Rayleigh–Taylor Instability

Experiments were performed to observe the growth of the turbulent, Rayleigh–Taylor unstable mixing layer generated between air and SF 6 , with an Atwood number of A=(ρ 2 –ρ 1 )/(ρ 2 +ρ 1 ) = 0.64, where ρ 1 and ρ 2 are the densities of air and SF 6 , respectively. A nonconstant acceleration with an average value of 2300g 0 , where g 0 is the acceleration due to gravity, was generated by interaction of the interface between the two gases with a rarefaction wave. Three-dimensional, multimode perturbations were generated on the diffuse interface, with a diffusion layer thickness of δ = 3.6 mm, using a membraneless vertical oscillation technique, and 20 experiments were performed to establish a statistical ensemble. The average perturbation from this ensemble was extracted and used as input for a numerical simulation using the Lawrence Livermore National Laboratory (LLNL) Miranda code. Good qualitative agreement between the experiment and simulation was observed, while quantitative agreement was best at early to intermediate times. Several methods were used to extract the turbulent growth constant α from experiments and simulations while accounting for time varying acceleration. Experimental, average bubble and spike asymptotic self-similar growth rate values range from α = 0.022 to α = 0.032 depending on the method used, and accounting for variable acceleration. Values found from the simulations range from α = 0.024 to α = 0.041. Here, the values of α measured in the experiments are lower than what are typically measured in the literature but are more in line with those found in recent simulations.

42 ENGINEERING↗
DOE OSTI2025

Analysis of runaway electron driven whistler wave instability experiments

Data acquired on the DIII-D tokamak were analyzed. The data are from experiments that were conducted to study an instability that is driven unstable by intense populations of electrons with MeV energies that are known as runaway electrons. The instability is a type of plasma wave called a whistler wave that occurs at frequencies above the ion cyclotron frequency but well below the electron cyclotron frequency. The waves were measured by magnetic fluctuation coils that are embedded in the DIII-D vacuum vessel wall. After upgrades to this diagnostic were completed, new experiments were conducted on July 13, 2020 in order to measure the toroidal mode number of the whistler waves and to extend the frequency of the detected waves. Through the use of mixers, instability between 600-700 MHz was detected. (The initial experiments only measured up to 200 MHz.) Analysis of the data was led by Hari Choudhury, a PhD student at Columbia University. Mr. Choudhury has submitted two papers for publication that include contributions by UC Irvine (UCI) Professor Heidbrink and his graduated PhD student Genevieve DeGrandchamp: “Detailed Characterization of Runaway Electron Driven Whistler Waves in Low-Density DIII-D Discharges” and “First Demonstration of Resonant Pitch-Angle Scattering of Relativistic Electrons by Externally-Launched Helicon Waves.” The first paper, which has been submitted to Physics of Plasmas, has significant contributions to both the data and the interpretation by UCI scientists. In contrast, UCI contributions to the second paper, which has been submitted to Physical Review Letters, are relatively minor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Positive Energy Functional for Massless Scalars in Rotating Black Hole Backgrounds of Maximal Ungauged Supergravity

Here, we outline a proof of the stability of a massless neutral scalar field ψ in the background of a wide class of four dimensional asymptotically flat rotating and “electrically charged” solutions of supergravity, and the low energy limit of string theory, known as STU metrics. Despite their complexity, we find it possible to circumvent the difficulties presented by the existence of ergo regions and the related phenomenon of superradiance in the original metrics by following a strategy due to Whiting, and passing to an auxiliary metric admitting an everywhere lightlike Killing field and constructing a scalar field ψ (related to a possible unstable mode ψ by a nonlocal transformation) which satisfies the massless wave equation with respect to the auxiliary metric. By contrast with the case for ψ, the associated energy density of ψ is not only conserved but is also non-negative.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions

The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.

Horiuchi, W. [Osaka Metropolitan University (Japan↗