Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distributed Energy Program”

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 1,189 records · Page 66

A Life Cycle Analysis of Refinery Fuel Products Air Emissions Using Refinery Emissions Inventory Data

With the transportation sector being a major source of air pollutant emissions, this study seeks to evaluate criteria pollutant emissions associated with the use of petroleum fuels in vehicles using the GREET® (Greenhouse gases, Regulated Emissions, and Energy use in Transportation) well-to-wheel model. The life cycle of petroleum fuels includes the major stages of crude recovery, crude transportation, crude refining, fuel distribution and fuel combustion (vehicle operation). This study focuses on quantifying refinery emissions, in terms of both refinery throughput and various refinery fuel products, by matching a refinery emissions inventory dataset with refinery operation information, at various levels (national, regional, facility, and process unit). The refinery emissions in terms refinery throughput (per refinery crude input, in g/mmbtu crude) are developed at a regional level for various Petroleum Administration for Defense Districts (PADD), and also aggregated at the national level. To allocate refinery emissions to refinery products, a methodology is developed using a refinery operation optimization Linear Programming (LP) model for individual refineries. The results from individual refinery facilities are then aggregated to obtain statistical national average (capacity weighted) emissions for each product. Because refineries vary by complexity of operation, configuration, crude quality, and product slate, the derived emissions demonstrate large variation for each product and per unit of crude input. Relative to previous estimates used in the GREET model, the present research shows reductions in refinery (and life cycle) emissions for most criteria pollutants. The reduction is especially significant for sulfur dioxide (SOx) emissions.

Allocation↗

Macrophysical properties of continental cumulus clouds from active and passive remote sensing

Cloud amount is an essential and extensively used macrophysical parameter of cumulus clouds. It is commonly defined as a cloud fraction (CF) from zenith-pointing ground-based active and passive remote sensing. However, conventional retrievals of CF from the remote sensing data with very narrow field-of-view (FOV) may not be representative of the surrounding area. Here we assess its representativeness using an integrated dataset collected at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site in Oklahoma, USA. For our assessment with focus on selected days with single-layer cumulus clouds (2005-2016), we include the narrow-FOV ARM Active Remotely Sensed Clouds Locations (ARSCL) and large-FOV Total Sky Imager (TSI) cloud products, the 915-MHz Radar Wind Profiler (RWP) measurements of wind speed and direction, and also high-resolution satellite images from Landsat and the Moderate Resolution Imaging Spectroradiometer (MODIS). We demonstrate that a root-mean-square difference (RMSD) between the 15-min averaged ARSCL cloud fraction (CF) and the 15-min averaged TSI fractional sky cover (FSC) is large (up to 0.3). We also discuss how the horizontal distribution of clouds can modify the obtained large RMSD using a new uniformity metric. The latter utilizes the spatial distribution of the FSC over the 100° FOV TSI images obtained with high temporal resolution (30 sec sampling). We demonstrate that cases with more uniform spatial distribution of FSC show better agreement between the narrow-FOV CF and large-FOV FSC, reducing the RMSD by up to a factor of 2.

ARM Active Remotely Sensed Cloud Locations (ARSCL)↗

Exploring baryon resonances with transition generalized parton distributions: status and perspectives

QCD gives rise to a rich spectrum of excited baryon states. Understanding their internal structure is important for many areas of nuclear physics, such as nuclear forces, dense matter, and neutrino-nucleus interactions. Generalized parton distributions (GPDs) are an established tool for characterizing the QCD structure of the ground-state nucleon. They are used to create 3D tomographic images of the quark/gluon structure and quantify the mechanical properties such as the distribution of mass, angular momentum, and forces in the system. Transition GPDs extend these concepts to N → N* transitions and can be used to characterize the 3D structure and mechanical properties of baryon resonances. They can be probed in high-momentum-transfer exclusive electroproduction processes with resonance transitions e + N → e' + M + N*, such as deeply-virtual Compton scattering (M = γ) or meson production (M = π, K, etc.), and in related photon/hadron-induced processes. This White Paper describes a research program aiming to explore baryon resonance structure with transition GPDs. This includes the properties and interpretation of the transition GPDs, theoretical methods for structures and processes, first experimental results from JLab 12 GeV, future measurements with existing and planned facilities (JLab detector and energy upgrades, COMPASS/AMBER, EIC, EicC, J-PARC, LHC ultraperipheral collisions), and the theoretical and experimental developments needed to realize this program.

