Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed generation”

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 775 records · Page 43

Optimal PV Inverter Control in Distribution Systems via Data-Driven Distributionally Robust Optimization

Distribution systems with high penetration of uncertain solar generation call for advanced control strategies of photovoltaics (PVs) inverters. This paper proposes a data-driven distributionally robust optimization (DDDRO) approach to optimally controlling the PV inverters to improve the system operation performance under solar power uncertainties. In the proposed DDDRO approach, a Wasserstein ball-based method is proposed to construct the distributional ambiguity set to model the uncertainties of PV generation through partial observations of historical data without knowing exact probability distributions. We further reformulate the computationally intractable DDDRO model to a mixed integer second order cone programming (MISOCP) problem. The effectiveness and out-of-sample performance of the proposed approach have been demonstrated on a modified IEEE 33-node system. We conduct a comparative study to compare the proposed method with traditional chance constrained programming (CCP). It shows that the proposed DDDRO approach can provide a less conservative yet robust solution to minimize the worse-case expectation of the total network loss while maintaining nodal voltages in a secure range.

Xue, Yaosuo↗

Measurement of double-parton scattering in inclusive production of four jets with low transverse momentum in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A measurement of inclusive four-jet production in proton-proton collisions at a center-of-mass energy of 13 TeV is presented. The transverse momenta of jets within |η| < 4.7 are required to exceed 35, 30, 25, and 20 GeV for the first-, second-, third-, and fourth-leading jet, respectively. Differential cross sections are measured as functions of the jet transverse momentum, jet pseudorapidity, and several other observables that describe the angular correlations between the jets. The measured distributions show sensitivity to different aspects of the underlying event, parton shower modeling, and matrix element calculations. In particular, the interplay between angular correlations caused by parton shower and double-parton scattering contributions is shown to be important. The double-parton scattering contribution is extracted by means of a template fit to the data, using distributions for single-parton scattering obtained from Monte Carlo event generators and a double-parton scattering distribution constructed from inclusive single-jet events in data. The effective double-parton scattering cross section is calculated and discussed in view of previous measurements and of its dependence on the models used to describe the single- parton scattering background.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Distributed Energy Resource-Cognizant Upgrade Paths to the Traditional Restoration Strategy of Utilities for Improved Load Restoration

Climate change has resulted in increasingly impactful and more frequent occurrences of extreme weather events. This trend poses a significant challenge for distribution utilities and system operators to ensure that there is uninterrupted power supply to critical loads in their networks under fault scenarios; however, currently, utilities deploying the automated fault location, isolation and restoration (FLISR) function in their advanced distribution management system (ADMS) do not take into account the available generation and load-modification capabilities of distributed energy resources present in the disconnected network due to an upstream isolated fault. This results in the network reconfiguration and restoration to result in sub-optimal load restoration. Therefore, this paper presents two approaches that can upgrade the existing FLISR capabilities of distribution utilities to significantly increase the restoration of critical loads. The performance of the proposed approaches is evaluated on a numerical model of a real distribution feeder in Georgia, USA.

DER↗

Distributed Energy Resource-Cognizant Upgrade Paths to the Traditional Restoration Strategy of Utilities for Improved Load Restoration: Preprint

Climate change has resulted in increasingly impactful and more frequent occurrences of extreme weather events. This trend poses a significant challenge for distribution utilities and system operators to ensure that there is uninterrupted power supply to critical loads in their networks under fault scenarios; however, currently, utilities deploying the automated fault location, isolation and restoration (FLISR) function in their advanced distribution management system (ADMS) do not take into account the available generation and load-modification capabilities of distributed energy resources present in the disconnected network due to an upstream isolated fault. This results in the network reconfiguration and restoration to result in sub-optimal load restoration. Therefore, this paper presents two approaches that can upgrade the existing FLISR capabilities of distribution utilities to significantly increase the restoration of critical loads. The performance of the proposed approaches is evaluated on a numerical model of a real distribution feeder in Georgia, USA.

DER↗

Enabling Grid-Aware Market Participation of Aggregate Flexible Resources

