Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid Connected”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 883 records · Page 49

SOC Synthetic Microstructure Bank

QUICK START: Start with property_library.html (can be found by typing the filename into the query box) and use the interactive table to filter, sort, and select a microstructure with the desired properties. Search for the alphabetic code to obtain the corresponding dataset. Full description: This is a bank of 1,970 unique 3-phase electrode microstructure files. When you account for reassigning phase IDs (e.g. declare that 1=Ni and 2=pore, instead of 1=pore and 2=Ni), it actually represents 5,910 unique electrode microstructures, each of which could be considered to be either an air or a fuel electrode (e.g. declare that the phase IDs correspond to pore, Ni, and YSZ; or that they correspond to pore, LSCF, and GDC; or whatever electron-conductor and ion-conductor combination is being studied). The voxel size is 50 nm and each electrode file contains a 4x4 grid of (12.5 micron)^3 sub-volumes. If placed together in a grid, they comprise a 50x50x12.5 micron electrode (note that the interfaces between sub-volumes will be sharp; this can be mitigated via simulating annealing/relaxation). The sub-volumes can also be used individually for a reasonably sized 12.5 micron cubic region-of-interest. A user can simply consider the voxel size to be a different value to rescale the volumes (and all of their morphological features, including particle size) as desired. These microstructures were generated using DREAM3D. The general procedure is outlined in https://doi.org/10.1016/j.jpowsour.2018.03.025 The file names are an alphabetic code having to do with the input parameters used in DREAM3D when they were generated. Most users would be best served by starting with the file property_library.html or property_library_subvols.html (which lists properties for each individual subvolume). These files contain a catalogue of the actual, measured properties of every microstructure in the database. Any combination of property values can be filtered and sorted until a desired electrode is found, at which point the user can find the file corresponding to that alphabetic code. The properties in the catalogue include connected TPB density, and for each phase: phase fraction, average particle size, polydispersity of particle size, tortuosity, and connected pair-wise interfacial area. They also include what fraction of each property is connected through to the interfaces of the volume. If the desired combination of properties is not found at first, remember that the phase IDs can be re-assigned arbitrarily, e.g. swapping 1s and 2s. In fact, the database was generated with this in mind so as not to generate redundant microstructures. If the database does not contain the desired property combinations, try to search for the other possible permutations of those properties with re-assigned phase IDs. Please cite https://doi.org/10.1149/10301.0909ecst for use. Please contact the maintainer, William K. Epting, for additional information or assistance.

3D microstructure↗

Achieving Cyber-Resilience for Power Systems using a Learning, Model-Assisted Blockchain Framework

The secure integration and management of distributed energy resources (DER) and power aggregators in the electric grid requires secure communications and a physics-aware Command and Control (C2) strategy. A Blockchain (BC)-based overlay network was developed to provide a security layer for the existing power grid network that mitigates risks in current and legacy network and C2 protocols. By integrating a Model-Assisted Machine Learning (MAML) framework with a Secure Blockchain Overlay Network (SBON) a defense-in-depth strategy was achieved. In our approach, the MAML framework leveraged a smart contract framework to gather network data and learn the dynamics of DER to develop detection strategies for attacks targeting sensors and actuators used by DER. The MAML framework learned dynamical systems models for individual DERs to detect sensor attacks. For DER we utilized a Digital Twin (DT) to accelerate the learning process for a model resistant to stealthy attacks. The project created DT for PV inverters and BESS. The DTs were coupled with a model-assisted, data-driven learning of DER behavior. Specifically, we evaluated architectures for model-based learning with model-free fine-tuning. Additionally, differential privacy techniques were used to obfuscate data, while still allowing the computation of attack detection results based on obfuscated data. The SBON developed leverages a private permissioned blockchain network orchestrated with the Hyperledger Fabric framework. To connect the cyber world, which orchestrates the blockchain fabric, and the physical world where the power network resides, we developed a system implementation to enable the secure interaction of the physical world and the abstracted blockchain.

