Search NASA⌕ Search

SEARCH · Search NASA

Results for “High Energy Density Physics”

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 253 records · Page 14

A new approach to include electron interaction effects in super transition array opacity theory

A method is presented for the improved calculation of super-shell partition functions which include the repulsive electron–electron interaction energy terms in the Boltzmann factor. Heretofore these interaction terms were approximately treated via the use of Feynman–Jensen inequalities. Further, such investigations are of particular interest for the super-transition-array approach of hot-plasma radiative opacity.

74 ATOMIC AND MOLECULAR PHYSICS↗

Towards the first plasma-electron screening experiment

The enhancement of fusion reaction rates in a thermonuclear plasma by electron screening of the Coulomb barrier is an important plasma-nuclear effect that is present in stellar models but has not been experimentally observed. Experiments using inertial confinement fusion (ICF) implosions may provide a unique opportunity to observe this important plasma-nuclear effect. Herein, we show that experiments at the National Ignition Facility (NIF) have reached the relevant physical regime, with respect to the density and temperature conditions, but the estimated impacts of plasma screening on nuclear reaction rates are currently too small and need to be increased to lower the expected measurement uncertainty. Detailed radiation hydrodynamics simulations show that practical target changes, like adding readily available high-Z gases, and significantly slowing the inflight implosion velocity, while maintaining inflight kinetic energy, might be able to push these conditions to those where plasma screening effects may be measurable. We also perform synthetic data exercises to help understand where the anticipated experimental uncertainties will become important. But challenges remain, such as the detectability of the reaction products, non-thermal plasma effects, species separation, and impacts of spatial and temporal gradients. This work lays the foundation for future efforts to develop an important platform capable of the first plasma electron screening observation.

inertial confinement fusion↗

Preparations for a European R&D roadmap for an inertial fusion demo reactor

A European consortium of 15 laboratories across nine nations have worked together under the EUROFusion Enabling Research grants for the past decade with three principle objectives. These are: (a) investigating obstacles to ignition on megaJoule-class laser facilities; (b) investigating novel alternative approaches to ignition, including basic studies for fast ignition (both electron and ion-driven), auxiliary heating, shock ignition, etc.; and (c) developing technologies that will be required in the future for a fusion reactor. This paperr presents a brief overview of these activities, along with new calculations relates the concept of auxiliary heating of inertial fusion targets, and provides possible future directions of research and development for the updated European Roadmap that is due at the end of 2020. This article is part of a discussion meeting issue ‘Prospects for high gain inertial fusion energy (part 2)’.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimizing Simulation Fidelity in Direct-Drive Inertial Confinement Fusion with Cassio

Recently, the National Ignition Facility (NIF) demonstrated that inertial confinement fusion (ICF) is capable to achieve thermonuclear (TN) ignition in the laboratory, making it a crucial method on the path to replicate the Sun’s power production mechanism on Earth. However, the physics governing the high-energy density environments is very complex and remains a challenge to fully understand and model. For example, dopants in the TN fuel are important diagnostic tools to extract the thermodynamic conditions of the plasma. However, if their concentration is chosen too high, they can significantly degrade the performance of an ICF capsule. In this study, we use the Los Alamos National Laboratory radiation-hydrodynamics code Cassio to model ICF implosions of capsules which contain deuterium fuel with high-Z dopants like Krypton and Argon from the high-Z campaign conducted 15 years ago. We focus on how chosen computational and physics parameters influence the implosion outcomes. By systematically changing the resolution of the computational mesh and the photon energies as well as modifying settings for the laser drive and TN fuel pre-heat effects, we assess the impact on experimentally measured performance metrics like neutron production from TN burn and x-ray emission during the implosion. Our results will help to improve the fidelity of simulations and guide future numerical studies and experimental designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dense plasma opacity via the multiple-scattering method

