Search NASA⌕ Search

SEARCH · Search NASA

Results for “Material point method”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Biofouling and Corrosion Study for a Novel Linear Guided Wave Energy Converter (Post Access Report)

The overall objective of this project was to examine the reliability and performance of antibiofouling coatings used for a wave energy converter (WEC) developed by E-Wave Technologies. The particular coatings were selected for their low toxicity and potential compatibility with aquaculture. The aim of this work was to 1) test coating solutions to prevent biofouling growth and saltwater corrosion on the static (paddle and attachment frame surface) components of the WEC that are submerged, 2) determine adhesion of the coatings to system components, and 3) assess the ease and effectiveness of biofouling cleaning to insure long term performance of the system. An analysis of commercial coatings was performed using methods to examine the prevention of biofouling and coating adhesion properties on two key materials of the WEC, which were 316L low carbon marine grade stainless steel (SS) and Ultra High Molecular Weight Polyethylene (PE). Three marine antifouling paints were selected based on their unique properties to test how different paint styles perform on different materials. The selected paints were ePaint Ecominder self-polishing paint with Zinc Omadine for slime control, Pettit ECO HRT Copper-Free ablative antifouling with Econea biocide, and Intersleek 1100SR foul release. Pacific Northwest National Laboratory (PNNL) prepared PE and SS substrates coated with the three paints and compared the performance against uncoated substrates when submerged in raw seawater for 3-, 6-, and 9-month (m) time periods. Results in adhesion testing indicated that Pettit and ePaint materials clearly bonded strongly to SS, but did not bond comparably well to PE. It was noted during adhesion testing that the Intersleek surfaces were especially difficult to test as the paint highly resists bonding to the epoxy adhesives used with the adherence testing platform. The wear rate of the coatings was not measured under this study; however, based on adhesion testing, coatings in the sliding regions of the Ewave device are expected to wear rapidly. Sandia National Laboratories (SNL) evaluated the adhesion of three different paints to PE and SS substrates which were exposed to a marine environment for time intervals of 0, 3, 6, and 9 months. From qualitative visual analysis of the 3 in2 coupons when pulled from the tank, the 3 in2 coupons generally only appeared to have biofouling consisting of filamentous algae or diatoms, which all have relatively low mass and can be easily wiped from the surface of coupons. Qualitative visual analysis indicated that ECO HRT and Unpainted were consistently worse than Intersleek and Ecominder at all time points. Results provide insight to aid with down-selection of commercial coatings under static conditions to support reliability of the WEC and potential maintenance schedules. This investigation was conducted using small coupon samples suspended in seawater and the development of testing rigs for dynamic component level testing is needed for future work. In addition, the potential toxicity of these commercial coatings on aquaculture has not been determined by this study. One recommendation is to conduct toxicity investigations at the Environmental Toxicity Laboratory at Oak Ridge National Laboratory.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

On-demand formation of Lewis bases for efficient and stable perovskite solar cells

In the fabrication of FAPbI 3 -based perovskite solar cells, Lewis bases play a crucial role in facilitating the formation of the desired photovoltaic α-phase. However, an inherent contradiction exists in their role: they must strongly bind to stabilize the intermediate δ-phase, yet weakly bind for rapid removal to enable phase transition and grain growth. To resolve this conflict, we introduced an on-demand Lewis base molecule formation strategy. This approach utilized Lewis-acid-containing organic salts as synthesis additives, which deprotonated to generate Lewis bases precisely when needed and could be reprotonated back to salts for rapid removal once their role is fulfilled. This method promoted the optimal crystallization of α-phase FAPbI 3 perovskite films, ensuring the uniform vertical distribution of A-site cations, larger grain sizes and fewer voids at buried interfaces. Perovskite solar cells incorporating semicarbazide hydrochloride achieved an efficiency of 26.1%, with a National Renewable Energy Laboratory-certified quasi-steady-state efficiency of 25.33%. These cells retained 96% of their initial efficiency after 1,000 h of operation at 85 °C under maximum power point tracking. Additionally, mini-modules with an aperture area of 11.52 cm 2 reached an efficiency of 21.47%. This strategy is broadly applicable to all Lewis-acid-containing organic salts with low acid dissociation constants and offers a universal approach to enhance the performance of perovskite solar cells and modules.

