Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Performance Trace”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Trace levels of PtRu on NiMo foam electrodes towards hydrogen evolution reaction

Electrochemical water splitting plays a critical role in high purity hydrogen production. To lower capital cost and energy consumption, efficient catalysts are required to boost the hydrogen evolution reaction (HER), especially in alkaline media. Although PtRu alloy materials are still considered the state-of-the-art catalyst for the HER under alkaline media, large-scale application is hindered by its scarcity and high cost. NiMo alloy catalysts can be potential Pt-free alternatives to drive the sluggish kinetics of the HER under alkaline media, but their HER activity is still not ideal due to the limitation of Volmer step. Herein, we investigated the impact of trace levels of PtRu catalysts on NiMo electrodes towards the HER by incorporating advantages of both PtRu and NiMo catalysts to minimize cost and maintain a promising HER activity. Trace levels of PtRu catalysts from 0.0173 to 0.2648 mg cm −2 were loaded on commercial NiMo foam electrodes by a facile electrodeposition technique utilizing ppm concentration levels of PtRu in electrolyte. The results show that 0.0173 mg cm −2 PtRu loading on NiMo electrode can significantly enhance HER activity compared to the pristine NiMo foam electrode (Tafel slopes drop from ∼170 to ∼60 mV dec −1 ) due to the remarkable decline of energy barriers towards water dissociation and hydrogen adsorption. Meanwhile, the rate determining step (RDS) switches from the Volmer step to the Heyrovsky step when the PtRu loading increases to 0.0675 mg cm −2 . Promising catalytic stability is observed at 100 mA cm −2 over the course of 50 h with a PtRu loading of 0.1198 mg cm −2 . In conclusion, this work demonstrates a potential strategy to decrease the cost and simultaneously maintain superior HER performance by integrating trace levels of PtRu catalyst with NiMo foam electrodes.

Electrochemical water splitting↗
DOE OSTI2023

A numerical study of the spectral radiative properties of packed bed with mixed bauxite and silica spheres

Bauxite and silica particles have gained increasing attention for applications in the field of concentrated solar power. In this work, a Monte Carlo ray-tracing simulation is performed to predict the radiative properties (absorptance, reflectance, and transmittance) of packed beds with mixed bauxite and silica spherical particles at wavelengths of 0.5, 2.6, 9.0, and 9.35 µm. These wavelengths are representative for the visible, near-infrared, and the mid-infrared regions that are important for solar and thermal radiation, respectively. A repeating unit-column approach is used to mathematically represent the particle bed. The effects of particle mixing ratio, volume fraction, and wavelength on the predicted radiative properties are examined. The obtained radiative properties are inputted to an inverse method to retrieve the effective absorption and scattering coefficients as well as the scattering albedo, which may be later used in a continuous-scale radiative heat transfer analysis. Furthermore, the independent scattering model is used to obtain the absorption coefficient and scattering albedo based on the absorption and scattering cross sections predicted by a Monte Carlo algorithm for a single particle. It is shown that the independent scattering model underpredicts the scattering coefficient for opaque particles but overpredicts the scattering coefficient for semitransparent particles for sufficiently high particle volume fractions. Here, the radiative properties calculated from the independent scattering model are compared to the full Monte Carlo simulation of the particle bed to examine the influence of particle mixing on dependent scattering.

42 ENGINEERING↗
DOE OSTI2026

Hydrogen Analysis by Gas Chromatography–Mass Spectrometry

