Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control of Distributed Energy Resources”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Potential Impact of Offshore Wind Energy on a Future Power System in the U.S. Northeast

This study aims to understand and quantify the potential impact of offshore wind on a future electricity system in the U.S. Northeast. In this analysis, a detailed representation of the Northeast power system is adopted, using a generation portfolio for 2024 paired with offshore wind nameplate capacities of 0 gigawatts (GW), 2 GW, and 7 GW. The analysis identifies points of offshore wind interconnection in the ISO-New England and New York Independent System Operator control areas and uses hourly wind profiles from the Wind Integration National Dataset Toolkit. Hourly simulations of the 2024 power system operations show an ability to accommodate the prescribed offshore wind capacities by adapting the system's generation dispatch. Curtailment levels of offshore wind range between 4% and 5%. Offshore wind generation displaces primarily natural-gas combined-cycle generation; however, requires increased flexibility from combined cycles through more frequent start-ups. The number of hours with transmission congestion increases because of offshore wind injection, with varying impact on a subregional level. Offshore wind's capacity credit was found to be 14.5%-28.3% and is lower than estimated in other large-scale power system studies, in part because this study considered a different weather year and a power system representation with higher shares of solar PV and onshore wind. The reliability contribution of resources is a topic area identified for more research and collaboration. The 7-GW scenario shows a reduction in locational marginal price of 11%, with production cost savings of up to 18% compared to the 0-GW scenario.

17 WIND ENERGY↗
DOE OSTI2025

A Comprehensive Study of the Impact of Inverter-Based Resource (IBR) Modeling and Control on Protection Relay Elements

Inverter-based resources (IBRs) exhibit fault responses that differ significantly from those of synchronous generators, which can challenge the reliable operation of many commonly used power system protection elements. The fault response of IBRs is primarily influenced by their control algorithms and configurations, but the impact of these controls on protection relays is not yet fully understood. This presentation provides a comprehensive study of how IBR modeling and controls affect transmission line protection. Key modeling and control aspects include the DC source, inverter model, power level control, current control, and current limiting. The study reveals that certain aspects - such as the type of DC source (battery, PV, or hybrid), inverter model (average vs. switching), and power level control methods (PQ dispatch vs. Vdc-Vac control for grid-following IBRs, and droop vs. VSM for grid-forming IBRs) - do not significantly affect relay response. However, faster control loops, such as current control and current limiting, do influence the relay behavior. Additionally, the study explores the effects of other factors, including momentary cessation, operating points, fast/slow current responses, and grid strength on relay performance. Finally, the study offers recommendations for both IBR and protection engineers to improve IBR fault response and enhance the reliability of protection systems.

14 SOLAR ENERGY↗
DOE OSTI2020

CyDER: A Cyber Physical Co-simulation Platform for Distributed Energy Resources in Smartgrids

The CyDER project aimed at developing an open-source, modular and scalable co-simulation platform for power grids with large shares of Distributed Energy Resources (DERs). The project partners are the Lawrence Berkeley National Lab (LBNL), Lawrence Livermore National Lab (LLNL), PG&E, SolarCity, and ChargePoint. The prime recipient is LBNL; SolarCity and ChargePoint were partners for the project’s first two years. Increased DER integration introduces a number of challenges in power grid operation including a more dynamic interaction between the transmission grid and distribution grids, and increased modeling complexity. Although specialized software exists to precisely model different components of the power system, it is far from trivial to integrate all various models and perform a holistic simulation. Instead of replicating all models in a common simulation program, a commonly accepted approach to tackle this model diversity is to couple third-party simulators and models through a co-simulation platform that coordinates information exchange among the various components. Following this line of research, this project’s objective was to develop a co-simulation platform based on a widely accepted industrial standard called Functional Mockup Interface (FMI). Within this process, the project developed models compliant with the FMI standard, called Functional Mockup Units (FMUs), and used them to perform various operational and planning power system analyses. Relying and building upon an industrial standard is the main differentiation of this project compared with previous or parallel efforts in the co-simulation area. Particular emphasis was put on delivering software utilities to facilitate setting up and running co-simulations by end-users. Furthermore, a strong aspect of this project is demonstrating that co-simulation techniques can be used to perform Hardware-in-the-Loop (HIL) simulations that couple software components (e.g., simulated models) with hardware components (e.g., real devices such PV systems and batteries). The long-term goal of CyDER project is to help establish FMI as a powerful standard for co-simulation and promote adoption by electric utilities and other interested stakeholders. The main accomplishments of the project include the development of several FMUs including distribution and transmission grid models, PV inverters with Volt/Var/Watt controllers, batteries, and predictive optimal controllers. Additionally, a unique software package was developed, called SimulatorToFMU, which is capable of exporting any Python-driven simulator or Python script as an FMU. This is an important contribution towards establishing FMI as one of the main co-simulation standards, because more and more third-party programs for sub-system modeling and simulation are delivered with Python APIs. The CyDER platform was used to perform PV hosting capacity analyses in real utility feeders with and without smart inverter controls, battery storage, and EV charging. Smart inverter controls include conventional Volt/Var/Watt controls for reactive power support and active power curtailment, but also predictive controls that optimize the charging and discharging profile of the battery connected on the DC side in order to minimize the customer’s economic benefit. Finally, an important result of this project is delivering an experimental setup that consists of residential-scale PV inverters with battery storage, a real-time grid simulator with an ideal voltage source as grid emulator, and micro Phasor Measurement Units (PMUs). All these components and additional software modules are coupled to one another using the FMI standard and can be co-simulated with the CyDER platform.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Utilities Perspective on Protection Challenges with High IBR Penetration

