Search NASA⌕ Search

SEARCH · Search NASA

Results for “granular systems”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

When and Where Lithium Plating Occurs, Its Correlation with Microstructure Heterogeneity, and the Mechanisms That Initiate and Self-Regulate Electrochemical Heterogeneity (A02-0444)

A microstructure scale electrochemical LIB model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. Model predicts active material particle surface roughness and size distribution (respectively, non-uniform curvature within and between particles) initiate in-plane heterogeneity, and that particle size heterogeneity at the separator interface controls the lithium plating preferential deposition ("Where"). These in-plane heterogeneities are then exacerbated by through-plane heterogeneities induced at fast charge as electrolyte depletion occurs and concentrates intercalation reaction near the anode-separator interface. Also, magnitude and occurrence of lithium plating is controlled by effective, or macroscale, microstructure parameters ("When"). As local states of charge start to diverge between nearby active material regions, overpotential differences induced by OCP difference kick in and contribute to reduce these SOC local heterogeneities. However, for staged materials such as graphite, with OCP profile alternating between plateaus and varying regions, this balancing mechanism is, respectively, inactive and active. This leads to a dynamic, non-monotonic, in-plane heterogeneity time evolution for state of charge and Faraday current density, for which their respective in-plane heterogeneity magnitude alternates. Such behavior has been modeled both for the whole electrode at the microstructure scale and at the particle scale. In-plane heterogeneities are usually considered to be detrimental, as they result in material non-uniform utilization (i.e., under and over stressed regions) and earlier degradations. However, this work provides a more granular approach as it discriminates between a harmful in-plane heterogeneity (non-uniform curvature) that triggers SOC in-plane heterogeneity, and a beneficial in-plane heterogeneity (Faraday current density) that contributes to reduce SOC in-plane heterogeneity. This work comprehensively explains the mechanisms that initiate, exacerbate, and regulate heterogeneity at the microstructure scale, while providing some design suggestions to reduce both in-plane and through-plane heterogeneities, as summarized in the graphical abstract.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Nanoparticle Tracking to Probe Transport in Porous Media

From the granular and fractured subsurface environment to highly engineered polymer membranes used in pharmaceutical purification, porous materials are ubiquitous in nature and industrial applications. In particular, porous media are used extensively in processes including water treatment, pharmaceutical sterilization, food/beverage processing and heterogeneous catalysis, where hindered mass transport is either essential to the process or a necessary, but undesirable limitation. Unfortunately, there are currently no universal models capable of predicting mass transport based on a description of the porous material because real porous materials are complex, and because many coupled dynamic mechanisms (e.g. adsorption, steric effects, hydrodynamic effects, electrostatic interactions, etc.) give rise to the observed macroscopic transport phenomena. While classical techniques, like nuclear magnetic resonance and dynamic light scattering, provide useful information about mass transport in porous media at the ensemble level, they provide limited insight into the microscopic mechanisms that give rise to complex phenomena such as anomalous diffusion, hindered pore-space accessibility, and unexpected retention under flow, among many others. To address this issue, we have developed refractive index matching imaging systems, combined with single-particle tracking methods, allowing the direct visualization of single particle motion within a variety of porous materials. Here, we summarize our recent efforts to advance the understanding of nanoparticle transport in porous media using single-particle tracking methods in both fundamental and applied scenarios. First, we describe the basic principles for two-dimensional and three-dimensional single-particle tracking in porous materials. Then, we describe concrete examples of nanoparticle transport in porous materials from two perspectives: (1) understanding fundamental elementary particle transport processes in porous media, including pore accessibility and cavity escape, which limit transport in porous media; (2) facilitating applications in industrial processes, e.g. by understanding the mechanisms of particle fouling and remobilization in filtration membranes. Finally, we provide an outlook of opportunities associated with investigating other types of mass transport in confined environments using single-particle tracking methods, including electrophoretic and self-propelled motion.

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

Characterizing the Influence of Organic Polymers on the Specific Reactivity of Particulate Remedial Amendments