The detection of hydrogen in complex gas mixtures is essential for many applications. Conventional approaches such as gas chromatography (GC) with thermal conductivity detection (TCD) and residual gas analyzers (RGAs) face significant limitations: TCD exhibits poor response when helium is used as a carrier gas, whereas RGAs lack chromatographic separation, preventing reliable quantification of hydrogen because of interference from other species. Traditional GC methods rely on dual-column configurations with packed or molecular-sieve porous layer open tubular (PLOT) columns to separate hydrogen from O2, N2, CO, CO2, CH4, and other hydrocarbons, increasing system complexity and limiting compatibility with mass spectrometry (MS) detection. In this work, we developed and validated a robust GC–MS method capable of directly detecting and quantifying hydrogen using electron ionization (EI) without dopants, reagent gases, or ion–molecule reaction schemes. By integrating a modified EI source and a cryogenically cooled single capillary column configuration, we achieved baseline separation of hydrogen from all major permanent gases and hydrocarbons in a refinery gas mixture using helium as the carrier gas. The method demonstrated excellent linearity, high sensitivity, and exceptional reproducibility. Adjustable sample-loop volumes and split ratios enabled optimization of peak shape and signal-to-noise performance for trace-level and percent-level hydrogen concentrations. Beyond hydrogen quantitation, the method provides simultaneous compositional profiling of other gases in a mixture in a single run, making it valuable for a wide range of tasks.

Lobodin, Vlad [ORNL]↗
DOE OSTI2022

The Ghost of Performance Reproducibility Past

The importance of ensemble computing is well established. However, executing ensembles at scale introduces interesting performance fluctuations that have not been well investigated. In this paper, we trace our experience uncovering performance fluctuations of ensemble applications (primarily constituting a workflow of GROMACS tasks), and unsuccessful attempts, so far, at trying to discern the underlying cause(s) of performance fluctuations. Is the failure to discern the causative or contributing factors a failure of capability? Or imagination? Do the fluctuations have their genesis in some inscrutable aspect of the system or software? Does it warrant a fundamental reassessment and rethinking of how we assume and conceptualize performance reproducibility? Answers to these questions are not straightforward, nor are they immediate or obvious. We conclude with a discussion about the performance of ensemble applications and ruminate over the implications for how we define and measure application performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fast HARDI Uncertainty Quantification and Visualization with Spherical Sampling

In this paper, we study uncertainty quantification and visualization of orientation distribution functions (ODF), which corresponds to the diffusion profile of high angular resolution diffusion imaging (HARDI) data. The shape inclusion probability (SIP) function is the state‐of‐the‐art method for capturing the uncertainty of ODF ensembles. The current method of computing the SIP function with a volumetric basis exhibits high computational and memory costs, which can be a bottleneck to integrating uncertainty into HARDI visualization techniques and tools. We propose a novel spherical sampling framework for faster computation of the SIP function with lower memory usage and increased accuracy. In particular, we propose direct extraction of SIP isosurfaces, which represent confidence intervals indicating spatial uncertainty of HARDI glyphs, by performing spherical sampling of ODFs. Our spherical sampling approach requires much less sampling than the state‐of‐the‐art volume sampling method, thus providing significantly enhanced performance, scalability, and the ability to perform implicit ray tracing. Our experiments demonstrate that the SIP isosurfaces extracted with our spherical sampling approach can achieve up to 8164× speedup, 37282× memory reduction, and 50.2% less SIP isosurface error compared to the classical volume sampling approach. We demonstrate the efficacy of our methods through experiments on synthetic and human‐brain HARDI datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Changing-look Active Galactic Nuclei Powered by Disk Tearing

Changing-look active galactic nuclei (CLAGN) feature order-of-magnitude variability in both the continuum and broad line luminosities on months-to-years-long timescales, and are currently unexplained. Simulations have demonstrated that rotating black holes (BHs) sometimes tear apart tilted accretion disks. These tearing events violently restructure the disk on timescales much shorter than a viscous timescale, hinting at a connection to CLAGN. Here, we show that disk tearing can power changing-look events. We report synthetic observations of an extremely high-resolution, three-dimensional, general-relativistic magnetohydrodynamic simulation of a geometrically thin, tilted accretion disk around a rapidly rotating, 10 8 M ⊙ BH. We perform ray-tracing calculations that follow the thermal disk emission to both a line-of-sight camera and to a distribution of cameras in a prescribed torus-like broad-line region (BLR). The continuum photoionizes the BLR, and we calculate the resulting spectrum. Both the continuum and line luminosities undergo order-of-magnitude swings on months-to-years-long timescales. We find shorter, week-long variability driven by the geometric precession of the inner disk and an intraday quasiperiodic oscillation driven by radial breathing of the inner disk. When the torn disk precesses, it causes asymmetric illumination of the BLR, driving time-evolving red-to-blue asymmetries of the broad emission lines that may be a smoking gun for disk tearing. We also make predictions for future photometric observations from ULTRASAT and Vera Rubin Observatory, given a Shakura–Sunyaev–like effective temperature profile, both of which may play an important role in detecting future changing-look events. Future work should consider nonthermal emission, especially from the time-variable corona.

