Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density estimation”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Approach to Evaluating Reorganization Energies of Interfacial Electrochemical Reactions

Reaction rate coefficients for electron-transfer processes at the electrode–electrolyte interface are commonly estimated by using the Butler–Volmer equation, but their values are inaccurate beyond a few tenths of volts of overpotential. The Marcus–Hush–Chidsey (MHC) formalism yields correct asymptotic behavior of the rate coefficients vs applied overpotential but has complex dependencies on the redox system’s intrinsic parameters, which can be difficult to model or measure. In this work, we bridge the two kinetics formalisms to estimate the reorganization energy, one of the important parameters for the MHC formalism, and investigate its dependence on other intrinsic parameters such as activation barriers, electronic coupling strength, and the density of states of the electrode surface. We examine the sensitivity of the reorganization energy to these parameters, establish some general relationships for accurately predicting rate coefficients using the MHC formalism over a wide range of applied overpotentials, and compare this approach to calculating MHC rate constants with other empirical approaches for the mechanisms of CO 2 reduction on different metal electrode surfaces.

Butler−Volmer↗
DOE OSTI2020

Visualizing Orbital Content of Electronic Bands in Anisotropic 2D Semiconducting ReSe 2

Many properties of layered materials change as they are thinned from their bulk forms down to single layers, with examples including indirect-to-direct band gap transition in 2H semiconducting transition metal dichalcogenides as well as thickness-dependent changes in the valence band structure in post-transition-metal monochalcogenides and black phosphorus. Here, we use angle-resolved photoemission spectroscopy to study the electronic band structure of monolayer ReSe 2 , a semiconductor with a distorted 1T structure and in-plane anisotropy. By changing the polarization of incoming photons, we demonstrate that for ReSe 2 , in contrast to the 2H materials, the out-of-plane transition metal dz2 and chalcogen pz orbitals do not contribute significantly to the top of the valence band, which explains the reported weak changes in the electronic structure of this compound as a function of layer number. We estimate a band gap of 1.7 eV in pristine ReSe 2 using scanning tunneling spectroscopy and explore the implications on the gap following surface doping with potassium. Additionally, a lower bound of 1.4 eV is estimated for the gap in the fully doped case, suggesting that doping-dependent many-body effects significantly affect the electronic properties of ReSe 2 . Our results, supported by density functional theory calculations, provide insight into the mechanisms behind polarization-dependent optical properties of rhenium dichalcogenides and highlight their place among two-dimensional crystals.

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

Ultraviolet spectroscopy of pressurized and supercritical carbon dioxide

Abstract Carbon dioxide (CO 2 ) is prevalent in planetary atmospheres and sees use in a variety of industrial applications. Despite its ubiquitous nature, its photochemistry remains poorly understood. In this work we explore the density dependence of pressurized and supercritical CO 2 electronic absorption spectra by vacuum ultraviolet spectroscopy over the wavelength range 1455-2000 Å. We show that the lowest absorption band transition energy is unaffected by a density increase up to and beyond the thermodynamic critical point (137 bar, 308 K). However, the diffuse vibrational structure inherent to the spectrum gradually decreases in magnitude. This effect cannot be explained solely by collisional broadening and/or dimerization. We suggest that at high densities close proximity of neighboring CO 2 molecules with a variety of orientations perturbs the multiple monomer electronic state potential energy surfaces, facilitating coupling between binding and dissociative states. We estimate a critical radius of ~4.1 Å necessary to cause such perturbations.

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

Testing general relativity on cosmological scales at redshift z ∼ 1.5 with quasar and CMB lensing