Experimental Nuclear Physics↗

Dynamic Probabilistic Risk Assessment Based Response Surface Approach for FLEX and Accident Tolerant Fuels for Medium Break LOCA Spectrum

After the Fukushima Daiichi Accident, the safety features such as accident tolerant fuel (ATF) and diverse and flexible coping strategies (FLEX) for existing nuclear fleets are being investigated by the US Department of Energy under the Light Water Reactor Sustainability Program. This research is being conducted to quantify the risk-benefit of these safety features. Dynamic probabilistic risk assessment (DPRA)-based response-surface approach has been presented to quantify the FLEX and ATF benefits by estimating the risk associated with each option. ATFs with multilayered silicon carbide (SiC), iron-chromium-aluminum, and chromium-coated zirconium cladding were considered in this study. While these ATF candidates perform better than the current zirconium cladding (Zr), they may introduce additional failure modes in some operating conditions. The fuel failure analysis modules (FAMs) were developed to investigate ATF performance. The dynamic risk assessments were performed using RAVEN, a DPRA tool, coupled with RELAP5 and FAMs. A cumulative distribution function-based index provided a mean of comparing the benefits of safety enhancements. For medium break loss of coolant accidents, FLEX operational timing window for each fuel type was estimated. Among these ATF candidates, SiC-type ATF was the most beneficial candidate for an increased safety margin than Zr-based fuel and was found to complement FLEX strategies in terms of risk and coping time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Unusual Multiwavelength Properties of the Gamma-Ray Source PMN J1603-4904

Context. We investigate the nature and classification of PMNJ1603-4904, a bright radio source close to the Galactic plane, which is associated with one of the brightest hard-spectrum gamma-ray sources detected by Fermi/LAT. It has previously been classified as a low peaked BL Lac object based on its broadband emission and the absence of optical emission lines. Optical measurements, however, suffer strongly from extinction and the absence of pronounced short-time gamma-ray variability over years of monitoring is unusual for a blazar. Aims. In this paper, we are combining new and archival multi-wavelength data of PMNJ1603-4904 in order to reconsider the classification and nature of this unusual gamma-ray source. Methods. For the first time, we study the radio morphology of PMNJ1603-4904 at 8.4GHz and 22.3GHz, and its spectral properties on milliarcsecond scales, based on VLBI observations from the TANAMI program. We combine the resulting images with multi-wavelength data in the radio, IR, optical/UV, X-ray, and gamma-ray regimes. Results. PMNJ1603-4904 shows a symmetric brightness distribution at 8.4 GHz on milliarcsecond scales, with the brightest, and most compact component in the center of the emission region. The morphology is reminiscent of a compact symmetric object (CSO). Such objects, thought to be young radio galaxies, have been predicted to produce gamma-ray emission but have not been detected as a class by the Fermi gamma-ray telescope so far. Sparse (u, v)-coverage at 22.3GHz prevents an unambiguous modeling of the source morphology at this higher frequency. Moreover, infrared measurements reveal an excess in the spectral energy distribution (SED), which can be modeled with a blackbody with a temperature of about 1600 K, and which is usually not present in blazar SEDs. Conclusions. The TANAMI VLBI data and the shape of the broadband SED challenge the current blazar classification of one of the brightest gamma-ray sources in the sky. PMNJ1603-4904 seems to be either a highly peculiar BL Lac object or a misaligned jet source. In the latter case, the intriguing VLBI structure opens room for a possible classification of PMNJ1603-4904 as a gamma-ray bright CSO.