accretion↗
DOE OSTI2020

Flexspander: augmenting expander networks in high-performance systems with optical bandwidth steering

Communication efficiency is one of the deciding factors in determining many of today’s high-performance computing (HPC) applications. Traditionally, HPC systems have been on static network topologies, making them inflexible to the variety of skewed traffic demands that may arise due to the spatial locality inherent in many applications. To handle traffic locality, researchers have proposed integrating optical circuit switches (OCSs) into the network architecture, which reconfigures the network topology to alter and dynamically adapt to the predicted traffic. In this paper, we present a novel reconfigurable network topology called Flexspander. Beyond offering a flexible interconnect, Flexspander also offers full flexibility in terms of construction and can be built with any arbitrary combination of commercial electrical packet switches and OCSs. We evaluate Flexspander performance through extensive simulations with multiple network traces, and our results show improved performance for Flexspander over currently proposed static and reconfigurable topologies in terms of the flow completion time.

Teh, Min Yee↗
DOE OSTI2023

Trace Element Sampling and Partitioning Modeling to Estimate Wastewater Composition and Treatment Performance at Coal Generators

The coal combusted at coal-fired power plants (CFPPs) contains several trace elements that impact human heath if released to the environment. Effective combinations of air pollution control devices (APCDs) and water pollution control devices (WPCDs) are critical to reducing these impacts. This work sampled pulverized CFPPs owned and operated by Louisville Gas & Electricity – Kentucky Utilities (LGE-KU) during baseload and cycling conditions. Next, an open-source and easy-to-implement trace element partitioning model was developed using publicly available datasets, literature studies of trace element partitioning, and sampling data to estimate trace element partitioning in APCDs to the gas, liquid, and solid phases exiting CFPPs. These estimates of the liquid phase trace element concentration in FGD wastewater were then used to estimate trace element behavior in WPCDs and evaluate treated wastewater effluent concentrations for compliance with the Effluent Limitation Guidelines (ELG) for the Steam Electric Generating Sector. Finally, cost estimates were developed for established and emerging wastewater treatment trains to identify the most cost-effective approaches for ELG compliance.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

I/O Bottleneck Detection and Tuning: Connecting the Dots using Interactive Log Analysis

Using parallel file systems efficiently is a tricky problem due to inter-dependencies among multiple layers of I/O software, including high-level I/O libraries (HDF5, netCDF, etc.), MPI-IO, POSIX, and file systems (GPFS, Lustre, etc.). Profiling tools such as Darshan collect traces to help understand the I/O performance behavior. However, there are significant gaps in analyzing the collected traces and then applying tuning options offered by various layers of I/O software. Seeking to connect the dots between I/O bottleneck detection and tuning, we propose DXT Explorer, an interactive log analysis tool. In this paper, we present a case study using our interactive log analysis tool to identify and apply various I/O optimizations. We report an evaluation of performance improvement achieved for four I/O kernels extracted from science applications.

Bez, Jean Luca↗
DOE OSTI2025

Analysis of Degradation and Electrochemical Behavior of Lithium–Oxygen Batteries under Lean Electrolyte Conditions

