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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

X-ray Pump-Probe Measurements using a Laser-Plasma Accelerator (Final Report)

This project funded one graduate student, Mario Balcazar, to perform experiments on high power laser facilities to study the application of X-rays generated by laser wakefield accelerated electrons for advanced radiography, in particular X-ray phase contrast imaging of hydrodynamic shocks/instabilities. The work made use of the DOE’s LaserNetUS facilities. We performed dynamic phase-contrast X-ray imaging of the shock generated by a 200 ps duration laser pulse (E = 1 J, I = 10 15 Wcm –2 ) with a 30 µm diameter liquid target, with unprecedented spatio-temporal resolution. This includes multi-shock generation within the liquid jet and plasma instability formation. Innovative electron-beam radiography was used to probe the laser-plasma interplay finding evidence of bilateral heating of the water followed by strong electric field generation. These measurements help explain some of the discrepancies between simulation and experiment and pave the way to better plasma diagnostic systems in HED and ICF physics experiments. This work was performed in collaboration with researchers from Michigan, LLNL, LBL and Sandia National Laboratories, Queen’s University Belfast and Imperial College.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. Here, the analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they become core-collapse supernovae (SNe). The resulting dense circumstellar medium may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomenon were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward-moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2–4 km s -1 convective motions seed the intrinsic convective instability associated with the high-entropy material deep in the envelope, enabling gas from deep within the envelope to escape and increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by ≈10× along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Margin to onset of nucleate boiling and flow instability studies for preliminary MITR design-demonstration element thermal-hydraulics

This study covers the onset of flow instability (OFI) preliminary results obtained from leveraging correlations, in addition to the preliminary thermal hydraulics results such as pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment (DDE) for the Massachusetts Institute of Technology Reactor (MITR). Current computational fluid dynamics (CFD) models in fluid structure interaction (FSI) have added the capability of assessing margins to onset of nucleate boiling (ONB). This study initiates the capability to model the margin to OFI and ONB presented for the MITR. Such study is supportive of the United States High Performance Research Reactor (USHPRR) program. Previous studies provided preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR DDE under conservative approximations for plate power distribution. In conclusion, this study focuses on providing insights into the OFI future research direction optimizing the transport of thermal energy, mass-flow rates, flow-channel geometries, and boundary conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

3D simulations of inertial confinement fusion implosions part 1: inline modeling of polarized cross beam energy transfer and subsequent drive anomalies on OMEGA and NIF

Inertial confinement fusion experiments are sensitive to cross-beam energy transfer (CBET), a laser-plasma instability that redistributes laser energy in the coronal plasma through self-generated ion acoustic wave (IAW) gratings. The detailed CBET coupling depends on the polarization state of the crossing wavefields. CBET itself can also scramble the beam polarizations by inducing ellipticity through the IAW grating, and rotating the seed polarization toward that of the pump. We develop a ray-based model that describes the polarized CBET coupling and that is compatible with the framework of 3D inline radiative hydrodynamics simulations. The model is implemented in the ASTER/IFRIIT code and verified against an academic test case and an offline polarized CBET post-processor. It is then applied to the detailed configuration of the distributed polarization rotator system on OMEGA, where results highlight how polarized CBET induces significant low modes in the collisional absorption source term. Finally, the modeling is applied to a simple indirect-drive configuration, comparing CBET calculations with 96 unpolarized or polarized beams with 24 unpolarized quads. It is shown that these cases produce similar power amplification per cone of beams grouped with similar polar angles. However, the 96 beam geometry itself is found to reduce azimuthal variations in quad power after the interaction and favors beams with larger polar angles within the cones, an effect that is amplified by the polarized CBET. Application of the model to inline calculations of OMEGA implosions are presented in a companion paper.

3D inline modeling↗
DOE OSTI2024

Systematic Local Simulations of Fast Neutrino Flavor Conversions with Scattering Effects

