Search NASASearch

SEARCH · Search NASA

Results for “surrogate 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.

930 records · Page 47

Selectivity mechanisms of ion intercalation in Prussian blue analogs

Prussian blue analogs (PBAs) are a family of materials with facile, reversible, and selective ion transport capability for various ions via electrochemical intercalation, owing to their vacancy structure. The large tunable compositional space of PBAs allows for manipulation of intercalation behavior and selectivity by controlling structural vacancy level through choice of transition metal centers and modifications to the synthesis process. However, a lack of understanding of the mechanisms of ion selectivity hinders the material’s design process. Here, for this work, we investigated the origins of ion selectivity using a model PBA, copper hexacyanoferrate, and focused on eight technologically and biologically prominent ions, for which we determined a sequence of selectivity: Rb + > K + > Na + > Ba 2+ > Sr 2+ ≈ Ca 2+ > Mg 2+ > Li + . We provide electrochemical, structural, and redox evidence of strong correlation between the ion identity, the dominant charge-compensating redox, and preferred occupancy site. Specifically, using synchrotron anomalous X-ray diffraction (AXRD), we reveal that monovalent ions exhibit significant association with the corner sites of the unit cell and iron redox, whereas divalent ions display affinity toward the center site with higher ratios of copper redox. Informed by selectivity results, we applied CuHCFe to Li purification and achieved 99.9% purity. Our findings demonstrate an approach to elucidating ion intercalation behavior in order to distinguish and manipulate material properties to optimize separation performance.

Prussian blue analog

Co-Location of Cellulosic Bioethanol and Alcohol-to-Jet (ATJ) Production Facilities for Targeted Scale-Up of Sustainable Aviation Fuel (SAF) Production

Achieving aerospace industry net-zero emissions by 2050 requires rapid scaling of sustainable aviation fuel (SAF) production. Leveraging existing infrastructure, proven technologies like Alcohol-to-Jet (ATJ), and low carbon intensity (CI) feedstocks (e.g., switchgrass and miscanthus) can support this transition and help achieve near-term emissions reduction targets. This study evaluates the implications of lignocellulosic ethanol biorefinery siting and integration with petroleum refineries to produce SAF across 1000 sites randomly sampled from areas suitable for perennial grasses in the U.S. rainfed region. To better understand the logistics of material transport and handoffs, we integrated models of biomass harvest, transport, ethanol, and ATJ production in a stochastic framework based on Monte Carlo simulations to characterize SAF minimum selling price (MSP) and carbon intensity (CI), considering site-specific parameters (e.g., feedstock production, transportation, taxes, incentives). The results indicate trade-offs between MSP and CI across locations, with median MSP ranging from 7.9 to 12.8 USD·gal −1 and CI from −9.7 to 39.4 gCO 2 e·MJ −1 . Despite high estimated decarbonization costs (580 USD·tonCO 2 e −1 ), our results indicate that site-specific deployment of ATJ with low-CI feedstocks can improve sustainability outcomes. The framework provides a systematic approach to assess cost and sustainability trade-offs across locations, considering the end-to-end supply chain and supporting an informed investment in SAF production.

09 BIOMASS FUELS

Laser Beam Welding Advancements for In-Space Servicing, Assembly, and Manufacturing

NASA is currently working to develop in-space servicing, assembly, and manufacturing (ISAM) capabilities for low Earth orbit and the lunar surface. One crucial technology for this effort is laser beam welding. Laser systems can perform joining, cleaning, cutting, and repair activities, which will enable the construction of large in-space structures that could not fit on a single launch vehicle, such as trusses for solar panels, radiators, or communications infrastructure. Multiple projects studying laser welding for space applications are currently underway at NASA Marshall Space Flight Center. One of these, the DIsk-Shaped Configurable and Modular vAcuum uNit (DISCMAN), is a compact vacuum chamber designed to support parameter development for laser welding in microgravity. It contains a rotating platen with weld samples made from aluminum, steel, and titanium, a high-power infrared laser, and integrated pumps for pulling vacuum inside the sample cartridge. The DISCMAN payload is planned to launch to the International Space Station, where welds will be performed under sustained microgravity inside the Bishop Airlock. Another effort underway at Marshall is the Lunar Assembly and Servicing by Autonomous Robotics (LASAR) initiative. This project uses a space-rated robotic arm equipped with a laser weld head, wire feeder, and multiple cameras to perform welds in a thermal vacuum chamber simulating the lunar surface environment. Some of these welds are done on snowflake joints, which are specially designed to slot together to join segments of trussN structures, allowing for the construction of tall surface infrastructure. DISCMAN, LASAR, and other projects are being carried out to advance the technological maturity of in-space laser beam welding, collect data to inform computational models, and learn reliable processes for creating weld joints in space. This work supports NASA’s greater goals to expand humanity’s presence in low Earth orbit, establish a permanent moon base, and eventually send crewed missions to Mars and beyond.

