Results for “load uncertainty” — Page 27

DOE OSTI2022

Measurement of the Effective Capture Cross Section of 238 Np in the High Flux Isotope Reactor

The cross sections of 237 Np and 238 Np are important for accurate modeling and simulation of 238 Pu in the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL). Uncertainties in these cross sections can impact the ability to predict and optimize the target design and loading for 238 Pu production targets. The effective capture cross section of 237 Np in the location of pneumatic tube 1 in HFIR was measured as a first step in the measurement of the 238 Np capture and fission cross sections. Here, we describe the flux measurements, 237 Np experiments, and data analysis of the 237 Np capture cross section in HFIR.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Cost Drivers for a Tokamak-Based Compact Pilot Plant

A comprehensive systems code that includes a range of physics and engineering considerations along with a simplified costing model has been utilized to evaluate the primary cost drivers for a compact tokamak pilot plant. The systems code has been benchmarked against several tokamak reactor designs and is utilized with sophisticated optimization algorithms to develop optimal solutions for a set of user-specified assumptions and design constraints. Here, in contrast to previous models that have focused on the cost of electricity as the key cost metric, this study uses the estimated capital cost of the facility. The analysis suggests that a pilot plant with the following features may offer potential for a cost-attractive pilot plant: A ~ 3, H 98y2 > 1.5, P net = 200 MW, 𝜏 𝑝⁢𝑢⁢𝑙⁢𝑠⁢𝑒 ~ 1 to 2 h utilizing rare-earth barium copper oxide (REBCO) magnet technology and a plug-bucked central solenoid/toroidal field (CS/TF) magnet support structure. While REBCO magnets offer some advantages relative to Nb 3 Sn magnets in all cases, the most gain is obtained when combined with the plug-bucked CS/TF bucking solution. Pulsed operation reduces capital cost requirements relative to steady-state operation, especially at low confinement. Cost sensitivity studies indicate that there are significant cost uncertainties associated with the achievable confinement quality, tritium breeding capability, attainable thermal efficiency, and achievable neutron wall loading, suggesting that these areas are the most critical areas in reducing the cost risk for a compact tokamak pilot plant. Further cost sensitivity studies indicate that the estimated cost is most sensitive to the underlying cost of the magnetic coils, providing further impetus to better establish cost-effective means for producing fusion magnets.

Compact pilot plant↗

DOE OSTI2021

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗

DOE OSTI2020

Neutron Spectrum and 10B(n,α) Measurements in ZPPR-15

Neutron spectra were measured at the core center with proton recoil spectrometers in ZPPR15A loading 66, ZPPR-15B loading 90 and ZPPR loading 186. In addition, neutron spectra were measured in the SSNA shield configuration in ZPPR-15B loading 121 and in the BCNA shield configuration in ZPPR-15B loading 126. The ZPPR-15 neutron spectrum measurement and data processing techniques were described and analyzed. Total uncertainties were derived for all five spectrum measurements. A back-to-back proportional counter was used to measure the 10 B(n,α) reaction rate at six locations in the simulated shield, radial reflector and outer core in the ZPPR-15B BCNA shield experiment. The measurements were described and analyzed. A total uncertainty was derived for the 10 B(n,α) reaction rate measurements. The ZPPR-15 experimental records were used to develop detailed as-built Monte Carlo models for the five neutron spectrum measurement configurations and for the six 10 B(n,α) measurement configurations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1991

Parametric entry corridors for lunar/Mars aerocapture missions

Parametric atmospheric entry corridor data are presented for Earth and Mars aerocapture. Parameter ranges were dictated by the range of mission designs currently envisioned as possibilities for the Human Exploration Initiative (HEI). This data, while not providing a means for exhaustive evaluation of aerocapture performance, should prove to be a useful aid for preliminary mission design and evaluation. Entry corridors are expressed as ranges of allowable vacuum periapse altitude of the planetary approach hyperbolic orbit, with chart provided for conversion to an approximate flight path angle corridor at entry interface (125 km altitude). The corridor boundaries are defined by open-loop aerocapture trajectories which satisfy boundary constraints while utilizing the full aerodynamic control capability of the vehicle (i.e., full lift-up or full lift-down). Parameters examined were limited to those of greatest importance from an aerocapture performance standpoint, including the approach orbit hyperbolic excess velocity, the vehicle lift to drag ratio, maximum aerodynamic load factor limit, and the apoapse of the target orbit. The impact of the atmospheric density bias uncertainties are also included. The corridor data is presented in graphical format, and examples of the utilization of these graphs for mission design and evaluation are included.

