Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid Optimization”

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 1,747 records · Page 97

eGridGPT: Trustworthy AI in the Control Room

This report outlines the synergy between human decision making and generative artificial intelligence (GenAI), where GenAI supports power system operators by analyzing procedures, suggesting actions, simulating scenarios with physics-based digital twins, and recommending optimal decisions. This report is the first research effort to apply large language models (LLMs), a type of GenAI, in the power grid control room. The authors describe the Electric Grid Generative Pretrained Transformer (eGridGPT), an LLM that virtually assists system operators. Developed with cybersecurity and regulatory requirements in mind, eGridGPT represents an opportunity to responsibly evolve control room technologies to meet the needs of a rapidly changing grid. As an innovative concept, eGridGPT seeks to spearhead productive discussions about the advanced technologies in the control room of the future amid the transition to clean energy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An Integrated Paradigm for the Management of Delivery Risk in Electricity Markets: From Batteries to Insurance and Beyond

If power systems transition to integrate higher amounts of variable renewable energy sources, storage technologies, and distributed energy resources (DERs), new risk management frameworks are necessary to ensure cost-effective and reliable power system operations. Projects funded by the Advanced Research Projects Agency-Energy (ARPA-E) Performance-based Energy Resource Feedback, Optimization, and Risk Management (PERFORM) program aim to contribute new risk management frameworks by developing methods to quantify and manage risk at grid asset and system levels. The National Renewable Energy Laboratory (NREL) led a PERFORM project in collaboration with the Johns Hopkins University, the Electric Power Research Institute (EPRI), kWh Analytics, Packetized Energy, and Imperial Consultants (ICON). The project addressed two challenges related to risk management in electricity markets: managing net load imbalances and flexibility from DERs. This final technical report presents a list of project accomplishments, activities, and outputs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Fourier analysis of the SOR iteration

The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor omega. It is shown that for the standard model problem of Poisson's equation on a rectangle, the optimal omega and corresponding convergence rate can be rigorously obtained by Fourier analysis. The trick is to tilt the space-time grid so that the SOR stencil becomes symmetrical. The tilted grid also gives insight into the relation between convergence rates of several variants.

Leveque, R. J.↗

Telerobotics - Display, control, and communication problems

An experimental telerobotics simulation is described suitable for studying human operator (HO) performance. Simple manipulator pick-and-place and tracking tasks allowed quantitative comparison of a number of calligraphic display viewing conditions. An enhanced perspective display was effective with a reference line from target to base, with or without a complex three-dimensional grid framing the view. This was true especially if geometrical display parameters such as azimuth and elevation were arranged to be near optimal. Quantitative comparisons were made possible, utilizing control performance measures such as root mean square error. There was a distinct preference for controlling the manipulator in end-effector Cartesian space for the primitive pick-and-place task, rather than controlling joint angles and then, via direct kinematis, the end-effector position. An introduced communication delay was found to produce decrease in performance. In considerable part, this difficulty could be compensated for by preview control information. The fact that neurological control of normal human movement contains a sampled data period of 0.2 s may relate to this robustness of HO control to delay.

Stark, Lawrence↗

Three-dimensional Euler time accurate simulations of fan rotor-stator interactions

A numerical method useful to describe unsteady 3-D flow fields within turbomachinery stages is presented. The method solves the compressible, time dependent, Euler conservation equations with a finite volume, flux splitting, total variation diminishing, approximately factored, implicit scheme. Multiblock composite gridding is used to partition the flow field into a specified arrangement of blocks with static and dynamic interfaces. The code is optimized to take full advantage of the processing power and speed of the Cray Y/MP supercomputer. The method is applied to the computation of the flow field within a single stage, axial flow fan, thus reproducing the unsteady 3-D rotor-stator interaction.

Boretti, A. A.↗

Multigrid with Overlapping Patches