Manufacturing

The Evolution of the Moon and the Terrestrial Planets

The thermal evolutions of the Moon, Mars, Venus, and Mercury are calculated theoretically starting from cosmochemical condensation models. An assortment of geological, geochemical, and geophysical data are used to constrain both the present-day temperatures and the thermal histories of the planets' interiors. Such data imply that the planets were heated during or shortly after formation and that all the terrestrial planets started their differentiations early in their history. Initial temperatures and core formation play the most important roles in the early differentiation. The size of the planet is the primary factor in determining its present-day thermal state. The Moon, smallest in size, is characterized as a differentiated body with a crust, a thick solid mantle, and an interior region which may be partially molten. It is presently cooling rapidly and is relatively inactive tectonically. Mercury, which probably has a large core, may have a 500-km-thick solid lithosphere and a partially molten core, if it is assumed that some heat sources exist in the core. If this is not the case, the planet's interior temperatures are everywhere below the melting curve for iron. The thermal evolution is dominated by the core separation and the high conductivity of iron that makes up the bulk of Mercury. Mars, intermediate in size, is assumed to have differentiated an Fe-FeS core. While the formation of an early crust is evident, large-scale melting and differentiation of the mantle silicates has occurred at least up until 1 b.y. ago. Present-day temperature profiles indicate moderate tectonic activity at the present time. Venus is characterized as a planet not unlike the Earth in many respects. Core formation has occurred probably during the first billion years after the formation. The present-day temperatures indicate a partially molten upper mantle overlain by a 100-km-thick lithosphere and a molten Fe-Ni core. We can expect that today Venus may have tectonic processes similar to the Earth's.

M Nafi Toksöz

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Hai Zhang

Seasonal variation of the surface cross-shelf exchange in the northern South China Sea: a Lagrangian perspective

Previous studies on cross-shelf exchange, predominantly adopted an Eulerian perspective, struggled to identify water sources and pathways. Using a high-resolution regional ocean modeling system (ROMS) and Lagrangian particle tracking, this study systematically investigates the seasonal variation and dynamics of surface cross-shelf exchange in the northern South China Sea (NSCS) from a Lagrangian perspective. Based on daily released 30-day drifter trajectories we identify the key pathways, source regions for surface cross-shelf exchange, revealing pronounced seasonal variability. Results show the surface cross-shelf exchange generally following monsoon-driven Ekman transport. However, local dynamics, especially topographic modulation, can outweigh the expected Ekman-driven transport, producing surface exchange patterns opposite to that predicted from the prevailing winds. Topographic effects vary across different scales. In the coastal waters of western Guangdong during winter, despite downwelling-favorable winds, the modulation of alongshore currents by island topography induces an eastward pressure gradient. This gradient, through geostrophic balance, drives offshore flows opposite to wind-driven Ekman onshore transport. Furthermore, the eastern widened shelf exhibits a distinct seasonal variation of cross-shelf exchange, with strong offshore transport (opposite to the direction of Ekman transport) in winter and exceedingly weak exchange in summer. Analysis of the underlying mechanisms reveals that this winter offshore transport is primarily attributed to geostrophic flows driven by surface pressure gradient that is jointly modulated by the Kuroshio intrusion and local widened shelf topography, and enhanced by cumulative submesoscale processes. In summer, a persistent and strong along-isobath jet acts as a dynamic barrier, effectively suppressing the exchange. These findings highlight the important role of topography in regulating surface material transport, and have important implications for predicting the advection and dispersion of plankton or oil spills over the continental shelf influenced by monsoon.

Hao, Dongliang

Experimental Studies of Open Rotor Installation Effects

Open rotor propulsion technologies offer an opportunity for reducing fuel burn. At the same time, great challenges arise from this radically different engine architecture in terms of aircraft system integration. The aeroacoustic effects of integration are one of those key challenges and, furthermore, can offer the possibility for total aircraft system noise reduction. The inter-relationship of the aerodynamic and aeroacoustic system integration effects is particularly important to enable future application. An extensive model scale test campaign was conducted to investigate a broad range of these open rotor installation effects for both a conventional and an unconventional airframe. The experimental campaign was conducted in the Boeing Low Speed Aeroacoustic Facility with specially designed modifications for efficient positioning of the open rotor rig relative to the airframe. The airframe was traversed remotely relative to the fixed open rotor rig providing for the investigation of a large number of installation positions. Eight positions around the main wing of the conventional airframe and eleven positions above the hybrid wing body airframe were documented. Other parameters investigated in the test campaign included forward flight simulation Mach number, angle of attack, rotor speed, and airframe control surface deployment. In addition, both airframes were modified for alternate configurations. The conventional airframe was configured with both a T and a U tail while the unconventional hybrid wing body airframe was configured with different vertical control surfaces. An extensive instrumentation package was deployed. The acoustics were documented with a fixed array of far field microphones, a traversing array of in-flow microphones, an out-of-flow traversing phased array, and unsteady pressure transducers mounted flush on the surface of the airframes. In addition, mean flow surveys were measured with an articulating arm traversing system. The flow field surveys were particularly valuable in documenting flow distortion effects for the various installation positions and those created by angle of attack. This presentation will report the key results obtained for open rotor installation effects and discuss future prospects with the perspective of these results.

