Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Interpenetrating Plasma Simulations

This project investigated multifluid models for simulating interpenetrating plasmas. These models solve the Euler equations for each plasma stream along with an appropriate model for electron dynamics. The project developed a novel code, Euclid, that simulates these models using high-order conservative spatial discretization methods and explicit and semi-implicit time integration methods. Euclid was used to solve several multifluid test cases representative of plasma dynamics in high-energy density physics experiments. Kinetic simulations of interpenetrating flows using the LOKI code revealed novel physical effects of the ion-ion streaming instability. A new multi-flow collision operator was developed for the LOKI code. A simple model of the ion-ion streaming instability was also developed. Several presentations were given on the project results; three journal articles have been published, one has been submitted, and two are planned.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2020

Time and space resolved diagnostics for plasma thermal-chemical instability of fuel oxidation in nanosecond plasma discharges

Here, an instability in a nanosecond pulsed dielectric barrier discharge plasma occurring in methane–oxygen–argon mixtures is experimentally observed and measured by 1D time-resolved in situ electric field measurements. This instability, which seems to be created by the positive feedback between plasma kinetics and plasma-assisted low temperature fuel oxidation, is studied using electric field induced second harmonic generation and direct ICCD imaging. The rapid formation of streamers from an originally uniform discharge appears to be caused by the chemical kinetics of plasma-assisted low temperature methane oxidation, which may be resulting in a new type of plasma instability: a thermal-chemical instability. The results also revealed that the occurrence of this possible thermal-chemical instability in a reactive flow drastically changes the plasma properties by forming multiple secondary discharges and possibly leads to micron-sized non-uniform electric distributions. Single shot uncalibrated measurements of the electric field of the micron sized streamers appears to show much greater strengths than the average electric field. Furthermore, one-dimensional data analysis shows the positive feedback loop between the streamers and the low temperature plasma assisted oxidation chemistry in the plasma thermal-chemical instability. The present finding advances the understanding plasma instability growth and provides a new way to control plasma uniformity in plasma-assisted combustion and plasma fuel reforming.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effect of triangularity on ion-temperature-gradient-driven turbulence

The linear and nonlinear properties of ion-temperature-gradient-driven turbulence with adiabatic electrons are modeled for axisymmetric configurations for a broad range of triangularities δ, both negative and positive. Peak linear growth rates decrease with negative δ but increase and shift toward a finite radial wavenumber kx with positive δ. The growth-rate spectrum broadens as a function of kx with negative δ and significantly narrows with positive δ. Here, the effect of triangularity on linear instability properties can be explained through its impact on magnetic polarization and curvature. Nonlinear heat flux is weakly dependent on triangularity for |δ|≤0.5, decreasing significantly with extreme δ, regardless of sign. Zonal modes play an important role in nonlinear saturation in the configurations studied, and artificially suppressing zonal modes increased nonlinear heat flux by a factor of about four for negative δ, increasing with positive δ by almost a factor of 20. Proxies for zonal-flow damping and drive suggest that zonal flows are enhanced with increasing positive δ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Structure formation in turbulence as an instability of effective quantum plasma

Structure formation in turbulence can be understood as an instability of “plasma” formed by fluctuations serving as effective particles. These “particles” are quantumlike in the sense that their wavelengths are non-negligible compared to the sizes of background coherent structures. The corresponding “kinetic equation” describes the Wigner matrix of the turbulent field, and the coherent structures serve as collective fields. This formalism is usually applied to manifestly quantumlike or scalar waves. Here, we show how to systematically extend it to more complex systems using compressible Navier–Stokes turbulence as an example. In this case, the fluctuation Hamiltonian is a five-dimensional matrix operator and diverse modulational modes are present. As an illustration, we calculate these modes for a sinusoidal shear flow and find two modulational instabilities. One of them is specific to supersonic flows, and the other one is a Kelvin–Helmholtz-type instability that is a generalization of the known zonostrophic instability. Our calculations are readily extendable to other types of turbulence, for example, magnetohydrodynamic turbulence in plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗
DOE OSTI2023

Fast particles in drift wave turbulence