Lithium-oxygen batteries have garnered much attention in the past decades due to their high energy density and active material abundancy. Decreasing the electrolyte volume is critical to achieving a high energy density for practical applications. Recent works have demonstrated a serious performance limitation at lower electrolyte fillings. Here, in this work, we study specifically the difference between flooded and lean electrolyte lithium-oxygen batteries by systematically tracing the source of performance degradation. We find that the key contributor to the decrease in performance under lean electrolyte conditions stems from the lithium metal anode rather than the cathode. Solid electrolyte interphase (SEI) compositionally appears to be the same from X-ray photoelectron spectroscopy but the coverage is less uniform, and impedance is much higher under lean electrolyte conditions. Such an effect likely produced the observed poorer cycle performance at 10 μL cm -2 (similar to 28% cathode's pore volume filled) vs 100 μL cm -2 (278% of cathode's pore volume) at a capacity of 1.0 mAh cm -2 (0.1 mAh cm -2 ) using an electrolyte composition of 1 M LiTFSI in TEGDME.

Córdoba, Daniel [Argonne National Laboratory (ANL)↗
DOE OSTI2022

Evaluation of digestively resistant or soluble fibers, short- and medium-chain fatty acids, trace minerals, and antibiotics in nonchallenged nursery pigs on performance, digestibility, and intestinal integrity

Abstract Three experiments (EXP) were conducted to determine the effect of feed additives on performance, intestinal integrity, gastrointestinal volatile fatty acids (VFA), and energy and nutrient digestion in nonchallenged nursery pigs. In EXP 1, 480 pigs (6.36-kg body weight, BW) were placed into 96 pens with 5 pigs/pen, and allotted to 1 of 10 dietary treatments: 1) negative control containing no feed additive (NC), 2) NC + 44 mg chlortetracycline and 38.5 mg tiamulin/kg diet (CTsb), 3) NC + 5% resistant potato starch (RSpo), 4) NC + 5% soluble corn fiber (SCF), 5) NC + 5% sugar beet pulp (SBP), 6) NC + 0.30% fatty acid mix (FAM), 7) NC + 0.10% phytogenic blend of essential oils and flavoring compounds (PHY), 8) NC + 50 mg Cu and 1,600 mg zinc oxide/kg diet (CuZn), 9) NC + 5% resistant corn starch (RScn), and 10) NC + 0.05% β-glucan (BG) for 28 d. There was no impact of dietary treatment on BW gain or feed intake (P ≥ 0.22). Pigs fed diets containing SCF, CTsb, and RSpo resulted in microbial community differences compared to pigs fed the NC (P < 0.05). In EXP 2, 48 barrows (12.8 kg BW) were selected at the end of EXP 1 and fed the same dietary treatments they had previously received: 1) NC, 2) NC + 5% RScn, 3) NC + 5% SCF, and 4) NC + FAM for 8 d. There was no effect of feeding diets containing RScn, SCF, or FAM on in vivo intestinal permeability (P ≤ 0.21). Ileal or colon pH, concentrations of VFA did not differ due to dietary treatment (P ≥ 0.36), but pigs fed diets containing FAM resulted in a greater butyric acid concentration in the cecum compared to pigs fed the NC (P ≤ 0.05). In EXP 3, 156 pigs (6.11 kg BW) were placed into 52 pens with 3 pigs/pen and allotted to 1 of 4 dietary treatments arranged in a factorial manner: 1) NC, 2) NC + 5% RSpo, 3) NC + 0.30% FAM, and 4) NC + 5% RSpo + 0.30% FAM for 24 d. Feeding pigs diets containing RSpo did not affect BW gain (P = 0.91) while pigs fed diets containing FAM grew improved BW gain (P = 0.09). Colonic butyric acid concentrations were greater in pigs fed diets containing RSpo (P = 0.03), while pigs fed diets containing FAM exhibited reduced total VFA concentrations (P = 0.11). The results indicate that supplementing diets with digestively resistant but fermentable fibers, short- and medium-chain fatty acids, or antibiotics do not have a consistent effect, positive or negative, on markers of intestinal integrity or barrier function, intestinal VFA patterns, ATTD of energy and nutrients, or on pig performance.

Agriculture↗
DOE OSTI2021

General relativistic MHD simulations of non-thermal flaring in Sagittarius A*