We investigate the dynamics of fast neutrino flavor conversions (FFCs) in the one-dimensional (1D) and zero-dimensional (0D) models, in which spatial advection is considered and ignored, respectively. In this study, we employ snapshots obtained by our self-consistent, realistic Boltzmann-neutrino-radiation-hydrodynamics simulations. We show that the FFC growth rate is considerably larger in the 1D model than in the 0D model, as expected from the previous linear analysis results. We find that the momentum space dimension does not significantly influence the neutrino transition probability in 1D models. On the other hand, in the 0D model without collisions, the FFC depends on the momentum space, and the azimuthal angle dependence breaks the periodicity of the FFC. Our study demonstrates that collisional instability can lead to further flavor conversions on a long timescale in 1D models after the asymptotic state of FFC has been reached. Such an effect should be taken into consideration when the fast and collisional flavor instabilities coexist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Reducing Richtmyer–Meshkov instability jet velocity via inverse design

In this work, we detail a novel application of inverse design and advanced manufacturing to rapidly develop and experimentally validate modifications to a shaped charge jet analog. The shaped charge jet analog comprises a copper liner, a high explosive (HE), and a silicone buffer. Here, we apply a genetic algorithm to determine an optimal buffer design that can be placed between the liner and the HE that results in the largest possible change in jet velocity. The use of a genetic algorithm allows for discoveries of unintuitive, complex, yet optimal buffer designs. Experiments using the optimal design verified the effectiveness of the buffer and validated the machine learning approach to hydrodynamic design optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electroconvection near an ion-selective surface with Butler–Volmer kinetics

We study the effects of interfacial kinetics on the electro-hydrodynamics of ion transport near an ion-selective surface using a combination of linear stability analysis and numerical simulation. The finite kinetics of the electrolyte–electrode interface affects the ion transfer and electroconvection in many ways. On a surface of fixed topography, such as a metal surface of slow and stable ion deposition or covered by a polymer membrane, the finite kinetics reduces the current in one-dimensional ion diffusion/migration, increases the critical voltage for the onset of the electroconvective instability, changes the dynamics of the electroconvection and the overlimiting current, and enhances the lateral ion diffusion within the interfacial layer. The first three effects are indirectly caused by the reaction kinetics and can be characterized by an effective voltage difference across the liquid electrolyte. In comparison, the last effect is controlled by a direct interplay between kinetics and nonlinear electroconvection. Scaling laws for ion transport and features of electroconvection are proposed. We also analyse the linear stability of a surface which evolves under ion deposition and find that the finite kinetics decreases the growth rate of both electroconvective and morphological instabilities and therefore modifies the wavenumber of the most unstable mode.

Mechanics↗
DOE OSTI2023

Radiative cooling effects on reverse shocks formed by magnetized supersonic plasma flows

We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental evidence of early-time saturation of the ion-Weibel instability in counterstreaming plasmas of CH, Al, and Cu

The collisionless ion-Weibel instability is a leading candidate mechanism for the formation of collisionless shocks in many astrophysical systems, where the typical distance between particle collisions is much larger than the system size. Multiple laboratory experiments aimed at studying this process utilize laser-driven (I≳10 15 W/cm 2 ), counterstreaming plasma flows (V≲ 2000 km/s) to create conditions unstable to Weibel-filamentation and growth. This technique intrinsically produces temporally varying plasma conditions at the midplane of the interaction where Weibel-driven B fields are generated and studied. Experiments discussed herein demonstrate robust formation of Weibel-driven B fields under multiple plasma conditions using CH, Al, and Cu plasmas. Linear theory based on benchmarked radiation-hydrodynamic FLASH calculations is compared with Fourier analyses of proton images taken ~5–6 linear growth times into the evolution. Our new analyses presented here indicate that the low-density, high-velocity plasma-conditions present during the first linear-growth time (~300–500 ps) sets the spectral characteristics of Weibel filaments during the entire evolution. It is shown that the dominant wavelength (~300 μm) at saturation persists well into the nonlinear phase, consistent with theory under these experimental conditions. However, estimates of B-field strength, while difficult to determine accurately due to the path-integrated nature of proton imaging, are shown to be in the ~10–30 T range, an order of magnitude above the expected saturation limit in homogenous plamas but consistent with enhanced B fields in the midplane due to temporally varying plasma conditions in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Exploration of cross-beam energy transfer mitigation constraints for designing an ignition-scale direct-drive inertial confinement fusion driver