14 SOLAR ENERGY↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2025

Calibration of RAFM Micromechanical Model for Creep Using Bayesian Optimization for Functional Output

A Bayesian optimization procedure is presented for calibrating a multimechanism micromechanical model for creep to experimental data of F82H steel. Reduced activation ferritic martensitic (RAFM) steels based on Fe(8–9)%Cr are the most promising candidates for some fusion reactor structures. Although there are indications that RAFM steel could be viable for fusion applications at temperatures up to 600°C, the maximum operating temperature will be determined by the creep properties of the structural material and the breeder material compatibility with the structural material. Due to the relative paucity of available creep data on F82H steel compared to other alloys such as Grade 91 steel, micromechanical models are sought for simulating creep based on relevant deformation mechanisms. As a point of departure, this work recalibrates a model form that was previously proposed for Grade 91 steel to match creep curves for F82H steel. Due to the large number of parameters (9) and cost of the nonlinear simulations, an automated approach for tuning the parameters is pursued using a recently developed Bayesian optimization for functional output (BOFO) framework (Huang et al., 2021, “Bayesian optimization of functional output in inverse problems,” Optim. Eng., 22, pp. 2553–2574). Incorporating extensions such as batch sequencing and weighted experimental load cases into BOFO, a reasonably small error between experimental and simulated creep curves at two load levels is achieved in a reasonable number of iterations. In conclusion, validation with an additional creep curve provides confidence in the fitted parameters obtained from the automated calibration procedure to describe the creep behavior of F82H steel.

42 ENGINEERING↗
DOE OSTI2024

Thermal and electric field driven rf breakdown precursor formation on metal surfaces

The phenomenon of electric breakdown poses serious challenges to the design of devices that operate in high electric field environments. Experimental evidence points toward breakdown events that are accompanied by elevated temperatures and dark current spikes, which is attributed to high-asperity nanostructure formation that enhances the local electric field and triggers a runaway process. However, the exact mechanistic origin of such nanostructures under typical macroscopic operational conditions of electric field and magnetic-field-mediated heating remains poorly understood. In this work, we simulate the evolution of a copper surface under the combined action of the electric fields and elevated temperatures. Using a mesoscale curvature-driven surface evolution model, we show how a copper surface can undergo a type of dynamical instability that naturally leads to the formation of sharp asperities in realistic experimental conditions. Exploring the combined effect of fields and temperature rise, we identify the critical regimes that allow for the formation of breakdown precursors. The results show that thermoelastic stresses, while not essential, can significantly lower the critical electric field required for runaway surface instability, which is consistent with experimental observations that thermal effects can increase breakdown rates. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Subsurface mechanical damage of fused silica glass during grinding by various sub-aperture tools with and without ultrasonics

