Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrodynamic 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 325 records · Page 18

Three-Dimensional Hydrodynamic Simulations of the Postimpact Proto-Earth

Astrophysical fluid configurations are susceptible to a variety of nonaxisymmetric instabilities under the combined effects of rotation, self-gravity, and thermal pressure. When strong enough, they can induce rapid transport of mass and angular momentum. Our own previous studies of nonaxisymmetric instabilities in model protostars and protostellar disks show that significant transport can occur on orbital timescales and that material can be ejected to large distances. In this contribution, we present three-dimensional simulations of the circumterrestrial debris belt that may have resulted from a giant impact. Our three-dimensional hydrodynamics code with self-gravity and artificial viscosity is fully second-order in space and time; the equations of hydrodynamics and the Poisson equation are solved on an Eulerian cylindrical grid. In the preliminary calculations presented here, we use a simplified EOS where the central proto-Earth is treated as an n = 1/2 polytropic fluid, surrounded by a more compressible, rapidly rotating, fluid disk that represents silicate vapor. Our initial disk parameters are generated from the endstate data of recent smoothed particle hydrodynamics giant-impact calculations. Ultimately, we wish to detennine under what conditions nonaxisymmetric instabilities grow in the postimpact disk and whether they facilitate the transport of material outside the proto-Earth's Roche Limit, leading to the formation of the Moon. In future work, we hope to include a more realistic EOS and the consequences of heating, cooling, and phase transitions.

Pickett, B. K.↗

Synchrotron radiography of Richtmyer–Meshkov instability driven by exploding wire arrays

We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine learning visualization tool for exploring parameterized hydrodynamics

We are interested in the computational study of shock hydrodynamics, i.e. problems involving compressible solids, liquids, and gases that undergo large deformation. These problems are dynamic and nonlinear and can exhibit complex instabilities. Due to advances in high performance computing it is possible to parameterize a hydrodynamic problem and perform a computational study yielding $\mathscr{O} (TB)$ of simulation state data. We present an interactive machine learning tool that can be used to compress, browse, and interpolate these large simulation datasets. This tool allows computational scientists and researchers to quickly visualize 'what-if' situations, perform sensitivity analyses, and optimize complex hydrodynamic experiments.

97 MATHEMATICS AND COMPUTING↗

A Shock Stabilization of the HLLC Riemann Solver for the Carbuncle Instability

The HLLC approximate Riemann solver improves upon the HLL Riemann solver by resolving contact discontinuities. This is a particularly desirable property for multi-material codes in which problems usually contain material interfaces. However, the HLLC solver is known to suffer from the carbuncle phenomenon, a numerical instability most apparent at grid-aligned shocks in multi-dimensional simulations. Many problems of interest, including high energy-density physics applications, require the accurate resolution of both material interfaces and hydrodynamic shocks. A variety of methods have been developed to cure this instability, with varying degrees of complexity. The objective of this work is to describe a simple approach to modify the HLLC Riemann solver and prevent the carbuncle instability. Here, the method is then demonstrated for assorted two-dimensional test problems known to exhibit the shock instability. The performance of the new solver is compared with that of the standard HLL and HLLC Riemann solvers.

97 MATHEMATICS AND COMPUTING↗

On spatial instability of a swirling jet - Theory and experiment

The spatial instability of a swirling jet is investigated both experimentally and theoretically. A hydrodynamic stability analysis is applied to an inviscid incompressible top-hat jet, with a swirl distribution of solid-body rotation and free vortex in and outside the vortex core, respectively. Both plane and helical instability modes are examined. It is found that the top-hat jet with swirl distribution of Rankine vortex type is unstable in all the modes studied. The higher the positive helicity, the less spatially unstable the jet behavior; the higher the negative helicity, the more spatially unstable this behavior becomes. A comparison is made between theoretical results and experimental data on a low-intensity swirling jet. The trend of the initial growth of the instability waves in the near field is captured.

Wu, C.↗

The acoustic characteristics of turbomachinery cavities

