Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion behavior”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Electrochemical mechanisms of leakage‐current in photovoltaic modules

Abstract This paper analyzes the mechanisms and pathways for leakage current flow observed in Si photovoltaic modules subjected to high temperature and humidity and a large voltage bias with respect to ground. The current inside of the frame is in the form of electron motion, but in the glass and polymer, it is in at least a large part attributable to the movement of ions. When the mode of current flow changes from electronic to ionic conduction, electrochemical reactions will take place at the interface. This can include reactions that produce volatile chemical species like H 2 , CO x , and O 2 , along with ionic species such as OH − and H 3 O + . Here, we see evidence of the importance of different charge carriers with different diffusion rates and the influence of electrochemical processes involved. The application of negative voltage to the cell circuit affects the resistivity of glass producing surfaces with poor conductivity but with some increases in the electrochemical potential producing complicated interactions that are important when the voltage is changed. In the polymer, there is the development of a space charge region and a chemical gradient providing two oppositional forces to current flow, which when released create a complicated discharge process. Here, we give a basic understanding of the charge/discharge of PV cells highlighting how the specific mechanisms are important in understanding some of the degradation processes in PV modules. We find that there is evidence of multiple significant charge carrier species with different diffusion time scales. The glass/polymer interface forms a depleted region of higher resistance after prolonged exposure to current. Charge also builds up at the polymer to cell antireflective coating interface and mostly flows to the gridlines to experience electrochemical reactions. These complexities result in non‐linear behavior where the apparent resistivities of the different layers change during charge/discharging processes, making the modeling of the current flow extremely difficult.

14 SOLAR ENERGY↗
DOE OSTI2020

Effect of Pressure and Thermal Cycling on Long-Term Oxidation in CO 2 and Supercritical CO 2

Concentrating solar power plant designers are interested in supercritical CO 2 (sCO 2 ) for the power block to achieve > 50% electrical efficiency at > 700 °C. The goal of this project was to develop a long-term (> 100 kh) lifetime model for sCO 2 compatibility using 10–15 kh laboratory exposures. Three Ni-based alloys (625, 282 and 740H) and an advanced austenitic stainless steel were evaluated here in long-term exposures at 700–800 °C using 500-h cycles in laboratory air, 0.1 MPa industrial grade (IG) CO 2 and 30 MPa supercritical IG CO 2 and using 10-h cycles in 0.1 MPa IG CO 2 and O 2 . Mass change data and quantification of the oxide scale thickness and depth of internal attack after 1000–10,000 h exposures at 750 °C indicate that these materials are compatible with the sCO 2 environments with modeling used to predict long-term behavior. Comparison of the 0.1 and 30 MPa 500-h cycle results did not show a significant effect of pressure on the reaction, and no significant internal carburization was observed under these conditions, even for the stainless steel, suggesting that chromia scales may be better C diffusion barriers than expected. For the Ni-based alloys, thermal cycling to simulate the solar duty cycle did not result in scale spallation after 15 kh in 10-h cycles or 4 kh in 1-h cycles at 750 °C. However, the stainless steel specimens formed an Fe-rich oxide after ~ 1500-h cumulative exposure time in both 1- and 10-h cycles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Empirical and mechanistic transient fission gas release model for high-burnup LOCA conditions

In response to the nuclear industry desire to extend burnup beyond current licensing practices, the US Nuclear Regulatory Commission (NRC) released a research information letter (RIL) that provides a basis for analyzing fuel fragmentation, relocation, and dispersal (FFRD) in light-water reactors. Of the five elements discussed, the most ambiguous is the significance of transient fission gas (FG) release (FGR) (tFGR) and its effects on fuel performance under loss-of-coolant accident conditions. In fresh fuel, FG migration and eventual release is primarily governed by diffusion-based mechanisms at higher temperatures (>1,000 °C). However, the mechanisms governing FGR changes as burnup increases. More recent research indicates that FGR increases as burnup increases, specifically under temperature transient conditions, and this release occurs at lower temperatures with a new release mechanism. This behavior has been attributed to microcracking in the fuel and is likely related to microstructure embrittlement with the presence of over pressurized FG bubbles. The NRC RIL outlines the complexity of the phenomenon and a need for a deeper understanding to adequately address FFRD for regulatory application. Here, therefore, this manuscript intends to summarize the publicly available tFGR data and discuss the observed dependencies (e.g., burnup, heating rate, sample geometry, terminal temperature). An empirical model has been developed and benchmarked against recently published experimental data. However, this empirical model is limited to conditions for which fitting data exist and, therefore, a high-level discussion is included to provide a roadmap for atomistically-informed multiscale modeling in conjunction with experimental data collection to develop a mechanistic tFGR model widely applicable to a broad range of nuclear fuel conditions at high burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High Temperature Oxidation Behavior of Aluminized Haynes 230