Increasing integration of distributed energy resources (DERs) within distribution feeders provides unprecedented flexibility at the distribution-transmission interconnection. With the new FERC 2222 order, DER aggregations are allowed to participate in energy market. To enable market participation, these virtual power plants need to provide their generation cost curves. This paper proposes efficient optimization formulations and solution approaches for the characterization of hourly as well as multi-time-step generation cost curves for a distribution system with high penetration of DERs. Network and DER constraints are taken into account when deriving these cost curves, and they enable active distribution systems to bid into the electricity market. The problems of deriving linear and quadratic cost curves are formulated as robust optimization problems and tractable reformulation/solution algorithm are developed to facilitate efficient calculations. The proposed formulations and solution algorithm are validated on a realistic test feeder with high penetration of flexible resources.

aggregated distributed energy resources↗

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗

Impact of High Penetration Distributed Energy Resources on the Bulk Electric System

The growth of distributed energy resources (DERs), mostly generation using intermittent renewable energy sources, along with the retirement of central generation (mostly conventional power plants using fossil-fuel-based resources) has implications for the steady-state and dynamic performance of the bulk electric system (BES). Among the DERs connected to the distribution system, photovoltaic (PV) systems are the most prevalent and have been installed at an increasing rate. Until recently, the penetration levels of DERs have not been high enough to create significant impacts on the reliability and security of power system operation. However, in recent years, the penetration levels of DERs have risen to a level that their impacts on the bulk electric system should be considered in a detailed fashion in planning and operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Review of Behind-the-Meter Solar Generation Modeling and Forecasting

Solar photovoltaic systems largely integrated within the distribution grid are operated 'behind-the-meter' and power generation cannot be directly monitored by most utilities. The increasing penetration of behind-the-meter solar photovoltaic systems can deter efficient network and market operations due to variability and uncertainty in net load, which is exacerbated by limited visibility and the difficulty in analyzing the hosting capacity. Risk introduced by behind-the-meter solar contributions may hinder reliable and secure grid operations due to biased system monitoring and forecasts. Accurate behind-the-meter estimations, together with capacity and specification forecasts, thus play a key role in balancing supply and demand and this article reviews the pertinent literature, identifying key characteristics and predictive methods for efficient behind-the-meter solar photovoltaic generation. Forecasting is central to methods herein. The fundamental characteristics of behind-the-meter solar forecasting, including which methods are applicable for scenario-driven use cases, are driven by the metrics most useful for system-wide performance evaluation. To this aim, the literature is reviewed with a focus on forecasting applications for aggregate, regional behind-the-meter generation useful to bulk system and utility operations. As distinguished from net load forecasting, subtleties in these coincident tasks are explored before concluding with recommendations for current practice and future implementations.

behind-the-meter↗

Predictive and Cooperative Voltage Control with Probabilistic Load and Solar Generation Forecasting

This paper proposes predictive cooperative voltage control method in a power system with high penetration of photovoltaic (PV) units. Cooperative distributed control of the reactive power output of PV inverters is coordinated with operation of voltage regulators (VRs) to maintain system voltages within an appropriate bandwidth. Probabilistic forecasting of the solar power generation and the loads is applied to estimate voltage changes which, in turn, are used to set the VR tap positions for preventing large voltage fluctuations with the lowest risk considering the voltage distribution estimation. The fine tuning of voltage adjustment is achieved by cooperative control of PV inverters to maintain a uniform voltage profile across the system. The proposed method is tested on the modified IEEE 123-node test feeder with high PV penetration using real insolation data and with constant loads replaced by several different load profiles. Simulation results demonstrate the effectiveness of the coordinated approach for voltage control with cooperative PV and predictive VR controls taking into account probabilistic load and solar power forecasts.

Cooperative Control↗

PacWave Grid Integration Study: Transient and Dynamic Conditions (Final Report)

