Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “apparent size”

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.

179 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Super-resolution surface slope metrology of x-ray mirrors

We present experimental, analytical, and numerical methods developed for reconstruction (deconvolution) of one-dimensional (1D) surface slope profiles over the spatial frequency range where the raw data are significantly perturbed due to the limited resolution of the measurement instrument. We characterize the spatial resolution properties of a profiler with the instrument's transfer function (ITF). To precisely measure the ITF, we apply a recently developed method utilizing test surfaces with 1D linear chirped height profiles of constant slope amplitude. Based on the results of the ITF calibration, we determine parameters of an analytical model for the ITF that is used in the original reconstruction software. Here, we treat surface slope metrology data obtained with the Optical Surface Measuring System (OSMS), using as a sensor an electronic autocollimator (AC) ELCOMAT-3000. The spatial resolution of the OSMS is limited by the size of the AC light-beam-collimating aperture. For the purposes of this investigation, the OSMS is equipped with a circular aperture with a diameter of 2.5 mm. This is a typical arrangement of most AC-based slope profilers developed for surface slope metrology of state-of-the-art x-ray mirrors. Using the example of surface slope metrology of two state-of-the-art elliptically shaped x-ray focusing mirrors, we demonstrate that the developed data reconstruction procedure allows us to significantly improve the accuracy of surface slope metrology with the OSMS over the spatial wavelength range from ~1.6 mm to 7 mm. Thus, the amplitude of the quasi-periodic error characteristic of the deterministic polishing process used appears to be higher by a factor of~2 than is apparent from the rough metrology data. Underestimation of the surface slope errors in this spatial wavelength range can lead to serious errors in the expected performance of x-ray mirrors in synchrotron beamlines, especially at modern light sources utilizing coherent x rays, where the perturbations can lead to increased speckle-like intensity variation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

In Situ Studies of Methanol Decomposition Over Cu(111) and Cu2O/Cu(111): Effects of Reactant Pressure, Surface Morphology, and Hot Spots of Active Sites

The dissociative adsorption of methanol was investigated on Cu(111) and ultrathin Cu 2 O films. We employed synchrotron-based Ambient Pressure X-ray Photoelectron Spectroscopy (AP-XPS) and Scanning Tunneling Microscopy (STM) to study the dynamics of gas–solid interactions, and calculations based on Density Functional Theory (DFT) were used to examine the reaction path. C 1s XPS spectra revealed that methanol underwent dissociative adsorption on plain Cu(111) to form methoxy (CH 3 O), formaldehyde (H 2 CO), and formate (HCOO) at a pressure range of 0.5–10 mTorr, with these species remaining on the surface after evacuation. This was accompanied by the appearance of a low coverage (~0.05 ML) of O ads in the O 1s which can be considered a highly active site for methanol activation. The high activity is apparent by a coverage of 0.8 ML of methoxy at room temperature. STM was unable to image these species at room temperature as they were highly mobile on metallic copper. In contrast, for CH 3 OH on Cu 2 O/Cu(111), STM showed clear hot spots for reaction and a complex array of adsorption structures. On the oxide substrate, there was decomposition of methanol to H 2 CO, CH 3 O, HCOO, and hydrocarbon species (CH x ) due to the subsequent interactions of methanol with lattice oxygen. Cu(111) remained entirely saturated with decomposition products under 10 mTorr of methanol (θ ≈ 0.97 ML), whereas the Cu 2 O overlayer was saturated at a much lower coverage (θ ≈ 0.30 ML). STM revealed rows and step edges of Cu 2 O decorated with decomposition products and metallic Cu islands ~5 nm in size. The difference in activity between Cu(111) and Cu 2 O/Cu(111) is attributed to the significant amount of O present on the oxide surface. Additionally, Density Functional theory (DFT) calculations described the XPS measurements well, showing a likely methanol dissociation to *CH 3 O and therefore a surface reduction. More importantly, the DFT results revealed that it was the chemisorbed oxygen on Cu 2 O/Cu(111) which oxidized the dissociated *CH 3 O to *HCOO and eventually CO 2 , while the reaction only involving upper oxygen on the Cu 2 O hexagonal ring led to the formation of H 2 CO.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Barrierless nucleation in glassy precursors drives zeolite formation

