SEARCH · Engineering Papers

Results for “Fixed grid”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2025

LLM-Based Adaptive Distribution Voltage Regulation Under Frequent Topology Changes: An In-Context MPC Framework

This paper proposes a large language model (LLM) based adaptive inverter control for distribution voltage regulation under frequent topology changes. We leverage the ability of the LLM to perform in-context learning and create a topology-adaptive surrogate model for power flow calculation. The surrogate model is then integrated with a long short-term memory-based load forecaster and a model predictive control (MPC) scheme to achieve the optimal inverter control that adapts to frequent topology changes. Unlike many existing works that assume fixed-topology grids or require the knowledge of all possible topologies when training a model, the proposed in-context MPC method tackles the distribution voltage control problem under various topologies and adapts to unknown topologies with limited data requirement for fine-tuning. The effectiveness of our method is demonstrated on a modified IEEE 123-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Optimal toolpath generation system and method for additively manufactured composite materials

The present disclosure relates to a system for optimizing an additive manufacturing (AM) process. The system may use a printing component for using a material to form a component in a layer-by-layer process. An electronic controller may control movement of one of the printing component or a substrate on which the component is being formed, in a manner to optimize a toolpath for the printing component as each layer of the component is formed. Optimization operations are performed using an optimization subsystem which enables manufacturability constraints, as well as optimized toolpaths for each layer of the component, to be defined using contours of a level-set function. The level-set function may be used to define the optimized toolpaths within a fixed, predetermined grid, and the optimized toolpaths then used to generate suitable code for controlling movement of the printing component relative to the substrate to create the part or structure in a layer-by-layer process, using the optimized toolpaths.

Lewicki, James↗

DOE OSTI2026

CHESS 2025: Spectrometer orthorectified at-sensor radiance from NEON AOP imaging spectroscopy surveys

This dataset provides Level 1 (L1) orthorectified at-sensor radiance derived from measurements collected by the Imaging Spectrometer-1 (NIS-1) onboard the NEON (National Ecological Observatory Network) Airborne Observation Platform (AOP) for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). NIS-1 captures light reflected from the Earth’s surface in 426 discrete wavelength bands as raw digital numbers (DNs; Level 0). These data are then calibrated to physical units (uW/cm²·sr·nm) following the processing steps described in the NEON Imaging Spectrometer Level 1B Calibrated Radiance Algorithm Theoretical Basis Document (ATBD; Gallery 2022). The data delivered here are the primary inputs for the surface reflectance product in “Custom surface reflectance, shade masks, and equivalent water thickness maps for the Colorado Headwaters Ecological Spectroscopy Study” (Carroll et al. 2026). For intertemporal comparison, the radiance data here are most directly relatable to the v2 radiance data in “NEON AOP Imaging Spectroscopy Survey of Upper East River Colorado Watersheds: Raw-Space Radiance and Observational Variable Dataset” (Goulden et al. 2018), to which the same processing methodology was applied. Together, the radiance and reflectance data enable users to exploit the unique reflection signatures of different surface objects for land cover classification, foliar trait mapping, plant vigor assessment, water content estimation, trace-element identification, and other scientific applications. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. Within each domain, data are delivered by flightline as orthorectified and calibrated hyperspectral rasters in Hierarchical Data Format version 5 (HDF5) format, with radiance values provided in uW/cm²·sr·nm on a fixed, uniform Universal Transverse Mercator (UTM) grid at 1 meter spatial resolution. The radiance rasters include all 426 NIS-1 spectral bands, along with associated quality-assurance (QA) and diagnostic and ancillary layers needed for atmospheric correction workflows. Orthorectified radiance is produced from pushbroom spectrometer observations by applying NEON’s radiometric calibration (including bad pixel masking, dark subtract, dark pedestal shift correction, electronic panel ghost correction, grating ghost correction, deblur correction and flat-fielding) and spectral calibration (using spectral response function band centers and full-width at half-maximum intensity), followed by geolocation and regridding to the fixed grid. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