Michael J Czech

Asymmetric pathways for lithium extraction and recovery based on the two-phase equilibrium of layered oxides

Electrochemical intercalation offers a promising platform for Li + extraction. However, only limited types of electrode materials have been investigated. The challenge to broaden and tailor materials for electrochemical intercalation-based Li + extraction lies in the lack of understanding of material’s response upon co-intercalation of multiple ions, therefore, paired process design to enable reversible Li + extraction and recovery. Here, we showcase the design of asymmetric ion pathways for Li + extraction and recovery for host material with complex Li + and Na + interaction using layered cobalt oxide as a model material. The two-phase equilibrium of Na 0.48 CoO 2 and Li 0.94 CoO 2 governs Li + selectivity when a high depth of intercalation is achieved (low vacancy level). We show that the relative rate between ion exchange and intercalation is critical to determine the ion pathways. The relationship can be quantitatively compared using the average pseudo ion exchange rate (C pseudoIX ) and the intercalation rate (C inter ). The ion pathways at the three regimes with C pseudoIX > C inter , C pseudoIX ~ C inter , and C pseudoIX < C inter are constructed. By selecting the optimized ion pathway and particle size, we demonstrate 9.7×10 4 Li + selectivity with 99% purity Li + recovery from an initial 1:1000 Li: Na molar ratio solution using 115 mAh/g specific capacity.

electrochemistry

X-Ray Polarization of Z-Type Neutron Star Low-Mass X-Ray Binaries II. Spectropolarimetric Analysis

IXPE has provided for the first time detailed energy- and time-resolved X-ray polarimetry of Z-type neutron star low-mass X-ray binaries (NS-LMXBs) as they move along their color-color diagrams (CCDs). These sources can reach the highest polarization observed for NS-LXMBs in the 2–8 keV range when they move along the horizontal branch. In a previous paper, we characterized the spectral state of a sample of Z-sources using the CCD and estimated the polarization with model-independent analysis. Here, we present detailed spectropolarimetric analysis for each source on each branch using data from IXPE, NICER, and NuSTAR. The continuum X-ray emission of all the sources is well described with a combination of thermal accretion disk emission plus a harder Comptonized component. In addition, reflection features, in particular the relativistically broadened Fe Kα line, are observed for our sources, except GX 5–1. For most of the sources and branches, the main contribution to the X-ray emission and polarization is due to Comptonization: moving from the horizontal branch (HB) to the normal branch (NB), the polarization degree (PD) in the 2–8 keV band varies from about 6% to 3–4%, while the PD is loosely constrained in the flaring branch (FB), due to the shorter exposures. These PD values are significantly higher than theoretical expectations for typical spreading or boundary layer configurations. The polarization of the disk is generally lower (below 3%) but still higher than predictions for an electron scattering-dominated, plane-parallel atmosphere above the disk observed at the corresponding inclination. Moreover, the polarization angle (PA) of the disk seems to be significantly misaligned and not perpendicular to that of Comptonization. We find no correlation between the polarization signal and the inclination, nor with the contribution of reflected photons throughout the Z-track.

accretion

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA

Conceptual design study of neutron detectors for safeguards measurement of an irradiated pebble

Nuclear material control and accounting (MC&A) of pebble-bed reactors (PBRs) is challenging because a PBR utilizes hundreds of thousands of identical, unmarked pebbles that are continuously recirculated through the core. To develop tools that enable the implementation of international safeguards, especially in the context of MC&A of spent pebbles, we designed and simulated three neutron detection concepts to determine fissile content in individual pebbles: a differential die-away (DDA) detector, a californium interrogation prompt neutron (CIPN) detector, and a passive neutron albedo reactivity (PNAR) detector using Monte Carlo calculations. Burnup calculations were performed on the spent pebbles from the PBMR-400 classic PBR. The varying neutron and gamma source terms, and isotopic compositions in the spent pebbles calculated at various burnup levels were used in the neutron detector models. DDA was found to be sensitive to the number of passes a pebble has had through the core and to the fissile content contained in a spent pebble. Optimization in the DDA design further increased the neutron count rates and thus reduced counting uncertainty. Meanwhile, passive neutron counting using the same detector body could distinguish pebbles with different numbers of passes, but its response was dominated by neutron-emitting actinides and was not sensitive to fissile content. On the other hand, the PNAR technique was not viable for a single pebble but performed reasonably for a 27-pebble array, which suggested potential use for verification measurements of containers filled with 27 or more spent pebbles.