Commercially available particulate amendments demonstrate high reactivity for effective treatment of water soluble organic and inorganic contaminants in laboratory studies; however, transport of these particles is constrained in the subsurface. In many field applications, particulate amendments are mixed with organic polymers to enhance mobility for direct push applications or stabilize suspensions for high mass loadings. As such, the interactions between particulate amendments, organic polymers and contaminant species need to be systematically investigated to properly understand mechanistic processes that facilitate predictive performance metrics for specific applications in situ . In this study, batch experiments were conducted to quantify the effects of organic polymers (xanthan gum, guar gum, and sodium alginate), polymer concentration (800 and 4,000 mg/L), and aging (up to 28 days) on chromate treatment rate and capacity by two classes of amendments: reductants [granular zero-valent iron (gZVI), micron-ZVI (mZVI), sulfur modified iron (SMI)], and an adsorbent (bismuth sub-nitrate). When particulate amendments were suspended in polymer solutions, reductants retained between 84–100% of the amendment treatment capacity. Conversely, the adsorbent maintained 63–97% relative treatment capacity of the no-polymer control. Polymer solutions had a more pronounced impact on the rate of chromate removal; first order rates of chemical reduction decreased by as much as 70% and adsorption by up to 81% relative to the no-polymer controls. Polymer–amendment aging experiments also showed decreased Cr(VI) treatment capacity; reductants decreased by as much as 24% and adsorption decreased by as much as 44% after 28 days of incubation. While polymer suspensions are needed to aid the injection of particulate amendments into the subsurface, the results from this study indicate potential losses of treatment capacity and a decrease in the rate of remedial performance due to the physical and chemical interactions between polymer suspensions and reactive particulate amendments. Simple batch systems provide baseline characterization of tripartite interactions for the removal of Cr(VI). Additional work is needed to quantify the full impact of polymers on remedial outcomes under site relevant conditions at field scale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Experimental demonstration of optical stochastic cooling

Particle accelerators and storage rings have been transformative instruments of discovery, and, for many applications, innovations in particle-beam cooling have been a principal driver of that success. Stochastic cooling (SC), one of the most important conceptual and technological advances in this area cools a beam through granular sampling and correction of its phase-space structure, thus bearing resemblance to a ‘Maxwell’s demon’. The extension of SC from the microwave regime up to optical frequencies and bandwidths has long been pursued, as it could increase the achievable cooling rates by three to four orders of magnitude and provide a powerful tool for future accelerators. First proposed nearly 30 years ago, optical stochastic cooling (OSC) replaces the conventional microwave elements of SC with optical-frequency analogues and is, in principle, compatible with any species of charged-particle beam. Here we describe a demonstration of OSC in a proof-of-principle experiment at the Fermi National Accelerator Laboratory’s Integrable Optics Test Accelerator. The experiment used 100-MeV electrons and a non-amplified configuration of OSC with a radiation wavelength of 950 nm, and achieved strong, simultaneous cooling of the beam in all degrees of freedom. This realization of SC at optical frequencies serves as a foundation for more advanced experiments with high-gain optical amplification, and advances opportunities for future operational OSC systems with potential benefit to a broad user community in the accelerator-based sciences.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Closed-form existence conditions for bandgap resonances in a finite periodic chain under general boundary conditions

Bragg scattering in periodic media generates bandgaps, frequency bands where waves attenuate rather than propagate. Yet, a finite periodic structure may exhibit resonance frequencies within these bandgaps. This is caused by boundary effects introduced by the truncation of the nominal infinite medium. Previous studies of discrete systems determined existence conditions for bandgap resonances, although the focus has been limited to mainly periodic chains with free–free boundaries. In this paper, we present closed-form existence conditions for bandgap resonances in discrete diatomic chains with general boundary conditions (free–free, free–fixed, fixed–free, or fixed–fixed), odd or even chain parity (contrasting or identical masses at the ends), and the possibility of attaching a unique component (mass and/or spring) at one or both ends. The derived conditions are consistent with those theoretically presented or experimentally observed in prior studies of structures that can be modeled as linear discrete diatomic chains with free–free boundary conditions. An intriguing case is a free–free chain with even parity and an arbitrary additional mass at one end of the chain. Introducing such an arbitrary mass underscores a transition among a set of distinct existence conditions, depending on the type of chain boundaries and parity. The proposed analysis is applicable to linear periodic chains in the form of lumped-parameter models, examined across the frequency spectrum, as well as continuous granular media models, or similar configurations, examined in the low-frequency regime.

Bastawrous, Mary V.↗
DOE OSTI2020

