Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid Independent”

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 91 records · Page 5

Renewable-battery hybrid power plants in congested electricity markets: Implications for plant configuration

Examining coupled renewable-battery power plants (“hybrids”) in congested areas provides insights into a future of increased wind and solar penetration. Our study focuses on two types of congested regions, Variable Renewable Energy (VRE)-rich Areas and Load Centers, and explores likely plant configuration choices for developers and transmission network planners. Here, this paper examines how hybrid value, comprising energy and capacity value, varies by plant configuration and congested region type considering factors such as storage duration, battery degradation, and ability to charge from the grid. We select plant locations from across the seven main U.S. independent system operators (ISOs). Hybrid value for each configuration is computed based on profit-maximizing plant operation given perfect foresight, according to observed wholesale power market real time prices from 2018 to 2021. In VRE-rich Areas, the median increase in energy value from extending storage duration from one to 4h is 29.4% for solar and 26.8% for wind, assuming low battery degradation costs and storage sized to 100% of the plant's nameplate generation capacity. Increasing storage duration beyond 4h does not substantially increase its value from energy markets, even in VRE-rich Areas. We find that solar hybrids reach a 90% capacity credit with 4h of storage, while wind hybrids require 8h of storage, based on the capacity factor of each hybrid during the top 100 net load hours.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Object-oriented analysis as a foundation for building climate storylines of compounding short-term drought and crop heat stress

Introduction: Crops are vulnerable to precipitation and heat extremes during late spring through summer. Methods: We analyzed for a north-central U.S. region short-term drought and agricultural heat stress during April-May-June-July. We used the 4-km Parameter Elevation Regression on Independent Slopes Model (PRISM) for observations, aggregated to a 25-km grid, and two 25-km Regional Climate Model version 4 (RegCM4) simulns used either GFDL- or MPI-GCM boundary conditions. We chose 1981-2000 as our contemporary time period, and 2041- 2060 as our scenario time period, which used the Representative Concentration Pathway 8.5 emissions scenario. We used object-oriented analysis to identify events of interest in observations and simulations by identifying objects in a space-time domain that meet specified criteria, such as exceeding a heat-stress temperature threshold. The event diagnosis allowed analysis of compound events, occurring when temperature and drought objects overlap. Results: Identified objects yielded events that can undermine agricultural productivity and which are thus relevant to decision makers, making them building blocks for possible climate storylines. The observations and simulations showed similar spatial distributions of event frequencies across the analysis region. However, the simulations attained this distribution by having fewer events that tend to cover larger areas compared to observed events, suggesting that the effective resolution of the simulations was coarser than their 25-km grids. Short-term drought frequency increased and heat-stress frequency decreased in transitioning to the scenario climate. When compounding occurred heat-stress events generally preceded the short-term drought events. The overlapping, compound events tended to be more extreme compared to non-overlapping events of either type. Discussion: The information yielded projected changes in these agriculturally motivated events. One prominent conditional behavior emerging from the work was that a heat-stress event should be a warning to watch for potential drought, as both could compound each other to more intense levels.

54 ENVIRONMENTAL SCIENCES↗

Multi-cell, triple pane, vacuum insulated glazing

A novel Triple Pane Vacuum Insulated Glazing (TPVIG) architecture consisting of three panes of glass separated by two offset grids is described. The grids are made of frit-coated fibers fused with the glass panes in a vacuum to create many independent hermetic cells. Unlike traditional VIGs, which are separated by support pillars, these cells significantly increase the redundancy in sealing, and the additional bonds formed by the grid lines reduce the stresses on the edge bonds thereby increasing VIG reliability. A numerical thermal analysis was conducted to determine the thermal performance of the VIG and suggested U-values below 0.5 W/m 2 -K (greater than R=11.4 hr-ft 2 -°F/Btu) can be obtained. A prototype of the VIG unit consisting of three glass panes (5.7 mm, 1.1 mm and 5.7 mm thick with low-e coating on two surfaces) was fabricated and tested for thermal performance using the hot box method. The experimental center of glass U-value for the VIG was found to be 0.42 W/m 2 K (R=13.2 hr-ft 2 -°F/Btu). Furthermore, this U-value is about four times lower than a typical Argon filled double pane window.

14 SOLAR ENERGY↗

Long-Duration Energy Storage Grid Integration-Valuation Framework and Incentive Gaps

