Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source term calculation”

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.

807 records · Page 45

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗

Potential State Regulatory Pathways to Facilitate Low-Carbon Fuels

States and the federal government are increasingly engaged in the challenges around decarbonizing the electric grid. In particular, regulators, consumers, stakeholders, and utilities recognize the need to carefully consider the role natural gas will play in a decarbonized future. A variety of technology and policy options to reduce greenhouse gas emissions associated with natural gas use are available, including energy efficiency programs, demand reduction tools, strategic electrification, and strategies to reduce emissions from natural gas production, transportation, and consumption. Low-carbon fuels – mainly renewable natural gas (RNG) and clean hydrogen – are being considered an important component of decarbonization goals. RNG and hydrogen may be able to meaningfully reduce emissions from processes independent of geologic natural gas, displacing emissions of methane, a powerful greenhouse gas. Although RNG and hydrogen are not cost-competitive today with geologic natural gas and are smaller in scale and potential than other decarbonization options, they can be explored as potential critical tools to decarbonize sectors that are difficult to electrify or shift off of natural gas entirely, such as air travel, industrial processes, maritime transport, long-distance trucking, space heating on cold days, and railroads (Nadel, 2022). The role of this report is to provide informational context for state utility regulators to understand the impacts of and challenges associated with broader integration of low-carbon fuels, followed by examples of state regulatory actions taken to date to facilitate the development of low-carbon fuels. Setting clear guidance to calculate the environmental benefits of low-carbon fuels and continuing federal and state investments in research and development to reduce costs relative to fossil fuels will be important steps to take to signal the desire to grow the market for these fuels. State public utility commissions may play a key role in setting regulatory frameworks for low-carbon fuels and ensuring that ratepayer funds, if utilized, are done so to further the public interest. This report is intended to summarize decisions that states have made to date on low-carbon fuels. In the spirit of understanding the current market and sharing information, this report provides success stories, and lessons learned across states as regulators implement varying strategies to achieve decarbonization objectives while maintaining their focus on affordability, safety, and reliability of the energy system. The report begins with an introduction of the role of natural gas in the U.S. economy (Section I) and background information on natural gas use, decarbonization, and low-carbon fuels (Section II). Next, the report describes the current market by discussing the scale of current production, emissions intensity, resource potential, and costs of low-carbon fuels compared to geologic natural gas (Section III). Following these sections, the report describes four strategies states have employed to facilitate low-carbon fuels: opening exploratory dockets, approving voluntary tariffs for customers, approving interconnection tariffs for producers, and considering portfolio-wide procurement targets (Section IV). This section lists states that have taken actions in each category, citing utility filings, commission decisions, stakeholder comments, and other relevant sources. Finally, the report concludes with suggested questions regulators may wish to consider regarding low-carbon fuels, in the interest of preparing to make decisions in the future (Section V). These questions include: Are there existing regulatory or technical barriers to voluntary purchases of low-carbon fuels? Can customers work with utilities to procure low-carbon fuels; are producers able to interconnect projects without significant barriers to entry? Should the infrastructure and/or commodity costs of low-carbon fuels be socialized among all ratepayers, or borne solely by the large commercial and industrial (C&I) customers currently driving the market? Should regulated natural gas and/or electric utilities own and operate low-carbon fuel production? How should regulators consider the unique decarbonization potential of low-carbon fuels, particularly for hard-to-abate sectors, in decision-making? Is additional direction or clarity from state policymakers needed? What no-regrets approaches can help facilitate both near-term RNG development and long-term development of hydrogen and other zero-carbon fuels? We collectively wish to express our gratitude to the U.S. Department of Energy, Office of Fossil Energy and Carbon Management, for supporting this report and other technical assistance resources for state regulators on natural gas topics. State regulators operate under a variety of policy environments, and states have vastly different types of energy resources, infrastructure, and customers. While there is no optimal regulatory, policy, or technological solution that will be successful in every state, state regulators can benefit by exchanging lessons learned with their peers across the country. We look forward to continued engagement with our fellow commissioners, commission staff, NARUC, the U.S. Department of Energy, and other stakeholders to develop sound regulation in the public interest.

03 NATURAL GAS↗

The Nature and Evolutionary History of GRO J1744-28