Solving boundary value problems with optimal efficiency requires adaptivity and multilevel techniques. Previously, an implementation of the AFACx algorithm is presented that is based on rectangular Cartesian grids. This implementation does not allow for the over]ap of grids that lie on the same level of refinement. We investigate the case in which these grids overlap. A standard technique for overlapping grids is the Schwarz algorithm. Some ways of using the Schwarz algorithm in a standard multigrid scheme are presented. Also, a problem that arises in some situations with non-aligned, overlapping grids is described. This situation comes up in a natural way when the Schwarz algorithm is used as a relaxation scheme within a multilevel algorithm. We identify the reason for the bad convergence and show that by more sophisticated interpolation the difficulties can be overcome. Then we present a multiplicative Schwarz algorithm for a large number of grids that has a high potential for parallelization. Finally we give some numerical results for the FACx algorithm with overlapping grids on each refinement level. The implementation of the described codes uses C++ and the array class libraries A++ and P++. Using the A++/P++ programming environment, it was possible to move from a serial code to a parallel code within a few days.

Berndt, Markus↗

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗

Gas Generator Feedline Orifice Sizing Methodology: Effects of Unsteadiness and Non-Axisymmetric Flow

Engine LH2 and LO2 gas generator feed assemblies were modeled with computational fluid dynamics (CFD) methods at 100% rated power level, using on-center square- and round-edge orifices. The purpose of the orifices is to regulate the flow of fuel and oxidizer to the gas generator, enabling optimal power supply to the turbine and pump assemblies. The unsteady Reynolds-Averaged Navier-Stokes equations were solved on unstructured grids at second-order spatial and temporal accuracy. The LO2 model was validated against published experimental data and semi-empirical relationships for thin-plate orifices over a range of Reynolds numbers. Predictions for the LO2 square- and round-edge orifices precisely match experiment and semi-empirical formulas, despite complex feedline geometry whereby a portion of the flow from the engine main feedlines travels at a right-angle through a smaller-diameter pipe containing the orifice. Predictions for LH2 square- and round-edge orifice designs match experiment and semi-empirical formulas to varying degrees depending on the semi-empirical formula being evaluated. LO2 mass flow rate through the square-edge orifice is predicted to be 25 percent less than the flow rate budgeted in the original engine balance, which was subsequently modified. LH2 mass flow rate through the square-edge orifice is predicted to be 5 percent greater than the flow rate budgeted in the engine balance. Since CFD predictions for LO2 and LH2 square-edge orifice pressure loss coefficients, K, both agree with published data, the equation for K has been used to define a procedure for orifice sizing.

Rothermel, Jeffry↗

Power system load flexibility forecasting

The example embodiments are directed to a system and method for forecasting load flexibility of a power grid. In one example, the method includes receiving temperature values associated with temperature set points of a plurality of loads that are included on a power grid, forecasting a flexibility of the plurality of loads using a polynomial-time mixed-integer non-linear programming (MINLP) optimization based on the received temperature values for the plurality of loads, and outputting information about the forecasted flexibility for display to a display device. The MINLP optimization performs the forecasting of the load flexibility on a fine-grained basis in comparison to conventional methods and is still fast enough that it can be computed in real-time.

Genc, Sahika↗

Summary: Techno-Economic Analysis of Solar Photovoltaics and Battery Energy Storage at a Vietnam Industrial Park

The Clean Energy Investment Accelerator conducted a case study analysis of battery energy storage system (BESS) feasibility for an industrial park in Vietnam using the National Renewable Energy Laboratory's (NREL's) REopt platform (a distributed energy modeling and optimization tool) to evaluate how BESS may reduce electricity costs, increase utilization of onsite renewable energy generation, and improve resilience to utility grid outages. This presentation summarizes the analysis and key takeaways.

batteries↗

Development of an advanced multi-scale subchannel model for CRUD risk assessments