The calculation of the optical properties of hot dense plasmas with a model that has self-consistent plasma physics is a grand challenge for high energy density science. Here we exploit a recently developed electronic structure model that uses multiple scattering theory to solve the Kohn-Sham density functional theory equations for dense plasmas. Here we calculate opacities in this regime, validate the method, and apply it to recent experimental measurements of opacity for Cr, Ni, and Fe. Good agreement is found in the quasicontinuum region for Cr and Ni, while the self-consistent plasma physics of the approach cannot explain the observed difference between models and the experiment for Fe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First-principles equation of state database for warm dense matter computation

We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF , B 4 C , BN , CH 4 , CH 2 , C 2 H 3 , CH , C 2 H , MgO , and MgSiO 3 . For all these materials, we provide the pressure and internal energy over a density-temperature range from ~ 0.5 to 50 g cm - 3 and from ~ 10 4 to 10 9 K, which are based on ~ 5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L - and K -shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. Furthermore, we derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H 2 O , H 2 O 2 , Al 2 O 3 , CO , and CO 2 to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comprehensive analytical model of the dynamic 𝑍 pinch

In this work we present an analytical 1D axisymmetric model describing the evolution of the dynamic 𝑍 pinch. This model is capable of predicting the trajectories of the imploding sheath's magnetic piston and preceding shock front, along with the velocity, pressure, density, and magnetic field profiles, for any time-dependent current, spatially varying initial density profile, and weak initial axial field. The implosion is divided into stages, with each stage described by a set of coupled ordinary differential equations derived from the ideal MHD equations. Comparisons with experimental data from the COBRA pulsed-power facility are quite promising and imply this model could prove useful in designing and analyzing future pulsed-power experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Average-atom calculations of bound-free and free-free cross sections in dense plasmas

Cross sections for photo-induced bound-free and free-free transitions in plasmas are evaluated in the average-atom approximation and applied to determine opacities of dense plasmas of light elements. Parameters characterizing the plasmas (chemical potential, average ionic charge, free electron density, bound and continuum wave functions, and occupation numbers) are obtained from the average-atom model. Lowest-order calculations of the free-free cross sections, which diverge in the low-frequency limit, are regularized by accounting for the finite electron-ion relaxation time. The resulting analysis provides the basis for average-atom studies of plasma opacities. Such studies are presented for dense lithium, beryllium, boron, and carbon. Applications are given to Rosseland mean opacities of dense hydrogen and deuterium plasmas and to comparisons of free-free to bound-free opacities in shock-compressed plasmas. Average-atom cross section and opacity calculations are extended to plasmas consisting of more than one ionic species, boron nitride, polystyrene, and a composite H, He, C plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimized recursion relation for the computation of partition functions in the superconfiguration approach

Partition functions of a canonical ensemble of non-interacting bound electrons are a key ingredient of the super-transition-array approach to the computation of radiative opacity. A few years ago, we published a robust and stable recursion relation for the calculation of such partition functions. In this work, we propose an optimization of the latter method and explain how to implement it in practice. The formalism relies on the evaluation of elementary symmetric polynomials, which opens the way to further improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A fast approximation to supershell partition functions

A formula for supershell partition functions, which play a major role in the Super Transition Array approach to radiative-opacity calculations, is derived as a functional of the distribution of energies within the supershell. It consists in an alternative expansion for an arbitrary number of electrons or holes which also allows for quick approximate evaluations with truncated number of terms in the expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast approximation to supershell partition functions: Explicit forms of the coefficients

In a previous work, we derived a formula for supershell partition functions, which are the cornerstone of the Super-Transition-Array approach to radiative-opacity calculations. The new expression takes the form of a functional of the distribution of energies within the supershell and allows for fast and accurate computations, truncating the number of terms in the expansion. The latter involves coefficients (denoted Γ $k$ ) for which we obtained a recursion relation. In the present short paper, as a complement of the previous one, we give an explicit formula for the coefficients Γ $k$ . Another recursion relation is also provided, as well as an alternative expression involving Bell polynomials. In conclusion, the connections with cycle indexes of permutation groups and with elementary symmetric polynomials are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the computation of moments in the Super-Transition-Arrays model for radiative opacity calculations

