Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pointing Performance Analysis and Requirements”

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 235 records · Page 13

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

2024 Buildings Technology Baseline: Dataset Documentation

The Buildings Technology Baseline is a curated and regularly updated dataset of current and projected performance, retail, and installed price data for all major building energy technologies needed to enable cost/benefit analyses. Building technology analyses require an up-to-date understanding of installation costs and cost-effectiveness of key building energy efficiency technologies. The dataset was assembled by Guidehouse during fiscal year 2024. Data was gathered from the 2024 National Residential Efficiency Measures Database (NREMDB), the 2023 Energy Information Administration Updated Buildings Sector Appliance and Equipment Costs and Efficiencies ("EIA Building Data Report"), DOE Lighting Market Model, the 2023 RSMeans database, and the 2020 Grid-Interactive Efficient Building Technology Cost, Performance, and Lifetime Characteristics ("GEB Data Report"), Lawrence Berkeley National Laboratory, various literature, as well as new data from online retailers, stakeholder interviews, and contractor databases in 2023 and 2024. The dataset has been reviewed by subject matter experts at NREL and DOE. The 2024 dataset release is intended to be a starting point for interested users to provide feedback. This database is not intended to provide specific cost estimates for a specific project. The cost estimates do not include any rebates or tax incentives that may be available for the measures. Rather, it is meant to help determine which measures may be more cost-effective. The National Renewable Energy Laboratory (NREL) makes every effort to ensure accuracy of the data; however, NREL does not assume any legal liability or responsibility for the accuracy or completeness of the information.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Cyber risk assessment and investment optimization using game theory and ML-based anomaly detection and mitigation for wide-area control in smart grids

The electric power grid is increasingly becoming susceptible to cyber attacks that exploit vulnerabilities in the smart grid control, information, and physical layers. Successful cyber attacks can have catastrophic impacts on the social and economic well-being of any nation all over the globe. It has, thus, become imperative to secure the smart grid against such adversarial actions to ensure stable, secure, and reliable operation of the grid. The existing research and industry practices prove to be inadequate in terms of providing pragmatic and effective defense methodologies and measures for long-term cybersecurity planning and real-time cybersecurity for grid operation. For example, existing works lack models that incorporate uncertain behavior of cyber-attackers and pragmatic defense measures for cyber risk assessment and cybersecurity investment optimization which often provide unreliable and strictly qualitative solutions to these problems. At the same time, with the growing number of cyber incidents in the grid, there still exists a need to develop attack-resilient algorithms for wide-area monitoring, protection, and control (WAMPAC) applications like the wide-area voltage control systems (WAVCS) for Flexible AC Transmissions Systems (FACTS) that lack in scalable and feasible solutions from the cybersecurity perspective. This dissertation proposes novel models and methodologies for: (1) Cybersecurity planning, and (2) Cybersecurity for system operation. The cybersecurity planning is achieved through cyber risk assessment and cybersecurity resource investment optimization for long-term cybersecurity of the grid using game theory and attack-defense trees. Cybersecurity for system operation consists of development of cyber anomaly detection and mitigation algorithms for flexible AC transmission system (FACTS) controller-based wide-area voltage control systems (WAVCS) using machine learning (ML), and software defined networking-based moving target defense network routing for achieving real-time cyber-physical security for grid operations. This is followed by hardware-in-the-loop (HIL) implementation and evaluation of these attack prevention, detection, and mitigation algorithms and methodologies showcasing their feasibility in a close to real-world environment. For cybersecurity planning, a novel approach involving a combination of game theory and attack defense trees (ADT) for optimal cybersecurity resource allocation in the smart grid is proposed. This methodology involves modeling of the cyber-physical smart grid substations as ADTs, defining attacker costs, defense costs, and attack probabilities for attack access points. Using game theoretical formulation, optimal defense strategies for the defender of the system to invest cybersecurity resources in the grid are obtained. Additionally, a game-theoretic framework is developed for quantitative cyber-physical risk assessment of the grid under a dynamically changing cyber threat space and uncertain behavior of cyber attackers which is further used to optimize investments in the smart grid's cybersecurity resources. The attacker, defender, and the smart grid system are modeled while incorporating attacker-stochasticity and federal guidelines for smart grid cybersecurity. This allows quantification of threat, vulnerabilities, and attack impact of the grid for quantitative risk assessment. The defender's budget to invest in the security resources in the grid is optimized based on the strategies leading to minimum system risk. The evaluation of the proposed solutions highlight the feasibility for practical implementation of these methodologies and algorithms in the smart grid, while taking the federal requirements and guidelines for smart grid security into consideration. For achieving cybersecurity for system operation, attack prevention, detection, and mitigation algorithms and methodologies are developed specifically for FACTS-based WAVCS. Anomaly detection and mitigation in the WAVCS are achieved using algorithms based on machine learning which involves offline training and testing of ML models with CPS datasets incorporating physics-based features that allow accurate distinction between system faults and cyber attacks. For attack prevention, a methodology based on software defined network (SDN)-based moving target defense (MTD) network routing is proposed that enables prevention of Denial of Service (DoS) type attacks on the smart grid communication system. Subsequently, these methodologies and algorithms are implemented and evaluated on an HIL testbed that allows for real-time attack prevention, detection, and mitigation of emulated cyber attacks on the WAVCS in a close to real-world environment. The results show highly accurate and efficient performance of the implemented algorithms and methodologies with the smart grid system operating within the NERC's system operation limits even in the presence of DoS and data integrity cyber attacks. This work opens up future research opportunities in other directions such as (1) Expanding cybersecurity planning methodologies to real-time cyber contingency analysis with different game formulations; and (2) Applying the cybersecurity for system operation algorithms to broader categories of wide-area control applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗

Techno-Economic Analysis of a Concentrating Solar Power Plant Using Redox-Active Metal Oxides as Heat Transfer Fluid and Storage Media