97 MATHEMATICS AND COMPUTING↗

ARM shortwave spectrometers to study the clear-cloud transition zone and mixing processes

The proposed research is a collaborative effort between NASA/Goddard Space Flight Center and the Hebrew University of Jerusalem. While the NASA team focused on analyzing ground-based hyper-spectral radiance observations to understand cloud edge properties and their connection to mixing processes, the Hebrew University team tackled the problem through cloud modeling activities. By approximating the shortwave spectra in the cloud-clear transition zone as a linear combination of purely clear and purely cloudy spectra we can characterize the variations of cloud optical thickness and cloud droplet effective radius in the transition zone. When applying this method to the measurements of a ground-based shortwave spectroradiometer at the ARM’s SGP site, representing continental conditions, and MAGIC field campaign between Log Angeles, California and Honolulu, Hawaii, representing maritime scenarios, we found that cloud optical depth consistently decreases in both cases, but droplet size decreases much more substantially for the continental regime, suggesting different mixing processes for the continental and maritime conditions. The investigation and measurements of radiation clouds were coupled with a unique cloud modeling. A novel spectral bin microphysics was developed and implemented to the System of Atmospheric Modeling (SAM). In order to resolve cloud transition zones with high spatial gradients of microphysical variables a unique high resolution (10 m) was used in simulations. The model calculates droplet size distributions in each grid point. The model output was transferred to the NASA/GSFC team for utilization in radiative calculations and testing of both radiative algorithm and model representation.

54 ENVIRONMENTAL SCIENCES↗

Blockchain for Fault-Tolerant Grid Operations Version 2.0

This report explores the potential of distributed ledger technology (DLT) as a transformative tool to enhance fault-tolerant operations in electrical distribution systems. Leveraging DLT's core attributes, including an immutable decentralized ledger, distributed consensus mechanisms, and state replication capabilities, this study focuses on three critical use cases. A central aspect of this research centers on the utilization of a consensus-driven ledger, providing actors within the system, such as distributed resources, with access to a reliable data repository. This empowers these actors to collaborate effectively and make informed decisions, all securely recorded on the blockchain. The first use case concentrates on data configuration, utilizing mathematical criteria---particularly, the chi-squared test for gross error detection---to identify trustworthy sensors for advanced decision-making. Building upon this foundation of trust, the second use case, topology identification, accurately determines circuit breaker states, unveiling the distribution network's topology. Ultimately, the third use case leverages this trust to execute switching actions, reconfiguring feeders and restoring power to disconnected customers after fault events. The concept of trust serves as a cornerstone in this approach, marking a departure from traditional fault location, isolation, and service restoration (FLISR) methods. Additionally, the blockchain-based architecture introduces decentralization, empowering disconnected areas to make autonomous decisions, even when communication with a central control center is disrupted. The primary contributions of this report are twofold: (1) a novel approach for evaluating distribution system voltage areas while preserving data ownership and (2) the implementation of interactions between distribution network areas using the actor model. Unlike the previous sequential approach for evaluating the area connection voltages, which required a radial network topology, this study's area model reduction enables a more versatile approach. The area model reduction addresses issues of prolonged data waiting times and multiple points of failure within the previous approach. Notably, the presented evaluation for the reduced network model area connection reveals a significant increase in the differences in voltage magnitudes. Simulation and evaluation of area agents across four distinct cases elucidate the area-level interaction behavior during a fault event. Simulations demonstrate that the proposed distributed FLISR (DFLISR) approach can successfully restore service to an affected area. Varying message delays and message loss probabilities in each simulation case underscore their impacts on restoration times, ranging from 3 min and 32 s to 6 min and 19 s. In contrast, power is not restored in an area in one of our simulation cases.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Upwind Moisture Controls on Interannual Variations of Precipitation and Vegetation in China's Drylands