GRO J1744-28 is the first known X-ray source to display bursts, periodic pulsations, and quasi-periodic oscillations. This source may thus provide crucial clues that will lead to an understanding of the differences in the nature of the X-ray variability from various accreting neutron stars. The orbital period is 11.8 days, and the measured mass function of 1.31 x 10(exp -4) solar mass is one of the smallest among all known binaries. If we assume that the donor star is a low-mass giant transferring matter through the inner Lagrange point, then we can show that its mass is lower than approximately 0.7 solar mass and probably closer to 0.25 solar mass. Higher mass, but unevolved, donor stars are shown to be implausible. We also demonstrate that the current He core mass of the donor star lies in the range of 0.20-0.25 solar mass. Thus, this system is most likely in the final stages of losing its hydrogen-rich envelope, with only a small amount of mass remaining in the envelope. If this picture is correct, then GRO J1744-28 may well represent the closest observational link that we have between the low-mass X-ray binaries and recycled binary pulsars in wide orbits. We have carried out a series of binary evolution calculations and explored, both systematically and via a novel Monte Carlo approach, the range of initial system parameters and input physics that can lead to the binary parameters of the present-day GRO J1744-28 system. The input parameters include both the initial total mass and the core mass of the donor star, the neutron-star mass, the strength of the magnetic braking, the mass-capture fraction, and the specifics of the core mass/radius relation for giants. Through these evolution calculations, we compute probability distributions for the current binary system parameters (i.e., the total mass, core mass, radius, luminosity, and K-band magnitude of the donor star, the neutron star mass, the orbital inclination angle, and the semimajor axis of the binary). Our calculations yield the following values for the GRO J1744-28 system parameters (with 95% confidence limits in parentheses): donor star mass: 0.24 solar mass (0.2-0.7 solar mass); He core mass of the donor star: 0.22 solar mass (0.20-0.25 solar mass); neutron-star mass: 1.7 solar mass (1.39-1.96 solar mass); orbital inclination angle: 18deg (7deg-22deg); semi- major axis: 64 lt-s (60-67 lt-s); radius of the donor star: 6.2 solar radius(6-9 solar radius); luminosity of donor star: 23 solar luminosity (15-49 solar luminosity), and long-term mass transfer rate at the current epoch: 5 x 10(exp -10)solar mass/yr (2 x 10(exp -10) to 5 x 10(exp -9)solar mass/yr). We deduce that the magnetic field of the underlying neutron star lies in the range of approximately 1.8 x 10(exp 11)G to approximately 7 x 10(exp 11)G, with a most probable value of 2.7 x 10(exp 11)G. This is evidently sufficiently strong to funnel the accretion flow onto the magnetic polar caps and suppress the thermonuclear flashes that would otherwise give rise to the type 1 X-ray bursts observed in most X-ray bursters. We present a simple paradigm for magnetic accreting neutron stars where X-ray pulsars, GRO J1744-28, the Rapid Burster, and the type 1 X-ray bursters may form a continuum of possible behaviors among accreting neutron stars, with the strength of the neutron-star magnetic field serving as a crucial parameter that determines the mode of X-ray variability from a given object.

Rappaport, S.↗

DIF3D-VARIANT 12.0: Updates and New Features

The DIF3D code has been a workhorse of fast reactor analysis work at Argonne National Laboratory for over 40 years. In 1995, a transport option called VARIANT was added to DIF3D to improve the flux solutions for fast reactor problems which we term DIF3D-VARIANT today. DIF3D-VARIANT performs nodal neutron transport calculations using P N or SP N theory in Cartesian and hexagonal two- and three-dimensional geometries. The limited computing capabilities of the time restricted DIF3D-VARIANT to use at most a 6 th order spatial approximation combined with a P3 flux approximation and P1 scattering kernel for a 33 group structure on most studied reactor problems. Computer capabilities have increased steadily since 1995 and today much larger space-angle-energy approximations are possible. This manuscript serves as an update to the theory section of the original DIF3D-VARIANT manual and details more than twenty years of changes made to DIF3D to make version 12 which was released on November 1 st , 2024. The primary focus of the initial work was to extend the space-angle approximations available in DIF3D-VARIANT such that the error due to transport approximations could be better understood. This work was started and completed in 2002 and marked the official version 10. Unfortunately, those higher order approximations could not be used at that time due to the memory constraints of the BPOINTER part of DIF3D (limited to 2 GB). In version 11, completed in 2012, BPOINTER was circumvented in DIF3D-VARIANT for the largest arrays by introducing a Fortran 90 module called LMA (Large Memory Array). This seamlessly replaces all of the functionality of the BPOINTER concept, but it allows 64 bit addressing for every array such that they can be larger than 2 GB. It is now common for DIF3D-VARIANT jobs to consume 50 GB of memory on modern workstations when using high order space-angle approximations and a large number of groups. Many improvements were made to version 11 from 2012 to 2022 when work to create version 12 started. For version 12, several parts of DIF3D were updated to improve performance and thread parallelism was introduced to further reduce the runtime. Numerous minor bugs were discovered in DIF3D-VARIANT as part of the process of creating the perturbation and sensitivity code PERSENT. All of these algorithmic problems were identified in the transition from version 10 to version 11 which prevented DIF3D-VARIANT from running efficiently and reliably. Firstly, the coarse mesh rebalance scheme would routinely diverge and a study detailed in this report demonstrates how it was also typically not effective. This is not a failure of the coarse mesh rebalance methodology, but a failure of its implementation in DIF3D-VARIANT for hexagonal geometries. The fission source extrapolation algorithm was also found to be unreliable on larger group structure problems, leading to divergence in some cases and a negligible improvement in performance overall. Finally, the “Omega” acceleration applied to the partial current solver routine of DIF3D-VARIANT was found to cause DIF3D-VARIANT to converge to the wrong answer. To resolve these issues, both the coarse mesh rebalance and fission source extrapolation were permanently disabled in version 11. The Tchebychev acceleration was put in as a temporary reliable alternative but it is generally inferior to coarse mesh rebalance or coarse mesh finite difference. For the Omega acceleration, the factor was restricted to guarantee that it would not cause follow-on errors in PERSENT. Due to limited funding to support maintenance and development of DIF3D in the last 10 years, no effort was spent since to resolve the outer iteration acceleration. Except for the threading work, all of the changes discussed in this manuscript refer to changes made between version 10 and version 11. Performance comparisons are done to demonstrate the improvements from version 9 to version 12. As will be demonstrated, the updated versi

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Major improvements on the longwave radiative interactions between surface and clouds in the Polar Regions in atmospheric global circulation model (GCM). Final Project Report