Zeolites are crystalline, microporous silicates widely used in catalysis and separations, yet the molecular mechanisms of their formation remain unresolved. Experiments indicate that hydrothermal synthesis of silica zeolites from clear solution proceeds through amorphous nanoaggregates that gradually develop zeolite order in an apparently continuous amorphous to crystal transformation. Here, we combine molecular simulations with advanced algorithms that identify zeolite order and computer vision to elucidate the pathway from clear solution to zeolite nanocrystal. We show that at conditions of hydrothermal synthesis of silica zeolites, the transformation of precursor aggregate into zeolite is not limited by nucleation barriers but by the slow dynamics of reorganization in the glassy precursor matrix. The negligible nucleation barriers result in spinodal-like crystallization that leads to a gradual formation of a mosaic of small crystallites that explain the seemingly continuous character of zeolite crystallization and the catalytic activity of X-ray amorphous, protozeolites and embryonic zeolites. We find that zeolite-like porosity and short-range order emerge early within glassy precursors, well before crystallinity is detected in transmission electron microscopy (TEM) images or X-ray diffraction. The nanoaggregate’s temperature-size phase diagram reveals a convergence of the zeolite–amorphous equilibrium and maximum crystallization rate at ~3 nm diameter nanoparticle diameters and ~200 °C. This convergence signals the termination of the first-order amorphous-to-zeolite transition. Our results provide a unifying framework for understanding nucleation of silica zeolites from solution and suggest that barrierless nucleation may govern the formation of other nanoparticle systems, including minerals and oxides synthesized far below their bulk melting points.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Compressed Expanded Natural Graphite (CENG) Processing for PCM Composites

The use of phase change materials (PCMs) in thermal energy storage applications has received considerable attention in recent decades. Organic PCMs are popular due to their high latent heat of fusion, noncorrosive properties, and relative stability over many charge and discharge cycles. A primary limitation of these materials is their low thermal conductivity. This has led researchers to develop various methods to increase thermal conductivity by seeding PCM with or impregnating them into more conductive materials. One method is to impregnate PCMs into compressed expanded natural graphite (CENG) matrices, which can improve thermal conductivity by a factor of 100. CENG matrices have received particular interest due to their low cost, high porosity, small (nano/micro) pore size, high pore density, high thermal conductivity, and ability to be compressed into many geometries. PCM/CENG matrix composites have been extensively studied; however, the effect that CENG processing has on PCM saturation and the overall matrix thermal conductivity has not been well investigated. This processing includes four major steps including graphite intercalation, thermal shock, compression, and PCM saturation. Intercalation involves soaking graphite flakes in sulfuric and/or nitric acid to intercalate the acids between the graphene layers. The graphite flakes are then subjected to a high-temperature thermal shock, during which the intercalated acid is gasified rapidly, pushing the graphene layers apart, resulting in accordion-shaped graphite "worms". The "worms" are then compacted to a desired bulk density and then soaked with molten PCM until fully saturated. The properties of the produced CENG matrix, and its ability to allow PCM permeation, are sensitive to the processing parameters, namely, the thermal shock temperature and exposure time, as well as the matrix apparent density or porosity. Here, we study the effect of the thermal shock conditions necessary to expand intercalated graphite flakes on PCM saturation and the expanded graphite's thermal conductivity and morphology. We found that the thermal shock temperature exhibits the greatest influence. At greater shock temperatures, SEM images showed that expanded graphite worms exhibited greater density of pores, thus increasing total surface area within the matrices. Increasing thermal shock temperature yielded greater overall PCM saturation, as well as an increased rate of saturation. Improvements in PCM saturation rate and overall saturation are obtained as the shock temperature is increased. Longer exposure to thermal shock also improves initial saturation rates and is beneficial if a shortened impregnation time is needed. Thermal shock conditions did not impact thermal conductivity; however, conductivity was largely affected by matrix porosity. A local maximum in axial thermal conductivity was observed at around 83% porosity, which is similar to that observed in previous studies.

