Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Burn”

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 955 records · Page 53

Violent mergers revisited: The origin of the fastest stars in the Galaxy

Binary systems of two carbon-oxygen white dwarfs are one of the most promising candidates for the progenitor systems of Type Ia supernovae. Violent mergers, where the primary white dwarf ignites when the secondary white dwarf smashes into it while being disrupted on its last orbit, were the first double degenerate merger scenario proposed that ignites dynamically. However, violent mergers likely contribute only a few percent to the total Type Ia supernova rate and do not yield normal Type Ia supernova light curves. Here we revisit the scenario, simulating a violent merger with better methods and, in particular, a more accurate treatment of the detonation. We find good agreement with previous simulations but with one critical difference: The secondary white dwarf being disrupted and accelerated towards the primary white dwarf and impacted by its explosion does not fully burn, and its core survives as a bound object. The explosion leaves behind a 0.16 M ⊙ star travelling 2800 km/s, making it an excellent (and so far the only) candidate to explain the origin of the fastest observed hypervelocity stars. We also show that before the explosion, 5 × 10 −3 M ⊙ of material predominantly consisting of helium, carbon, and oxygen had already been ejected at velocities above 1000 km/s. Finally, we argue that if a violent merger made the hypervelocity stars D6-1 and D6-3 and violent mergers require the most massive primary white dwarfs in binaries of two carbon-oxygen white dwarfs, there has to be a much larger population of white dwarf mergers with slightly lower mass primary white dwarfs. Because this population likely represents ≫10% of the Type Ia supernovae rate, it can essentially only give rise to normal Type Ia supernovae.

Astronomy and AstroPhysics↗

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗

Density determination of the thermonuclear fuel region in inertial confinement fusion implosions

Understanding of the thermonuclear burn in an inertial confinement fusion implosion requires knowledge of the local deuterium–tritium (DT) fuel density. Neutron imaging of the core now provides this previously unavailable information. Two types of neutron images are required. The first is an image of the primary 14-MeV neutrons produced by the D +T fusion reaction. The second is an image of the 14-MeV neutrons that leave the implosion hot spot and are downscattered to lower energy by elastic and inelastic collisions in the fuel. These neutrons are measured by gating the detector to record the 6–12 MeV neutrons. Using the reconstructed primary image as a nonuniform source, a set of linear equations is derived that describes the contribution of each voxel of the DT fuel region to a pixel in the downscattered image. Using the measured intensity of the 14-MeV neutrons and downscattered images, the set of equations is solved for the density distribution in the fuel region. The method is validated against test problems and simulations of high-yield implosions. The calculated DT density distribution from one experiment is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantitative metrics for evaluating thermonuclear design codes and physics models applied to the National Ignition Campaign

In this study, physics models and design codes for hot dense plasmas undergoing thermonuclear burn are evaluated objectively using statistical metrics that compare the difference between calculations and data relative to the experimental uncertainties. The analysis is applied to the National Ignition Campaign (NIC) because it is relevant, comprehensive, and well documented. The statistics confirm that a key process afflicting NIC performance is mix driven by hydrodynamic instabilities as approximated here using the KL model [G. Dimonte and R. Tipton, Phys. Fluids 18, 085101 (2006)]. New physics models are also presented for instability-driven magnetic fields [B. Srinivasan et al., Phys. Rev. Lett. 108, 165002 (2012)] and the Coulomb logarithm for electron–ion thermal relaxation [G. Dimonte and J. Daligault, Phys. Rev. Lett. 101, 135001 (2008)]. The plasma-generated magnetic fields improve the agreement between code and data in a statistically significant manner by reducing the electron thermal conduction in the hot-spot via the Hall term. The Coulomb logarithm is presented in a numerically practical form that incorporates recent theoretical advances. However, it does not improve the agreement between the code and data because the thermal relaxation is so fast in non-ignited NIC plasmas that the typical 20% error does not change the tight coupling between electrons and ions. Even with these improved physics models, the one-dimensional simulations presented here are not able to describe the high-convergence NIC implosions in a statistically acceptable manner. Of the heroic multi-dimensional simulations of Clark et al. [Phys. Plasmas 23, 056302 (2016)], only those in three-dimensions satisfy the statistical acceptance criterion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502