1. Summary of Work and Accomplishments in Year 1 (01/15/2015-01/14/2016). We developed new ice cloud parameterizations based on state-of-the-art MODIS Collection 6 (MC6) ice particle model. Features of the model include the consideration of surface roughness of ice crystals and spectral consistency in cloud property retrievals. Using more than 14,000 particle size distributions (PSDs) from in-situ measurements and the latest ice particle single-scattering property library, our parameterizations provide band-averaged bulk-scattering properties. According to statistics of fittings, the parameterization schemes represent the band-averaged bulk-scattering properties very well. Compared to the ice particle model developed by Fu (1998), flux simulations with the MC6 particle model result in larger upward flux at the top of the atmosphere (TOA) and smaller at the surface for typical cirrus clouds. In the assessment of the effect of longwave light scattering, our results also support the conclusions of previous studies. The error in flux calculations reaches about 10 W/m 2 when a cirrus cloud is thin and its particle size is small. 2. Summary of Work and Accomplishments in Year 2 (01/15/2016-01/14/2017) Clouds are a major modulator of the global radiation budgets. However, representations of clouds in general circulation models (GCMs) neglect many important processes, particularly those related to ice clouds. One major uncertainty source is the radiative transfer scheme in GCMs. To reduce the computational burden in most GCMs, absorption is the only radiative transfer process considered in longwave spectral bands. Recently, a state-of-the-art MC6 ice model has been available, which provides spectral consistency in ice cloud retrievals. With an advanced light scattering library for ice particles, we performed sensitivity studies of the potential impacts on the climate including light scattering for MC6 model and other 11 ice cloud models. Specifically, we modified RRTMG_LW (the GCM version of Longwave Rapid Radiative Transfer Model), a radiative transfer module widely used in GCMs and numerical prediction models, to identify flux simulation differences with and without considering light scattering processes. The simulations with light scattering come from a rigorous radiative solver, DISORT (Discrete Ordinates Radiative Transfer Program for a Multi-Layered Plane-Parallel Medium), are used as benchmarks. The results show that the weighted annual mean biases of RRTMG_LW for the upward flux at the top of the atmosphere, the downward flux at the surface, and the net flux into the atmosphere are about 0.8 ± 0.3, -0.1 ± 0.05, and -0.7 ± 0.3 W m-2, respectively. According to the different the ice crystal shapes, the weighted annual mean bias for heating rate is about -0.006 ± 0.02 to 0.04 K/day. 3. Summary of Work and Accomplishments in Year 3 (01/15/2017-01/14/2018) Since absorption dominates optical properties of clouds in the longwave (LW) spectrum, most of general circulation models (GCMs) only take absorption properties of clouds into account in the radiative transfer modules in order to reduce computational costs. To quantify the biases of excluding LW scattering when clouds exist, we simulated fluxes and heating rates by using satellite data from Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO), CloudSat, Clouds and the Earth’s Radiant Energy System (CERES) and Moderate Resolution Imaging Spectrometer (MODIS) merged products (CCCM) in 2010. In annual global average, neglecting LW scattering causes overestimation by 2.6 W/m 2 for the top-of-atmosphere (TOA) upward fluxes and underestimation by 1.2 W/m 2 for surface downward fluxes. In particular, regional extreme biases are approximately 12 W/m 2 for TOA upward fluxes and are approximately -3.6 W/m 2 for surface downward fluxes, or are approximately 5% of the global averaged outgoing longwave radiation (OLR) and 1% of the global averaged surface downward fluxes, respectively. In terms of heating rates, on average, LW scattering heats the whole atmosphere column by 0.0045 K/day, heats the cloud layers by 0.0420 K/day, heats the surface by 0.028 K/day, and cools the tropopause by 0.018 K/day.

58 GEOSCIENCES↗

Laser Sounder for Measuring Atmospheric CO2 Concentrations: Progress Toward Ascends