The subsurface mechanical damage (SSD) depth after grinding fused silica glass with a comprehensive set of sub-aperture fixed abrasive grinding tools [cup, wheel, belt, pad, and rotary face mill (with and without ultrasonics)] and process parameters has been statistically measured using the taper wedge technique and evaluated. Consistent with a previously reported grinding model [J. Non-Cryst. Solids 352, 5601–5617 (2006) Crossref , Materials Science and Technology of Optical Fabrication (Wiley & Sons, 2018)], based on the sliding indentation fracture by sliding particles or asperities where the normal load per particle determines the depth of the fracture (and ultimately the overall SSD depth), the dominant factor controlling SSD depth was found to be the abrasive size regardless of the tool type and process conditions. Compared to full aperture grinding methods, the overall SSD depth was higher using the sub-aperture tools, likely due to the higher effective pressure and higher load per particle distribution. Here, in addition to abrasive size, a significant reduction in SSD depth was achieved by: (1) reducing the load distribution on the abrasive particles via increase in contact area and/or decrease in mechanical loading; (2) using a more compliant host tool medium; and (3) in what we believe is a more novel way, using ultrasonics. Combining low abrasive size, larger contact area, and a compliant host, the 6 µm diamond in a resin matrix (Trizact) on a foam pad led to very low SSD depth (~ 4.6 µm), relatively fast grinding rate (186 mm 3 /h), and little or no figure degradation. This grinding tool/process is an attractive choice for final grind, resulting in significantly reduced polish out (i.e., “grey out”) times. With the rotary face mill tool, the use of ultrasonics consistently led to a SSD depth reduction (ranging from 17%–34%). A new fracture mechanics-based model, to the best of our knowledge, where the relevant normal load is parallel to the feed direction, has been developed to explain how ultrasonics leads to lower SSD depth. The key factors, supported by finite element stress analysis and load measurements, are (1) the initiation of fractures at higher z heights during the tool’s ultrasonic vertical oscillations, thus propagating less deep into workpiece; (2) reduction in load (and therefore reduction in fracture propagation distance) due to smaller tool-workpiece feed direction contact area (again caused by higher heights relative to depth of cut); (3) upward movement of the tool during oscillation leads to fracturing toward the surface instead into the depth; and finally (4) at tool’s lowest point of oscillation cycle, there may not be enough time for the fracture to propagate to its full length.

Optics and optical instruments↗
DOE OSTI2025

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Spin Waves in Dirac Semimetal Ca 0.6⁢ Sr 0.4⁢ MnSb 2 Investigated with Neutrons by the Diffraction Method

The tunability of Dirac semimetals with antiferromagnetic Mn layers is of great interest. The observed sign change of interlayer magnetic coupling between CaMnBi 2 and SrMnBi 2 suggests ionic substitution as a potential tuning mechanism. If so, novel behavior near the compensation point could be expected. Here, to explore this, we study a mixed-cation analog, Ca 0.6⁢ Sr 0.4⁢ MnSb 2 , where Bi is replaced by Sb. Conventional inelastic neutron scattering is impractical due to the small crystal size (m ≈ 0.28 g) available for compositional studies; however, we find that using a neutron diffractometer with a wide-angle area detector we can obtain a good quality spin-wave signal, which is shaped by energy-momentum conservation and retains spectroscopic information even without direct energy analysis. Spin-wave modeling reveals an interlayer coupling quantitatively similar to SrMnSb 2 , indicating it is not directly tuned by ionic size and that the sign change in Bi-based compounds likely arises from the observed change in lattice symmetry. Beyond this key insight, our results present an efficient method for parametric and compositional studies of spin dynamics in small crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The corrosion behavior of borosilicate glass in the presence of cementitious waste forms

Borosilicate glasses are widely used for radioactive waste disposal due to their ability to incorporate a variety of contaminants and radionuclides while exhibiting high durability in various disposal scenarios. This research evaluated the dissolution of borosilicate glass using both single-pass-flow-through (ASTM C1662-18) and product consistency test (ASTM C1285-21) methods with different solutions, including a cementitious-contacted water (called grout-contacted, GC, from this point) and solutions with varying levels of dissolved cementitious species such as Si, Ca, Al. Here, the results indicated that the presence of Ca plays a crucial role in suppressing glass corrosion, as evidenced by the slower normalized dissolution rates, which were one order of magnitude lower for boron and two orders of magnitude lower for rhenium, observed in both Ca-amended and GC solutions compared to the pH 12 buffer solution. This effect is attributed to the formation of a dense, low-porosity, and strongly bonded calcium silicate hydrate (CSH) layer on the glass surface, which implies that a glass corrosion process is influenced by ion exchange involving alkali ions Na + , K + , Ca 2+ , and hydrogen-containing species. A small number of glass particles treated in the GC solution showed minor corrosion pits in the form of shallow craters with an average diameter of approximately 500 μm. This observation is correlated with a significant reduction, 2000 to 3000 times lower, in the cumulative volume of glass pores, indicating that smaller pore voids were “sealed” in the presence of Ca 2+ ions, likely attributed to the formation of CSH precipitation or other corrosion products such as calcium carbonate saturated from the grout solution. These findings suggest that the presence of dissolved Ca in the GC solution can slow down the dissolution of borosilicate glass, contrary to the expected trend of higher dissolution rates resulting from exposure to high alkaline and thus higher pH solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermal conductivity suppression in uranium-doped thorium dioxide due to phonon-spin interactions