ABSTRACT We test general relativity (GR) at the effective redshift $\bar{z} \sim 1.5$ by estimating the statistic EG, a probe of gravity, on cosmological scales $19 - 190\, h^{-1}{\rm Mpc}$. This is the highest redshift and largest scale estimation of EG so far. We use the quasar sample with redshifts 0.8 < z < 2.2 from Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 16 as the large-scale structure (LSS) tracer, for which the angular power spectrum $C_\ell ^{qq}$ and the redshift-space distortion parameter β are estimated. By cross-correlating with the Planck 2018 cosmic microwave background (CMB) lensing map, we detect the angular cross-power spectrum $C_\ell ^{\kappa q}$ signal at $12\, \sigma$ significance. Both jackknife resampling and simulations are used to estimate the covariance matrix (CM) of EG at five bins covering different scales, with the later preferred for its better constraints on the covariances. We find EG estimates agree with the GR prediction at $1\, \sigma$ level over all these scales. With the CM estimated with 300 simulations, we report a best-fitting scale-averaged estimate of $E_G(\bar{z})=0.30\pm 0.05$, which is in line with the GR prediction $E_G^{\rm GR}(\bar{z})=0.33$ with Planck 2018 CMB + BAO matter density fraction Ωm = 0.31. The statistical errors of EG with future LSS surveys at similar redshifts will be reduced by an order of magnitude, which makes it possible to constrain modified gravity models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Deep-field analytical calibration

The next generation of imaging surveys, including the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), Euclid, and the Nancy Grace Roman Space Telescope, will provide unprecedented constraints on cosmology using weak gravitational lensing. To fully exploit this statistical power, shear measurement methods must achieve sub- per cent accuracy while mitigating systematic biases from noise, the point-spread function (PSF), blending, and shear-dependent detection. The analytical calibration framework (AnaCal) has demonstrated such accuracy but requires adding noise to images, reducing effective depth. We introduce Deep-Field Analytical Calibration (DEEP-FIELD AnaCal), an extension of AnaCal that uses deep-field images to compute shear responses while preserving the statistical power of wide-field data. We validate DEEP-FIELD AnaCal on isolated and blended galaxy image simulations with LSST-like conditions, finding it meets the stringent requirement of multiplicative bias $|m| < 3\times 10^{-3}$ at 99.7 per cent confidence. Compared to standard AnaCal applied to wide-field images, DEEP-FIELD AnaCal increases the effective galaxy number density from 17 to 30 arcmin$^{-2}$ for simulated 10-yr LSST data. With deep fields $10\times$ longer than the wide field, we find pixel noise variance in shear estimation is reduced by 30 per cent and overall uncertainty by $\sim 25~{{\ \rm per\ cent}}$. Finally, using the LSST Deep Drilling Fields strategy, we assess sample variance and find an equivalent calibration uncertainty of $\lesssim 0.3~{{\ \rm per\ cent}}$. These results demonstrate that DEEP-FIELD AnaCal offers a promising path to achieve the required shear calibration for upcoming weak lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Multimodal fission from self-consistent calculations

When multiple fission modes coexist in a given nucleus, distinct fragment yield distributions appear. Multimodal fission has been observed in a number of fissioning nuclei spanning the nuclear chart, and this phenomenon is expected to affect the nuclear abundances synthesized during the rapid neutron-capture process (𝑟-process). In this study, we generalize the previously proposed hybrid model for fission-fragment yield distributions to predict competing fission modes and estimate the resulting yield distributions. Here, our framework allows for a comprehensive large-scale calculation of fission-fragment yields suited for 𝑟-process nuclear network studies. Nuclear density functional theory is employed to obtain the potential energy and collective inertia tensor on a multidimensional collective space defined by mass multipole moments. Fission pathways and their relative probabilities are determined using the nudged elastic band method. Based on this information, mass and charge fission yields are predicted using the recently developed hybrid model. Fission properties of fermium isotopes are calculated in the axial quadrupole-octupole collective space for three energy density functionals (EDFs). Disagreement between the EDFs appears when multiple fission modes are present. Within our framework, the UNEDF⁢1 HFB EDF agrees best with experimental data. Calculations in the axial quadrupole-octupole-hexadecapole collective space improve the agreement with the experiment for SkM*. We also discuss the sensitivity of fission predictions on the choice of EDF for several superheavy nuclei. Fission-fragment yield predictions for nuclei with multiple fission modes are sensitive to the underlying EDF. For large-scale calculations in which a minimal number of collective coordinates is considered, UNEDF⁢1 HFB provides the best description of experimental data, though the sensitivity motivates robust quantification of the uncertainties of the theoretical model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Reliability comparisons between additively manufactured and conventional SiC–Si ceramic composites

