Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiation Hydrodynamics”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Possible octupole deformation of 208 Pb and the ultracentral v 2 to v 3 puzzle

Recent measurements have established the sensitivity of ultracentral heavy-ion collisions to the deformation parameters of nonspherical nuclei. In the case of 129 Xe collisions, a quadrupole deformation of the nuclear profile led to an enhancement of elliptic flow in the most central collisions. In 208 Pb collisions a discrepancy exists in similar centralities, where either elliptic flow is overpredicted or triangular flow is underpredicted by hydrodynamic models; this is known as the v 2 -to-v 3 puzzle in ultracentral collisions. Motivated by low-energy nuclear structure calculations, we consider the possibility that 208Pb nuclei could have a pear-shape deformation (octupole), which has the effect of increasing triangular flow in central PbPb collisions. Using the recent data from ALICE and ATLAS, we reexamine the v 2 -to-v 3 puzzle in ultracentral collisions, including new constraints from recent measurements of the triangular cumulant ratio v 3 {4}/v 3 {2} and comparing two different hydrodynamic models. We find that while an octupole deformation would slightly improve the ratio between v 2 and v 3 , it is at the expense of a significantly worse triangular flow cumulant ratio. In fact, the latter observable prefers no octupole deformation, with β 3 ≲ 0.0375 for 208 Pb, and is therefore consistent with the expectation for a doubly-magic nucleus even at top collider energies. Furthermore, the v 2 -to-v 3 puzzle remains a challenge for hydrodynamic models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Observational Signatures of Planets in Protoplanetary Disks: Temperature Structures in Spiral Arms

High-resolution imaging of protoplanetary disks has unveiled a rich diversity of spiral structure, some of which may arise from disk-planet interaction. Using 3D hydrodynamics with β cooling to a vertically stratified background, as well as radiative-transfer modeling, we investigate the temperature rise in planet-driven spirals. In rapidly cooling disks, the temperature rise is dominated by a contribution from stellar irradiation, 0.3%–3% inside the planet radius but always < 0.5% outside. When cooling time equals or exceeds dynamical time, however, this is overwhelmed by hydrodynamic PdV work, which introduces a ~10%–20% perturbation within a factor of ~2 from the planet’s orbital radius. We devise an empirical fit of the spiral amplitude Δ(T)=(M{sub p}/M{sub th}){sup c}(f{sub PdV}e{sup -t{sub a}{sub r}{sub m}/t{sub c}}+f{sub rad}) to take into account both effects. Where cooling is slow, we find also that temperature perturbations from buoyancy spirals—a strictly 3D, nonisothermal phenomenon—become nearly as strong as those from Lindblad spirals, which are amenable to 2D and isothermal studies. Our findings may help explain observed thermal features in disks like TW Hydrae and CQ Tauri, and underscore that 3D effects have a qualitatively important effect on disk structure.

42 ENGINEERING↗
DOE OSTI2022

Kinematic Evidence for an Embedded Planet in the IM Lupi Disk

Abstract We test the hypothesis that an embedded giant planet in the IM Lupi protostellar disk can produce velocity kinks seen in CO line observations as well as the spiral arms seen in scattered light and continuum emission. We inject planets into 3D hydrodynamics simulations of IM Lupi, generating synthetic observations using Monte Carlo radiative transfer. We find that an embedded planet of 2–3 M Jup can reproduce non-Keplerian velocity perturbations, or “kinks”, in the 12 CO J = 2–1 channel maps. Such a planet can also explain the spiral arms seen in 1.25 mm dust continuum emission and 1.6 μ m scattered-light images. We show that the wake of the planet can be traced in the observed peak velocity map, which appears to closely follow the morphology expected from our simulations and from analytic models of planet–disk interaction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Observational signatures of circumbinary discs – I. Kinematics