Ling, Lisa M.↗

NASA NTRS1992

A probabilistic fracture mechanics approach for structural reliability assessment of space flight systems

A probabilistic fracture mechanics approach for predicting the failure life distribution due to subcritical crack growth is presented. A state-of-the-art crack propagation method is used in a Monte Carlo simulation to generate a distribution of failure lives. The crack growth failure model expresses failure life as a function of stochastic parameters including environment, loads, material properties, geometry, and model specification errors. A stochastic crack growth rate model that considers the uncertainties due to scatter in the data and mode misspecification is proposed. The rationale for choosing a particular type of probability distribution for each stochastic input parameter and for specifying the distribution parameters is presented. The approach is demonstrated through a probabilistic crack growth failure analysis of a welded tube in the Space Shuttle Main Engine. A discussion of the results from this application of the methodology is given.

Sutharshana, S.↗

NASA NTRS2007

Risk Assessment of Bone Fracture During Space Exploration Missions to the Moon and Mars

The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions for specific activities or events during a particular space mission. The key elements of the BFxRM are the mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population and unknowns associated with the effects of the space environment. Consequently, parameter distributions were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission scenarios. The model predicts an increase in the probability of fracture as the mission length increases and fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and resource planning in order to mitigate the risk of bone fracture in space.

Lewandowski, Beth E.↗

NASA NTRS2008

Risk Assessment of Bone Fracture During Space Exploration Missions to the Moon and Mars

The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions for specific activities or events during a particular space mission. The key elements of the BFxRM are the mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population and unknowns associated with the effects of the space environment. Consequently, parameter distributions were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission scenarios. The model predicts an increase in the probability of fracture as the mission length increases and fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and resource planning in order to mitigate the risk of bone fracture in space.

Lewandowski, Beth E.↗

NASA NTRS2012

Recent Progress on Deep Blue Aerosol Algorithm as Applied TO MODIS, SEA WIFS, and VIIRS, and Their Intercomparisons with Ground Based and Other Satellite Measurements

The impact of natural and anthropogenic sources of aerosols has gained increasing attention from scientific communities in recent years. Indeed, tropospheric aerosols not only perturb radiative energy balance by interacting with solar and terrestrial radiation, but also by changing cloud properties and lifetime. Furthermore, these anthropogenic and natural air particles, once generated over the source regions, can be transported out of the boundary layer into the free troposphere and can travel thousands of kilometers across oceans and continents resulting in important biogeochemical impacts on the ecosystem. With the launch of SeaWiFS in 1997, Terra/MODIS in 1999, and Aqua/MODIS in 2002, high quality comprehensive aerosol climatology is becoming feasible for the first time. As a result of these unprecedented data records, studies of the radiative and biogeochemical effects due to tropospheric aerosols are now possible. In this talk, we will demonstrate how this newly available SeaWiFS/MODIS aerosol climatology can provide an important piece of puzzles in reducing the uncertainty of estimated climatic forcing due to aerosols. We will start with the global distribution of aerosol loading and their variabilities over both land and ocean on short- and long-term temporal scales observed over the last decade. The recent progress made in Deep Blue aerosol algorithm on improving accuracy of these Sea WiFS / MODIS aerosol products in particular over land will be discussed. The impacts on quantifying physical and optical processes of aerosols over source regions of adding the Deep Blue products of aerosol properties over bright-reflecting surfaces into Sea WiFS / MODIS as well as VIIRS data suite will also be addressed. We will also show the intercomparison results of SeaWiFS/MODIS retrieved aerosol optical thickness with data from ground based AERONET sunphotometers over land and ocean as well as with other satellite measurements. The trends observed in global aerosol loadings of both natural and anthropogenic sources based upon more than a decade of combined MODIS/SeaWiFS data (1997-2011) will be discussed. We will also address the importance of various key issues such as differences in spatial-temporal sampling rates and observation time between different satellite measurements could potentially impact these intercomparisons results, especially for using the monthly mean data, and thus on estimates of long-term aerosol trends.

Hsu, N. Christina↗