AREVA NP has developed a process, based upon standard licensing methods and specialized applications, for the assessment of the risk of crud induced corrosion and crud induced axial off-set. This streamlined process allows for the evaluation of the crud risk for full core subchannel by subchannel models, by applying local conditions on each fuel rod/subchannel, with axial nodes spanning several centimeters. AREVA NP has applied these tools successfully to units from various NSSS vendors. To date, there have been no plant performance issues that have been outside the performance as predicted by the AREVA NP risk assessment tools. The study presented in this paper utilized the core simulator and thermal-hydraulic packages of this tool set to evaluate the effects of changing various fuel cycle design (e.g. radial zoning patterns) and fuel assembly design parameters (e.g. fuel rod diameter, spacer grid placement) on thermal-hydraulic key decision variables (KDVs) of interest. The study may be valuable as a first step towards optimizing the fuel cycle design, optimizing the fuel lattice design, or evaluating potential operating changes such as a power up-rate or modification of the cycle length for a given PWR plant. The various geometric inputs to the AREVA NP risk assessment process were tested by individually adjusting a single parameter and observing the effect. Multi-parameter studies are also included to assess the impact of optimized fuel cycle designs. The sensitivity of the model, in terms of thermal-hydraulic results, is presented. Provided a relationship between the thermal-hydraulic results to both crud deposition thickness and lithium uptake in the zirconium cladding oxide layer, approximate impacts to the risk of crud induced corrosion and crud induced axial off-set can be estimated. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Techno-Economic Assessment of Data Center Load Demand Powered by Small Modular Reactors and Distributed Energy Resources

The rapid increase in data center energy demand, driven by AI and large-scale data processing, poses significant challenges to global energy infrastructure. Data centers require substantial and reliable energy for continuous operations and high-performance computing. Current electrical grids face issues such as transmission bottlenecks and aging infrastructure, making it difficult to meet these demands. Integrating inverter-based-resources (IBRs) like solar and wind presents both opportunities and challenges due to their intermittent nature. Small Modular Reactors (SMRs) offer a promising solution with their enhanced safety, modularity, reliability, and scalability, providing consistent base load power ideal for data center operations. This study presents a comprehensive techno-economic assessment of powering data center load demand using a combination of SMRs and IBRs with grid-connected and islanded mode. This study utilized Idaho National Laboratory’s (INL) HPC data center hourly load profiles and Xendee microgrid optimization platform to conduct the analysis. In this configuration, SMRs serves as the primary base load power source, consistently providing a steady supply of electricity necessary to meet the minimum load demand of the data center with support from the IBRs. Key performance indicators such as Levelized Cost of Electricity (LCOE), Net Present Value (NPV) has been calculated to assess the economic feasibility. The findings from this research will underscore the strategic benefits of integrating SMR plant with DERs – particularly for critical infrastructure load such as data centers.

14 - SOLAR ENERGY↗

Impact of Blockchain Delay on Grid-Tied Solar Inverter Performance

This paper investigates the impact of the delay resulting from a blockchain, a promising security measure, for a hierarchical control system of inverters connected to the grid. The blockchain communication network is designed at the secondary control layer for resilience against cyberattacks. To represent the latency in the communication channel, a model is developed based on the complexity of the blockchain framework. Taking this model into account, this work evaluates the plant’s performance subject to communication delays, introduced by the blockchain, among the hierarchical control agents. In addition, this article considers an optimal model-based control strategy that performs the system’s internal control loop. The work shows that the blockchain’s delay size influences the convergence of the power supplied by the inverter to the reference at the point of common coupling. In the results section, real-time simulations on OPAL-RT are performed to test the resilience of two parallel inverters with increasing blockchain complexity.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Optimal Energy Storage System and Smart Switch Placement in Dynamic Microgrids With Applications to Marine Energy Integration