High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Empirical correlations between the function of entropy ( Z S ) and net artificial viscous work in a shock physics hydrocode

Entropy is a state variable that may be obtained from any thermodynamically complete equation of state (EOS). However, hydrocode calculations that output the entropy often contain numerical errors; this is not because of the EOS, but rather the solution techniques that are used in hydrocodes (especially Eulerian) such as convection, remapping, and artificial viscosity. Here, in this work, empirical correlations are investigated to reduce the errors in entropy without altering the solution techniques for the conservation of mass, momentum, and energy. Specifically, these correlations are developed for the function of entropy Z S , and they depend upon the net artificial viscous work, as determined via Sandia National Laboratories’ shock physics hydrocode CTH. These results are a continuation of a prior effort to implement the entropy-based CREST reactive burn model in CTH, and they are presented here to stimulate further interest from the shock physics community. Future work is planned to study higher-dimensional shock waves, shock wave interactions, and possible ties between the empirical correlations and a physical law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hotspot conditions achieved in inertial confinement fusion experiments on the National Ignition Facility

We describe the overall performance of the major indirect-drive inertial confinement fusion campaigns executed at the National Ignition Facility. With respect to the proximity to ignition, we can describe the performance of current experiments both in terms of no-burn ignition metrics (metrics based on the hydrodynamic performance of targets in the absence of alpha-particle heating) and in terms of the thermodynamic properties of the hotspot and dense fuel at stagnation—in particular, the hotspot pressure, temperature, and areal density. We describe a simple 1D isobaric model to derive these quantities from experimental observables and examine where current experiments lie with respect to the conditions required for ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A model for radiative heating of a high-Z pusher

Several approaches to Inertial Confinement Fusion (ICF), including double-shell, pushered-single-shell, and the Revolver designs, have fuel surrounded by pushers made from high-Z materials. An advantage of these designs is that radiation emitted by the hot fuel will be absorbed and re-radiated into the fuel to reduce cooling. This process is referred to as radiation-trapping, and it lowers the fuel temperature required for ignition. To elucidate the physics of radiation trapping, a simplified model comprising a set of coupled differential equations has been developed to model the arbitrary time-dependent wall temperature from the Hammer and Rosen solution [J. H. Hammer and M. D. Rosen, Phys. Plasmas 10, 1829 (2003)] to a Marshak wave. The derivation of the model and a set of analytical power-law solutions will be presented. A discussion on numerical implementation of the differential equations into a model for burn in ICF capsules is also included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The rate of development of atomic mixing and temperature equilibration in inertial confinement fusion implosions

The MARBLE project is a novel inertial confinement fusion platform for studying the development of atomic mixing and temperature equilibration in inertial confinement fusion implosions and their impact on thermonuclear burn. Experiments involve the laser-driven implosion of capsules filled with deuterated engineered foams whose pores are filled with a gaseous mixture of hydrogen and tritium. By varying the size of the foam pores, we can study the timescale of the development of atomic mix relative to the development of thermal equilibrium between species. In contrast, previous separated reactant experiments have only provided information on the total amount of mix mass. Additionally, we report on the series of MARBLE experiments [first reported in Haines et al., Nat. Commun. 11, 544 (2020)] performed on the University of Rochester's OMEGA laser facility and detailed and highly resolved three-dimensional radiation-hydrodynamic simulations of the implosions. In both the experimental and simulation results, we observe that the reactants do not achieve thermal equilibrium during the course of the implosion except in atomically mixed regions—i.e., that atomic mixing develops faster than thermal equilibration between species. The results suggest that ion temperature variations in the mixture are at least as important as reactant concentration variations for determining the fusion reaction rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing laser preheat for laser-driven magnetized liner inertial fusion using soft x-ray emission