This study aims to incorporate the effects of fast particles into our present fluid model for tokamak transport. The parameter ε f = ω / ω f, where ω is the mode frequency and ω f is the typical frequency of the fast particles, which enters as a factor in front of the fast particle response. Thus, for trapped fast particles, where ω f = ω pres the precession frequency of the fast particles, this parameter is of order 10 – 2 for drift waves, and thus, the fast particle response can be neglected. However, ε f will be of order 1 for fast particle modes such as in the fishbone instability. An important turbulence property, affecting both these limits, is resonance broadening. Effects of resonance broadening have recently been considered for fast particle instabilities, often coupled directly to the linear growth rate, while we here consider the original Dupree formulation where the turbulence directly drives a nonlinear frequency shift. Resonance broadening has a general tendency to counteract dissipative wave particle resonances. This has been observed for fast particle instabilities. Here, there is a resonant external source for the fast particles, so the instability survives if this source is dominant over the resonance broadening. For drift waves, however, external sources are not resonant since ε f << 1. Furthermore, the resonance broadening is able to remove the dissipative wave particle resonance completely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗
DOE OSTI2023

Large Eddy Simulation of Rotationally Induced Ingress and Egress around an Axial Seal between Rotor and Stator Disks

In gas turbines, the hot gas exiting the combustor can have temperatures as high as 2000 °C, and some of this hot gas enter into the space between the stator and rotor disks (wheelspace). Since the entering hot gas could damage the disks, its ingestion must be minimized. This is carried out by rim seals and by introducing a cooler flow from the compressor (sealing flow) into the wheelspace. Ingress and egress into rim seals are driven by the stator vanes, the rotor and its rotation, and the rotor blades. This study focuses on the ingress and egress driven by the rotor and its rotation. This is carried out by performing wall-resolved large eddy simulation (LES) around an axial seal in a rotor–stator configuration without vanes and blades. Results obtained show the mechanisms by which the rotor and its rotation induce ingress, egress, and flow trajectories. Kelvin–Helmholtz instability was found to create a wavy shear layer and displacement thickness that produces alternating regions of high and low pressures around the rotor side of the seal. Vortex shedding on the backward-facing side of the seal and its impingement on the rotor side of the seal also produces alternating regions of high and low pressures. The locations of the alternating regions of high and low pressures were found to be statistically stationary and to cause ingress to start on the rotor side of the seal. Vortex shedding and recirculating flow in the seal clearance also cause ingress by entrainment. With the effects of the rotor and its rotation on ingress and egress isolated, this study enables the effects of stator vanes and rotor blades to be assessed.

42 ENGINEERING↗
DOE OSTI2025

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Flow Z-pinch plasma production on the FuZE experiment

The Fusion Z Pinch Experiment (FuZE) investigates sheared-flow stabilization of classic m = 0 and m = 1 instabilities in Z pinches with an embedded axial flow. FuZE consists of a 100 cm coaxial plasma accelerator, where neutral gas is ionized and accelerated in a pulsed electrical discharge, followed by a 50 cm assembly region, where pinches are formed. Maintaining the pinch requires continuous plasma injection provided by a deflagration mode in the coaxial accelerator. Two discharge modes, with and without deflagration, are investigated on FuZE. Pinch formation is observed with deflagration only. Plasma velocities in the assembly region are found to match the E × B velocity estimated in the accelerator based on a 1D circuit model, indicating that a 1D MHD approximation may offer a valid description of the plasma in the accelerator channel. The velocity of magnetic field propagation is found to agree with the snowplow model based on momentum conservation, and the lifetime of the pinch is shown to be in agreement with constraints imposed by mass conservation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nonlinear harmonics coupled by parallel wave propagations in a time-dependent plasma flow

In a time-dependent flow, nonlinear harmonics can be excited by coupling between linear waves and flow-induced harmonic waves. Examining the dispersion relations and selection rules for the coupling, we investigate nonlinearly coupled harmonics for waves propagating along the magnetic field line in a magnetized plasma, as well as waves in an unmagnetized plasma. The coupled harmonics in a plasma flow are described by analytic dispersion relations and selection rules. This nonlinear coupling is corroborated by the particle-in-cell simulation. Here, the coupled-harmonics model describes a mechanism for the excitation of nonlinear harmonics from linear waves in a time-dependent flow. The spectral analysis of the dispersion relation provides a useful way to evaluate the spatiotemporal behavior of a plasma flow.

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