Utilities have seen rapid increase in Solar, Wind, and Battery Energy Storage Resources that interconnect to their electric system through Inverters. Inverter-Based Resources (IBRs) have fault current characteristics that are unlike the fault current response of traditional rotating-machine-based generators, which is well known and repeatable. IBR’s non-traditional fault current behavior is due to the IBR control scheme, which is configured to provide a clean AC output but also protect the inverter’s sensitive power electronics devices from damage, one source of which is overcurrent. This results in low fault current magnitude, low or no negative sequence current injection, the variability of sequence component currents, the variability of voltage with respect to current angles, and the lack of inertia. The control scheme also results in a fault current response that can vary between manufactures and between models of the same manufacturer. High penetration of IBRs can adversely affect the protection schemes applied in areas with high penetration of IBRs. With the proliferation of IBRs, utilities are finding out that conventional protection schemes are not adequately equipped to protect the electric systems. This is mainly because the existing protection elements and practices have been designed based on the fault current response of conventional rotating machines. In several cases, the available literature does not provide any clear solution for the issues when the protection scheme does not operate properly near IBRs. This presentation identifies various protection challenges due to IBRs that industry is facing, from the utility perspective. Instead of facing on one issue, we are looking broadly on all the challenges that system protection has experienced with high penetration of IBRs. Based on the IBR response from various utilities during real fault events and gathering perspective from different utility SMEs via questionnaire, the presentation summarizes on gathered data and internal experiences.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Characterizing Current THD’s Dependency on Solar Irradiance and Supraharmonics Profiling for a Grid-Tied Photovoltaic Power Plant

The rapidly increasing distributed energy resources (DERs) in power systems are now getting interconnected to set community grid structures, where power quality will be a major concern. The grid-to-grid (G2G) bidirectional power transfer among the distribution microgrid will not be considered commercially feasible unless the upstream harmonics are under the limits. The aggregation of such harmonics, measured as total harmonic distortion (THD), is feared to be beyond tolerable limits with the progression of rooftop grid-tied PV-like installations. Hence, this THD needs to be characterized with DER generation end variables. In this work, the photovoltaic (PV) DERs’ dependency on environment variables such as irradiance was profiled in the context of generating and injecting harmonics into the grid. A mathematical model of a grid-tied three-phase PV DER was developed as part of this correlation characterization, matching the fundamental unit structure of a 1.4 MW solar canopy located on the Florida International University (FIU) Miami campus. To determine the qualitative association with produced THD patterns, the model was evaluated with various irradiance settings. A real-time digital simulation (RTDS) platform was used to verify it. Following this confirmation, sets of data from power quality meters at the point of common coupling and FIU field sensors were utilized to validate further the correlation model. The results showed that the grid current’s THD exhibited a high correlation with the irradiance profile and its variation over time. The early morning and late afternoon periods of the day, associated with a low irradiance, constantly had higher harmonics generated from the PV DER. The midday THD was rather rational with partial shadings, hence a geolocation-dependent factor. These findings were verified by an RTDS and validated by real field data. In quantifying the THD injected by a single DER at a high-frequency (2–150 kHz) supraharmonics (SH) level, a 3% peak increment in magnitude was observed from the high-fixed to the low-fixed irradiance profile. The correlation characteristics depicted that the hybrid microgrid suffered from a daytime-dependent harmonic insertion from the grid-tied DER. This is a global problem unless specific measures are taken to mitigate the harmonics. The electrically notorious higher-frequency SH was found to increase proportionally. The G2G power transfer can be limited because of the higher THD in the early morning and late afternoon, which will also worsen because the numbers of grid-tied PV DERs (i.e., rooftop solar and industrial solar) are likely to increase rapidly soon. The community grid structure can thus have a controlled harmonics filtration setup purposefully designed to address the findings of this work, which also fall within the scope of our future research.