In this work, impact of low level of uranium (U) atom substitution on thermal conductivity of thorium dioxide (ThO 2 ) is investigated. ThO 2 is an electronic insulator with a wide optical band-gap and no unpaired electrons whose thermal transport is governed by phonons. U-substitution introduces unpaired f-electrons resulting in paramagnetic behavior of U-ThO 2 at room temperature, which significantly suppresses its thermal conductivity. A single crystal of U-ThO 2 with graded composition of U is grown using a hydrothermal synthesis method, and thermal conductivity measurements are performed in regions with uniform composition of U at levels of 0%, 6%, 9% and 16%. Measured thermal conductivity profiles over 77–300 K temperature range are analyzed using an analytical expression for phonon-mediated thermal transport based on Klemens-Callaway model. Temperature dependent thermal conductivity is found to deviate significantly from the Rayleigh scattering trend expected for a simple substitutional point defect with a small perturbation to mass and interatomic forces. With the resonant scattering term, observed large suppression of thermal conductivity at low temperatures can be closely reproduced. Additionally, the extracted phonon-spin coupling constants imply a nonlinear relation of phonon-spin interaction intensity with respect to U doping percentage. Our study reveals how phonon-spin scattering contributed by unpaired f-electrons in U atoms influences thermal transport in the U-ThO 2 system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2025

Innovative use of CaO in combination with amino acid salt to convert CO 2 as CaCO 3 nanoparticles under mild pH and low temperature

The increasing demand for sustainable CO 2 management has driven the development of innovative methods that can convert point source CO 2 into value-added products. In this study, CaO in combination with amino acid salt was used to convert CO 2 into CaCO 3 nanoparticles. Different from the conventional method where CO 2 diffusion presents a major challenge in reacting with CaO to form CaCO 3 , amino acid salt solvent was applied to absorb CO 2 first and then rapidly reacted with CaO to form CaCO 3 nanoparticles (∼50 nm) at a low temperature (e.g., 60 °C). Our experiments showed that at a glycine (Gly)/NaOH ratio of 2:1 or 3:1, the solution pH values during the CO 2 absorption and conversion were about 8–9 at 60 °C, while at a ratio of 1:1, the solution pH values were about 9–11; without Gly, the solution pH values were about 12. Gly-optimized solvent substantially reduced corrosion risk to reactors. In addition, the use of amino acid (i.e., Gly) led to much smaller CaCO 3 particles, distinctly different chemical phases, and fundamentally different chemical reactions. Moreover, in the presence of Gly, the solution pH was completely reversed and the solution was regenerated for cyclic use when CaO was added. The solvent was recyclable and reusable, highlighting the cost-effectiveness and sustainability of this approach. The Gly-modulated CaCO 3 nanoparticles may have significant potential for industrial applications in the biomedicine, construction, plastics, and rubber industries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Aboveground Rather Than Belowground Productivity Drives Variability in Miscanthus × giganteus Net Primary Productivity