CIPN

A Panchromatic View of the Restless SN 2009ip Reveals the Explosive Ejection of a Massive Star Envelope

The double explosion of SN 2009ip in 2012 raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN 2009ip during its remarkable rebrightenings. High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the Very Large Array, Swift, Fermi, Hubble Space Telescope, and XMM) constrain SN 2009ip to be a low energy (E (is) approximately 10(exp 50) erg for an ejecta mass approximately 0.5 solar mass) and asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at approximately 5 × 10( exp 14) cm with M approximately 0.1 solar mass, ejected by the precursor outburst approximately 40 days before the major explosion. We interpret the NIR excess of emission as signature of material located further out, the origin of which has to be connected with documented mass-loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic-ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, which later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the massive luminous progenitor star survived. Irrespective of whether the explosion was terminal, SN 2009ip brought to light the existence of new channels for sustained episodic mass loss, the physical origin of which has yet to be identified.

view

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer

Analysis of Impact Induced Damage and its Effect on Structural Integrity of Space Flight Composite Overwrapped Pressure Vessels

The objective of this research work has been to provide analytical background and support to the ongoing experimental program at NASA, White Sands Test Facility, involving testing composite overwrapped pressure vessels (COPV) for impact damage and cyclic pressurization. Preliminary theoretical basis, including the governing equations for a shallow shell subjected to internal pressure, has been established. Effects of the Griffith type cracks on the structural integrity of the cylindrical vessel were evaluated by methods of Fracture Mechanics. The results indicate that the effective mass of the pressure vessel is an important factor influencing the response to impact events. We also have found that the material properties of the target, contained in the constitutive equations of the composite attached to the Aluminum liner, dominate the impact event in the low velocity range, the material properties become less important, while the target mass distribution and the impactor mass become more significant as the velocity of the impactor increases. Therefore, at high-velocity impact it is not only the kinetic energy of the impactor but also its mass which has a significant effect on the dynamics of the event, and consequently on the induced damage. This work also suggests a methodology for an assessment of the rate of loading effects on the degradation of the material toughness associated with a high-velocity impact where the rate effects become significant. To model the rate dependence of the material response a viscoelastic-plastic constitutive equations were assumed, and on this basis predictions are made regarding the rate dependent material resistance curve. Other dynamic phenomena associated with the impact event have been treated in the framework of the Computational Mechanics using the courtesy of Prof. P. Guebelle and his graduate student at University of Illinois at Urbana-Champaign who have an access to a super-fast computer located on their campus. Finally, the guidelines for a follow-up research program are provided in the body of this report. They address three major areas: theoretical research, numerical studies, and further experimental work.

Michael P Wnuk

Optical transparency of LiF (100) shock compressed to ∼360 GPa

High-purity [100] lithium fluoride (LiF) is the most widely used optical window in dynamic compression experiments due to its wide bandgap and well-characterized mechanical response. Recent plate-impact experiments established the [100] LiF Hugoniot to ∼230 GPa and demonstrated shock-induced melting onset at 182 GPa with complete melting by 195 GPa; theoretical models predict LiF optical transparency to nearly 900 GPa. To experimentally examine the optical transparency of [100] LiF at higher pressures and in the liquid state, laser-driven shock experiments were performed at peak stresses ranging from 223 to 363 GPa. Optical response was examined by measuring the particle velocity histories at the Kapton/LiF interface using laser interferometry at 532 and 1550 nm wavelengths; in-material particle velocities were obtained using established refractive-index corrections. Continuous photonic Doppler velocimetry fringes were observed across the entire stress range, demonstrating that LiF remains transparent to 1550 nm light throughout the multi-megabar regime investigated. At 532 nm, fringe visibility depended on the reflector coating: aluminum mirrors provided signals to ∼235 GPa, while gold mirrors extended this limit to 270.5 GPa, indicating that the shorter-wavelength response is likely sensitive to experimental configuration rather than to the loss of LiF transparency. Continued optical transparency to at least 360 GPa indicates that shock-melted LiF does not display bandgap closure over the stress range explored. Furthermore, these results provide direct experimental constraints on the high-pressure optical response and establish LiF (100) as a robust optical window material for laser-driven dynamic compression experiments approaching 400 GPa.

Renganathan, P. [Argonne National Laboratory (ANL)

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

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato