Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiation Hydrodynamics Modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Radiation hydrodynamics modeling of the highest compression inertial confinement fusion ignition experiment from the National Ignition Campaign

The recently completed National Ignition Campaign (NIC) on the National Ignition Facility (NIF) showed significant discrepancies between post-shot simulations of implosion performance and experimentally measured performance, particularly in thermonuclear yield. This discrepancy between simulation and observation persisted despite concerted efforts to include all of the known sources of performance degradation within a reasonable two-dimensional (2-D), and even three-dimensional (3-D), simulation model, e.g., using measured surface imperfections and radiation drives adjusted to reproduce observed implosion trajectories [Clark et al., Phys. Plasmas 20, 056318 (2013)]. Since the completion of the NIC, several effects have been identified that could explain these discrepancies and that were omitted in previous simulations. In particular, there is now clear evidence for larger than anticipated long-wavelength radiation drive asymmetries and a larger than expected perturbation seeded by the capsule support tent. This paper describes an updated suite of one-dimensional (1-D), 2-D, and 3-D simulations that include the current best understanding of these effects identified since the NIC, as applied to a specific NIC shot. The relative importance of each effect on the experimental observables is compared. In combination, these effects reduce the simulated-to-measured yield ratio from 125:1 in 1-D to 1.5:1 in 3-D, as compared to 15:1 in the best 2-D simulations published previously. While the agreement with the experimental data remains imperfect, the comparison to the data is significantly improved and suggests that the largest sources for the previous discrepancies between simulation and experiment are now being included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

RADIATION HYDRODYNAMICS MODELS OF THE INNER RIM IN PROTOPLANETARY DISKS

Many stars host planets orbiting within a few astronomical units (AU). The occurrence rate and distributions of masses and orbits vary greatly with the host star’s mass. These close planets’ origins are a mystery that motivates investigating protoplanetary disks’ central regions. A key factor governing the conditions near the star is the silicate sublimation front, which largely determines where the starlight is absorbed, and which is often called the inner rim. We present the first radiation hydrodynamical modeling of the sublimation front in the disks around the young intermediate-mass stars called Herbig Ae stars. The models are axisymmetric and include starlight heating; silicate grains sublimating and condensing to equilibrium at the local, time-dependent temperature and density; and accretion stresses parameterizing the results of MHD magnetorotational turbulence models. The results compare well with radiation hydrostatic solutions and prove to be dynamically stable. Passing the model disks into Monte Carlo radiative transfer calculations, we show that the models satisfy observational constraints on the inner rim’s location. A small optically thin halo of hot dust naturally arises between the inner rim and the star. The inner rim has a substantial radial extent, corresponding to several disk scale heights. While the front’s overall position varies with the stellar luminosity, its radial extent depends on the mass accretion rate. A pressure maximum develops near the location of thermal ionization at temperatures of about 1000 K. The pressure maximum is capable of halting solid pebbles’ radial drift and concentrating them in a zone where temperatures are sufficiently high for annealing to form crystalline silicates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Radiation hydrodynamics modelling of kilonovae with SNEC

ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

DATA-DRIVEN RADIATIVE HYDRODYNAMIC MODELING OF THE 2014 MARCH 29 X1.0 SOLAR FLARE

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph ( IRIS ) in Mg ii h and k, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DST/IBIS) in H α 6563 Å and Ca ii 8542 Å, and the Reuven Ramaty High Energy Solar Spectroscope Imager ( RHESSI ) in hard X-rays. We constructed a multithreaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated Mg ii line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km s{sup −1}) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of ≈800 km below the transition region can turn the simulated Mg ii line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A New Two-fluid Radiation-hydrodynamical Model for X-Ray Pulsar Accretion Columns

Previous research centered on the hydrodynamics in X-ray pulsar accretion columns has largely focused on the single-fluid model, in which the super-Eddington luminosity inside the column decelerates the flow to rest at the stellar surface. This type of model has been relatively successful in describing the overall properties of the accretion flows, but it does not account for the possible dynamical effect of the gas pressure. On the other hand, the most successful radiative transport models for pulsars generally do not include a rigorous treatment of the dynamical structure of the column, instead assuming an ad hoc velocity profile. In this paper, we explore the structure of X-ray pulsar accretion columns using a new, self-consistent, “two-fluid” model, which incorporates the dynamical effect of the gas and radiation pressures, the dipole variation of the magnetic field, the thermodynamic effect of all of the relevant coupling and cooling processes, and a rigorous set of physical boundary conditions. The model has six free parameters, which we vary in order to approximately fit the phase-averaged spectra in Her X-1, Cen X-3, and LMC X-4. In this paper, we focus on the dynamical results, which shed new light on the surface magnetic field strength, the inclination of the magnetic field axis relative to the rotation axis, the relative importance of gas and radiation pressures, and the radial variation of the ion, electron, and inverse-Compton temperatures. The results obtained for the X-ray spectra are presented in a separate paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Integrated modeling of cryogenic layered highfoot experiments at the NIF

Integrated radiation hydrodynamic modeling in two dimensions, including the hohlraum and capsule, of layered cryogenic HighFoot Deuterium-Tritium (DT) implosions on the NIF successfully predicts important data trends. The model consists of a semi-empirical fit to low mode asymmetries and radiation drive multipliers to match shock trajectories, one dimensional inflight radiography, and time of peak neutron production. Application of the model across the HighFoot shot series, over a range of powers, laser energies, laser wavelengths, and target thicknesses predicts the neutron yield to within a factor of two for most shots. The Deuterium-Deuterium ion temperatures and the DT down scattered ratios, ratio of (10–12)/(13–15) MeV neutrons, roughly agree with data at peak fuel velocities <340 km/s and deviate at higher peak velocities, potentially due to flows and neutron scattering differences stemming from 3D or capsule support tent effects. These calculations show a significant amount alpha heating, 1–2.5× for shots where the experimental yield is within a factor of two, which has been achieved by increasing the fuel kinetic energy. This level of alpha heating is consistent with a dynamic hot spot model that is matched to experimental data and as determined from scaling of the yield with peak fuel velocity. Lastly, these calculations also show that low mode asymmetries become more important as the fuel velocity is increased, and that improving these low mode asymmetries can result in an increase in the yield by a factor of several.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Radiation transport and energetics of laser-driven half-hohlraums at the National Ignition Facility

Experiments that characterize and develop a high energy-density half-hohlraum platform for use in benchmarking radiation hydrodynamics models have been conducted at the National Ignition Facility (NIF). Results from the experiments are used to quantitatively compare with simulations of the radiation transported through an evolving plasma density structure, colloquially known as an N-wave. A half-hohlraum is heated by 80 NIF beams to a temperature of 240 eV. This creates a subsonic diffusive Marshak wave, which propagates into a high atomic number Ta{sub 2}O{sub 5} aerogel. The subsequent radiation transport through the aerogel and through slots cut into the aerogel layer is investigated. We describe a set of experiments that test the hohlraum performance and report on a range of x-ray measurements that absolutely quantify the energetics and radiation partition inside the target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Investigation of Relationship between High-energy X-Ray Sources and Photospheric and Helioseismic Impacts of X1.8 Solar Flare of 2012 October 23

The X-class solar flare of 2012 October 23 generated continuum photospheric emission and a strong helioseismic wave (“sunquake”) that points to an intensive energy release in the dense part of the solar atmosphere. We study properties of the energy release with high temporal and spatial resolutions, using photospheric data from the Helioseismic Magnetic Imager (HMI) on board Solar Dynamics Observatory , and hard X-ray observations made by RHESSI . For this analysis we use level-1 HMI data (filtergrams), obtained by scanning the Fe i line (6731 Å) with the time cadence of ∼3.6 s and spatial resolution of ∼0.″5 per pixel. It is found that the photospheric disturbances caused by the flare spatially coincide with the region of hard X-ray emission but are delayed by ≲4 s. This delay is consistent with predictions of the flare hydrodynamics RADYN models. However, the models fail to explain the magnitude of variations observed by the HMI. The data indicate that the photospheric impact and helioseismic wave might be caused by the electron energy flux, which is substantially higher than that in the current flare radiative hydrodynamic models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Modeling Hydrodynamic Instabilities and Mix in National Ignition Facility Hohlraums