The next generation of space-based, active remote sensing instruments for measurement of tropospheric CO2 promises a capability to quantify global carbon sources and sinks at regional scales. Active (laser) methods will extend CO2 measurement coverage in time, space, and perhaps precision such that the underlying mechanisms for carbon exchange at the surface can be understood with .sufficient detail to confidently project the future of carbon-climate interaction and the influence of remediative policy actions. The recent Decadal Survey for Earth Science by the US National Research Council has recommended such a mission called the Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) for launch in 2013-2016. We have been developing a laser technique for measurement of tropospheric CO2 for a number of years. Our immediate goal is to develop and demonstrate the method and instrument technology that will permit measurements of the CO2 column abundance over a horizontal path and from aircraft at the few-ppmv level. Our longer-term goal is to demonstrate the required capabilities of the technique, develop a space mission approach, and design the instrument for an ASCENDS-type mission. Our approach is to use a dual channel laser absorption spectrometer (i.e., differential absorption in altimeter mode), which continuously measures from a near-polar circular orbit. We use several co-aligned tunable fiber laser transmitters allowing simultaneous measurement of the absorption from a CO2 line in the 1570 nm band, O2 extinction in the oxygen A-band (near 765 nm), and aerosol backscatter in the same measurement path. We measure the energy of the laser echoes at nadir reflected from land and water surfaces, day and night. The lasers have spectral widths much narrower than the gas absorption lines and are turned on and off the selected CO2 and O2 lines at kHz rates. The gas extinction and column densities for the CO2 and O2 gases are estimated from the ratio of the on and off-line singnals via the DIAL technique. We used pulsed laser signals, photon counting detectors, and time gating to isolate the laser returns from the surface, and to reject photons scattered from thin clouds and aerosols. High signal-to-noise ratios are required and the CO2 estimates can be sensitive to small drifts or other errors in the instrument, so the absorption estimates need to be quite stable for hours. We have constructed a breadboard version of the CO2 sensor that uses a low power fiber laser and a 20 cm diameter telescope. We have used it to make measurements of CO2 absorption in the laboratory and over 200-m to 2-km long open horizontal paths. These have been done in several sessions extending over multiple days, which allows us to assess the measurement stability and to compare absorption variations to readings from an external in situ CO2 sensor. We have also calculated characteristics of the technique for space including its expected measurement performance for different modulation types, and have performed an initial space mission accommodation study. We sill describe these results in the presentation.

Abshire, J. B.↗

Collaboration with NSTX in Calculations of Radiofrequency and Neutral Beam Heating and Current Drive

The overall objective of this project was to advance the understanding of heating and current drive necessary for toroidal magnetic confinement fusion energy reactors, and to provide accurate models for these sources to aid in understanding magneto-hydrodynamic, transport and boundary physics issues. Increased understanding is provided by computer simulation using the underlying physics laws, and will lead to greater reliability and increased efficiency of a power producing fusion reactor. The proposed work entailed applications of existing CompX Fokker-Planck, ray tracing, full-wave microwave, and particle diffusion calculator codes to the National Spherical Torus Experiment at Princeton Plasma Physics Laboratory. Special features of NSTX include that large ion orbits are emphasized in its geometry and thus offer a major opportunity within fusion research for focused exploration of their ubiquitous effects, including comparison with simulation codes. The CompX CQL3D Fokker-Planck code was be further adapted to the finite-orbit-width physics. Comparison with experiment ensured that the computer codes have sufficient accuracy to accurately model the experiments. Nonthermal distributions of energetic ions are produced in experiments-- and modeled with the codes-- exhibiting the signatures of auxiliary heating, transport, and radiation processes within the plasma. Using physics-based modeling, in conjunction with experimental data, nonthermal distributions can provide keys to validation and quantification of auxiliary heating physics. Nonthermal distributions also hold particular information on radial transport processes, not obscured by the smoothing Maxwellization effects of Coulomb collisions. Generally, radial transport is not well understood but is a key to efficient fusion energy production. By comparison with Fokker-Planck modeling, the velocity dependence of radial transport coefficients becomes experimentally accessible, particularly in cases of well-modeled sources such as neutral beams, fast wave power, or the toroidal electric field. It is clear from modeling of fast-wave interactions with ions that finite-orbit-width effects play a significant role. For example, a different set of ions is computed to intersect a given cyclotron resonance depending on whether or not radial orbit drifts are taken into account. Also, radial transport, associated bootstrap current, and alpha channeling depend on radial drifts. Ion cyclotron heating above the first harmonic interacts preferentially with higher velocity ions, giving large radial drifts. The CQL3D Fokker-Planck code, maintained by the proposers, is a physics based model used extensively within the US fusion community to compare with and interpret experimental observations. At the beginning of this grant, bounce-averages of the rf quasilinear and other operators have been performed neglecting the radial drifts which give finite orbit widths to the particle orbits. The proposed work included a major upgrade to CQL3D to include finite-orbit-width effects be accounted for by performing the bounce-averages of rf and neutral beam fast ions along guiding center orbits. This CQL3D enhancement provided microwave and neutral beam contributions to neoclassical radial transport coefficients for nonthermal distributions, and substantially increased the comprehensiveness of the ion model. Additional work validated Fokker-Planck and ray tracing modeling of neutral beam and high harmonic fast wave absorption and current drive against the NSTX experiment, as had been very successful for electron cyclotron current drive in the DIII-D tokamak at General Atomics and lower hybrid waves on C-Mod at MIT. The validation work was in conjunction with staff and graduate students at Princeton Plasma Physics Laboratory. This DOE grant supported fusion energy research, a potential long-term solution to the world's energy needs. The research grant under consideration, DE-SC0006614, provided research support at the 0.3 Full-Time-Equivalent level for application of computer models to aid in understanding and projecting efficacy of heating and current drive sources in the National Spherical Torus Experiment, a tokamak variant, at the Princeton Plasma Physics Laboratory. The NSTX experiment explores the physics of very tight aspect ratio, almost spherical tokamaks, aiming at producing steady-state fusion plasmas. The current drive is an integral part of the steady-state concept, maintaining the magnetic geometry in the steady-state tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CO2 Capture Strategies via Mineralization with Industrial Waste Brines