The effect of reactive air aluminizing (RAA) on the high-temperature oxidation behavior of the H230 superalloy is studied at 850 °C for 600 hours. Mass gain values are measured and microstructure is studied by XRD, SEM, and EDX. Post oxidation mass gain for the aluminized alloy is reduced compared to wrought alloy. A three-zone microstructure of diffusion coating is observed in RAA H230. Formation and growth of alumina scale during RAA process and oxidation is attributed to the complementary roles of ß-NiAl + ?'-Ni3Al phases in the composite coating.

Rashidi, S↗
DOE OSTI2023

Interfacial microstructure and mechanical properties of rotary inertia friction welded dissimilar 422 martensitic stainless steel to 4140 low alloy steel joints

In this work, dissimilar rotary inertia friction welds between 422 martensitic stainless steel and 4140 martensitic low-alloy steel were made to fabricate prototype heavy-duty diesel engine pistons. The influence of the inertia friction welding process and post weld heat treatment (PWHT) temperature on the interfacial microstructure evolutions and corresponding effects on mechanical properties of the 422/4140 welds were evaluated in detail. Carbon diffused from the 4140 side to the 422 side during PWHT at 650 °C for 1.5 h, causing the formation of a hard carbide-rich layer on the 422 side, and a softer but discontinuous C-depleted layer the 4140 side. PWHT at 700 °C for 1.5 h greatly accelerated C diffusion across the interface relative to 650 °C, resulting in a thicker hard carbide-rich layer and a relatively thick and continuous layer of coarse C-depleted grains (ferrite) on the 4140 side. In addition, the PWHT temperature greatly influenced the tensile properties and fracture behavior of the welds, with the 650 °C PWHT-ed samples failing predominately in a ductile manner in the 4140 heat affected zone during tensile testing. Conversely, the 700 °C PWHT specimens exhibited a strength reduction compared with the 650 °C PWHT specimens because of additional coarsening of the interfacial ferrite layer and softening of the base materials during PWHT, with brittle fracture between the hard and soft layers the predominate failure mechanism. Based on the findings, a reduced PWHT temperature and/or time, minimizing the hardness differential of the base metals, and pre-heating the 422 steel prior to welding are the potential pathways to achieve a more optimal balance between desirable tempering and stress relief of the weld microstructure and undesirable C migration across the weld interface, and to reduce the strength mismatch across the weld.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Synthesis and Properties of the Helium Clathrate and Defect Perovskite [He 2– x $^\square$ x ][CaNb]F 6

The defect double perovskite [He 2– x $^\square$ x ][CaNb]F 6 , with helium on its A-site, can be prepared by the insertion of helium into ReO 3 -type CaNbF 6 at high pressure. Upon cooling from 300 to 100 K under 0.4 GPa helium, ~60% of the A-sites become occupied. Helium uptake was quantified by both neutron powder diffraction and gas insertion and release measurements. After the conversion of gauge pressure to fugacity, the uptake of helium by CaNbF 6 can be described by a Langmuir isotherm. The enthalpy of absorption for helium in [He 2– x $^\square$ x ][CaNb]F 6 is estimated to be ~+3(1) kJ mol –1 , implying that its formation is entropically favored. Helium is able to diffuse through the material on a time scale of minutes at temperatures down to ~150 K but is trapped at 100 K and below. The insertion of helium into CaNbF 6 reduces the magnitude of its negative thermal expansion, increases the bulk modulus, and modifies its phase behavior. On compressing pristine CaNbF 6 , at 50 and 100 K, a cubic (Fm$\overline{3}$m) to rhombohedral (R$\overline{3}$) phase transition was observed at <0.20 GPa. However, a helium-containing sample remained cubic at 0.4 GPa and 50 K. CaNbF 6 , compressed in helium at room temperature, remained cubic to >3.7 GPa, the limit of our X-ray diffraction measurements, in contrast to prior reports that upon compression in a nonpenetrating medium, a phase transition is detected at ~0.4 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Evaluation of Common Thermoplastic Polymers in High-Pressure Cycling Hydrogen Under Ambient and Cold Environments as Applicable to the Hydrogen Infrastructure

