Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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 217 records · Page 12

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Morphology of shocked lateral outflows in colliding hydrodynamic flows

Supersonic interacting flows occurring in phenomena, such as protostellar jets, give rise to strong shocks and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. Here, we introduce variations in the flow parameters of density, velocity, and cross-sectional radius of the colliding flows in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows being deflected from the horizontal when the flows are of differing cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using neutrons and x rays to measure plasma conditions in a solid sphere of deuterated polyethylene compressed to densities of 35 g/cc at temperatures of 2 keV and pressures of 40 Gbar

This paper describes an experiment that shock compresses the center of a solid deuterated polyethylene sphere, CD 2 , to densities of 35 g/cc and temperatures of 2 keV with corresponding pressure of 40 Gbar. The design employs a strong spherically converging shock launched through a solid ball of material using a Hohlraum radiation drive. As the shock coalesces at the center it produces a hot spot that we characterize by measuring the x-ray self-emission and 2.45 MeV neutrons emitted. Two-dimensional images and time-resolved measurements of the x rays emitted determine the size and time duration of the hot spot, leading to an estimated 2k eV electron temperature. The neutron time of flight spectrometer measures an average ion temperature of 1.06 ± 0.15 keV and neutron yield of 7.0 (±0.5) × 10 9 DD neutrons. Our new distribution function tool enables us to create a forward model of the experimental data based on 1D radiation-hydrodynamic simulations, leading to a better understanding of the plasma conditions that produce the measured neutrons and x rays. Furthermore, our simulations indicate that the x rays are produced in a short-lived hot-dense core over tens of picoseconds, whereas the neutron emission continues for about 200 ps, as the hot core starts to expand, thereby leading to a lower mean temperature of the plasma during neutron production. This finding is in agreement with the experimental data, and we therefore conclude that the forward-modeling is a useful tool for inferring the conditions of the hot spot in a laser-driven implosion during burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The colliding planar shocks platform to study warm dense matter at the National Ignition Facility

We have developed an experimental platform at the National Ignition Facility that employs colliding planar shocks to produce warm dense matter with uniform conditions and enable high-precision equation of state measurements. The platform uses simultaneous x-ray Thomson scattering and x-ray radiography to measure the density, electron temperature, and ionization state in warm dense matter. The experimental platform is designed to create a large volume of uniform plasma (approximately 700×700×150μm 3 ) at pressures approaching 100 Mbar and minimize the distribution of plasma conditions in the x-ray scattering volume, significantly improving the precision of the measurements. Here, in this study, we present the experimental design of the platform and compare hydrodynamic simulations to x-ray radiography data from initial experiments studying hydrocarbons, producing uniform densities within ±25% of the average probed condition. We show that the platform creates a homogeneous plasma that can be characterized using x-ray Thomson scattering. Thus, the new platform enables accurate measurements of plasma conditions necessary to test models for the equation of state and ionization potential depression in the warm dense matter regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Al–W gradient density materials—Processing and dynamic ramp compression

Materials with high-density gradients are desired for controlling loading paths in dynamic compression, important for studying material properties in extreme conditions and inertial confinement fusion. The large density difference between Al and W makes them ideal choices for producing gradient density materials, but their extremely different melting temperatures make them challenging to fabricate simultaneously. We report a method for producing Al–W porosity-free materials with a fourfold increase in density (2.7–11 g/cm 3 ) across the composition range, from Al-rich to W-rich, without intermetallic phase formation. This was achieved by understanding the aluminum-dominated densification behavior and examining the influence of pressure and temperature on the densification of Al–W composites. Dynamic compression experiments conducted with the Al–W gradient density material produced shock ramp compressions as expected based on the designed composition, and the performed hydrodynamics simulations showed excellent agreement with experimental results. The results demonstrate that current activated pressure-assisted densification allows for the easy and rapid fabrication of gradient density materials with significant density gradients and tailored compositions, facilitating precise control of the loading paths. These materials have the potential to create customized pressure drives for advancing the fields of material science in extreme environments and dynamic compression.

Alloys↗

Sensitivity of an integrated experiment to uncertainty in the high explosive equations of state