Quantifying the carbon (C) uptake of Miscanthus × giganteus ( M × g ) in both aboveground and belowground structures (e.g., net primary productivity (NPP)) and differences among methodological approaches is crucial. Our objectives were to directly measure Mxg NPP and evaluate the effects of nitrogen application, location, and belowground biomass sampling methods. We hypothesize that increased nitrogen application increases the overall NPP of M × g and that quantifying rhizome biomass using excavations will produce the lowest variability between replicates. We collected biomass from mature M × g stands from three locations in Iowa with three nitrogen application rates and one site in Illinois. We destructively sampled at two time points, when rhizome mass is anticipated to be at a minimum (initial) and anticipated to be at its maximum (peak). Biomass was collected from 1 × 1 m quadrats in which one in-clump and one beside-clump cores were collected and then excavated to 30 cm depth to extract all rhizomes. We found that aboveground M × g NPP ranged from 15.4 Mg DM ha –1 year –1 to 36.4 Mg DM ha –1 year–1 and belowground M × g NPP ranged from 4.4 Mg DM ha –1 year –1 to 19.6 Mg DM ha –1 year –1 . M × g NPP varied across sites, fertilization, and calculation assumptions. Aboveground NPP (yield) was on average 68.7% of the total NPP. Root-to-shoot ratios at peak biomass decreased with nitrogen application rate, from an average of 1.9 for 0 N plots to 0.89 for 224 N fertilized plots. There was more variation in core data than from excavations; however, when in-clump and beside-clump cores were averaged together, core and excavation averages were not different. Overall, these results show that the range of mature M × g NPP is driven by aboveground productivity, influenced by nitrogen application and site. Our results provide useful data to constrain agro-ecosystem models and provide crucial insights for future perennial belowground sampling.