Large coal-fired power plants (>500 MW) account for 30% of global CO2 emissions, and long-term management of this CO2 to is urgently needed mitigate global temperature increases. Sequestration of CO2 within stable mineral carbonates (e.g., CaCO3) represents an attractive emission reduction strategy because it offers a leakage-free alternative to geological storage of CO2 in an environmentally friendly form. We have previously described a mineralization process in which divalent cations are sourced from various waste streams (e.g., produced water and brackish water) and alkalinity is induced via regenerable ion-exchange materials (Bustillos et. al. Frontiers in Energy Research. 2020, 8, 352). In our process, aqueous carbonate-bearing streams with pH > 8 are produced by contacting fresh water and carbon dioxide with various ion-exchange materials (e.g., Na form zeolites or ion exchange resins). These streams are mixed with produced water containing varying concentrations (~0.01 – 1.0 M) of Ca2+ leading to the precipitation of solid calcium carbonate (PCC). This process has the advantages of using regenerable solids in a simple and continuous process to increase the pH of water by ion exchange instead of relying on the consumption of costly and unsustainable sources of alkalinity (e.g., sodium hydroxide). While once-through column experiments showed the above benefits, the same were yet to established in a steady-state process with recycle streams. In this work, we set up a process simulation to quantify the energy requirements and CO2 emissions associated with the process and seek optimal produced water compositions and CO2 concentrations (5 – 20 vol%). The process simulation was set up in ASPEN Plus using eRNTL as the thermodynamic property method and sequential modular strategy. Ion exchange alkaline solution was simulated using sodium hydroxide and validated against the experimental data obtained from once-through kinetic experiments. Nanofiltration and reverse osmosis membrane steps were also implemented for the separation of divalent cations and production of fresh water and a regeneration stream following mineralization. Sensitivity analysis was carried out using a range of produced water compositions (0.01 – 1.0 M Ca2+, 0.001 – 0.15 M Mg2+, 0.5 – 3.5 M Na+ and 0.0004 – 0.002 M Fe2+) according to the United States Geological Survey (USGS) database. Calcium carbonate yields increased with increasing CO2 concentrations and were maximized using produced water compositions with larger Ca2+ concentrations. Maximum calcium carbonate yields produced at 5 vol%, 12 vol% and 20 vol% CO2 were 2.3 mmol/L, 5.5 mmol/L, and 9.3 mmol/L, respectively, with the formation of brucite (a magnesium hydroxide phase, Mg(OH)2) and goethite (an iron hydroxide phase, FeOOH) as the primary contaminant phases (99% calcite, 0.6% brucite, 0.4% goethite), which agree with phases detected by XRD experimentally. These results indicate high purity calcium carbonate can be precipitated using industrial waste streams. Consequentially, energy consumption and net CO2 emissions were minimized where precipitated calcium carbonate was maximized for all produced water compositions and CO2 concentrations. Minimum energy consumptions were 0.21 kWh/ton CO2 processed, with 98% of the energy input required coming from the membrane filtration steps. Produced water compositions with large Na+ concentrations (> 0.5 M) were effective at reducing energy consumptions due to faster regeneration time of ion exchange materials. Additionally, calculated net CO2 emissions were negative for the process and ranged from -0.02 kg/ton CO2 to -0.15 kg/ton CO2 processed, indicating a low emission process. We will also present techno-economic assessment showing the economic benefits of the current process as an alternative to the addition of stoichiometric bases to induce alkalinity for the precipitation of CaCO3.

Simonetti, Dante↗

R and DE Robotic Sensor Intern Team: Remote 3D Mapping Overview