compressed expanded natural graphite↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
DOE OSTI2023

Modern temperatures in central–north Greenland warmest in past millennium

The Greenland Ice Sheet has a central role in the global climate system owing to its size, radiative effects and freshwater storage, and as a potential tipping point. Weather stations show that the coastal regions are warming, but the imprint of global warming in the central part of the ice sheet is unclear, owing to missing long-term observations. Current ice-core-based temperature reconstructions are ambiguous with respect to isolating global warming signatures from natural variability, because they are too noisy and do not include the most recent decades. By systematically redrilling ice cores, we created a high-quality reconstruction of central and north Greenland temperatures from ad 1000 until 2011. Here we show that the warming in the recent reconstructed decade exceeds the range of the pre-industrial temperature variability in the past millennium with virtual certainty (P < 0.001) and is on average 1.5 ± 0.4 degrees Celsius (1 standard error) warmer than the twentieth century. Our findings suggest that these exceptional temperatures arise from the superposition of natural variability with a long-term warming trend, apparent since ad 1800. The disproportionate warming is accompanied by enhanced Greenland meltwater run-off, implying that anthropogenic influence has also arrived in central and north Greenland, which might further accelerate the overall Greenland mass loss.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Trading water for carbon in the future: Effects of elevated CO 2 and warming on leaf hydraulic traits in a semiarid grassland

Abstract The effects of climate change on plants and ecosystems are mediated by plant hydraulic traits, including interspecific and intraspecific variability of trait phenotypes. Yet, integrative and realistic studies of hydraulic traits and climate change are rare. In a semiarid grassland, we assessed the response of several plant hydraulic traits to elevated CO 2 (+200 ppm) and warming (+1.5 to 3°C; day to night). For leaves of five dominant species (three graminoids and two forbs), and in replicated plots exposed to 7 years of elevated CO 2 , warming, or ambient climate, we measured: stomatal density and size, xylem vessel size, turgor loss point, and water potential (pre‐dawn). Interspecific differences in hydraulic traits were larger than intraspecific shifts induced by elevated CO 2 and/or warming. Effects of elevated CO 2 were greater than effects of warming, and interactions between treatments were weak or not detected. The forbs showed little phenotypic plasticity. The graminoids had leaf water potentials and turgor loss points that were 10% to 50% less negative under elevated CO 2 ; thus, climate change might cause these species to adjust their drought resistance strategy away from tolerance and toward avoidance. The C4 grass also reduced allocation of leaf area to stomata under elevated CO 2 , which helps explain observations of higher soil moisture. The shifts in hydraulic traits under elevated CO 2 were not, however, simply due to higher soil moisture. Integration of our results with others' indicates that common species in this grassland are more likely to adjust stomatal aperture in response to near‐term climate change, rather than anatomical traits; this contrasts with apparent effects of changing CO 2 on plant anatomy over evolutionary time. Future studies should assess how plant responses to drought may be constrained by the apparent shift from tolerance (via low turgor loss point) to avoidance (via stomatal regulation and/or access to deeper soil moisture).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

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

Tandem Catalytic Antioxidant Nanoparticles Comprising Cerium Carbonate and Photoactive Metal Oxides

While photoactive metal oxides such as TiO 2 find widespread use in paints and coatings as well as cosmetics and suncare products, they also generate reactive oxygen species (ROS), which degrade materials and are associated with human-health pathologies. Here, we demonstrate a robust and potent catalytic antioxidant consisting of earth-abundant cerium carbonate nanoparticles and micron-size Ce 2 (CO 3 ) 3 . 8H 2 O, which are characterized by powder X-ray diffraction and scanning nanobeam electron diffraction. When dispersed with photoactive metal oxides, these cerium carbonate catalysts decrease the photodecomposition rate of organic dyes and commercial pigment colorants in aqueous media by up to 820-fold, as well as in acrylic coatings. X-ray photoelectron spectroscopy and kinetic experiments support the same tandem catalysis mechanism of photoprotection when using both micron-size Ce 2 (CO 3 ) 3 . 8H 2 O and cerium carbonate nanoparticles. This mechanism involves ROS disproportionation (catalyzed by cerium carbonate) and H 2 O 2 decomposition (partially catalyzed by TiO 2 ) pathways, both of which cerium carbonate also catalyzes on its own, crudely mimicking the function of the cascade system of superoxide dismutase and catalase enzymes. When cerium carbonate nanoparticles were dispersed at 2 wt % in polymethylmethacrylate, the transparency of the polymer film was preserved and the photo-oxidative degradation of the polymer was prevented following UV irradiation at 254 nm, which otherwise resulted in the loss of optical properties and hydroxylation as characterized by ATR-FTIR spectroscopy, in the control polymer lacking cerium carbonate. Similar observations were made regarding color preservation in paint films comprising dye and insoluble commercial colorant pigments. The material chemistry associated with this photoprotection catalysis is subtle and emphasizes the importance of both Ce(III) and carbonate together, as both CePO 4 and Na 2 CO 3 are inactive. In conclusion, this emphasis is also apparent in comparisons of photoprotection catalysis with previously reported cerium oxide nanoparticles, which are significantly less active compared with Ce 2 (CO 3 ) 3 . 8H 2 O under the same conditions.

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

Intrinsic Kinetic Limitations in Substituted Lithium-Layered Transition-Metal Oxide Electrodes

Substituted Li-layered transition-metal oxide (LTMO) electrodes such as Li x Ni y Mn z Co 1-y-z O 2 (NMC) and Li x Ni y Co 1-y-z Al z O 2 (NCA) show reduced first cycle Coulombic efficiency (90-87% under standard cycling conditions) in comparison with the archetypal Li x CoO 2 (LCO; similar to 98% efficiency). Focusing on Li x Ni 0.8 Co 0.15 Al 0.0.5 O 2 as a model compound, we use operando synchrotron X-ray diffraction (XRD) and nuclear magnetic resonance (NMR) spectroscopy to demonstrate that the apparent first-cycle capacity loss is a kinetic effect linked to limited Li mobility at x > 0.88, with near full capacity recovered during a potentiostatic hold following the galvanostatic charge- discharge cycle. This kinetic capacity loss, unlike many capacity losses in LTMOs, is independent of the cutoff voltage during delithiation and it is a reversible process. The kinetic limitation manifests not only as the kinetic capacity loss during discharge but as a subtle bimodal compositional distribution early in charge and, also, a dramatic increase of the charge-discharge voltage hysteresis at x > 0.88. 7 Li NMR measurements indicate that the kinetic limitation reflects limited Li transport at x > 0.86. Electrochemical measurements on a wider range of LTMOs including Li x (Ni,Fe) y Co 1-y O 2 suggest that 5% substitution is sufficient to induce the kinetic limitation and that the effect is not limited to Ni substitution. In this paper, we outline how, in addition to a reduction in the number of Li vacancies and shrinkage of the Li-layer size, the intrinsic charge storage mechanism (two-phase vs solid-solution) and localization of charge give rise to additional kinetic barriers in NCA and nonmetallic LTMOs in general.

25 ENERGY STORAGE↗
DOE OSTI2025

Evolving Dark Sector and the Dark Dimension Scenario