This report describes the results of PSCAD simulations that were performed in 2020 to assess the impacts of PacWave generation on CLPUD's 12.47-kV distribution and 69-kV subtransmission systems. PacWave South (PacWave) is a wave energy test facility planned by Oregon State University. PacWave shore facilities will be located south of Seal Rock on the Oregon coast. PacWave is expected to be operational in 2022 and will connect up to 10 MW of generation to CLPUD's Seal Rock distribution feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Electric field distribution in a non-self-sustained RF discharge with ionization generated by Ns discharge pulses

Abstract Electric field in a capacitively coupled, non-self-sustained RF discharge nitrogen plasma with external ionization generated by high-voltage ns pulses has been measured by ps Electric field induced second harmonic generation (EFISH). The measurements are made both in the bulk of the plasma and in the sheaths, using the absolute calibration by the Laplacian field between two plane electrodes. The results are compared with the kinetic modeling calculations. The RF electric field in the sheaths is significantly higher compared to that in the plasma, due to the displacement of the electrons by the drift oscillations and the resultant plasma self-shielding. However, the kinetic modeling predictions indicate that the electron impact ionization in the sheaths is largely ineffective, due to the low electron density. The reduction of the electric field in the plasma by the self-shielding in the sheaths is moderate, such that the energy is coupled to the plasma by the below-breakdown RF field. The peak RF field in the plasma is in the range of 15–25 Td, indicating the efficient vibrational excitation of N 2 by electron impact. The modeling predictions suggest that the targeted vibrational excitation of molecular species in a non-self-sustained RF discharge with external ionization is scalable to high pressures, electron densities, and discharge powers. The present approach can be extended to the vibrational excitation of other molecular species where the vibrational relaxation is relatively slow, such as CO, CO 2 , and H 2 .

Physics↗

Life Cycle Emissions Factors for Electricity Generation Technologies

This dataset consists of a table containing the distribution of literature estimates of greenhouse gas emissions for the following electricity generation and storage technologies: biopower, coal, concentrating solar power, geothermal, hydrogen storage, hydropower, lithium-ion battery storage, natural gas, nuclear, ocean, oil, photovoltaic, pumped-storage hydropower, and wind. Quartile estimates of life cycle emissions factors in units of grams of carbon dioxide equivalent per kilowatt hour of generation (g CO2e/kWh) are provided for the following life cycle stages: one-time upstream, ongoing combustion, ongoing non-combustion, one-time downstream, and total. Literature estimates were compiled by the LCA Harmonization study and subsequent updates, as detailed in the factsheet which accompanies this dataset, https://www.nlr.gov/docs/fy21osti/80580.pdf .

01 COAL, LIGNITE, AND PEAT↗

Selective Depolymerization for Sculpting Polymethacrylate Molecular Weight Distributions

Chain-end reactivation of polymethacrylates generated by reversible-deactivation radical polymerization (RDRP) has emerged as a powerful tool for triggering depolymerization at significantly milder temperatures than those traditionally employed. In this study, we demonstrate how the facile depolymerization of poly(butyl methacrylate) (PBMA) can be leveraged to selectively skew the molecular weight distribution (MWD) and predictably alter the viscoelastic properties of blended PBMA mixtures. By mixing polymers with thermally active chain ends with polymers of different molecular weights and inactive chain ends, the MWD of the blends can be skewed to be high or low by selective depolymerization. This approach leads to the counterintuitive principle of the “destructive strengthening” of a material. As a result, we demonstrate, as a proof of concept, the encryption of information within polymer mixtures by linking Morse code with the MWDs before and after selective depolymerization, allowing for the encoding of data within blends of synthetic macromolecules.

36 MATERIALS SCIENCE↗

State of Common Grid Services Definitions

This document is prepared as part of the Department of Energy’s Grid Modernization Laboratory Consortium (GMLC) 2.5.2 project, whose goal is to develop and socialize a common set of grid service definitions relevant to grid-related interactions with distributed energy resources (DER: responsive generation, storage, and loads), and to advance the concept and requirements of the Energy Services Interface (ESI) to the point of launching related interface standards and guides that can be implemented in communication protocols and business process definitions. The notion of “grid services” is integral to the definition of an ESI because a key principle of the ESI is that it permits coordination between grid operators and DER facilities in a way that is service-oriented, with an understanding of performance expectations. This document reviews the current state of grid service definitions, including those actively used in the market today as well as new services that have been proposed for future implementation. The document describes grid services used in transmission as well as distribution systems. In defining grid services, this document also distinguishes between two fundamental concepts: an “operational objective” and a “grid service,” which describes a generator’s or customer’s expected physical performance in delivering power to or consuming power from the grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Development of physics-consistent conditional diffusion model to overcome data scarcity in critical heat flux