14 SOLAR ENERGY↗
DOE OSTI2024

DOE 10658 Project Sensitivity-Driven Wide Area Protection (SWAP) Coordination Tools for IBRs

Utilities have seen a rapid increase in solar, wind, and battery energy storage systems that interconnect to their electric system through inverters. Inverter-Based Resources (IBRs) have fault current characteristics that are unlike the fault current response of traditional rotating-machine-based generators, which is well-known and repeatable. IBR’s non-traditional fault current behavior is due to the IBR control scheme, which is configured to provide a clean AC output but also protects the inverter’s sensitive power electronics devices from damage from conditions such as overcurrent. This results in a low fault current magnitude, low or no negative sequence current injection, the variability of sequence component currents, the variability of voltage with respect to current angles, and the lack of inertia. The control scheme also results in a fault current response that can vary between manufacturers and between models of the same manufacturer. High penetration of IBRs can adversely affect the protection schemes applied in areas with high penetration of IBRs. With the proliferation of IBRs, utilities are finding out that conventional protection schemes are not adequately equipped to protect the electric systems. This is mainly because the existing protection elements and practices have been designed based on the fault current response of conventional rotating machines. In several cases, the available literature does not provide any clear solution for the issues when the protection scheme does not operate properly near IBRs. This report identifies various protection challenges due to IBRs that industry is facing, from the utility perspective. Instead of focusing on one issue, this report looks broadly at all the challenges that system protection has experienced with high penetration of IBRs. Based on the IBR response from various utilities during real fault events and gathering perspectives from different utility SMEs via questionnaire, the report summarizes the gathered data and internal experiences.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Power-Hardware-in-the-Loop Experiments of a Microgrid With a Grid-Forming Battery Inverter: Preprint

Microgrids continue to proliferate, and they are transitioning away from using conventional generating resources to increasingly relying on inverter-based resources (IBRs) as the voltage and frequency leaders. It is crucial to evaluate the capability of IBRs to provide microgrid stability and resilience. Hardware-in-the-loop (HIL) experiments were conducted to de-risk the field deployment of the San Diego Gas & Electric Company Borrego Springs Microgrid, where a battery inverter was upgraded with grid-forming (GFM) capability to serve as the island leader. This paper presents the HIL experimental results from an HIL test bed that uses a power-hardware-in-the-loop (PHIL) interface with a power inductor that was previously developed for PHIL simulations of microgrids where the inverters need to switch modes, i.e., between grid-following and GFM as the microgrid transitions between grid-connected and islanded operation. This paper presents more details on the interface and HIL simulation results of the planned islanding and load steps in islanded operation to show the effectiveness of the inverters in managing the voltage and frequency.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Simple, Secure, Internet Delivery of MOOSE-based Applications