Dryland precipitation depends on upwind and local moisture sources via moisture recycling. How upwind moisture variations affect interannual variations of downwind precipitation and vegetation in China's drylands remains unclear. We used high-resolution moisture tracking data sets and found terrestrial moisture (93%) was the dominant moisture source for China's drylands, especially from drylands themselves (46%). In most dryland grids, we observed strong correlations between precipitation and upwind moisture sources from 2003 to 2022 (median r = 0.55), with a more significant effect in drier areas. These demonstrated the upwind moisture control on interannual variations of dryland precipitation, in which internal moisture from drylands exceeds the influence of external terrestrial sources. The upwind moisture variations, especially the recycled moisture of drylands, propagate to influence downwind vegetation greenness in precipitation-sensitive dryland areas. Our findings revealed that upwind moisture variations induced by climate or land-cover changes have important implications for water and food security in China's drylands.

54 ENVIRONMENTAL SCIENCES↗

Assessing Customer Experience and Business Models around Price-to-Device Communication and Smart Control Pathways in CalFlexHub

California is facing three major challenges in electrical grid operation: renewable overgeneration, steep evening ramping, and growing peak demand. The state has identified dynamic retail price response as a key strategy evidenced by CPUC’s Dynamic Rates proceeding and CEC’s Load Management Standards. Furthermore, the CEC launched a $16M “California Load Flexibility Research and Deployment Hub (CalFlexHub)” administered by Berkeley Lab to accelerate price-response flexible load technologies in buildings and EV charging. There are more than 16 laboratory and field demonstration projects in CalFlexHub, each demonstrating innovative automated price-response technologies. CalFlexHub tests various pathways through which hourly price signals and triggered control commands are communicated to load-flexible devices such as smart thermostats, heat pumps, water heaters, and EVs. We identified seven unique communication and control pathways, which involve combinations of third-party cloud, device OEM’s cloud, building central gateway, and local controller in between the price server and the load-flexible devices. It is important for utilities and policy makers to understand the long-term implications of each pathway in designing future programs and creating related policies and mandates for market transformation. We propose an evaluation framework including the following aspects: ● Functionality: connectivity and uptime, resilience, and optimization; ● Customer experience: simplicity in setup, troubleshooting support, continuity, customer choice, first cost, and ongoing cost; ● Business model and scalability: advance interoperability, holistic solution, bridge unique gap, customer base, and value streams and pricing structures. In this paper, we identify emerging business models associated with each communication pathway and discuss their positive features and challenges from the above aspects.

Liu, Jingjing↗

Molecular-Resolution Electron Imaging of Defects and Dynamics at the Ice-Water Interface

Water crystallization into hexagonal ice (type I h ) is one of the most critical processes relevant to the Earth’s environment and human activities. However, despite recent breakthroughs in imaging non-equilibrium condensed ice structures, the ice-water interface has never been imaged at a molecular resolution. This is primarily due to the low stability of the hydrogen bonds in ice under high-resolution microscopy conditions and a lack of methods to prepare compatible samples. Here, this presentation describes the first molecular-resolution imaging of ice crystallized from liquid water and the ice-water interface using high-resolution transmission electron microscopy (HRTEM). By encapsulating deionized (DI) water between two amorphous carbon (a-C) TEM grids and subsequently freezing it with liquid N 2 on a cryo sample stage, we generated two types of ice: non-equilibrium, condensed ice from the atmosphere and encapsulated ice from the DI water (Fig. 1A). Condensed ice usually shows irregular, spherulitic shapes (Fig. 1B). Selected area electron diffraction (SAED) shows that they are a mixture of cubic and hexagonal crystals (Fig. 1C). On the contrary, encapsulated ice forms thin films that contain large-area single-crystalline regions of hexagonal ice oriented along the [0001] zone axis (Fig. 1D). Differential electron energy-loss spectroscopy (EELS) confirmed the high purity of the encapsulated ice free from organic contaminations that are common in other encapsulation methods for HRTEM such as graphene liquid cells. These single-crystalline areas are robust under the electron beam up to ~100 e/Å 2 s. Aberration-corrected HRTEM imaging in these areas achieved a line resolution of ~1.3 Å (Fig. 1E and F). This platform allows us to study near-equilibrium ice structures and dynamics at an unprecedented spatial resolution (Fig. 2). For example, we discovered subdomain-rich regions near the defective crystal edges despite the structure appearing single-crystalline according to diffraction criteria. These subdomains connect via low-angle grain boundaries with flat energy landscapes as a function of tilt angles (according to simulations), showing the high tolerance of ice to defect structures. When we tuned the sample temperature and electron flux rate, we observed radiolysis-controlled bubble generation and dissolution in ice single crystals near a steady state of bubble dynamics. Furthermore, rich beam-induced melting and recrystallization dynamics were observed at the ice-water interface with lattice resolution. These data represent the first observation of the ice-water phase transformation at the sub-nanometer level. In summary, the methods developed in this work enabled molecular-resolution observations of ice and the ice-water interface and shed light on the microstructures and phase transformation pathways. Finer control on the temperature, electron irradiation profile, and imaging detector could eventually lead to real-time observation of ice nucleation in water and address long-standing questions in the nucleation pathways.

74 ATOMIC AND MOLECULAR PHYSICS↗

Safe Operations at Roadway Junctions: Intelligent Roadway Infrastructure as Functional Interlocking

Automated vehicle (AV) technology is quickly maturing, and the corresponding infrastructure systems that evaluate traffic and communicate to vehicles requires sophisticated sensing and perception technologies, referred to as intelligent roadway infrastructure (IRI), to complement emerging AV capabilities. IRI provides signals to vehicles, indicating right-of-way for vehicles and communicating to approaching AVs that no other vehicle is failing to yield. This capability, denoted as safety-affirmative signaling, provides a green light or a green arrow as appropriate and affirms through communication links to connected vehicles when it is safe to proceed. About 36% of collisions occur at intersections, with most occurring upon left turns (22.2%) or crossing over (12.6%), and only a small percentage (1.2%) while turning right at an intersection. Of all intersection crashes about half (52.5%) of those vehicles were traveling through a signalized intersection 2. Safety-affirmative signaling would guarantee safety of AV fleet vehicles, by providing the interlocking principle, a term from automated train control that only allows progression through a railway intersection after affirming no opportunity for a crash exists. IRI through safety-affirmative signaling would bring performance and safety to complex roadway intersections where AV transit fleet service is most needed, as well as safety benefits to traditional, non-automated vehicles and vulnerable road users. The implementation of IRI has functional, programmatic, and technical challenges. Research work performed at the National Renewable Energy Laboratory (NREL) in an integrative approach encapsulating these themes, and termed infrastructure perception and control (IPC) is motivated by improved performance (travel time), improved safety (reduced collisions), and improved energy efficiency (less fuel burned and minimized production of greenhouse gases). IPC is intended not only for roadway and intersection applications but also in extension to inform complementary buildings and grid systems to enable better co-management, as vehicles and their charging needs become increasingly integrated into the built environment. The NREL IPC project presents an open-source framework, architecture, and supporting technology to implement IRI, addressing critical issues such as fusion of data, reliability, standardization of data interfaces, and confidence of detection. The framework is informed by previous experience in U.S. Department of Defense research technology, specifically in the use of radar to detect, identify, and track aerial threats. These principles combined with multi-sensor fusion provides for a complete digital twin with known and measurable confidence and accuracy from which safety-affirmative signaling can be developed and deployed.

ADVANCED PROPULSION SYSTEMS,MATHEMATICS AND COMPUT↗

Transmission Value in 2023: Market Data Shows the Value of Transmission Remained High in Certain Locations Despite Overall Low Wholesale Electricity Prices

