Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

2D Magnetohydrodynamic Simulations of the Electrothermal Instability in Metallic Liners

The Virginia Tech (VT) Plasma Dynamics Laboratory Computational (PDCL) and Lawrence Livermore National Laboratory (LLNL) are performing two dimensional (2D) simulations of the electrothermal instability (ETI) using the LLNL multi-physics code Ares. Ares is a multi-physics arbitrary–Lagrangian-Eulerian (ALE) code developed by LLNL and is of particular use in studying magnetohydrodynamic (MHD) instabilities like the ETI due to its resistive MHD, magnetic diffusion, and radiative-hydrodynamics packages. Among its capabilities, it has the ability to model material strength, perform adaptive mesh refinement (AMR), and incorporate a wide variety of equations of state models and conductivity models. The 2D Ares simulation model created by VT-LLNL for studying the development and growth of the electrothermal instability has been configured with initial conditions based on the Mykonos Electrothermal Instability II (METI-II) experiments described by this grant and conducted by team members at the University of Nevada (UNR), the University of New Mexico (UNM), and Sandia National Laboratories. Previously, preliminary 2D Ares simulations of the ETI had been run to approximately 80ns. The rods in these simulations were initiated with sinusoidal perturbations at a similar order of magnitude to those measured on the aluminum rods used for the Mykonos experiment. This model has been improved by increasing the spatial resolution of the simulations and running the simulations further in time. In addition to the simulation run-times extending, the preliminary sinusoidal perturbation has been replaced with a perturbation derived from amplitude measurements by the experimental team, thereby correlating the simulation inputs better to the experimental runs. These new runs are capable of reaching 120ns of simulated time for the uncoated cases and to 200ns the 41 μm coated cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Preshot Report for OMEGA Double Cylinders (CylDRT 22B)

The primary goals for the May 2022 shot day at OMEGA are to: image an imploding inner cylinder as proof-of-principle for double cylinder experiments, measure the growth of pre-seeded perturbations on the inner cylinder, and measure the axial non-uniformity of the implosion with a sidelighter that provides a transverse view of the target. We will compare these experimental results to xRAGE (2D-only) and FLASH (2D and 3D) radiation-hydrodynamics calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulation of Direct-Drive Hybrid Using Two Opposed Beams for Inertial Fusion Energy

Xcimer Energy is working to deliver high levels of laser light at the costs required for applications in fusion energy. The architecture of their integrated laser system is unique, in that it naturally supports a target which is illuminated by two opposed beams, and in a manner that can readily complement a reactor. This INFUSE proposal was submitted to investigate a target concept termed the Direct-Drive Hybrid or DDH, and make special use of radiation hydrodynamic simulations and expertise available to the Laboratory for Laser Energetics. The associated capabilities are unique to the field, and take great advantage of investments and advances made by the DOE and DOD over many decades. The primary goals were to stand up the DDH design in calculations for the first time, study and refine the concept, and provide the scientific basis for understanding and projecting proposals by Xcimer Energy. The main accomplishments have been summarized below.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Viscosity Measurements in Extreme Conditions

Radiation-hydrodynamics simulations are critical to the NNSA complex for understanding high energy density physics and for addressing problems in national security science. However, current rad-hydro codes typically do not account for viscosity, which generally leads to large uncertainties in problems involving mixing of materials under shock loading. This report details progress made towards experimental measurements of fluid viscosity at high pressures and temperatures achieved by shock wave compression techniques, as well as complementary material model development and hydrodynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dataset from 3D FLASH Computation of National Ignition Facility Shot N210201-001

This data comes from a 3D FLASH [Fryxell et al. ApJS 131, 273 (2000); Tzeferacos et al. HEDP 17, 24-31(2015)] computation of a cylindrical implosion experiment performed at the National Ignition Facility (shot number N210201-001). FLASH is an open-access radiation-hydrodynamics code being used to simulate high-energy-density systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Determining the Structure of Rotating Massive Stellar Cores with Gravitational Waves

The gravitational wave (GW) signal resulting from stellar core collapse encodes a wealth of information about the physical parameters of the progenitor star and the resulting core-collapse supernova (CCSN). We present a novel approach to constrain CCSN progenitor properties at collapse using two of the most detectable parts of the GW signal: the core-bounce signal and evolution of the dominant frequency mode from the protoneutron star. We focus on the period after core bounce but before explosion and investigate the predictive power of GWs from rotating CCSNe to constrain properties of the progenitor star. We analyze 34 2D and four 3D neutrino-radiation-hydrodynamic simulations of stellar core collapse in progenitors of varied initial mass and rotation rate. Extending previous work, we verify the compactness of the progenitor at collapse to correlate with the early ramp-up slope, and in rotating cases, also with the core angular momentum. Combining this information with the bounce signal, we present a new analysis method to constrain the pre-collapse core compactness of the progenitor. Because these GW features occur less than a second after core bounce, this analysis could allow astronomers to predict electromagnetic properties of a resulting CCSN even before shock breakout.

Astronomy & Astrophysics↗
DOE OSTI2021

Three-dimensional Core-collapse Supernova Simulations with 160 Isotopic Species Evolved to Shock Breakout