Traditionally, hydrodynamics simulations are performed with a single equation-of-state (EOS) to describe each material. These EOSs typically have a physics-informed functional form with adjustable parameters that are calibrated in order to replicate small-scale data. However, because the calibration data have uncertainty and there are typically inherent degeneracies in fitting the EOS, there are actually multiple EOSs that might be consistent with calibration data. In this work, we perform uncertainty quantification (UQ) for the reactant and product equations of state for the high explosive PBX 9501 to yield an ensemble of EOSs that match the uncertain small-scale calibration data. We then simulate an experiment of an explosively formed penetrator repeatedly with different EOSs to both validate the UQ analysis and determine the effects of EOS uncertainty on the prediction of quantities of interest in the experiment. In general, we find good agreement between the simulation predictions and the experimental measurements, and we identify an EOS variable that contributes most directly to the spread in the predictions as the EOSs are varied.

36 MATERIALS SCIENCE↗

Demonstration of low-mode shape control in indirect-drive double shell implosions at the NIF

Double shell inertial confinement fusion is a concept for achieving robust thermonuclear burn that uses dense metal shells to compress deuterium-tritium (DT) fuel to fusion conditions. Double shell implosions are typically indirectly driven and involve a target that consists of a low-Z ablator, a foam layer, and a high-Z pusher surrounding the DT fuel. The goal of the campaign is to achieve a volumetric burn as radiation losses from the DT fuel are trapped by the opaque high-Z shell. The overall performance of double shell implosions relies on the efficient collisional transfer of kinetic energy between layers. The efficiency of this transfer (and therefore the overall performance of a given implosion) is degraded by the presence of low-mode asymmetries. P2 asymmetries are often observed in spatially resolved 2D radiographs of nominal double shell implosions. This work discusses three such experiments: one with an oblate P2 asymmetry, one with a prolate P2 asymmetry, and one with an approximate spherical symmetry. After performing a shape analysis of the oblate and prolate implosions to quantify asymmetries, these experimental results were compared with the results of hydrodynamic simulations for the two experiments. Differences between the experiment and simulation were then used to design an approximately spherical implosion by altering the incident laser cone fraction. Radiographs from the experiment that implemented the modified cone fraction show evidence of an implosion that is approximately spherical until bang time. This design is intended to serve as a point design for future studies that will seek to optimize various aspects of the double shell target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First graded metal pushered single shell capsule implosions on the National Ignition Facility

Graded metal pushered single shell (PSS) capsules are predicted to be a viable alternative to low-Z capsule indirect drive inertial confinement fusion (ICF) implosions for achieving high fusion yields [MacLaren et al., Phys. Plasmas 28, 122710 (2021)]. The first experiments with Be/Cr-graded metal PSS capsules indicate that the implementation of the principle design feature, the graded density inner metal layer, has succeeded in producing a stable implosion with performance in agreement with predictions. With 50% Cr concentration in the pusher, PSS capsules have greater than ~2× higher shell densities during stagnation for enhanced core confinement and radiation trapping at ~35% lower shell implosion velocities than low-Z capsules. High-energy >30 keV inflight shell radiography recorded 215 km/s implosion velocities and show that implosion Legendre mode P 2 asymmetry can be tuned via inner-to-outer beam wavelength separation, similar to other implosions. Shell radiographs and neutron core images show similar P 2 asymmetry, suggesting no symmetry swings between peak implosion velocity and stagnation times. Despite the modest implosion velocities, gas-filled deuterium–tritium capsule implosions generate 10 15 neutron yields at relatively modest core ion temperatures of 2.75 keV, indicating that in spite of the high-density inner layer, the implosions have been stabilized by the design density gradient. When compared with hydrodynamic simulations, the measured yield-over-simulated is 35% due to fuel–pusher mix and other perturbations such as the capsule fill tube. Finally, simple analytical scalings of hot spot pressure and neutron yield show that PSS implosions reach similar performance at lower implosion velocities and higher shell densities to low-Z ICF capsules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Hot Circumsingle Disks Drive Binary Black Hole Mergers in Active Galactic Nucleus Disks

