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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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↗
DOE OSTI2025

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↗
DOE OSTI2021

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↗
NASA NTRS2003

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↗
DOE OSTI2018

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↗
DOE OSTI2019

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↗
DOE OSTI2021

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↗
NASA NTRS2000

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.↗
NASA NTRS2004

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.↗
NASA NTRS2012

Analysis of Vegetation Index Variations and the Asian Monsoon Climate

Vegetation growth depends on local climate. Significant anthropogenic land cover and land use change activities over Asia have changed vegetation distribution as well. On the other hand, vegetation is one of the important land surface variables that influence the Asian Monsoon variability through controlling atmospheric energy and water vapor conditions. In this presentation, the mean and variations of vegetation index of last decade at regional scale resolution (5km and higher) from MODIS have been analyzed. Results indicate that the vegetation index has been reduced significantly during last decade over fast urbanization areas in east China, such as Yangtze River Delta, where local surface temperatures were increased significantly in term of urban heat Island. The relationship between vegetation Index and climate (surface temperature, precipitation) over a grassland in northern Asia and over a woody savannas in southeast Asia are studied. In supporting Monsoon Asian Integrated Regional Study (MAIRS) program, the data in this study have been integrated into Giovanni, the online visualization and analysis system at NASA GES DISC. Most images in this presentation are generated from Giovanni system.

Shen, Sunhung↗
DOE OSTI2019

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 sCO 2 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 sCO 2 expander and will guide more robust designs tailored towards sCO 2 cycle applications. A DGS was tested in 1MWe sCO 2 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↗
DOE OSTI2024

Flexible dynamic boundary microgrid operation considering network and load unbalances

Flexible microgrids with dynamic boundaries have recently been introduced in the literature. With the ability to reconfigure the topology of the microgrids dynamically through remotely controlled switches, flexible microgrids with dynamic boundaries can further improve the resiliency and energy efficiency of microgrids with distributed energy resources (DERs). This paper focuses on the optimal operation considering one of the predominant characteristics of microgrids and distribution systems – unbalanced networks and loads. In existing literature, balanced modeling of microgrids is more common due to its attractive simplicity. The three-phase power unbalance has not been considered as a constraint on the generation units in a microgrid. Further, negative sequence constraints have also been neglected. In this article, we propose a set of constraints that is specifically related to the capabilities of inverter interfaced resources to supply unbalanced current/power when the microgrid is islanded from the main distribution grid. We incorporate the new set of constraints into two optimization formulations leveraging two convex relaxations of the three-phase power flow equations: mixed-integer linear programming (MILP) and mixed-integer semidefinite programming (MISDP) that optimize the dispatch of controllable switches and DERs in the microgrid. The algorithms are then extended to networked microgrids with grid-forming sources. We test the algorithms on a realistic community microgrid model in Puerto Rico as well as standardized IEEE distribution test feeders. The testing results demonstrate the performance of the proposed algorithms. The MILP is fast and scalable, and the MISDP enforces the negative sequence voltage constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1988

Galileo Earth-Venus trajectory correction maneuver design

As a result of delays in the Space Shuttle program and a re-evaluation of safety concerns which resulted in cancellation of the Shuttle Centaur Upper Stage program, the Galileo mission to Jupiter has undergone a substantial redesign effort. Use of the Inertial Upper Stage (IUS) for injection has necessitated a complicated Venus-Earth-Earth gravity assist trajectory in order to reach Jupiter. One flyby of Venus and two flybys of Earth are required to compensate for the lower injection energy available from the IUS. As the spacecraft was not originally designed for the environment within 1 AU of the Sun, modifications to the spacecraft and operating procedures have been developed. This paper describes some of these changes and how they have affected the design and implementation of trajectory correction maneuvers on the Earth-Venus leg of the mission. A strategy of biasing the IUS target at injection to control the statistical velocity distribution of a subsequent trajectory correction maneuver is discussed.

Wilson, M. G.↗
DOE OSTI2024

HiFiPDV2 v. 2

SAND2024-01061O The HiFiPDV-2 software is similar to other photonic Doppler velocimetry (PDV) analysis codes in that it converts signals from intensity-time space to frequency-time space using a short-time Fourier transform (STFT). It then analyzes peaks of the resulting power spectral density information to extract the velocity-time history. The unique aspects of HiFiPDV and HiFiPDV-2 are that this reduction process is repeated for an array of reasonable STFT input variables. The resulting population distribution of output velocity-time histories is evaluated to determine the most likely velocity history and its associated uncertainty. This uncertainty is defined as the systematic uncertainty of the PDV signal as it is a function of input values to the reduction process. The HiFiPDV-2 program brings new functionality and efficiency to these calculations—most notably significant reductions in memory usage and calculation times, the ability to evaluate frequency-shifted signals, and the ability to filter out baseline signals. Developers recommend that users have at least 2 GB of RAM per CPU core in their system. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Voorhees, Travis↗
DOE OSTI2018