String theory naturally leads to the expectation that dark energy is not stable, and may be evolving as captured by the Swampland de Sitter conjectures. Moreover, motivated by the distance conjecture, a unification of dark sector has been proposed, where the smallness of dark energy leads to one extra dimension of micron size with dark matter being the Kaluza--Klein graviton excitations in this extra dimension. We consider the natural possibility that the radius of the dark dimension varies as the dark energy decreases, leading to the variation of the dark matter mass. This correlates the decrease of the dark energy with the variation of the dark matter mass as they depend on the variations of a scalar field $ϕ$ controlling the radius of the extra dimension. A simple realization of this idea for small range of $ϕ$ is adequately captured by choosing a potential which is locally of the form $V=V_0\ {\rm exp}(-cϕ)$ and dark matter mass $m_{\rm DM}=m_0\ {\rm exp}(-c' ϕ)$ where the sign of $ϕ$ is chosen such that $c'\geq 0$ while we have two choices for the sign of $c$ depending on whether the dark dimension expands or shrinks when the dark energy dominates. We find excellent agreement with recent experimental data from DESI DR2 combined with SN measurements (from DES, Union3 or Pantheon+) and reproduces the same significance as CPL parametrization with the added benefit of providing a natural explanation for the apparent phantom behavior ($w<-1$) reported by DESI and DES based on a physical model. DESI and SN datasets independently favor non-zero values of $c'$ and $c$, respectively, both lying within the expected $\mathcal{O}(1)$ range suggested by the Swampland criteria. Moreover, our best fit value $c'\simeq 0.05 \pm 0.01$ is remarkably consistent with the experimental upper bound of $c'\lesssim 0.2$ demanded by the lack of detection of fifth force in the dark sector.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2023

Heavy Versus Light Lanthanide Selectivity for Graphene Oxide Films is Concentration Dependent

Rare earths are important materials in various technologies such as catalysis and optoelectronics. Graphene oxide (GO) is a promising material for separation applications, including the isolation of lanthanides from complex mixtures. Previous works using fatty acid monolayers have demonstrated preferential heavy versus light lanthanide adsorption, which has been attributed to differences in lanthanide ion size. In this work, we used interfacial X-ray fluorescence measurements to reveal that GO thin films at the air/water interface have no lanthanide selectivity for dilute subphases. However, at high subphase concentrations, ~8 times more Lu is adsorbed than La. By comparing the GO results with an ideal monolayer with a carboxylic acid headgroup, arachidic acid (AA), we demonstrate that the number of Lu ions adsorbed to GO is significantly higher than the number expected to compensate for the surface charge. Vibrational sum frequency generation (SFG) spectroscopy results on both GO thin films and AA monolayers reveal a red-shifted SFG signal in the OH region, which we attribute to partial dehydration of the adsorbed ions and carboxylic acid headgroups. Liquid surface X-ray reflectivity data show that the GO thin film structure does not significantly change between the very dilute and concentrated subphases. We speculate that the functional groups of both GO and AA facilitate cation dehydration, which is essential for ion adsorption. Heavy lanthanide Lu has stronger ion–ion correlations that can overcome the electrostatic repulsion between cations at higher concentrations compared to light lanthanide La, meaning GO and AA can exhibit apparent overcharge with Lu. Lastly, the layered structure of the GO films and reactive chemical nature of GO itself can accommodate ion adsorption.

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

Cavitation-resistant intergranular precipitates enhance creep performance of θ'-strengthened Al-Cu based alloys

Tensile and compressive creep properties of a quaternary Al-Cu-Mn-Zr (ACMZ) alloy and its commercial counterpart (Al-Cu-Mn-Zr with Ni, Co and Sb additions, RR350) are investigated at 300°C. At low stresses up to 30 MPa where diffusional creep dominates, creep resistance is the same in tension and compression and RR350 deforms more slowly than ACMZ, consistent with RR350 alloy's larger linear fraction of intergranular precipitates (Al7Cu2(NiFe) and Al9FeNi for RR350 vs. θ-Al2Cu for ACMZ) and a reduced fraction of precipitate-free zones near grain boundaries. At stresses between 30 and 80 MPa, dislocation creep with a stress exponent n ~ 3 becomes rate-limiting in compression, which is expected to be controlled by θ' precipitates within the grain bulk. By contrast, in tension, enhanced creep rate and higher apparent stress exponents are measured, consistent with cavitation at intergranular precipitates becoming increasingly dominant as the stress increases. In the dislocation creep regime, RR350 alloy is again more creep resistant than ACMZ alloy, which is related to three mechanisms (i) a reduced fraction of softer precipitate-free zones, (ii) more effective load transfer to intergranular precipitates, and (iii) reduced cavitation. A model for cavitation is applied to calculate tensile creep rates from compressive creep rates and the model successfully predicts the improved tensile creep resistance of the RR350 alloy. Overall, the present investigation underscores the importance of intergranular grain boundary precipitates, in addition to strengthening θ' precipitates, in enhancing the creep resistance of Al-Cu alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Demonstrate new plasticity models for doped UO 2 that capture dislocation mechanisms

In light water reactors, fuel vendors are investigating the use of dopants to modify the properties of UO 2 pellets, with the goal of improving pellet-cladding mechanical interactions during operation. Dopants are expected to ‘soften’ the pellets; that is, the doped pellets have higher plastic deformation than conventional UO 2 . This leads to a reduction in the severity of mechanical pellet-cladding interactions, helping to reduce the hoop strain on the cladding. By minimizing the strain exerted by the pellet on the cladding, it is anticipated that cladding performance under accident conditions can be enhanced (i.e., lowering the risk of burst during a LOCA). Dopants such as chromium (Cr) promote grain growth during pellet fabrication, leading to larger grains; therefore, understanding the link between chemistry, microstructure and mechanical deformation (enhanced creep rates) behavior of UO 2 is critical to helping operators further substantiate the benefits of doping UO 2 . Historically, the nuclear energy industry has relied on empirical models to make assessments of performance. Compared to empirical models, mechanistic physics-based models provide benefits, such as, fewer data points for validation and better extrapolation where experimental data is scarce or non-existent. In this report, Bayesian inference techniques have been applied to a previously developed lower length-scale-informed diffusional creep model. The objective is to i) infer lower-length-scale parameter distributions from available experiment and then ii) determine the uncertainties in the measurable quantity (in this case creep rates) after propagating the inferred lower length scale parameter uncertainties. The approach requires many evaluations of the model, which becomes computationally insurmountable; therefore, a neural-network model is trained to data obtained by sampling the full model over the most important parameters. This neural-network is then used in the Bayesian inference approach to determine probability distributions in the parameter values that represent the uncertainty in the model given what is known from the experiments (posterior). A significant reduction compared to conservative initial (prior) uncertainties is achieved through inference against the experimental data, demonstrating the efficacy of this approach. Furthermore, by accounting for uncertainties in the experimental conditions and sample non-stoichiometry, it is possible to resolve apparent discrepancies in experimental measurements within a self-consistent grain boundary (Coble) creep model that is sensitive to chemistry. This work has been written up and submitted to Nuclear Technology for a special issue on accelerated fuel qualification (AFQ). This uncertainty quantification (UQ) work not only improves the diffusional model, while accounting for uncertainty, but also establishes a framework which can readily be applied to the mechanistic models of dislocation deformation developed in this study. The most likely values from the Bayesian analysis are incorporated into our UO 2 diffusional creep model and a lower length scale-informed irradiation UO 2 creep mechanistic model to generate a dataset. This dataset has been provided to our INL collaborators for training an artificial neural network surrogate model, which will be implemented in the BISON fuel performance code to assess how the results differ from those currently obtained using a fully empirical model and that of using the nominal (uncalibrated) atomic scale parameters in our mechanistic model. Plastic deformation (creep and glide) in UO 2 is a complex phenomenon, governed by multiple underlying processes such as local defect concentrations, applied stresses, and microstructural characteristics. Consequently, there is a need for a meso-scale model with polycrystalline resolution capable of extrapolating to large grain sizes applicable to doped UO 2 , where data is limited and the model can help bridge the knowledge gap. By integrating atomistic data into the polycrystal LApx code, it becomes possible to predict dislocation climb and glide plasticity that simple analytical models cannot accurately represent. The application of atomic-scale data within LApx demonstrated the importance of climb and glide mechanisms in reproducing high-stress UO 2 behavior. Behaviors such as this are crucial to capture and implement in BISON, as parts of the fuel pellet can reach temperatures where glide can occur before pellet cracking. This model which captures dislocation based mechanisms for UO 2 is then used to stand up the doped model accounting for larger grain sizes. It was found that larger grain sizes can lead to enhanced deformation rates in the glide regime, and therefore can help with the pellet cladding mechanical interaction. Therefore if the fuel pellet reaches conditions (stress/temperature) where glide is active, the enhanced creep rates for larger grains in the glide regime (doped UO 2 ) can help with pellet cladding mechanical interactions. Plastic deformation in UO 2 involves multiple mechanisms, including diffusional creep, dislocation climb, and glide. This milestone contains two parts: (1) UQ of a pre-existing lower length scale informed mechanistic diffusional creep model, and (2) development of a new LApx based model for dislocation-mediated creep mechanisms in UO 2 , with application to large-grain doped UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Petrophysical and geomechanical properties of gas hydrate-bearing sediments recovered from Alaska North Slope 2018 Hydrate-01 Stratigraphic Test Well