Given these challenges and current modeling gaps on Long Duration Energy Storage (LDES), enhancing the structure and design of existing planning, operations, and organized wholesale markets can better characterize the value of LDES to the power system. To more thoroughly assess the gaps and barriers to LDES investment and readiness for integration into a future grid, we conducted stakeholder outreach through an online survey, interviews with individual independent system operators/regional transmission organizations, and a literature review. Based on this assessment, we identified a set of opportunities for LDES focused development, including a framework to quantify the contributions of LDES on resource adequacy, reliability, and resiliency. Specifically, we identify the potential demand for and benefits of an open-source, LDES-centric evaluation framework that can guide future planning, operations, market design, and policy reforms.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimal local truncation error method for 3-D elasticity interface problems

The paper deals with a new effective numerical technique on unfitted Cartesian meshes for simulations of heterogeneous elastic materials. Here, we develop the optimal local truncation error method (OLTEM) with 27- point stencils (similar to those for linear finite elements) for the 3-D time-independent elasticity equations with irregular interfaces. Only displacement unknowns at each internal Cartesian grid point are used. The interface conditions are added to the expression for the local truncation error and do not change the width of the stencils. The unknown stencil coefficients are calculated by the minimization of the local truncation error of the stencil equations and yield the optimal second order of accuracy for OLTEM with the 27-point stencils on unfitted Cartesian meshes. A new post-processing procedure for accurate stress calculations has been developed. Similar to basic computations it uses OLTEM with the 27-point stencils and the elasticity equations. The post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing numerical techniques (e.g., with finite elements). Numerical experiments show that at an accuracy of 0.1% for stresses, OLTEM with the new post-processing procedure significantly (by 10 5 -10 9 times) reduces the number of degrees of freedom compared to linear finite elements. OLTEM with the 27-point stencils yields even more accurate results than high-order finite elements with wider stencils.

42 ENGINEERING↗

Generator Interconnection Cost Analysis in the Midcontinent Independent System Operator (MISO) territory

Electric transmission system operators (ISOs, RTOs, or utilities) require new large generators seeking to connect to the grid to undergo a series of impact studies before they can be built. This process establishes what new transmission equipment or upgrades may be needed before a project can connect to the system and assigns the costs of that equipment. Berkeley Lab has collected interconnection cost data from interconnection studies for the Midcontinent Independent System Operator (MISO), representing nearly 50% of all projects requesting interconnection from 2010 to 2020. Project-level cost summary data are available for download on this page. We find: -Average interconnection costs have grown as the number of interconnection requests have escalated -Projects that have completed all required interconnection studies have the lowest cost compared to applicants still actively working through the interconnection process or those that have withdrawn. -Broader network upgrade costs are the primary driver of recent cost increase. -Potential interconnection costs for wind, storage, and solar are larger than for natural gas -Larger generators have greater interconnection costs in absolute terms, but economies of scale exist on a per kW basis. -Interconnection costs vary by location Berkeley Lab will publish a series of short analytical papers of generator interconnection costs to the transmission system for MISO, PJM, SPP, ISO-NE and NYISO, which you can find at https://emp.lbl.gov/interconnection_costs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Optimizing power system restoration with damaged communications

Utility procedures for power system blackstart and restoration typically assume that energization decisions can be reliably communicated across the grid. In reality, the communications and control network would likely also be affected in power outages, such as those caused by extreme weather events or cyber-attacks. This paper studies the effect of damage to the power system communications and control infrastructure on restoration operations following a blackout. We model the communications infrastructure as a graph, overlaying the power grid, and imposing the requirement that every energized element in the power grid be observable from a control center. We expand on a specialized branch-and-bound algorithm from the literature to optimize the restoration process and devise an initialization heuristic and a rounding heuristic to improve solution speed. We perform numerical experiments on synthetic systems for Illinois and Texas with outages based on a solar flare or hurricane. We compare the results of our specialized branch-and-bound algorithm to the results from (i) the initialization heuristic alone, (ii) a variation of this heuristic that we use as a baseline, and (iii) the restoration optimization for the power system without communications constraints. Here, we find that damage to the communications infrastructure significantly increases the time required to re-energize the grid. Moreover, by simultaneously optimizing communications repairs and grid energization decisions, we are able to re-energize the grid significantly faster than if communications repairs and energization decisions were made independently or with partial coordination, motivating improvements to current industry practice.

97 MATHEMATICS AND COMPUTING↗

An embedded 3D fracture modeling approach for simulating fracture-dominated fluid flow and heat transfer in geothermal reservoirs