Thermoplastic polymers are required to perform unfailingly under stringent conditions of changing pressures (35 MPa to 70 MPa) and temperatures (-40°C to +85°C) in storage and fueling operations. Although diffusivity in glassy polymers appear to be an order of magnitude lower than elastomers, polymer microstructural attributes such as degree of crystallinity and presence of polar and non-polar groups on the main chain with and without branching can play a substantial role in influencing their behaviors in hydrogen environments. In the work described here, the effect of high-pressure hydrogen cycling on PEEK, PTFE, PA11, HDPE, POM (and different commercial grades of these polymers) under ambient and cold (-40°C) conditions is addressed. Ex-situ characterization for polymer changes included density measurements for impact of hydrogen retention, hardness changes using nanoindentation, storage modulus and glass transition changes using dynamic mechanical and thermal analysis (DMTA), degree of crystallization with Differential Scanning Calorimetry (DSC) and tensile testing following ASTM D412. In other evaluations, solid-state nuclear magnetic resonance (NMR), attenuated total reflectance Fourier transform infra-red spectroscopy (ATRFTIR), X-ray diffraction (XRD), and X-ray CT results are presented. The overall impact of these evaluations is to establish a technical basis for the behaviors of common thermoplastics in cycling hydrogen environments under ambient and cold temperatures while establishing the relationship between polymer structure-based properties and hydrogen transport effects.

Menon, Nalini C.↗
DOE OSTI2021

The Corrected Distortion model for Lagrangian spray simulation of transcritical fuel injection

In this work, we present a detailed implementation and validation of the droplet modeling framework proposed by Dahms and Oefelein (2016) into the engine commercial CFD software CONVERGE using the User Defined Function (UDF) interface. The model accounts for the nonlinear deformation and oscillation experienced by liquid spray droplet injected into high pressure and temperature. Lagrangian spray simulations of Engine Combustion Network (ECN) Spray A are performed. Model validation against standard experimental measurements of liquid velocity, vapor mixture fraction is conducted. To perform more rigorous model validation, new experimental measurements based on Diffused Back Illumination (DBI) are introduced. The new measurements are processed for Projected Liquid Volume (PLV), which offers as close as possible one-to-one model validation for liquid penetration while offering new insights into the spray physics. Comparison with a One-D model based on adiabatic mixing theory by Siebers (1999) and Desantes et al. (2007) are also conducted. Through these model validation exercises, it is shown that the new framework improves liquid-phase penetration predictions, following a tendency for enhanced evaporation, compared to the standard approach for both Reynolds Average Navier Stokes (RANS) and Large Eddy Simulation (LES). Additionally, at the liquid length, maximum mixture fraction values predicted by the new approach are in good agreement those of an adiabatic mixing model. Qualitative analysis of the spray behaviors during the early stage of the injection process reveals that the proposed framework predicts significant increase in droplet evaporation rate with lower drop drag compared to the current standard approach.

42 ENGINEERING↗
DOE OSTI2022

Behavior of B in passivating gels formed on International Simple Glass in acid and basic pH

International Simple Glass altered for years in silica saturated conditions develops a passivating gel, which retains some mobile elements such as B and Ca in its inner region. It has been suggested that the retention of these elements impacts the dissolution rate of the glass. However, the elements' retention mechanism, speciation, and diffusion properties are still unknown. This study finds that simultaneous presence of B and Ca in the solution can dramatically slow down glass alteration. However, splitting the B and Ca into two different solutions with other glass elements did not give similar protection against glass alteration. These experiments along with other experiments conducted in basic or acid pH, in which isotopically tagged B can diffuse in gels revealed that (i) water diffusion is not dramatically affected by the retention of B and Ca in passivating gels, and (ii) slow diffusion and high retention of B is observed only in gels developed in basic pH during the tracing experiments. Finally, this last mechanism therefore appears as one of the many mechanisms controlling the glass dissolution rate in the residual regime.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fluid and hybrid simulations of the ionization instabilities in Hall thruster

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. Further, the predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coherent band-edge oscillations and dynamic longitudinal-optical phonon mode splitting as evidence for polarons in perovskites

The coherence of collective modes, such as phonons and polarons, and their modulation of electronic states is long sought in complex systems, which is a crosscutting issue in photovoltaics and quantum electronics. In photovoltaic cells and lasers based on metal halide perovskites, the presence of polarons, i.e., photocarriers dressed by the macroscopic motion of charged lattice, assisted by terahertz (THz) longitudinal-optical (LO) phonons, has been intensely studied yet is still debated. This may be key for explaining the remarkable properties of the perovskite materials, e.g., defect tolerance, long charge lifetimes, and diffusion lengths. In this work we use the intense single-cycle THz pulse with peak electric field up to E THz = 1000 kV/cm to drive coherent polaronic band-edge oscillations at room temperature in CH 3 NH 3 PbI 3 (MAPbI 3 ). We reveal the oscillatory behavior is dominated by a specific quantized lattice vibration mode at $ω$ LO ~ 4 THz, which is both dipole and momentum forbidden. THz-driven coherent polaron dynamics exhibits distinguishing features: room temperature coherent oscillations at ω LO longer than 1 ps in both single crystals and thin films, mode-selective modulation of different band-edge states assisted by electron-phonon interaction, and dynamic mode splitting at low temperature due to entropy and anharmonicity of organic cations. Our results demonstrate intense THz-driven coherent band-edge modulation is a powerful probe of electron-lattice coupling phenomena and polaronic quantum control in perovskites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonlinear ion transport mediated by induced charge in ultrathin nanoporous membranes