DOE OSTI2020

Asteroseismic modelling of solar-type stars: a deeper look at the treatment of initial helium abundance

ABSTRACT Detailed understanding of stellar physics is essential towards a robust determination of stellar properties (e.g. radius, mass, and age). Among the vital input physics used in the modelling of solar-type stars which remain poorly constrained, is the initial helium abundance. To this end, when constructing stellar model grids, the initial helium abundance is estimated either (i) by using the semi-empirical helium-to-heavy element enrichment ratio, (ΔY/ΔZ), anchored to the standard big bang nucleosynthesis value, or (ii) by setting the initial helium abundance as a free variable. Adopting 35 low-mass, solar-type stars with multiyear Kepler photometry from the asteroseismic ‘LEGACY’ sample, we explore the systematic uncertainties on the inferred stellar parameters (i.e. radius, mass, and age) arising from the treatment of the initial helium abundance in stellar model grids. The stellar masses and radii derived from grids with free initial helium abundance are lower compared to those from grids based on a fixed ΔY/ΔZ ratio. We find the systematic uncertainties on mean density, radius, mass, and age arising from grids which employ a fixed value of ΔY/ΔZ and those with free initial helium abundance to be ∼ 0.9 per cent, ∼ 2 per cent, ∼ 5 per cent, and ∼ 29 per cent, respectively. We report that the systematic uncertainties on the inferred masses and radii arising from the treatment of initial helium abundance in stellar grids lie within the expected accuracy limits of ESA’s PLATO, although this is not the case for the age.

Nsamba, Benard↗

DOE OSTI2022

Stringent σ8 constraints from small-scale galaxy clustering using a hybrid MCMC + emulator framework

ABSTRACT We present a novel simulation-based hybrid emulator approach that maximally derives cosmological and Halo Occupation Distribution (HOD) information from non-linear galaxy clustering, with sufficient precision for DESI Year 1 (Y1) analysis. Our hybrid approach first samples the HOD space on a fixed cosmological simulation grid to constrain the high-likelihood region of cosmology + HOD parameter space, and then constructs the emulator within this constrained region. This approach significantly reduces the parameter volume emulated over, thus achieving much smaller emulator errors with fixed number of training points. We demonstrate that this combined with state-of-the-art simulations result in tight emulator errors comparable to expected DESI Y1 LRG sample variance. We leverage the new abacussummit simulations and apply our hybrid approach to CMASS non-linear galaxy clustering data. We infer constraints on σ8 = 0.762 ± 0.024 and fσ8(zeff = 0.52) = 0.444 ± 0.016, the tightest among contemporary galaxy clustering studies. We also demonstrate that our fσ8 constraint is robust against secondary biases and other HOD model choices, a critical first step towards showcasing the robust cosmology information accessible in non-linear scales. We speculate that the additional statistical power of DESI Y1 should tighten the growth rate constraints by at least another 50–60 ${{\ \rm per\ cent}}$, significantly elucidating any potential tension with Planck. We also address the ‘lensing is low’ tension, which we find to be in the same direction as a potential tension in fσ8. We show that the combined effect of a lower fσ8 and environment-based bias accounts for approximately $50{{\ \rm per\ cent}}$ of the discrepancy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Laser diagnostics to characterize the in-flame growth of platinum nanoparticles manufactured by the reactive spray deposition technology

