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

Fluctuating hydrodynamics and the Rayleigh–Plateau instability

The Rayleigh–Plateau instability occurs when surface tension makes a fluid column become unstable to small perturbations. At nanometer scales, thermal fluctuations are comparable to interfacial energy densities. Consequently, at these scales, thermal fluctuations play a significant role in the dynamics of the instability. These microscopic effects have previously been investigated numerically using particle-based simulations, such as molecular dynamics (MD), and stochastic partial differential equation–based hydrodynamic models, such as stochastic lubrication theory. In this paper, we present an incompressible fluctuating hydrodynamics model with a diffuse-interface formulation for binary fluid mixtures designed for the study of stochastic interfacial phenomena. An efficient numerical algorithm is outlined and validated in numerical simulations of stable equilibrium interfaces. We present results from simulations of the Rayleigh–Plateau instability for long cylinders pinching into droplets for Ohnesorge numbers of Oh = 0.5 and 5.0. Both stochastic and perturbed deterministic simulations are analyzed and ensemble results show significant differences in the temporal evolution of the minimum radius near pinching. Short cylinders, with lengths less than their circumference, were also investigated. As previously observed in MD simulations, we find that thermal fluctuations cause these to pinch in cases where a perturbed cylinder would be stable deterministically. Finally, we show that the fluctuating hydrodynamics model can be applied to study a broader range of surface tension–driven phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Hybrid target design for imprint mitigation in direct-drive inertial confinement fusion

A target design for mitigating the Rayleigh-Taylor instability is proposed for use in high energy density and direct-drive inertial confinement fusion experiments. In this scheme, a thin gold membrane is offset from the main target by several-hundred microns. A strong picket on the drive beams is incident upon this membrane to produce x rays which generate the initial shock through the target. The main drive follows shortly thereafter, passing through the ablated shell and directly driving the main target. The efficacy of this scheme is demonstrated through experiments performed at the OMEGA EP facility, showing a reduction of the Rayleigh-Taylor instability growth which scales exponentially with frequency, suppressing development by at least a factor of 5 for all wavelengths below 100 μm. Furthermore, this results in a delay in the time of target perforation by ~40%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stemflow Hydrodynamics

Stemflow hydrodynamics is the study of water movement along the exterior surface area of plants. Its primary goal is to describe water velocity and water depth along the stem surface area. Its significance in enriching the rhizosphere with water and nutrients is not in dispute. Yet, the hydrodynamics of stemflow have been entirely overlooked. This review seeks to fill this knowledge gap by drawing from thin film theories to seek outcomes at the tree scale. The depth‐averaged conservation equations of water and solute mass are derived at a point. These equations are then supplemented with the conservation of momentum that is required to describe water velocities or relations between water velocities and water depth. Relevant forces pertinent to momentum conservation are covered and include body forces (gravitational effects), surface forces (wall friction), line forces (surface tension), and inertial effects. The inclusion of surface tension opens new vistas into the richness and complexity of stemflow hydrodynamics. Flow instabilities such as fingering, pinching of water columns into droplets, accumulation of water within fissures due to surface tension and their sudden release are prime examples that link observed spatial patterns of stemflow fronts and morphological characteristics of the bark. Aggregating these effects at the tree‐ and storm‐ scales are featured using published experiments. The review discusses outstanding challenges pertaining to stemflow hydrodynamics, the use of dynamic similarity and 3D printing to enable the interplay between field studies and controlled laboratory experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The rate of development of atomic mixing and temperature equilibration in inertial confinement fusion implosions

The MARBLE project is a novel inertial confinement fusion platform for studying the development of atomic mixing and temperature equilibration in inertial confinement fusion implosions and their impact on thermonuclear burn. Experiments involve the laser-driven implosion of capsules filled with deuterated engineered foams whose pores are filled with a gaseous mixture of hydrogen and tritium. By varying the size of the foam pores, we can study the timescale of the development of atomic mix relative to the development of thermal equilibrium between species. In contrast, previous separated reactant experiments have only provided information on the total amount of mix mass. Additionally, we report on the series of MARBLE experiments [first reported in Haines et al., Nat. Commun. 11, 544 (2020)] performed on the University of Rochester's OMEGA laser facility and detailed and highly resolved three-dimensional radiation-hydrodynamic simulations of the implosions. In both the experimental and simulation results, we observe that the reactants do not achieve thermal equilibrium during the course of the implosion except in atomically mixed regions—i.e., that atomic mixing develops faster than thermal equilibration between species. The results suggest that ion temperature variations in the mixture are at least as important as reactant concentration variations for determining the fusion reaction rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Exploring capabilities of Micro-Fabricated 3D-printed capsules for studying effects of material mix on thermonuclear burn