In 2023 additional electricity transmission would have provided the most value for links that crossed between grid interconnection regions in the United States (the Western Interconnection, the Eastern Interconnection, the Texas Interconnection) or crossed between system operator regions within the same interconnection. Many multi-interconnection or multi-region links had values of greater than $\$20$/MWh, or up to $\$175$ million/yr per 1 GW expanded transmission (subject to limits to the depth of the market at each side of the link). In contrast, many links within regions, or between regions in the northeast, had relatively low values in 2023, following the overall decline in wholesale electricity prices in 2023 compared with 2021-2022. The most valuable link in 2023, at $\$61$ /MWh, was between Texas and the Southwest. Multiple events in 2023 (high natural gas prices in the western U.S., and high summer temperatures in Texas and the Southwest) were observed to have driven this high value. Of particular note, high prices in Texas occurred at a largely distinct set of hours from high prices in the Southwest, helping to drive up the value of transmission in total and demonstrating significant value to both regions. This example demonstrates the unique value of transmission (compared to other solutions, such as building local generation resources) in delivering benefits to multiple regions given its ability to connect areas of the country that inevitably face differing circumstances.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

AGGREGATE: dAta-driven modelinG preservinG contRollable dEr for outaGe mAnagemenT and rEsiliency (Final Report)

The AGGREGATE project team successfully developed and validated various modules for outage management. Brief summaries of each module are provided to showcase their strength for outage management and restoration for a distribution system with a high penetration of connected distribution energy resources (DERs). In recent years, inverter-based DERs have been widely deployed in distribution system. A most of behind-the-meter (BTM) solar power generation is not visible to the utility. The data-driven DER and load estimation modules are using machine learning (ML) and artificial intelligence (AI) to manage this issue, which provides an opportunity for distribution system operators (DSOs) to operate systems and make decisions in real-time for a distribution system with a high penetration of DERs deployed. Also, the estimated DER and true load can be further leveraged in network aggregation and cold-load pick up estimation for reducing the computing complexity and providing for fast restoration. After load demand and DER power generations have been estimated, the information will support topology and state estimation (SE). The topology estimation module demonstrated the viability of mixed integer linear programming (MILP) formulation to estimate the most likely operational radial topology and outage sections using power flow measurements, historical/estimated load and DERs data and smart meter ping measurements. Formulation includes continuous (power flow, load and DERs data) and binary measurements (smart meter ping measurements) in a single formulation. Errors in continuous data and binary data are modeled as normal distribution and Bernoulli distribution, respectively. In the future distribution grid, the power injection from controllable DERs will be essential for efficient and resilient grid operation. However, determining the optimal DER injections and restoration actions is dependent on knowledge of the system states. State estimation (SE), already the cornerstone of transmission energy management systems, will become commonplace in distribution management systems as more measurements become available from deployment of automated metering infrastructure (AMI). Observability analysis is the first step in SE, as it determines the sufficiency of the available measurements for accurately estimating the current system states. A new type of pseudo-measurement called a Correlational Measurement (CM) is introduced in this module, to enhance the observability of the system to enable more accurate SE. CMs encapsulate knowledge of correlation between demand patterns for similar classes of loads as well as injection patterns for same-technology renewable DERs. During grid contingency scenarios, DERs have been traditionally disconnected, without any fault ride-through capabilities. However, with new regulations and better technology, it is feasible for these resources to contribute to the grid’s restoration after an adverse event and hence enhance resilience. The controllability module proposes a two-step restoration scheme for the power system restoration process by leveraging additional degrees of freedom in power electronics interfaced DERs for mitigating voltage problems. In a resilience mode without the utility system, the distribution grid relies on DERs to serve critical load. In such a severe event with multiple faults on the distribution feeders, actuation of various protective devices (PDs) divides the distribution system into electrical islands. The undetected actuated PDs due to fault current contributions from DERs can delay the restoration process, thereby reducing the system resilience. The Advanced Outage Management (AOM) and the Advanced Feeder Restoration (AFR) modules developed in this project provide improved system resilience with multiple DERs. AOM identifies the faulted sections and actuated PDs in a distribution system with DERs by incorporating smart meter data. The most credible outage scenario including fault locations, PD actuations, and fault indicator (FI) failures is identified by a set of binary integer linear programming incorporating hypotheses. The AFR module serves to restore a distribution system with available energy resources taking into consideration the availability of utility sources and DERs. By partitioning the system into islands, critical load will be served with the available generation resources within islands based on the solution of a MILP. When the utility systems become available, the optimal path will be determined by a spanning tree search algorithm that reconnects these islands back to substations and restores the remaining load. The transmission and distribution (T&D) co-simulation module was used to validate the effect of a control action performed on the distribution side assets as it propagates to the transmission side. This ensures that the control action performed results in a feasible operating point on both the transmission and the distribution system. In addition to validation, the team used the T&D co-simulation module to demonstrate how distribution system assets can be used to mitigate issues on the transmission system. Specifically, the team demonstrated that appropriate switching operations on the distribution side can alleviate the line overload condition on the transmission side without causing new operational constraint violations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Air quality impacts from the development of unconventional oil and gas well pads: Air toxics and other volatile organic compounds