NASA NTRSDate not supplied

Experimental Methods Using Force Application of A Single Boom for A 500-M²-Class Solar Sail

Solar sailing missions rely on deployable systems for large area-to-mass ratios once in space, while still being small enough for launcher envelopes in the stowed configuration. Many of these deployable systems feature booms that are flattened and subsequently coiled onto a spool / hub. As part of a collaborative deployable space structures research effort between the National Aeronautics and Space Administration (NASA) and German Aerospace Center (DLR), a boom deployment system for a future 500 m² solar sail has been developed. To achieve the respective solar sail size goal, 16.5-m-long booms produced by NASA were integrated into a DLR-designed deployer mechanism. This considerable size, as well as the lightweight construction of the booms and respective deployable systems makes testing on the ground a significant challenge. Some systems for gravity compensation as well as vertical testing to minimize the influence of gravity have been used in the past. However, an uncertainty factor remains towards the behaviour in the space environment. The focus of this paper is the load application testing of the integrated boom-deployment mechanism system under microgravity condition, as well as a comparison to vertical testing under gravity. Testing in microgravity was performed during a parabolic flight test inside an aircraft and it included stowage and full deployment of a single boom along the longitudinal axis of the aircraft. The load application was split into two categories: static testing, which induced a linear force ramp to the static boom for either compression or compression-bending; and dynamic testing, which applied a constant force to a boom during extension provided by the deployer mechanism. Both types of tests were performed multiple times at two distinct lengths of the boom, fully deployed (12.76 m) and ~28.6 % deployed (3.65 m). More parameters that are vital to the test philosophy are the angles of attack in the force application and the highest force applied. The data acquisition used for the applied load and deflection measurements of the boom is also presented.

Martin Richter↗

NASA NTRSDate not supplied

Experimental Methods Using Force Application of A Single Boom for A 500-M²-Class Solar Sail

Solar sailing missions rely on deployable systems for large area-to-mass ratios once in space, while still being small enough for launcher envelopes in the stowed configuration. Many of these deployable systems feature booms that are flattened and subsequently coiled onto a spool / hub. As part of a collaborative deployable space structures research effort between the National Aeronautics and Space Administration (NASA) and German Aerospace Center (DLR), a boom deployment system for a future 500 m² solar sail has been developed. To achieve the respective solar sail size goal, 16.5-m-long booms produced by NASA were integrated into a DLR-designed deployer mechanism. This considerable size, as well as the lightweight construction of the booms and respective deployable systems makes testing on the ground a significant challenge. Some systems for gravity compensation as well as vertical testing to minimize the influence of gravity have been used in the past. However, an uncertainty factor remains towards the behaviour in the space environment. The focus of this paper is the load application testing of the integrated boom-deployment mechanism system under microgravity condition, as well as a comparison to vertical testing under gravity. Testing in microgravity was performed during a parabolic flight test inside an aircraft and it included stowage and full deployment of a single boom along the longitudinal axis of the aircraft. The load application was split into two categories: static testing, which induced a linear force ramp to the static boom for either compression or compression-bending; and dynamic testing, which applied a constant force to a boom during extension provided by the deployer mechanism. Both types of tests were performed multiple times at two distinct lengths of the boom, fully deployed (12.76 m) and ~28.6 % deployed (3.65 m). More parameters that are vital to the test philosophy are the angles of attack in the force application and the highest force applied. The data acquisition used for the applied load and deflection measurements of the boom is also presented.

Martin Richter↗

DOE OSTI2023

Impacts of Dispatch Strategies and Forecast Errors on the Economics of Behind-the-Meter PV-Battery Systems

To assess the economic value of batteries in hybrid PV-battery systems, one must create a dispatch profile for the battery. Many analyses of battery value assume perfect forecasts of PV generation and load, determining an upper limit on the value of the battery. Prior work that accounts for forecast uncertainty often does so in the context of a single dispatch algorithm, which does not provide a baseline for comparison. Furthermore, when multiple dispatch algorithms are assessed with uncertainty, the benefits considered are for diesel generation in a microgrid, not retail rate savings. This work addresses the gaps in the literature by comparing the performance of both heuristic and optimal dispatch algorithms for retail rate savings under forecast uncertainty, and provides comparisons of the robustness of these algorithms and their associated estimates of economic value. We find that using a perfect forecast can overestimate the value of hybrid PV-battery systems between 1% and 8% compared to the reality of using a day-ahead forecast, depending on the dispatch algorithm used. Thus, accounting for forecast uncertainty in system design and analysis will significantly improve the accuracy of modeled system values.