In this paper, we describe an efficient modeling approach, named embedded discrete fracture method (EDFM), for incorporating arbitrary 3D, discrete fractures, such as hydraulic fractures or faults, into modeling fracture-dominated fluid flow and heat transfer in fractured geothermal reservoirs. This technique allows 3D discrete fractures to be discretized independently from surrounding rock volume and inserted explicitly into a primary fracture/matrix grid, generated without including 3D discrete fractures in prior. An effective computational algorithm is developed to discretize these 3D discrete fractures and construct local connections between 3D fractures and fracture/matrix grid blocks representing the surrounding rock volume. The constructed gridding information on 3D fractures is then added to the primary grid. This embedded fracture modeling approach can be directly implemented into a developed geothermal reservoir simulator via the integral finite difference (IFD) method or with TOUGH2 technology. This embedded fracture modeling approach is very promising and computationally efficient to handle realistic 3D discrete fractures with complicated geometries, connections, and spatial distributions. Compared with other fracture modeling approaches, it avoids cumbersome 3D unstructured, local refining procedures, and increases computational efficiency by simplifying Jacobian matrix size and sparsity, while maintaining enough accuracy. Several numeral simulations are presented to demonstrate the utility and robustness of the proposed technique. Our numerical experiments show that this approach captures all the key patterns about fluid flow and heat transfer dominated by fractures in these cases. Thus, this approach is readily available to the simulation of fractured geothermal reservoirs with both artificial and natural fractures.

15 GEOTHERMAL ENERGY↗

A Tutorial on the Control of Floating Offshore Wind Turbines: Stability Challenges and Opportunities for Power Capture

Climate change is a serious threat facing humanity. The United States (U.S.) and many other countries are increasing the amount of electrical power generated from renewable energy sources in an effort to combat climate change and ensure energy independence. The U.S. has set goals to achieve a 100% decarbonized electric grid by 2035 and a net-zero emissions economy by no later than 2050. Renewable energy currently accounts for about 20% of the U.S. power grid. According to the U.S. Department of Energy and the National Renewable Energy Laboratory (NREL), in the U.S. in 2021, wind and solar photovoltaic generation supplied 9.1% and 4% of total electricity generation, respectively, and the latest publicly available data show that hydropower represented 6.6% of all electricity generated in the U.S. in 2019. As wind farms have been built in many of the best wind resource areas on land, the U.S. and many other countries are turning to offshore wind for further growth of wind power capacity. Further, the U.S. has committed to deploy 30 GW of offshore wind by 2030, a significant increase from the 0.04 GW of installed offshore wind in the country as of the end of 2023. Many other countries have also established ambitious goals or plans to increase the amount of installed offshore wind power.

17 WIND ENERGY↗

Deploying E3’s RESERVE Tool to Enable Advanced Operation of Clean Grids

Energy and Environmental Economics, Inc. (E3) developed an open-source machine learning model, RESERVE, for deriving ancillary services timeseries in deeply decarbonized electricity grids. E3 used a bespoke PLEXOS production simulation model of the California Independent System Operator’s (CAISO) balancing area to validate RESERVE’s ability to enable production cost, greenhouse gas emissions (GHG), and renewable energy curtailment savings. These savings were modeled by comparing PLEXOS cases with RESERVE’s outputs to PLEXOS cases with CAISO’s incumbent reserve product in the Western Energy Imbalance Market (EIM)’s 15-minute market. E3 also tested cases with solar operating flexibly to provide reserves. E3 found that, in a 2030 modeling year, using RESERVE and flexible solar enabled significant production cost, GHG and curtailment savings versus the incumbent CAISO method in cases with low penetrations of lithium-ion batteries. However, with the full 14 gigawatts (about 30% of peak CAISO demand) of 4-hour lithium-ion batteries that are expected to be installed by 2030, these savings approach zero due to batteries saturating ancillary services markets. E3 also found significant savings under a 2019 benchmarking year.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Electron bunch dynamics and emission in particle-in-cell simulations of relativistic laser–solid interactions: On density artifacts, collisions, and numerical dispersion

Sub-optical-cycle dynamics of dense electron bunches in relativistic-intensity laser–solid interactions lead to the emission of high-order harmonics and attosecond light pulses. The capacity of particle-in-cell simulations to accurately model these dynamics is essential for the prediction of emission properties because the attosecond pulse intensity depends on the electron density distribution at the time of emission and on the temporal distribution of individual electron Lorentz-factors in an emitting electron bunch. Here, we show that in one-dimensional collisionless simulations, the peak density of the emitting electron bunch increases with the increase in the spatial resolution of the simulation grid. When collisions are added to the model, the peak electron density becomes independent of the spatial resolution. Collisions are shown to increase the spread of the peaks of Lorentz-factors of emitting electrons in time, especially in the regimes far from optimum generation conditions, thus leading to lower intensities of attosecond pulses as compared to those obtained in collisionless simulations.

