Search NASA⌕ Search

SEARCH · Search NASA

Results for “Rayleigh-Taylor 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 163 records · Page 9

Granular instability in fluidized beds: A small-over-large instability

For different particle diameters arranged in a gas fluidized bed, a new gravitational instability is predicted. The first order stability analysis yields conditions for growth of the most unstable wavenumber both as a function of particle diameter and dimensionless Bond number. Gravity acts to distribute the greater fluid mass associated with smaller particles, while surface tension stabilizes the entire bed against accelerations. The phenomenon, called the 'Small-over-Large' (SL) instability, evolves somewhat akin to a Rayleigh-Taylor mechanism and in addition to fundamental interest, carries implications for predicting solute redistribution and gravitational alignments in particle beds and chemical separation columns.

Noever, David A.↗

Interstellar clouds in high-speed, supersonic flows: Two-dimensional simulations

We present a series of gasdynamical simulations of the interaction of a dense, cool interstellar cloud with a high-speed, supersonic wind that confines and accelerates the embedded cloud. Our goal is to attempt to determine if such clouds can survive various potentially disruptive instabilities, that occur at their peripheries, long enough to be accelerated to speeds which are comparable to the wind velocity. These simulations are performed using two-dimensional, Eulerian gas dynamics on both an axisymmetric (about the cloud axis) and 'slab' geometric grid. The spatial and temporal resolutions of the simulations are varied over a wide range to investigate the effects of small-scale instabilities on the overall acceleration of clouds and the development of large-scale, disruptive instabilities. Also, we study the effects of wind/cloud Mach number variations by changing the wind speed constant at about 12 km/s (which corresponds to a cloud temperature of 10,000 K). The current simulations track the evolution of clouds as they are accelerated to speeds approximately 4-5 times greater than their internal sound speeds. Furthermore, the models with the highest resolution were extended far beyond quasi-linear Rayleigh-Taylor growth times reaching 6-7 Rayleigh-Taylor growth times for the largest scale instabilities before being terminated because of the accumulation of errors at the rear grid boundary.

Schiano, A. V. R.↗

Observations and 3D MHD Modeling of Small-scale Solar Prominence Oscillations

Small scaleoscillations were recently observed in prominences with high resolution Hinode/SOT and IRIS SJI in unprecedented detailand identified as fast magnetosonic waves. The high-resolution observations provide evidence of flows of cool material, fluid instabilities such as Kelvin-Helmholtz and Rayleigh-Taylor, and the spectroscopic information provides evidence of Doppler shifts, heating, and cooling of the prominence material. Recently nonlinear fast magnetosonic waves were observed withperiods in the range of 5–11 minutes and wavelengths in the plane of the sky (POS) of ~2000kmand flows, and the flows in narrow threads withPOS speed ∼16–46 km/s. Thenonlinear fast magnetosonic waves and flows were modeled using 2.5D MHD (Ofman & Kucera 2020). Here, we extend the studyof the nonlinear fast magnetosonic waves using more realistic 3D MHD model. The model includes the 3D structure of the density and the magnetic fieldas well as mode coupling anddemonstrate the propagating density compressions associated with the nonlinear fast magnetosonicwaves, as well as the guided modesin the pillar. The 3DMHD modelingimprove understanding ofthe POSeffects in limb prominence oscillations observations.The results of the modeling are useful for coronal seismology of the prominence structure.

Leon Ofman↗

On the study of hydrodynamic instabilities in the presence of background magnetic fields in high-energy-density plasmas

Blast-wave-driven hydrodynamic instabilities are studied in the presence of a background B-field through experiments and simulations in the high-energy-density (HED) physics regime. In experiments conducted at the Laboratoire pour l’utilisation des lasers intenses (LULI), a laser-driven shock-tube platform was used to generate a hydrodynamically unstable interface with a prescribed sinusoidal surface perturbation, and short-pulse x-ray radiography was used to characterize the instability growth with and without a 10-T B-field. The LULI experiments were modeled in FLASH using resistive and ideal magnetohydrodynamics (MHD), and comparing the experiments and simulations suggests that the Spitzer model implemented in FLASH is necessary and sufficient for modeling these planar systems. These results suggest insufficient amplification of the seed B-field, due to resistive diffusion, to alter the hydrodynamic behavior. Although the ideal-MHD simulations did not represent the experiments accurately, they suggest that similar HED systems with dynamic plasma-β (=2μ0ρv2/B2) values of less than ∼100 can reduce the growth of blast-wave-driven Rayleigh–Taylor instabilities. These findings validate the resistive-MHD FLASH modeling that is being used to design future experiments for studying B-field effects in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deceleration-stage Rayleigh–Taylor growth in a background magnetic field studied in cylindrical and Cartesian geometries