The Department of Energy (DOE) nuclear complex includes many hazardous and challenging environments that must be inspected remotely and routinely for structural integrity evaluation and long-term planning purposes. The H-Canyon Air Exhaust (HCAEX) Tunnel biannual inspection is an example of such an environment. Presently inspections are performed by teleoperated robotic crawlers provide low resolution video feedback to the structural engineers. Advanced sensors such as Lidar and panospheric hi-res cameras have the potential to provide additional valuable data in future inspections. UT-Austin Lidar work in support of an H-Canyon deployment began with the DOE EM sponsored H-Canyon Integrated Research Project. The R and DE Robotic Sensor intern team includes two of the UT-Austin Masters students working on Lidar development to facilitate the pulling of university developed solutions to the field. The goal of this project is to develop a suite of sensors that could be used to gather advanced data in future HCAEX tunnel inspections. Of particular interest is the ability to map the tunnel walls at an accuracy that will allow for the comparison of two scans performed at different time periods to estimate the concrete erosion rate. Additionally, the development of the 3D point cloud in the GPS denied environment is the first step in being able to collect and overlay other 3D data such as radiation mapping. Project specific goals include: - Design and develop 3D mapping sensor suite. Merge Lidar point cloud with high res camera images Incorporate Inertial Measurement Unit (IMU) to assist with positioning. - Develop sensor software and user GUI Sensor data collection and integration software User GUI for reach back system. - Test solution on crawler in tunnel test bed Use crawler base similar to H-Canyon planned crawler Design and build test bed simulating tunnel conditions. To facilitate rapid sensor integration, software coding and user GUI development, the Robot Operating System (ROS) was employed as the development platform. ROS is a set of open source software libraries and tools that assist with the building of robot applications. ROS includes drivers for many sensors, state-of-the-art algorithms, and powerful developer tools. Sensors selected for the summer intern solution included assessment as to whether they possessed ROS drivers and associated sample code. Additionally, the ROS 3D Robot Visualizer (rviz) tool has been employed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Measurement of the U 235 ( n , f ) prompt fission neutron spectrum from 10 keV to 10 MeV induced by neutrons of energy from 1 MeV to 20 MeV

The characterization of fission-driven nuclear systems primarily relies on calculations of neutron-induced chain reactions, and these calculations require evaluated nuclear data as input. Calculation accuracy heavily depends on input nuclear data evaluation accuracy, and thus high precision on the experimental input to the nuclear data evaluation is essential for fundamental quantities like the energy spectrum of neutrons emitted from neutron-induced fission (i.e., the prompt fission neutron spectrum, PFNS). Despite decades of measurement efforts, prior to the measurements described in this work there were only three literature data sets for the 235 U(n,f) PFNS at incident neutron energies above 1.0 MeV considered reliable for inclusion in nuclear data evaluations and no reliable data sets above 3.0 MeV incident neutron energy. In this work we report on new measurements of the 235 U(n,f) PFNS spanning a grid of 1.0–20.0 MeV in incident neutron energy and 0.01–10.0 MeV in outgoing (PFNS) neutron energy. These measurements were carried out at the Weapons Neutron Research facility at the Los Alamos Neutron Science Center and used a multifoil parallel-plate avalanche counter target with both a Li-glass and a liquid scintillator detector array in separate experiments to span the quoted outgoing neutron energy ranges. The PFNS results are shown in terms of the energy spectra themselves as well as the average PFNS energy $(\langle{E}\rangle)$ and ratios of $\langle{E}\rangle$ at forward and backward angles. Here, the results are compared with literature data and selected nuclear data evaluations. Generally, the data agree with the ENDF/B-VIII.0 evaluation below 5.0-MeV incident neutron energy and more closely with the JEFF-3.3 evaluation above 5.0 MeV, though no evaluations considered for comparison in this work agree with the data across all of the incident and outgoing neutron energies shown, especially in regions where the third-chance fission process becomes available. Additionally, we show a ratio of the present PFNS results for 235 U(n, f) with a recent and highly correlated experiment to measure the 239 Pu(n, f) PFNS at the same experimental facility and with nearly identical equipment and analysis procedures. Many observations reported in this work are the first of their kind and represent significant advancements for knowledge of the 235 U(n, f) PFNS.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Self-Powered Autonomous Sensing System for Arctic Ocean using a Frequency-multiplied Cylindrical Triboelectric Nanogenerator