We present results for a one-dimensional quasi-steady-state thermodynamic model developed for a 111.7 MW e concentrating solar power (CSP) system using a redox-active metal oxide as the heat storage media and heat transfer agent integrated with a combined cycle air Brayton power block. In the energy charging and discharging processes, the metal oxide CaAl 0.2 Mn 0.8 O 2.9-δ (CAM28) undergoes a reversible, high temperature redox cycle including an endothermic oxygen-releasing reaction and exothermic oxygen-incorporation reaction. Concentrated solar radiation heats the redox-active oxide particles under partial vacuum to drive the reduction extent deeper for increased energy density at a fixed temperature, thereby increasing storage capacity while limiting the required on sun temperature. Direct counter-current contact of the reduced particles with compressed air from the Brayton compressor releases stored chemical and sensible energy, heating the air to 1,200°C at the turbine inlet while cooling and reoxidizing the particles. The cool oxidized particles recirculate through the solar receiver subsystem for another cycle of heating and reduction (oxygen release). We applied the techno-economic model to 1) size components, 2) examine intraday operation with varying solar insolation, 3) estimate annual performance characteristics over a simulated year, 4) estimate the levelized cost of electricity (LCOE), and 5) perform sensitivity analyses to evaluate factors that affect performance and cost. Simulations use hourly solar radiation data from Barstow, California to assess the performance of a 111.7 MW e system with solar multiples (SMs) varying from 1.2 to 2.4 and storage capacities of 6–14 h. The baseline system with 6 h storage and SM of 1.8 has a capacity factor of 54.2%, an increase from 32.3% capacity factor with no storage, and an average annual energy efficiency of 20.6%. Calculations show a system with an output of 710 GWh e net electricity per year, 12 h storage, and SM of 2.4 to have an installed cost of $329 million, and an LCOE of 5.98 ¢/kWh e . This value meets the U.S. Department of Energy’s SunShot 2020 target of 6.0 ¢/kWh e ( U. S Department of Energy, 2012 ), but falls just shy of the 5.0 ¢/kWh e 2030 CSP target for dispatchable electricity ( U. S Department of Energy, 2017 ). The cost and performance results are minimally sensitive to most design parameters. However, a one-point change in the weighted annual cost of capital from 8 to 7% (better understood as a 12.5% change) translates directly to an 11% decrease (0.66 ¢/kWhe) in the LCOE.

Gorman, Brandon T.↗

Accelerated Aging Humidity Chamber for Nuclear Grade HEPA Filter Media

A unique accelerated aging humidity chamber for simulating the natural aging process of nuclear grade high efficiency particulate air (HEPA) filter media is being designed, implemented, and tested. The nuclear industry currently lacks information regarding the shelf life and service life of HEPA filter media. The Department of Energy recommends disposing of HEPA filters 10 years after the manufacture date for filters operated under dry conditions. Studies have shown the tensile strength and water repellency of HEPA filter media diminish with age. The chamber's principal function is to sustain an elevated relative humidity condition for an extended duration specified by the Arrhenius aging model. To accomplish this, various design parameters were established: (1) consistently achieve uniform humidity conditions within the chamber to equally expose media sheets, (2) ensure the chamber is sufficiently sealed to mitigate the uncontrolled infiltration of air, (3) confidently secure the media sheets to avert the introduction of unnecessary stressing/creasing, (4) fabricate the chamber from clear polycarbonate material to allow inspection of the media sheets, (5) provide a safe and ergonomic design for personnel, (6) ensure the chamber design is replicable and relatively easy to fabricate/assemble. A humidity source, fan, PID controller, humidity and temperature transmitters, and data loggers are required for proper operation and control of the chamber. The humidity source must accommodate target values up to 95% relative humidity throughout the chamber. The three humidity and temperature transmitters are strategically spaced throughout the chamber to ensure uniform conditions. The 10 cubic feet per minute (CFM) fan selected intends to accomplish the desired air change rate inside the chamber of approximately 20 air changes per hour (ACH). The 2 plenums on the upper and lower portions of the chamber utilize a perforated design to uniformly distribute air, promote air mixing, and control the air velocity entering the chamber. Industrial blueprint hanging clamps are being retrofitted to be seated within the chamber drawers and firmly secure the nuclear grade filter media. The aged media will be evaluated using autopsy methodology, and the results obtained intend to help clarify the useful life of nuclear grade HEPA filters. Nuclear grade high efficiency particulate air (HEPA) filters are defined as disposable, extended-media, dry-type filters with a rigid casing enclosing the full depth of the pleats and have a minimum particle removal efficiency of 99.97%. [1] - HEPA filters are constantly exposed to stressors and subsequently have been shown to degrade as the filter ages. Properties such as tensile strength and water repellency, along with others, are negatively impacted. [2] - Due to the scarcity of naturally aged HEPA filters available, it is crucial to develop accelerated aging methods which effectively mimic the natural degradation effects. - After accelerated aging, the filters are subject to testing in the Axial Flow Large Scale Test Stand (ALSTS) and other analysis techniques such as the same qualification tests performed directly after being manufactured. - Quality data is necessary for future decisions regarding the lifetime of nuclear grade HEPA filters. - The objective is to develop a prototype accelerated aging humidity chamber for nuclear grade HEPA filter media. The chamber should expose the media to commonly occurring stressors to expedite the degradation process caused by aging. PID yields reliable and reproducible results: - CFD model correlates well with actual performance; - Further analysis is needed to complete characterization; - Future work for the chamber is still to be completed. The data shown in Figure 5 and Table 1 indicate the system is functioning as designed. As can be determined from the system response evaluation in Table 1, the PID controller is yielding reliable and reproducible results at each set point. As the set point increases, the time constant and 20% settling time increase accordingly. From the results in Table 1, it stands to reason that the CFD model correlates well with the actual performance of the system. Each system response parameter for the 90%RH* CFD model is slightly lower than the parameters collected from the physical model tested at an equivalent set point. This is to be expected considering the inlet in the CFD model is a constant supply of water vapor, whereas the physical model uses a PID controller to monitor the amount of water vapor needed and reacts accordingly. More analysis is still to be done on the prototype chamber. A pressure test will be conducted in order to determine the hourly leak rate of the chamber. This test will involve constant air flow into the chamber until the desired pressure is reached. A uniformity test will also be conducted in order to ensure each filter medium is being exposed to equivalent relative humidity levels. This test will involve an array of humidity meters strategically oriented inside the chamber to record any possible gradients. Future work for the chamber includes dehumidification abilities, tests while fully loaded with filters, and temperature control. Once completed, the autopsy team at Institute of Clean Energy Technology will conduct necessary accelerated aging studies as needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Preliminary Design of Reactor Cavity Cooling System for a Horizontal Compact HTGR

The Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed MIGHTR in 3 years and support its commercialization as a safe and low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. This report documents the preliminary design study of the RCCS for the HC-HTGR. It includes the establishment of the design requirements, a high-level design study by initial scoping calculations, and preliminary performance calculations of the HC-HTGR RCCS design. Design requirements for the HC-HTGR RCCS have been established to guide preliminary design activities and scoping performance calculations. Initial scoping calculations including estimation of the water inventory, estimation of HVAC thermal capability, and a parametric study on loop dimensions by standalone RCCS analysis. Based on scoping calculation results, a set of baseline dimensions of the HC-HTGR RCCS was derived. A water panel modeling approach was investigated to explore various potential design options for the water panel under consideration for the HC-HTGR RCCS using RELAP5-3D. A test case study was performed to assess the prediction capability of two modeling approaches. The results were compared with CFD simulations conducted in constant RPV temperature and heat flux boundary conditions. It confirms the capability of the RELAP5-3D modeling approach to include all important heat transfer mechanisms expected in the HC-HTGR RCCS operation conditions. Then, a reference RELAP5-3D model for the 1/8 th of a compartment of the preliminary design of the HC-HTGR RCCS was developed. A preliminary performance analysis was conducted to evaluate single-phase natural circulation performance with different top tank temperature values and panel conduction performance in various operation conditions. From single-phase natural circulation performance analysis, the system operation mode was investigated in normal operating and limiting design conditions. It showed operation mode in a subcooled state with a proper top tank water cooling system. Parasitic heat loss by both internal air flow and RCCS was estimated, showing it satisfies maintaining below target maximum heat loss of the HC-HTGR RCCS. From the panel conduction performance analysis, two candidate materials for the riser tube such as carbon steel and stainless steel were compared in the thermal performance of HC-HTGR RCCS. From a single water panel test compared with CFD simulation results, it was confirmed that the current capability of the RELAP5-3D modeling approach for the water panel predicts the thermal conduction of two different materials of the water panel. Then, system-level thermal performance analysis was performed for 1/8 th of the compartment of the preliminary HC-HTGR RCCS design. It was first observed that the current preliminary HC-HTGR RCCS design had minimal impact on the overall thermal performance of the water panel by changing pipe material from carbon steel to stainless steel. From Argonne’s effort on the ongoing water-based NSTF testing program, several considerations other than the thermal performance point of view were addressed to be considered in selecting pipe materials.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Measurements and computational analysis of the natural decay of 176 Lu

Background: Mainly because of its long half-life and despite its scientific relevance, spectroscopic measurements of 176 Lu forbidden β decays are very limited and lack formulation of shape factors. A direct precise measurement of its Q value is also presently unreported. In addition, the description of forbidden decays provides interesting challenges for nuclear theory. The comparison of precise experimental results with theoretical calculations for these decays can help to test underlying models and can aid the interpretation of data from other experiments. Purpose: Perform the first precision measurements of 176 Lu β-decay spectra and attempt the observation of its electron capture decays, as well as perform the first precision direct measurement of the 176 Lu β-decay Q value. Compare the shape of the precisely determined experimental β spectra to theoretical calculations, and compare the end point energy to that obtained from an independent Q value measurement. Method: Here, the 176 Lu β-decay spectra measurements and the search for electron capture decays were performed with an experimental setup that employed lutetium-containing scintillator crystals and a NaI(Tl) spectrometer for coincidence counting. The β decay Q value was determined via high-precision Penning trap mass spectrometry (PTMS) with the LEBIT facility at the National Superconducting Cyclotron Laboratory. The β-spectrum calculations were performed within the Fermi theory formalism with nuclear structure effects calculated using a shell model approach. Results: Both β transitions of 176 Lu were experimentally observed and corresponding shape factors formulated in their entire energy ranges. The search for electron capture decay branches led to an experimental upper limit of 6.3×10 –6 relative to its β decays. The 176 Lu β-decay and electron capture Q values were measured using PTMS to be 1193.0(6) and 108.9(8) keV, respectively. This enabled precise β end point energies of 596.2(6) and 195.3(6) keV to be determined for the primary and secondary β decays, respectively. The conserved vector current hypothesis was applied to calculate the relativistic vector matrix elements. The β-spectrum shape was shown to significantly depend on the Coulomb displacement energy and on the value of the axial vector coupling constant gA, which was extracted according to different assumptions. Conclusion: The implemented self-scintillation method has provided unmatched observations of 176 Lu, independently validated by the first direct measurements of its β-decay Q value by Penning trap mass spectrometry. Theoretical study of the main β transition led to the extraction of very different effective gA and log10f values, showing that a high-precision description of this transition would require a realistic nuclear structure with nucleus deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GeoMelt{sup R} In-Container Vitrification (ICV){sup TM} for Fukushima Daiichi Water Treatment Secondary Wastes - 20212