Experiments have identified the Rayleigh–Taylor (RT) instability as one of the greatest obstacles to achieving inertial confinement fusion. Consequently, mitigation strategies to reduce RT growth and fuel–ablator mixing in the hotspot during the deceleration phase of the implosion are of great interest. In this work, the effect of seed magnetic fields on deceleration-phase RT growth are studied in planar and cylindrical geometries under conditions relevant to the National Ignition Facility (NIF) and Omega experiments. The magnetohydrodynamic (MHD) and resistive-MHD capabilities of the FLASH code are used to model imploding cylinders and planar blast-wave-driven targets. Realistic target and laser parameters are presented that suggest the occurrence of morphological differences in late-time RT evolution in the cylindrical NIF case and a measurable difference in spike height of single-mode growth in the planar NIF case. The results of this study indicate the need for target designs to utilize an RT-unstable foam–foam interface in order to achieve sufficient magnetic field amplification to alter RT evolution. Benchmarked FLASH simulations are used to study these magnetic field effects in both resistive and ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rayleigh-Taylor modes in constant-density incompressible fluids accelerated by radiation pressure

The paper examines the behavior of linear perturbations in an incompressible fluid undergoing acceleration by radiation pressure, with reference to processes occurring in quasars, supernovae, and planetary nebulae. It is shown that, contrary to prior expectation, fluids accelerated by radiation pressure, are not always unstable to Rayleigh-Taylor modes. Some are, in fact, unstable, but the nature of the instability is qualitatively different.

Krolik, J. H.↗

Stability of fireballs and gamma-ray bursts

Fireballs are an essential part of any cosmological gamma-ray burst. We derive a stability criterion for fireballs and show that fireballs are Rayleigh-Taylor unstable in any region in which the entropy decreases outward. The instability begins to operate when the fireball becomes matter dominated. Among the possible implications of the instability are the following: (1) Conversion of a fraction of the radiation energy to a convective energy expressed in the motion of bubbles relative to each other. (2) Penetration of fast bubbles through slower ones and creation of high-gamma regimes which are essential for efficient conversion of the energy to gamma-rays. (3) Formation of rapid time variation (of the scale of the bubbles) in the observed gamma-rays.

Waxman, Eli↗

Post-Shot Report for OMEGA Double Cylinders (CylDRT 22B)

The direct-drive double cylinder experimental platform is a high-energy-density (HED) science platform designed to image an imploding cylindrical target. The target consists of a directly-driven outer cylinder and a shock-and-collision driven inner cylinder. The purpose of this platform is to study hydrodynamic instability growth on the inner cylinder, the outer surface of which is classically Rayleigh-Taylor unstable during the acceleration phase. We present results from recent experiments at the OMEGA laser facility. These experiments were designed as a proof-of-principle for the platform, using the same cylinder exterior dimensions and direct-drive beam configuration as previous single cylinder experiments. In these experiments, three sets of targets were fielded: no machined perturbations (smooth), a sinusoidal mode-10 perturbation on the outer surface of the inner cylinder, and a mode-20 perturbation on the outer surface of the inner cylinder. The primary diagnostic was a gated x-ray framing camera which imaged the backlit inner cylinder on-axis for sixteen frames over a time window of 1 ns for each shot. A second side-lit radiograph captured one image per shot, diagnosing axial uniformity. In this report we include an overview of the results from both the backlighter and the sidelighter diagnostics. We discuss at length the experimental analysis process. Finally, we present the results of the smooth target implosion trajectory and compare them to post-shot simulations. We see favorable agreement between simulation and experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The conductive propagation of nuclear flames. I - Degenerate C + O and O + Ne + Mg white dwarfs

The paper determines the physical properties - speed, width, and density structure - of conductive burning fronts in degenerate carbon-oxygen (C + O) and oxygen-neon-magnesium (O + Ne + Mg) compositions for a grid of initial densities and compositions. The dependence of the physical properties of the flame on the assumed values of nuclear reaction rates, the nuclear reaction network employed, the thermal conductivity, and the choice of coordinate system are investigated. The occurrence of accretion-induced collapse of a white dwarf is found to be critically dependent on the velocity of the nuclear conductive burning front and the growth rate of hydrodynamic instabilities. Treating the expanding area of the turbulent burning region as a fractal whose tile size is identical to the minimum unstable Rayleigh-Taylor wavelength, it is found, for all reasonable values of the fractal dimension, that for initial C + O or O + Ne + Mg densities above about 9 x 10 exp 9 g/cu cm the white dwarf should collapse to a neutron star.