Structural Tolerance Factor Approach to Defect-Resistant I 2 -II-IV-X 4 Semiconductor Design

Recent work on quaternary semiconductors Cu 2 BaSn(S,Se) 4 and Ag 2 BaSnSe 4 for photovoltaic and thermoelectric applications, respectively, has shown the promise of exploring the broader family of defect-resistant I 2 -II-IV-X 4 materials (where I, II, and IV refer to the formal oxidation state of the metal cations and X is a chalcogen anion) with tetrahedrally coordinated I/IV cations and larger II cations (i.e., Sr, Ba, Pb, and Eu) for optoelectronic and energy-related applications. Chemical dissimilarity among the II and I/IV atoms represents an important design motivation because it presents a barrier to antisite formation, which otherwise may act as electronically harmful defects. We herein show how all 31 experimentally reported I 2 -II-IV-X 4 examples (with large II cations and tetrahedrally coordinated smaller I/IV cations), which form within five crystal structure types, are structurally linked. Based on these structural similarities, we derive a set of tolerance factors that serve as descriptors for phase stability within this family. Despite common usage in the well-studied perovskite system, Shannon ionic radii are found to be insufficient for predicting metal–chalcogen bond lengths, pointing to the need for experimentally derived correction factors as part of an empirically driven learning approach to structure prediction. Here, we use the tolerance factors as a predictive tool and demonstrate that four new I 2 -II-IV-X 4 compounds, Ag 2 BaSiS 4 , Ag 2 PbSiS 4 , Cu 2 PbGeS 4 , and Cu 2 SrSiS 4 , can be synthesized in correctly predicted phases. One of these compounds, Ag 2 PbSiS 4 , shows potentially promising optoelectronic properties for photovoltaic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The mechanical response of micron-sized molecular crystals

Abstract Microstructures and corresponding properties of compacted powders ultimately depend on the mechanical response of individual particles. In principle, computational simulations can predict the results of powder compaction processes, but the selection of appropriate models for both particle–particle interactions and particle deformations across all relevant length scales remain nontrivial tasks, especially in material systems lacking detailed mechanical property information. The work presented here addresses these issues by conducting uniaxial compressions in situ inside of a scanning electron microscope to characterize the mechanical response of individual micron-sized particles of a molecular crystal, hexanitrohexaazaisowurtzitane (CL-20). This experimental approach enabled the collection of quantitative force and displacement data alongside simultaneous imaging to capture morphology changes. The results reveal information about elastic deformation, yield, plastic deformation, creep, and fracture phenomena. Accordingly, this work demonstrates a generalizable approach for assessing the mechanical response of individual micron-sized molecular crystal particles and utilizing those responses in particle-level models. Graphic abstract

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermal Evolution of Internal Strain in Doped PbTe

Recent improvements in the efficiency of heat-to-electricity energy conversion in lead chalcogenide thermoelectrics involve reducing the thermal conductivity by incorporating large amounts of internal strain. The extent to which typical lead chalcogenide processing techniques (such as doping, ball milling, and densification) increase internal strain and dislocation density must be quantified to improve materials design. In this study, neutron powder diffraction is leveraged to evaluate the internal strain introduced by ball milling in doped and undoped powders. Doping with Na and/or Eu increases internal strain beyond ball milling alone, with the greatest increase from combining the two dopants. Strain recovery occurs in each powder above 400 K but can be suppressed by co-doping, indicating a strong dopant–dislocation interaction in this system. Therefore, high-temperature processing of PbTe powders should be avoided if high internal strain is desired. Low-temperature densification and/or rapid pressing techniques may be key to maintaining internal strain in the final pressed pellet. The diffraction peak asymmetry and correlated elastic softening measured in pressed PbTe pellets in past studies were not observed in the precursor powders measured here, suggesting that measurements of the Debye temperature on final pressed pellets are required to examine the influence of defect-induced elastic softening on thermal conductivity. This work provides key guidance for defect engineering to maximize internal strain and thermoelectric performance in PbTe thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Optimal binning of correlated measurements