Sgr A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics code H-AMR to perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing code bhoss to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ~60 h light curves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the light curves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetized accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms per cent amplitudes, $\gtrsim 150{{\ \rm per\ cent}}$ both in the NIR and the X-rays, with changes in the accretion rate driving the 230 GHz flux variability, in agreement with Sgr A* observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

High-Burnup BWR LOCA Burst Analysis Using High-Fidelity Multiphysics Simulations

The US nuclear industry is looking to improve on the operating economics of the current fleet of light-water reactors (LWRs). One way of achieving this is by operating fuel to higher burnup. In pressurized water reactors (PWRs), relaxing the current burnup limit will allow for cycle length extensions and power uprates; in boiling water reactors (BWRs) it may allow for improved fuel utilization and reduced feed assemblies, as well as more efficient power uprates and increased capacity factors that will support the Administration’s Executive Order to facilitate 5 GW of power uprates at existing nuclear facilities. However, one of the key limitations to operating fuel to higher burnup is the risk of fuel fragmentation, relocation, and dispersal (FFRD). Recognizing the high interest in extending burnup limits, the US Nuclear Regulatory Commission (NRC) has issued Draft Regulatory Guide DG-1434, which defines an approach that would be acceptable to the NRC for addressing FFRD risk. The approach defined will require better understanding of the phenomena leading to FFRD as well as best-estimate simulation methods to understand FFRD risk in high-burnup cores. The Nuclear Energy Advanced Modeling and Simulation program is supporting the FFRD industry challenge problem through development of state-of-the-art, high-fidelity modeling and simulation LWR analysis capabilities; namely, the BISON fuel performance code and the VERA core simulator software. These tools, along with the US NRC TRACE system analysis code, have been utilized for analysis of FFRD risk in both PWR and BWR cores in recent years. The work documented in this report addresses the lack of high-fidelity research for BWRs and builds on a previous activity where the framework has been applied to Cycles 16 through 18 of Limerick Unit 1, a BWR/4, with introduction of 8 high-burnup lead use assemblies (HBLUAs) that were representative of the 8 HBLUAs loaded into Limerick Unit 2 in 2021. VERA was used in this previous activity to model rod-by-rod depletion in these cycles, and its solution was used to initialize a TRACE simulation of a large-break loss-of-coolant accident (LBLOCA) at the end of Cycle 18. In the work documented in this report, the TRACE model was improved by refining the core mesh and utilizing a new feature that allows for capturing the full 3D VERA power distribution in the model. This allows for a more detailed solution for setting BISON boundary conditions. Furthermore, the solutions from VERA and TRACE were used to set up and perform BISON simulations of about 1,000 rods sampled from the core, including all burnup levels. Utilizing two cladding burst models, it was shown that no fuel rods were predicted to burst during the postulated LBLOCA transient. Additionally, a sensitivity study was performed by artificially increasing linear heat rate during the postulated LBLOCA to identify parameters that correlate with rod burst susceptibility. Burnup, fission gas release, and hoop strain were all found to be positively correlated with rod burst susceptibility. Small-break loss-of-coolant accident (SBLOCA) analyses were also performed; these analyses predicted cladding temperature increases that were bounded by the LBLOCA cladding temperatures for all small break sizes studied for this plant. However, future refinements to the plant response assumptions during the SBLOCA could impact the predicted cladding response. Finally, a benchmark study was performed between CTF and TRACE for LOCA conditions to better qualify CTF for BWR LOCA modeling.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Coupled Decay Heat and Thermal Hydraulic Capability for Loss-of-Coolant Accident Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle for longer durations. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than has traditionally been allowed in these reactors. Thus, there are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al.. To investigate these concerns, advanced modeling and simulation capabilities are being leveraged to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to more accurately model these phenomena; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs for the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required. Prior LOCA simulations have relied on the TRACE systems code, which can model the core and primary loop during accident conditions. TRACE includes many models for various aspects of the primary loop, but two sets of models are important for this report. First, TRACE uses a lumped-fuel approach for modeling the core. This approximates the ~50,000 fuel rods in the core with a much smaller number of rods. The rods can be lumped in various ways as determined by the user. For example, one lumped rod may be used to represent all rods in an assembly, sometimes with an additional rod representing the hottest fuel rod. However, due to runtime constraints and complexity of modeling, a more common approach is to group several assemblies or larger regions of the core into single lumped rods. These lumping schemes apply not only to fuel rods but to flow channels as well. Second, TRACE has several different models for treating decay heat, ranging from pregenerated decay heat curves based on an ANSI/ANS-5.1 standard (hereinafter abbreviated simply as ANSI) to explicit time-dependent heat inputs from the user. None of these models account for differences in isotopics between different rods, which is an approximation the work in this report seeks to eliminate. This report focuses on the implementation of coupled decay heat capabilities in the Virtual Environment for Reactor Applications (VERA) code suite to address a gap identified in previous LOCA simulations. This constitutes an improvement for both the lumped-fuel and decay heat models in TRACE. VERA has been developed to perform high-fidelity, whole-core multiphysics simulations for LWRs. Previously, during the Consortium for Advanced Simulation of LWRs (CASL) program, the emphasis was on providing accurate steady-state analysis—with a secondary focus on reactivity insertion accident (RIA) analysis—to address operational challenges in the nuclear energy industry. Under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, these capabilities are being extended to a broader range of transient analyses with the goal of quantifying the risk of fuel failures such as FFRD. To properly model such conditions with VERA, decay heat calculations have been integrated with the multiphysics to enable rod-by-rod thermal hydraulic (TH) conditions to be driven by the decay heat in long-running accidents such as LOCAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Polycrystalline SnSe with a thermoelectric figure of merit greater than the single crystal

