Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mathematical Models”

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.

544 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Recent Progress in the Electroweak Structure of Light Nuclei Using Quantum Monte Carlo Methods

Nuclei will play a prominent role in searches for physics beyond the Standard Model as the active material in experiments. In order to reliably interpret new physics signals, one needs an accurate model of the underlying nuclear dynamics. In this review, we discuss recent progress made with quantum Monte Carlo approaches for calculating the electroweak structure of light nuclei. We place particular emphasis on recent β decay, muon capture, neutrinoless double β decay, and electron scattering results.

Physics↗
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
DOE OSTI2026

Feedback, physics, and forecasts: The emerging paradigm of machine learning-driven battery research

Machine learning (ML) is reshaping how we understand, predict, and optimize electrochemical systems. In batteries, ML accelerates discovery across chemistry, design, and operation by transforming massive experimental and simulated datasets into predictive, interpretable models. This review consolidates a decade of progress in ML-driven battery innovation, from early-cycle feature extraction to operando image analysis and physics-informed modeling. We categorize approaches by data domain and physical fidelity, emphasizing interpretable ML for diagnostics, reinforcement learning for control, and multi-objective optimization for lifetime extension strategies. Additionally, we demonstrate how integrated models accelerate discovery, reduce testing time, and guide sustainable design. Economic analyses furthermore illustrate how these advances can lower cost per cycle and improve circularity. Together, these developments chart a path toward self-optimizing, sustainable battery technologies.

artificial intelligence↗
DOE OSTI2026

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
DOE OSTI2024

Understanding and Predicting the Spatially Resolved Adsorption Properties of Nanoporous Materials

Using knowledge from statistical thermodynamics and crystallography, we develop an image–image translation model, called SorbIIT, that uses three-dimensional grids of adsorbate–adsorbent interaction energies as input to predict the spatially resolved loading surface of nanoporous materials over a broad range of temperatures and pressures. SorbIIT consists of a closed-form differential model for loading-surface prediction and a U-Net to generate spatial differential distributions from the energy grids. SorbIIT is trained using the energy grids and adsorbate distributions (obtained from high-throughput simulations) of 50 synthesized and 70 hypothetical zeolites and applied for predicting the adsorption of carbon dioxide, hydrogen sulfide, n-butane, 2-methylpropane, krypton, and xenon in other zeolites from 256 to 400 K. In conclusion, employing a quadratic isotherm model for the local differentiation, SorbIIT yields mean R 2 values of 0.998 for total adsorption and 0.6904 for local adsorption with a resolution of 0.2 Å, and a value of 0.721 for the structural similarity of the local loading distribution.

