Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid computing”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computational and numerical analysis of AC optimal power flow formulations on large-scale power grids

We report that alternating current optimal power flow (AC-OPF) is a fundamental tool in electric utilities to determine optimal operation of the various resources. Typically, the AC-OPF problem uses power balance formulation containing voltages and power equations. Yet, there is no comprehensive comparison of the different AC-OPF formulations, especially for large-scale networks. This paper presents a detailed comparative evaluation of different formulations of the AC-OPF problem on networks ranging from 9-bus to 25,000 buses. Three different formulations: 1) power balance with polar voltages, 2) power balance with Cartesian voltages, and 3) current balance with Cartesian voltages are discussed in detail by comparing their characteristics, and numerical and computational performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

A Privacy-Aware Federated Learning Framework for Distributed Energy Resource Analytics in Constrained Environments

To be resilient against extreme weather events, the rural communities in Puerto Rico are leveraging distributed energy resources (DER). However, computing frameworks sup-porting the grid in critical decision-making are still largely centralized. Sensitive consumer data are transmitted over the Internet or cellular networks to a secondary or tertiary node. It guarantees better situational awareness at the cost of a wider attack surface, jeopardizing user privacy, as more DER come online. Cloud, Edge, and Fog computing all require data aggregation at some level. This paper introduces a privacy-aware federated learning framework that leverages the Fog model by pushing analytics all the way to the DER and load assets. These local models train on individual asset data and transmit only learned parameters (such as weights) over secure communications to a global decision-maker. By abstracting personally identifiable consumer data without impacting decision optimality, this framework better aligns with distributed power generation paradigm.

Sundararajan, Aditya↗
DOE OSTI2024

Sienna Modeling Framework [Slides]

NREL's Sienna modeling framework effectively builds, solves, and analyzes the scheduling problems and dynamic simulations of quasi-static infrastructure systems. It uses a modular framework to answer different questions about future energy systems, fundamentally advancing the nation's ability to model individual and integrated infrastructure systems at a range of spatial and temporal scales. This presentation will include NREL power grid researcher Clayton Barrows.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

CFD-Informed Spacer Grid Models in the Subchannel Code CTF

Data from computational fluid dynamics (CFD) simulations can be utilized to improve the predictive capabilities of the subchannel code CTF. Past work incorporating CFD results into CTF has included accounting for directed crossflow and turbulent mixing between adjacent subchannels. Spacer grids impact fluid flow in various ways, and this work focuses on four spacer grid effects, namely: 1. Directed crossflow; 2. Turbulent mixing; 3. Enhanced downstream heat transfer; 4. Pressure drop correlation. Heat transfer augmentation occurs partially due to the disturbances in the developed flow and from the additional turbulence generated by the presence of obstructions in the fluid flow's path. The boundary layer disruption is caused in part by the newly-directed flow just downstream of the spacer grid vanes. This effect lessens as the distance from the grid increases and the flow returns to its developed state, closer to what would be expected in a bare subchannel. Current plans of this work are to develop a method for incorporating the effects of spacer grids on the heat transfer calculation within CTF. Proposed techniques include the addition of an enhancement factor to the local heat transfer coefficient. The inclusion of CFD data from properly-designed simulations will allow CTF to better predict parameters such as fluid temperatures near spacer grids. Improved pressure loss coefficients for each grid can also be generated using CFD simulations. The heat transfer modification has been studied, with other authors noting correlations which are dependent on grid properties such as the blockage ratio. Parameters considered in this study are the turbulent kinetic energy (TKE) and vorticity. There are currently no models to calculate or store either of these parameters within CTF. When accounting for spacer grid effects in CTF, it is important to not double-count the effects towards two different CTF models. This requires carefully designing CFD simulations around capturing specific phenomena from the grids and being able to apply these CFD results to all appropriate flow conditions. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Fast-Tracking Technology Solutions to Boost Cybersecurity in Virtual Power Plant and Aggregator Operations

This presentation was presented at the Aggregation and Grid Security Workshop - held on June 17-18, 2025, at NREL in Golden, Colorado. The goal of the two-day workshop was to address the critical cybersecurity challenges for the future electric grid. This presentation overviews NREL's work in evaluating cybersecurity technology solutions for virtual power plant and aggregator operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Multi-level optimization with the koopman operator for data-driven, domain-aware, and dynamic system security

Cyber-Physical Systems (CPSs) like the power grid are critically important but also increasingly vulnerable; ensuring reliable system operation in the face of disruptions is becoming more and more challenging. Multi-Level Optimization (MLO) is a powerful way to model adversarial interactions, which naturally makes it applicable to studying CPS security. However, MLO typically does not address underlying system dynamics, and incorporating nonlinear dynamics is generally infeasible. In this paper, we show how to combine MLO with the Koopman Operator (KO) to remedy this. The KO maps nonlinear dynamics to a lifted space in which those dynamics are linear, thus making it ideal for use with MLO. Moreover, the structure of the KO also provides convenient ways to incorporate domain knowledge into the data-driven process of learning the KO representation of a given system. Here we then demonstrate the use of MLO-KO on a small example problem taken from the power grid domain, discuss the scalability and computational cost of MLO-KO, and identify future research directions for this work.

