Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow instability”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Nature and the nonlinear evolution of electrostatic waves associated with the AMPTE solar wind releases

The nonlinear evolution of the electrostatic waves observed during the AMPTE (Active Magnetosphere Particle Tracer Explorers) solar wind releases is investigated. Previous linear studies indicate two distinct sets of instabilities may be responsible for the generation. One set consists of ion-acoustic type instabilities which are insensitive to the presence of a background magnetic field, while the other group corresponds to the modified two stream instabilities and requires the solar wind flow to be across the ambient magnetic field. To establish which set of instabilities is more viable a detailed linear Vlasov theory has been conducted by numerically solving the full electromagnetic dispersion relation. In addition both the plasma wave as well as the magnetic field measurements by the IRM (Ion Release Module) spacecraft were used to correlate the frequency and the power of the observed waves with the magnitude and the direction of the solar wind magnetic field. The results of these analyses indicate that the ion-acoustic type instabilities have growth rates that are an order of magnitude or more larger than those of the modified two stream instabilities. Results show that both the solar wind protons and the released ions may be heated and accelerated in the directions oblique to the solar wind flow velocity.

Omidi, N.↗
NASA NTRS1975

Further study of the atmospheric gravity waves over the eastern seaboard on 18 March 1969

Washington, D.C., microbarograph records for Mar. 18, 1969 reveal gravity-wave-associated pressure oscillations which appear to be directly related to upper tropospheric wave structure observed at the same time with a Wallops Island 10-cm wavelength radar. The consistency between the two sets of data provides observational support for the hypothesis that shear instability in the upper tropospheric flow is the mechanism responsible for the generation of such waves. The comparison presented suggests that microbarograph arrays might be useful adjuncts to future radar studies of upper tropospheric wave dynamics, supplying such wave parameters as phase velocity and wavelength in favorable cases. The radar data pertinent to this event show an apparent vertical wave phase variation, permitting a very approximate and somewhat uncertain estimate of the wave-associated vertical flux of horizontal momentum, which is found to be about 4 dyn per sq cm. Wave generation by shear instability in the upper tropospheric air flow and the resulting vertical momentum transport may be an important element of the global atmospheric momentum budget.

Hooke, W. H.↗
NASA NTRS2010

Development and Breakdown of Goertler Vortices in High Speed Boundary Layers

The nonlinear development of G rtler instability over a concave surface gives rise to a highly distorted stationary flow in the boundary layer that has strong velocity gradients in both spanwise and wall-normal directions. This distorted flow is susceptible to strong, high frequency secondary instability that leads to the onset of transition. For high Mach number flows, the boundary layer is also subject to the second mode instability. The nonlinear development of G rtler vortices and the ensuing growth and breakdown of secondary instability, the G rtler vortex interactions with second mode instabilities as well as oblique second mode interactions are examined in the context of both internal and external hypersonic configurations using nonlinear parabolized stability equations, 2-D eigenvalue analysis and direct numerical simulation. For G rtler vortex development inside the Purdue Mach 6 Ludwieg tube wind tunnel, multiple families of unstable secondary eigenmodes are identified and their linear and nonlinear evolution is examined. The computation of secondary instability is continued past the onset of transition to elucidate the physical mechanisms underlying the laminar breakdown process. Nonlinear breakdown scenarios associated with transition over a Mach 6 compression cone configuration are also explored.

Li, Fei↗
DOE OSTI2022

Local Lagrangian reduced-order modeling for the Rayleigh-Taylor instability by solution manifold decomposition