Application packaging and distribution are the final steps for delivering software to end-users; both are frequently neglected when creating scientific software. Commercial businesses rely on electronic distribution systems that have rendered disk drives obsolete. Still, national laboratories continue to rely heavily on removable media to distribute and limit access to controlled applications. With increasing concerns of unauthorized copying of sensitive applications, a modern distribution system that utilizes cryptographically secure communication and authentication protocols has been developed. This new distribution system will secure the chain of custody for nuclear software while simultaneously simplifying access to these tools. This report summarizes four primary advancements made toward the secure distribution of Nuclear Energy Advanced Modeling and Simulation (NEAMS)-developed, Multiphysics Object Oriented Simulation Environment (MOOSE)-based applications: application installation, package distribution, automated package building, and distribution of documentation. NEAMS is currently developing more than ten separate applications based on the open-source MOOSE Framework. Distribution of these applications has primarily been accomplished by distributing source code, with end-users compiling the applications themselves. This work created a mechanism where MOOSE applications can be installed in a similar way to any other software. This allows both administrators and end-users simplified access to runnable executables. With this new installation capability, it was then possible to rethink distribution. A new, secure capability for delivering MOOSE-based applications over the internet has been created. This system requires unique cryptographic tokens for authentication, greatly securing the custody chain for software. Once granted access, installation of any NEAMS code can be accomplished with these terminal commands: "conda install ncrc" "ncrc install ncrc-bison." After these two commands (and authenticating) the BISON application will be securely down- loaded from Idaho National Laboratory (INL)’s servers, installed, and ready to use. To enable this new distribution capability to be successful, the open-source Continuous Integration, Verification, Enhancement, and Testing (CIVET) Continuous Integration (CI) capability was augmented to add Continuous Delivery (CD). CD enables the automated building and packaging of MOOSE-based applications as they are modified by development teams, ensuring that our customers can obtain up-to-date versions of the software at any time. The need for instruction on how to use these applications was addressed through modifications to the MOOSE documentation system. The MooseDocs capability, which enables robust documentation of MOOSE-based applications, has been extended to allow both for the installation of documentation and the packaging of documentation with installed applications. Together, these enhancements form the core of a new, secure distribution mechanism for nuclear simulation tools. In concert with the Nuclear Computational Resource Center (NCRC), NEAMS- developed applications will now be straightforward to obtain securely.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Power-Hardware-in-the-Loop Experiments of a Microgrid with a Grid-Forming Battery Inverter

Microgrids continue to proliferate, and they are transitioning away from using conventional generating resources to increasingly relying on inverter-based resources (IBRs) as the voltage and frequency leaders. It is crucial to evaluate the capability of IBRs to provide microgrid stability and resilience. Hardware-in-the-loop (HIL) experiments were conducted to de-risk the field deployment of the San Diego Gas & Electric Company Borrego Springs Microgrid, where a battery inverter was upgraded with grid-forming (GFM) capability to serve as the island leader. This paper presents the HIL experimental results from an HIL test bed that uses a power-hardware-in-the-loop (PHIL) interface with a power inductor that was previously developed for PHIL simulations of microgrids where the inverters need to switch modes, i.e., between grid-following and GFM as the microgrid transitions between grid-connected and islanded operation. This paper presents more details on the interface and HIL simulation results of the planned islanding and load steps in islanded operation to show the effectiveness of the inverters in managing the voltage and frequency.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2014

Lunar Flashlight: Mapping Lunar Surface Volatiles Using a Cubesat

Water ice and other volatiles may be located in the Moon's polar regions, with sufficient quantities for in situ extraction and utilization by future human and robotic missions. Evidence from orbiting spacecraft and the LCROSS impactor suggests the presence of surface and/or nearsurface volatiles, including water ice. These deposits are of interest to human exploration to understand their potential for use by astronauts. Understanding the composition, quantity, distribution, and form of water/H species and other volatiles associated with lunar cold traps is identified as a NASA Strategic Knowledge Gap (SKG) for Human Exploration. These polar volatile deposits could also reveal important information about the delivery of water to the Earth- Moon system, so are of scientific interest. The scientific exploration of the lunar polar regions was one of the key recommendations of the Planetary Science Decadal Survey. In order to address NASA's SKGs, the Advanced Exploration Systems (AES) program selected three lowcost 6-U CubeSat missions for launch as secondary payloads on the first test flight (EM1) of the Space Launch System (SLS) scheduled for 2017. The Lunar Flashlight mission was selected as one of these missions, specifically to address the SKG associated with lunar volatiles. Development of the Lunar Flashlight CubeSat concept leverages JPL's Interplanetary Nano- Spacecraft Pathfinder In Relevant Environment (INSPIRE) mission, MSFC's intimate knowledge of the Space Launch System and EM-1 mission, small business development of solar sail and electric propulsion hardware, and JPL experience with specialized miniature sensors. The goal of Lunar Flashlight is to determine the presence or absence of exposed water ice and its physical state, and map its concentration at the kilometer scale within the permanently shadowed regions of the lunar south pole. After being ejected in cislunar space by SLS, Lunar Flashlight deploys its solar panels and solar sail and maneuvers into a low-energy transfer to lunar orbit. The solar sail and attitude control system work to bring the satellite into an elliptical polar orbit spiraling down to a perilune of 30-10 km above the south pole for data collection. Lunar Flashlight uses its solar sail to shine reflected sunlight into permanently shadowed regions, measuring surface albedo with a four-filter point spectrometer at 1.1, 1.5 1.9, and 2.0 microns. Water ice will be distinguished from dry regolith from these measurements in two ways: 1) spatial variations in absolute reflectance (water ice is much brighter in the continuum channels), and 2) reflectance ratios between absorption and continuum channels. Derived reflectance and reflectance ratios will be mapped onto the lunar surface in order to distinguish the composition of the PSRs from that of the sunlit terrain. Lunar Flashlight enables a low-cost path to in-situ resource utilization (ISRU) by identifying operationally useful deposits (if there are any), which is a game-changing capability for expanded human exploration.