Experimental measurements are commonly represented on a discrete grid, requiring a balance between granularity and statistical noise. Two strategies have traditionally been used to improve such representations: selecting an appropriate bin width to control discretization error and applying kernel-based smoothing to suppress fluctuations. Despite their shared goal, these approaches have largely developed independently, without a unified statistical description of how discretization and correlation jointly determine measurement precision. Here, we extend the discussion of optimal interval averaging to a correlation-aware setting by Gaussian process regression, which explicitly accounts for correlations among neighboring bins. Starting from first principles, we derive the mean-squared error of discretized measurements and obtain closed-form asymptotic expressions for the optimal bin width and correlation length. When recast in reduced variables, the theory reveals distinct universal scaling laws governing the error in the correlation-free and correlation-controlled regimes. Characterized by intrinsically smooth intensity profiles and counting-based statistics, neutron scattering measurements are well suited for demonstrating the enhanced error contraction enabled by inter-bin correlations. We show that such improvement is achievable over the experimentally accessible Q-range and across multiple instruments and material systems. These results show that explicitly accounting for correlations systematically reshapes the limits of precision in discretized, noise-limited measurements. More broadly, the framework provides a transferable statistical foundation for optimizing data representation, inference, and experimental design across the physical and data sciences.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗
DOE OSTI2025

A Typology of Decision-Making Tasks for Visualization

Despite decision-making being a vital goal of data visualization, little work has been done to differentiate decision-making tasks within the field. While visualization task taxonomies and typologies exist, they often focus on more granular analytical tasks that are too low-level to describe large complex decisions, which can make it difficult to reason about and design decision-support tools. In this paper, we contribute a typology of decision-making tasks that were iteratively refined from a list of design goals distilled from a literature review. Our typology is concise and consists of only three tasks: CHOOSE, ACTIVATE, and CREATE. Although decision types originating in other disciplines exist, we provide definitions for these tasks that are suitable for the visualization community. Our proposed typology offers two benefits. First, the ability to compose and hierarchically organize the tasks enables flexible and clear descriptions of decisions with varying levels of complexities. Second, the typology encourages productive discourse between visualization designers and domain experts by abstracting the intricacies of data, thereby promoting clarity and rigorous analysis of decision-making processes. We demonstrate the benefits of our typology through four case studies, and present an evaluation of the typology from semi-structured interviews with experienced members of the visualization community who have contributed to developing or publishing decision support systems for domain experts. Our interviewees used our typology to delineate the decision-making processes supported by their systems, demonstrating its descriptive capacity and effectiveness. Finally, we present preliminary findings on the usefulness of our typology for visualization design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Overview of advanced plasma-facing materials testing for Fusion Pilot Plants at DIII-D

Characterization and testing of advanced plasma-facing materials (PFMs) for Fusion Pilot Plants (FPP) is being conducted at the DIII-D National Fusion Facility through the ongoing two-year FPP Candidate Materials Thrust. Year one tested 17 novel materials utilizing the Divertor Materials Evaluation System (DiMES), with samples analyzed pre- and post-experiment via SEM, EDS, and confocal microscopy. Repeatable reference discharges were developed to ensure uniformity between experiments, including a new strike-point rastering scenario to provide more uniform heat/particle flux across DiMES during ELMing H-mode discharges. Various sample geometries and temperatures were used to achieve FPP-relevant conditions, including samples angled 10° towards the incident plasma flux and pre-heating up to 500 °C. The first exposure of liquid lithium (Li) capillary porous structures in a tokamak demonstrated uniform emission of Li vapor and suppression of Li droplets in H-mode when preheated to 350 °C. Dispersoid-strengthened W with 1 wt% TaC, TiC, and ZrC exposed to H-mode showed cracking and dispersoid ejection for all varieties except TiC, providing a clear down-selection. Ultra-high temperature ceramic materials TiB 2 and ZrB 2 showed minimal degradation under L-mode exposure. Silicon carbide (SiC) fiber composites showed arcing along edges, while CVD SiC remained pristine. Atmospheric plasma-sprayed W and SiC coatings endured H-mode exposure without macroscopic delamination; SiC exhibited granular ejection, while W showed increased outgassing. Additional W-based alloys were stress tested in H-mode, including Ni-based W heavy alloys, W f SiC f /W composites, W multi-principle element alloys, and functionally-graded W/SiC, to varying degrees of success.

DIII-D↗
DOE OSTI2024

Synthesis of Three Isoelemental MXenes and Their Structure–Property Relationships

