Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
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 OSTI2021

Microstructure evolution of gadolinium doped cerium oxide under large thermal gradients

In this report the effects of large thermal gradient annealing on the microstructure of 10 mol% gadolinium doped ceria (GDC) were investigated. GDC powder was prepared by solvent deficient method and sintered at 1650 °C for 10 h to achieve dense ceramics with ~8 μm grain size. The densified GDC samples were subsequently annealed using a 60 W infrared laser at over 2100 °C for 1 h under a thermal gradient equivalent to ~0.3–0.5 °C/μm. The post-annealed samples at 2150 °C for 1 h exhibit grains with average length and width of 37 and 28 μm, respectively. Electron backscattered diffraction (EBSD) analysis revealed that the post-annealed sample at 2150 °C consists of grains oriented close to five principal directions (<4 3 10>, <0 0 1> and <13 1 14> on [0 0 1], and <7 6 20> and <7 2 7> on [0 1 0]) within a tolerance angle of ±10°, whereas the grains of the pre-annealed sample are randomly oriented. Gadolinium diffuses 20–30 μm away from the irradiated surface, with the measured composition of regions deeper than 30 μm, Ce 0.86 Gd 0.14 O 1.93 , is close to that of the pre-annealed sample, Ce 0.87 Gd 0.13 O 1.94 . Enhancement of total conductivity of the post-annealed GDC (1.1 × 10 -3 S cm -1 at 500 °C, and 2.1 × 10 -2 S cm -1 at 700 °C) is observed when compared to the pre-annealed GDC (3.1 × 10 -5 S cm -1 at 500 °C, and 1.7 × 10 -3 S cm -1 at 700 °C), and points to the decrease in the grain boundary (GB) resistivity. This could be attributed to both the change in GB area and grain alignment.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impact of high-dose gamma-ray irradiation on electrical characteristics of N-polar and Ga-polar GaN p–n diodes

We investigate the impact of high-dose gamma-ray irradiation on the electrical performance of Ga-polar and N-polar GaN-based p-n diodes grown by metalorganic chemical vapor deposition. We compare the current density-voltage (J-V), capacitance-voltage (C-V), and circular transfer length method (CTLM) characteristics of the p-n diodes fabricated on Ga-polar and N-polar orientations before and after irradiation. The relative turn-on voltage increases for the Ga-polar diodes with increasing irradiation dose, while it increases initially and then starts to decrease for the N-polar diodes. The p-contact total resistance increases for Ga-polar and decreases for N-polar samples, which we attribute to the formation of point defects and additional Mg activation after irradiation. The J-V characteristics of most of the tested diodes recovered over time, suggesting the changes in the J-V characteristics are temporary and potentially due to metastable occupancy of traps after irradiation. X-ray photoelectron spectroscopy (XPS) and photoluminescence (PL) measurements reveal the existence of different types of initial defects and surface electronic states on Ga-polar and N-polar samples. Gallium vacancies (V Ga ) are dominant defects in Ga-polar samples, while nitrogen vacancies (V N ) are dominant in N-polar samples. The presence of a higher concentration of surface states on Ga-polar surfaces compared to N-polar was confirmed by calculating the band bending and the corresponding screening effect due to opposite polarization bound charge and ionized acceptors at the surface. The difference in surface stoichiometry in these two orientations is responsible for the different behavior in electrical characteristics after gamma-ray interactions.

36 MATERIALS SCIENCE↗
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 OSTI2022

Control over emergence and alignment of transient blisters in thermo-responsive gels using hierarchically patterned substrates