batteries↗

DOE OSTI2023

Impacts of Dispatch Strategies and Forecast Errors on the Economics of Behind-the-Meter PV-Battery Systems: Preprint

To assess the economic value of batteries in hybrid PV-battery systems, one must create a dispatch profile for the battery. Many analyses of battery value assume perfect forecasts of PV generation and load, determining an upper limit on the value of the battery. Prior work that accounts for forecast uncertainty often does so in the context of a single dispatch algorithm, which does not provide a baseline for comparison. Furthermore, when multiple dispatch algorithms are assessed with uncertainty, the benefits considered are for diesel generation in a microgrid, not retail rate savings. This work addresses the gaps in the literature by comparing the performance of both heuristic and optimal dispatch algorithms for retail rate savings under forecast uncertainty, and provides comparisons of the robustness of these algorithms and their associated estimates of economic value. We find that using a perfect forecast can overestimate the value of hybrid PV-battery systems between 1% and 8% compared to the reality of using a day-ahead forecast, depending on the dispatch algorithm used. Thus, accounting for forecast uncertainty in system design and analysis will significantly improve the accuracy of modeled system values.

batteries↗

DOE OSTI2020

Summary of Expedited Results from Samples Supporting Tank Closure Cesium Removal (TCCR) Batch 3

Savannah River Remediation (SRR) is currently operating the Tank Closure Cesium Removal (TCCR) process to remove 137 Cs from tank waste supernate using an ion exchange process. As part of that process, Savannah River National Laboratory (SRNL) receives and analyzes samples in support of the qualification of each batch to be processed. SRNL recently received supernate samples retrieved from Tank 10H as well as in-tank batch contact samples for characterization in support of qualifying Batch 3 for processing through the TCCR unit. SRNL received and characterized a set of dip samples collected from Tank 10H (one surface and one variable depth sample). No solids were observed in either sample and therefore the samples were combined, and a suite of analyses were performed on the composite sample. The density of the combined sample was 1.174 g/mL (0.18 %RSD), which was consistent with the measured sodium concentration of 3.54 M (0.08 %RSD). The 137Cs activity in the sample was 4.75E+07 dpm/mL (4.28 %RSD) and the total Cs concentration was calculated to be 1.54 mg/L (1.16E-05 M). The in-tank batch contact samples consist of 0.1 g of crystalline silicotitanate (CST) contained within a teabag device. Duplicate samples were submerged in Tank 10H supernate for a period of about 13 days, after which time they were retrieved and transferred to SRNL for analysis. At SRNL the CST was rinsed to remove excess salt solution and the CST was then air dried before being digested for analysis. Results of the analysis indicated a 137 Cs loading of 3.66E+10 ± 1.40E+09 dpm/g or 16.5 ± 0.63 Ci/kg CST . This value represents a bounding upper limit as it includes the addition of two sigma uncertainty from replicate analysis of the individual teabag samples as well as the addition of the small amount of 137 Cs activity measured in the rinse solutions. The above values are based on the air-dried mass of CST. Correcting to the true dry mass using a F-factor of 0.8191 results in a maximum loading of 20.2 ± 0.77 Ci/kg CST .

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2014

Advanced Modeling and Uncertainty Quantification for Flight Dynamics; Interim Results and Challenges