Abstract Binary black hole (BBH) mergers, particularly those with component masses in the pair-instability gap, may be produced by hierarchical mergers in the disks surrounding Active Galactic Nuclei (AGNs). While the interaction of an embedded BBH with an AGN disk is typically assumed to facilitate a merger, recent high-resolution hydrodynamical simulations challenge this assumption. However, these simulations often have simplified treatments for gas thermodynamics. In this work, we model the possible consequence of various feedback from an embedded BBH with a simple model that maintains an enhanced temperature profile around each binary component. We show that when the minidisks around each BH become hotter than the background by a factor of three, the BBH orbital evolution switches from expansion to contraction. By analyzing the gravitational torque profile, we find that this change in direction is driven by a weakening of the minidisk spirals and their positive torque on the binary. Our results highlight the important role of thermodynamics around BBHs and its effect on their orbital evolution, suggesting that AGN disks could be efficient factories for BBH mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

Three-dimensional hot-spot x-ray emission tomography from cryogenic deuterium–tritium direct-drive implosions on OMEGA

A three-dimensional model of the hot-spot x-ray emission has been developed and applied to the study of low-mode drive asymmetries in direct-drive inertial confinement fusion implosions on OMEGA with cryogenic deuterium–tritium targets. Here, the steady-state model assumes an optically thin plasma and the data from four x-ray diagnostics along quasi-orthogonal lines of sight are used to obtain a tomographic reconstruction of the hot spot. A quantitative analysis of the hot-spot shape is achieved by projecting the x-ray emission into the diagnostic planes and comparing this projection to the measurements. The model was validated with radiation-hydrodynamic simulations assuming a mode-2 laser illumination perturbation resulting in an elliptically shaped hot spot, which was accurately reconstructed by the model using synthetic x-ray images. This technique was applied to experimental data from implosions in polar-direct-drive illumination geometry with a deliberate laser-drive asymmetry, and the hot-spot emission was reconstructed using spherical-harmonic modes of up to ℓ = 3. A 10% stronger drive on the equator relative to that on the poles resulted in a prolate-shaped hot spot at stagnation with a large negative A 2,0 coefficient of A 2,0 = –0.47 ± 0.03, directly connecting the modal contribution of the hot-spot shape with the modal contribution in laser-drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tidal Stripping of a White Dwarf by an Intermediate-mass Black Hole

During the inspiralling of a white dwarf (WD) into an intermediate-mass black hole (~10 2-5 M ⊙ ), both gravitational waves (GWs) and electromagnetic (EM) radiation are emitted. Once the eccentric orbit's pericenter radius approaches the tidal radius, the WD would be tidally stripped upon each pericenter passage. The accretion of this stripped mass would produce EM radiation. It is suspected that the recently discovered new types of transients, namely the quasiperiodic eruptions and the fast ultraluminous x-ray bursts, might originate from such systems. Modeling these flares requires a prediction of the amount of stripped mass from the WD and the details of the mass supply to the accretion disk. We run hydrodynamical simulations to study the orbital parameter dependence of the stripped mass. We find that our results match the analytical estimate that the stripped mass is proportional to z 5/2 , where z is the excess depth by which the WD overfills its instantaneous Roche lobe at the pericenter. The corresponding fallback rate of the stripped mass is calculated, which may be useful in interpreting the individual flaring light curve in candidate EM sources. We further calculate the long-term mass-loss evolution of a WD during its inspiral and the detectability of the GW and EM signals. The EM signal from the mass-loss stage can be easily detected: the limiting distance is ~320(M h /10 4 M ⊙ ) Mpc for the Einstein Probe. The GW signal, for space-borne detectors such as Laser Interferometer Space Antenna or TianQin, can be detected only within the Local Supercluster (~33 Mpc).

79 ASTRONOMY AND ASTROPHYSICS↗

Evidence of restricted heat transport in National Ignition Facility Hohlraums

We present experimental evidence of restricted electron thermal conduction in the high-Z coronal plasma regions of laser-driven Hohlraums on the National Ignition Facility. Four separate measurements, three of which are direct observations of Hohlraum dynamics, corroborate this finding. (1) The velocity of the coronal plasma ablated and heated by the outer-cone laser beams is determined by time-dependent imaging of the gold plasma plume, or “bubble.” The velocities of the incoming plume (perpendicular to the Hohlraum axis) are consistent with high-fidelity 2D radiation-hydrodynamic simulations using flux-limited thermal electron conduction with a flux multiplier f = 0.03. Simulations using f = 0.15, which is very nearly classical Spitzer–Härm transport, predict plume velocities slower than measured. (2) Specific features in time-resolved images of the Hohlraum wall at an angle of 19° are also more consistent with f = 0.03 simulations compared to f = 0.15. (3) Spectroscopic tracers were added to the Hohlraum wall in the outer-beam bubble region. The ratios of hydrogen-like to helium-like line emission are sensitive to the electron temperature of the bubble. The hydrogen-like to helium-like ratios extracted from the time-integrated spectra of manganese and cobalt tracers from two observation angles are consistent with f = 0.03 and not with f = 0.15. (4) The time of peak capsule emission, or “bang time,” an integrated measurement, is also more consistent with f = 0.03 than with f = 0.15. Overall, while these findings do not identify the causes of restricted thermal conduction in Hohlraums, they motivate future experiments to test specific hypotheses and focus on model development in the regions of the plasma exhibiting restricted transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Galaxy cluster matter profiles - I. Self-similarity, mass calibration, and observable-mass relation validation employing cluster mass posteriors