Sun, Yangzesheng [Univ. of Minnesota, Minneapolis,↗
DOE OSTI2025

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗
DOE OSTI2026

Getting Warmer: IceCube Nears Freeze Out

IceCube has recently detected a diffuse population of high-energy neutrinos arising from the Milky Way. We use this high-significance detection to place new limits on dark matter (DM) annihilation to neutrinos with two complementary approaches. The first method uses the background-subtracted Galactic longitude distribution of shower events to place a conservative bound on the DM annihilation cross section that does not rely on any assumed Galactic cosmic ray emission model; the resulting limits on the velocity-averaged annihilation cross section improve upon existing bounds by factors of a few. The second method uses the template-dependent neutrino energy spectra from the Inner Galaxy, inferred under different Galactic cosmic ray emission models. This complementary approach shows that the inferred Galactic neutrino intensities are already sensitive to DM contributions near the thermal-relic benchmark for a range of TeV-scale DM masses, though this comparison is more model-dependent. Our results demonstrate that measurements of diffuse Galactic neutrino emission can be used as a powerful probe of DM annihilation into neutrinos. Future observations with IceCube-Gen2 and KM3NeT will substantially extend this sensitivity, potentially allowing a decisive test of the thermal freeze-out mechanism with Galactic neutrino observations.

Mukhopadhyay, Mainak [Fermilab; Chicago U., KICP; ↗
DOE OSTI2026

EXERGETIC: De-Risking Next-Generation Resilient Geothermal Hybrids via At-Scale Evaluation Using Virtual Emulation Digital Twin Environment for Efficient Operation

The DOE-GTO-funded project, award number 5.1.2.12, entitled "EXERGETIC - De-risking Next Generation Resilient Geothermal Hybrids via at-Scale Evaluation Using a Virtual Emulation Digital Twin Environment for Efficient Operation," advances the solution to these challenges by developing and validating a geothermal co-emulation environment implemented at the National Laboratory of the Rockies (NLR)'s Advanced Research on Integrated Energy Systems (ARIES) platform. This framework enables the de-risking of next-generation geothermal and geothermal hybrid systems through high-fidelity modeling, real-time digital emulation, advanced control strategies, and techno-economic assessment. The project focused on geothermal hybrid configurations that integrate geothermal power plants with concentrated solar power and underground thermal energy storage, enabling enhanced efficiency, flexibility, and grid support capabilities. The main goal of this project was the development of a geothermal digital co-emulation environment to demonstrate the technical and economic value of geothermal hybrid systems and their contribution to grid stability and flexibility. The EXERGETIC framework combined physics-based models, controls, and real assets at ARIES, including digital real-time simulators (DRTS), a 20-MW-scale controllable grid interface (CGI), and a 2-MW conventional generator. Detailed transient models were developed for the key subsystems of a hybrid geothermal plant, including parabolic trough solar collectors, reservoir thermal energy storage (RTES), and a binary Organic Rankine Cycle (ORC) power plant. The ORC model explicitly captured thermal inertia and off-design operation and integrated control strategies to dynamically respond to electric load profiles. The models were validated against published experimental and numerical studies, demonstrating strong agreement and confirming the accuracy and robustness of the modeling approach. The resulting digital twin represents geothermal-solar-storage systems at multiple scales (1 MW to 100 MW) and enables realistic emulation of grid-connected operation. The control architecture allows the geothermal resource to provide stable baseload generation, while solar and stored thermal energy supply flexible, dispatchable support during periods of high demand or variable grid conditions. A key contribution of the EXERGETIC project is the demonstration that geothermal hybrid systems can be designed to be active grid assets rather than passive baseload generators. Using the ARIES platform, the digital twin was evaluated under multiple grid scenarios, including load following, voltage support at the distribution level, and frequency response at the transmission level. Results show that hybrid geothermal systems can respond effectively to dynamic grid conditions, providing inertia-like behavior, primary frequency support, and voltage regulation through coordinated control. In addition to the performance and grid services capability analysis of geothermal and hybrid geothermal systems, the EXERGETIC project also focused on scalability and techno-economic analysis of geothermal hybrid plants. In particular, for the scalability analysis, machine-learning (ML)-based surrogate models were trained using data generated from the geothermal digital twin under different grid-connected scenarios and plant capacities. These ML models demonstrated strong interpolation and extrapolation capabilities across plant sizes, accurately reproducing both steady-state and transient responses with very low errors. Regarding the techno-economic analysis, plant performance results were integrated with cost models for hybrid geothermal systems, and the levelized cost of electricity (LCOE) was used as the main economic metric to evaluate system performance across a range of system capacities, solar shares, solar multiples, and storage durations. Results indicate that economies of scale significantly reduce geothermal LCOE as plant capacity increases, with large-scale systems (25-100 MW) achieving substantially lower costs than small plants. Hybridization with solar thermal energy and storage further improves economic performance by increasing capacity utilization and enabling flexible dispatch. In addition, thermal storage plays a critical role in reducing LCOE by maximizing geothermal, solar, and stored energy resources. In summary, the results from this project demonstrate that geothermal hybrid systems represent a promising alternative for increasing the energy conversion efficiency of geothermal technologies, contributing to the preservation of geothermal resources, and supporting the transition of geothermal plants from traditional baseload resources into flexible, resilient, and cost-competitive energy conversion technologies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dark Energy Survey year 6 results: Magnification modeling and its impact on galaxy clustering and galaxy-galaxy lensing cosmology

Gravitational lensing magnification alters the observed spatial distribution of galaxies and must be accounted for to prevent biases in cosmological probes of the large-scale structure. We investigate its effects on the Dark Energy Survey Year 6 galaxy clustering and galaxy-galaxy lensing analyses using the fiducial lens (position tracer) sample M ag L im++. Magnification bias is parameterized by a coefficient that describes the response of the number of selected objects per unlensed area element to a change in the lensing convergence. We quantify this coefficient using the BALROG synthetic source injection catalog to account for the complexity of the selection function, and compare these results with simplified estimates. The resulting values of the magnification coefficients for each redshift bin are [3.16 ± 0.08, 2.76 ± 0.21, 4.09 ± 0.15, 4.42 ± 0.16, 4.90 ± 0.29, 4.83 ± 0.25]. Relative to Year 3, this analysis provides more precise and accurate magnification bias estimates through a larger BALROG area and reweighting to better match the data properties. Here, the cosmological results are robust when tested against various magnification parameter prior choices and also when adding cross-clustering between lens redshift bins. Neglecting magnification, however, introduces significant systematic shifts: relative to the fiducial analysis with Gaussian priors centered on the BALROG -derived estimates, we observe shifts of 1.37σ in S 8 and -0.84σ in Ω m (with cosmic shear included: -0.61σ in S 8 and -0.71σ in Ω m ), in agreement with findings from simulated data, demonstrating that magnification must be modeled to avoid biases. Freeing the magnification bias in lens bin 2 leads to unphysical negative values, further justifying its exclusion from the fiducial Year 6 analysis.

Cosmological parameters↗
DOE OSTI2024

Collective excitations and low-energy ionization signatures of relativistic particles in silicon detectors

Abstract Solid-state detectors with a low energy threshold have several applications, including searches of non-relativistic halo dark-matter particles with sub-GeV masses. When searching for relativistic, beyond-the-Standard-Model particles with enhanced cross sections for small energy transfers, a small detector with a low energy threshold may have better sensitivity than a larger detector with a higher energy threshold. In this paper, we calculate the low-energy ionization spectrum from high-velocity particles scattering in a dielectric material. We consider the full material response including the excitation of bulk plasmons. We generalize the energy-loss function to relativistic kinematics, and benchmark existing tools used for halo dark-matter scattering against electron energy-loss spectroscopy data. Compared to calculations commonly used in the literature, such as the Photo-Absorption-Ionization model or the free-electron model, including collective effects shifts the recoil ionization spectrum towards higher energies, typically peaking around 4–6 electron-hole pairs. We apply our results to the three benchmark examples: millicharged particles produced in a beam, neutrinos with a magnetic dipole moment produced in a reactor, and upscattered dark-matter particles. Our results show that the proper inclusion of collective effects typically enhances a detector’s sensitivity to these particles, since detector backgrounds, such as dark counts, peak at lower energies.

Physics↗
DOE OSTI2025

Field Validation of Thermoelectric Generation System at Holcim Cement Plant in Alpena, Michigan

Executive Summary Project Background The Industrial Technology Validation (ITV) program aims to identify and demonstrate the performance of new, emerging, and underutilized energy-saving technologies in the industrial sector to help inform decisions to help accelerate their commercialization and deployment, as well as to help make industries more competitive. This ITV demonstration evaluated a thermoelectric generation (TEG) technology at a cement plant, aiming to reduce energy demand in the cement industry. A median cement plant consumes 5.73 million British thermal units per ton of clinker production (resulting in 0.838 metric tons of carbon dioxide [CO₂] emissions per ton of clinker) (Boyd and Zhang 2011, EPA 2021), equivalent to approximately 6.9 trillion British thermal units (TBtu) per year in energy consumption at a cement plant producing 3,300 tons of clinker per day.¹ Collaborating with Holcim, Advanced Thermovoltaic Systems (ATS) developed and deployed a pilot-scale thermoelectric power system to efficiently capture and convert waste heat to electricity. The system leverages the Seebeck effect to convert temperature differences on two sides of semiconductor cartridges into electrical power (ScienceDirect, n.d.). This generation is realized with minimal moving parts compared to existing waste-heat-to-generation solutions and allows capture from heat sources with temperatures as low as 150°C. This project aimed to validate a scalable solution applicable for capturing medium-temperature waste heat, including ambient losses from other high-temperature processes, and high-temperature sources less suitable for other waste-heat-to-power solutions. By recovering this otherwise wasted heat, this project intends to validate improvements to overall process efficiency through reduction in purchased electricity, thereby reducing operational costs while enhancing resiliency and competitiveness. Description and Scope This study evaluated the performance of a TEG system from ATS as a solution to convert waste heat into useful power at a Holcim cement plant in Alpena, Michigan. This plant is a fully integrated cement plant that has been operating since 1907. The facility operates continuously (24/7/365) with approximately 250 employees and five long dry kilns, yielding a total production capacity of 7,852 tons of cement per day (EPA 2023). Currently, the Alpena plant uses waste heat boilers to convert waste heat from the exhaust of each kiln into steam, which drives steam turbine generators. The ATS TEG is being evaluated for its potential to supplement the steam turbines by capturing the remaining lower grade heat. This technology is also being considered for other Holcim plants where steam turbines are not a viable option. ATS installed a pilot-scale TEG unit with an array of 582 individual thermoelectric semiconductor cartridges, of which 573 were operational. The cartridges are sandwiched between 48 hot plates and 49 cold plates. Each cartridge is designed to generate 20 watts (W) of gross power at a hot-side temperature of 240°C and cold-side temperature of 20°C. As such, the total gross generation capacity of the installed system is 11.5 kilowatts (kW) at design conditions. The system configuration for the evaluation was designed to prioritize convenience of installation and minimize disruption to production at the site, while ensuring that the heat required can be obtained for evaluating the TEG system at various operational conditions. To accomplish this, a portion of the steam supplied to Alpena’s steam turbine generation system was diverted to be used as the heat source for the TEG system, while water was supplied to the cold side of the system from nearby Lake Huron. This configuration was designed for the evaluation of the pilot-scale system to assess the performance at different conditions. A commercial-scale system will likely vary from the pilot system depending on typical configurations, including both scale and application. Future commercial applications of the ATS system would involve integrating the system into the exhaust from kiln preheaters, clinker coolers, or radiant heat capture from kiln shells for the heat source. For the cold source, a range of cooling solutions can be considered, including a mechanical cooling system, depending on the location and the application. To increase the generation capacity for commercial applications, the technology provider is working toward developing a commercial-scale TEG system, which would combine multiple TEG units (each similar in design to the pilot system) together. The scope of this evaluation includes the pilot-scale TEG system and all impacted equipment including pumps, controllers, and power handling equipment. Study Objectives The evaluation's goal was to assess the potential of the ATS TEG system to generate useful electrical power by capturing waste heat from cement production kilns. The objectives of this study are to evaluate and verify the following claims made by ATS regarding the pilot-scale system installed at the Holcim Alpena plant. The following design parameters and claims are also outlined in Table ES- 1 and Table ES- 2: • Gross Power: The thermoelectric system converts heat into power to create gross power, the total measured power generated by the system. The 573 active cartridge pilot-scale system is expected to generate 11.5 kW of gross power at the designed hot-side temperature of 240°C and cold-side temperature of 20°C. Power production is dependent on the temperature difference between the heat source (ultimately from the waste heat) and cold temperature supply source. • Net Power: The net power is the total usable power provided to the site by the TEG system after deducting parasitic power loads from the gross generated power. Supplementary equipment is required to operate the TEG system including pumps, controllers, and, in certain anticipated applications, mechanical cooling, which introduce parasitic loads to system operation. After deducting the parasitic loads from the gross power generation, ATS anticipates achieving a net power generation of 7.5 kW from the pilot-scale system. • Thermal Efficiency: The thermal efficiency is the percent of the total heat transferred to the TEG system that is converted to gross power. Historically, TEGs have a thermal efficiency of 2%–5% (DOE 2008). Prior industrial-scale TEG systems, such as the E1 TEG offered by Alphabet Energy, operated at an efficiency of 2.5% (Lamonica, 2014). ATS anticipates achieving an average efficiency of 4.8% or higher in converting heat energy to usable electricity. • Cartridge Performance: The TEG system comprises 573 active individual semiconductor cartridges, each of which generates a portion of the total power. Cartridge optimization and selection is an important design consideration for potential future TEG system design performance. Therefore, understanding the distribution of gross power and efficiency within the pilot system is vital to understanding what is achievable. At a design hot-side temperature of 240°C and cold-side temperature of 20°C, ATS anticipates a cartridge performance of 20 W of gross power per cartridge at an efficiency of 4.8% per cartridge. In addition to evaluating the claimed performance of the TEG pilot-scale unit, the study estimated the potential annual impacts of a scaled-up commercial system used to capture kiln waste heat over annual operations. The evaluation estimated the gross and net annual electric generation achievable by capturing heat from the two proposed tap-in points: the kiln exhaust and the clinker cooler exhaust; see Section 2.1 for details. Two use cases were examined: • Holcim Alpena: The Holcim Alpena site consists of long dry kilns with superheater boilers, which differs from the rest of Holcim’s cement plant portfolio and results in lower waste heat temperatures. The study estimates gross and net annual generation using the superheater boiler exhaust and clinker cooler exhaust, based on 2023 operational data. • Typical Installation: Common cement plants have preheater kilns with higher exhaust temperatures than Holcim Alpena across a range of production rates. The study estimates gross and net annual generation using the preheater exhaust and clinker cooler exhaust, with a sensitivity analysis to account for the typical range of preheater exhaust temperatures, clinker cooler exhaust temperatures, and clinker production rates. Methodology The evaluation methodology followed a measurement and verification (M&V) strategy based on the International Performance Measurement and Verification Protocol Option B through comprehensive measurements and analyses of the affected systems. Evaluation data was collected from March 9 to March 11, 2024, the test period of the pilot TEG system. During the test period, in coordination with the ITV team, the ATS team adjusted system operations to capture the range of variability expected for each of the variables pertinent to performance of the system. The methodology consisted of two parts: evaluating the performance of the pilot unit's TEG system and estimating the annual TEG impact in terms of gross and net power based on a given waste heat profile. First, the evaluation of the thermoelectric generation performance of the pilot unit relative to the claims was performed by analyzing the collected test data. Gross power of the pilot TEG system was directly measured. Net power was determined by deducting the measured parasitic power from the gross power. The gross power generation was compared to heat transferred to the system by the working fluid (which was heated by steam generated from the kiln waste heat) to calculate the thermal efficiency achieved by the system. Performance of individual semiconductor cartridges within the pilot array was also assessed in terms of measured gross cartridge power and calculated cartridge thermal efficiency. The second part of the evaluation estimated the annual TEG impacts in terms of gross power and net power (calculated from the difference between gross power and parasitic power). This analysis comprised development of mathematical regression models for gross power and parasitic power, with assessment of each model’s goodness-of-fit characteristics to ensure satisfaction of statistical requirements. The models predicted the gross power generation, the parasitic load based on the temperature difference between the hot working fluid and the cold-side fluid (cold water from Lake Huron) entering the system, the volumetric flow rate of the cold-side fluid at the inlet, and the volumetric flow rate of the hot working fluid at the inlet. The annual impact analysis considered a theoretical commercial-scale system sized to capture the available waste heat at a cement plant, consisting of linked pilot-scale units that receive heat from a theoretical gas-to-working-fluid heat exchanger. To estimate annual impacts at the Alpena plant, the gross power and parasitic power regression models were applied to the arrays in the theoretical commercial-scale system. The heat supplied to the unit was calculated based on the kiln run time, annual production, kiln exhaust waste heat, and clinker cooler waste heat derived from 2023 Holcim Alpena kiln operational data. Net power impacts were calculated by deducting the resulting parasitic power from the estimated gross power. Inputs for the model were generated from a combination of hourly data, assumed design considerations for TEG system scale-up from the pilot-scale unit, and assumptions regarding TEG system operations. This analysis was then used as the basis for estimating annual impacts of typical TEG installation at cement plants, by applying sensitivity analyses to key kiln operational characteristics including kiln preheater exhaust temperatures, cooler clinker exhaust temperatures, and plant daily production rates across a range of expected values. Project Results/Findings Table ES- 2 and Table ES- 2 provide a summary of the operating conditions and evaluation results compared to the stated claims from the technology provider. Key takeaways include: • Gross Power: The peak gross power achieved during the testing period was 10.0 kW, compared to the 11.5 kW expected for 573 active cartridges. The claimed gross power was associated with a target hot side of 240°C; however, the system only received a maximum hot-side mean plate temperature of 212°C during the testing period. • Net Power: The pilot-scale unit exceeded the claims for net power, achieving a peak of 7.7 kW net compared to a claim of 7.5 kW. One factor contributing to the higher achieved net power is the relatively high water pressure available through Lake Huron. The pilot TEG system did not require cold-side pumps during the test, whereas most installations would. This reduced the parasitic loads on the system, ultimately contributing to higher net power relative to the gross power. • Thermal Efficiency: The pilot-scale unit outperformed the claimed efficiency, achieving a peak system efficiency of 5.0% thermal efficiency compared to the stated 4.8%. • Cartridge Performance: To compare cartridge performance against claims, the study focused on the third day of testing, which aimed for conditions closest to the design specifications, with a hot side of 240°C and cold-side exit temperature of 6.4°–30°C. On this day, the mean gross power observed in the cartridges within the TEG array was 18.1 W/cartridge, and the peak performance was 34.7 W/cartridge. The estimated mean cartridge efficiency was 5.2%, and the estimated efficiency at peak gross cartridge power was 10%. The regression models developed for gross power generation and parasitic loads were used to estimate the generation impact for given heat input to the TEG from the working fluid (captured from the waste heat) and from the cold loop (Lake Huron) on an hourly basis for a year of operation. Based on this analysis, installation of a commercial-scale TEG system at the Holcim cement plant in Alpena, Michigan, with a waste heat exchanger of 0.85 effectiveness, would generate up to 391 kW of net power, translating to between 920,000 and 1,800,000 kilowatt-hours (kWh) in net electricity per year. Based on typical grid emissions for Alpena, this would avoid estimated net emissions by 752 metric tons of CO₂ annually.² The sensitivity analysis estimated that typical TEG system installations at cement plants could generate an average of 56–1,040 kW of net power, or between 488,000 and 9,110,000 kWh of net energy. This generation potential is most significantly affected by plant production rates and also influenced by preheater and clinker cooler exhaust temperatures. Applying the national average emission rate, typical commercial-scale installations at Holcim plants are projected to avoid between 182 and 3,401 metric tons of CO₂ annually per site. Table ES- 3 shows a summary of the estimated annual impacts.³ While parasitic loads are significant and vary by application, this analysis assumed the use of heating loop pumps and access to Lake Huron as a cold sink. This setup assumed no need for cooling loop pumps due to the available water pressure at the test site. Applications that require cooling towers or additional equipment are likely to experience higher parasitic loads. Therefore, the study’s estimates are most applicable to scenarios with similar parasitic load configurations—namely, access to a high-pressure cold sink. Applicability to other locations may be limited, as differing conditions could necessitate additional pumps and cooling systems, potentially impacting performance significantly.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Phonon screening and dissociation of excitons at finite temperatures from first principles

The properties of excitons, or correlated electron–hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron–hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio GW -Bethe–Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and SrTiO 3 . We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron–hole pairs.

Science & Technology - Other Topics↗
DOE OSTI2024

A novel digital lifecycle for Material‐Process‐Microstructure‐Performance relationships of thermoplastic olefins foams manufactured via supercritical fluid assisted foam injection molding

Abstract This research significantly enhances the applicability of thermoplastic olefins (TPOs) in the automotive industry using supercritical N 2 as a physical foaming agent, effectively addressing the limitations of traditional chemical agents. It merges experimental results with simulations to establish detailed material‐process‐microstructure‐performance (MP2) relationships, targeting 5–20% weight reductions. This innovative approach labeled digital lifecycle (DLC) helps accurately predict tensile, flexural, and impact properties based on the foam microstructure, along with experimentally demonstrating improved paintability. The study combines process simulations with finite element models to develop a comprehensive digital model for accurately predicting mechanical properties. Our findings demonstrate a strong correlation between simulated and experimental data, with about a 5% error across various weight reduction targets, marking significant improvements over existing analytical models. This research highlights the efficacy of physical foaming agents in TPO enhancement and emphasizes the importance of integrating experimental and simulation methods to capture the underlying foaming mechanism to establish material‐process‐microstructure‐performance (MP2) relationships. Highlights Establishes a material‐process‐microstructure‐performance (MP2) for TPO foams Sustainably produces TPO foams using supercritical (ScF) N 2 with 20% lightweighting Shows enhanced paintability for TPO foam improved surface aesthetics Digital lifecycle (DLC) that predicts both foam microstructure and properties DLC maps process effects & microstructure onto FEA mesh for precise prediction

Engineering↗
DOE OSTI2026

Upper ocean mixing, surface heat fluxes, and heat content variability in the upper 150 m during Hurricane Laura (2020)

Observations (from thermistor sensors and a hurricane glider) and model-based temperature were used to answer questions including the following: 1) What physical and environmental conditions explain the rapid intensification of Hurricane Laura? 2) How did the pre-existing warm mixed layer at Stone mooring (StM) influence the degree of cooling in the mixed layer? 3) How did vertical mixing, surface heat fluxes, and advective processes collectively shape the mixed layer heat evolution at StM? Observations measured by thermistor sensors showed that Hurricane Laura induced a cooling of -1.2oC in the mixed layer when Hurricane Laura’s wind speed increased to 47 m/s on 26 August. This contrasts with the -1.04oC mixed layer temperature change estimated from model-based temperature. At StM, the presence of a 31oC warm mixed layer and elevated heat content (60-80kJ/cm2) effectively preconditioned the upper ocean ahead of Hurricane Laura’s passage. After Hurricane Laura passed, turbulent mixing associated with the storm transported the heat anomalies downward into subsurface layers. In addition, this study shows that a surface heat flux of 5.04 kJ/cm2 supported the intensity of Hurricane Laura as it traversed over the StM. Using a mixed layer heat balance model, this study shows that entrainment flux and surface flux mostly contributed to the observed change in the mixed layer temperature.

Tsai, Senam↗