The Japanese government is supporting development work implemented jointly by Veolia subsidiaries Kurion Japan, K.K., Veolia Nuclear Solutions (VNS), Inc., and Veolia Nuclear Solutions Federal Services, LLC for treating radioactive waste generated from Fukushima Daiichi Nuclear Power Station (NPS) water treatment using the GeoMelt{sup TM} In-Container Vitrification (ICV){sup TM} technology. The initial work consisted of glass formulation and engineering-scale testing which was completed in 2018, in the frame of an IRID (International Research Institute for Nuclear Decommissioning) program as part of a project subsidized by Japan's Ministry of Economy, Trade and Industry (METI). The Fukushima Daiichi NPS Mid- and Long-Term Road-map requires investigation of methods to stabilize solid wastes (and to immobilize radioisotopes in the wastes) generated as a result of emergency response and decommissioning activities. Cooling water treatment has resulted in a significant amount of solid and slurry secondary wastes (mostly adsorbents and ion-exchange materials) which will require processing at some point. GeoMelt{sup R} ICV{sup TM} is a joule-heated melter technology which uses a refractory-lined single-use container combining the melter and disposal container. There is no pouring required nor concerns with refractory corrosion which allows the process to accommodate a wide range of waste chemistries and high waste loadings. The testing described here consisted of three engineering-scale melts, each processing between 212 kg and 240 kg of waste simulants, glass formers, and non-radioactive cesium (Cs) and strontium (Sr) tracers. Continuous isokinetic stack sampling of off-gas emissions was performed for each test in order to calculate Cs and Sr retention in the glass wasteform. Single-pass retention of Cs in the final glass wasteform ranged from 91.46 to 99.30%, and single-pass retention of Sr ranged from 99.76 to 100%. Planned particulate recycle will increase these retention levels. Melt 1 processed a mixture of KUR-EH (a zeolite-based ion-exchange material), simulated Advanced Liquid Processing System (ALPS) Carbonate and Iron Slurries, and glass additives. Melt 2 processed a mixture of KUR-EH, KUR-TSG (a titanate-based adsorbent), and glass additives. Melt 3 processed a mixture of KUR-EH, simulated barium sulfate/iron ferrocyanide sludge (AREVA sludge), and glass additives. Waste loadings for these melts ranged from 70 weight percent (wt%) to 82 wt%. Vitrification produces a waste form much denser than the stored water treatment secondary waste wastes, resulting in significant volume reduction. Volume reductions for the three tests ranged from 74 to 79 vol%. Vitrification produces a chemically durable wasteform. Pacific Northwest National Laboratory (PNNL) tested three glass samples from each engineering-scale melt) by the Materials Characterization Center 1 (MCC-1) test, an international standard leach test of the chemical durability of nuclear waste glasses. PNNL also obtained one U.S. reference glass (EA Glass) and two Japan reference glasses P0798) and tested these under the same MCC-1 conditions (90 deg. C, 10 m-1, DIW, and 7, 14, 28-day) as the GeoMelt{sup R} ICV{sup TM} glasses. The GeoMelt{sup R} ICV{sup TM} glasses exhibited lower total normalized releases and 14- to 28-day normalized release rates than the three reference glasses. These results suggest that the GeoMelt{sup R} ICV{sup TM} glasses have durabilities on par with high-level waste glasses under standard test conditions. Post-melt process sampling and analysis indicated no Cs migration into the melter refractory materials and very little deposition of Cs or Sr onto the melter hood or off-gas piping internals. The results of the testing indicated good Cs retention in the glass, high volume reduction and waste loadings, and excellent chemical durability. These factors are important to minimize treatment costs and to protect workers and the environment. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology for pin lattice models is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Technical Performance of Refractory Liners for Molten Chloride Salt Thermal Energy Storage Systems

A chloride-based molten-salt system that uses a ternary blend of MgCl2/KCl/NaCl is investigated to provide higher temperature thermal energy storage capability. Despite higher thermal stability, molten chlorides present several unique challenges, including the design of internal refractory-ceramic liners to prevent the corrosion and thermal stress of alloy tank shells. This work discusses issues and potential solutions related to containment of molten chloride salt, specifically the optimization of the refractory material at the molten salt interface (hot face). The down-selected hot face candidate limits permeation of salt through the material and forms a highly stable secondary surface phase in equilibrium with the molten salt. A mortar is created using the corrosion resistant hot face brick. Brick and mortar composite are subjected to mechanical stress/strain analysis, in order to calculate composite material properties and better inform thermomechanical models. The U.S. Department of Energy Generation 3 (DOE Gen3) program seeks to develop higher efficiency CSP plants that can provide cost-competitive, flexible power in the U.S. electric grid. The proposed Gen3 Liquid Pathway CSP plant closely resembles the configuration of current nitrate salt power towers with two-tank storage (Gen2). The differences between Gen2 and Gen3 include the types of compatible materials used in salt storage tank construction. Stainless steel loses strength at Gen3 temperatures, and although nickel superalloys would be capable of withstanding sustained high temperatures, these materials are prohibitively expensive at scale. Uninsulated tank shells also pose a significant risk as common steels are highly susceptible to chemical attack from molten chloride salt. To address these concerns, refractory-ceramic based containment materials are proposed to line the inside of the hot and cold storage tanks. In doing so, stainless or carbon steel shells may be used in construction depending on the level of insulation provided. The composition of the internal liner requires careful consideration to maximize the efficacy of multiple parameters including corrosion resistance, strength at operating temperature, durability, and cost. This is particularly true for the material at the interface with the salt, known as the "hot face", which is responsible for protecting the insulating layers between the tank shell and the hot face brick layer. The molten salt in this system is superheated over 300 °C above its freeze temperature. Therefore, unlike other industrial processes which use refractory-lined vessels, it is not expected that a freeze plane will develop in the hot face. Therefore, the hot face must be designed to withstand chemical corrosion and inhibit permeation of molten salt into the insulating layers. A down selection was performed to identify a hot face candidate best equipped to maintain thermal, mechanical, and chemical integrity when exposed to molten salt over extended periods of time. Long-duration chemical capability experiments were conducted with the down selected hot face refractory fully immersed in molten chloride salt for up to 3000 hours. The average salt penetration does not exceed 100 microns. When extrapolated to 20 and 30 years of continuous exposure, the expected salt penetration depth is approximately 2.0 mm and 2.9 mm, respectively. A magnesium-rich secondary phases develops at the salt/refractory interface. X-ray diffraction identifies the material as forsterite (Mg2SiO4), which is reported to form synthetically in molten chloride salt solutions. These results suggest the selected hot face will adequately inhibit salt permeation. While there is optimism that the hot face brick will inhibit salt permeation, mortar joints are typically the weakest point of a refractory brick liner. From a thermochemical perspective, differing thermal expansion coefficients may result in the mortar and brick to grow independent of each other, creating gaps through which molten salt can penetrate. To address this issue, NREL has developed an in-house mortar composed of the down selected hot face brick that has been shown to be compatible with the salt. Compressive stress/strain analyses have been performed on the brick/ mortar composites to generate stress/strain curves. Modulus of elasticity and Poisson's ratio of the composite may be calculated from the stress/strain curves, in order to provide more representative data to finite element mechanical models for accurate approximation of stress on the tank shell and the amount of thermal expansion expected within the tank liner.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Macro-physical Properties of Shallow Cumulus from Integrated ARM Observations (Final Report)