The MXene family has rapidly expanded since its discovery in 2011 to include nearly 50 unique MXenes, not accounting for solid solutions and diverse surface terminations. However, a question raised since their discovery has been: What is the effect of n? In other words, how does the number of layers affect the MXene properties? To date, no direct study of the impact of n has been conducted due to the lack of isoelemental MXene compositions spanning more than two n values. Herein, we report on a system of three MXenes with identical M-site chemistries, (Mo 2/3 V 1/3 ) n+1 C n T x (n = 1, 2, and 3), allowing for the study of MXene structure–property relationships across n, for the first time. Chemical analysis of the samples shows complete and partial ordering of the M-elements in the n = 2 and 3 samples, respectively. We show that sample stability gradually evolves as n is increased from 1 to 3, while electronic and electrochemical properties exhibit more significant changes in going from n = 1 to 2 than from n = 2 to 3.

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

Customer-Focused Key Performance Indicators for Electric Vehicle Charging

To systematically improve the public charging experience, EV charging industry stakeholders need to define and measure it precisely. Many stakeholders currently measure aspects of the charging experience, but they typically employ metrics that are either operational in nature, such as charger uptime and mean time between failures, or composite customer satisfaction indices. To improve the customer experience most effectively, the industry needs metrics that define the charging experience from the perspective of the customer, not business operations. Furthermore, industry practitioners need granular metrics to know what specific aspects of the charging experience need improvement. This report defines such customer-focused metrics, called key performance indicators (KPIs).

33 - ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

CHESS 2025: Waveform LiDAR data from NEON AOP surveys

This dataset provides Level 1 (L1) full-waveform light detection and ranging (LiDAR) data collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). These data were acquired to enable characterization of vegetation structure and other three-dimensional features of the land surface, and to evaluate structural changes that may have occurred between a prior LiDAR acquisition in 2018 and the 2025 overflight. Waveform LiDAR data can provide more detailed information about objects on the ground than discrete point clouds typically do, and they are often used for granular target segmentation and characterization of subcanopy vegetation. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. LiDAR data were acquired using the Optech Galaxy Prime Airborne LiDAR Terrain Mapper onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). These are the primary waveform LiDAR data delivered by NEON and are provided per flightline in compressed Pulsewaves format, an open-source binary file standard. A Pulsewaves object comprises a two files: a pulse (.pls) file, which stores the geographic origin, outgoing vector, and metadata for every laser pulse emitted by the scanner, and a wave file (.wvs), which stores the sequential amplitude samples of the outgoing pulse and the returning signals. The files are published here in their compressed forms (.plz, .wvz). All waveform data were processed following the theoretical workflow described in the NEON L0-to-L1 Waveform LiDAR Algorithm Theoretical Basis Document (Krause and Goulden 2022a); however, the Pulsewaves output format differs from a legacy format described in that document. Waveform amplitude samples are recorded at 1 nanosecond intervals. All coordinates are provided in meters. Horizontal coordinates are referenced in Universal Transverse Mercator (UTM) zone 13N and the World Geodetic System (WGS) 1984 ensemble datum. Elevations are referenced to Geoid12A. Waveform data for the UPTA survey area were collected without incident and the published records are complete. However, both the ALMO and CRBU collections experienced issues that resulted in incomplete data for those areas. On collection day 2018-06-16 a hardware failure caused the waveform digitizer to lose data from the eastern edge of the ALMO site (Figure 22). The waveform data for flightlines 2–20 could not be extracted from the digitizer, and the data proved unrecoverable. As a result, a portion of the site does not have coverage with waveform data. Although no hardware failure was observed during collection over the CRBU area, final waveform files generated by vendor software contained only ~25% of the expected number of return pulses. After discovery, NEON initiated troubleshooting with the vendor. The root cause of the data ablation had not been identified at the time of publication. Additional data will be published in an update to this package if further recovery proves successful. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2024

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An Open-Source Framework for the Computational Analysis and Design of Autothermal Chemical Processes