Knowledge of petrophysical and geomechanical properties of gas hydrate-bearing sediments are essential for predicting reservoir response due to gas production from gas hydrate reservoirs. That information will be critical parameters for designing production well completion such as specification of depressurization pump, water storage tank, and mesh size of sand screen. In December 2018, Stratigraphic Test Well Hydrate-01 was drilled in the western part of the Prudhoe Bay Unit, Alaska North Slope as part of the technical planning effort for a future long-term production test by collaborative team of DOE/NETL, USGS, and MH21-S (Boswell et al., 2020, Collett et al., 2020, Okinaka et al., 2020). Data set of logging-while-drilling (LWD) were acquired (Haines et al., 2020, Suzuki et al., 2019) and core sampling depth was determined on-site.Side-wall pressure coring was conducted to recover gas hydrate-bearing sediments from two reservoir sections named Unit-B and Unit-D. A total of 34 cores were successfully recovered by 5 runs of a wire-line deployed pressure corer (CoreVault® System - Halliburton). Core analysis scheme of this project are shown in Figure 1. All cores were quenched in liquid nitrogen while at high pressure in the laboratory of Stratum Reservoir, LLC. at Anchorage (Figure 1, a)). And the cores were removed from the pressure corer autoclave with temperature support of dry ice and stored under liquid nitrogen at atmospheric pressure. 19 damaged cores were processed to index property measurements directly such as grain size, grain density. 4 of another 17 cores were depressurized and trimmed for making a plug to analyze petrophysical properties of host sediments. Unsteady-state permeability test was conducted to obtain relative water permeability to gas and core scale NMR T2 distribution measurement was performed for evaluating pore size distribution at Houston (Figure 1, b)). Remained high quality 13 cores were preserved with gas hydrate for advanced laboratory analysis. National Institute of Advanced Industrial Science and Technology, as a part of the Japanese National Hydrate Research Program (MH21-S, funded by Ministry of Economy, Trade and Industry), received the samples at Sapporo, Japan for advanced core analysis. High-resolution X-ray CT was used to analyze the quality of the samples, which showed undisturbed lithological layers. Cores were lathed into cylindrical shape and distributed for multi property measurements (Figure 1, c)).<p>As a result, sediment from Unit-D is categorized as silty sand at ~37% porosity with ~80% gas hydrate saturation. An average hydration number n = 6.16 was measured by Raman spectroscopy. An average intrinsic permeability of ~400 mD and in situ effective permeability (with hydrate) on the order of ~10 mD. The Unit B recovered cores consisted of well sorted sand at ~40% porosity with ~95% gas hydrate saturation. An average intrinsic permeability of ~1 Darcy and in situ effective permeability on the order of ~30 mD was measured for the Unit B cores. There was a small permeability reduction due to porosity loss with increasing effective stress that simulated consolidation behavior along with depressurization in the highly permeable sandy sediment. The apparent minimum change in porosity and permeability may be caused by the low compressibility of quartz sand grains in the recovered cores. XRD and thermal conductivity analysis also suggested high quartz content. Triaxial compression tests established internal friction angles based on the Mohr-Coulomb's failure criterion, which are 40° for hydrate-bearing sediment and 29.8° for hydrate free sediment.</p>