In the Super-Transition-Array statistical method for the computation of radiative opacity of hot dense matter, the moments of the absorption or emission features involve partition functions with reduced degeneracies, occurring through the calculation of averages of products of subshell populations. Here, in the present work, we discuss several aspects of the computation of such peculiar partition functions, insisting on the precautions that must be taken in order to avoid numerical difficulties. In a previous work, we derived a formula for supershell partition functions, which takes the form of a functional of the distribution of energies within the supershell and allows for fast and accurate computations, truncating the number of terms in the expansion. The latter involves coefficients for which we obtained a recursion relation and an explicit formula. We show that such an expansion can be combined with the recurrence relation for shifted partition functions. We also propose, neglecting the effect of fine structure as a first step, a positive-definite formula for the Super-Transition-Array moments of any order, providing an insight into the asymmetry and sharpness of the latter. The corresponding formulas are free of alternating sums. Several ways to speed up the calculations are also presented.

74 ATOMIC AND MOLECULAR PHYSICS↗

Inertial-confinement fusion-plasma-based cross-calibration of the deuterium-tritium γ-to-neutron branching ratio

The deuterium-tritium (D-T) γ -to-neutron branching ratio [ H 3 ( d , γ ) He 5 / H 3 ( d , n ) He 4 ] has been determined previously under inertial-confinement fusion (ICF) conditions and in beam-target based experiments. In the former case, neutron-induced backgrounds are mitigated compared to the latter due to the short pulse nature of ICF implosions and the use of gas Cherenkov γ -ray detectors. Furthermore, an added benefit of ICF based measurements is the ability to achieve lower center-of-mass energies as compared to accelerators. Previous ICF based experiments however report a large uncertainty in the D-T γ -to-neutron branching ratio of ≈ 48 % , which arises from the necessity of an absolute detector calibration and/or a cross-calibration against the D- He 3 γ -to-proton branching ratio. A more precise value for the branching ratio based on data taken at the OMEGA laser facility is reported here, which relies on a cross-calibration against the better known C 12 neutron inelastic scattering cross section. A D-T branching ratio value of ( 4.6 ± 0.6 ) × 10 - 5 is determined by this method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Azimuthal magnetic field distribution in gas-puff Z -pinch implosions with and without external magnetic stabilization

An experimental study of the magnetic field distribution in gas-puff Z pinches with and without a preembedded axial magnetic field (B z0 ) is presented. Spatially resolved, time-gated spectroscopic measurements were made at the Weizmann Institute of Science on a 300 kA, 1.6 μs rise time pulsed-power driver. The radial distribution of the azimuthal magnetic field, B θ , during the implosion, with and without a preembedded axial magnetic field of B z0 = 0.26 T, was measured using Zeeman polarization spectroscopy. The spectroscopic measurements of B θ were consistent with the corresponding values of B θ inferred from current measurements made with a B-dot probe. One-dimensional magnetohydrodynamic simulations, performed with the code TRAC-II, showed agreement with the experimentally measured implosion trajectory, and qualitatively reproduced the experimentally measured radial B θ profiles during the implosion when B z0 = 0.26 T was applied. Simulation results of the radial profile of B θ without a preembedded axial magnetic field did not qualitatively match experimental results due to magneto-Rayleigh-Taylor (MRT) instabilities. Furthermore, our analysis emphasizes the importance of MRT instability mitigation when studying the magnetic field and current distributions in Z pinches. Discrepancies of the simulation results with experiment are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of ionic disorder on the principal shock Hugoniot

Here, the effect of ionic disorder on the principal Hugoniot is investigated using multiple scattering theory to very high pressure (Gbar). Calculations using molecular dynamics to simulate ionic disorder are compared to those with a fixed crystal lattice, for both carbon and aluminum. For the range of conditions considered here we find that ionic disorder has a relatively minor influence. It is most important at the onset of shell ionization and we find that, at higher pressures, the subtle effect of the ionic environment is overwhelmed by the larger number of ionized electrons with higher thermal energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