As part of the NASA Vehicle Systems Safety Technologies (VSST), Assuring Safe and Effective Aircraft Control Under Hazardous Conditions (Technical Challenge #3), an effort is underway within Boeing Research and Technology (BR&T) to address Advanced Modeling and Uncertainty Quantification for Flight Dynamics (VSST1-7). The scope of the effort is to develop and evaluate advanced multidisciplinary flight dynamics modeling techniques, including integrated uncertainties, to facilitate higher fidelity response characterization of current and future aircraft configurations approaching and during loss-of-control conditions. This approach is to incorporate multiple flight dynamics modeling methods for aerodynamics, structures, and propulsion, including experimental, computational, and analytical. Also to be included are techniques for data integration and uncertainty characterization and quantification. This research shall introduce new and updated multidisciplinary modeling and simulation technologies designed to improve the ability to characterize airplane response in off-nominal flight conditions. The research shall also introduce new techniques for uncertainty modeling that will provide a unified database model comprised of multiple sources, as well as an uncertainty bounds database for each data source such that a full vehicle uncertainty analysis is possible even when approaching or beyond Loss of Control boundaries. Methodologies developed as part of this research shall be instrumental in predicting and mitigating loss of control precursors and events directly linked to causal and contributing factors, such as stall, failures, damage, or icing. The tasks will include utilizing the BR&T Water Tunnel to collect static and dynamic data to be compared to the GTM extended WT database, characterizing flight dynamics in off-nominal conditions, developing tools for structural load estimation under dynamic conditions, devising methods for integrating various modeling elements into a real-time simulation capability, generating techniques for uncertainty modeling that draw data from multiple modeling sources, and providing a unified database model that includes nominal plus increments for each flight condition. This paper presents status of testing in the BR&T water tunnel and analysis of the resulting data and efforts to characterize these data using alternative modeling methods. Program challenges and issues are also presented.

Hyde, David C.↗

NASA NTRS2012

Uncertainty Due to Unsteady Fluid/Structure Interaction for the Ares I Vehicle Traversing the Transonic Regime

Rapid reduced-order numerical models are being investigated as candidates to simulate the dynamics of a flexible launch vehicle during atmospheric ascent. There has also been the extension of these new approaches to include gust response. These methods are used to perform aeroelastic and gust response analyses at isolated Mach numbers. Such models require a method to time march through a succession of ascent Mach numbers. An approach is presented for interpolating reduced-order models of the unsteady aerodynamics at successive Mach numbers. The transonic Mach number range is considered here since launch vehicles can suffer the highest dynamic loads through this range. Realistic simulations of the flexible vehicle behavior as it traverses this Mach number range are presented. The response of the vehicle due to gusts is computed. Uncertainties in root mean square and maximum bending moment and crew module accelerations are presented due to assumed probability distributions in design parameters, ascent flight conditions, gusts. The primary focus is on the uncertainty introduced by modeling fidelity. It is found that an unsteady reduced order model produces larger excursions in the root mean square loading and accelerations than does a quasi-steady reduced order model.

Bartels, Robert E.↗

DOE OSTI2024

Solar, Wind, and Load Forecasting Dataset for MISO, NYISO, and SPP Balancing Areas

The Performance-based Energy Resource Feedback, Optimization, and Risk Management (PERFORM) program is an initiative intended to foster "a fundamental shift in grid management rooted in an understanding of asset risk and system risk" (ARPA-E 2020). Launched by the Advanced Research Projects Agency-Energy (ARPA-E), the program supports efforts to incorporate uncertainty in electric power decision making. In support of PERFORM, the National Renewable Energy Laboratory (NREL) has produced a set of time-coincident forecasts of solar, wind, and load profiles. As part of Phase I of the PERFORM effort, NREL created a dataset that consists of one year of time-coincident load, wind, and solar actuals and probabilistic forecasts based on data from the Electric Reliability Council of Texas (ERCOT) (Bryce et al. 2023). In Phase II, NREL developed similar datasets for three other U.S. Independent System Operators (ISO): the Midcontinent Independent System Operator (MISO), the New York Independent System Operator (NYISO), and the Southwest Power Pool (SPP).

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Robust Scheduling of Networked Microgrids for Economics and Resilience Improvement

The benefits of networked microgrids in terms of economics and resilience are investigated and validated in this work. Considering the stochastic unintentional islanding conditions and conventional forecast errors of both renewable generation and loads, a two-stage adaptive robust optimization is proposed to minimize the total operating cost of networked microgrids in the worst scenario of the modeled uncertainties. By coordinating the dispatch of distributed energy resources (DERs) and responsive demand among networked microgrids, the total operating cost is minimized, which includes the start-up and shut-down cost of distributed generators (DGs), the operation and maintenance (O&M) cost of DGs, the cost of buying/selling power from/to the utility grid, the degradation cost of energy storage systems (ESSs), and the cost associated with load shedding. The proposed optimization is solved with the column and constraint generation (C&CG) algorithm. The results of case studies demonstrate the advantages of networked microgrids over independent microgrids in terms of reducing total operating cost and improving the resilience of power supply.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