Cohen, B. A.↗
DOE OSTI2025

Real-World Cyber Security Demonstration for Networked Electric Drives

In this article, we present the design and implementation of a cyber-physical security testbed for networked electric drive systems, aimed at conducting real-world security demonstrations. To our knowledge, this is one of the first security testbeds for networked electric drives, seamlessly integrating the domains of power electronics and computer science, and cybersecurity. By doing so, the testbed offers a comprehensive platform to explore and understand the intricate and often complex interactions between cyber and physical systems. The core of our testbed consists of four electric machine drives, meticulously configured to emulate small-scale but realistic information technology (IT) and operational technology (OT) networks. This setup both provides a controlled environment for simulating a wide array of cyber-attacks, and mirrors potential real-world attack scenarios with a high degree of fidelity. The testbed serves as an invaluable resource for the study of cyber-physical security, offering a practical and dynamic platform for testing and validating cybersecurity measures in the context of networked electric drive systems. As a concrete example of the testbed's capabilities, we have developed and implemented a Python-based script designed to execute step-stone attacks over a wireless local area network (WLAN). This script leverages a sequence of target IP addresses, simulating a real-world attack vector that could be exploited by adversaries. To counteract such threats, we demonstrate the efficacy of our developed cyber-attack detection algorithms, which are integral to our testbed's security framework. Furthermore, the testbed incorporates a real-time visualization system using InfluxDB and Grafana, providing a dynamic and interactive representation of networked electric drives and their associated security monitoring mechanisms. This visualization component not only enhances the testbed's usability but also offers insightful, real-time data for researchers and practitioners, thereby facilitating a deeper understanding of cyber-physical security dynamics in networked electric drive systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Multi-Timescale Integrated Dynamics and Scheduling for Solar (MIDAS-Solar) (Final Technical Report)

Solar photovoltaic (PV) installations have experienced unprecedented growth in the United States. PV will become not only an energy producer but also a necessary provider of ancillary services at multiple timescales. Conventional methods to simulate power systems operations - such as long-term production simulation (which typically considers schedules from hours to minutes by using an optimization framework) and short-term transient studies (which simulate dynamics from seconds to sub-seconds using state variables and differential equations) - are not sufficient for studying the multiple-timescale variation of solar generation and its impact on system reliability. Long-term system economics and short-term system dynamics are highly coupled, particularly when the penetration level of renewable generation is extremely high, because the uncertainty and variability of solar generation will impact both power system steady-state and dynamic performance. This project helps meet and exceed the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Solar Energy Technologies Office goal of systems integration by directly addressing this stability and reliability challenge for power grid planning and operation. We have developed a temporally comprehensive, closed-loop simulation model, named Multi-timescale Integrated Dynamics and Scheduling (MIDAS), that seamlessly simulates power system operations from economic scheduling (day-ahead to hours) to dynamic response analysis (seconds to sub-seconds). For schedules with very high levels of inverter-based resources (IBRs), up to and including 100%, the stability of grid controls has been evaluated through electromagnetic transient (EMT) simulations and power-hardware-in-the-loop (PHIL) simulations of key transient events at key schedule points. Specifically, MIDAS provides: 1) a closed-loop simulation framework for simulating timescales from economic scheduling to dynamic stability analysis; 2) machine learning-based stability assessment; 3) EMT modeling and analysis for large-scale power systems; 4) MIDAS PHIL test bed. We worked with Hawaii Electric Companies to apply the MIDAS study framework to a Maui grid study. The entire island's transmission system was modeled in detail - from a yearly scheduling model, to a second-level frequency dynamic model, down to a sub-second-scale EMT model to address critical stability issues. The project demonstrated how MIDAS can help system planners and operators assess system reliability and stability while the power grid is marching toward a high-renewable, high-IBR future. In this Maui grid study, we found that 100% instantaneous IBR operation is achievable in EMT simulation and PHIL testing, and grid planners and operators might need new analysis/simulation tools to assess grid reliability and stability in the scheduling stage. MIDAS will bring Maui and other systems closer to 100% clean and stable energy futures. (In this study, we examined transient stability. Other topics necessary for 100% IBR operation, such as protection and resource adequacy, were not examined.)