Radiation hydrodynamics modeling of hohlraum experiments on the National Ignition Facility has shown a chronic overprediction of hohlraum performance compared to experiment. Depending on the particular experiment, up to several hundred kilojoules of laser energy must be artificially subtracted from the simulations in order to match experimentally measured quantities like the time of peak neutron production (bang time). Fine-scale hydrodynamic mixing of the high-Z hohlraum wall with the low-Z hohlraum gas, an effect omitted in current hohlraum simulations, has long been hypothesized as possible source of some of this degradation. A preliminary computational assessment using a spherical hohlraum model has been made. It shows that the wall-gas interface is both unstable during certain phases of the experiment and has the potential for turbulent mixing of the gas and wall material, depending on the hohlraum gas fill and pulse shape. The same simulations also indicate that mixing at this interface can delay the simulated bang time by up to 140 ps with high gas fills but appears negligible for low gas fills.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

National Opacity Program: Analysis of Opacity-relevant X-ray Emission Spectra (Milestone Report ID 7121)

We have produced high energy density iron plasmas at temperatures above 1 keV and electron densities exceeding 1023 cm -3 (~ 1 g/cm 3 ) using the Orion laser at the Atomic Weapons Establishment. These plasmas were created by irradiating 50 µm diameter layered targets with frequency doubled (λ = 527 nm), 1 ps laser pulses focused to a 100 µm diameter producing an irradiance of ~ 2 x 10 18 W/cm 2 . The buried layer targets consist of 160 nm iron sulfide (FeS), 60 nm potassium chloride (KCl), and 15 nm carbon. The combined layers are tamped on both sides with 3 µm of parylene-N. The x-ray emission from the plasma was measured using two time-resolved, and four time-integrated Bragg crystal spectrometers, as well as one time-integrated imaging system. One time-resolved x-ray spectrometer measured emission from L-shell transitions in highly charged iron, the other from K-shell transitions in helium-like S 14+ and hydrogen-like S 15+ . The time-integrated spectrometers are intensity-calibrated and measured emission from K-shell transitions in sulfur, potassium, chlorine, and both K-shell and L-shell transitions in iron. The density and temperature of the plasma were determined by modeling the x-ray spectra using different spectral and hydrodynamic modeling packages. A brief overview of the uncertainties associated with the measurements and models are presented. We also give an overview of our 1-D HYDRA-DCA radiation-hydrodynamics model and improvements for future work. Our results aid in assessing experimental uncertainties associated with plasma uniformity and with x-ray emission employed as diagnostics in opacity experiments at temperatures and densities not achievable elsewhere and represent a significant step in creating and diagnosing plasmas near LTE. These results are summarized as part of the completion requirements for milestone 7121

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

X-RAY OBSERVATIONS OF BOW SHOCKS AROUND RUNAWAY O STARS. THE CASE OF ζ OPH AND BD+43°3654

Non-thermal radiation has been predicted within bow shocks around runaway stars by recent theoretical works. We present X-ray observations toward the runaway stars ζ Oph by Chandra and Suzaku and of BD+43°3654 by XMM-Newton to search for the presence of non-thermal X-ray emission. We found no evidence of non-thermal emission spatially coincident with the bow shocks; nonetheless, diffuse emission was detected in the vicinity of ζ Oph. After a careful analysis of its spectral characteristics, we conclude that this emission has a thermal nature with a plasma temperature of T ≈ 2 × 10{sup 6} K. The cometary shape of this emission seems to be in line with recent predictions of radiation-hydrodynamic models of runaway stars. The case of BD+43°3654 is puzzling, as non-thermal emission has been reported in a previous work for this source.

79 ASTRONOMY AND ASTROPHYSICS↗