The introduction of a high-precision, 3D-printing method, known as two-photon polymerization (2PP), has created a new means of producing novel targets to study inertial confinement fusion (ICF). Here, this work aims to explore if 2PP fabricated capsules can meet ICF requirements such as size, uniformity, deuteration degree, and ability to maintain gas fill. Two types of capsules with 3D-printed lattices were evaluated on the Omega-60 Laser Facility: (1) carbon-deuterium-oxygen (CDO) lattice with a H 2 gas fill and (2) carbon-hydrogen-oxygen (CHO) lattice with a D 2 gas fill. Experimental results and numerical simulations, which assumed complete mixing of capsule materials, reasonably agreed for both target types. A further study into the effects of capsule preheat or 2PP geometry on material mix is underway. Our work demonstrates that the unique capabilities of customizable, 3D-printed capsules present promising opportunities for further investigation into more sophisticated mix and burn experiments in ICF.

2PP↗
DOE OSTI2021

Understanding the evolution of liquid and solid microjets from grooved Sn and Cu samples using radiography

Experiments were performed on grooved Sn and Cu samples to study the temporal evolution of microjets. Jets were generated by the impact of gun-launched flyer plates against the back of grooved targets made from either Cu or Sn (groove depth of ~250 μm). The Hugoniot states in the various Sn targets encompassed conditions where solid phases are maintained throughout (7 and 16 GPa) and also conditions where melting occurs upon the release of compression (25 and 34 GPa); the transition occurs near a Hugoniot pressure of 23 GPa. Cu targets at 27 and 56 GPa provide comparisons in which the jets move at similar speeds but remain solid. In all cases, the spatial distribution of mass within the microjets was measured using high-speed synchrotron radiography. The result is a time history of the jet thickness profile from which quantities like total jet mass and jet velocity can be derived. In both the solid and liquid states, we generally observe that an increase in the shock strength leads to an increase in jet mass. However, this trend breaks down for Hugoniot states near the transition from continuously solid to melted-on-release. This is evidenced by the observation that there was no difference in the rate of mass flow in Sn jets at 16 and 25 GPa, while similar pressure jumps on either side of this range caused substantial changes in the jet mass. This contrasts with the behavior of smaller polishing defects that were present on the same samples (~1 μm deep). From these, no ejecta mass was detected below the melt boundary, but obvious microjets were generated once melting occurred. This indicates that crossing the bulk melt-on-release threshold can alternately promote or inhibit the flow of mass into microjets based on the amplitude of the initial perturbation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗
DOE OSTI2022

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal (rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multiphase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(\sigma _{\rho /\langle \rho \rangle }^{2\Gamma }-1)^{1/2}(1-\sigma _{\rho /\langle \rho \rangle }^{-2})^{-1/2}\mathcal {M}^{-1}\Gamma ^{-1/2}$, with the density dispersion σρ/<ρ>, the turbulent Mach number $\mathcal {M}$, and the polytropic exponent Γ. Purely solenoidal and purely compressive driving correspond to b ~ 1/3 and b ~ 1, respectively. As a result, using simulations in which a shock is driven into a multiphase medium with structures of different sizes and Γ < 1, we find b ~ 1 for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Localized mix-induced radiative cooling in a capsule implosion at the National Ignition Facility

In this work, w present direct measurements of electron temperature variations within an inertially confined deuterium-tritium plasma caused by localized mix of higher-Z materials into the central hot spot. The data are derived from newly developed differentially filtered penumbral imaging of the bremsstrahlung continuum emission. Our analysis reveals distinct localized emitting features in the stagnated hot-spot plasma, and we infer spatial variations in the electron temperature: the mixed region is 660±130eV colder than the surrounding hot-spot plasma at 3.26±0.11keV. Our analysis of the energy flow shows that we measure approximately steady-state conditions where the radiative losses in the mix region are balanced by heat conduction from the surrounding hot deuterium-tritium plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Anomalous Absorption by the Two-Plasmon Decay Instability

Radiation-hydrodynamic simulations of directly driven fusion experiments at the Omega Laser Facility predict absorption accurately when targets are driven at low overlapped laser intensity. Discrepancies appear at increased intensity, however, with higher-than-expected laser absorption on target. Strong correlations with signatures of the two-plasmon decay (TPD) instability—including half-harmonic and hard-x-ray emission—indicate that TPD is responsible for this anomalous absorption. Scattered light data suggest that up to ≈ 30 % of the laser power reaching quarter-critical density can be absorbed locally when the TPD threshold is exceeded. A scaling of absorption versus TPD threshold parameter was empirically determined and validated using the laser–plasma simulation environment code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multidimensional Radiation Hydrodynamics Simulations of Pulsational Pair-instability Supernovae

Stars with masses of 80–130 M ⊙ can encounter pulsational pair-instability at the end of their lives, which triggers consecutive episodes of explosive burning that eject multiple massive shells. Collisions between these shells produce bright transients known as pulsational pair-instability supernovae (PPI SNe) that may explain some extreme supernovae. In this paper, we present the first 2D and 3D radiation hydrodynamics simulations of PPI SNe with the CASTRO code. Radiative cooling causes the collided shells to evolve into thin, dense structures with hot spots that can enhance the peak luminosity of the SN by factors of 2–3. The light curve peaks at 1.9–2.1 × 10 43 erg s -1 for 50 days and then plateaus at 2–3 × 10 42 erg s -1 for 200 days, depending on the viewing angle. The presence of 12 C and 16 O and the absence of 28 Si and 56 Fe in its spectra can uniquely identify this transient as a PPI SN in follow-up observations. Our models suggest that multidimensional radiation hydrodynamics is required to model the evolution and light curves of all shell-collision SNe, such as Type IIne, not just PPI SNe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Further developments on observations of the Taylor instability in a dusty plasma

We previously reported the observation of a hydrodynamic-type interchange instability occurring sporadically at the curved interface between two dust clouds of different densities, in a three-dimensional dusty plasma [Pacha et al., Phys. Plasmas 19, 014501 (2012)]. The mechanism for the instability was not conclusively identified. Further analysis of the video images has been performed, showing that the instability is associated with the appearance of large amplitude dust acoustic waves or shocks. Evidence of a dust acoustic shock at the boundary, which could drive a Richtmyer–Meshkov-type instability, will be presented. We also discuss the possibility, suggested by Avinash and Sen [Phys. Plasmas 22, 083707 (2015)], that the observed instability is due to a Rayleigh–Taylor instability driven by a component of the gravitational acceleration, which is perpendicular to the curved boundary.

Pacha, K. A.↗
DOE OSTI2020

Supernova Hydrodynamics: A Lab-scale Study of the Blast-driven Instability Using High-speed Diagnostics

A novel experimental approach to study the blast-driven instability at a nondiffuse, gaseous interface with a density gradient is presented. Under Euler similarity, this approach enables study of dissipative-scale hydrodynamics relevant to many astrophysical and laboratory high energy density phenomena in a well-resolved manner. The instability is initiated by passing a Taylor–Sedov blast wave originating from a controlled detonation through a perturbed and stably stratified interface between two gases. The facility and driving blast wave are characterized to obtain repeatable conditions and capture large ensembles of time-resolved Mie scattering imaging that show consistent hydrodynamic development. We analyze the instability evolution between different gas pairs to demonstrate the wide range of development and turbulent behavior that may occur between different supernova layers. The mean evolution of the hydrodynamic instability is compared to a buoyancy–drag model that is frequently used to estimate perturbation growth in supernova mixing research. We propose a time delay to this model in order to reproduce the measured instability behavior and demonstrate model robustness in handling flows driven by a time-varying acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advancing our understanding of instabilities in high-energy-density systems through the study of benchmark datasets with advanced modeling tools

Numerous types of pulsed power driven inertial confinement fusion (ICF) and high energy density (HED) systems rely on implosion stability to achieve desired temperatures, pressures, and densities. Sandia National Laboratories Pulsed Power Sciences Center’s main ICF platform, Magnetized Liner Inertial Fusion (MagLIF), suffers from implosion instabilities which limit attainable fuel conditions and can compromise fuel confinement. This Truman Fellowship research primarily focused on computationally exploring (a) methods for improving our understanding of hydrodynamic and magnetohydrodynamic instabilities that form during cylindrical liner implosions, (b) methods for mitigating implosion instabilities, particularly those that degrade performance of MagLIF targets, and (c) novel MagLIF target designs intended to improve target performance primarily via enhanced implosion stability. Several multi-dimensional computational tools were used, including the magnetohydrodynamics code ALEGRA, the radiation-magnetohydrodynamics code HYDRA, and the magnetohydrodynamics code KRAKEN. This research succeeded in executing and analyzing simulations of automagnetizing liner implosions, shockless MagLIF implosions, dynamic screw pinch driven cylindrical liner implosions, and cylindrically convergent HED instability studies. The methods and tools explored and developed in this Truman Fellowship research have been published in several peer-reviewed journal articles and will serve as useful contributions to the fields of pulsed power science and engineering, particularly pertaining to pulsed power ICF and HED science.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