Unconventional oil and natural gas development (UOGD) has expanded rapidly across the United States in recent decades and raised concerns about associated air quality impacts. While significant effort has been made to quantify methane emissions, relatively few observations have been made of Volatile Organic Compounds (VOCs), especially during drilling and completion of new wells. Extensive air monitoring during development of several large, multi-well pads in Broomfield, Colorado, in the Denver-Julesburg Basin, provides a novel opportunity to examine changes in local air toxics and other VOC concentrations during well drilling and completions and production. These operations offer an especially useful case to study as several management practices were implemented to reduce emissions (e.g., electrified, grid-powered drill rigs and closed loop fluid handling systems to reduce truck traffic and limit fluid handling on the pad). With simultaneous measurements of methane and 50 VOCs from October 2018 to December 2022 at as many as 19 sites near well pads, in adjacent neighborhoods, and at a more distant reference location, we identify impacts from each phase of well development and production. Use of weekly, time-integrated canisters, a Proton Transfer Reaction Mass Spectrometer (PTR-MS), continuous photoionization detectors (PID) to trigger canister collection upon detection of VOC-rich plumes, and an instrumented vehicle, provided a powerful suite of measurements to characterize both transient plumes and longer-term changes in air quality. Prior to the start of well development, VOC gradients were small across Broomfield. Once drilling commenced, concentrations of oil and gas (O&G) related VOCs, including alkanes and aromatics, increased around active well pads. Concentration increases were clearly apparent during certain operations, including drilling, coil tubing/millout operations, and production tubing installation. Emissions of C 8 –C 10 n-alkanes during drilling operations highlighted the importance of VOC emissions from synthetic drilling mud chosen to reduce odor impacts. More than 90 samples were collected of transient plumes. Using composition measurements, meteorological data, and information about well pad activities, these plumes were connected with specific UOGD operations including drilling, flowback, and production equipment maintenance. The chemical signatures of these plumes differed by operation type (e.g., C 8 –C 10 n-alkanes constituted a larger fraction of measured VOCs in drilling-related plumes). Concentrations of individual, oil and gas-related VOCs in these plumes were often several orders of magnitude higher than in background air, with maximum ethane and benzene concentrations of 79,600 and 819 ppbv, respectively. Because these plumes typically impact a monitoring site for just several minutes, they are easily missed by slower-responding instruments. Study measurements highlight future emission mitigation opportunities during UOGD operations, including better control of emissions from shakers that separate drill cuttings from drilling mud, production separator maintenance operations, and periodic emptying of sand cans during flowback operations.

54 ENVIRONMENTAL SCIENCES↗

Towards an IPv6-only WLCG: More successes in reducing IPv4