Multiwavelength↗

Global transpiration data from sap flow measurements: the SAPFLUXNET database

Plant transpiration links physiological responses of vegetation to water supply and demand with hydrological, energy, and carbon budgets at the land–atmosphere interface. However, despite being the main land evaporative flux at the global scale, transpiration and its response to environmental drivers are currently not well constrained by observations. Here we introduce the first global compilation of whole-plant transpiration data from sap flow measurements (SAPFLUXNET, https://sapfluxnet.creaf.cat/, last access: 8 June 2021). We harmonized and quality-controlled individual datasets supplied by contributors worldwide in a semi-automatic data workflow implemented in the R programming language. Datasets include sub-daily time series of sap flow and hydrometeorological drivers for one or more growing seasons, as well as metadata on the stand characteristics, plant attributes, and technical details of the measurements. SAPFLUXNET contains 202 globally distributed datasets with sap flow time series for 2714 plants, mostly trees, of 174 species. SAPFLUXNET has a broad bioclimatic coverage, with woodland/shrubland and temperate forest biomes especially well represented (80 % of the datasets). The measurements cover a wide variety of stand structural characteristics and plant sizes. The datasets encompass the period between 1995 and 2018, with 50 % of the datasets being at least 3 years long. Accompanying radiation and vapour pressure deficit data are available for most of the datasets, while on-site soil water content is available for 56 % of the datasets. Many datasets contain data for species that make up 90 % or more of the total stand basal area, allowing the estimation of stand transpiration in diverse ecological settings. SAPFLUXNET adds to existing plant trait datasets, ecosystem flux networks, and remote sensing products to help increase our understanding of plant water use, plant responses to drought, and ecohydrological processes. SAPFLUXNET version 0.1.5 is freely available from the Zenodo repository (https://doi.org/10.5281/zenodo.3971689; Poyatos et al., 2020a). The “sapfluxnetr” R package – designed to access, visualize, and process SAPFLUXNET data – is available from CRAN.

58 GEOSCIENCES↗

SelenITA: A Dual Point Lunar Mission to Characterize the Near Surface Dust and Electromagnetic Plasma Environment

SelenITA is a dual CubeSat mission that will provide the first multi-point dust, plasma, and magnetic field measurements in lunar orbit. This mission will advance the understanding of the electromagnetic space environment at the Moon in support of the Artemis program, exploration, and the geosciences. Here we present the science rationale motivating the mission. The candidate mission science objectives include: (1) constrain the origins of crustal magnetic fields; (2) determine the nature of plasma interactions with crustal magnetic fields; (3) characterize plasma waves and turbulence at the Moon; (4) characterize the lunar surface electric potential in varying plasma environments; (5) constrain the composition, thermal state, and structure of the lunar upper mantle and crust; (6) characterize the ionizing radiation in lunar orbit; and lastly, (7) determine the density of the dust exosphere as a function of latitude, longitude, and altitude, including the lunar polar space environment. The measurement requirements include simultaneous two-point observations of the 3-component vector magnetic field, plasma distribution (flux, energy, density, temperature), and single-point observations of energetic particles (protons, electrons, gamma rays), and dust. These measurements are important because it helps us understand how future astronauts, robots, and space hardware will live and work on the lunar surface.

space weather↗

SelenITA: A Dual Point Lunar Mission to Characterize the Near Surface Dust and Electromagnetic Plasma Environment

SelenITA is a dual CubeSat mission that will provide the first multi-point dust, plasma, and magnetic field measurements in lunar orbit. This mission will advance the understanding of the electromagnetic space environment at the Moon in support of the Artemis program, exploration, and the geosciences. Here we present the science rationale motivating the mission. The candidate mission science objectives include: (1) constrain the origins of crustal magnetic fields; (2) determine the nature of plasma interactions with crustal magnetic fields; (3) characterize plasma waves and turbulence at the Moon; (4) characterize the lunar surface electric potential in varying plasma environments; (5) constrain the composition, thermal state, and structure of the lunar upper mantle and crust; (6) characterize the ionizing radiation in lunar orbit; and lastly, (7) determine the density of the dust exosphere as a function of latitude, longitude, and altitude, including the lunar polar space environment. The measurement requirements include simultaneous two-point observations of the 3-component vector magnetic field, plasma distribution (flux, energy, density, temperature), and single-point observations of energetic particles (protons, electrons, gamma rays), and dust. These measurements are important because it helps us understand how future astronauts, robots, and space hardware will live and work on the lunar surface.

space weather↗

An approximate solution for interlaminar stresses in laminated composites: Applied mechanics program

An approximate solution for interlaminar stresses in finite width, laminated composites subjected to uniform extensional, and bending loads is presented. The solution is based upon the principle of minimum complementary energy and an assumed, statically admissible stress state, derived by considering local material mismatch effects and global equilibrium requirements. The stresses in each layer are approximated by polynomial functions of the thickness coordinate, multiplied by combinations of exponential functions of the in-plane coordinate, expressed in terms of fourteen unknown decay parameters. Imposing the stationary condition of the laminate complementary energy with respect to the unknown variables yields a system of fourteen non-linear algebraic equations for the parameters. Newton's method is implemented to solve this system. Once the parameters are known, the stresses can be easily determined at any point in the laminate. Results are presented for through-thickness and interlaminar stress distributions for angle-ply, cross-ply (symmetric and unsymmetric laminates), and quasi-isotropic laminates subjected to uniform extension and bending. It is shown that the solution compares well with existing finite element solutions and represents an improved approximate solution for interlaminar stresses, primarily at interfaces where global equilibrium is satisfied by the in-plane stresses, but large local mismatch in properties requires the presence of interlaminar stresses.

Rose, Cheryl A.↗

Optimization of Energy Storage System Economics and Controls by Incorporating Battery Degradation Costs in REopt

The use of stationary electrochemical energy storage systems utilizing lithium-ion batteries has increased rapidly as the production scale and price for lithium-ion batteries has decreased. These energy storage systems are crucial for maintaining grid resiliency, especially for grids operating with high penetration of renewable energy generation assets or for with a variety of distributed energy generation and storage systems. One challenging factor for the development of battery energy storage systems is estimating the proper sizing, in terms of both power and energy, that minimizes total costs over the lifetime of the systems; this calculation is difficult in simple cases, where a battery is costed independently, but is extremely challenging when building loads and electrical generation by photovoltaic resources are also considered. REopt is a techoeconomic optimization tool developed by NREL to address these challenges. Previously, battery degradation has been priced by simply assuming a 10-year replacement schedule for battery systems. However, this does not account for varying degradation trends observed across real-world batteries, or allow for batteries to be operated in a degradation-aware manner that optimizes battery dispatch based on operating costs. This work incorporates a battery life model into REopt. This battery life model is simple, so that it may be solvable within the constrains of a mixed-integer linear optimization problem, but is fit to accelerated aging data recorded in the lab. To achieve the best possible accuracy for lifetime estimates given these constraints, parameters for the battery life model in REopt are estimated by fitting 20-year simulations of battery life after identifying state-space battery degradation model from accelerated aging data. Comparisons of battery life predicted in REopt and from the state-space battery degradation model to ensure validity of lifetime estimates made by REopt. Battery life and cost is optimized by controlling three decision to minimize system life cost: battery sizing, daily state-of-charge, and daily energy-throughput. The cost of battery degradation as a function of these control variables is then estimated assuming two possible maintenance strategies: replacement, where the entire battery system is replaced if cell reach an end-of-life capacity threshold; and augmentation, which establishes a fund to pay for continual purchase of new batteries to maintain the initial energy capacity of the system. These two strategies offer conservative (for replacement) and optimistic (for augmentation) bounds for total system cost. The degradation cost incurred by these strategies is then used to control battery dispatch decisions, operating the battery in a degradation-aware manner that maximizes battery lifetime while also providing energy when favorable. Because the mixed-integer linear program has perfect foresight of future energy needs, batteries with degradation costs are always operated using 'just-in-time' charging, which is unrealistic, as no energy is left in the storage system to perform other energy services or to serve as emergency back-up power. To combat this, an inequality constraint on the average annual state-of-charge is imposed, and the sensitivity of system cost to average stored energy, e.g., the cost of system resiliency, can be quantified. Analysis of results has several conclusions, for instance, oversizing of battery storage systems is not a cost burden when battery storage is an optimal solution, as any additional battery capacity can simply be utilized to avoid costs of purchasing energy from a utility.

battery↗

Key insights from US Department of Energy Better Plants workforce development bootcamps (2022–2025)

This study examines the effectiveness of the US Department of Energy’s Better Plants Program Bootcamps, which are designed to enhance participants’ technical skills in improving energy efficiency and optimizing operations in manufacturing facilities. Through the analysis of survey data collected from 529 participants across 9 bootcamps, the research investigates the motivations, benefits, and demographic trends of attendees. The findings reveal that skill acquisition and improvement are primary drivers for participation, with key benefits including hands-on training on diagnostic equipment and software tools, networking opportunities, and access to technical resources. The analysis shows strong participation from sectors characterized by high energy consumption and employment, such as chemical and transportation equipment manufacturing. Over 50% of participants have job titles that include “EHS” or “Energy” showing their key roles in leading energy efficiency and energy management efforts in manufacturing. Furthermore, the analysis highlights the distribution of participants across managerial, engineering, and technical roles, revealing a higher representation of managers and engineers. This observation suggests a need for targeted outreach to engage technicians, equipment operators, maintenance staff, and floor workers to ensure comprehensive workforce development. The post-bootcamp survey showed that the participants highly valued the opportunities for peer learning and idea exchange, and the benefits they gained from them. This research contributes to the advancement of manufacturing education by demonstrating the efficacy of specialized training in addressing critical industry challenges and fostering a more competent and empowered workforce.

Energy efficiency↗

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Constellation-X Observatory

The Constellation X-ray mission will revolutionize X-ray spectroscopy and will provide a key component of NASA's '!Beyond Einstein" program. Over the past 35 years, X-ray astronomy has grown from infancy to near-maturity and is now an indispensable astrophysical tool. Constellation-X will transform this field from the exploratory spectroscopic era of Chandra and XMM to an era of fully comprehensive spectroscopic studies. The mission will consist of four 1.6-meter X-ray telescopes orbiting the Earth/Sun system. By providing a factor of 100 increase in throughput over previous missions, as well as high spectral resolution and a broad energy bandpass, Constellation-X will obtain high quality spectra for all classes of X-ray sources over a wide range of luminosity and redshift. It will study the warping of space and time near black holes, determine the distribution of ordinary matter, dark matter and dark energy, and probe detailed physical processes that occur at temperatures, densities and pressures far beyond those achievable in Earth-bound laboratories.

Weaver, Kimberly↗

Physics based Reliability Models for Supercritical CO 2 Turbomachinery Components

CSP plants using sCO 2 power cycle can potentially achieve high thermal conversion efficiency at low capital cost due to compact turbo-machinery and other components. An sCO 2 expander and improved heat exchanger is expected to provide a major stepping stone for achieving CSP power at $0.06/kW-hr LCOE, energy conversion efficiency > 50%, and total power block cost < $1,200/kW installed. Hybrid gas bearings (HGB) and Dry gas seals (DGS) are two key components that are critical for robust and efficient operation of future sCO 2 power cycles. From a system perspective, HGBs can provide substantial benefits (~3 to 5%) to the modular CSP operation, while use of DGSs limits the leakage losses as well as the windage losses in the generator and can increase overall generator efficiency by ~ 8%. The life limiting mechanisms of HGBs and DGSs in high pressure, high temperature sCO2 environment were not well understood and hence this program aims at developing physics-based performance prediction and life prediction models for these two components. The output of the performance models serves as input to the life models. The performance models were validated using experimental data and using data from literature. DGS and HGB models have shed light on some risks associated with its operation in sCO2 expander and will guide more robust designs tailored towards sCO 2 cycle applications. A DGS was tested in 1MWe sCO2 turbo expander to study its performance and validate the performance models. To understand the effect of high pressures, high temperatures and sCO 2 environment on crack initiation, crack propagation and LCF/HCF life of these components, a novel experimental setup was developed. The LCF test data was used to calibrate and validate the life prediction models. The effect of sCO 2 corrosion on the life of Ni-based super-alloys is quantified using long duration experiments. Advanced microstructure and spectroscopic analyses conducted shed light on some key differences between IN718 and IN617 in terms of oxidation morphology, chemical species diffusion and trapping, the formation of protective corrosion resistant layers and changes in surface properties. A novel crack initiation model was developed which takes into account strain distribution, temperature, hold times and activation energy at the gas-metal interface and predicts the number of cycles to initiate the crack in the Inconel alloy. An advanced 3D fracture mechanics-based crack propagation model was also developed which takes into account the transient loading mission of the component and predict how fast the crack will grow and what will be its final shape and size. High energy X-Ray tomography measurements were conducted to study and validate 3D morphology of crack growth inside the specimen. These two models along with chemical kinetics models obtained from corrosion tests can be combined to predict the total life of the component under given long term operating conditions. A Bayesian Hybrid Modeling (BHM) probabilistic framework was developed to quantify uncertainty in physics-based performance models and in crack initiation and crack propagation models, to calibrate the models and to validate the models with statistical confidence. Advanced coupled physics multi-phase flow models were developed to study condensation, boiling and phase change behavior of CO 2 under various sCO 2 turbomachinery and heat exchanger conditions. A novel experimental setup was also designed to observe phase change of CO 2 during its transition across liquid-gas-supercritical state lines.

14 SOLAR ENERGY↗

High Energy Storage Capacity Low Cost Iron Flow Battery

A new flow battery was proposed that utilizes low cost materials: iron as the only active element, cheap aqueous electrolytes, and inexpensive separators. During charging, ferrous iron (Fe 2+ ) is oxidized to ferric iron (Fe 3+ ) at the positive electrode while it (Fe 2+ ) is reduced to form iron metal (Fe 0 ) at the negative electrode. Because iron is plated at the negative electrode during charging, conventional electrode structures couple the energy storage capacity and the power rating of the battery. In order to decouple the energy and power ratings and regain the economic advantages of a flow battery, a slurry electrode design was proposed for use as the negative electrode. The slurry electrode is made by flowing electrically conductive particles in an electrolyte containing the dissolved iron species. On charging, iron is plated onto the particles. The particles can then carry the iron metal out of the cell to be stored in external reservoirs. This project looked to develop a slurry for the negative electrode of the all-iron flow battery that can operate at 200 mA/cm 2 with <100 mV overpotential, be pumped with parasitic energy costs < 10%, and promote plating on the particles. A scalable and stackable cell design of 1100 cm 2 was developed. The project was initially under contract in mid-January 2013. The project was converted from a one-year proof of concept seedling to an additional two-year program, then extended with an 18-month plus-up, and finally with an additional 6-month approved cost-overrun. During the seeding project, the concept of plating on carbon particles in a slurry was demonstrated. However, this was accomplished with high-cost multi-walled carbon nanotubes as the slurry material. The following two-year program discovered several low-cost carbon blacks that could also be used in this application. These findings were supported by an analysis that demonstrated the effectiveness of the carbon that takes into account not just the effective electronic conductivity of the slurry, but also the surface area per volume of the slurry. It was demonstrated that stable, flowable slurries could be created with these carbons, even when the carbons where fully loaded with iron during charging. During the two-year program, laboratory-scale 50 cm 2 single-cells demonstrated charge-discharge cycling at current densities of 150 mA/cm 2 and energy efficiencies of >50%. Meanwhile, cost modeling continued to show promise for this energy storage approach. Following the conversion period, the overall objective of the nominally 18-month plus up period was to close the technical gap of a 25X scale-up and reduce risk to enable commercial investment required to develop a slurry iron flow battery product. A specific objective was to demonstrate 150 cycles of performance exceeding 70% energy efficiency with a 10-cell stack of 1156 cm 2 cells operating at 100 mA/cm 2 . The 10-cell short stack was designed to have a power delivery rating of approximately 1 kW, which represents a size commensurate with the requirements of our commercialization partner at that time. A complete battery support system for the full-size cell was designed, fabricated, and tested. This includes pumps (multiple types were considered, and constant volume hose pumps were found to be acceptable for performance and cost), electrolyte rebalance reactor, and sensors and control systems. The support system proved to be well-behaved and worked well. A new low-cost membrane was developed for this system based on a microporous polymer support with a layer of hydrogel polymer. The membrane had low resistivity, and significantly reduced pressure driven fluid flow across the membrane (as compared to the un-coated microporous support). This membrane approach was scaled-up in conjunction with a commercial toll coating company, and a large amount of membrane was manufactured to support large cell and stack testing. A new concept for an electrolyte rebalance reactor to react hydrogen gas from the negative electrolyte (from plating inefficiency) with excess ferric ion in the positive electrolyte was developed. This design was very simple, low cost, and effective. The full scale 1156 cm 2 single-cell was designed, modeled, fabricated and tested. Modeling showed that uniform flow-distribution should be obtained with the designed flow distributor. The design took into consideration the need to stack the cells, and components were fabricated for a ten-cell stack. Testing of the cell for charge and discharge demonstrated current densities of 50 mA/cm 2 with efficiencies of 50%. However, optimization and long cycling could not be accomplished because the slurry electrode would plug within the cell and block the flow. This would shut down the cell after several charge-discharge cycles. Various approaches were pursued to address this problem with some progress. However, there was not adequate time to find solutions to the problem although several promising approaches were identified for future work. The final deliverable for this proposed 18-month effort was to be a 1 kW, 6 kWh slurry flow battery, with complete balance of plant, ready for testing at an ARPA-E ‘Charges’ site with an estimated capital cost less than $\$ $45/kWh for a six-hour system (excluding power electronics). The deliverable milestone was not achieved for reasons mentioned above related to the slurry plugging. At the completion of this project, it was determined that the cost targets are achievable. The cost estimate for a 1 MW system with a one-hour capacity was estimated at $\$ $140/kWh (including power electronics) and for a six-hour system was estimated at $\$ $40/kWh (including power electronics). These costs include the battery, the balance of plant, the power conditioning and a concrete pad. This program resulted in three patents including the slurry-based iron flow battery, the rebalance reactor, and the membrane separator. At the time of this report, slurry-based iron flow battery patents have been issued in the US, China, Japan, Korea, and Europe. The other patent applications were under various stages of examination. A licensee was identified and the company initiated a scale-up and development program to take this technology to the marketplace, and CWRU staff have transferred their know-how and knowledge base to the company. This technology will enhance the economic and energy security of the United States by enabling the use of intermittent renewable energy technology, such as wind and solar power, as dispatchable resources. This report summarizes the results of this research and development program.

25 ENERGY STORAGE↗

Collaborative Decision Approach for Electricity Pricing-demand Response Stackelberg Game

Demand response programs are considered as a valuable resource in smart grids that provide several advantages of load shifting, peak load reduction, mediating intermittency of renewable energy integration, etc. Flexible price-based incentives have been recognized as a critical strategy in motivating and compensating consumers' load adjustment actions for successful implementation of demand response. Game theoretical approaches, especially Stackelberg games are popularly adopted to model the relationship between electricity price and customers' demand response and solved by the classical centralized backward induction (BI) method. However, the BI method generally requires convexity of the follower's model for necessary optimality conditions, and the computational time of any centralized approach increases sharply with larger problem instances. In this paper, the Stackelberg game of electricity pricing-demand response between a distribution system operator (DSO) and load aggregators (LAs) is decomposed based on a collaborative optimization (CO) framework, where each LA is treated as a discipline with its own domain constraints (e.g. building temperature control), while the DSO at the system level tries to reduce the solution discrepancy and guide the searching towards optimality. Several groups of comparison experiments have demonstrated the effectiveness of the proposed collaborative decision approach in solving the demand response game.

Chen, Yang↗

Photon-Counting Microlaser Rangers, Transponders, and Altimeters

Unlike current manned systems, NASA's next generation SLR2000 Satellite Laser Ranging (SLR) station is fully autonomous. eye-safe, relatively compact and inexpensive. and, during daytime tracking operates at signal-to-noise ratios several orders of magnitude below unity. Tiny, passively Q-switched microlasers generate ultra-short pulses with output energies on the order of 100 micron-J at few kHz rates to achieve mm-level ranging precision to satellite altitudes of 20,000 km. Special ranging receivers, combined with Poisson statistical analysis of the received photon distribution, enable the system to rapidly and reliably identify and extract the single photon laser echoes from the solar background. The enhanced rate of return, combined with a uniform signal strength, can actually drive down both systematic and random range errors. The new SLR2000 technology has already spawned exciting new applications. Compact microlaser altimeters, capable of mapping the surface of a planet or other celestial body at multikilohertz rates, is one such application, and a high altitude, airborne version is currently being developed under NASA's Instrument Incubator Program. Interplanetary microlaser transponders would be capable of performing decimeter ranging or subnanosecond time transfer to spacecraft throughout the inner Solar System. resulting in improved knowledge of planetary motions and liberations and enhanced General Relativity experiments.

Degnan, John J.↗

Linearized Aeroelastic Solver Applied to the Flutter Prediction of Real Configurations

A fast-running unsteady aerodynamics code, LINFLUX, was previously developed for predicting turbomachinery flutter. This linearized code, based on a frequency domain method, models the effects of steady blade loading through a nonlinear steady flow field. The LINFLUX code, which is 6 to 7 times faster than the corresponding nonlinear time domain code, is suitable for use in the initial design phase. Earlier, this code was verified through application to a research fan, and it was shown that the predictions of work per cycle and flutter compared well with those from a nonlinear time-marching aeroelastic code, TURBO-AE. Now, the LINFLUX code has been applied to real configurations: fans developed under the Energy Efficient Engine (E-cubed) Program and the Quiet Aircraft Technology (QAT) project. The LINFLUX code starts with a steady nonlinear aerodynamic flow field and solves the unsteady linearized Euler equations to calculate the unsteady aerodynamic forces on the turbomachinery blades. First, a steady aerodynamic solution is computed for given operating conditions using the nonlinear unsteady aerodynamic code TURBO-AE. A blade vibration analysis is done to determine the frequencies and mode shapes of the vibrating blades, and an interface code is used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor is used to interpolate the mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh. Then, LINFLUX is used to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. Finally, a post-processor uses the unsteady pressures to calculate the generalized aerodynamic forces, eigenvalues, an esponse amplitudes. The eigenvalues determine the flutter frequency and damping. Results of flutter calculations from the LINFLUX code are presented for (1) the E-cubed fan developed under the E-cubed program and (2) the Quiet High Speed Fan (QHSF) developed under the Quiet Aircraft Technology project. The results are compared with those obtained from the TURBO-AE code. A graph of the work done per vibration cycle for the first vibration mode of the E-cubed fan is shown. It can be seen that the LINFLUX results show a very good comparison with TURBO-AE results over the entire range of interblade phase angle. The work done per vibration cycle for the first vibration mode of the QHSF fan is shown. Once again, the LINFLUX results compare very well with the results from the TURBOAE code.

Reddy, Tondapu S.↗