Reactive Spray Deposition Technology (RSDT) is an atmospheric pressure, flame-based advanced manufacturing method used to fabricate Membrane Electrode Assemblies (MEAs) for proton exchange fuel cells and water electrolyzers. RSDT combines the flame synthesis of catalyst nanoparticles and their deposition onto the membrane in one step. The properties of the synthesized nanoparticles, such as their Size Distribution Function (SDF), determine the manufactured electrodes’ performance, which can be evaluated when the deposition is complete via ex-situ characterization of the products. Efforts to improve RSDT for manufacturing state-of-the-art MEAs can be significantly enhanced with integrated laser diagnostics that enable in situ measurement of the synthesized catalyst nanoparticles. This paper reports the implementation of laser diagnostics in an RSDT facility and evaluates their potential to assist the manufacturing process. Laser Light Scattering (LS) and Laser-Induced Incandescence (LII) measurements are performed in the oxygen-rich zone of the flame at various distances from the flame fuel jet nozzle downstream of its luminous region. The measurements quantitatively track the evolution in size and volume fraction of platinum-based nanoparticles in two flames that yield different catalytic properties in the manufactured electrodes. The profiles of the measured nanoparticles’ volume fraction along the flame axis can be estimated a priori so that the average nanoparticle sizes can be measured in quasi-real time via LS. Nanoparticles experience an extremely slow growth rate while being convected at distances from 150 nm to 300 nm from the fuel nozzle. Complementary characterizations of the synthesized nanoparticles are performed ex-situ via High-Angle Annular Darkfield Scanning Transmission Electron Microscopy (HAADF-STEM) image analyses of samples collected on grids at a fixed distance from the fuel jet nozzle. Comparing the results from laser and microscopy techniques not only cross-validates the findings but also provides the parameters to infer the bimodal SDF in-situ and yields evidence that the coagulation efficiency of the synthesized nanoparticles has extremely low values in the investigated zone of the flames.

Rayleigh Light Scattering (LS)↗

DOE OSTI2021

Long-Run Marginal CO 2 e Emission Rates for End-Use Electricity Consumption in the State of Washington [Slides]

This analysis contains an estimate of the long-run marginal emission rate for the electric sector in the state of Washington. The long-run marginal emission rate is an estimate of the rate of emissions that would be either induced or avoided by a long-term (i.e., more than several years) change in electrical demand. The metric explicitly takes into account both the underlying evolution of the electric grid, as well as the potential for an incremental change in electrical demand to influence the structural evolution of the grid (i.e., the building and retiring of capital assets, such as generators and transmission lines). It is therefore distinct from the more-commonly-known short-run marginal, which also identifies a marginal emission rate but treats the grid assets as fixed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Long-run Marginal CO2 Emission Rates Workbooks for 2020 Standard Scenarios Cambium Data

This dataset has been superseded by a new set of workbooks that can be found here: https://data.nlr.gov/submissions/183 These workbooks contain modeled estimates of long-run marginal CO2 emission rates (LRMER) for the contiguous United States. The LRMER is an estimate of the rate of emissions that would be either induced or avoided by a long-term (i.e., more than several years) change in electrical demand. It incorporates both the projected changes to the electric grid, as well as the potential for an incremental change in electrical demand to influence the structural evolution of the grid (i.e., the building and retiring of capital assets, such as generators and transmission lines). It is therefore distinct from the more-commonly-known short-run marginal, which treats grid assets as fixed. In addition to year-over-year data, the Levelized LRMER worksheet within each workbook is set up to produce a levelized long-run marginal emission rate based on user-provided inputs. These levelized LRMER values are intended for analysts to use when estimating the emissions induced (or avoided) by a long-term change in end-use electricity demand. Three future scenarios are provided in separate workbooks: A Mid-case (i.e., business-as-usual), and two scenarios with relatively higher or lower renewable energy costs. For more details on these scenarios, see the Standard Scenarios 2020 Report ( https://www.nlr.gov/docs/fy21osti/77442.pdf ). For more data underlying each scenario, see the Standard Scenarios 2020 project (Cambium data) at https://cambium.nlr.gov/ . This data was produced as part of the Cambium project. For more details about the methodology, see the Cambium Documentation: Version 2020 ( https://www.nlr.gov/docs/fy21osti/78239.pdf ). This data is planned to be updated annually. Information on the latest versions can be found at https://www.nlr.gov/analysis/cambium.html .

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Long-run Marginal Emission Rates for Electricity - Workbooks for 2021 Cambium Data