The Worldwide Large Hadron Collider Computing Grid (WLCG) community’s deployment of dual-stack IPv6/IPv4 on its worldwide storage infrastructure has been very successful. Dual-stack is not, however, a viable longterm solution; the HEPiX IPv6 Working Group has focused on studying where and why IPv4 is still being used, and how to flip such traffic to IPv6. The agreed end goal is to turn IPv4 off and run IPv6-only over the wide-area network to simplify both operations and security management.This paper reports our work since the CHEP2023 conference. Firstly, we present our campaign to deploy IPv6 on CPU services and Worker Nodes, with a deadline of end of June 2024. Then, the WLCG Data Challenge (DC24) performed in February 2024 was an excellent opportunity to observe the percentage of data transfers carried by IPv6. We observed the predominance of IPv6 in data transfers during DC24 and were able to understand yet more reasons for the use of IPv4 and areas for remedial action.The paper ends with the working group’s plans for moving WLCG to “IPv6- only”. One aspect of this is the possible automated use of IPv6-only clients configured with a customer-side translator, or CLAT, together with a deployment of NAT64 using what is often known as “IPv6-Mostly”, enabling IPv6-only sites to connect to non-WLCG IPv4-only services.

Attebury, Garhan [U. Nebraska, Lincoln]↗

Data-driven method for electric vehicle charging demand analysis: Case study in Virginia

Electric vehicle (EV) adoption in the U.S. will be accelerated by the historic $7.5 billion public investments in EV charging infrastructure. Careful analysis of EV charging demands plays a vital role in understanding the energy requirements, power grid impact, and smart charging management opportunities of EVs. To this end, this paper develops a data-driven trip-chaining-based modeling framework including five steps: Trip data acquisition and preprocessing, EV adoption modeling, travel itinerary synthesis, EV charging demand simulation and EV load profile generation. The developed analysis framework was demonstrated using real-world data for one region in Virginia, U.S. The results show that the proposed modeling framework can work effectively. For the study region in 2040, the predicted number of plug-in EVs is 470,114, resulting in a weekly charging demand of 38,078,127 kWh (55% home, 9% work, and 36% public) in September and 45,920,358 kWh (61% home, 9% work, and 30% public) in February.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Memory forensic analysis of a programmable logic controller in industrial control systems

In industrial control systems (ICS), programmable logic controllers (PLCs) are used to automate physical processes such as nuclear plants and power grid stations, and are often subject to cyber attacks. As in conventional IT domain, the memory analysis of the PLCs can help answer important forensic questions about the attack, such as the presence of malicious firmware, injection of modified control logic (the program running on the PLC), and manipulation of I/O devices (e.g., sensors and actuators). Unlike conventional IT domain, PLCs have heterogeneous hardware architecture, proprietary firmware and control software, making it challenging to employ a unified framework for their memory forensics. For merely extracting artifacts of forensic importance, reverse-engineering the firmware is a tedious task, and the effort needs to be repeated for every PLC model. As a community, a step-wise approach to tackle this challenge is to analyze the memory of specific PLCs, and subsequently find a generic framework applicable to all PLCs. Our work is a step forward in this direction. By following a methodology that focuses on the functional layer of PLCs instead of reverse engineering the firmware, we analyze the digital forensic artifacts available in a common PLC, Allen-Bradley ControlLogix 1756-L61. Before diving into the memory dump, we analyze the PLC control software to create a list of important artifacts that are sure to exist in the PLC memory dump. The approach employs a setup where PLC control software RSLogix-5000 is connected to the PLC, and the memory dump can be obtained as and when needed. We create test cases that sequentially highlight each category of artifacts, followed by an examination of the resultant impact on memory. After attaining the listed artifacts, we employ conventional string and known data searches to extract interesting information present in this PLC's memory. The memory analysis profile, presented as a Python library and shared with the community, can help a forensic investigator to readily extract forensic artifacts from the same model's controller. The adopted approach may help researchers in creating memory profile of other PLCs, and ultimately formulating a generic PLC memory analysis framework.

Rais, Muhammad Haris↗

Evaluation of Nominal Energy Storage at Existing Hydropower Reservoirs in the US