We present a study of the weak lensing inferred matter profiles ΔΣ(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel’dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles ΔΣ(R/R200c)/(ρcritR200c) with a lower dispersion than the unscaled ΔΣ(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M200|ζ̂, λ̂, z), which describes the individual cluster masses given their tSZE (ζ̂) and optical (λ̂, z) observables. This method enables simultaneous constraints on richness λ-mass and tSZE detection significance ζ-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPT×DES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPT×DES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the λ-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical ΛCDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.Key words: gravitational lensing: weak / galaxies: clusters: general / large-scale structure of Universe

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental characterization of hot-electron emission and shock dynamics in the context of the shock ignition approach to inertial confinement fusion

We report on planar target experiments conducted on the OMEGA-EP laser facility performed in the context of the shock ignition (SI) approach to inertial confinement fusion. The experiment aimed at characterizing the propagation of strong shock in matter and the generation of hot electrons (HEs), with laser parameters relevant to SI (1-ns UV laser beams with I ~10 16 W/cm 2 ). Time-resolved radiographs of the propagating shock front were performed in order to study the hydrodynamic evolution. The hot-electron source was characterized in terms of Maxwellian temperature, T h , and laser to hot-electron energy conversion efficiency η using data from different x-ray spectrometers. The post-processing of these data gives a range of the possible values for T h and η [i.e., T h [keV] ϵ (20, 50) and η ϵ (2%, 13%)]. These values are used as input in hydrodynamic simulations to reproduce the results obtained in radiographs, thus constraining the range for the HE measurements. According to this procedure, we found that the laser converts ~10% ± 4% of energy into hot electrons with T h = 27 ± 8 keV. The paper shows how the coupling of different diagnostics and numerical tools is required to sufficiently constrain the problem, solving the large ambiguity coming from the post-processing of spectrometers data. The effect of the hot electrons on the shock dynamics is then discussed, showing an increase in the pressure around the shock front. Furthermore, the low temperature found in this experiment without pre-compression laser pulses could be advantageous for the SI scheme, but the high conversion efficiency may lead to an increase in the shell adiabat, with detrimental effects on the implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyα Forest

The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyα forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit–accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible IGM thermal histories by varying the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen at redshifts 2.2 ≤ z ≤ 5.0 and with the helium Lyα opacity at redshifts 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyα forest data to break model degeneracies, we find that the IGM experienced two main reheating events over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z R ≈ 6.0 with a first IGM temperature peak of T 0 ≃ 1.3 × 10 4 K and is followed by the reionization of He ii that completes by z R ≈ 3.0 and yields a second temperature peak of T 0 ≃ 1.4 × 10 4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

A platform for non-local thermodynamic equilibrium atomic physics research using the buried layer target approach at a kJ level laser facility

Here we present a design and first use of a kJ level laser facility for research of non-local thermodynamic equilibrium atomic physics using the buried layer target method. The target design included a metal layer buried inside a plastic tamper with thicknesses tailored to the expected laser intensities. The target was illuminated from each side by two laser beams with intensities of 0.5–5 × 10 14 W/cm 2 . The advanced diagnostic suite included static and time-resolved imagers and spectrometers with various spectral resolutions. A 3D printed dual elliptically curved spectrometer is presented, and its results are compared to a traditional crystal spectrometer. Experimental results and radiation hydrodynamic simulations demonstrate that the target achieved the desired thermodynamic conditions of n e ≈ 10 21 –10 22 cm -3 and T e ≈ 1–2 keV.

3D printing↗