Internal fluid flows are subject not only to self-sustained oscillations of the purely hydrodynamic type but also to the coupling of the instability with the acoustic mode of the surrounding cavity. This situation is common to turbomachinery, since flow instabilities are confined within a flow path where the acoustic wavelength is typically smaller than the dimensions of the cavity and flow speeds are low enough to allow resonances. When acoustic coupling occurs, the fluctuations can become so severe in amplitude that it may induce structural failure of engine components. The potential for catastrophic failure makes identifying flow-induced noise and vibration sources a priority. In view of the complexity of these types of flows, this report was written with the purpose of presenting many of the methods used to compute frequencies for self-sustained oscillations. The report also presents the engineering formulae needed to calculate the acoustic resonant modes for ducts and cavities. Although the report is not a replacement for more complex numerical or experimental modeling techniques, it is intended to be used on general types of flow configurations that are known to produce self-sustained oscillations. This report provides a complete collection of these models under one cover.

Lucas, M. J.↗

Reconstructing Richtmyer–Meshkov instabilities from noisy radiographs using low dimensional features and attention-based neural networks

We develop an ML-based approach for density reconstruction based on transformer neural networks. This approach is demonstrated in the setting of ICF-like double shell hydrodynamic simulations wherein the parameters related to material properties and initial conditions are varied. The new method can robustly recover the complex topologies given by the Richtmyer-Meshkoff instability (RMI) from a sequence of hydrodynamic features derived from radiographic images corrupted with blur, scatter, and noise. A noise model is developed to characterize errors in extracting features from synthetic radiographs of the simulated density field. The key component of the network is a transformer encoder that acts on a sequence of features extracted from noisy radiographs. This encoder includes numerous self-attention layers that act to learn temporal dependencies in the input sequences and increase the expressiveness of the model. This approach is shown to exhibit an excellent ability to accurately recover the RMI growth rates, despite the gas-metal interface being greatly obscured by radiographic noise. Our approach can be applied in a broad array of fields involving shock physics and material science.

47 OTHER INSTRUMENTATION↗

[Kelvin-Helmholtz instability in protostellar jets]

NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.

Stone, James↗

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modern developments in shear flow control with swirl

Passive and active control of swirling turbulent jets is experimentally investigated. Initial swirl distribution is shown to dominate the free jet evolution in the passive mode. Vortex breakdown, a manifestation of high intensity swirl, was achieved at below critical swirl number (S = 0.48) by reducing the vortex core diameter. The response of a swirling turbulent jet to single frequency, plane wave acoustic excitation was shown to depend strongly on the swirl number, excitation Strouhal number, amplitude of the excitation wave, and core turbulence in a low speed cold jet. A 10 percent reduction of the mean centerline velocity at x/D = 9.0 (and a corresponding increase in the shear layer momentum thickness) was achieved by large amplitude internal plane wave acoustic excitation. Helical instability waves of negative azimuthal wave numbers exhibit larger amplification rates than the plane waves in swirling free jets, according to hydrodynamic stability theory. Consequently, an active swirling shear layer control is proposed to include the generation of helical instability waves of arbitrary helicity and the promotion of modal interaction, through multifrequency forcing.

Farokhi, Saeed↗

Modern developments in shear flow control with swirl

Passive and active control of swirling turbulent jets is experimentally investigated. Initial swirl distribution is shown to dominate the free jet evolution in the passive mode. Vortex breakdown, a manifestation of high intensity swirl, was achieved at below critical swirl number (S = 0.48) by reducing the vortex core diameter. The response of a swirling turbulent jet to single frequency, plane wave acoustic excitation was shown to depend strongly on the swirl number, excitation Strouhal number, amplitude of the excitation wave, and core turbulence in a low speed cold jet. A 10 percent reduction of the mean centerline velocity at x/D = 9.0 (and a corresponding increase in the shear layer momentum thickness) was achieved by large amplitude internal plane wave acoustic excitation. Helical instability waves of negative azimuthal wave numbers exhibit larger amplification rates than the plane waves in swirling free jets, according to hydrodynamic stability theory. Consequently, an active swirling shear layer control is proposed to include the generation of helical instability waves of arbitrary helicity and the promotion of modal interaction, through multifrequency forcing.

Farokhi, S.↗

Hydroscaling indirect-drive implosions on the National Ignition Facility