Hartman, Theodore [Univ. of Illinois at Urbana-Cha↗
DOE OSTI2026

A systematic analytical framework for multi-source municipal solid waste characterization for energy recovery

Advancing municipal solid waste (MSW) management from disposal-oriented practices toward circular, value-driven systems requires standardized methodologies capable of identifying material composition and resource recoverable potential at the point of generation. Despite extensive research, MSW characterization remains fragmented due to inconsistences in sampling methodologies, waste sorting categories, and temporal coverage across previous studies which limit cross-site comparability, reproducibility, and constrain the reliable evaluation of potential resource recovery pathways. This lack of consistency has hindered the development of a unified framework for MSW characterization and resource assessment. This study introduces a standardized, field-validated protocol for MSW sampling and composition analysis that ensures consistent, traceable data across diverse waste sources. The protocol integrates randomized spatial sampling, systematic material sorting, and controlled subsampling for multi-site and multi-season field campaigns. Validation included MSW collection from residential, grocery, restaurant, and school MSW streams across five U.S. states, including Maryland, Idaho, Virginia, Ohio, and Mississippi, to demonstrate the protocol’s ability to identify source-based composition patterns relevant to resource recovery applications. Grocery and restaurant streams were dominated by food waste and high-moisture organics, while school waste contained higher paper content and residential waste showed greater heterogeneity. Aggregation into energy-relevant fractions highlighted practical recovery pathways via anaerobic digestion or gasification, supporting data-driven planning, policy, and circular economy strategies for sustainable waste management across waste sources.

09 BIOMASS FUELS↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Woody Feedstock 2022 State of Technology Report

The U.S. Department of Energy promotes production of advanced liquid transportation fuels from lignocellulosic biomass by funding fundamental and applied research that advances the state of technology (SOT). As part of its involvement in this mission, Idaho National Laboratory completes an annual SOT report for n th -plant and 1 st -plant woody biomass feedstock logistics. The purpose of the SOT is to provide the status of feedstock supply system technology development for woody biomass to biofuels relative to technical targets and cost goals from specific design cases, based on data and experimental results. Conventional feedstock supply systems need to be modified to meet the demands of conversion pathways, specifically to have the ability to adjust the quality of the raw biomass materials. Advanced systems incorporate innovative methods of material handling, preprocessing and supply chain configuration. In advanced designs, variability of the raw biomass can be reduced to produce feedstocks of a uniform format, moving toward biomass commoditization. Against this backdrop, the 2022 Woody SOT for low-ash woody feedstocks utilizes feedstock fractionation by incorporating technologies that can separate the biomass into its anatomical fractions (wood, bark, needle, and extrinsic ash) to reduce impurities and attempt to maximize the retention of usable fractions that satisfy downstream quality considerations. By using a series of air classification steps, this strategy can reduce the extrinsic ash in forest residues, separate out a majority of the incoming needles (which can be supplied to alternate markets), and maximize the retention of whitewood in the usable fraction. The fractionated forest residues are then mixed with clean-pine chips in a 50-50 blend to prepare the feedstock for the desired conversion pathway. The n th -plant analysis estimated the delivered cost for the feedstock at $\$$69.23/dry ton (2016$\$$) which represents a $\$$6.64/dry ton decrease compared to the cost estimate of the 2021 Woody SOT supply system for low-ash woody feedstocks. The quality requirements in the 2022 Woody SOT were identical to those of the 2021 Woody SOT at = 1.00 wt % ash and = 50.51 wt% carbon. The cost savings derive primarily from reductions in dry matter losses during air classification. The GHG emissions for the n th -plant analysis were estimated at 178.39 kg CO2e/dry ton compared to 178.71 kg CO 2 e/dry ton in the 2021 Woody SOT, a decrease of 0.32 kg CO2e/dry ton. The small change stems from an increase in emissions attributed to preprocessing and slightly larger savings in emissions from transportation. In the 1 st -plant analysis of the 2022 Woody SOT system, the average throughput was estimated to be approximately 2,128 dry tons/day or 96.51% of the name plate capacity. During the simulation the daily throughput ranged from 1,090 dry tons/day to 2,200 dry tons/day, or 49.43% to 99.75% of the daily nameplate capacity. After the year of operation 722,403 tons of processed feedstock were produced in total without regard to quality considerations (99.64% of the annual nameplate capacity). The variability in throughput was primarily caused by equipment failures in the system. Regular failures, downtime caused by routine maintenance per manufacturer guidelines, contributed to a majority 62.50% of failures and 62.60% of downtime. Failures due to wear were the other cause of disruption within the system, impacting the rotary shear and orbital screen and accounting for 37.50% of the failures and 37.40% of the total downtime. Ultimately the system was on stream for 87.84% during the simulation period, which is only 2.16 percentage points below the nth-plant assumption for on-stream time. The production cost of the system averaged $\$$71.66/dry ton. The costs ranged from a minimum of $\$$71.23/dry ton to a maximum of $\$$2,115.30/dry ton. When dry matter losses (disposed low-quality fractions as well as other losses such as in grinders) were considered the costs increased to an average of $\$$75.11/dry ton with a minimum of $\$$74.69/dry ton and a maximum of $\$$2,136.86/dry ton...

09 BIOMASS FUELS↗
DOE OSTI2026

Pixelated plastic scintillator array manufacturing using fast-, photo-curable resin

Pixelated plastic scintillator arrays can serve as high efficiency and high resolution neutron imaging detectors. Manufacturing these arrays is intensive in both time and labor. This article presents a fabrication method based on additive manufacturing for two-dimensional plastic organic scintillator arrays using a custom-built automated assembly machine and a custom photocurable resin that has significant non-aromatic acrylate oligomer content. The process involves two main stages: fully autonomous production of one-dimensional layered arrays, followed by semi-autonomous cutting and stacking to form two-dimensional pixel arrays. One-dimensional arrays were manufactured at a rate of around 4 layers per hour with minimal defects and tight dimensional tolerances, while two-dimensional arrays up to 7 x 7 pixels and 70 mm in length were completed in approximately 3.5 hours. Final arrays exhibited dimensional deviations of less than 0.5 mm. Two-dimensional arrays read out by a multi-anode photomultiplier tube demonstrated per-pixel position resolution and pulse-shape discrimination, enabling gamma–neutron interaction separation in mixed radiation environments.

36 MATERIALS SCIENCE↗