Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Equation-of-state measured via x-ray phase contrast imaging for Epon 828/DEA epoxy

Epoxies are a broad class of polymer materials often used as adhesive, structural or binding materials. Epon 828 is an epoxy resin that can be polymerized with a variety of curing agents with the choice of curing agent potentially having an effect on the resulting epoxy polymer’s material properties. In this study, the dynamic behavior of Epon 828 epoxy resin cured with diethanolamine (DEA) is investigated through a series of tamped Richtmyer-Meshkov instability (RMI) experiments measured with x-ray phase-contrast imaging. The measured shock and particle velocities are combined with data in the literature to calibrate Mie-Grüneisen equations-of-state (EOS) for portions and combinations of the collective dataset. The calibrated Mie-Grüneisen EOS are validated against particle velocity profiles extracted from published literature using the Eulerian hydrocode CTH. Here, the Mie-Grüneisen EOS fit to only the tamped RMI experimental data presented here most closely follows the particle velocity profile in the published literature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Extended magnetohydrodynamics simulations of thin-foil Z-pinch implosions with comparison to experiments

Cylindrical foil liners, with foil thicknesses on the order of 400 nm, are often used in university-scale Z-pinch experiments (~1 MA in 100 ns) to study physics relevant to inertial confinement fusion efforts on larger-scale facilities (e.g., the MagLIF effort on the 25-MA Z facility at Sandia National Laboratories). The use of ultrathin foil liners typically requires a central support rod to maintain the structural integrity of the liner target assembly prior to implosion. The radius of this support rod sets a limit on the maximum convergence ratio achievable for the implosion. In recent experiments with a support rod and a preimposed axial magnetic field, helical instability structures in the imploding foil plasma were found to persist as the foil plasma stagnated on the rod and subsequently expanded away from the rod. We have now used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) to study these experiments. The results suggest that it is the support rod which is responsible for the helical structures persisting beyond stagnation. Furthermore, we find that as the radius of the support rod decreases (i.e., as the convergence ratio increases), the integrity and persistence of the helical modes diminish. In the limit with no support rod, we find that the structure of the final stagnation column is governed by the structure of the central precursor plasma column. Furthermore, these simulation results and their comparisons to experiment are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗
DOE OSTI2021

Liner implosion experiments driven by a dynamic screw pinch

This paper expands upon recent experimental results [Campbell et al., Phys. Rev. Lett. 125, 035001 (2020)], where thin-foil liner implosions were driven by a dynamic screw pinch (DSP) and found to have magneto-Rayleigh–Taylor instability (MRTI) amplitudes up to three times smaller than in implosions driven by a standard z-pinch (SZP). Here, the expanded discussion presented herein includes: (1) a detailed comparison of the MRTI growth measured in the experiment with that calculated from theory; (2) measurements of axial magnetic field injection into the liner interior prior to the implosion, as well as the subsequent compression of this field during the implosion; (3) an in-depth description of how the helical geometry of the DSP can result in earlier implosion and stagnation times relative to the SZP; and (4) particle-in-cell simulations showing different electron drift behavior in the anode–cathode gap of the DSP relative to the SZP, and how this difference may be related to the different current waveforms recorded during the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A mechanism for reduced compression in indirectly driven layered capsule implosions

High-yield implosions on the National Ignition Facility rely on maintaining low entropy in the deuterium–tritium fuel, quantified by its adiabat, in order to efficiently couple energy to the hot spot through high compression of the fuel layer. We present very-high-resolution xRAGE simulation results that study the impacts of interfacial mixing and the jetting of materials due to surface defects, defects on internal interfaces, voids, and engineering features on fuel layer compression. Defects and voids are typically neglected in implosion simulations due to their small size and three-dimensional geometry. Our results showed that supersonic jets of material arise through weak spots in the shell at peak implosion velocity that prevent uniform compression of the fuel layer even when they do not introduce contaminant into the hot spot. This occurs despite maintaining low fuel entropy, since the formation of the weak spots involves nonradial displacement of fuel mass. In contrast, simulations show that fuel–ablator mixing due to interfacial instabilities has a much smaller impact on compression. We show that defects on interior interfaces of plastic capsules decrease compression by 15% to 25% and interfacial mixing between the ablator and fuel decreases compression by less than 1% for implosions with plastic or high-density carbon (HDC) ablators. For low adiabat implosions, the impact of jetting seeded by the support tent can also decrease the compression by 25%. We demonstrate that the inclusion of interior defects in simulations can explain the inferred compression in two fielded plastic capsule implosions and that the inclusion of voids, for which available characterization has large uncertainties, in simulations of HDC capsule implosions has a qualitatively consistent impact. This mechanism offers a potential explanation for persistently overestimated fuel compression in design simulations of layered implosions on the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids

For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. [Formula: see text] is the radially directed acceleration at the interface. ρ i denotes the density, μ i is the viscosity, and R i is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio [Formula: see text], density ratio [Formula: see text], spherical harmonic mode n, [Formula: see text], [Formula: see text], and the dimensionless growth rate [Formula: see text], where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ~ 1, low order harmonics can become the most unstable only for s >> 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Validating methods for modeling composition gradients in planar shock experiments

An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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↗
DOE OSTI2023

Non-Richtmyer–Meshkov instability ejecta production based on shallow bubble collapse

The study of shock-driven ejecta production has focused on Richtmyer–Meshkov instability (RMI) growth from geometric features of the material surface. Extensive study of this mechanism under both single- and multiple-shock conditions has found that the ejected mass tends to be closely associated with the shocked surface phase, and its temperature is not dramatically greater than the hydrodynamic shock temperature of the bulk. In this work, we propose and demonstrate a new ejecta production mechanism that can occur under multiple-shock conditions based on the collapse of bubbles near the free surface of the material. This mechanism produces ejected mass that is much greater in quantity than observed in the RMI case. The particles are much hotter than predicted by the shock Hugoniot state, and the ejected mass does not appear to be strongly dependent upon initial surface finish. The ejecta source extends into the material with no clear remaining free surface. We name this mechanism Shallow Bubble Collapse (SBC) and discuss the conditions under which it activates. We demonstrate resolved modeling methods that enable the calculation, design, and study of SBC as a mechanism and perform a series of experiments to compare with the models. Under some multiple-shock conditions, SBC ejection produces ten times more ejected mass than RMI growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Linear and Nonlinear Kelvin–Helmholtz Instability and Magnetohydrodynamic Wave Emission in Sheared Astrophysical Plasma Flows

The evolution of the Kelvin–Helmholtz instability (KHI) and magnetohydrodynamic (MHD) wave emission is investigated at shear-flow boundaries of magnetized plasmas. While MHD wave emission has been suggested to be only possible during the nonlinear stages, we find that there is also significant wave emission during the KHI's linear stages. These emitted MHD waves may have stronger impacts than KHI surface waves since they can act to transport energy away from the local region of the shear flow. The removal of energy from the shear-flow region, instead of just the local redistribution considered in previous studies, and its propagation away from the interface could have major implications for the evolution of astrophysical objects characterized by fast plasma flow shears.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Code-to-code comparison between FLASH and HYDRA in gas-puff Z-pinch modeling

The numerical modeling of gas-puff Z pinches involves the nonlinear coupling of a broad range of complex, multi-physics phenomena that makes such simulations challenging. The challenge is further compounded by nonlinear instabilities that can impact the dynamics of imploding gas-puff Z pinches, such as the magneto Rayleigh–Taylor instability (MRTI). If the growth rate and amplitude of the latter is comparable to the relevant timescales and properties of the imploding plasma, the MRTI can dramatically alter implosion dynamics, dictate pinch stability, and govern the plasma properties achievable in pulsed-power-driven laboratory experiments. National Laboratories and academic teams have developed numerical tools that can accurately model Z-pinch configurations and provide reliable design capabilities that can guide experimental choices and assist in interpreting experimental results. Most such tools, however, are not broadly available. Here, we present newly developed Z-pinch simulation capabilities of the publicly available FLASH code, applied in the study of MRTI growth and dynamical effects in gas-puff implosions. To verify the new implementations, we perform a comparison of FLASH gas-puff implosion simulations with previously published calculations with the HYDRA code from Lawrence Livermore National Laboratory, which have been validated with experimental data from the CESZAR pulsed-power driver at the University of California, San Diego. The experiments involved double- and triple-nozzle configurations, in an experimental attempt to stabilize the pinch to the MRTI. The code-to-code comparison shows similar results between the FLASH and HYDRA simulations, supporting the use of FLASH in the modeling of future gas-puff Z-pinch experiments at CESZAR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
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 OSTI2021