Abstract An additively manufactured reaction‐bonded silicon carbide ceramic composite is fabricated using a bimodal powder feedstock and the binder‐jetting printing technique. On fabricating, the ceramic is investigated to report its composition, mechanical, and thermal properties at room temperature and high temperatures (up to 750). The ceramic has a density of 2.76 g/cc, and shows a hardness of 21.74 GPa, and a flexure strength of 207.5 MPa at room temperature. The mechanical property of flexure strength is used to estimate its failure probability under applied stress, using a Weibull distribution. The mechanical properties and failure probabilities are compared with those of conventionally fabricated reaction‐bonded silicon carbide ceramics reported in the literature. These ceramics were fabricated using methods such as slip‐casting, compacting, and tape‐casting. The comparison is used to elucidate the advantages and areas of improvement of the present additive manufacturing technique for reaction‐bonded ceramics.

14 SOLAR ENERGY↗
DOE OSTI2026

Galaxy Clusters Selected via the Sunyaev-Zel'dovich Effect in 5 year data from the SPT-3G Main Survey

We report a new galaxy cluster catalog, selected using the thermal Sunyaev-Zel'dovich (SZ) effect, from 5 years of observations of the SPT-3G Main field. Drawn from arcminute-resolution data with white noise levels of 3.2, 2.5, and 8.9 $μ$K-arcmin at 95, 150, and 220 GHz, respectively, the sample consists of 8,892 cluster candidates detected above significance $ξ=4$, with an expected purity of $>82\%$ (4,480 at $ξ\ge5$ with purity $>99\%$). Using optical and infrared data we have confirmed 7,190 candidates as clusters. The sample spans a mass range $7.9 \times 10^{13}M_\odot/h_{70}< M_\textrm{500c} < 1.6 \times 10^{15}M_\odot/h_{70}$ with a median mass of $1.65 \times 10^{14}M_\odot/h_{70}$, and a redshift range of $0.0371$ and 271 at $z>1.5$. Compared to previous SZ cluster samples from South Pole Telescope and Atacama Cosmology Telescope data, the SPT-3G sample is highly consistent in mass and redshift but is significantly deeper, with per-cluster detection signal-to-noise 2-4 times higher and a cluster density of 4.5 confirmed clusters/deg$^2$. We cross match with eRASS1 cluster and point source catalogs, finding 1,279 and 1,319 matches, respectively. The SPT and eROSITA cluster mass estimates are in relatively good agreement. We perform a series of validation checks using both internal data splits and comparisons to external samples. These tests show increasing correlated (dusty) emission with redshift, with a $\sim17\times$ larger 220 GHz temperature increment for clusters at $z\sim1.5$ than $z\sim0.25$, but only weak evidence for correlated synchrotron emission. Finally, a number of clusters are flagged as candidate strong gravitational lenses.