We present five morphological and kinematic criteria to aid in asserting the binary nature of a protoplanetary disc, based on 3D hydrodynamical simulations of circumbinary discs post-processed with Monte Carlo radiative transfer. We find that circumbinary discs may be identified by (i) a central cavity, (ii) spiral arms both in and outside of their central cavities, (iii) non-localized perturbations in their iso-velocity curves, (iv) asymmetry between the lines of maximum speed of the blue shifted and redshifted wings, and (v) asymmetry between the area of the blue shifted and redshifted wings. In this study, we provide quantitative metrics for the last two criteria that can be used, in conjunction with the morphological criteria, to signal whether a protoplanetary disc is likely to be a circumbinary disc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks

Reacting astrophysical flows can be challenging to model, because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach, based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting, that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double-detonation problem with a moderately sized astrophysical nuclear reaction network and explore the time step size and reaction network tolerances, to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

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

Nonperturbative heavy-flavor transport approach for hot QCD matter

The heavy charm and bottom quarks are unique probes of the transport properties of the quark-gluon plasma (QGP) and its hadronization in high-energy nuclear collisions. A key challenge in this context is to embed the interactions of the heavy quarks in the expanding medium compatible with the strong-coupling nature of the QGP, and thus to unravel the underlying microscopic mechanisms. In the present work we progress toward this goal by combining recent T -matrix interactions for elastic scattering with an effective transport implementation of gluon radiation, and apply these in a Langevin framework in a viscous hydrodynamic evolution. Hadronization of heavy quarks is evaluated using a modern recombination model with 4-momentum conservation, supplemented with fragmentation constrained by data in proton-proton collisions. Deploying this approach to charm-hadron observables in Pb-Pb collisions at the LHC yields fair agreement with experiment while also identifying areas of further systematic improvement of the simulations and its current input.

Heavy-flavor transport↗
DOE OSTI2020

Characterization of an imploding cylindrical plasma for electron transport studies using x-ray emission spectroscopy

We report here on the characterization of the conditions of an imploding cylindrical plasma by time-resolved x-ray emission spectroscopy. Knowledge about this implosion platform can be applied to studies of particle transport for inertial confinement fusion schemes or to astrophysical plasmas. A cylindrical Cl-doped CH foam within a tube of solid CH was irradiated by 36 beams (I total ~5×10 14 W/cm 2 , 1.5ns square pulse, and E total ~ 16.2kJ) of the OMEGA-60 laser to radially compress the CH toward the axis. The analysis of the time-resolved spectra showed that the compression can be described by four distinct phases, each presenting different plasma conditions. First the ablation of the cylinder is dominant; second, the foam is heated and induces a significant jump in emission intensities; third, the temperature and density of the foam reaches a maximum; and finally, the plasma expands. Ranges for the plasma temperature were inferred with the atomic physics code SCRAM (Spectroscopic Collisional-Radiative Atomic Model) and the experimental data have been compared to hydrodynamic simulations performed with the 2D code FLASH, which showed a similar implosion dynamic over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Revealing the Galaxy–Halo Connection through Machine Learning

Abstract Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. The complexity of these physics can prove difficult to disentangle and obfuscate how mass-dependent trends in the galaxy population originate. Here, we train a machine-learning method called Explainable Boosting Machines (EBMs) to infer how the stellar mass and star formation rate of nearly 6 million galaxies simulated by the Cosmic Reionization on Computers project depend on the physical properties of halo mass, the peak circular velocity of the galaxy during its formation history v peak , cosmic environment, and redshift. The resulting EBM models reveal the relative importance of these properties in setting galaxy stellar mass and star formation rate, with v peak providing the most dominant contribution. Environmental properties provide substantial improvements for modeling the stellar mass and star formation rate in only ≲10% of the simulated galaxies. We also provide alternative formulations of EBM models that enable low-resolution simulations, which cannot track the interior structure of dark matter halos, to predict the stellar mass and star formation rate of galaxies computed by high-resolution simulations with detailed baryonic physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A line-binned treatment of opacities for the spectra and light curves from neutron star mergers

ABSTRACT The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy r-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound–bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion–opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and thus eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion–opacity formalism. The agreement between the line-binned and expansion–opacity results is explained as arising from the similarity in their opacities in the limit of low optical depth, where radiation transport is important in the ejecta. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semirelativistic approaches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