Long-term planning and operation of hydropower reservoirs require an understanding of both water and energy storage. As energy storage needs of the evolving grid increase, we must account for the water and energy storage potential of these reservoirs. Given the limitations of current data on existing hydropower, we compile statistics related to storage volume and hydraulic head from publicly available data sets and examine differences in descriptions of US hydropower storage. Assembled characteristics are used to calculate nominal energy storage capacity, a simple measure of potential to generate power from a given volume of water, not factoring in detailed constraints. Inventory-based estimates of energy storage are calculated at 2,075 dams, which helps put the potential for US hydropower to support energy storage in context with similar evaluations in other regions and with other energy storage technologies. The national energy storage capacity ranges between 34.5 and 45.1 TWh depending on the information used, with 52% of energy storage located at the 10 largest reservoirs in the US. Energy storage capacities are also calculated at 236 dams with historical volume and elevation data. Finally, reservoir inflows provide context for the storage volumes and sensitivities to hydrologic variability. Larger reservoirs with greater storage volume to inflow ratios are concentrated in the Western US, but the majority of hydropower reservoirs store less than the annual inflow. We address several infrastructure and water resource informatics challenges and highlight remaining issues, including representing seasonal or shorter variability in water volumes and representing connected hydropower facilities.

13 HYDRO ENERGY↗

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Despite decades of practice, finite-size errors in many widely used electronic structure theories for periodic systems remain poorly understood. For periodic systems using a general Monkhorst-Pack grid, there has been no comprehensive and rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order Møller-Plesset perturbation theory (MP2), which are the simplest wavefunction based method, and the simplest post-Hartree-Fock method, respectively. Such calculations can be viewed as a multi-dimensional integral discretized with certain trapezoidal rules. Due to the Coulomb singularity, the integrand has many points of discontinuity in general, and standard error analysis based on the Euler-Maclaurin formula gives overly pessimistic results. The lack of analytic understanding of finite-size errors also impedes the development of effective finite-size correction schemes. We propose a unified analysis to obtain sharp convergence rates of finite-size errors for the periodic HF and MP2 theories. Our main technical advancement is a generalization of the result of Lyness [Math. Comp. 30 (1976), pp. 1–23] for obtaining sharp convergence rates of the trapezoidal rule for a class of non-smooth integrands. Our result is applicable to three-dimensional bulk systems as well as low dimensional systems (such as nanowires and 2D materials). Our unified analysis also allows us to prove the effectiveness of the Madelung-constant correction to the Fock exchange energy, and the effectiveness of a recently proposed staggered mesh method for periodic MP2 calculations (see X. Xing, X. Li, and L. Lin [J. Chem. Theory Comput. 17 (2021), pp. 4733–4745]). In conclusion, our analysis connects the effectiveness of the staggered mesh method with integrands with removable singularities, and suggests a new staggered mesh method for reducing finite-size errors of periodic HF calculations.

97 MATHEMATICS AND COMPUTING↗

A Randomization-Based, Zero-Trust Cyberattack Detection Method for Hierarchical Systems

This paper demonstrates a novel randomization-based approach for verifying power system control signals with application to detecting cyberattacks. We consider fully connected hierarchical systems containing multiple local agents and a global "trust" agent. The global agent uses a time-varying randomized assignment scheme to identify corrupt network links based on principles of zero trust and majority rule. To evaluate the performance of this detection approach, we implement our algorithm in MATLAB and run it against nearly 43 million unique attack scenarios spanning a range of system sizes. For each scenario, the algorithm determines whether the identified corruptions satisfy a set of validity constraints reflecting network topology and uses that result to say whether the recovered state value for one or more local agents is malicious. We compare the algorithm's determination to the true state of the system to assess performance and find that classification accuracy converges to 100% as system size increases, suggesting that the validity constraints become more difficult to satisfy for larger systems. We further explore the scenarios that evade detection to understand practical implications for employing this detection approach.

cybersecurity↗