42 ENGINEERING↗
DOE OSTI2022

Model Predictive Control for a Grid-interactive Efficient Thermal Storage-integrated Heat Pump System

Building heating and cooling systems can be used to overcome the mismatch between the intermittent supply of renewable power and the fluctuating demand for electricity. A novel underground thermal energy storage integrated with a dual-source heat pump has been proposed to mitigate the mismatch while meeting the thermal demand of buildings efficiently. Conventional thermostat control with heuristic rules cannot provide intelligent decisions to maximize the thermal efficiency and flexibility of the proposed system. Advanced control strategies like model predictive control (MPC) have provided a new paradigm for grid-interactive efficient building operation with the advancement of computation and sensing. This study developed an MPC for the proposed system to provide grid service for Demand Side Management and minimize the operating cost of building owners. A control-oriented dynamic model of the proposed system has been developed. Given an objective function and proper constraints, an optimization problem is formulated to determine the optimal control strategy of the system. Dynamic Programming is adopted to solve the optimization problem. A rule-based control (RBC) is also developed to achieve similar goals. Short-term simulations are conducted to compare the system performance resulting from the two controls. The simulation results indicate that the MPC performs more intelligently than the RBC in charging thermal energy storage and selecting heat pump sources by taking advantage of the predicted cooling demands of the building and the performance of the integrated system. As a result, the MPC could save energy and reduce operating costs compared with the RBC. A case study shows that, for a 3-day operation, the MPC saves 36.9% energy and reduces 38.5% operating cost compared with the RBC.

Shi, Liang↗
DOE OSTI2023

Metrics for Evaluating Grid Service Provision from Communities of Grid-interactive and Efficient Buildings and other DER

Tapping into the flexibility of aggregations of Grid-interactive and Efficiency Buildings (GEBs) represent a large opportunity to cost-effectively improve operations in future low-carbon power grids. The US Department of Energy's Connected Communities program seeks to demonstrate collections of these GEBs at 10 different sites across the country. However, consistent metrics for evaluation across all projects are needed to build confidence in the approach. This paper presents the metrics relating to Grid Service Provision that will be computed at the 10 demonstrations. The metrics proposed quantity the magnitude of service offers, the consistency and quality of services provided by the community, and the individual DER contributions to the community-level service. The goal of this work is to present these metrics and describe some of their intended insights. The hope is that these insights will be useful and build confidence for grid operators, regulators, and aggregators and practitioners as they look to deploy these resources in grids of the future.

MacDonald, Jason S↗
DOE OSTI2018

A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm

The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. A novel, implicit, particle-in-cell based scheme for the hybrid model is derived for fully 3D electromagnetic problems with multiple ion species, which features global mass, momentum and energy conservation. The scheme includes sub-cycling and orbit-averaging for the ions, with cell-centered finite differences and implicit midpoint time advance. To reduce discrete particle noise, the scheme allows arbitrary-order shape functions for the particle-mesh interpolations and the application of conservative binomial smoothing. The algorithm is verified for a number of test problems to demonstrate the correctness of the implementation, the unique conservation properties, and the favorable stability properties of the new scheme. In particular, there is no indication of unstable growth of the finite-grid instability for a population of cold ions drifting through a uniform spatial mesh, in a set-up where several commonly used non-conservative schemes are highly unstable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A block-tridiagonal solver with two-level parallelization for finite element-spectral codes

Two-level parallelization is introduced to solve a massive block-tridiagonal matrix system. One-level is used for distributing blocks whose size is as large as the number of block rows due to the spectral basis, and the other level is used for parallelizing in the block row dimension. The reason for the added parallelization dimension is to retard the saturation of the scaling due to communication overhead and inefficiencies in the single-level parallelization only distributing blocks. As a technique for parallelizing the tridiagonal matrix, the combined method of “Partitioned Thomas method” and “Cyclic Odd–Even Reduction” is implemented in an MPI-Fortran90 based finite element-spectral code (TORIC) that calculates the propagation of electromagnetic waves in a tokamak. The two-level parallel solver using thousands of processors shows more than 5 times improved computation speed with the optimized processor grid compared to the single-level parallel solver under the same conditions. Three-dimensional RF field reconstructions in a tokamak are reflected as examples of the physics simulations that have been enabled by this algorithmic advance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Interleaved Cuk Converter Wave Energy System With Advanced Control and Grid Support Functions for Seamless Integration

This paper presents an innovative wave energy conversion system that integrates an interleaved Cuk converter with advanced nonlinear control for seamless grid integration. The system efficiently extracts power on the DC side using the interleaved Cuk converter, while an inverter manages power transfer to the grid or load on the AC side. A nonlinear control architecture, based on the Lyapunov energy function, ensures stable and optimal operation under varying wave conditions, effectively addressing the challenges of variability and unpredictability inherent in wave energy. The proposed system also incorporates features to enhance power quality and minimize losses, making it a robust solution for renewable energy integration. The system's effectiveness in harnessing wave energy and achieving smooth grid integration is validated through comprehensive computer simulations in MATLAB/Simulink, with case study results demonstrating its capabilities.