We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino radiation hydrodynamic simulations of the neutrino-driven phase performed with the Chimera code. We consider a 9.6 M ⊙ zero-metallicity progenitor starting from both 2D and 3D Chimera models and a 10 M ⊙ solar-metallicity progenitor starting from a 2D Chimera model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh–Taylor (R-T) "bullets" determines how the metal-rich ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum 56 Ni velocities of ~1950 and ~1750 km s -1 at shock breakout from 2D and 3D 9.6 M ⊙ Chimera models, respectively, due to the bullets' ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher-velocity structures is suppressed when the 2D Chimera model and 3D FLASH model meshes are aligned. The development of faster-growing spherical-bubble structures, as opposed to the slower-growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10 M ⊙ model, which achieves maximum 56 Ni velocities of ~2500 km s -1 at shock breakout.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Characterizing the Directionality of Gravitational Wave Emission from Matter Motions within Core-collapse Supernovae

Abstract We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T /∣ W ∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

Astronomy & Astrophysics↗
DOE OSTI2024

How Population III Supernovae Determined the Properties of the First Galaxies

Abstract Massive Population III stars can die as energetic supernovae that enrich the early Universe with metals and determine the properties of the first galaxies. With masses of about 10 9 M ⊙ atz≳ 10, these galaxies are believed to be the ancestors of the Milky Way. This paper investigates the impact of Population III supernova remnants (SNRs) from both Salpeter-like and top-heavy initial mass functions (IMFs) on the formation of first galaxies with high-resolution radiation-hydrodynamical simulations with theENZOcode. Our findings indicate that SNRs from a top-heavy Population III IMF produce more metals, leading to more efficient gas cooling and earlier Population II star formation in the first galaxies. From a few hundred to a few thousand Population II stars can form in the central regions of these galaxies. These stars have metallicities of 10 −3 –10 −2 ,Z ⊙ , greater than those of extremely metal-poor (EMP) stars. Their mass function follows a power-law distribution with dN ( M * ) / dM * ∝ M * α , whereM * is stellar mass, andα= 2.66–5.83 and is steeper for a top-heavy IMF. We thus find that EMP stars were not typical of most primitive galaxies.

Astronomy & Astrophysics↗
DOE OSTI2025

The BTSbot-nearby Discovery of SN 2024jlf: Rapid, Autonomous Follow-up Probes Interaction in an 18.5 Mpc Type IIP Supernova

We present observations of the Type IIP supernova (SN) SN 2024jlf, including spectroscopy beginning just 0.7 days (∼17 hr) after first light. Rapid follow-up was enabled by the new BTSbot-nearby program, which involves autonomously triggering target-of-opportunity requests for new transients in Zwicky Transient Facility data that are coincident with nearby (D < 60 Mpc) galaxies and identified by the BTSbot machine learning model. Early photometry and nondetections shortly prior to first light show that SN 2024jlf initially brightened by >4 mag day −1 , quicker than ∼90% of Type II SNe. Early spectra reveal weak flash ionization features: narrow, short-lived (1.3 < τ[days] < 1.8) emission lines of Hα, He II , and C IV . Assuming a wind velocity of v w = 50 km s −1 , these properties indicate that the red supergiant progenitor exhibited enhanced mass loss in the last year before explosion. We constrain the mass-loss rate to $1{0}^{-4}\lt \dot{M}\,[{M}_{\odot }\,{\mathrm{yr}}^{-1}]\lt 1{0}^{-3}$ by matching observations to model grids from two independent radiative hydrodynamics codes. BTSbot-nearby automation minimizes spectroscopic follow-up latency, enabling the observation of ephemeral early-time phenomena exhibited by transients.

core-collapse supernovae↗
DOE OSTI2025

Measurement and Modeling of Electron Temperature in Laboratory Photoionized Plasmas Relevant to Astrophysics

The electron temperature of photoionized plasmas characterizes the thermalization of photoelectrons, impacts the charge-state distribution, emissivity, and opacity through atomic recombination processes, and is needed to perform detailed comparisons with theory predictions. We discuss temperature measurements in laboratory photoionized plasmas and a comparison with model calculations done with several theory approximations and codes. These include a radiation-hydrodynamics simulation and two nonequilibrium heating and ionization models that tracked the evolution of the internal energy of the electrons. Furthermore, for the same physics model and X-ray flux time history, calculations were performed assuming steady-state or time-dependent conditions. The time history of steady-state results correlates with that of the X-ray flux, while that of the time-dependent cases does not, and it is qualitatively and quantitatively different from the steady-state case. Steady-state results significantly overestimated temperature measurements, while time-dependent results produced better approximations, which suggests the importance of transient effects in the experiment and also the need for time-resolved measurements.

Laboratory astrophysics↗
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
DOE OSTI2022

Design of inertial fusion implosions reaching the burning plasma regime