Timmes, F. X.↗

Integrity of the Plasma Magnetic Nozzle

This report examines the physics governing certain aspects of plasma propellant flow through a magnetic nozzle, specifically the integrity of the interface between the plasma and the nozzle s magnetic field. The injection of 100s of eV plasma into a magnetic flux nozzle that converts thermal energy into directed thrust is fundamental to enabling 10 000s of seconds specific impulse and 10s of kW/kg specific power piloted interplanetary propulsion. An expression for the initial thickness of the interface is derived and found to be approx.10(exp -2) m. An algorithm is reviewed and applied to compare classical resistivity to gradient-driven microturbulent (anomalous) resistivity, in terms of the spatial rate and time integral of resistive interface broadening, which can then be related to the geometry of the nozzle. An algorithm characterizing plasma temperature, density, and velocity dependencies is derived and found to be comparable to classical resistivity at local plasma temperatures of approx. 200 eV. Macroscopic flute-mode instabilities in regions of "adverse magnetic curvature" are discussed; a growth rate formula is derived and found to be one to two e-foldings of the most unstable Rayleigh-Taylor (RT) mode. After establishing the necessity of incorporating the Hall effect into Ohm s law (allowing full Hall current to flow and concomitant plasma rotation), a critical nozzle length expression is derived in which the interface thickness is limited to about 1 ion gyroradius.

Gerwin, Richard A.↗

The Impact of Anomalous Resistivity in Vacuum Contaminant Plasmas on the Electrothermal Instability

This manuscript presents an assessment of the electrothermal instability (ETI) in the presence of anomalous resistivity (AR) in vacuum contaminant plasmas (VCP) when applied to a magnetized liner inertial fusion (MagLIF)-like load.Pulsed-power driven dielectrically coated metallic liners, like in MagLIF, experience the current-driven electrothermal instability which occurs when a material’s resistivity changes with temperature and is subject to ohmic heating. Large scale pulsed-power facilities that use magnetically insulated transmission lines (MITL) have been shown to generate low-density plasma which enters the target chamber and coalesces around the load. The low-density high-temperature vacuum contaminant plasmas (VCP) can parasitically divert current from the load through causing a short in the anode-cathode gap inside the target chamber. Resistive magnetohydrodynamic (MHD) simulations of these VCP experience unphysical runaway ohmic heating due to under predicting the resistivity by using a purely collisional resistivity model.AR provides a physics-based way to address this runaway heating through increasing the resistivity in a proportional way with the drift speed. In this work, 1D simulations probe the effect that AR in VCP has on the magnetic diffusion rate, and 2D simulations show how this effect manifests in the nonlinear striation form of the ETI for a MagLIF-like load. Beryllium and aluminum dielectrically coated liners are used for the 1D and 2D simulations in this work. The1D simulations show that a VCP causes a delay in the current delivery to the load by upwards of 8 ns at 60 ns into atypical current pulse for an 1 × 10 7 A scale pulsed-power accelerator. The 2D simulations show the delay observed in the 1D simulations is reduced substantially (4 ns) by 2D turbulence that disrupts the VCP layer early in time. Here, the 2Dsimulated ETI growth varies across AR models, more so for the beryllium liner than the aluminum liner because the beryllium liner shows an enhanced rate of penetration for the magnetic diffusion wave in comparison to the aluminum liner. The 2D simulations show the bulk dielectric thickness varies across AR models with the Davidson AR model beingthe largest and the Buneman AR model being the smallest, and in connection with the thickness the Rayleigh-Taylor bubble-spike distances varies correspondingly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Convective instabilities in SN 1987A

Following Bandiera (1984), it is shown that the relevant criterion to determine the stability of a blast wave, propagating through the layers of a massive star in a supernova explosion, is the Schwarzschild (or Ledoux) criterion rather than the Rayleigh-Taylor criterion. Both criteria coincide only in the incompressible limit. Results of a linear stability analysis are presented for a one-dimensional (spherical) explosion in a realistic model for the progenitor of SN 1987A. When applying the Schwarzschild criterion, unstable regions get extended considerably. Convection is found to develop behind the shock, with a characteristic growth rate corresponding to a time scale much smaller than the shock traversal time. This ensures that efficient mixing will take place. Since the entire ejected mass is found to be convectively unstable, Ni can be transported outward, even into the hydrogen envelope, while hydrogen can be mixed deep into the helium core.