HYDRO ENERGY,TIDAL AND WAVE POWER↗
DOE OSTI2024

Impact of Reordering on the LU Factorization Performance of Bordered Block-Diagonal Sparse Matrix

Power engineers rely on computer-based simulation tools to assess grid performance and ensure security. At the core of these tools are solvers for sparse linear equations. When transformed into a bordered block-diagonal (BBD) structure, part of the sparse linear equation solving can be parallelized. This work focuses on using the Schur-complement-based method for LU factorization on BBD matrices, specifically, Jacobian matrices from large-scale systems. Our findings show that the natural ordering method outperforms the default ordering method in computational performance for each block of the BBD matrix. This observation is validated using synthetic 25k-bus and 70k-bus cases, showing a speedup of up to 38% when using natural ordering without permutation. Additionally, the impact of the number of partitions is studied, and the result shows that computational performance improves with more, smaller partitions in the BBD matrices.

BBD matrix↗
DOE OSTI2020

A CIERRA Gridded Product Climatology for LIS/OTD [Slides]

As a stepping-stone towards developing improved GLM-CIERRA grids, I am testing pixel matching techniques on the LIS/OTD data; The goal is to maintain accuracy while vastly improving computational efficiency; The end result will be standardized grids at the nominal resolution of each instrument (5 km for LIS, 10 km for OTD) that take into account the pixel geometry for each event that comprises a given flash

58 GEOSCIENCES↗
DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗
DOE OSTI2023

A Gradient Based Subgrid-Scale Parameterization for Ocean Mesoscale Eddies

Mesoscale eddies play an important role in transport of heat and biogeochemical tracers in the global ocean circulation. Resolving these energetic eddies, however, is challenging in ocean general circulation models (OGCM) because it requires a horizontal grid spacing of ≲1/8° that is computationally expensive. As a result, we are required to parameterize mesoscale eddy effects on large-scale ocean flows. In this work, we introduce a new subgrid-scale (SGS) model that is developed based on a Taylor series expansion of resolved variables to parameterize subgrid mesoscale eddy transports and momentum fluxes in OGCM. We have performed an a priori study to evaluate the performance of our new gradient model using high-resolution ocean simulations. Our results show that the gradient model well predicts the actual SGS thickness fluxes in the zonal and meridional directions in coarse-resolution simulations with the grid spacing ≳1/4°. The unresolved kinetic energy at the ocean surface is also skillfully estimated. More importantly, unlike current mesoscale eddy parameterizations, which are mainly developed based on an assumption of flat bottom ocean, our new SGS model can capture the structure of unresolved standing meanders at the ocean surface. We have also developed a dynamic procedure for setting in non-dimensional parameters in our new parameterization through a non-ad hoc and tuning-free method. Overall, this work suggests that implementing the gradient model in OGCM can improve the model accuracy with an affordable computational cost in eddy-permitting and non-eddying simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A Physical Model Enhanced Data Driven Method for High-Resolution Residential Load Profile Generation

Residential buildings account for significant energy consumption, creating opportunities to offer grid services. As electric utilities seek to implement effective system operation strategies, understanding residential energy consumption patterns becomes essential; However, the time intervals of load profiles measured by utilities' smart meters are typically from 15 minutes to 60 minutes. The low-resolution data make it hard to extract appliance-level load information, which is critical for providing grid services. This paper presents a load profile generator designed to produce synthetic load profiles for residential buildings that emphasizes the importance of accurate representations of realistic energy consumption patterns. The generator takes realistic low-resolution residential load measurements and weather data as inputs, producing 1-minute interval profiles that match the characteristics of the original profiles. Further, this generator can be used to populate load profiles in areas where actual measurements are limited to improve the ability of utilities to analyze their distribution systems. By providing more high-resolution residential building load profiles, this tool supports electric utilities to enhance their residential building load control strategies and improve overall grid stability.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

A framework for WRF to WRF-IBM grid nesting to enable multiscale simulations

With advances in computational power, mesoscale models, such as the Weather Research and Forecasting (WRF) model, are often pushed to higher resolutions. As the model’s horizontal resolution is refined, the maximum resolved terrain slope will increase. Because WRF uses a terrain-following coordinate, this increase in resolved terrain slopes introduces additional grid skewness. At high resolutions and over complex terrain, this grid skewness can introduce large numerical errors that require methods, such as the immersed boundary method, to keep the model accurate and stable. Our implementation of the immersed boundary method in the WRF model, WRF-IBM, has proven effective at microscale simulations over complex terrain. WRF-IBM uses a non-conforming grid that extends beneath the model’s terrain. Boundary conditions at the immersed boundary, the terrain, are enforced by introducing a body force term to the governing equations at points directly beneath the immersed boundary. Nesting between a WRF parent grid and a WRF-IBM child grid requires a new framework for initialization and forcing of the child WRF-IBM grid. This framework will enable concurrent multi-scale simulations within the WRF model, improving the accuracy of high-resolution simulations and enabling simulations across a wide range of scales.

54 ENVIRONMENTAL SCIENCES↗