Ultrathin membranes with nanoporous conduits show promise for ionic separations and desalination applications, but the mechanisms underlying the nonlinear ionic transport observed in these systems are not well understood. Here, in this work, we demonstrate how induced charge at membrane interfaces can lead to nonlinear ionic transport and voltage-dependent conductance through such channels. The application of an electric field on a polarizable membrane leads to induced charges at the membrane interfaces. The induced charges in turn are screened by diffuse charges in the electrolyte, which are acted upon by the electric field. For extremely thin membranes, the induced charge effect can be significant even for moderate applied voltages commonly used in experiments. We apply a continuum Poisson-Nernst-Planck model to characterize the current-voltage behavior of ultrathin membranes over a wide parameter space. The predictions of the model are compared to recent experiments on graphene and MoS 2 membranes in an electric field. We expect the role of induced charge to be especially pronounced in the limit of atomically thin membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chaotic Dynamics Driven by Particle-Core Interactions

High-intensity beams in modern linacs are frequently encircled by diffuse halos, which drive sustained particle losses and result in gradual degradation of accelerating structures. In large part, the growth of halos is facilitated by internal space-charge forces within the beams, and detailed characterization of this process constitutes an active area of ongoing research. A partial understanding of dynamics that ensue within space-charge dominated beams is presented by the particle-core interaction paradigm – a mathematical model wherein single particle dynamics, subject to the collective potential of the core, are treated as a proxy for the broader behavior of the beam. In this work, we investigate the conditions for the onset of large-scale chaos within the framework of this model, and demonstrate that the propensity towards stochastic evolution is strongly dependent upon the charge distribution of the beam. In particular, we show that while particle motion within a uniformly charged beam is dominantly regular, rapid deterministic chaos readily arises within space-charge dominated Gaussian beams. Importantly, we find that for sufficiently high values of the beam’s space charge and beam pulsation amplitude, enhanced chaotic mixing between the core and the halo can lead to an enhanced radial diffusion of charged particles. We explain our results from analytic grounds by demonstrating that chaotic motion is driven by the intersection of two principal resonances of the system, and derive the relevant overlap conditions. Additionally, our analysis illuminates a close connection between the mathematical formulation of the particle-core interaction model and the Andoyer family of integrable Hamiltonians

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Response of 11 B enriched ZrB 2 ultra-high temperature ceramic to neutron irradiation at elevated temperatures

ZrB 2 , an ultra-high temperature ceramic (UHTC) is being considered for use in fusion reactor first-wall structures, yet its response to irradiation remains poorly understood. This study employed scanning/transmission electron microscopy (S/TEM), synchrotron X-ray diffraction (XRD), finite element calculations, and thermal property measurements to thoroughly investigate the neutron-irradiation effects on 11 B-enriched ZrB 2 . Neutron irradiations were conducted at 220 °C and 620 °C, with a neutron fluence of 2.2 × 10 25 neutron/m 2 (energy > 0.1 MeV), resulting in 3.9 dpa and 4200 appm He. The study revealed the unusual prevalence of prism loops and a > c anisotropic lattice swelling, likely linked to the low c/a ratio of ZrB 2 , leading to grain boundary microcracking. Reducing the grain sizes was effective in reducing intergranular cracking and macroscopic swelling. The observation of cavities in ZrB 2 irradiated at 620 °C, as opposed to 220 °C, prompts questions about the temperature at which vacancies in ZrB 2 become mobile, and the role of neutron absorption by 10 B in elevating irradiation temperatures. Isotopic enrichment in 11 B proves to be a viable strategy for mitigating helium production in transition-metal diborides, which is a critical consideration for nuclear applications. Irradiation-induced defects reduce the thermal diffusivity and conductivity of ZrB 2 by a factor of 4–9, which has important implications for its role as a plasma-facing material in fusion reactors that drive high heat fluxes through first-wall materials. Here, this comprehensive study lays the foundation for understanding ZrB 2 behavior under neutron irradiation and highlights important phenomena to consider for various material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