100% Renewables↗
DOE OSTI2020

Demand Response Optimization and Management System for Real-TIme (DROMS-RT)

To design and demonstrate DROMS-RT, a highly distributed Demand Response Optimization and Management System for Real-Time (DROMS-RT) power flow control to support large scale integration of distributed renewable generation into the grid. AutoGrid developed a novel control and communications platform to allow highly dispatchable demand response (DR) services in time frames suitable for providing ancillary services to the transmission grid. These services will be substantially less expensive and more efficient than other forms of ancillary services options currently available to manage the intermittency associated with large-scale renewable integration. DROMSRT successfully leveraged Automated Demand Response (ADR) by fundamentally re-thinking the architecture of the DR platform from the ground up and by developing innovative new technologies in a number of areas related to DR. DROMS-RT leveraged the low-cost, open, interoperable DR signaling technology, OpenADR, and low-cost, internet-protocol based telemetry solutions to reduce the cost of hardware. This allowed DROMS-RT to provide dynamic price signals to millions of OpenADR clients. Statistically rigorous signal processing techniques were developed to reliably detect even small load reductions in the presence of noisy baseline profiles. Novel forecasting engines based on modern online machine learning algorithms enabled accurate individualized forecasts for customer loads in the presence of dynamic pricing signals, and a real-time decision engines enabled continuous optimization and optimal dispatch of DR resources across a large portfolio of heterogeneous loads that respond at varying time-scales. Moreover, the real-time optimization conducted by the decision engine can utilize grid physics to maximize load reduction at the transmission system in addition to the distribution sites, for more efficient grid operation. Finally, the Software-as-a-Service (SaaS) availability of the DROMS-RT platform has reduced the cost of deployment and enable participation of small commercial and residential customers in DR who otherwise would not be able to do so.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

GT Flex: A Coordinated Multi-Building Pilot Study

Buildings are a significant and untapped resource for providing utility electric grid services. Recent studies have estimated that buildings could reduce the peak demand on the electric grid in the U.S. by almost 25% through effective combinations of energy efficiency (EE) measures and load flexibility strategies (Langevin et al. 2021). The U.S. Department of Energy (DOE) has established a goal to triple energy efficiency and demand flexibility in both residential and commercial sections by 2030 compared to 2020 levels (Satchwell et al. 2021). Such findings place buildings alongside electric vehicles, photovoltaics, electric batteries, and other distributed energy resources (DERs) as primary technologies needed for supporting high renewable energy generation grids. Coordinating and optimizing multiple buildings and other DERs is more beneficial and valuable when compared with individual buildings and DERs operating as siloed resources, uncoordinated with others (Olgyay et al. 2020). A pilot study at the Georgia Institute of Technology (GIT) was conducted to evaluate value propositions of a multi-building scale project seeking carbon reduction, energy efficiency and grid-interactive capabilities, by demonstrating the means by which stakeholders can determine the technical and financial merits of grid-interactivity and energy efficiency technologies coordinated across multiple assets. The study focused on analyzing technical feasibility of deploying thermal load flexibility strategies at the multi-building scale, coordinated to not exceed existing infrastructure constraints at the pilot site. Results show that campus can provide 3-3.5 MW of potential load shed over a 4-hour event window through coordinated dispatch of thermal cooling load flexibility without exceeding existing infrastructure capacities. Under future high renewable scenarios, this thermal flexibility resource is also valuable when coordinated to reduce curtailment of intermittent renewables. Economic analyses were performed to effectively communicate various value propositions of grid-interactive and efficient building (GEB) thermal flexibility strategies. Load flexibility presents a financial value proposition to the campus today. By conducting rationalized, coordinated dispatch in response to real time price (RTP) fluctuations, the campus can benefit materially from daily price arbitrage. The RTP signal acts as an aggregating mechanism between the utility and customer to call on demand flexibility resources, with a large portion of the benefit deriving from a relatively small number of days. Realizing and maximizing this benefit with thermal load flexibility requires careful attention to the timing of pricing signals and parameterization of dispatch to overcome efficiency penalties. Grid value and signals are expected to evolve over time, and thermal load flexibility shows potential to adapt dispatch logic to support intermittent renewable generation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Rapid Load Transition for Integrated Solid Oxide Fuel Cell - Gas Turbine (SOFC-GT) Energy Systems: A Demonstration of the Potential for Grid Response

