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.

956 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗
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↗
NASA NTRS1977

Heat Flow and Thermal History of the Moon

An analysis was made of current heat flow data and thermal models of lunar evolution which satisfy the diverse information that has accumulated on internal processes.

Ye A Lyubimova↗
NASA NTRSDate not supplied

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗
NASA NTRSDate not supplied

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

vaned diffuser↗
NASA NTRSDate not supplied

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

centrifugal compressor↗
NASA NTRS2026

LAI Assimilation Schedules to Constrain Uncertain Cultivars and Soils in CERES-Maize

In anticipation of large-domain crop model applications where precise local configuration and calibration is not possible, we describe benefits and potential drawbacks of employing a crop pest module to achieve leaf area index (LAI) assimilation into a high performing CERES-Maize crop model configuration at a field experiment site in Perry, Iowa. Simulation experiments explore the use of MODIS satellite-derived LAI to constrain and adjust LAI to counter imprecise cultivar and soil configurations often occurring in the absence of high-quality local information. Simulations using single-day, window, and continuous LAI replacement across 5 cultivars and 2 soil calibration approaches for nine corn rotation years from 2004 to 2020 led to different yield outcomes and reverberations throughout the field environment. Evaluating variance and mean bias, results indicate minimal interventions in early vegetative and grain-filling stages were more beneficial than use of continuous LAI adjustments, as they minimized disruptions to the internal resource balances governing plant stresses and grain production. LAI adjustment was particularly helpful in constraining growth related to uncertain thermal unit requirements and leaf tip appearance rates (P1 and PHINT cultivar parameters, respectively). Findings underscore the need to assimilate additional state variables to ensure internal biophysical coherence. This approach shows promise for applications spanning wider domains with prediction time pressure where detailed configuration, more complex assimilation methods, or recalibration of crop model parameters may not be practical.

Phenology↗
NASA NTRS2026

A Summary of Test and Analysis Results from a Second Lift+Cruise Full-Scale Drop Test

The realization of advanced air mobility markets is enabling new forms of transportation to take shape in the United States and around the world. Though currently in development, as these markets mature, new types of vertical take-off and landing (VTOL) vehicles have been undergoing development for use. There are many factors which must be addressed prior to these types of vehicles becoming viable alternative forms of transportation in these markets. These factors include incorporation into the existing airspaces, the logistics of operating in urban environments, along with numerous factors associated with safety and reliability. To address some of the safety aspects associated with the development of these new types of vehicles, NASA has been conducting research into the performance of an example electric VTOL (eVTOL) aircraft as a part of the Revolutionary Vertical Lift Technology (RVLT) project. Over the course of this research, many aspects including the development of energy absorbing components, the evaluation of seating systems, the development of advanced finite element material model systems and the acquisition of full-scale vehicle impact data were investigated. The report will discuss aspects related to the acquisition of full-scale vehicle data which occurred in the form of a full-scale impact test conducted in the Summer of 2025. This test was on a NASA designed Lift+Cruise composite cabin test article and represented a partial capstone in the entirety of previous eVTOL research conducted for the project. In this test, a variety of experiments were included in order to investigate the effect of a full-scale environment on the experiment results. In parallel, the development of a computational impact model to simulate the full-scale test will be discussed in this report. A model of the Lift+Cruise test article was developed utilizing data collected from previous sub- and full-scale test data and then simulated in the current test environment. The model development, its use in pre-test predictions, and its use in post-test correlation will all be presented. This report will present the test data acquired from the Lift+Cruise test and document several of the results obtained. One intended result is to determine the effect of a complex full-scale crash impact on the identification of occupant injury risk within seat and vehicle designs. A second intended result is to determine whether high-fidelity models can be used with some confidence in the prediction of test events and can allow for additional test cases to be simulated without the need of having to conduct additional tests. The overall goal of the test is to provide the community with data that can be used for design, development or certification efforts, along with providing data on what an example eVTOL crash incident could entail.

energy storage systems↗
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↗
DOE OSTI2024

In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back

Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.

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

Value of Geothermal Energy Storage for Supply-Side and Demand-Side Applications

This report presents the results of a study examining the value potential for geothermal energy storage (GES), a long-duration energy storage resource that stores thermal and/or geomechanical energy in the subsurface. GES could benefit the overall U.S. power system by temporally shifting electricity generation (supply-side) or meeting building heating and cooling load (demand-side). This report analyzes supply-side and demand-side opportunities independently because of differences in applications and models. Currently there is significant uncertainty about the development costs for GES, with only a limited number of demonstration plants for electric energy storage and building heating and cooling storage developments. In this report, we estimate the value of supply-side and demand-side GES to the bulk power system in the contiguous United States. Because of the significant uncertainty about GES development costs, this analysis does not consider GES deployment costs but instead focuses on the value of GES to the U.S. electricity system. The estimated values of GES provide reference points for economically competitive commercial cost targets. Supply-side GES is modeled as part of an enhanced geothermal system (EGS) generation plant in NREL's Regional Energy Deployment System (ReEDS) capacity expansion model (Ho et al. 2021). In contrast to conventional geothermal plants, which generate constant power, EGS plants have unique features that may allow for in-reservoir energy storage for flexible generation. Demand-side GES for heating and cooling, including seasonal hot and cold storage and short-duration heat pump storage, is incorporated into a price-taker model using Cambium electricity marginal cost projections. To establish an upper bound for the value of GES, analysis focused on favorable scenarios for storage with high generation from zero marginal cost, variable renewable energy resources. High penetrations of variable renewable energy generation can increase hourly electricity price variability, which increases the value of temporal energy arbitrage for storage technologies like GES.

15 GEOTHERMAL ENERGY↗
NASA NTRS2026

Ice in Convective Storm Environments: Assessing Impacts of Microphysics on Simulated Reflectivities

Satellite radar retrievals of clouds and precipitation rely on assumptions about the microphysical properties of hydrometeors such as particle size distributions (PSDs) and massdimensional relationships, which have been shown to vary in distinct environments. A deeper understanding of linkages between mass-dimensional relationships, ice crystal shape, and radar reflectivities across cloud and precipitation regimes is necessary for developing new retrievals for upcoming satellite radar missions, including the INvestigation of Convective UpdraftS (INCUS), as well as future missions that will feature space-borne radars observing clouds, convection, and precipitation. This study investigates ice microphysical characteristics through aircraft in-situ and remote sensing observations from two convective events during the Mid-latitude Continental Convective Clouds Experiment (MC3E) in the Southern Great Plains, and tests implications for ice particle habit in forward modeled radar observations. During MC3E, the Ka/Ku- band High Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) was onboard the NASA ER-2 aircraft, while in-situ measurements were taken onboard the University of North Dakota Citation. In-situ PSDs and particle shape information from optical array probe data are used to forward model radar reflectivities from particle scattering databases and are compared against HIWRAP reflectivities. Various particle combinations can be used to find agreement with HIWRAP observations, and the forward modeled reflectivities were sensitive to particle type. Varying the riming on the aggregate and graupel particles resulted in larger than 10 dB differences in reflectivities. The in-situ and radar observations imply the occurrence of aggregation, size sorting and lofting particles.

Julia A. Shates↗
DOE OSTI2025

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials↗