Fair-weather shallow cumuli (ShCu) play an important role in many climate-related processes. Irregular geometry of ShCu and their strong temporal and spatial variability make it challenging to observe ShCu holistically and to represent them correctly in climate models. To improve ShCu parameterizations, information on both vertically and horizontally resolved cloud properties is required. Commonly, the vertically resolved cloud properties are provided by zenith pointing lidar-radar observations with a very narrow field of view (FOV). Thus, these “pencil-beam” properties may not be representative of a larger surrounding area. Limited number of areal-averaged cloud properties, such as fractional sky cover (FSC), are offered typically by wide-FOV observations. The main goal of our project was to integrate advantages of the narrow-FOV (vertical structure of clouds) and wide-FOV (spatial arrangement of clouds) observations for an improved characterization of single-layer ShCu observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site for an 18-yr period (2000-2017). There are four major accomplishments of our project, First, an updated operational cloud classification for days with ShCu has been suggested and evaluated through a detailed comparison with the manually curated records. Our classification extends successfully the latest ARM cloud type Value Added Product (VAP) based on the Active Remote Sensing of Clouds (ARSCL) cloud product by incorporating both cloud fraction (CF) provided by narrow-FOV ceilometer data and FSC from wide-FOV images offered by a Total Sky Imager (TSI). Moreover, our classification allows one to identify impact of instrumentation changes at the SGP site, namely the transition to KAZRARSCL with the updated cloud radar, on the identification of periods with single-layer ShCu. Second, a new approach that resolves cloud area distributions for a given region (up to 4x4 km 2 ) has been suggested and cloud equivalent diameters (CEDs) have been estimated for the first time. These estimations have been performed over a wide range of cloud sizes (about 0.01–3.5 km) with high temporal resolution (30s) using wide-FOV TSI images and cloud base height (CBH) provided by complementary narrow-FOV lidar measurements. Our simple and computationally inexpensive approach offers a previously unavailable dataset for process studies in the convective boundary layer and evaluation of ShCu parameterizations in cloud-resolving models. Third, a long-term integrated record of ShCu macrophysical properties has been developed. The developed record represents the longest available compilation of events with ShCu and includes (i) a novel visualization of the spatial variability in cloud cover both along- and across-wind directions, (ii) updated estimates of narrow-FOV CF and wide-FOV FSC, (iii) updated narrow-FOV CBH, and (iv) complementary data, such as wind speed and direction from the 915-MHz Radar Wind Profiler (RWP) data. The developed record has been used successfully to assess conventional observational estimates of cloud cover and their sensitivity to the following two factors: (i) instrument-dependent cloud detection and data merging criteria and (ii) FOV configuration. Fourth, co-variability of the ShCu macrophysical properties and environmental parameters has been analyzed for a 3-yr period (2016-2018). Our initial analysis includes diurnal changes of FSCs obtained for clouds with small, moderate and large CEDs and several environmental parameters, such as lifted condensation level (LCL) and mixed layer height (zi). Preliminary results of our analysis suggest that the horizontal extent of ShCu is controlled substantially by the sign and magnitude of difference between these two parameters (zi-LCL): the CED tends to grow with increase of this difference (zi exceeds LCL). We have initiated relationships between the ShCu and key atmospheric parameters that control both the development and evolution of ShCu using our new data product, which combines effectively the advantages of narrow-FOV data offered by zenith pointing cloud radars and lidars and wide-FOV TSI images. While the latest instrumentation at the ARM sites may address these challenging relationships in the future, we believe that the historical ARM data at the SGP site has not yet been fully utilized. Overall, our data product can be used by researchers working on a wide range of climate-related projects. These projects may include (i) a comprehensive evaluation of outputs from the Large-Eddy Simulation (LES) and single-column models for their future improvement, (ii) the representativeness of “short-period” results obtained from the previous model and observational studies and (iii) the planning of future field campaigns with focus on improved understanding of the diurnal cycle of cumulus convection.

54 ENVIRONMENTAL SCIENCES↗

Low-Cost Heliostat for High-Flux Small-Area Receivers (Final Technical Report)

This project analyzed a two-stage heliostat concept consisting of a tracking stage and a concentrating stage. The tracking stage uses mirrors mounted on a common drive that move to track the sun. The concentrating stage consists of stationary mirrors that each have a unique angle to direct rays towards a small-area, high-flux, point-focused receiver. By splitting the collection and concentrating process into two stages, multiple small, inexpensive mirrors can share a structure and be controlled by a single drive in the tracking stage. The project effort developed modeling techniques that were specifically relevant to this two-stage heliostat concept. Both field-level and unit-level models were developed. The field-level model does not explicitly consider unit-level losses which are predicted by the unit-level model and then integrated into the field-level model through a correlation referred to as an efficiency modifier. This approach is referred to as the two-model approach; the development and demonstration of this two-model approach for a multi-stage heliostat technology is a key outcome of this work. The field-level model is used to design a field that hits a specific design day power given a set of heliostat design parameters. An oversized field is simulated and then heliostat units are removed based on their annual energy production in order to generate the highest performing field. The field reduction procedure fits a smooth curve fit to annual energy production as a function of position in the field which has the effect of reducing the noise that is otherwise caused by the Monte Carlo ray tracing technique. This approach is referred to as the annual energy fit method and substantially reduces computational run time for a given field level modeling accuracy. The annual energy fit approach enables the selection of a properly sized, high-performing field using orders of magnitude fewer rays than would otherwise be possible and the development of this approach is a second key outcome of this work. These models are used within a genetic optimization algorithm in order to optimize the geometric parameters associated with a heliostat in order to achieve the lowest cost per unit of collected design day power. The cost modeling that underlies the optimization is a simple, scaling type analysis backed up by a much more detailed Design for Manufacture and Assembly (DFMA) analysis. Although the figure of merit used for optimization was not cost per mirror area, this metric is reasonable to use as a means of comparison. The optimally designed 500 kW design has a tracking mirror specific cost of $181.85/m 2 , which is significantly larger than the target value and also larger than the current state of the art. The cost of the torque-tube type linkages contributed substantially to the overall cost. Based on this observation, potentially attractive alternative design configuration utilizing a capstan type actuation system should be investigated. Finally, NREL compared the performance of the two-stage heliostat to the performance of a focused and different sized flat conventional heliostats and showed that, as expected, additional losses versus the convention heliostat caused by a worse cosine efficiency, two stages of reflection, and interstage interactions. The two-stage heliostat requires around 75% more reflective area than a flat 1x1 meter conventional heliostat (similar to a focused heliostat) and 40% more than a flat 2x2 meter conventional heliostat.