The compression of direct-drive inertial confinement fusion (ICF) targets is strongly impacted by cross-beam energy transfer (CBET), a laser-plasma instability that limits ablation pressure by redirecting laser energy outward and that is projected to be mitigated by laser bandwidth. Here, we explore various CBET mitigation constraints to guide the design of future ICF facilities. First, we find that the flat, Gaussian, and Lorentzian spectral shapes have similar CBET mitigation properties, and a flat shape with nine spectral lines is a good surrogate for what can be obtained with other spectral shapes. Then, we conduct a comprehensive study across energy scales and ignition designs. 3D hydrodynamic simulations are used to derive an analytical model for the expected CBET mitigation as a function of laser and plasma parameters. From this model, we study the bandwidth requirements of conventional and shock ignition designs across four different energy scales and find that they require between 0.5 and 3±0.2% relative bandwidth. Further, best mitigation is achieved when the beam radius over critical radius R b /R c is kept low during the drive while the plasma temperature is kept high. In a steady state, we find that the bandwidth required to mitigate 85% of CBET scales as (R b /R c ) 2.15 L n –0.58 I 0.7 , where L n is the density scale length, and I the laser intensity. Finally, we find that the chamber beam port layout does not influence CBET mitigation. In the case of a driver using many monochromatic beamlets, we find that ~10 beamlets per port is required, with diminishing returns above ~20.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimental and computational study of thin-layer Rayleigh–Taylor instability development during deceleration with and without an externally applied magnetic field

The importance of mitigating the Rayleigh–Taylor instability (RTI) in inertial confinement fusion (ICF) is critical to successfully achieve high gain fusion yield. Consequently, understanding the seed mechanisms of RTI and the potential evolution of RTI in ICF relevant conditions is crucial. Single feature perturbations consistently demonstrate non-linear RTI evolution, for which an experimental platform on OMEGA-EP is developed. Manufacturing defects introduced into the target design require exploration of unanticipated changes to RTI development and an identification of targets that will still render quantifiable physics results. Consequently, it is presented that the inherent 3D nature of experimental targets necessitates 3D modeling for accurate design work and predictive modeling of experimental targets, especially when high resolution imaging diagnostics, like Fresnel Zone Plates, are utilized. A study of the morphology of the RTI evolution due to changing initial conditions and the presence of an externally applied magnetic field are also explored. Experimental data show thin-layer RTI morphology comparable to resistive magneto hydrodynamic 3D results. A discussion on the impacts of an externally applied magnetic field makes the case for continued efforts to observe a magnetic field's impact on RTI morphology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design of high-pressure iron Rayleigh–Taylor strength experiments for the National Ignition Facility

Iron is an important metal, scientifically and technologically. It is a common metal on Earth, forming the main constituent of the planet's inner core, where it is believed to be in solid state at high pressure and high temperature. It is also the main component of many important structural materials used in quasistatic and dynamic conditions. Laser-driven Rayleigh–Taylor instability provides a means of probing material strength at high pressure and high temperature. The unavoidable phase transition in iron at relatively low pressure induces microstructural changes that ultimately affect its strength in this extreme regime. This inevitable progression can make it difficult to design experiments and understand their results. Here, we address this challenge with the introduction of a new approach: a direct-drive design for Rayleigh–Taylor strength experiments capable of reaching up to 400 GPa over a broad range of temperatures. We use 1D and 2D hydrodynamic simulations to optimize target components and laser pulse shape to induce the phase transition and compress the iron to high pressure and high temperature. At the simulated pressure–temperature state of 350 GPa and 4000 K, we predict a ripple growth factor of 3–10 depending on the strength with minimal sensitivity to the equation of state model used. The growth factor is the primary observable, and the measured value will be compared to simulations to enable the extraction of the strength under these conditions. These experiments conducted at high-energy laser facilities will provide a unique way to study an important metal.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hydrodynamic fluctuations near a Hopf bifurcation: Stochastic onset of vortex shedding behind a circular cylinder