The Rayleigh-Taylor instability is a classical hydrodynamic instability of great interest in various disciplines of science and engineering, including astrophysics, atmospheric sciences and climate, geophysics, and fusion energy. Analytical methods cannot be applied to explain the long-time behavior of the Rayleigh-Taylor instability, and therefore, numerical simulation of the full problem is required. However, in order to capture the growth of amplitude of perturbations accurately, both the spatial and temporal discretizations need to be extremely fine for traditional numerical methods, and long-time simulation may become prohibitively expensive. In this paper, we propose efficient reduced order model techniques to accelerate the simulation of the Rayleigh-Taylor instability in compressible gas dynamics. Here, we introduce a general framework for decomposing the solution manifold to construct the temporal domain partition and temporally-local reduced order model construction with varying Atwood number. We propose two practical approaches in this framework, namely decomposition by physical time and by penetration distance. Numerical results are presented to examine the performance of the proposed approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Numerical simulations of onset and growth of Rayleigh–Taylor instability involving solids in converging geometry

Numerical simulation of Rayleigh–Taylor instabilities involving solids has been carried out in a spherically converging geometry. A solid shell is driven toward the center by high-pressure gas outside. Low-pressure gas inside the shell is then compressed, causing the solid shell to slow down and eventually bounce back. Here, simulation results show the absence of the onset condition unlike in the planar geometry equivalent. Results also show that disturbance growth seems to be strongly associated with the deformation direction instead of the acceleration direction. For example, when the interface is compressed in the cross-stream direction, the disturbance at the gas–solid interface grows even under the Rayleigh–Taylor stable situation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Nonlinear error field response in the presence of plasma rotation and real frequencies due to favorable curvature

This works describes the interaction of an phenomenon known to arise in tokamak fusion experiments with the magnetic perturbations that can arise at the edge of such experiments. We present nonlinear NIMROD resistive MHD simulations of the response of a rotating plasma to an error field when the plasma has weakly damped linear tearing modes (TMs), stabilized by a pressure gradient and favorable curvature. The favorable curvature leads to the Glasser effect: the occurrence of real frequencies and stabilization with positive stability index Δ'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs’ birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Direct measurements of the 3D plasma velocity in single-helical-axis RFP plasmas

The first local velocity measurements of helical equilibrium plasmas in the Reversed Field Pinch (RFP) Single Helical Axis (SHAx) state using a Charge Exchange Recombination Spectroscopy (CHERS) diagnostic are presented. Measurements show strong axisymmetric and non-axisymmetric flow, with n ¼ 5 components of flow related to the ðm; nÞ¼ð1; 5Þ dominant magnetic mode on the order of the axisymmetric flow in certain regions of the plasma, as well as significant n > 5 flow. Flow measurements are compared with NIMROD simulations of visco-resistive, single-fluid MHD in toroidal and cylindrical geometries with limited axial periodicity. Both measurements and the simulation with toroidal geometry show stronger inboard flows relative to the outboard flows, which is attributed to the toroidal geometry of the device. In the experiment, the n ¼ 5 component of flow is phase shifted from the reconnection-like flow pattern observed in the single-fluid simulations, possibly due to decoupling of the ion and electron fluids over much of the plasma. Finally, the strength of the helical angular flow shear relative to the critical shear necessary to disrupt nonlinear coupling between tearing modes is calculated around the helical magnetic axis. The shear in the measured flow is on the order of the theoretical critical threshold needed to nonlinearly decouple modes, but the measurement uncertainty in the gradient of the flow is large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Suppression of ITG turbulence due to spectral shift during biasing induced H-mode on HBT-EP

Investigations of biasing induced H-mode discharges on HBT-EP show that the edge turbulence is consistent with the ion temperature gradient) mode and have allowed for the controlled observation of the effect of applied flow shear on the turbulence. Measurements of the radial wavenumber spectrum of floating potentials at the edge show that the turbulence intensity decreases with increasing shift in the spectrum average ⟨ k r ⟩ when increasing amounts of bias probe voltage (and increasing amounts of flow shear) is applied. This is in agreement with the spectral shift model [Staebler et al., Phys. Rev. Lett. 110, 055003 2013] for turbulence suppression via sheared flow. A shift in the wavenumber spectrum occurs at applied electrode voltages and currents below the threshold needed for an L–H transition, and a dithering transition is obtained when biasing near the threshold. Suppression of blob-filament turbulence in the scrape-off layer (SOL) precedes the L–H transition, with the SOL turbulence remaining low throughout the dithering phase, despite the modulation of turbulence levels in the nearby edge. This demonstrates that the SOL turbulence “decouples” from the edge turbulence. The spectral shift in the measured radial wavenumber is corroborated by the direct measurement of eddy tilt angle using a novel time delay analysis technique first developed for Doppler reflectometry [Pinzon et al., Plasma Phys. Controlled Fusion 61, 105009 (2019)] but adapted here for floating potential measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toward the core-edge coupling of delta-f and total-f gyrokinetic models