Tracking the Sun: Installed Price Trends for Distributed Photovoltaic Systems in the United States - 2018 Edition

Berkeley Lab’s Tracking the Sun report series is dedicated to summarizing installed prices and other trends among grid-connected, distributed solar photovoltaic (PV) systems in the United States. The present report, the 11th edition in the series, focuses on systems installed through year-end 2017, with preliminary trends for the first half of 2018. As in years past, the primary emphasis is on describing changes in installed prices over time and variation in pricing across projects based on location, project ownership, system design, and other attributes. New to this year, however, is an expanded discussion of other project characteristics in the large underlying data sample. Future editions will include more of such material, beyond the report’s traditional focus on installed pricing. The trends described in this report derive primarily from project-level data reported to state agencies and utilities that administer PV incentive programs, solar renewable energy credit (SREC) registration systems, or interconnection processes. In total, data were collected and cleaned for more than 1.3 million individual PV systems, representing 81% of U.S. residential and non-residential PV systems installed through 2017. The analysis of installed pricing trends is based on a subset of roughly 770,000 systems with available installed price data. Key findings from this year’s report include the following: -Installed Prices Continued to Fall through 2017 and into 2018. National median installed prices in 2017 were $3.7 per watt (W) for residential systems (a $0.2/W or 6% decline from the prior year), $3.1/W for small non-residential systems (a $0.4/W or 11% decline), and $2.2/W (a $0.1/W or 5% decline) for large non-residential systems. Data for the first half of 2018 show an additional drop of $0.1/W for residential and small non-residential systems, and effectively no change for large non-residential systems. These recent trends are generally consistent with the pace of price declines since 2014, and mark a slowing from the years immediately preceding (2009-2013) when prices fell by roughly $1/W per year. -Installed Prices Vary Widely Across Individual Projects. Among residential systems installed in 2017, 20% were priced below $3.0/W (the 20th percentile value) while another 20% were above the 80th percentile at $4.5/W. Non-residential systems also exhibit wide spreads, albeit shifted downward, from $2.4/W to $4.1/W for small non-residential and from $1.8/W to $2.8/W for large non-residential projects. These pricing spreads have persisted over time and are both reflective of the heterogeneity in PV systems and markets, as well as suggestive of the potential for further price declines. -Pronounced Pricing Disparities Exist Across States. State-level median installed prices in 2017 ranged from $2.6/W to $4.5/W for residential systems, from $2.2/W to $4.0/W for small non-residential systems, and from $2.1/W to $2.4/W for large non-residential systems. Three of the largest state markets (California, Massachusetts, and New York) are relatively high-priced, pulling overall U.S. median prices upward. These cross-state pricing differences reflect both idiosyncratic features of particular states as well as more-fundamental differences related to market and policy conditions. -Installed Pricing Variation Also Reflects Differences in System Design, Installer, and Customer Characteristics (among other factors). Significant pricing differences were also observed across system sizes within each customer segment (with median prices ranging from $3.2/W to $4.5/W across residential systems of varying size, and from $2.1/W to $3.7/W across non-residential systems), between third-party owned (TPO) and host-owned residential systems ($3.3/W vs. $3.8/W), across residential installers ($2.1/W to $9.6/W for host-owned systems and $1.1/W to $5.1/W for TPO systems), between residential new construction and retrofit systems in California ($2.3/W vs. $3.9/W), between large non-residential systems at commercial sites and tax-exempt customer sites ($2.1/W vs. $2.6/W), and between systems with standard vs. premium efficiency modules ($3.6/W vs. $4.2/W among residential systems). A public version of the underlying dataset can be downloaded at https://emp.lbl.gov/tracking-the-sun.

14 SOLAR ENERGY↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2023

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2000

Theoretical Technology Research for ISTP/SOLARMAX

During the last decade, we have been developing theoretical tools to support the scientific objectives of the International Solar Terrestrial Physics (ISTP) program. Results from our mission-oriented theory program have contributed significantly to the development of predictive capabilities by using real upstream solar wind conditions as input to our models and forecasting events observed downstream near Earth. We also developed the capability to unravel the complex information contained in ion velocity distribution functions measured near the Earth to determine their origin and energization process. During solar maximum, solar flares and coronal mass ejections (CMEs) dominate the sun's activity. It is now widely accepted that the impact of CMEs (or magnetic clouds) with the Earth's magnetosphere is the cause of most magnetic storms during solar maximum. One important aspect of a CME is the occurrence of solar energetic particle (SEP) events. During these events, protons, electrons, and heavy ions of solar origin are accelerated to very high energies by shock waves driven out from the sun. We carried out a series of large-scale kinetic (LSK) simulations to model the effect of SEPs on the near-Earth environment and the accessibility of these high-energy particles to the inner magnetosphere. We present the results of these studies.

Ashour-Abdalla, Maha↗