Surfaces with tunable microscale textures are vital in a large variety of technological applications, including heat transfer, antifouling and adhesion. To facilitate such broad-scale use, there is a need to create surfaces that undergo reconfigurable changes in topology and thus, enable switchable functionality. To date, there is a relative dearth of methods for engineering surfaces that can be actuated to change topography over a range of length scales, and hence, form tunable hierarchically structured layers. Combining modeling and experiments, we design a geometrically patterned, thermo-responsive poly (N-isopropylacrylamide) gel film that undergoes controllable hierarchical changes in topology with changes in temperature. At the bottom, the film is covalently bound to a solid, curved substrate; at the top, the film encompasses longitudinal rectangular ridges that are oriented perpendicular to the underlying cylindrical curves. At temperatures below lower critical solution temperature (LCST), the swollen gel exhibits 3D variations in polymer density and thickness defined by the gel's top and bottom topography. As the temperature rises above LCST, the interplay between the upper ridges and lower curves in the gel drives non-uniform, directional solvent transport, the nucleation and propagation of a phase-separated higher-density skin layer, and the resulting pressure buildup within the film. These different, interacting kinetic processes lead to an instability, which produces transient microscopic blisters in the film. Through simulations, we show how tuning the width of the ridges modifies the propagation of a skin layer and creates localized pressure build-up points, which enables control over the emergence, distribution, and alignment of the microscopic blisters. Futhermore, we provide a simple argument to predict the size of such microscopic features. Experiments confirm our predictions and further highlight how our computational model enables the rational design of topographical transitions in these tunable films. The development of actuatable, hierarchically structured films provides new routes for achieving switchable functionality in actuators, drug release systems and adhesives.

42 ENGINEERING↗
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 OSTI2023

Migration of Ga vacancies and interstitials in β – Ga 2 O 3

Pathways and energy barriers for the migration of Ga vacancies (V Ga ) and Ga interstitials (Ga i ) in β–Ga 2 O 3 are explored using hybrid functional calculations and the nudged elastic band method. Considering β–Ga 2 O 3 as primarily being an n-type semiconductor, we focus on defect charge states relevant under such conditions: $V^{3 –}_{Ga}$, Ga$^{3+}_{i}$, and Ga$^{+}_{i}$. Notably, we describe a mechanism by which V Ga can transform between its different split configurations. In all cases, the intermediate state consists of a vacancy split between three Ga sites—a three-split vacancy—which enables passage over a significantly lower energy barrier. This is because it avoids the unfavorable simple vacancy at the tetrahedral Ga site. The proposed mechanism lowers the overall barrier for $V^{3–}_{Ga}$ diffusion along the [001] crystal direction from 1.73 to 0.97 eV, whereas the 2.08 eV barrier for the [100] and [010] directions is unaffected. For Ga$^{3+}_{i}$, we obtain similar overall migration barriers of 0.72, 0.80, and 1.02 eV for the [010], [001], and [100] directions, respectively. Ga$^{+}_{i}$ exhibits a strong preference for diffusion within the large eight-sided channel; the overall migration barrier is 0.92 eV for the [010] direction, and 2.16 eV for the [001] and [100] directions. The limiting step for the two latter directions involves ionization of Ga$^{+}_{i}$ followed by a jump to an adjacent large eight-sided channel as Ga$^{3+}_{i}$. Finally, our results are discussed in light of experimental observations of thermally activated recovery processes in irradiated material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chiral Topological Surface States on a Finite Square Photonic Crystal Bounded by Air

Chiral, topologically protected, photonic surface states can be found at the boundary between gyrotropic photonic crystals where a changing magnetic field induces different topology across the interface. Typically, photonic crystals with either a suitable band structure on both sides of the interface to provide a band gap and evanescent decay of the surface states away from the interface, or an outer layer with engineered material properties is required. In this paper, we show the emergence of topological, unidirectional surface states at the termination of finite gyrotropic photonic crystals with a simple square lattice and C 4 rotational symmetry bounded by a vacuum, eliminating the need for an outside layer to enable chiral surface modes. Here, we start from an infinite, time-reversal-symmetry-breaking photonic crystal with a band gap associated with bands with nonzero Chern numbers, different from all-zero Chern numbers in air. We then modify the photonic crystal to move this band gap below the light line, while maintaining the Chern-number discontinuities. Band-structure calculations for a supercell approximating a photonic crystal finite in the direction normal to the surface demonstrate the existence, dispersion, and chirality of the surface mode. Extensive direct scattering calculations for a point source and spatial Fourier analysis further reveal a unidirectional free-space topological surface state, which propagates counterclockwise around the surface of a finite photonic crystal, providing a nearly foolproof way to cross-check the surface-mode band structure unaffected by backscattering from local defects. Additionally, scattering simulations allow an independent characterization of the state dispersion and unveil the robustness of the topological plasmonic mode propagation around the 90° bends of the structure, being due to only radiation leakage. In contrast to buried topological surface states, the observed surface modes at the photonic crystal–air interface have the advantage of being accessible to the outside world, allowing one to take advantage of the defect-tolerant backscattering-free surface modes to engineer emission from photonic crystal surfaces into arbitrary free-space beam shapes and directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2020