Physics↗

Review of the November 11, 1905 Earthquake: Seismic Intensity Data, Magnitude, and Location Report

The macroseismic location and magnitude of the Shoshone 1905 earthquake has been based on the application of objective and quantitative approaches that produce estimates of earthquake source location and magnitude, and associated uncertainties. The performance of all three approaches applied depends critically on the characteristics of the seismic intensity field. Regional intensity attenuation relations were derived based on a calibration set of earthquakes and then, the procedure was validated using three independent very recent events. The macroseismic epicenter location was determined by either a grid search procedure, the spatial distribution of the sites experiencing the largest intensities, or a combination of both. The macroseismic magnitude estimates are based on the derived regional intensity attenuation relations. There is consistency in the location and magnitude of the November 1905 earthquake as obtained by the three methods. The overlapping region of the 68% bootstrap resampling solutions suggests a reliable epicenter location within 41.550N ? 42.060N and 113.478W ? 112.320W. This region encloses the original Oaks (1992) epicenter location. Macroseismic magnitude estimates range between 5.87 and 6.05.

58 GEOSCIENCES↗

Bias Corrected NOAA HRRR Wind Resource Data for Grid Integration Applications

To address the need for regularly updated wind resource data, NREL has processed the High-Resolution Rapid Refresh (HRRR) outputs for use in grid integration modeling. The HRRR is an hourly-updated operational forecast product produced by the National Oceanic and Atmospheric Administration (NOAA) (Dowell et al., 2022). Several barriers have prevented the HRRR's widespread proliferation in the wind energy industry: missing timesteps (prior to 2019), challenging file format for wind energy analysis, limited vertical height resolution, and negative bias versus legacy WIND Toolkit data (2007-2013). NREL has applied re-gridding, interpolation, and bias-correction to the native HRRR data to overcome these limitations. This results in the now-publicly-available bias corrected and interpolated HRRR (BC-HRRR) dataset for weather years 2015 to 2023. Bias correction is necessary for wind resource consistency across weather years to be used simultaneously in planning-focused grid integration studies alongside the original WIND Toolkit data. We show that quantile mapping with the WIND Toolkit as a historical baseline is an effective method for bias correcting the interpolated HRRR data: the BC-HRRR has reduced mean bias versus comparable gridded wind resource datasets (+0.12 m/s versus Vortex) and has very low mean bias versus ground measurement stations (+0.01 m/s) (Buster et al., 2024). BC-HRRR's consistency with the legacy WIND Toolkit allows NREL to extend grid integration analysis to 15+ weather years of wind data with low-overhead extensibility to future years as they are made available by NOAA. As with historical datasets like the WIND Toolkit, BC-HRRR is intended for use in grid integration modeling (e.g., capacity expansion, production cost, and resource adequacy modeling) both independently and alongside the legacy WIND Toolkit.

Array↗

Scaling and performance portability of the particle-in-cell scheme for plasma physics applications through mini-apps targeting exascale architectures

We perform a scaling and performance portability study of the particle-in-cell scheme for plasma physics applications through a set of mini-apps we name "Alpine", which can make use of exascale computing capabilities. The mini-apps are based on Independent Parallel Particle Layer, a framework that is designed around performance portable and dimension independent particles and fields. We benchmark the simulations with varying parameters such as grid resolutions (5123 to 20483) and number of simulation particles (109 to 1011) with the following mini-apps: weak and strong Landau damping, bump-on-tail and two-stream instabilities, and the dynamics of an electron bunch in a charge-neutral Penning trap. We show strong and weak scaling and analyze the performance of different components on several pre-exascale architectures such as Piz-Daint, Cori, Summit and Perlmutter. While the scaling and portability study helps identify the performance critical components of the particle-in-cell scheme in the current state-of-the-art computing architectures, the mini-apps by themselves can be used to develop new algorithms and optimize their high performance implementations targeting exascale architectures.

Muralikrishnan, Sriramkrishnan↗

The impact of forest-controlled snow variability on late-season streamflow varies by climatic region and forest structure