Bleem, L. E. [Argonne (main); Chicago U., KICP; Ch↗
DOE OSTI2022

Updated Godiva-IV Benchmark Preview

A note of errata prepended to the Godiva-IV delayed-critical benchmark (HEU-MET-FAST- 086) identifies two corrections to be made to the model: The glory hole in the spindle should be made larger, and the height of the safety block should be made smaller (and therefore its density made larger). In addition, the safety block at its full-in position is closer to the inner subassembly plate than was modeled in the benchmark. These changes have been made to HEU-MET-FAST-086 Case 4 in order to estimate the effect on $\kappa$ eff and on the neutron flux spectrum. Using smaller separation, a smaller safety block, and a larger glory hole caused keff to increase by 453 ± 1 pcm from the benchmark. The latest nuclear data, ENDF/B-VIII.0, have also been used, causing $\kappa$ eff to increase another 37 ± 1 pcm. Flux spectra were compared in a modeled fission foil (near the center of Godiva-IV in the glory hole) and at three external (point) detectors. Within the fission foil, using ENDF/B-VIII.0 induced changes in the flux spectrum similar in size to the changes due to using smaller separation, a smaller safety block, and a larger glory hole. At the detectors, using ENDF/B- VIII.0 induced changes in the flux spectrum much larger than those due to changing the model. In other words, the corrections to the benchmark model cause a large increase in $\kappa$ eff , but the changes to the neutron flux spectrum are small compared to those caused by using the latest nuclear data. This study presents a preview of results expected during the reevaluation of the Godiva IV benchmark, but it is not a substitute for the full reevaluation.A note of errata prepended to the Godiva-IV delayed-critical benchmark (HEU-MET-FAST- 086) identifies two corrections to be made to the model: The glory hole in the spindle should be made larger, and the height of the safety block should be made smaller (and therefore its density made larger). In addition, the safety block at its full-in position is closer to the inner subassembly plate than was modeled in the benchmark. These changes have been made to HEU-MET-FAST-086 Case 4 in order to estimate the effect on $\kappa$ eff and on the neutron flux spectrum. Using smaller separation, a smaller safety block, and a larger glory hole caused $\kappa$ eff to increase by 453 ± 1 pcm from the benchmark. The latest nuclear data, ENDF/B-VIII.0, have also been used, causing $\kappa$ eff to increase another 37 ± 1 pcm. Flux spectra were compared in a modeled fission foil (near the center of Godiva-IV in the glory hole) and at three external (point) detectors. Within the fission foil, using ENDF/B-VIII.0 induced changes in the flux spectrum similar in size to the changes due to using smaller separation, a smaller safety block, and a larger glory hole. At the detectors, using ENDF/B- VIII.0 induced changes in the flux spectrum much larger than those due to changing the model. In other words, the corrections to the benchmark model cause a large increase in $\kappa$ eff , but the changes to the neutron flux spectrum are small compared to those caused by using the latest nuclear data. This study presents a preview of results expected during the reevaluation of the Godiva IV benchmark, but it is not a substitute for the full reevaluation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Theoretical Insights on the Fragmentation of Cellulosic Oligomers to Form Hydroxyacetone and Hydroxyacetaldehyde

The presence of heavy unknown oligomeric sugar products in bio-oil is evidenced in experimental results reported in the literature. In this paper, we study the fragmentation reactions yielding acetol and glycolaldehyde from oligomeric sugars following previous work on dehydration reactions to propose structures of these oligomers. Acetol and glycolaldehyde are primary products of cellulose fast pyrolysis but the fragmentation reaction mechanism of these compounds from oligomers merits further study. The density functional theory (DFT) approach was employed to study this reaction. Results revealed that acetol and glycolaldehyde fragments are favorably removed from the non-reducing end based on their thermodynamic stabilities. Theoretical FTIR and NMR spectra were calculated to aid in understanding the structures of the oligomeric sugars. Also, the thermodynamics and physical properties of these compounds were estimated using the Group Contribution Method (GCM). In conclusion, these properties are essential in the design of processing technologies and the upgrading of products.

10 SYNTHETIC FUELS↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

Increased adaxial stomatal density is associated with greater mesophyll surface area exposed to intercellular air spaces and mesophyll conductance in diverse C 4 grasses

Summary Mesophyll conductance ( g m ) is the diffusion of CO 2 from intercellular air spaces ( IAS ) to the first site of carboxylation in the mesophyll cells. In C 3 species, g m is influenced by diverse leaf structural and anatomical traits; however, little is known about traits affecting g m in C 4 species. To address this knowledge gap, we used online oxygen isotope discrimination measurements to estimate g m and microscopy techniques to measure leaf structural and anatomical traits potentially related to g m in 18 C 4 grasses. In this study, g m scaled positively with photosynthesis and intrinsic water‐use efficiency ( TE i ), but not with stomatal conductance. Also, g m was not determined by a single trait but was positively correlated with adaxial stomatal densities ( SD ada ), stomatal ratio ( SR ), mesophyll surface area exposed to IAS ( S mes ) and leaf thickness. However, g m was not related to abaxial stomatal densities ( SD aba ) and mesophyll cell wall thickness ( T CW ). Our study suggests that greater SD ada and SR increased g m by increasing S mes and creating additional parallel pathways for CO 2 diffusion inside mesophyll cells. Thus, SD ada , SR and S mes are important determinants of C 4 ‐ g m and could be the target traits selected or modified for achieving greater g m and TE i in C 4 species.

Pathare, Varsha S.↗
DOE OSTI2025

Investigation of boron powder flow rates on real-time wall

The limit of boron flow rates for real-time conditioning of the first walls has been systematically investigated in the Experimental Advanced Superconducting Tokamak (EAST) with a full metal wall. Initially, solid boron injection demonstrated effective control over carbon impurities and deuterium recycling on the basis of pre-discharge boronization. A minimum flow rate, identified between 1.0 mg/s and 2.0 mg/s, was necessary for actively improving wall conditions under specific plasma operating scenarios, with this effect progressively enhancing as boron flow rates increased. Additionally, a maximum flow rate, estimated between 3.5 mg/s and 8.0 mg/s, was identified for these plasma conditions. When boron flow rates exceeded this maximum, boron-induced fueling effects influenced the plasma line-averaged density, and at excessively high flow rates, plasma disruption was observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Designing a boron nitride polyethylene composite for shielding neutrons

Neutrons are encountered in many different fields, including condensed matter physics, space exploration, nuclear power, and healthcare. Neutrons interacting with a biological target produce secondary charged particles that are damaging to human health. The most effective way to shield neutrons is to slow them to thermal energies and then capture the thermalized neutrons. These factors lead us to consider potential materials solutions for neutron shields that maximize the protection of humans while minimizing the shield mass and adapt well to modern additive manufacturing techniques. Using hexagonal boron nitride (hBN) as a capture medium and high-density polyethylene (HDPE) as a thermalization medium, we aim to design the optimal internal structure of h 10 BN/HDPE composites by minimizing the effective dose, which is a measure of the estimated radiation damage exposure for a human. Through Monte Carlo simulations in Geant4, we find that the optimal structure reduces the effective dose up to a factor of 72 over aluminum (Al) and up to a factor of 4 over HDPE; this is a significant improvement in shielding effectiveness that could dramatically reduce the radiation exposure of occupational workers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Smallest remnants of early matter domination

The evolution of the universe prior to Big Bang Nucleosynthesis could have gone through a phase of early matter domination which enhanced the growth of small-scale dark matter structure. If this period was long enough, self-gravitating objects formed prior to reheating. Here, we study the evolution of these dense early halos through reheating. At the end of early matter domination, the early halos undergo rapid expansion and eventually eject their matter. We find that this process washes out structure on scales much larger than naively expected from the size of the original halos. We compute the density profiles of the early halo remnants and use them to construct late-time power spectra that include these non-linear effects. We evolve the resulting power spectrum to estimate the properties of microhalos that would form after matter-radiation equality. Surprisingly, cosmologies with a short period of early matter domination lead to an earlier onset of microhalo formation compared to those with a long period. In either case, dark matter structure formation begins much earlier than in the standard cosmology, with most dark matter bound in microhalos in the late universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Disruption thermal load mitigation with shattered pellet injection on the Joint European Torus (JET)

Disruption mitigation remains a critical, unresolved challenge for ITER. To aid in addressing this challenge, a shattered pellet injection (SPI) system was installed on JET and experiments conducted at a range of thermal energy fractions and stored energies in excess of 7 MJ. The primary goals of these experiments were to investigate the efficacy of the SPI on JET and the ability of the plasma to assimilate multiple pellets. Single pellet injections produced a saturation in total radiated energy (Wrad) with increasing injected neon content, suggesting total radiation of stored thermal energy. Further increases in injected neon quantities resulted in reduced cooling times and current quench (CQ) durations, indicating higher impurity assimilation. No significant variation in CQ duration or W rad was observed when varying the deuterium content at fixed neon quantities. Additionally, higher assimilation, inferred by shorter CQ durations, was measured when a mechanical punch was used to launch the pellets and this was attributed to a lower pellet velocity leading to higher solid content in the pellet plume and larger fragments penetrating deeper into the plasma. Radiation asymmetries averaged over the cooling time were inferred from Emis3D and ranged from 1.6 to 1.9. Asymmetries averaged over the entire disruption sequence were found to increase at higher thermal energy fractions. The radiated energy fractions decreased with increasing thermal energy fractions but this trend was eliminated when toroidal asymmetries were accounted for with Emis3D. Pure deuterium pellets were able to produce cooling times of up to 75 ms with a gradual loss in thermal stored energy of up to 80%. Experiments with multiple pellet injection indicated W rad can be increased through pellet superposition and density can be increased with an additional D2 injection without a reduction in W rad . KPRAD modelling accurately reproduced the cooling times and the CQ duration at high thermal energies. Assimilation estimates from KPRAD indicated CQ rates scale strongly whilst W rad scales weakly and saturates with assimilated neon content. Comparable W rad can be achieved with lower assimilated neon quantities as longer cooling times are attained. Thus reduced neon content can be preferential in a thermal load mitigation scheme as it may reduce radiation asymmetries and prevent flash melting.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Molecule-surface interaction from van der Waals-corrected semilocal density functionals: The example of thiophene on transition-metal surfaces

Semilocal density functional approximations are widely used. None of them can capture the long-range van der Waals (vdW) attraction between separated subsystems, but they differ remarkably in the extent to which they capture intermediate-range vdW effects responsible for equilibrium bonds between neighboring small closed-shell subsystems. The local density approximation (LDA) often overestimates this effect, while the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) underestimates it. The strongly constrained and appropriately normed (SCAN) meta-GGA often estimates it well. All of these semilocal functionals require an additive nonlocal correction such as the revised Vydrov-Van Voorhis 2010 (rVV10) to capture the long-range part. This work reports adsorption energies and the corresponding geometry of aromatic thiophene (C 4 H 4 S) bound to transition-metal surfaces. The adsorption process requires a genuine interplay of covalent and weak binding and requires a simultaneously accurate description of surface and adsorption energies with the correct prediction of the adsorption site. All these quantities must come from well-balanced short- and long-range correlation effects for a universally applicable method for weak interactions with chemical accuracy. Here, our methods indicate that the correct interplay is not present in any combination of recent meta-GGA's and rVV10. The simple short-range damping of the vdW correction scheme that is practically successful in the combination of GGA's and vdW approximations is less transferable in SCAN+rVV10 or in the revised version, revSCAN+rVV10. In addition, we present accurate random-phase-approximation-quality adsorption energies from a model based on the one of Zaremba and Kohn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Validation of Wichner Predictive Model for Chronic Oxidation by Moisture of Nuclear Graphite