The coupling of core reduced delta-f and edge total-f gyrokinetic models could enable a significant speed-up of the whole device modeling simulations. It is thus interesting to compare these models in the core where they will be coupled. The main difference between these models is the presence of the zeroth order term on the right-hand side of the total-f gyrokinetic equation. This term is responsible for driving a radial electric field that is also associated with Geodesic acoustic mode-like oscillations. To investigate the coupling between these models, the subsequent large-scale gyrokinetic 3D turbulence simulations have been performed starting from an already saturated gyrokinetic axisymmetric equilibria (collisions are not included). This work has also been extended to couple different axisymmetric and turbulent models available in XGC. It is shown that the radial electric field and its drive have to be taken care of consistently while coupling different models together.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Non-ideal MHD growth of current interchange tearing modes at plasma edge and response to externally-imposed flow

Current interchange tearing modes (CITMs) have been proposed to explain an intermittent eruption of electric current toward a divertor at an edge region of a tokamak device. Extended magnetohydrodynamic (MHD) simulations combined with the Scrape-Off-Layer diffusion model, which has been developed by the authors [Miura et al., Phys. Plasmas 24, 092111 (2017)], are carried out in order to study the growth of a CITM under influences of two-fluid and gyro-viscous effects, or non-ideal MHD effects. Numerical simulations both with and without an externally imposed azimuthal flow show that the growth of a CITM is observed for a relatively small flow velocity, and that a CITM can be a candidate mechanism for the intermittent eruption of the current in a tokamak, whether the growth is under influences of non-ideal MHD effects or not. Here, the growth can be suppressed when an externally imposed azimuthal flow is sufficiently large to cause a finite radial displacement because such a displacement prevents a transition from an interchange to a tearing mode that is the key process of a CITM. Furthermore, it is also found that a stripe pattern similar to so-called streamers is formed in the course of the CITM growth. The basic nature of this streamer-like structure is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measuring higher-order moments of neutron-time-of-flight data for cryogenic inertial confinement fusion implosions on OMEGA

Ion temperatures serve as an important diagnostic for inertial confinement fusion (ICF) implosions. In direct-drive ICF experiments on OMEGA, neutron-time-of-flight (nTOF) data are used to infer the ion temperature of the fusing plasma produced in the implosion experiment. The analysis of the nTOF data requires an assumption about the shape of the underlying source signal. Since the source nTOF signal is a near-replica of the neutron energy spectrum, an ideal Gaussian shape, corresponding to the neutron energy spectrum of a uniform temperature plasma, is routinely employed. However, spatial and temporal variations of the ion temperature in the plasma give rise to higher-order moments, which were first described by Munro [Nucl. Fusion 56, 036001 (2016)]. In this work, we show a simpler alternative analysis to derive moments of the neutron energy spectrum for a plasma with variations in ion temperature. We also present a revised analysis of measured nTOF signals that uses a model with an additional degree of freedom to take into account the effect of ion temperature variations on the shape of the spectrum. Compared to presently used nTOF analysis, the revised analysis yields on average ≈2× more accurate fits to the data and up to 15% higher ion temperatures for cryogenic experiments. As a result, we quantify the ion temperature inflation caused by radially symmetric fluid flows, which are present even in a symmetric implosion, and which serve as a lower bound on the ion temperature inflation in real implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