One of the last remaining milestones in fusion research before reaching ignition is creating a burning plasma state, where alpha particles from deuterium-tritium (DT) fusion reactions redeposit their energy as the dominant source of heating in the plasma. The indirect-drive inertial confinement fusion approach at the National Ignition Facility (NIF) uses a laser-generated radiation cavity (hohlraum) to spherically implode DT fuel to high temperatures and densities in a central ”hot spot”. Here, we deliver more energy to the hot spot than ever before, while maintaining the extreme pressures required for inertial confinement, by increasing the size of the implosion compared to previous experiments. We develop more efficient hohlraums, to drive these larger implosions within NIF’s current laser energy and power capability and control symmetry by moving energy between laser beams and by changing the shape of the hohlraum. These designs resulted in record fusion powers of 1.5 petawatts, greater than the input power of the laser, and 170 kJ of fusion energy. Radiation hydrodynamics simulations show alpha particle heating as the dominant term in the hot spot energy balance, e.g. a burning plasma state. This work is expected to motivate future studies of burning plasmas and improve predictive capability by providing a benchmark for modeling used to understand the proximity to ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

On the mass of the Cepheid X Cygni

Using a new non-Lagrangian hydrodynamic radiation transport code, models of the Cepheid X Cygni have been calculated and compared with new observations. Adopting an effective temperature of 5300 K, the luminosity was varied for stellar masses of 6, 8, and 9 solar masses in order to obtain a period of 16.4d. The results indicate that in order to find reasonable agreement between the light curves as calculated and observed, the mass used in the models must be nearly the evolutionary mass, not some 60% less as found in earlier work for Cepheids with periods from 7d to 10d. In addition, the new observations support the physical reality of the so-called 'artificial viscosity dip'.

Davis, C. G.↗
NASA NTRS1986

Grain processes in massive star formation

Observational evidence suggests that stars greater than 100 M(solar) exist in the Galaxy and Large Magellanic Cloud (LMC), however classical star formation theory predicts stellar mass limits of only approx. 60 M(solar). A protostellar accretion flow consists of inflowing gas and dust. Grains are destroyed as they are near the central protostar creating a dust shell or cocoon. Radiation pressure acting on the grain can halt the inflow of material thereby limiting the amount of mass accumulated by the protostar. We first consider rather general constraints on the initial grain to gas ratio and mass accretion rates that permit inflow. We further constrain these results by constructing a numerical model. Radiative deceleration of grains and grain destruction processes are explicitly accounted for in an iterative solution of the radiation-hydrodynamic equations. Findings seem to suggest that star formation by spherical accretion requires rather extreme preconditioning of the grain and gas environment.

Wolfire, M. G.↗
NASA NTRS1992

A fast operator perturbation method for the solution of the special relativistic equation of radiative transfer in spherical symmetry

A fast method for the solution of the radiative transfer equation in rapidly moving spherical media, based on an approximate Lambda-operator iteration, is described. The method uses the short characteristic method and a tridiagonal approximate Lambda-operator to achieve fast convergence. The convergence properties and the CPU time requirements of the method are discussed for the test problem of a two-level atom with background continuum absorption and Thomson scattering. Details of the actual implementation for fast vector and parallel computers are given. The method is accurate and fast enough to be incorporated in radiation-hydrodynamic calculations.

Hauschildt, P. H.↗
NASA NTRS1992

Non-LTE radiating acoustic shocks and Ca II K2V bright points

We present, for the first time, a self-consistent solution of the time-dependent 1D equations of non-LTE radiation hydrodynamics in solar chromospheric conditions. The vertical propagation of sinusoidal acoustic waves with periods of 30, 180, and 300 s is calculated. We find that departures from LTE and ionization recombination determine the temperature profiles of the shocks that develop. In LTE almost all the thermal energy goes into ionization, so the temperature rise is very small. In non-LTE, the finite transition rates delay the ionization to behind the shock front. The compression thus goes into thermal energy at the shock front leading to a high temperature amplitude. Further behind the shock front, the delayed ionization removes energy from the thermal pool, which reduces the temperature, producing a temperature spike. The 180 s waves reproduce the observed temporal changes in the calcium K line profiles quite well. The observed wing brightening pattern, the violet/red peak asymmetry and the observed line center behavior are all well reproduced. The short-period waves and the 5 minute period waves fail especially in reproducing the observed behavior of the wings.

Carlsson, Mats↗
NASA NTRS1992

An interpretation of the multipeaked structure in X-ray bursts

Some X-ray bursts from sources such as 4U/MXB 1636 - 53, 4U 1608 - 52, and possibly MXB 1727 - 34, display a multipeaked structure in the bolometric flux that is difficult to reconcile with the standard thermonuclear flash model or with the effects due to photospheric expansion and contraction. Understanding this phenomenon is important primarily because it might be due to a variable Thomson optical depth along a line of sight threading the burst-induced accretion disk corona. We explore the conditions under which these burst profiles might be observed, in light of the somewhat surprising recent results from the radiative-hydrodynamic simulation of the coronal structure. We find that double- (and triple-) peaked bursts are produced at inclination angles in the range 40-65 deg. Interestingly, the orbital information for 4U/MXB 1636 - 53 and 4U 1608 - 52 suggests i is less than 70 deg, in these sources. The orbital inclination of MXB 1728 - 34 is not yet known.

Melia, Fulvio↗