The main project object was to develop software tools for simulating non-equilibrium autothermal processes, improving the prospects for identifying and designing such systems. The project demonstrates the use of these tools to simulate autothermal pyrolysis, a process recently developed at the pilot-plant scale at Iowa State University. In such process, instead of externally heating a reactor to pro-vide the enthalpy of pyrolysis, sufficient oxygen in the form of air is introduced into the reactor to support partial oxidation of reactants and products with the exothermic energy released supporting endothermic pyrolysis reactions. A fluidized bed is used to assure good mixing of biomass and oxidant and provide an isothermal reaction environment. The amount of oxygen required depends upon the kind of biomass being pyrolyzed and parasitic heat losses from the reactor. For example, for woody biomass pyrolyzed under conditions that simulate adiabatic operation, equivalence ratios can be as low as 0.06, compared to 0.20 or higher for autothermal gasifiers. By removing the heat transfer bottleneck of conventional pyrolysis, operation in autothermal mode allowed a significant increase in reactor throughput process, approaching five times the throughput of the conventionally operated pyrolyzer. Different simulation strategies were considered and developed: a zero-dimensional chemistry model was used to verify the applicability of kinetic schemes to predict biomass fast pyrolysis in autothermal conditions. Conventional chemical reactor models such as the plug flow reactor and the partially stirred reactor were used to investigate the role of mixing in the fluidized bed pyrolyzer and to establish the impact of mixing time on the gas-phase reactions. A comprehensive multiphase computational fluid dynamics (mCFD) framework, including polydisperse granular phase modeling and detailed chemical kinetics was formulated and used to model the experimental setup for autothermal biomass fast pyrolysis at ISU. Multiphase CFD was also used to investigate the role of biomass feed positioning on the mixing of biomass in the pyrolizer. Finally, a reduced order model (ROM), suitable to be implemented in process simulators was obtained. Both the mCFD and the ROM were validated against experiments.

09 BIOMASS FUELS↗
DOE OSTI2026

Hierarchical semi-Markov models with duration-aware dynamics for activity sequences

Residential electricity demand at granular scales is driven by what people do and for how long. Accurately forecasting this demand for applications like microgrid management and demand response therefore requires generative models for activities that can produce realistic daily activity sequences, capturing both the timing and duration of human behavior. This paper develops a generative model of human activity sequences using nationally representative time-use diaries at a 10-min resolution. We use this model to quantify which demographic factors are most critical for improving predictive performance. We propose a hierarchical semi-Markov framework that addresses two key modeling challenges. First, a time-inhomogeneous Markov router learns the patterns of “which activity comes next.” Second, a semi-Markov hazard component explicitly models activity durations, capturing “how long” activities realistically last. To ensure statistical stability when data are sparse, the model pools information across related demographic groups and time blocks. The entire framework is trained and evaluated using survey design weights to ensure our findings are representative of the U.S. population. On a held-out test set, we demonstrate that explicitly modeling durations with the hazard component provides a substantial and statistically significant improvement over purely Markovian models. Furthermore, our analysis reveals a clear hierarchy of demographic factors: Sex, Day-Type, and Household Size provide the largest predictive gains, while Region and Season, though important for energy calculations, contribute little to predicting the activity sequence itself. The result is an interpretable and robust generator of synthetic activity traces, providing a high-fidelity foundation for downstream energy systems modeling.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Characterization and Understanding of Thermal Transport in Stationary and Moving Particle Beds for Concentrating Solar Power

Concentrating solar power (CSP) coupled with thermal energy storage (TES) is being considered as an appealing solution to deliver stable, dispatchable, and inexpensive electricity generation from renewable solar energy. The heat from the HTF is stored in a TES system that is significantly less expensive and more scalable than electrochemical storage, enabling electricity generation when the sunlight is unavailable. The general concentrating solar thermal technologies are being pushed to higher temperature for higher efficiency and lower cost in a diverse range application in electricity generation (in CSP), industrial heating, and thermochemical fuel productions. Inert and abundant solid particles are becoming increasingly more important for these emerging applications. Heat transfer in both stationary and moving particle beds have been studied for many decades, revealing a great deal of insights on the mechanisms. Yet, there are still several unresolved problems from both fundamental understanding and technological applications points of view. In this Chapter, we will attempt to provide an up-to-date picture of thermal transport in stationary and moving particle beds. In both cases, we will review historical and contemporary developments in the theoretical, computational, and experimental aspects of heat transfer and related behaviors in particle beds. Here, we will discuss both the advancements and limitations of the existing studies and point to future developments towards deeper fundamental understandings and better engineering applications.

14 SOLAR ENERGY↗