When do waves drive plasma flows?

Flows and rotation, particularly E×B rotation, are critical to improving plasma performance, and waves are a primary tool of plasma control. Thus, it is paramount to understand under what conditions waves can drive E×B flows in plasmas. In this didactic review, an invited paper accompanying the 2023 Marshall N. Rosenbluth Doctoral Thesis Award, this question is answered in the context of momentum-conserving quasilinear theory. There are two primary frameworks for momentum-conserving quasilinear theories that can handle both resonant and nonresonant particles: Eulerian averaging theories and oscillation-center Hamiltonian theories. There are also two different paradigmatic wave problems: plane-wave initial value problems, and steady-state boundary value problems. Here, it is shown that each of these frameworks “naturally” works better with a different problem type. By using these theories, one finds a great difference in the behavior of time- vs space-dependent waves. A time-evolving plane wave can only drive flow if the electromagnetic momentum of the wave, given by the Poynting flux, changes. This result precludes flow drive by any planar electrostatic wave. In contrast, a steady-state spatially evolving wave can drive flow whenever there is divergence in the flux of Minkowski momentum, a completely different physical quantity. This review aims to provide a high-level, intuitive understanding of the very different behaviors observed for these two types of problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear coupling of tearing modes in reversed field pinch plasmas with stepped pressure profiles

A theory of the three-wave coupling of triplets of tearing modes in toroidal pinches [i.e., either reversed field pinches (RFPs) or tokamaks] was proposed by R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)]. However, this theory only applies to toroidal pinches with negligible equilibrium plasma pressure gradients. Such a limitation is particularly inappropriate to RFPs. This paper generalizes the analysis of R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)] in order to take the equilibrium pressure gradient into account. However, for the sake of simplicity, a stepped pressure profile, rather than a continuous profile, is employed. In the limit in which the number of steps becomes very large, the results obtained from the generalized theory are presumably equivalent to those that would have been achieved using a continuous pressure profile. The generalized theory is used to investigate the formation of the characteristic toroidally localized pattern of phase-locked m = 1 and m = 0 tearing modes in RFP plasmas that is known as the “slinky” pattern. The incorporation of the equilibrium plasma pressure into the analysis is found to be of crucial importance when determining the properties of the pattern. This is the case because the plasma pressure controls the number of unstable m = 1 and m = 0 tearing modes, and also significantly affects the strength of three-wave coupling, as well as the phase relation between the phase-locked m = 1 and m = 0 modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Mode transitions and spoke structures in E × B Penning discharge

Two-dimensional particle-in-cell simulations in the (radial-azimuthal) plane perpendicular to the axial direction of a cylindrical $E$ x $B$ Penning discharge are presented. The low-pressure discharge is self-consistently supported by plasma ionization from the electron beam injected axially, along the direction of the external magnetic field. It is shown that with the increasing strength of the external magnetic field, the discharge undergoes a sequence of transitions between several azimuthal modes. Azimuthal m > 1 spiral arm structures are excited at low magnetic field values as plasma confinement improves and the radial density profile becomes peaked. With a larger field, spiral arms with m > 1 are replaced by the m = 1 spoke mode, most clearly seen in plasma density. A transition from spiral arms to the spoke regime occurs when the plasma potential in the center changes from weakly positive (or zero) to negative. Further increase in the magnetic field results in a well-developed m = 1 spoke mode with additional small-scale higher-frequency m > 1 structures inside and around the spoke. It is shown that while ionization and collisions affect some characteristics of the observed fluctuations, the basic features of the spoke and m > 1 spiral structure remained similar without ionization. The role of energy conservation in small-scale high-frequency modes and spoke dynamics is discussed. It is demonstrated that in regimes with the m = 1 spoke mode, additional m = 4 harmonics of the ion and electron fluxes to the wall appear due to the square boundary. The frequency of the m = 1 mode is weakly affected by the geometry of the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