Here, we investigate hydrodynamic fluctuations in the flow past a circular cylinder near the critical Reynolds number Re c for the onset of vortex shedding. Starting from the fluctuating Navier-Stokes equations, we perform a perturbation expansion around Re c to derive analytical expressions for the statistics of the fluctuating lift force. Molecular-level simulations using the direct simulation Monte Carlo method support the theoretical predictions of the lift power spectrum and amplitude distribution. Notably, we have been able to collect sufficient statistics at distances Re ⁡/ Re c – 1 = O ⁡(10 –3 ) from the instability that confirm the appearance of non-Gaussian fluctuations, and we observe that they are associated with intermittent vortex shedding. These results emphasize how unavoidable thermal-noise-induced fluctuations become dramatically amplified in the vicinity of oscillatory flow instabilities and that their onset is fundamentally stochastic.

42 ENGINEERING↗
DOE OSTI2020

Ring patterns generated by an expanding colloidal meniscus

The drop-and-dry is a common technique allowing for creation of periodic nanoparticle (NP) structures for sensing, photonics, catalysis, etc. However, the reproducibility and scalability of this approach for fabrication of NP-based structures faces serious challenges due to the complexity of the simple, at first glance, evaporation process. In this work we study the effect of the spatial confinement on the NP self-assembly under slow solvent evaporation, when the air-liquid-substrate contact line (CL) expands from the center towards the walls of a cylindrical cell, forming a toroid. Using in situ video monitoring of the stick-slip CL motion, we find regular hydrodynamic perturbations in the meniscus, and reveal fine details of the formation of quasiperiodic rings of close packed NP layers. We report that drying of the toroidal NP droplet has a number of important differences from drying of the classical hemispherical colloidal drops. In toroidal drops we observe linear-in-time average meniscus motion, in contrast to the hemispherical drops where the meniscus moves as a square root of time. While both droplet geometries produce NP ring patterns, the ring width for the toroidal drop decreases with increasing ring radius, while it decreases with decreasing the radius of the hemispherical drop. We suggest that free ligands are the main cause of the Marangoni instabilities driving the periodic vorticity in the meniscus. Finally, we show that the usually ignored contact line tension may yield a considerable contribution to the CL pinning causing the CL slip-stick motion and the ring formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Laser pulse-length dependent ablation and shock generation in silicon at 5 × 10 14 W/ cm 2 intensities

The effect of laser pulse duration on energy coupling into a planar silicon target is investigated in experiments at the OMEGA-EP facility by varying the laser pulse length τ —spanning 3 orders of magnitude from 100 ps to 10 ns—while maintaining a constant peak laser intensity, I 0 = 5 × 10 14 W/ cm 2 . In theoretical models, the ablation pressure primarily scales for a given material with laser intensity and wavelength, which are all fixed variables here, allowing us to explore the specific role of laser pulse duration. Two-dimensional radiation-hydrodynamics simulations benchmarked with optical probing of the expanding plasma show that the pulse duration is critical for the ablation pressure to reach a steady state. Moreover, the pulse duration impacts shock decay and multiple wave effects, which strongly dictate the evolving shock profile that propagates within the laser-shocked target as ultimately measured by rear-surface diagnostics. The shock velocities inferred from the theoretical model, after considering shock decay, impedance matching, and shock Hugoniot, are found to be in good agreement with velocimetry measurements. However, discrepancies are observed with simulations for the shorter (0.1 ns) and longer (10 ns) pulse durations, which are respectively attributed to unaccounted contributions of kinetic absorption mechanisms and instabilities in simulations. Published by the American Physical Society 2024

Bailly-Grandvaux, M. (ORCID:0000000175294013)↗
DOE OSTI2022

Measurements of the imploding plasma sheath in triple-nozzle gas-puff z pinches

Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