Here this paper studies a dynamic microgrid (DMG) planning problem that places energy storage systems (ESSs) and smart switches (SSWs) optimally in the system. We apply the proposed methodology to applications concerning marine renewable energy (MRE). MRE is an emerging clean energy resource with enormous capacity but volatile and intermittent energy output profiles. Innovative grid-integration technologies designed to enhance the reliability of an MRE-integrated system are needed. However, there are still limited studies in this regard. Existing works have shown the promising prospect of using a dynamic microgrid (DMG) operational concept to accommodate renewable resources in distribution systems, but they usually assume fixed ESS and SSW installations. To further improve the operational flexibility of DMGs, we propose a DMG planning methodology that optimally places ESSs and SSWs so that a DMG with MRE is warranted with proper resource adequacy and topological flexibility in both the contingency and normal operations. We use realistic case studies based on a real-world distribution network and the U.S. Department of Energy's MRE dataset to verify the value and validity of the proposed work.

25 ENERGY STORAGE↗

Optimizing the physical design and layout of a resilient wind, solar, and storage hybrid power plant

We report as wind and solar technologies improve and their costs decrease, the share of power produced by these sources will increase. As the market penetration increases, these power sources will need to provide grid services, such as dispatchability, in addition to providing energy. One way to reduce variability, provide higher quality power to the grid, and address local grid stability issues is through colocating wind and solar power plants. In addition to operating reliably during normal operating conditions, in scenarios with high penetrations of renewable generation, it is important that these hybrid plants can withstand production disruptions and continue to supply power despite prolonged resource reduction, extreme weather events, or other disruptions. In this paper, we present a methodology to optimize a wind-solar-battery hybrid power plant down to the component level that is resilient against production disruptions and that can continually produce some minimum required power. We introduce the models and assumptions we used to simulate a hybrid power plant as well as the design variable parameterization and specific methods we used to optimize the plant. We demonstrate the performance of our method by comparing a plant optimized for different objectives, generation outage durations, minimum power requirements, and power purchase agreements. Although the plant design is sensitive to model parameters and various other assumptions, our results demonstrate some of the optimal designs that occur in different scenarios and what one should expect when designing a hybrid wind-solar-storage power plant.

14 SOLAR ENERGY↗

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗

Open V2X at ESIF: June 2015 - February 2017 (Final Report)

The University of Delaware (UD), and their industrial partners, have carried out testing, documentation, and standards development at ESIF. These activities are intended to help meet the objective of NREL's INTEGRATE program, to enable clean energy technologies to increase the hosting capacity of the grid, by providing grid services in a holistic manner using an open source, interoperable platform. In particular, the 'Open V2X at ESIF' project worked to facilitate standardization and interoperability for Electric Vehicle technologies, including battery electric vehicles and plug-in hybrid vehicles (hence EV) to provide grid services. This was done by interconnecting and testing EV and EVSE (charging station) hardware from different manufacturers, by developing standards and test protocols for standards, and by documenting new protocols to be proposed as standards. Contractor-supplied equipment used in the testing included Vehicle Smart Links (VSLs) installed in an EV and School Bus, single-phase and three-phase charging stations, and an optimizing aggregator (software running on a server, hence just 'Aggregator'). The Aggregator software was developed by U of Delaware and is now licensed to Nuvve Corp for commercial development, and licensed to NRG/EVgo for demonstrations. The Aggregator monitors capacity and operations of the EVs and EVSEs, dispatches them in response to requests, and reports their operations and metered power flow. In operation, as vehicles move, plug-in and unplug, the aggregator tracks the total capacity available within a defined jurisdiction or RTO (Regional Transmission Organization, synonym for ISO or TSO), and dispatches part or all of that capacity based on market or direct dispatch signals from a grid entity. The VSL is the data acquisition and control module in the resource itself (EVSE or EV), with rich communication between a single Optimizing Aggregator and many VSL (one in each resource). All requests and dispatch events are also logged by the aggregator to make the provision of grid power auditable. For these tests, the dispatch was typically a recorded signal, the VSL was in the EV or EVSE, and the recording of experimental performance was often carried out by the normal logging functions of the aggregator.

24 POWER TRANSMISSION AND DISTRIBUTION↗