Benz, Willy↗

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Star-forming instabilities of a decelerating plane-parallel slab of finite thickness

In order to study how supernova blast waves might catalyze star formation, the stability of a slab of decelerating gas of finite thickness is explored. Previous analyses are extended by applying shock-like boundary conditions to the leading edge of the slab and by studying the effects of arbitrary deceleration. Contrary to some earlier claims, it is found that blast waves can indeed accelerate the rate of star formation in the interstellar medium by compressing the interstellar gas and enabling it to cool. Also, it is demonstrated that in an incompressible fluid, the symmetric and antisymmetric modes in the case of zero acceleration transform continuously into Rayleigh-Taylor and gravity-wave modes as acceleration grows more important. The dynamical overstability of an isothermal shock wave and its implications for star formation are discussed.

Voit, G. Mark↗

Magnetic-Field Effects in Unstable High-Energy-Density Plasmas

The influence of magnetic fields in high-energy-density (HED) systems is a widely unexplored field of research that can play an important role in magnetized inertial confinement fusion (ICF) and in astrophysical systems. Hydrodynamic instabilities are presently believed to be the dominant degradation mechanism in ICF implosions. Therefore, understanding the influence of an imposed B-field on RT growth is important for future mitigation strategies as well as fundamental physics of ICF implosions that utilize B-fields for enhanced performance. The Crab Nebula, one of the most observed objects in our universe, contains elongated spikes of material created by hydrodynamic instabilities, however these spikes do not break-up and become turbulent as they evolve late in time. One hypothesis for this behavior is that B-fields around these spikes prevent vortex generation typical in these types of systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mixing and Transport in the Solar Nebula

Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.

Boss, Alan P.↗

Amplification of laser imprint in the presence of strong, externally imposed, target-normal magnetic fields

Here, an experiment was performed on the OMEGA EP Laser System to investigate the effects of magnetic fields on laser imprint. A 30 µ⁢m thick CH target was driven by a single beam delivering 3.0 kJ of UV energy in a 5 ns square pulse without smoothing by spectral dispersion. The Rayleigh-Taylor amplification of this beam's imprint on the target surface was monitored using face-on gated x-ray radiography from a Gd backlighter. Magnetic fields of up to 45 T were applied normal to the target surface. Analysis of the resulting radiographs shows a 60 ± 13% increase in the spectrally resolved surface perturbation amplitudes, consistent at all times and for all unsaturated frequencies. This consistency indicates that the increase in perturbation amplitudes was due to a change in the initial amplitudes rather than Rayleigh-Taylor growth. This is supported by the trajectory of individual modes and the inferred time-averaged perturbation growth rates. Laser imprint is therefore inferred to have increased due to strong magnetic fields that remain in the conduction zone of the target, suppressing off-axis electron motion and limiting the effects of thermal smoothing.

direct drive↗

Scalar mixing in a Kelvin-Helmholtz shear layer and implications for Reynolds-averaged Navier-Stokes modeling of mixing layers

Large-eddy simulation of a temporally evolving Kelvin-Helmholtz (KH) mixing layer is performed with the tenth-order compact difference code miranda to examine the steady-state behavior of a passive scalar in a shear-driven mixing layer. It is shown that the integral behavior of scalar variance in a KH mixing layer behaves similarly to the integral behavior of scalar variance in a Rayleigh-Taylor (RT) mixing layer, and mixedness of the simulated KH shear layer tends towards a value of about 0.8. It is further shown that if the k-L-a-V Reynolds-averaged Navier-Stokes (RANS) model [B. E. Morgan et al., Phys. Rev. E 98, 033111 (2018)], calibrated to reproduce steady-state mixing in an RT layer, is applied to simulate a KH mixing layer, the RANS model will significantly overpredict the magnitude of scalar variance in the KH layer. A straightforward addition to the k-L-a-V model is then suggested, and self-similarity analysis is applied to determine constraints on model coefficients. Furthermore, it is shown that with the addition of a buoyancy production term in the model equation for scalar variance, it becomes possible to eliminate the model deficiency and match steady-state mixedness in simulations of both RT and KH mixing layers with a single model calibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