An autonomous sensing system for collecting environmental observations in the Arctic region is critical for the estimation and prediction of climate change. Ocean waves are a great source of energy for these sensing systems but there has been limited research done on small-scale energy harvesting applications in the Arctic Ocean (subsea or surface). The available wave energy in the Arctic Ocean is lower than the typical ocean wave energy due to the low wave frequency, height, and operating months. Although the available wave energy depends on the specific location in the Arctic Ocean, the wave height decreases everywhere during the winter (Jan-March) [1]. Our target location in this work is the Beaufort and Chukchi Seas, which are ice covered for about 195 days per year leaving only a few months (June-November) for wave energy harvesting [2]. During these few months, the average wave frequency is 0.2 Hz and the most common wave frequency is around 0.15 Hz. Many energy harvesting methods are unsuitable for use in the Arctic Ocean because of to the cold temperature, low wave frequency, and low wave height. Triboelectric nanogenerators (TENG) are one of the few energy harvesting methods that excel in these conditions. Zhong et al. [3] designed a stacked pendulum-structured TENG for low-frequency ocean wave energy and reported the device generated a peak power density of 11.2 W/m3 under the low wave frequency of 0.2 Hz. The autonomous sensing system we developed for the Arctic Ocean, the Arctic-TENG, is based on a frequency-multiplied cylindrical TENG (FMC-TENG) which is an optimized TENG configuration for Arctic Ocean conditions due to the high power density under low-frequency wave conditions [4]. Figure 1 shows an FMC-TENG with multiple pairs of free-standing triboelectric-layer mode materials (Aluminum and FEP). The mass and magnet attached to the rotor store gravitational potential energy which is released as kinetic energy when the potential energy overcomes the repulsive magnetic force generated by the opposing magnet attached to the stator. This force unbalance triggers a sudden rotation and swinging motion of the mass which increases the angular velocity of the system and therefore enhances the output power of the TENG system. The main components of the Arctic-TENG are a 3D-printed rotor and stator, electric and dielectric material adhesive tape, and bearings. All these materials were cold-soaked and tested at -40 °C in a chest freezer. Both the 3D-printed parts and adhesives were confirmed to have a minimal effect from the cold temperature. Multiple bearings were tested and the starting torque of each one was compared both at room temperature and -40 °C. The bearing with the lowest starting torque at -40 °C was selected for use in the Arctic-TENG. The Arctic-TENG was tested using an out-of-water motor-driven wave simulator (Figure 2). The wave simulator allows for controlled testing at a wave height of 0.2 m and frequencies between 0.1 Hz to 0.5 Hz. This wave simulator was used for both room temperature and -40 °C testing. The Arctic-TENG generated significantly more power at -40 °C compared to room temperature for each frequency tested. The system stored energy in a supercapacitor via a power management circuit, and the amount of energy stored per day was calculated at different wave frequencies. Based on the conditions at the proposed deployment location (days of non-ice-covered ocean and wave frequency), the amount of stored energy per year from Arctic-TENG was calculated to be enough energy for two transmissions every day. Durability testing was completed on the Arctic-TENG at room temperature to determine the lifetime. The rotor of the Arctic-TENG was attached to a DC motor and spun continuously for several millions of cycles without degradation of the electrical output, demonstrating the feasibility for long-term operation.

Jung, Hyunjun↗

Non–linear bonding trends in maleonitrile-1,2–dithiolate complexes of the transuranium actinides

The trivalent actinides are produced in the nuclear fuel cycle during power production and provide the largest long-term radiation dose in used nuclear fuel. It is ideal for these elements to be removed from used nuclear fuel for disposal and a necessity for fuel recycling. A key challenge to this is the similarity of chemical behavior of the trivalent actinides to the lanthanides that are also present as fission products in used fuel. Thus far, some of the most effective separations of actinides from lanthanides utilise chelating agents containing sulfur moieties such as dithiophosphinates that selectively bind to actinide ions because of a greater bond covalency relative to lanthanide ions. Typically, greater differences between actinide and lanthanide ions are observable the more ligands and chelators bonds have a covalent character. Here, a series of complexes of the trivalent actinides Np(III) through Cf(III) (excluding Bk(III)) with maleonitrile-1,2-dithiolate (mnt 2– ) are synthesized along with their lanthanide counterparts (La(III) – Nd(III), Sm(III) – Gd(III), Dy(III)), in order to characterize the nature of chemical bonds with these metal ions and a polarizable, non-innocent, sulfur-donor ligand. The metal-sulfur bonds in these complexes trend shorter than measured for lanthanides with equivalent ionic radii. However, particularly large deviations are observed in the neptunium and plutonium complexes in both structure and bonding, resulting in a nonlinear bond length trendline for the actinide series. Density Functional Theory (DFT) calculations with Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Order (NBO) analyses indicate that for the neptunium and plutonium complexes, the presence of increased 5f-orbital participation, energy degeneracy of the metal and ligand orbitals, and the structure packing result in shortened M–S bonds. The stabilization of the energy of the 5f-orbitals and the decrease in f-contribution to bonding orbitals in the later actinides results in structural properties more similar to the lanthanide complexes.

07 ISOTOPE AND RADIATION SOURCES↗

Multi-Sourced Collaboration for the Production and Refining of Rare Elements and Critical Metals (Final Technical Report)

