Electrohydrodynamic Rayleigh-Taylor bulk instability
Stability criterion for initially static, stratified fluid subject to electrical stress
Melcher, J. R.↗
SEARCH · Search NASA
Results for “TAYLOR 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.
Bound oscillations on thin magnetic flux tubes - Convective instability and umbral oscillations
The possibility that 'tube waves' can be trapped on slender solar magnetic flux tubes is investigated. For rigid isothermal flux tubes, it is found that the flux tube geometry can by itself lead to waves which are trapped on the part of the tube that expands with height. Some geometries lead to trapped modes with eigenperiods near 180 s, if parameters appropriate to sunspot umbrae are chosen. It is possible that the umbral oscillations are a manifestation of such trapped waves, if sunspot umbrae consist of an assembly of slender flux tubes, as in the spaghetti model of Parker (1979). For flux tubes which have a constant ratio of Alfven speed to sound speed, it is found that it is primarily the variation of temperature with height which determines whether trapped waves can exist. Certain temperature profiles lead to disturbances for which omega squared is less than zero, corresponding to convective instability or Rayleigh-Taylor instability.
Hollweg, J. V.↗
The nonlinear breakup of a magnetic layer - Instability to interchange modes
Motivated by considerations of the solar toroidal magnetic field, the behavior of a layer of uniform magnetic field embedded in a convectively stable atmosphere is studied. Since the field can support extra mass, such a configuration is top-heavy and thus instabilities of the Rayleigh-Taylor type can occur. For both static and rotating basic states, the evolution of the interchange modes (no bending of the field lines) is followed by integrating numerically the nonlinear compressible MHD equations. The initial Rayleigh-Taylor instability of the magnetic field gives rise to strong shearing motions, thereby exciting secondary Kelvin-Helmholtz instabilities which wrap the gas into regions of intense vorticity. The subsequent motions are determined primarily by the strong interactions between vortices which are responsible for the rapid disruption of the magnetic layer.
Cattaneo, F.↗
Competing mechanisms at vibrated interfaces of density-contrast fluids
Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.
Direct numerical simulations↗
An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa
Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.
97 MATHEMATICS AND COMPUTING↗
Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures
Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.
Environmental Sciences & Ecology↗
Stabilization of Liner Implosions via a Dynamic Screw Pinch
Magnetically driven implosions are susceptible to magnetohydrodynamic instabilities, including the magneto-Rayleigh-Taylor instability (MRTI). To reduce MRTI growth in solid-metal liner implo-sions, the use of a dynamic screw pinch (DSP) has been proposed [P. F. Schmit et al., Phys. Rev. Lett. 117, 205001 (2016)]. In a DSP con guration, a helical return-current structure surrounds the liner, resulting in a helical magnetic eld that drives the implosion. Here, we present the rst experimental tests of a solid-metal liner implosion driven by a DSP. Using the 1-MA, 100{200-ns COBRA pulsed-power driver, we tested three DSP cases (with peak axial magnetic elds of 2 T, 14 T, and 20 T) and a standard z-pinch (SZP) case (with a straight return-current structure and thus zero axial eld). The liners had an initial radius of 3.2 mm and were made from 650-nm-thick aluminum foil. Images collected during the experiments reveal that helical MRTI modes developed in the DSP cases, while non-helical (azimuthally symmetric) MRTI modes developed in the SZP case. Additionally, the MRTI amplitudes for the 14-T and 20-T DSP cases were smaller than in the SZP case. Speci cally, when the liner had imploded to half of its initial radius, the MRTI amplitudes for the SZP case and for the 14-T and 20-T DSP cases were, respectively, 1.1-0.3 mm, 0.7-0.2 mm, and 0.3-0.1 mm. Relative to the SZP, the stabilization obtained using the DSP agrees reasonably well with theoretical estimates.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Ionospheric Instabilities and Shock Capturing Algorithms
We describe the basic physics of the gradient-drift instability and the collisional Rayleigh-Taylor instability, the two plasma instabilities which dominate the production of small-scale structure in the earth's ionosphere. We then focus our attention on the nonlinear evolution of these instabilities for two particular cases: 1) the recursive bifurcation and "freezing" of artificial ion clouds released in the ionosphere; and 2) equatorial spread F. Finally, we discuss the numerical challenges faced in attempting to numerically simulate these instabilities, and the role played by shock-capturing algorithms therein.
Zalesak, Steven T.↗
Radio emission from young supernova remnants - Effects of an inhomogeneous circumstellar medium
The evolution of young supernova remnants has been modeled using a one-dimensional hydrodynamics code. Rayleigh-Taylor instabilities and the growth of Rayleigh-Taylor fingers have been included in the code. Turbulent dynamo amplification of magnetic fields and both turbulent and shock acceleration of relativistic electrons have been included macroscopically. From this, the distribution of synchrotron luminosity in the remnant has been calculated. It is found that the radio morphology of model remnants expanding into a homogeneous medium does not agree with observations. Expansion into a circumstellar medium with many small cloudlets does produce radio shells which agree with observations. It is suggested that supernova remnants reflect the interaction of ejected matter with a cloudy circumstellar medium.
Dickel, John R.↗
Material effects on dynamics in triple-nozzle gas-puff Z pinches
Here, the gas-puff Z pinch has a long history with myriad applications as an efficient neutron or x-ray source. Its simplicity as a load configuration makes it suitable for studying fundamental plasma physics phenomena such as instabilities and energy transport. For example, the implosion of cylindrical shells onto a fusion fuel are inherently susceptible to instability growth on their external surfaces; if such instabilities are unmitigated, then the consequences in terms of degraded performance can be substantial. Similarly, mitigating heat transport from a hot fuel to its colder surrounding container can make fusion conditions more easily achievable. Here we have conducted a systematic study of triple-nozzle (outer liner, inner liner, fuel) gas puffs using two-dimensional (2D) magnetohydrodynamic simulations to investigate the effect of load material on the relevant dynamics. Analogous to past studies on spherical blast waves and converging shock waves, a trend emerges linking increased radiative cooling, lower adiabatic index, and increased magneto-Rayleigh-Taylor instability growth. Notably, our results suggest that, for the present configuration, Ar radiates less than both Ne and Kr during the early stages of the implosion while mass is being swept up and perturbations begin to seed instability growth. Consequently, pinches with Ar on the outer surface exhibit more stable 2D behavior. Here we also present a parameter scan of thermonuclear neutron yield, Y, as a function of peak current, I pk and dopant concentration with Ne or Ar, depending on the inner liner material. Above 6 MA, our results suggest Y ∝ $I$$^{5}_{pk}$ and even substantial mixing (10% by volume) of Ne into the fuel does not drastically reduce yield, suggesting an Ar/Ne/fuel configuration may reliably achieve DD thermonuclear yields of 10 13 – 10 14 /cm in the 10–20 MA range.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Techniques for improved statistical convergence in quantification of eddy diffusivity moments
While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪(100) simulations rather than 𝒪(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.
general physics↗
Mass and momentum transport in the Tilted Rocket Rig experiment
Here, high-fidelity large eddy simulations of the inclined Rayleigh-Taylor instability in the Tilted Rocket Rig experimental configuration [Smeeton and Youngs, AWE Report No. O 35/87 (1987)] are performed using a tenth-order compact finite difference code. These simulations are analyzed for spatial distributions of turbulent kinetic energy, turbulent mass flux velocity, species mass fraction flux, species mass fraction variance, and Favre-averaged Reynolds stresses at two time instances, t = 37 ms and t = 52 ms. Additionally, the vertical distribution of the components of the unclosed budget equations over the center of the domain for these quantities are examined. The dominant terms of these budget equations are further decomposed to examine the principal contributions to these terms in each axis. Notably, the principal contribution to the horizontal turbulent mass flux velocity budget is found to be from a term which is commonly neglected in many Reynolds-averaged Navier-Stokes models.
42 ENGINEERING↗
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↗
Stereo Imaging Velocimetry of Mixing Driven by Buoyancy Induced Flow Fields
Mixing of two fluids generated by steady and particularly g-jitter acceleration is fundamental towards the understanding of transport phenomena in a microgravity environment. We propose to carry out flight and ground-based experiments to quantify flow fields due to g-jitter type of accelerations using Stereo Imaging Velocimetry (SIV), and measure the concentration field using laser fluorescence. The understanding of the effects of g-jitter on transport phenomena is of great practical interest to the microgravity community and impacts the design of experiments for the Space Shuttle as well as the International Space Station. The aim of our proposed research is to provide quantitative data to the community on the effects of g-jitter on flow fields due to mixing induced by buoyancy forces. The fundamental phenomenon of mixing occurs in a broad range of materials processing encompassing the growth of opto-electronic materials and semiconductors, (by directional freezing and physical vapor transport), to solution and protein crystal growth. In materials processing of these systems, crystal homogeneity, which is affected by the solutal field distribution, is one of the major issues. The understanding of fluid mixing driven by buoyancy forces, besides its importance as a topic in fundamental science, can contribute towards the understanding of how solutal fields behave under various body forces. The body forces of interest are steady acceleration and g-jitter acceleration as in a Space Shuttle environment or the International Space Station. Since control of the body force is important, the flight experiment will be carried out on a tunable microgravity vibration isolation mount, which will permit us to precisely input the desired forcing function to simulate a range of body forces. To that end, we propose to design a flight experiment that can only be carried out under microgravity conditions to fully exploit the effects of various body forces on fluid mixing. Recent flight experiments, by the P.I. through collaboration with the Canadian Space Agency (STS-85, August 1997), aimed at determining the stability of the interface between two miscible liquids inside an enclosure show that a long liquid column (5 cm) under microgravity isolation conditions can be stable, i.e. the interface remains sharp and vertical over a short time scale; thus transport occurs by molecular mass diffusion. On the other hand, when the two liquids were excited from a controlled vibration source (Microgravity Vibration Isolation Mount) two to four mode large amplitude quasi-stationary waves were observed. The data was limited to CCD recording of the dynamics of the interface between the two fluids. We propose to carry out flight experiments to quantify the dynamics of the flow field using Stereo Imaging Velocimetry and measure the concentration field using laser fluorescence. The results will serve as a basis to understand effects of g-jitter on transport phenomena, in this case mass diffusion. As the measurement of the kinematics of the flow field will shed light on the instability mechanism. The research will allow measurement of the flow field in microgravity environment to prove two hypotheses: (1) Maxwell's hypothesis: finite convection always exists in diffusing systems, and (2) Quasi-stationary waves inside a bounded enclosure in a microgravity environment is generated by Kelvin-Helmholtz instability; resonance of the interface which produces incipient mixing is due to Rayleigh-Taylor instability. The first hypothesis can be used as a benchmark experiment to illustrate diffusive mixing. The second hypothesis will lead to the understanding of g-jitter effects on buoyancy driven flow fields which occur in many situations involving materials processing, and other basic fluid physics phenomena. In addition, the second hypothesis will also provide insight in how Rayleigh-Taylor and Kelvin-Helmholtz instabilities propagate concentration fronts during mixing. Measurement of the flow field using SIV is important because it is the flow field which causes instability at the interface between the two fluids. Mixing driven by buoyancy induced flow fields will be addressed both experimentally and computationally. The experimental effort will address the kinematics of mixing: stretching, transport and chaos. Quantification of the mechanisms of mixing will consists of measuring the flow field using the SIV system at Glenn and capturing the dynamics of the interface, to measure mass transport, using a CCD camera. These experiments will be carried out within the framework of Earth's gravity and g-jitter microgravity acceleration as in a Space Shuttle environment or the International Space Station. The g-jitter will be induced and controlled using a tunable vibration isolation platform to isolate against vibration as well as input periodic and random vibration to the system. The parametric range of the microgravity experiment will be extended from the experiments on STS-85 to investigate higher mode quasi-stationary waves (8 to 12), as well as resonance regions which leads to chaos and turbulence. Ground-based experiments will focus on effects of vibration on stably stratified fluid layers in order to scale for possible scenarios in a microgravity environment. These vibrations will be subjected perpendicular to the concentration field on the ground since the parallel case can only be carried out in a microgravity environment. The concept of dynamical similarity will be applied to tune the experiments as closely as possible to a Space Shuttle environment or the International Space Station. The computational effort will take advantage of the Computational Laboratory at Glenn to corroborate the experimental findings with predictions of the dynamics of the flow field using the codes FLUENT (finite difference based) and FIDAP (finite element based). We will investigate two important cases, single-fluid model to address dilute systems with negligible jump in viscosity and the more general two-fluid model which accounts for finite jump in viscosity. Apart from its microgravity relevance, this experiment is well suited to study dynamics in nonlinear systems.
Duval, W. M. B.↗
Kinetic Simulations of the Kruskal–Schwarzschild Instability in Accelerating Striped Outflows: Dynamics and Energy Dissipation
Astrophysical relativistic outflows are launched as Poynting-flux dominated, yet the mechanism governing efficient magnetic dissipation, which powers the observed emission, is still poorly understood. We study magnetic energy dissipation in relativistic "striped" jets, which host current sheets separating magnetically dominated regions with opposite field polarity. The effective gravity force g in the rest frame of accelerating jets drives the Kruskal–Schwarzschild instability (KSI), a magnetic analog of the Rayleigh–Taylor instability. By means of 2D and 3D particle-in-cell simulations, we study the linear and nonlinear evolution of the KSI. The linear stage is well described by linear stability analysis. The nonlinear stages of the KSI generate thin (skin-depth-thick) current layers, with length comparable to the dominant KSI wavelength. There, the relativistic drift-kink mode and the tearing mode drive efficient magnetic dissipation. The dissipation rate can be cast as an increase in the effective width Δ eff of the dissipative region, which follows dΔ eff /dt ≃ 0.05 $\sqrt{Δ_{eff} g}$. Our results have important implications for the location of the dissipation region in gamma-ray burst and active galactic nuclei jets.
79 ASTRONOMY AND ASTROPHYSICS↗
Final Report for Tech-X Corporations’s contribution to Enhancing Understanding of High Energy Density Plasmas from Wire Array and Solid Liner Implosions Using Fluid Modeling with Kinetic Closures
Recent results from experiments and simulation of magnetically driven pulsed power liners have explored the role of the early-time electrothermal instability in the evolution of the magneto-Rayleigh-Taylor instability. Our focus will be on understanding the development of such instabilities and the potential stabilization mechanisms via electron thermal conduction, viscosity, and large magnetic fields which we expect could play a significant role in supporting the success of the MagLIF program. Experiments have shown that studies of high-energy density plasmas from wire-array implosions require physics modeling that goes well beyond simple models such as ideal magnetohydrodynamics. The goal of this work is to provide increased understanding of these experiments by employing simulations with a multi-fluid extended-MHD model which uses kinetic closures for thermal conduction, resistivity and viscosity. We will use codes easily available to the wider research community, including university students, with a secondary goal of providing the community with well-benchmarked tools capable of advanced modeling of high-energy-density plasmas.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Stability of an oscillated fluid with a uniform density gradient
Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.
Jacqmin, David↗
Collisional interchange instability. I - Numerical simulations of intermediate-scale irregularities
Numerical simulations of the generalized Rayleigh-Taylor instability are presented. The model and simulations are applicable to bottomside and topside spread F, unstable barium cloud dynamics, and collisional interchange instability in general. The principal result is that the evolution of the instability tends to an anisotropic state consisting of nearly sinusoidal (quasi-periodic) variation along the effective electric field, and shocklike structures propagating perpendicular to E(eff) along the extrema of the quasi-periodic structures. The spectral properties of the nonlinear state are analyzed using one-dimensional power spectra calculated along spatial trajectories for selected angles to E(eff). In this way a direct comparison to in situ probe data can be made. The inherent anisotropy of the nonlinear state is reflected in major qualitative differences between the spectra taken parallel to and perpendicular to E(eff). The fundamental finding of the present work is that anisotropy in interchange dynamics is much greater than had been previously reported. This strong anisotropy can explain much of the spectral and spatial structural characteristics of both bottomside and topside spread F. In a companion paper a comparison of the simulation results to various in situ data sets is given.
Zargham, S.↗