Previous studies have documented how forests influence snow at fine spatial scales, but none have documented the influence that existing forest-snow variability has on streamflow. Here to test how much forest-controlled snow variability influences streamflow, a tiling parameterization based on classifications from high-resolution (1–3 m) vegetation maps was incorporated into the Distributed Hydrology and Soil Vegetation Model (DHSVM). Within each grid cell (90–150 m), the tiling parameterization simulated forest-snow variability with four independently evolving snowpacks. Each tile had unique radiation conditions to represent conditions underneath the canopy, in exposed areas, and along north- and south-facing forest edges. This tiled parameterization was used to test where and when detailed forest-snow modelling should be considered further and where and when the impacts are too small to be worth the effort. To test this, tiled model simulations of streamflow were compared to non-tiled model simulations in the Sierra Nevada, CA, the Jemez Mountains, NM, and the Eastern Cascades, WA. In Tuolumne, CA, the tiled model simulated little difference in grid cell average SWE, and late-season streamflow decreased by only 3%–4% compared to the non-tiled model. In Jemez, NM, the tiled model decreased late-season streamflow by 18% due to increased sublimation. In Chiwawa, WA, the tiled model increased late-season streamflow by 15% due to high shortwave radiation attenuation and less longwave radiation enhancement from the forest. Furthermore, within the Chiwawa, a substantial silvicultural practice was synthetically implemented to increase the north-facing edge's fractional area. This silvicultural experiment, which used the same fractional forest area in all simulations increased late-season streamflow by 35% compared to tiled model simulations that did not represent forest edges. In conclusion, representing forest-SWE variability had an effect on late-season streamflow in some watersheds but not in others based on the fractional area of the forest edges, forest characteristics, and climate conditions.

54 ENVIRONMENTAL SCIENCES↗

Co-simulation Framework for Community-scale Building-grid Integration [SWR-21-75]

Distributed energy resources (DERs), including rooftop solar, energy storage, and flexible loads, are gaining popularity as costs decline and as building owners and utilities realize their benefits. DERs can improve distribution system efficiency, help prevent the need for expensive grid upgrades, and increase the resilience of local communities. However, they can also cause difficulties in grid operations and can require controls to achieve their benefits. To address this challenge, NREL researchers have developed a community-scale solution that assesses the impacts of DERs and their control strategies on a distribution system. The framework has been shown to reduce solar photovoltaic (PV) curtailment to 0%, mitigate the adverse impact of solar variability on the distribution voltage, and provide up to 5-day critical load support during emergency events. Utilizing 5 different modules representing the feeder, buildings, home energy management systems, an aggregator, and a utility controller, NREL expects this simulation technology to play a critical role in the continued integration of DERs. According to the Energy Information Administration (EIA), solar curtailments accounted for 94% of the total energy curtailed in the California Independent System Operator (CAISO) in 2020. By enabling Independent System Operators (ISOs) and utility operators to bring solar curtailments to 0%, the electrical grid can become less dependent on fossil-fueled power generation sources. NREL's co-simulation framework contains five major components: Distribution Feeder Model: describes the distribution feeder topology using OpenDSS, including the locations of all DERs. Residential Building Model: simulates a large number of buildings at a high resolution using OCHRETM. The model is equipped to control equipment based on signals from an external module. The model includes major household appliances such as HVAC and a water heater, non-dispatchable load models, a distributed PV system, and a home battery system. Home Energy Management System: optimizes the controls for the devices in a home using foreseeTM. The control can adjust based on the user preferences including cost, comfort, and convenience. In hierarchical control scenarios, where the houses follow signals from an aggregator, the home energy management system provides a flexibility band with a range of power and follows the dispatch signals received from aggregator. Community-Level Aggregator: solves for optimal energy dispatch based on the flexibility bands received from each home and the grid service signal received from the utility controller. Utility-Level Controller: provides grid signals for voltage control using Distributed Energy Resources (DERs), such as solar systems, in the community.

Balamurugan, Sivasathya Pradha↗

Alternative mixed integer linear programming optimization for joint job scheduling and data allocation in grid computing

This paper presents a novel approach to the joint optimization of job scheduling and data allocation in grid computing environments. We formulate this joint optimization problem as a mixed integer quadratically constrained program. To tackle the nonlinearity in the constraint, we alternatively fix a subset of decision variables and optimize the remaining ones via Mixed Integer Linear Programming (MILP). We solve the MILP problem at each iteration via an off-the-shelf MILP solver. Our experimental results show that our method significantly outperforms existing heuristic methods, employing either independent optimization or joint optimization strategies. We have also verified the generalization ability of our method over grid environments with various sizes and its high robustness to the algorithm setting.

97 MATHEMATICS AND COMPUTING↗

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