This document reports the results of multi-year extensive research on chronic oxidation of nuclear graphite by traces of moisture that might be present in the helium coolant of high-temperature gas-cooled reactors (HTGRs). The research focuses on four grades of nuclear graphite that are considered potential candidates for component fabrication in advanced nuclear reactors. The four grades are (1) PCEA (medium-fine grain, extruded, manufactured by GrafTech Inc., USA); (2) NBG-17 (medium-fine grain, vibrationally molded, manufactured by SGL Carbon, Germany/France); (3) IG-110 (microfine grain, isostatically molded, manufactured by Toyo Tanso, Japan); and (4) 2114 (microfine grain, isostatically pressed, manufactured by Mersen, USA). The experimental part of this research comprised three objectives that were consistently pursued for each graphite material: (1) accelerated oxidation tests in well-controlled conditions (temperature, gas composition) and analysis of chemical kinetic models and of their parameters for each graphite grade; (2) measurements of water vapor transport rates from H 2 O/He mixtures through thin graphite slabs for estimation of the water effective diffusivity specific to each material; (3) controlled oxidation of graphite specimens in H 2 O/He mixtures and characterization of the effects of oxidation temperatures on the density profiles of oxidized layers. These three groups of experimental results were interpreted using the model developed by Wichner et al. for predicting the penetration depth of oxidants (moisture or oxygen) as a function of oxidation conditions. Basically, starting from the kinetic and transport properties of each graphite grade, the ultimate goal was predicting the stable density profile in the oxidized layer and its dependence on the oxidation temperature. The profiles obtained this way were then compared against direct observations of oxidized layers in specimens obtained in well-controlled oxidation conditions. The general agreement between predictions and observations validated the general applicability of the Wichner model. The results obtained in this project were presented at annual Department of Energy–Advanced Reactor Technologies program reviews and at specialist meetings and conferences and were published in Oak Ridge National Laboratory (ORNL) reports and in professional journals. The results published in 2018 had an immediate impact in 2019, when a group from three US universities used ORNL kinetic data as input information in a complex simulation of chronic oxidation by moisture using the same graphite grades during long-time operation in a modular high-temperature gas-cooled reactor. Their conclusion, in line with the conclusions obtained at the end of this project, was that chronic graphite oxidation by moisture during normal operation conditions of an HTGR will be limited to a thin surface layer of graphite blocks at the reactor bottom, where the highest temperatures are. That will not affect the integrity of graphite components and reactor safety, provided the impurity concentration in the helium coolant is rigorously controlled and maintained at the design specifications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