14 SOLAR ENERGY↗

Herbaceous Feedstock 2022 State of Technology Report

The U.S. Department of Energy promotes production of advanced liquid transportation fuels from lignocellulosic biomass by funding fundamental and applied research that advances the state of technology (SOT). As part of its involvement in this mission, Idaho National Laboratory completes an annual SOT report for nth-plant and 1st-plant herbaceous biomass feedstock logistics. The purpose of the SOT is to provide the status of feedstock supply system technology development for herbaceous biomass to biofuels relative to technical targets and cost goals from specific design cases, based on data and experimental results. Although conventional feedstock supply systems form the backbone of the emerging biofuels industry, they have limitations that restrict widespread implementation on a national scale. To meet the demands of the future industry, the feedstock supply system must shift from the conventional system to what has been termed “advanced” supply systems. In advanced designs, a distributed network of aggregation and processing centers, termed “depots,” are employed near the points of biomass production (i.e., the field or forest) to reduce feedstock variability and produce feedstocks of a uniform format, moving toward biomass commoditization. The 2022 Herbaceous SOT is part of a vision of achieving an implemented advanced feedstock supply system, which produces a stable, tradable commodity at the decentralized distributed depot. It utilizes feedstock fractionation by incorporating technologies that can separate the biomass into its anatomical fractions (leaves, husks, stems and cobs) to reduce impurities and produce fractions that satisfy downstream quality considerations. By using a series of air classification steps, this strategy can reduce the extrinsic ash in corn stover and produce enriched tissue fractions that can be blended to a conversion specification or converted individually in optimized biochemical conversion campaigns. Additionally, a majority of the leaves (which do not meet the quality specification) are separated out early and can be supplied to alternate markets. The 2022 Herbaceous SOT incorporates an advanced biomass fractionation and processing system to produce pellets enriched tissues from three-pass corn stover. The resulting enriched pellets are delivered to the biorefinery individually where they can be blended to a specification or converted in campaigns where the conditions are optimized for each tissue. Unused fractions can be sent to a a midstream market or to a different conversion process that is better suited to their properties to offset the cost of the delivered feedstock. The main benefits from the proposed system can be summarized as: (1) $6.86/dry ton (2016$) lower cost for the air classification due to elimination of the requirement to discard the high ash lights fraction; (2) $1.56/dry ton lower delivered cost by selling the unsuitable leaf fraction into the feed market as a midstream co-product (assuming a selling price that is 11% higher than their cost of production); (3) 0.98% increase in carbohydrate content (from 60.16% to 61.14%); and (4) 0.97% decrease in ash content (from 6.00% to 5.03%) compared to the 2021 Herbaceous SOT. Overall, the 2022 nth-plant Herbaceous SOT predicts a modeled delivered feedstock cost of $78.64/dry ton (2016$) if it is assumed that the enriched leaf fraction is sold at its production cost; this is a slight increase of $0.43/dry ton increase from the 2021 Herbaceous SOT nth-Supply case cost. The increased cost derived from a $0.38/dry ton increase in transportation and handling cost to procure more biomass (to replace the enriched leaf fraction that was not delivered to the biorefinery. The total preprocessing cost was $0.27/dry ton higher than the 2021 result because of updates to energy consumption, purchasing price and dry matter loss data for the rotary shear ($3.00/dry ton increase) and the pelleting mill ($4.52/dry ton increase). The data utilized were generated in pilot-scale tests in the Biomass Feedstock National User Facility (BFNUF) at INL and at Forest Concepts, including tests for rotary shear and pelleting of the air classified fractions. A greenhouse gas emissions analysis was performed by Argonne National Laboratory using the most up to date version of the Greenhouse Gases, Regulated Emissions, and Energy use in Transportation model (GREET®). The analysis showed an increase of 17.34 kg CO2e/dry ton from the 2021 SOT (67.71 kg CO2e/ton in the 2021 Herbaceous SOT to 85.05 kg CO2e/ton in the 2022 Herbaceous SOT). The net increase is primarily attributed to increased energy consumption in pelleting mill.

09 BIOMASS FUELS↗

Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (Final Report)

This document presents final technical findings for the project Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (FE0031782). The project is part of the U.S. Department of Energy Oil and Gas Program to develop and advance technologies that can significantly improve the recovery efficiencies of unconventional oil and gas resources. The overall objective of this project is to improve the ultimate recovery from unconventional oil and gas (UOG) resources in the United States by developing a method for using unrefined natural gas liquids (NGLs) as treatment fluids to improve hydrocarbon production. Horizontal extended-lateral drilling coupled with high volume hydraulic fracturing has significantly increased production from UOG resources in the U.S. However, the recovery efficiency is low compared to the estimated oil and gas in place. Recent data indicate that less than 10% of the oil in the liquid-rich UOG reservoirs is produced. Alternative completion methods using NGLs could increase production (Battelle, 2016; Wan et al., 2013; Wan, 2013; Downey et al, 2021); however, field validation tests are needed to develop an approach that is economical, efficient, and compatible in the UOG setting to advance towards commercial deployment. This project aims to develop and field test a method to improve recovery of oil resources in UOG shale plays by using Y-Grade NGLs, or a similar combination of NGLs, as treatment fluids. Refined NGLs have been used as a hydraulic treatment fluid in UOG plays for decades and are shown to be particularly effective because their miscibility with oil allows oil to flow more freely; however, the use of Y-Grade (unrefined) NGLs has not been studied. The use of Y-Grade NGLs would be advantageous over refined NGLs because Y-Grade NGLs do not require infrastructure or investment in refining and are already being produced from many UOG reservoirs. The concept was tested and monitored in the field at a commercial well site owned by project partner Hopco, Ltd. The project team, which consists of multiple oil and gas operators, Linde Gas North America LLC (Linde) and the Ohio Division of Geological Survey (ODGS), has extensive experience with oil and gas production in the Appalachian basin and the ability to work together quickly to solve technical issues and research needs. A key part of the proposed work was the use of existing wells for field testing and monitoring. A total of four wells (three vertical and one horizontal) were available for this project. One of the vertical wells was utilized as the test well for the NGL treatment test. A nearby vertical well was used for microseismic monitoring. The remaining vertical well and the horizontal well provided a baseline for typical UOG production in the oil window of the Utica/Point Pleasant (UPP). Major technical tasks of the project include characterization of the geotechnical properties of the UPP with an emphasis on the field site; design and planning for the NGLs testing; field testing and monitoring; analysis and integration of field data; and economic and resource/reserve assessment. The Shoman monitoring well and Doughty NGL treatment well were successfully plugged back during September-October 2020 in preparation for treatment and monitoring. A nitrogen diagnostic fracture injection test (DFIT) was completed on the Doughty well on July 22, 2021, consisting of 133,000 scf (91 Bbl.) of nitrogen. A nitrogen foam frac was completed in the Utica-Point Pleasant interval on August 17, 2021 with funding from outside sources. A microseismic monitoring array was installed in the Shoman well and monitored microseismic activity during the Doughty well frac job. Y-Grade NGL injection commenced on August 26, 2021. A total of 215 Bbl. was injected but the job was shut down due to a small leak on the suction hose on the pump truck. The Y-Grade treatment resumed on August 27, 2021 and an additional 726 Bbl. of Y-Grade NGL was injected at a well head pressure of 3850 psi. Total volume of injected Y-Grade over the two days of injection was 941 Bbl. The well was shut-in for 17 days following injection to allow the Y-Grade NGLs to soak on the formation. Y-Grade treatment flow back commenced on 9/13/2021 on a weekly basis. Production data, including surface pressures, oil, nitrogen, natural gas, and flow times was measured and recorded. Periodic gas samples were collected and analyzed to determine composition of flowback gas. As of July 2022, the treatment well had produced 726 Bbl of oil and 2,888 mcf gas. In August 2022, tubing and packer in the well were removed and a pump was installed to enhance oil recovery. Currently, the operator is producing the well about 2 days a week for a few hours. The performance of the NGL treatment test was evaluated based on reservoir simulations of the treatment process, processing of well testing data, analysis of micro seismic monitoring data, and production data analysis. This analysis suggested that oil production in the small test would continue through 2025. Upscaling the treatment to a horizontal Utica Point-Pleasant well would allow more oil recovery, but the process would involve more investment, services, and operational support. An economic analysis was conducted for scenarios aimed at upscaling the NGL treatment process for more typical horizontal UPP wells in the Appalachian Basin.

02 PETROLEUM↗

Enhancement of Optical Efficiency of CSP Mirrors for Reducing O&M Cost via Near-Continuous Operation of Self-Cleaning Electrodynamic Screens (EDS). Final Report

Over the past decade, techno-economical advancements in solar energy systems, particularly the improvement in the conversion efficiency of PV modules made with mono-crystalline silicon solar cells from 12% to 20%, as well as the cost reduction in manufacturing by a factor of about 10%, have made it possible to achieve a levelized cost of electricity (LCOE) in PV plants that is comparable to, or less than, the cost of deriving electricity from fossil fuels. Operating PV plants in mid-latitude sun-belt regions, where the solar irradiance level is highest, provides a high annual energy-yield (kWh/kWp) due to two factors: (1) the availability of predictable high solar irradiance throughout the year with the fewest interruptions in solar flux from clouds and rain, and (2) the increased conversion efficiency of crystalline solar cells, as recombination loss has decreased with increased intensity of the sunlight that illuminates the silicon solar cells. Semi-arid and desert regions, however, are plagued by high atmospheric dust concentrations and frequent sand storms. The deposition of a layer of dust on the optical surfaces of solar collectors such as PV modules and concentrating mirrors reduces the transmission efficiency of sunlight that actually reaches the solar cells or receivers, resulting in high energy-yield soiling loss. There are two major cost components to operating a solar plant: (1) installation costs, and (2) operation-and-maintenance (O&M) costs. There is no fuel cost; hence operating a solar plant in a semi-arid or desert region provides high returns on investment if soiling losses are mitigated via efficient cleaning methods and optimized cleaning frequency. If solar collectors are not cleaned, the accumulation of dust layers on solar collectors may cause the operation of such plants in arid regions to become economically unviable. Washing solar collectors with water and detergent, as is most commonly done now, is an efficient method for cleaning. The conventional approach in utility-scale solar plants is to use a large truck with a water tank and pump system for spraying deionized water on the surface of the solar collectors. Robotic cleaning with brushes, used for many solar plants, requires lesser water for cleaning. The water consumed using semi-automated cleaning of PV modules in utility-scale solar plants is approximately 2 liters/m 2 per cleaning cycle. The total optical surface area of the solar collectors in 1 TW-scale solar installation will be more than 3 × 109 m 2 ; hence an enormous amount of water be needed for cleaning. There simply is not enough fresh water in the sun-belt areas of the world for predicted cleaning needs. In solar power plants, the estimated cost of cleaning solar collectors includes expenses related to the equipment used, labor, cost of transportation of water, the energy required for cleaning, plus ancillary costs whereas the cost of water is not considered. The water used is obtained from sources located close to plant sites, unmindful of the environmental and societal impact as the power plants are oftentimes located in regions that face severe drought. This practice is very similar to the cost calculations in deriving the levelized cost of electricity (LCOE) in conventional power plants based on burning fossil fuel such as coal and gas, while disregarding the cost of climate change and health effects. Unless a water-free or low-water cleaning method is established, the expansion of solar plants may lose public support in areas suffering long intervals of drought. The goal of this research project has been the development and application of the Electrodynamic Screen (EDS) as a means for a water-free, scalable cleaning process applicable to solar-power installations, including rooftop applications. We describe here the development of an EDS film-based cleaning process as an emerging method for use on PV modules, parabolic troughs, and heliostats. This report aims to show the feasibility of integrating or retrofitting EDS films onto the optical surfaces of solar collectors (both PV and CSP) while maintaining high transmission or reflection efficiency. The cleaning action provided by the EDS film is an active method to remove dust deposits by electrodynamic force. Current lab-scale prototype EDS films, retrofitted onto solar panels and mirrors, have shown to be capable of maintaining optical transmission or specular-reflection efficiencies higher than 90% of initial values under clean conditions. The optical surfaces of solar collectors laminated with EDS films can remove more than 90% of deposited dust when the EDS is activated for less than two minutes. As an electrodynamic dust removal process, the EDS film-based method is designed primarily for the removal of dust in solar installations located in semi-arid and desert areas, where the atmosphere is often dry and dusty and rainfall is infrequent. While EDS film application minimizes water consumption and facilitates cleaning as frequently as needed, it has limitations in removing contaminants such as soot, organic pollutants deposited as fine films on the surface, and bird droppings. The dust removal efficiency of the EDS is maximum at relative humidity RH is 40 to 50% and decreases at when RH > 65%. Many solar plant sites undergo diurnal and seasonal cycles of high ambient, early afternoon temperatures and an RH that reaches the dew point early in the morning. These variations in atmospheric conditions do not limit the operation of the EDS film. This report presents a brief review of the progress and the potential of EDS film technology for mitigating the impact of dust on solar collectors via water-free cleaning, as well as current technical challenges regarding efficiency and durability. Our experimental data on the performance of EDS films show that: (1) the dust-removal efficiency (DRE) can reach levels higher than 90%, (2) the specular reflectivity (SR) of EDS film-laminated second-surface mirrors reach levels in excess of 90%, (3) the specular reflectivity restoration (SRR) can exceed 90%, (4) the output-power restoration (OPR) of PV modules can exceed 95%, and (5) the optical transmission efficiency (TE) of the EDS films can be greater than 90%. Working with Sandia National Laboratories (NM), Corning Research and Development Corporation (NY), Eastman Kodak (NY), Tomark-Worthen Industries (NH), and EDS Chile SPA (Chile), we have produced EDS film-laminated PV modules and demonstrated their self-cleaning functions without requiring water. We have demonstrated that the operational range of EDS films will cover the expected ambient temperature of solar fields at RH cycling varying from 20 to 95% as long as the EDS films are activated in the RH range 20 to 50%. (Typical solar-field climates in deserts and semi-arid lands often reach near dew point in coastal areas.) Our experiments on the application of hydrophobic-fluorinated nanoparticle coatings on EDS film surfaces show that the EDS operational range can be extended to higher RH levels that approach the dew point. At Eastman Kodak, as one of our industrial partners, we were able to establish a process for manufacturing EDS films using flexographic printing of the electrodes onto transparent polymer films. This process utilizes an existing manufacturing line at Eastman Kodak that allows fabrication of medium-scale EDS films (26 cm × 30 cm). The manufacturing process has the capacity to produce EDS films at high production speeds. The medium-scale EDS films that have been printed at Kodak were evaluated in the lab at Boston University. The EDS films produced at Kodak were laminated at Tomark-Worthen using an industrial scale vacuum laminator to produce EDS film stacks, which can be affixed onto the optical surfaces of PV modules or concentrating mirrors. The EDS film stack consists of the EDS films that have Willow® Glass (WG) which has a thickness of 100 μm as the front surface. The WG sheets obtained from Corning Research and Development Corporation are customized in size and shape to cover the active area of the EDS films. The back surface of the EDS film stack is integrated onto the optical surface of the solar panel or mirror using optically clear adhesive (OCA) films or silicone adhesives. The OCA films (thickness 25 μm) are produced by 3M. The EDS film stacks have the architecture: WG/OCA/EDS Film/OCA/ over PV module or solar mirror. The power-supply units needed for activating the electrodes of the EDS film were designed and produced at Boston University. These power supply units provide three-phase, 1.2 kV voltage pulses at a very low current (micro-ampere) level and at a low frequency (≈ 5 Hz). The voltage pulses are applied to the electrodes in a sequence such that the train of pulses resembles a unidirectional traveling wave of electrical field on the surface of the EDS film. The dust particles on the surface become charged electrostatically and are levitated by the Coulomb force. The lateral sweeping action of the traveling electric field created by the three-phase voltages pulses then sweeps the dust off the surface. The energy consumed by the EDS electrodes is less than 0.2 Wh/m 2 /cleaning cycle, enabling energy-efficient restoration of output power (OPR) of PV module or specular reflectivity restoration (SRR) for solar mirrors. The EDS system consists of (1) an EDS film stack laminated onto the solar collectors, (2) connection of the EDS film to its power supply unit and (3) Interconnection of the power supply to PV modules or solar mirrors. Design, construction and assessment of field-testing units that have EDS stack laminated PV modules for evaluating the performance of the EDS films in solar fields is being carried out at BU. Our progress under this project, aimed at the advancement of EDS film technology, has reached DOE Technology Readiness Level (TRL) 6. Based on the extensive laboratory evaluations and limited field trials, as well as contacts with potential users, we believe that the technology has reached its commercial stage. We are conducting a cost analysis using the National Renewable Energy Laboratory (NREL) System Advisory Model (SAM) and are preparing for field trials of EDS films in different solar fields in the US, Chile, India and in the Middle East. A brief description of EDS film performance, construction and testing of the field-test unit, autonomous operation of the field-test unit for evaluating EDS performance in increasing energy yield, associated revenue savings, and water conservation is presented.

14 SOLAR ENERGY↗