Detrimental effects and mitigation of the joint feature in double shell implosion simulations

Double shell capsules provide an attractive option in inertial confinement fusion experiments due to their potential for achieving a low-convergence, robust burn. However, these designs suffer from symmetry degradation and accompanying reduced fuel confinement due to the currently necessary joint between the two hemispheres of the outer shell. Here, the gap widens as a result of the excess ablation pressure produced by x rays that penetrate the joint during the drive phase, and this perturbation grows and imprints onto the inner shell during the collision. xRAGE Eulerian radiation-hydrodynamic simulations predict significant reductions in deuterium–tritium fusion yields compared to joint-less simulations when the depth of the outer joint is increased, whereas the performance is less sensitive to the depth of the inner gap. Here we examine the technique of plating the insides of the outer gap with a high-Z material to mitigate the impact of this feature. Gold-plating in quantities comparable to or exceeding the “missing” outer shell mass shows promise toward restoring both implosion symmetry and yield closer to the joint-less levels, and synthetic diagnostics suggest that high-energy x-ray radiographs can capture this shape retention of the inner and outer shells in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Entropy generation from hydrodynamic mixing in inertial confinement fusion indirect-drive targets

The increase in entropy from the physical mixing of two adjacent materials in inertial confinement fusion (ICF) implosions and gas-filled hohlraums is analytically assessed. An idealized model of entropy generation from the mixing of identical ideal-gas particles across a material interface in the presence of pressure and temperature gradients is applied. Physically, mix-driven entropy generation refers to the work done by the gases in expanding into a larger common volume from atomic mixing under the condition of no internal energy change, or work needed to restore the initial unmixed state. Furthermore, the effect of a mix-generated entropy increase is analytically shown to lead to less compression of the composite ICF fluid under adiabatic conditions. The amount of entropy generation is estimated to be ~10 J for a Rayleigh–Taylor-induced micrometer-scale annular mixing layer between the solid deuterium–tritium fuel and (undoped) high-density carbon pusher of an imploding capsule at the National Ignition Facility (NIF). This level of entropy generation is consistent with lower-than-expected fuel compressions measured on the NIF [Hurricane et al., Phys. Plasmas 26, 052704 (2019)]. The degree of entropy increase from mixing of high-Z hohlraum wall material and low-Z, moderate- to high-density gas fills is estimated to lead to ~100 kJ of heat generation for NIF-scale experiments [Moody et al., Phys. Plasmas 21, 056317 (2014)]. This value represents a significant fraction of the inferred missing x-ray drive energy based on observed delays in capsule implosion times compared with mainline simulations [Jones et al., Phys. Plasmas 19, 056315 (2012)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2021

A pushered capsule implosion as an alternate approach to the ignition regime for inertial confinement fusion

We report in inertial confinement fusion, the threshold for ignition is a highly dynamic quantity as the sources and sinks of power in the hot spot can vary rapidly. In this article, we consider the ignition condition as a race between heating and disassembly rates and make use of a prior solution to the fusion hot-spot thermodynamics to develop a Lawson-like ignition criteria for pressure × confinement time (p-τ) vs temperature. Low-Z capsule designs reach the temperature for this threshold using as much of the shell as feasible as ablator but then are limited in τ by low stagnated mass. An alternate approach, the pushered single shell (PSS) design [D. D.-M. Ho, S. MacLaren, and Y. Wang, “High-yield implosions via radiation trapping and high rho-R,” paper presented at the 60th Annual Meeting of the APS Division of Plasma Physics, 2018], introduces a dense inner layer of Mo-Be alloy that is smoothly graded outward to pure Be, increasing the confinement time at stagnation and lowering the temperature requirement at the ignition threshold. Here, we describe a PSS ignition design for the National Ignition Facility and use the theory as well as simulations to compare it with the low-Z capsule approach. Additionally, we show how an adjustment to the design is used to anticipate the effects of mixing at the fuel–ablator interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