Laser heating of a gas in a cylindrical liner held in by a thin foil window is a critical process in the MagLIF (magnetized liner inertial fusion) concept. Window burn-through and gas heating for OMEGA scale MagLIF cylinders as a function of time have been determined using spectrally integrated soft x-ray diagnostics. Window laser absorption is classified in terms of the emitted x-rays from the window plasma as a function of laser energy, and shows that laser energy absorbed is weakly dependent on incident intensity. Radiation-hydrodynamic simulations overestimate the amount of laser energy absorbed by the window as evidenced by the increase in x-ray radiation across several photon energy bands compared to experiments. Gas temperatures inferred from soft x-ray emission from the front 1 mm of the liner are shown to evolve in time in a similar manner to simulation predictions. In conclusion, soft x-ray emission from the gas within the region of the liner that is normally imploded is shown to meet the 100 eV requirements set by the initial point design for laser-driven MagLIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experiments to explore the influence of pulse shaping at the National Ignition Facility

The shaping of the drive pulse in time is a key tool in the design of fusion experiments that use inertia to confine burning plasmas. Here, it is directly related to the adiabat and compressibility of the DT fuel, and the characteristics of the laser and target that are needed to ignite. With this in mind, we have performed experiments at the National Ignition Facility that test small changes in the shape of the pulse. In contrast to theory, we find implosions at lower adiabats can have reduced yield and areal density. We discuss implications to performance and the mechanism(s) that could be responsible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-uniform longitudinal current density induced power saturation in GaAs-based high power diode lasers

The output power of modern 975 nm GaAs-based broad area diode lasers is limited by increasing carrier and photon losses at high bias. Here, we use experiment and one-dimensional calculations on these devices to reveal that higher current densities (and hence higher local recombination rates and higher losses) arise near the front facet due to spatial hole burning and that the non-uniformity is strongly affected by laser geometry, which is more severe for longer resonators and less severe for higher front facet reflectivity. Specifically, we use devices with a segmented p-contact to directly measure the current distribution along the resonator and compare this with laser simulation. Devices with a 6 mm resonator show 29% more current at the front than back, twice as large as the 15% current non-uniformity in devices with a 3 mm resonator. In contrast, increased front facet reflectivity (20% rather than 0.8%) is shown to almost halve the current non-uniformity from 29% to 18% in devices with a 6 mm resonator and reduces power saturation. Although the magnitude of current non-uniformity in experiment and theory is broadly consistent, in experiment, an additional divergence is seen in current flow (and hence recombination rate) near the facets, and earlier power saturation occurs. Finally, we discuss the possible saturation mechanisms that are not included in the simulation.

47 OTHER INSTRUMENTATION↗

Deficiencies in compression and yield in x-ray-driven implosions

This paper analyzes x-ray–driven implosions that are designed to be less sensitive to 2-D and 3-D effects in hohlraum and capsule physics. Key performance metrics including the burn-averaged ion temperature, hot-spot areal density, and fusion yield are found to agree with simulations where the design adiabat (internal pressure) is multiplied by a factor of 1.4. Furthermore, these results motivate the development of a simple model for interpreting experimental data, which is then used to quantify how improvements in compression could help achieve ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental validation of shock propagation through a foam with engineered macro-pores

The engineered macro-pore foam provides a new way to study thermonuclear burn physics by utilizing capsules containing deuterated (D) foam and filling tritium (T) gas in the engineered macro-pores. The implosion of a thermonuclear capsule filled with an engineered macro-pore foam will be complex due to the interaction of a shock wave with the engineered macro-pores. It is our goal to quantify how substantially complex foam structures affect the shape of shock and bulk shock speed. A cylinder-shape shock tube experiment has been designed and performed at the Omega Laser Facility. In order to examine how a foam structure will affect shock propagation, we performed several tests varying (1) engineered macro-pore size, (2) average foam density, and (3) with/without neopentane (C 5 H 12 ) gas. X-ray radiographic data indicate that shock speed through engineered macro-pore foams depends strongly on average foam density and less on pore size. In this work, experimental shock propagation data helped guide two numerical simulation approaches: (1) a 2D simulation with homogenizing foams rather than explicitly simulating engineered macro-pores and (2) a 2D toroidal-pore approximation adopting a toroidal-tube geometry to model engineered macro-pores.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of asymmetric explosive density evolution in the deflagration-to-detonation transition for porous explosives