Deep generative modeling provides a powerful pathway to overcome data scarcity in energy-related applications where experimental data are often limited. By learning the underlying probability distribution of the training dataset, deep generative models, such as the diffusion model, can generate high-fidelity synthetic samples that statistically resemble the training data. Such synthetic data generation can significantly enrich the size and diversity of the available training data, and more importantly, improve the robustness of downstream machine learning models in predictive tasks. The objective of this paper is to investigate the effectiveness of diffusion models for overcoming data scarcity in nuclear energy applications. By leveraging a public dataset on critical heat flux which covers a wide range of commercial nuclear reactor operational conditions, we developed a diffusion model that can generate an arbitrary amount of synthetic samples. Since a vanilla diffusion model can only generate samples randomly, we also developed a conditional diffusion model capable of generating targeted critical heat flux data under user-specified thermal-hydraulic conditions. The performance of the diffusion model was evaluated based on its ability to capture empirical feature distributions and pair-wise correlations, as well as to maintain physical consistency. The results showed that both the diffusion model and conditional diffusion model can successfully generate realistic and physics-consistent critical heat flux data. Furthermore, uncertainty quantification results demonstrate that the conditional diffusion model is highly effective in augmenting critical heat flux data while maintaining acceptable levels of uncertainty.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Cybersecurity Considerations for Grid-Connected Batteries with Hardware Demonstrations

The share of renewable and distributed energy resources (DERs), like wind turbines, solar photovoltaics and grid-connected batteries, interconnected to the electric grid is rapidly increasing due to reduced costs, rising efficiency, and regulatory requirements aimed at incentivizing a lower-carbon electricity system. These distributed energy resources differ from traditional generation in many ways including the use of many smaller devices connected primarily (but not exclusively) to the distribution network, rather than few larger devices connected to the transmission network. DERs being installed today often include modern communication hardware like cellular modems and WiFi connectivity and, in addition, the inverters used to connect these resources to the grid are gaining increasingly complex capabilities, like providing voltage and frequency support or supporting microgrids. To perform these new functions safely, communications to the device and more complex controls are required. The distributed nature of DER devices combined with their network connectivity and complex controls interfaces present a larger potential attack surface for adversaries looking to create instability in power systems. To address this area of concern, the steps of a cyberattack on DERs have been studied, including the security of industrial protocols, the misuse of the DER interface, and the physical impacts. These different steps have not previously been tied together in practice and not specifically studied for grid-connected storage devices. In this work, we focus on grid-connected batteries. We explore the potential impacts of a cyberattack on a battery to power system stability, to the battery hardware, and on economics for various stakeholders. We then use real hardware to demonstrate end-to-end attack paths exist when security features are disabled or misconfigured. Our experimental focus is on control interface security and protocol security, with the initial assumption that an adversary has gained access to the network to which the device is connected. We provide real examples of the effectiveness of certain defenses. This work can be used to help utilities and other grid-connected battery owners and operators evaluate the severity of different threats and the effectiveness of defense strategies so they can effectively deploy and protect grid-connected storage devices.

25 ENERGY STORAGE↗

Weibel instability beyond bi-Maxwellian anisotropy

The shape of the anisotropic velocity distribution function, beyond the realm of strict Maxwellians can play a significant role in determining the evolution of the Weibel instability dictating the dynamics of self-generated magnetic fields. For non-Maxwellian distribution functions, we show that the direction of the maximum growth rate wave vector changes with shape. In this work, we investigate different laser-plasma interaction model distributions which show that their Weibel generated magnetic fields may require closer scrutiny beyond the second moment (temperature) anisotropy ratio characterization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