These workbooks contain modeled estimates of long-run marginal emission rates (LRMER) for the contiguous United States. The LRMER is an estimate of the rate of emissions that would be either induced or avoided by a long-term (i.e., more than several years) change in electrical demand. It incorporates both the projected changes to the electric grid, as well as the potential for an incremental change in electrical demand to influence the structural evolution of the grid (i.e., the building and retiring of capital assets, such as generators and transmission lines). It is therefore distinct from the more-commonly-known short-run marginal, which treats grid assets as fixed. The Levelized LRMER worksheet within each workbook is set up to produce a levelized long-run marginal emission rate based on user-provided inputs. These levelized LRMER values are intended for analysts to use when estimating the emissions induced (or avoided) by a long-term change in end-use electricity demand. There are two workbooks that supply the data at two different geographic resolutions: states and GEA regions (20 regions that are similar to, but not exactly the same as, the US EPA's eGRID regions). For more data underlying these emissions factors, see the Cambium 2021 project at https://cambium.nrel.gov/. For more details on the inputs into the scenarios available in the workbooks, see the Standard Scenarios 2021 Report (https://www.nrel.gov/docs/fy22osti/80641.pdf). This data was produced as part of the Cambium project. For more details about the methodology, see the Cambium Documentation: Version 2021 (https://www.nrel.gov/docs/fy22osti/81611.pdf). This data is planned to be updated annually. Information on the latest versions can be found at https://www.nrel.gov/analysis/cambium.html.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Long-run Marginal Emission Rates for Electricity - Workbooks for 2022 Cambium Data

These workbooks contain modeled estimates of long-run marginal emission rates (LRMER) for the contiguous United States. A LRMER is an estimate of the rate of emissions that would be either induced or avoided by a change in electric demand, taking into account how the change could influence both the operation as well as the structure of the grid (i.e., the building and retiring of capital assets, such as generators and transmission lines). It is therefore distinct from the more-commonly-known short-run marginal, which treat grid assets as fixed. Long-run marginal emissions rates are generally appropriate to use when trying to comprehensively estimate the impact of a long-lived (i.e., more than several years) intervention. There are two workbooks that supply the data at two different geographic resolutions: states and GEA regions (20 regions that are similar to, but not exactly the same as, the US EPA's eGRID regions). For more data underlying these emissions factors, see the Cambium 2022 project at https://scenarioviewer.nrel.gov/. For more details on input assumptions and methodology see the associated report (Cambium 2022 Scenario Descriptions and Documentation, https://www.nrel.gov/docs/fy23osti/84916.pdf). This data is planned to be updated annually. Information on the latest versions can be found at https://www.nrel.gov/analysis/cambium.html.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2021

A Kinetic Monte Carlo Approach for Simulating Cascading Transmission Line Failure

In this work, cascading transmission line failures are studied through a dynamical model of the power system operating under fixed conditions. The power grid is modeled as a stochastic dynamical system where first-principles electromechanical dynamics are excited by small Gaussian disturbances in demand and generation around a specified operating point. In this context, a single line failure is interpreted in a large deviation context as a first escape event across a surface in phase space defined by line security constraints. The resulting system of stochastic differential equations admits a transverse decomposition of the drift, which leads to considerable simplification in evaluating the quasipotential (rate function) and, consequently, computation of exit rates. Tractable expressions for the rate of transmission line failure in a restricted network are derived from large deviation theory arguments and validated against numerical simulations. Extensions to realistic settings are considered, and individual line failure models are aggregated into a Markov model of cascading failure inspired by chemical kinetics. Cascades are generated by traversing a graph composed of weighted edges representing transitions to degraded network topologies. Numerical results indicate that the Markov model can produce cascades with qualitative power-law properties similar to those observed in empirical cascades.

Kinetic Monte Carlo, rare events, Cascading Failur↗

DOE OSTI2025

Designing Yield Grids for Airborne Nuclear Explosions to Maintain Decrements in Rayleigh Waveform Amplitudes

This document details a method for selecting grid points that sample the yields of airborne, low yield nuclear explosions that constrains the relative size of seismic waveforms sourced by explosions that are parameterized by this grid. In detail, we impose this constraint so that the vertical component of ground displacement that is sourced by explosions with a fixed height of burst (HoB) at grid point k change by a fixed factor of a when compared to the waveform amplitudes sourced by explosions at grid point k + 1, as observed at a fixed standoff distance. This gridding thereby preserves amplitude ratios along the grid so that they are agnostic to k. We assert that this method is useful in tests that apply algorithms against synthetic or augmented datasets for which the performance of the algorithm must be measured against source size, and waveform amplitudes must change uniformly (e.g., Berg et al, 2025).

58 GEOSCIENCES↗

DOE OSTI2021

Method and system for determining a geographical location of an unmanned aerial vehicle flown in proximity to a power grid

A method for determining geolocation of a UAV near a power grid includes detecting, via a transceiver, a carrier signal transmitted from a first grid node to identify the node's fixed geolocation. A response signal may be transmitted from a second grid node in response to the carrier signal to identify a fixed geolocation of the second grid node, or the UAV may process the carrier signal. A processor determines time-of-flight of the carrier signal, e.g., using the response signal, and derives the UAV's geolocation using the time-of-flight. Determining time-of-flight may include referencing a lookup table indexed by time-of-arrival at the transceiver of the modulated carrier and response signals. A timestamp may indicate time-of-transmission of the carrier and response signals, respectively. Deriving geolocation may include subtracting time-of-transmission of the response signal from that of the carrier signal. A system includes the processor and transceiver.

Moore, Andrew J.↗

DOE OSTI2021

Interpolated Sounding and Gridded Sounding Value-Added Products

Standard Atmospheric Radiation Measurement (ARM) Climate Research Facility sounding files provide atmospheric state data in one dimension of increasing time and height per sonde launch. Many applications require a quick estimate of the atmospheric state at higher time resolution. The INTERPOLATEDSONDE (i.e., Interpolated Sounding) Value-Added Product (VAP) transforms sounding data into continuous daily files on a fixed time-height grid, at 1-minute time resolution, on 332 levels, from the surface up to a limit of approximately 40 km. The grid extends that high so the full height of soundings can be captured; however, most soundings terminate at an altitude between 25 and 30 km, above which no data is provided. Between soundings, the VAP linearly interpolates atmospheric state variables in time for each height level. In addition, INTERPOLATEDSONDE provides relative humidity scaled to microwave radiometer (MWR) observations.The INTERPOLATEDSONDE VAP, a continuous time-height grid of relative humidity-corrected sounding data, is intended to provide input to higher-order products, such as the Merged Soundings (MERGESONDE; Troyan 2012) VAP, which extends INTERPOLATEDSONDE by incorporating model data. The INTERPOLATEDSONDE VAP also is used to correct gaseous attenuation of radar reflectivity in products such as the KAZRCOR VAP.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Unleashing the Frequency: Multi-Megawatt Demonstration of 100% Renewable Power Systems with Decentralized Communication-Less Control Scheme

Power systems, which range in size from small microgrids to island systems to large regional grids, are typically managed by a central controller that requires complex communication methods and can be unreliable and pose cyber security risks in certain applications, especially when controlling a larger number of nodes. We propose an inherently robust, scalable method of integration using multiple energy storage systems and distributed energy resources, which does not require any means of dedicated communication. This method moves beyond the paradigm of controlling grid frequency at a fixed value (e.g., 60 Hz), instead allowing the frequency to fluctuate within certain limits (e.g., 59.6-60.4 Hz). With a greater operating range, the frequency can carry necessary information from energy storage systems to highly variable distributed energy resources like photovoltaics, wind, hydro, etc.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Improved Understanding of Multicentury Greenland Ice Sheet Response to Strong Warming in the Coupled CESM2‐CISM2 With Regional Grid Refinement

The simulation of ice sheet‐climate interactions, such as surface mass balance fluxes, is sensitive to model grid resolution. Here we simulate the multi‐century evolution of the Greenland Ice Sheet (GrIS) and its interaction with the climate using the Community Earth System Model version 2.2 (CESM2.2) including an interactive GrIS component (the Community Ice Sheet Model v2.1 [CISM2.1]) under an idealized warming scenario (atmospheric CO 2 increases by 1% yr -1 until quadrupling the pre‐industrial level and then is held fixed). A variable‐resolution (VR) grid with 1/4° regional refinement over the broader Arctic and 1° resolution elsewhere is applied to the atmosphere and land components, and the results are compared with conventional 1° lat‐lon grid simulations to investigate the impact of grid refinement. Compared with the 1° runs, the VR run features a slower rate of surface melt, especially over the western and northern GrIS, where the ice surface slopes gently toward the periphery. This difference pattern originates primarily from higher snow albedo and, thus, weaker albedo feedback in the VR run. The VR grid better captures the CISM ice sheet topography by reducing elevation discrepancies between CAM and CISM and is, therefore, less reliant on the downscaling algorithm, which is known to underestimate albedo gradients. The sea level rise contribution from the GrIS in the VR run is 53 mm by year 150 and 831 mm by year 350, approximately 40% and 20% less than that of the 1° runs, respectively.

Earth System Model↗

DOE OSTI2021

Adaptive critic design-based reinforcement learning approach in controlling virtual inertia-based grid-connected inverters

In this report, an adaptive critic design (ACD) approach is proposed to control the phase and voltage of a grid-connected virtual synchronous generator (VSG). The penetration of fast responding inertia-less power converters significantly affect the stability of the power system, especially weak systems such as micro grids. The concept of virtual inertia addresses this concern by virtually emulating the behavior of a synchronous generator. However, the conventional VSG is designed based on two conditions: (i) fixed operating point and (ii) inductive grid connections. The performance of VSGs in low-voltage semi-resistive microgrids is far from optimal. To overcome the aforementioned concerns, a heuristic dynamic programing (HDP) approach is proposed to optimally control grid-connected VSGs. The neural-network-based inherence of the HDP enables the proposed technique to adapt to any impedance angle. The HDP controller includes two subnetworks: (i) the action network that controls the system optimally and (ii) the critic network, which evaluates the effectiveness of the action network. The simulation and experimental results are provided to evaluate the effectiveness of the proposed technique. As shown, the HDP-based approach illustrates a better performance in comparison with the conventional PI-based VSG in various operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Tri-Level Linear Programming Model for Automatic Load Shedding Using Spectral Clustering

Traditional load shedding schemes can be inadequate in grids with high renewable penetration, leading to unstable events and unnecessary grid islanding. Although for both manual and automatic operating modes load shedding areas have been predefined by grid operators, they have remained fixed, and may be sub-optimal due to dynamic operating conditions. In this work, a distributed tri-level linear programming model for automatic load shedding to avoid system islanding is presented. Preventing islanding is preferred because it reduces the need for additional load shedding besides the disconnection of transmission lines between islands. This is crucial as maintaining the local generation-demand balance is necessary to preserve frequency stability. Furthermore, uneven distribution of generation resources among islands can lead to increased load shedding, causing economic and reliability challenges. This issue is further compounded in modern power systems heavily dependent on non-dispatchable resources like wind and solar. The upper-level model uses complex power flow measurements to determine the system areas to shed load depending on actual operating conditions using a spectral clustering approach. The mid-level model estimates the area system state, while the lower-level model determines the locations and load values to be shed. The solution is practical and promising for real-world applications.

Baquedano-Aguilar, Mario D.↗