Rapid load transition is an essential requirement for integrated energy systems to maintain grid resilience as more renewable resources are added to the grid. Integrated solid oxide fuel cell - gas turbine (SOFC-GT) systems can provide high efficiency and low carbon emissions over a broad range of turndown. These hybrids also have the potential to enable rapid grid response. The challenge has been to demonstrate effective control strategies to manage load transitions. In the present study, a load transition of ~50% was achieved in 10 s using a novel but simple strategy. Power demand on the SOFC and the GT were ramped down concurrently. During this transition, the SOFC anode fuel was manipulated to maintain SOFC fuel utilization while the cathode inlet air flow and temperature were also manipulated to thermally protect the SOFC. This study was conducted using the Hybrid Performance (Hyper) facility at the National Energy Technology Laboratory (NETL) in a co-simulation environment with the Idaho National Laboratory (INL)'s grid-simulation. The load ramping strategy was tested using a hardware-based cyber-physical simulation methodology. The results demonstrate a high-fidelity representation of SOFC-GT hybrid dynamics and validation of the control strategy. Thermal and electrochemical transients indicated that the SOFC was well protected during rapid load turndown without violating operability constraints. This demonstration revealed the non-linear nature of tightly coupled SOFC-GT system components, especially the non-linear response of SOFC cathode air flow and inlet temperature controls. These results highlight the needs and challenges in developing adaptive automatic controls for autonomous rapid load transitions. This work demonstrates that SOFC-GT hybrids are a viable option to provide the fast-ramping characteristics essential to accommodate high levels of variable renewable power while maintaining grid resilience, reliability, and environmental performance. The results also demonstrate the utility of co-simulation in advancing the tightly-coupled integrated energy systems needed to meet goals for zero-carbon power generation.

DIRECT ENERGY CONVERSION,POWER TRANSMISSION AND DI↗
DOE OSTI2023

Investigation on Tuning Power–Frequency Droop for Improved Grid–Forming Inverter and Synchronous Generator Transient Load Sharing

Poor transient load sharing has been observed during the parallel operation of synchronous generators (SGs) and droop-controlled grid-forming (GFM) inverter-based resources (IBRs) in islanded microgrids. This can result in overcurrent protections, causing the inverters to trip, which leads to the collapse of voltage and frequency. In this article, small-signal analysis supported by electromagnetic transients (EMTs) simulation is performed on a detailed model of a microgrid containing IBRs and an SG to identify underdamped modes and their sensitivities to the power–frequency (P-F) droop parameter. Further EMT simulation of a more complex microgrid model with multiple GFM IBRs and SGs is performed to determine whether trends identified from the modal analysis extend to cases where the number of SGs in operation may vary. The objective is to investigate the effectiveness of tuning this parameter to reduce the likelihood of overcurrents in response to load and generation disturbances. The work uncovers both fast and slow modes of concern. Primarily, the findings show that increasing the P-F droop value improves the damping of the slower mode while degrading the damping of the faster mode, leading to an upper limit on the damping improvement of the GFM IBR output current. In EMT simulations of the more complex microgrid, the greatest damping improvement still exhibits significant maximum overshoot. Furthermore, the optimal value of the P-F droop parameter is sensitive to the number of SGs in operation. The conclusion is that tuning the P-F droop parameter is not an effective means to avoid the overcurrent tripping of GFM IBRs during response to large load or generation disturbances when operating in parallel to SGs in islanded microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-zirconium (Nb-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