The project objective was to develop a feasible and cost-effective method for recovering rare earth elements (REEs) and critical materials (CMs) from coal and coal byproducts, resulting in high-purity individually separated REEs and CMs. The targeted REEs included Y, Pr, Nd, Gd, Dy, and Sm, with a purity of over 99.5%, while the CMs included Co, Mn, Ga, Sr, Li, Ni, Zn, and Ge, with a purity of over 90%. The project aimed to design a prototype facility capable of producing 1-3 tonnes/day of high-purity REO mixes. The work was divided into four designated circuits: 1) REE extraction and concentration, 2) REE separation and purification, 3) RE metal production, and 4) CM production. To achieve these goals, the project involved 11 tasks, including technology reviews, research, process flow diagram development, mass balance estimation, and preliminary technical-economic analysis. The project team included researchers from the University of Kentucky, University of Alabama and Virginia Tech as well as process specialists from Argonne National Laboratory. MP Materials provided technical support regarding rare earth markets and processing while Alliance Coal performed resource assessment. The project included a market analysis for Nd/Pr, Tb, Dy, Gd, Y, Co, Mn, Li, Sr, Ga, Ni, Zn, and Ge. These analyses provided insights into the supply and demand trends as well as historic and future projections of market price relative to purity requirements for these elements. Two coal resources were selected for the project: the West Kentucky No. 13 (Baker) Seam and an undisclosed lignite resource in the Illinois coal basin. The estimated quantities of REEs in these resources were calculated based on production samples and drilling data. It was estimated that there is adequate supply for an operation producing one metric ton daily of higher purity mixed rare earth oxides (MREO) for approximately 20 years at a site located in western Kentucky. In Circuit 1, project data was obtained from a pilot heap leach and REE concentration facility. It was concluded that the existing circuit, which generated a MREO concentrate, two types of CM mixed products, and Li- and Sr-containing waters, would be suitable feed for circuits 2-4. Data from the first-of-its-kind coal coarse refuse heap leach pilot pad played a crucial role in estimating reliable elemental concentrations of the pregnant leaching solution (PLS). The average total REE concentration in the PLS was found to be 28.6 ppm. In Circuit 2, several concepts were explored including a novel process referred to as solvent-assisted chromatography (SAC). This concept involved a novel columnar reactor that incorporated multiple mixer/settlers, thereby enabling the operation of counter-flowing aqueous and organic phases. Unfortunately, due to project time constraints, a complete fundamental modeling analysis could not be completed to fully evaluate the technology. Molten salt electrowinning was considered as an alternative for circuit 3 following circuit 2 purification circuit utilizing the novel SAC process. A mass and energy balance of Nd reduction to metal in a fluoride containing molten salt electrolyte was conducted. Comparisons were made with the current state of Asian molten salt electrorefining, and potential improvements in siphoning rare earth metals (REM) from the reactor were presented. A cost estimate was performed for the production of 1 tonne per day, which yielded a total of $2.29 million for the nine electrowinning (EW) cells required. The selected option for circuits 2 and 3 was a plasma distillation process, which initially separates rare earth elements (REEs) from other elements. This is followed by selective electrowinning in various ionic liquids. The selection was made on the basis of thermodynamic modeling and experimental data previously published by a project partner. The combination offers an innovative approach to integrated refining and RE metal production. For Circuit 4, an extensive literature review was conducted for the processing of the CMs. The ultimate decision was to utilize a combined plasma and ionic liquid process as well to produce individual high-purity concentrates of Zn, Ni, Co, Mn, and Mg. A separate flowsheet for Li and Sr was recommended, which would yield carbonates of these elements. Due to the lack of suitable experimental data at this time, a process recommendation could not be provided but several methods have been proposed for consideration. Lastly, a techno-economic analysis (TEA) was conducted to assess the effectiveness of the proposed process for further investigation. The TEA results revealed a capital expense (CapEx) of $737 million and an annual operational expense (OpEx) of $220 million. Due to the selected elements, the hypothetical heap leach pad can produce 1 metric tonne per day of REO equivalent, but a conscious decision was made to only treat targeted REEs, resulting in the production of 0.4 metric tonne of REM. An estimated annual revenue of $90.87 million was projected based on standard market pricing information provided by the funding agency. During the TEA, ten different modules were evaluated for costing purposes. The precipitation circuit was identified as the largest single operational expense, followed by the Mg/Mn process due to the amount of treated metal. In terms of capital expenditures, the heap leach process incurred the highest cost, followed by the Mg/Mn process. The scalability of the plasma process is a crucial consideration since the reactors cannot be scaled beyond the largest demonstrated size due to their reliance on surface area of the slag and vapor phase. The purity estimate for the REEs are generally 98%±2% to produce a metal. The purity level being lower than the project objective was due to the lack of specific experimental data needed to tighten the tolerance of the estimates. Based on literature and previous experience, the CMs are estimated as follows; Ga (95%+, metal), Sr (95%+, carbonate), Li (95%+, carbonate), Ni (98%±2%, metal), Zn (95%+, metal sponge), Ge (95%+, metal), Co (98%±2%, metal), and Mn (98%±2%, metal).

01 COAL, LIGNITE, AND PEAT↗