Abstract Thermoelectric materials generate electric energy from waste heat, with conversion efficiency governed by the dimensionless figure of merit, ZT. Single-crystal tin selenide (SnSe) was discovered to exhibit a high ZT of roughly 2.2–2.6 at 913 K, but more practical and deployable polycrystal versions of the same compound suffer from much poorer overall ZT, thereby thwarting prospects for cost-effective lead-free thermoelectrics. The poor polycrystal bulk performance is attributed to traces of tin oxides covering the surface of SnSe powders, which increases thermal conductivity, reduces electrical conductivity and thereby reduces ZT. Here, we report that hole-doped SnSe polycrystalline samples with reagents carefully purified and tin oxides removed exhibit an ZT of roughly 3.1 at 783 K. Its lattice thermal conductivity is ultralow at roughly 0.07 W m –1 K –1 at 783 K, lower than the single crystals. The path to ultrahigh thermoelectric performance in polycrystalline samples is the proper removal of the deleterious thermally conductive oxides from the surface of SnSe grains. These results could open an era of high-performance practical thermoelectrics from this high-performance material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive “full-wave” simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. Furthermore, the numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves on previously derived analytical approximations of the MGO integral.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Characterizing the Impact of Heliostat Size, Focus Method, and Optical Error on Concentrating Solar Tower Systems

Concentrating Solar Power (CSP) tower systems use heliostat fields to direct solar energy to a central receiver. The efficiency of the heliostat field is crucial for achieving thermal generation targets and is also the largest cost component, making up 30-40% of total project costs. Developers face a tradeoff: reducing heliostat costs may decrease optical performance and increase land use, while improving performance could raise manufacturing costs (through higher precision). The optimal design depends on land, labor, and manufacturing costs. Additionally, inherent performance limitations exist - optical errors result in larger solar images as heliostats move farther from the receiver, increasing spillage and attenuation losses. Ultimately, there exists a point where receiver design power can no longer be achieved at design conditions for a given solar field area. This point depends on key heliostat field design parameters including field layout, heliostats size, focus method, and heliostat optical error. In this work, we explore and characterize the performance trade-offs between heliostat size, focus method, and optical error using Monte Carlo ray-tracing simulations using NREL's high performance computing resources. The goal of this fundamental investigation is to understand practical heliostat design performance impacts that can be used to produce more cost-effective heliostat fields while still achieving plant thermal generation targets.

14 SOLAR ENERGY↗