Accelerated Aging Humidity Chamber for Nuclear Grade HEPA Filter Media

A unique accelerated aging humidity chamber for simulating the natural aging process of nuclear grade high efficiency particulate air (HEPA) filter media is being designed, implemented, and tested. The nuclear industry currently lacks information regarding the shelf life and service life of HEPA filter media. The Department of Energy recommends disposing of HEPA filters 10 years after the manufacture date for filters operated under dry conditions. Studies have shown the tensile strength and water repellency of HEPA filter media diminish with age. The chamber's principal function is to sustain an elevated relative humidity condition for an extended duration specified by the Arrhenius aging model. To accomplish this, various design parameters were established: (1) consistently achieve uniform humidity conditions within the chamber to equally expose media sheets, (2) ensure the chamber is sufficiently sealed to mitigate the uncontrolled infiltration of air, (3) confidently secure the media sheets to avert the introduction of unnecessary stressing/creasing, (4) fabricate the chamber from clear polycarbonate material to allow inspection of the media sheets, (5) provide a safe and ergonomic design for personnel, (6) ensure the chamber design is replicable and relatively easy to fabricate/assemble. A humidity source, fan, PID controller, humidity and temperature transmitters, and data loggers are required for proper operation and control of the chamber. The humidity source must accommodate target values up to 95% relative humidity throughout the chamber. The three humidity and temperature transmitters are strategically spaced throughout the chamber to ensure uniform conditions. The 10 cubic feet per minute (CFM) fan selected intends to accomplish the desired air change rate inside the chamber of approximately 20 air changes per hour (ACH). The 2 plenums on the upper and lower portions of the chamber utilize a perforated design to uniformly distribute air, promote air mixing, and control the air velocity entering the chamber. Industrial blueprint hanging clamps are being retrofitted to be seated within the chamber drawers and firmly secure the nuclear grade filter media. The aged media will be evaluated using autopsy methodology, and the results obtained intend to help clarify the useful life of nuclear grade HEPA filters. Nuclear grade high efficiency particulate air (HEPA) filters are defined as disposable, extended-media, dry-type filters with a rigid casing enclosing the full depth of the pleats and have a minimum particle removal efficiency of 99.97%. [1] - HEPA filters are constantly exposed to stressors and subsequently have been shown to degrade as the filter ages. Properties such as tensile strength and water repellency, along with others, are negatively impacted. [2] - Due to the scarcity of naturally aged HEPA filters available, it is crucial to develop accelerated aging methods which effectively mimic the natural degradation effects. - After accelerated aging, the filters are subject to testing in the Axial Flow Large Scale Test Stand (ALSTS) and other analysis techniques such as the same qualification tests performed directly after being manufactured. - Quality data is necessary for future decisions regarding the lifetime of nuclear grade HEPA filters. - The objective is to develop a prototype accelerated aging humidity chamber for nuclear grade HEPA filter media. The chamber should expose the media to commonly occurring stressors to expedite the degradation process caused by aging. PID yields reliable and reproducible results: - CFD model correlates well with actual performance; - Further analysis is needed to complete characterization; - Future work for the chamber is still to be completed. The data shown in Figure 5 and Table 1 indicate the system is functioning as designed. As can be determined from the system response evaluation in Table 1, the PID controller is yielding reliable and reproducible results at each set point. As the set point increases, the time constant and 20% settling time increase accordingly. From the results in Table 1, it stands to reason that the CFD model correlates well with the actual performance of the system. Each system response parameter for the 90%RH* CFD model is slightly lower than the parameters collected from the physical model tested at an equivalent set point. This is to be expected considering the inlet in the CFD model is a constant supply of water vapor, whereas the physical model uses a PID controller to monitor the amount of water vapor needed and reacts accordingly. More analysis is still to be done on the prototype chamber. A pressure test will be conducted in order to determine the hourly leak rate of the chamber. This test will involve constant air flow into the chamber until the desired pressure is reached. A uniformity test will also be conducted in order to ensure each filter medium is being exposed to equivalent relative humidity levels. This test will involve an array of humidity meters strategically oriented inside the chamber to record any possible gradients. Future work for the chamber includes dehumidification abilities, tests while fully loaded with filters, and temperature control. Once completed, the autopsy team at Institute of Clean Energy Technology will conduct necessary accelerated aging studies as needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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↗