A goal of the laser-based National Ignition Facility (NIF) is to increase the liberated fusion energy “yield” in inertial confinement fusion experiments well past the ignition threshold and the input laser energy. One method of increasing the yield, hydrodynamic scaling of current experiments, does not rely on improving compression or implosion velocity, but rather increases the scale of the implosion to increase hotspot areal density and confinement time. Indirect-drive (Hohlraum driven) implosions carried out at two target sizes, 12.5% apart, have validated hydroscaling expectations. Moreover, extending comparisons to the best-performing implosions at five different capsule sizes shows that their performance also agrees well with hydroscaling expectations even though not direct hydroscales of one another. In the future, by switching to a reduced loss Hohlraum geometry, simulations indicate that we can drive 20% larger-scale implosions within the current power and energy limitations on the NIF. Finally, at the demonstrated compression and velocity of these smaller-scale implosions, these 1.2× hydroscaled implosions should put us well past the ignition threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The stability of cosmic-ray-dominated shocks - A secondary instability

Previous studies in one spatial dimension showed that the precursors of cosmic-ray-dominated shocks are unstable against traveling acoustic disturbances. Here we report that a secondary, Rayleigh-Taylor type instability can exist in cosmic-ray-dominated media influenced by the acoustic instability. Using the local WKB approximation, the growth rate of the secondary instability is shown to be comparable to that of the one-dimensional acoustic instability itself in the cases we have considered. The nonlinear development has been followed numerically with a two-dimensional PPM hydrodynamics code that also incorporates the two-fluid cosmic-ray energy equation. We show that the secondary instability may cause the precursor and postshock flows to become highly turbulent.

Ryu, Dongsu↗

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Postexplosion hydrodynamics of supernovae in red supergiants

Shock propagation, mixing, and clumping are studied in the explosion of red supergiants as Type II supernovae using a two-dimensional smooth particle hydrodynamic (SPH) code. We show that extensive Rayleigh-Talor instabilities develop in the ejecta in the wake of the reverse shock wave. In all cases, the shell structure of the progenitor is obliterated to leave a clumpy, well-mixed supernova remnant. However, the occurrence of mass loss during the lifetime of the progenitor can significantly reduce the amount of mixing. These results are independent of the Type II supernova explosion mechanism.

Herant, Marc↗

The importance of laser wavelength for driving inertial confinement fusion targets. I. Basic physics

We reinvestigate the role that laser wavelength plays in driving inertial confinement fusion (ICF) targets. Different assumptions underlie previous analytic frameworks that provide predictions for wavelength scaling of many important target parameters. These are explored and compared to radiation-hydrodynamics simulations of laser-driven targets. We are particularly interested here in lasers with wavelengths between 0.193 μm [wavelength of the Argon Flouride (ArF) excimer laser] and 0.527 μm (the frequency-doubled glass Nd:glass laser). Short-wavelength drivers have significant advantages for directly driven ICF targets, which are summarized here. We show that constraints such as providing a certain pressure or avoiding laser-plasma instability thresholds allow shorter laser wavelengths to provide energy savings, pressure enhancements, and/or higher hydrodynamic efficiencies. We also consider potential disadvantages, such as increased laser imprint or exposure to the Landau–Darrieus instability. These are shown to be either minor and/or can be easily remediated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrodynamic Dryout in Two-Phase Flows: Observations of Low Bond Number Systems

Dryout occurs readily in certain slug and annular two-phase flows for systems that exhibit partial wetting. The mechanism for the ultimate rupture of the film is attributed to van der Waals forces, but the pace towards rupture is quickened by the surface tension instability (Rayleigh-type) of the annular film left by the advancing slug and by the many perturbations of the free surface present in the Re(sub g) approximately 0(10(exp 3)), Re(sub l) approximately 0(10(exp 4)), and Ca approximately 0(10(exp -1) flows. Results from low-gravity experiments using three different test fluids are presented and discussed. For the range of tests conducted, the effect of increasing viscosity is shown to eliminate the film rupture while the decrease of surface tension via a surfactant additive is shown to dramatically enhance it. Laboratory measurements using capillary tubes are presented which reveal the sensitivity of the dryout phenomena to particulate and surfactant contamination. Rom such observations, dryout due to the hydrodynamic-van der Waals instability can be expected in a certain range of flow parameters in the absence of heat transfer. The addition of heat transfer may only exacerbate the problem by producing thermal transport lines replete with "hot spots." A caution to this effect is issued to future space systems designers concerning the use of partially wetting working fluids.

Weislogel, Mark M.↗

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