From Hydrodynamics to Jet Quenching, Coalescence, and Hadron Cascade: A Coupled Approach to Solving the $R_{AA}Ⓧv_{2}$ Puzzle

Hydrodynamics and jet quenching are responsible for the elliptic flow $v_{2}$ and suppression of large transverse momentum ($p_{T}$) hadrons, respectively, two of the most important phenomena leading to the discovery of a strongly coupled quark-gluon plasma in high-energy heavy-ion collisions. A consistent description of the hadron suppression factor $R_{AA}$ and $v_{2}$, especially at intermediate $p_{T}$, however, remains a challenge. We solve this long-standing $R_{AA}Ⓧv_{2}$ puzzle by including quark coalescence for hadronization and final state hadron cascade in the coupled linear Boltzmann transport-hydro model that combines concurrent jet transport and hydrodynamic evolution of the bulk medium. We illustrate that quark coalescence and hadron cascade, two keys to solving the puzzle, also lead to a splitting of $v_{2}$ for pions, kaons, and protons in the intermediate $p_{T}$ region. We demonstrate for the first time that experimental data on $R_{AA}, v_{2}$, and their hadron flavor dependence from low to intermediate and high $p_{T}$ in high-energy heavy-ion collisions can be understood within this coupled framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Feasibility of measuring the speed of sound of the quark-gluon plasma from the multiplicity and mean 𝑝 𝑇 of ultracentral heavy-ion collisions

The mean transverse momentum ⟨𝑝 𝑇 ⟩ of hadrons has been observed experimentally and in numerical simulations to have a power-law dependence on the hadronic multiplicity 𝑁 in ultracentral relativistic heavy-ion collisions: ⟨𝑝 𝑇 ⟩∝𝑁 𝑏 UC . It has been put forward that this exponent 𝑏 UC is the speed of sound of quark-gluon plasma measured at a temperature determined from ⟨𝑝 𝑇 ⟩. We study step by step the connection between (i) the energy and entropy of hydrodynamic simulations and (ii) experimentally measurable observables. We show that an argument based on energy and entropy should yield an exponent equal to the pressure over energy density 𝑃/ɛ, rather than the speed of sound 𝑐$_s^2$; however, we also observe that ⟨𝑝 𝑇 ⟩ and 𝑁 are not sufficiently accurate proxies for the energy and entropy to make this possible in practice. From simulations, we find that the exponent 𝑏 UC is significantly different whether the “effective volume” is strictly constant or not, a condition that cannot be enforced experimentally. Additional tests using a modified equation of state find that the exponent 𝑏 UC exhibits a variable degree of correlations with the speed of sound and with 𝑃/ɛ, but is not an accurate measurement of either quantity in general.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Mixing in ICF implosions on the National Ignition Facility caused by the fill-tube

The micrometer-scale tube that fills capsules with thermonuclear fuel in inertial confinement fusion experiments at the National Ignition Facility is also one of the implosion's main degradation sources. It seeds a perturbation that injects the ablator material into the center, radiating away some of the hot-spot energy. This paper discusses how the perturbation arises in experiments using high-density carbon ablators and how the ablator mix interacts once it enters the hot-spot. Both modeling and experiments show an in-flight areal-density perturbation and localized x-ray emission at stagnation from the fill-tube. Simulations suggest that the fill-tube is degrading an otherwise 1D implosion by ~2×, but when other degradation sources are present, the yield reduction is closer to 20%. Characteristics of the fill-tube assembly, such as the through-hole size and the glue mass, alter the dynamics and magnitude of the degradation. Overall, these aspects point the way toward improvements in the design, some of which (smaller diameter fill-tube) have already shown improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Mapping the Planetary Wake in HD 163296 with Kinematics

Abstract We map the planetary wake associated with the embedded protoplanet creating the CO kink in the disk of HD 163296. We show that the wake can be traced by a series of correlated perturbations in the peak velocity map. The sign change of the perturbations across the disk’s major axis confirms that the wake induces predominantly radial motion, as predicted by models of planet–disk interaction. These results provide the first direct confirmation of planet wakes generated by Lindblad resonances. Mapping the wake provides a constraint on the disk aspect ratio, which is required to measure the mass of the planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