We report deflagration-to-detonation transition phenomena in cylindrical columns of porous explosives based on octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX). The HMX powder and LX-14 (95.5 wt. % HMX and 4.5 wt. % estane polyurethane) prills were observed with high-speed cameras and flash x rays under conditions that allow for precise scaling and registration of the recorded images. These experimental results test the conventional assumption of axisymmetric, effective one-dimensional burning. Polycarbonate confinement tubes allowed this direct imaging with steel spheres employed as radio-opaque fiducials. Spheres embedded in the explosive column revealed the displacement of the explosive prior to deflagration. X-ray and fast camera images show that a dense plug of compact explosive develops ahead of the deflagration front. The internal fiducials register both internal changes in the plug and changes in the position between successive x-ray images. Tantalum witness strips placed on the inner wall of the explosive channel revealed radial expansion of the tube bore at successive x-ray exposures, indicating the internal pressure of about 100 MPa just before detonation. Images indicate that the deflagration is not only asymmetric but also involves newly observed transient phenomena just before detonation. Furthermore, we typically see a brief but continuous increase in deflagration speed and the formation of a deflagration channel that bypasses the plug just prior to detonation. We hypothesize that these phenomena play an important role in the deflagration-to-detonation transition itself.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New methodology for measuring electron density perturbations caused by plasma coherent modes using profile reflectometry: Magnitudes and radial profiles in DIII-D

New capabilities of fast-sweep frequency-modulated profile reflectometry are explored to measure electron density ne perturbation magnitudes and radial profiles due to plasma coherent modes in DIII-D. The first approach is based on the frequency analysis of phase perturbations associated with high frequency (~MHz) Alfvén eigenmodes (AEs). The measurement of ~5.5 MHz fast-ion-driven global Alfvén eigenmodes (GAEs) is demonstrated in a neutral beam-heated DIII-D plasma. The GAE induced a broad radial distribution of phase perturbations in the profile reflectometer data. Analysis of these data determined the effective cutoff location displacement and the estimated ne fluctuation profile. In the second approach, high resolution ne profiles are used directly to determine the radial structure of ne perturbations due to a neo-classical tearing mode. These new measurements broaden the application of profile reflectometry and advance the development of AE spectroscopy as a tool for non-invasive diagnosis of fast-ion-driven modes in DIII-D and burning plasmas such as ITER.

47 OTHER INSTRUMENTATION↗

Proof-of-concept of a neutron time-of-flight ellipsoidal detector

The time-resolved measurement of neutrons emitted from nuclear implosions at inertial confinement fusion facilities is used to characterize the fusing plasma. Several significant quantities are routinely measured by neutron time-of-flight (nToF) detectors in these experiments. Current nToF detectors use scintillators as well as solid-state Cherenkov radiators. The latter has an inherently faster time response and can provide a co-registered γ-ray measurement as well as improved precision in the bulk hot-spot velocity. This work discusses a nToF ellipsoidal detector that also utilizes a solid-state Cherenkov radiator. The detector has the potential to achieve a fast instrument response function allowing for characterization of the γ-ray burn history as well as the ability to field the detector closer to the fusion source. Proof-of-concept testing of the nToF ellipsoidal detector has been conducted at the National Ignition Facility using commercial optics. A time-resolved neutron signal has been measured from the diagnostic. Preliminary simulations corroborate the results.

Jeet, J. (ORCID:0000000291091563)↗