Yoneda, Jun↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Controlling Structural, Electronic, and Energy Flow Dynamics of Catalytic Processes through Tailored Nanostructures

MoS 2 (molybdenum disulfide) is a highly-versatile catalyst material for support of numerous reactions from hydrodesulfurization and denitrogenation to the focus of this renewal proposal: hydrogenation of CO/CO 2 towards (higher) alcohols. At the same time, MoS 2 is a non-toxic, environmentally-benign and rather inert material- which under ambient conditions for some time has served as a lubricant and, more recently, as next-generation electronic material. The apparent contrast between inertness and stability in ambient, and catalytic activity under reactive conditions is puzzling and calls for a synergistic theoretical and experimental investigation with the long-term objective of enabling the rational design of MoS 2 -based catalysts for alcohol-formation reactions by providing a microscopic understanding of the environmental factors that determines site activity and selectivity. Our research project seeks answers to the questions (a) what conformation does MoS 2 adopt under reaction conditions (as opposed to that under ultrahigh vacuum and low temperatures)?; (b) what reaction pathways exist on such a material?; (c) how can the local environment of the active sites be manipulated so as to make MoS 2 an efficient catalyst for production of higher alcohol from syngas? In particular, research strategies will explore how the basal plane composed of sulfur atoms can be activated so as to exhibit a reactivity of its own, by addressing three research targets and building on extensive preliminary and enabling work: (1) vacancies and vacancy aggregates on the basal plane; (2) non-local catalyst transformation through alkali doping, hydrogenation and phase transition; (3) fabrication of a metal-nanoparticle-activated MoS 2 system, in which particle anchoring, reactive sites and pathways as well as selectivity are controlled by design, as an example of predictive development of a catalyst material,. All strategies are directed to improve the efficacy of the key reactive sites and selectivity of chemical pathways by design, to replace the inefficient methodology of trial and error in catalyst development. This research project represents a synergistic combination of computational guidance, foundational surface-science-based experiments and validation under reactive conditions that aims at transformative new insights into the working of MoS 2 -based hydrogenation-catalysts, a topic squarely at the center of the interest of DOE BES. Guided and led by Talat Rahman, a computational physicist, this project will apply density functional theory to understand structure, reaction pathways and chemical potentials associated with MoS 2 -based CO/CO 2 hydrogenation, augmented by kinetic Monte Carlo methods for reaction rates and prefactors as well as ab-initio molecular dynamics for evaluation of thermal stability. Complementary experimental input and validation will originate from co-PI Ludwig Bartels, a physical chemist and materials scientist, whose group focuses on local imaging and preparation of MoS 2 materials, from co-PI Peter Dowben, an experimental physicist, whose group is expert in the spectroscopy of occupied and unoccupied electronic states, and from senior collaborator Michael White of Brookhaven National Laboratory, whose group generates high-resolution electronic and activity information on size-selected well-defined metal chalcogenide clusters. This collaborative effort will enable a comprehensive understanding of the correlation of structural integrity and catalytic activity of MoS 2 in forms ranging from extended films to individual particles with known geometries and binding sites. Alcohol formation from syngas is a rapidly emerging application that has great potential through facile, economic, and decentralized biomass gasification. CO 2 activation is one of the most pressing concerns of our time: increasing CO 2 levels in the atmosphere change the climate and expose the globe to environmental transformations with the potential for enormous economic and societal impact. We will investigate CO/CO 2 hydrogenation via an interdisciplinary research collaboration with established synergy – one that, in accordance with the mission of the DOE, involves accredited Hispanic-Serving Institutions and that, through student exchange with international collaborators and National Labs, will directly benefit a broad spectrum of communities and generate human resources in sciences essential to their future.

2D materials↗