Engineering topics

Yang, Lin H.

Publications and source records attributed to Yang, Lin H..

3 records

DOE OSTI2021

Mechanical and Thermal Properties for Uranium and U–6Nb Alloy from First-Principles Theory

Elasticity, lattice dynamics, and thermal expansion for uranium and U–6Nb alloy (elastic moduli) are calculated from density functional theory that is extended to include orbital polarization (DFT+OP). Introducing 12.5 at.% of niobium, substitutionally, in uranium softens all the cii elastic moduli, resulting in a significantly softer shear modulus (G). Combined with a nearly invariant bulk modulus (B), the quotient B/G increases dramatically for U–6Nb, suggesting a more ductile material. Lattice dynamics from a harmonic model coupled with a DFT+OP electronic structure is applied for α uranium, and the obtained phonon density of states compares well with inelastic neutron-scattering measurements. The Debye temperature associated with the lattice dynamics falls within the range of experimentally observed Debye temperatures and it also validates our quasi-harmonic (QH) phonon model. The QH Debye–Grüneisen phonon method is combined with a DFT+OP electronic structure and used to explore the anisotropic thermal expansion in α uranium. The anomalous negative thermal expansion (contraction) of the b lattice parameter of the α-phase orthorhombic cell is relatively well reproduced from a free-energy model consisting of QH-phonon and DFT+OP electronic structure contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Role of opacity at the 9 keV back lighter energy used in measuring the equation of state of boron at pressures up to a Gbar

Many experiments have been conducted at the National Ignition Facility to measure the equation of state Hugoniot of plastic, boron, and diamond at extreme pressures up to a Gbar. The “Gbar” design employs a strong spherically converging shock launched through a solid ball of material using a hohlraum radiation drive. The shock front conditions are characterized using X-ray radiography, typically at energies near 9 keV. Here we examine how the opacity of boron at 9 keV changes at high pressures and temperatures. Understanding this is vital to unfolding the density in the shock front as pressures exceed 100 Mbar. We compare opacity calculations from a number of methods including the legacy XSN opacity tables, super transition array (STA), average atom (AVAT), and detailed configuration accounting (DCA) methods. We examine how the changing opacity is correlated with the K-shell occupation calculated using the opacity methods as well as other electronic structure methods such as MECCA, Purgatorio, and an all-electron Kohn-Sham density functional theory calculation that uses optimized all-electron norm-conserving Vanderbilt pseudopotentials (ONCV). We also examine the relative contribution of free-free opacity to the overall opacity, which is dominated by bound-free absorption.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