Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid Analysis”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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↗
DOE OSTI2020

Dynamic Optimal Dispatch of Energy Systems with Intermittent Renewables and Damage Model

With the increasing penetration of intermittent renewable energy sources into the grid, there is a growing need for process systems-based strategies that integrate dispatchable and variable energy systems for supplying the demand while maintaining grid reliability. The proposed framework corresponds to a dynamic mixed-integer linear programming optimization approach that integrates coal-fired and natural gas-fired power plants, NaS batteries for energy storage, and solar/wind energy to supply the demand. This optimization approach considers an economic goal and constraints to provide power balance while maintaining the overall damage of the natural gas combined cycle (NGCC) power plant drum under a maximum stress as well as avoiding the overheating of the NGCC superheater and reheater. Renewable curtailment levels are also retained at minimum levels. Case studies are analyzed considering different loads and renewable penetration levels. The results show that the demand was met for all cases. Grid flexibility was mostly provided by the NGCC, while the batteries were used sparingly. In addition, considering a CO 2 equivalent analysis, the environmental performance was intrinsically connected to grid flexibility and the level of renewable penetration. Stress analysis results reinforced the necessity for an equipment health-related constraint.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Perturbation-Based Diagnosis of False Data Injection Attack Using Distributed Energy Resources

Modern smart grid relies on various sensor measurements for its operational control. In a successful false data injection attack, the attacker manipulates the measurements from the grid sensors such that undetected errors are introduced into the estimates of the system parameters leading to catastrophic situations. This paper proposes a novel perturbation based false data injection attack detection mechanism that utilizes inverter based distributed energy resources (DERs) to create low magnitude perturbation signal in the distribution system voltage that is inconsequential to the normal grid operation. Two voltage sensitivity analysis based algorithms are designed to identify the optimal set of DERs that can create the voltage perturbation signal of desired magnitude. An analytical method of voltage sensitivity analysis is used to compute the magnitude of voltage perturbation signal at each node in a computationally efficient manner. Then, a detection mechanism is developed that checks for the presence of the perturbation sequence in each sensor measurement. A sensor measurement is deemed authentic if the voltage perturbation signal is present in the data. In case of sensor malfunction or cyber-attack, the perturbation signal will not be present in the measurement data. Performance of the proposed attack detection mechanism is validated via simulation of the IEEE 69 bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

FLOW BATTERY COST AND COMMERCIALIZATION TOOL (FlowBaCC) v.1.0

The Flow Battery Cost and Commercialization Tool (FlowBaCC) is a python-based framework used to project the costs and commercialization timelines for redox flow batteries. FlowBaCC combines learning curves for component costs and performance with adoption curves for market diffusion to project future system costs. The tool uses the bottom-up, spreadsheet-based Battery Performance and Cost Model for Flow Batteries (BatPaC-Flow) as the engine to translate learning curve inputs into full system costs. FlowBaCC provides capital costs and levelized costs of energy storage. The tool enables scenario analysis, sensitivity analysis, and uncertainty quantification (via grid and Monte Carlo methods).

Fu, Xiaoxu [Argonne National Laboratory (ANL), Arg↗
DOE OSTI2020

Development of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an Integrated Energy System (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. A new RAVEN (Risk Analysis Virtual ENvironment) plugin for grid and capacity optimization (HERON) has been developed for optimizing the production variables of the IES given the demand profile. Currently, only the limits that affect the production variables and their corresponding time rates of change are considered (explicit constraints). However, other variables are additionally subject to constraints, but the associated limits are not accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response effects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. The Reference Governor (RG) algorithm is traditionally used for enforcing state and control constraints by modifying the set-point trajectories supplied to the feedback regulators. In our application, the RG is coupled within an iterative loop with the HERON-power dispatcher to generate optimal trajectories that ensure the operational constraints are met. A data-driven procedure to derive a representation of the dynamics of the controlled system was developed. First, the variables that represented the state of the system are selected (PCA-based approach), and then state-space representation matrices are derived from the collected measurements (DMDc algorithm). A preliminary version of the developed workflow based on Linear Time Invariant matrices was assessed by adopting a two-unit test case. More sophisticated versions of this workflow foreseeing the on-line derivation of system matrices will be deployed in FY 2021. Finally, a “plug-and-play” library of controllers and state observers was developed in Dymola. Some aspects of the current configuration of the IES unit components, e.g., the encapsulation of the control schemes into dedicated blocks, are consistent with the “plug-and-play” philosophy. Other features, e.g., the system buses collecting the input and the output variables, are not. For this reason, once listed and described the limits of the current configuration, necessary modifications to the plant model interface are presented. As a test-case, the interfaces of the SES model in the RAVEN-based HYBRID framework were reworked accordingly, and two different control schemes were applied to the same plant model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Enhancing Convection Parameterization for Next Generation E3SM (Final Technical Report)

This technical report summarizes the achievements during the project period. The general goal of the project is to enhance the representation of atmospheric deep convection in the DOE global climate model E3SM. By implementing a stochastic convection scheme into E3SM, we improved the simulation of precipitation intensity and frequency. As a result, it has a profound impact on the simulation of other fields in the model, including aerosol scavenging by precipitation and aerosol radiative effects. We also implemented a convective cloud microphysics scheme into E3SMv2 to represent more accurately the interaction between convection and grid-scale clouds. An extensive analysis of model simulations as part of the Next Generation Development (NGD) convection scheme assessment project found that including the convective microphysics scheme improved the simulation of Madden-Julian oscillations (MJO) in the model. To make convective parameterization schemes more scale-aware as the GCM resolution increases, we analyzed the output from cloud-resolving model (CRM) simulations of both organized and unorganized convection by coarse-graining the CRM data to different GCM resolutions. Using a novel approach to modify the convective closure, we demonstrated that the convection scheme used in E3SM (ZM scheme) can be made scale-aware.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Data Validation for Hosting Capacity Analyses [Slides]

The National Renewable Energy Laboratory (NREL), in partnership with the Interstate Renewable Energy Council (IREC), has recently released a report which identifies a suite of best practices for producing trusted, validated hosting capacity analysis results reflecting real-world grid conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Modeled Results of Four Residential Energy Efficiency Measure Packages for Deriving Advanced Building Construction Research Targets

The Advanced Building Construction (ABC) Initiative from the U.S. Department of Energy Building Technologies Office is working to accelerate industrialized construction innovations for decarbonizing buildings. To inform performance and cost targets for research under the ABC Initiative, this analysis used the ResStock™ tool to evaluate the energy savings, utility bill impacts, and carbon emissions impacts of four simulated upgrade packages with specific target performance levels on a large sample of residential dwelling units (approximately 550,000) representative of the U.S. housing stock.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Deploying Adversarial Attacks in Super-Resolution Models

Reliable super-resolution methods are crucial for applications like remote sensing, grid resilience and disaster impact analysis, and standoff biometrics. These methods infuse additional high-frequency information into reconstructions, allowing for better contextualization and image intelligence. However, super-resolution models can also introduce hallucinations or other unseen vulnerabilities that could be exploited by an adversary. This is further compounded by the prominence of deep learning in these models, as models are often blindly applied on out-of-distribution images. In this work, we implement adversarial attacks in common open-source super-resolution models and examine their impact on reconstructions and downstream classification tasks. We find that an adversarially trained super-resolution model can produce high-quality reconstructions that degrade downstream classifications. Moreover, these attacks do not require access to low-resolution imagery or class labels at inference time. These results demonstrate the vulnerability of super-resolution methods to malicious actors and motivates the development of a detector for super-resolution adversarial attacks. Further exploration of adversarial attacks in this domain is required to ensure trustworthiness and robustness of super-resolution models for national security applications.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Estimating and Calibrating DER Model Parameters Using Levenberg–Marquardt Algorithm in Renewable Rich Power Grid

The proliferation of inverter-based distributed energy resources (IBDERs) has increased the number of control variables and dynamic interactions, leading to new grid control challenges. For stability analysis and designing appropriate protection controls, it is important that IBDER models are accurate. This paper focuses on the accurate estimation and parameter calibration of DER_A, a recently proposed aggregated IBDER model. In particular, we focus on the parameters of the reactive power–voltage regulation module. We formulate the problem of parameter tuning as a non-linear least square minimization problem and solve it using the Levenberg–Marquardt (LM) method. The LM method is primarily chosen due to its flexibility in adaptively selecting between the steepest descent and Gauss–Newton methods through a damping parameter. The LM approach is used to minimize the error between the actual measurements and the estimated response of the model. Further, the computational challenges posed by the numerical calculation of the Jacobian are tackled using a quasi-Newton root-finding approach. The proposed method is validated on a real feeder model in the northeastern part of the United States. The feeder is modeled in OpenDSS and the measurements thus obtained are fed to the DER_A model for calibration. The simulation results indicate that our approach is able to successfully calibrate the relevant model parameters quickly and with high accuracy, with a total sum of square error of 3.57 × 10 –7 .

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ~$\frac{1}{3}$ of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2% for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

Worldwide cloud cover model

Classifying worldwide cloudiness into homogeneous regions, using a satellite data set containing day IR, night IR, incoming, and absorbed solar radiation measurements on a 2.5-degree latitude-longitude grid is considered. Methods of analysis are presented.

Smith, O. E.↗
NASA NTRS1981

Low-level circulations induced by an upper-level jet in a fine mesh numerical model

Attention is given to varying the initial conditions in the Penn State three-dimensional numerical model in the development of the low-level jet over the Ohio Valley during the 10 May 1973 severe weather event. The results of two forecasts are summarized: (1) the control forecast which interpolated the National Meteorological Center's operational analysis to the mesoscale model grid; and (2) modifying the control forecast by inserting artificial soundings and surface observations. Enhancements are made on the basis of subsynoptic scale analyses by Uccellini and surface data. The results indicate that the low-level jet was in part induced by pressure changes associated with deceleration and horizontal divergence of air in the exit region of the upper-level jet. The enhancement of the initial moisture field using surface reports of clouds and precipitation has a substantial effect on precipitation.

Burkhart, R. P.↗
NASA NTRS1987

Aerodynamic design modification of a hypersonic wind tunnel nozzle by CSCM with high order accuracy

An improved version of the CSCM implicit Navier-Stokes solver with a flexible data structure has been applied, in conjunction with a sophisticated patched grid system, to the aerodynamic analysis of hypersonic axisymmetric contoured nozzles of the NASA-Ames 3.5-foot hypersonic wind tunnel. The results obtained are found to be in good agreement with experiment. Based on the analysis, a newly designed throat contour for the Mach 14 nozzle has been recommended.

Yang, J. Y.↗
NASA NTRS1992

Three-dimensional time-marching aeroelastic analyses using an unstructured-grid Euler method

Modifications to a three dimensional, implicit, upwind, unstructured-grid Euler code for aeroelastic analysis of complete aircraft configurations are described. The modifications involve the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. The paper presents a detailed description of the time marching aeroelastic procedure and presents comparisons with experimental data to provide an assessment of the capability. Flutter results are shown for an isolated 45 degree swept-back wing and a supersonic transport configuration with a fuselage, clipped delta wing, and two identical rearward-mounted nacelles. Comparisons between computed and experimental flutter characteristics show good agreement, giving confidence in the accuracy of the aeroelastic capability that was developed.

Rausch, Russ D.↗
NASA NTRS1992

Three-dimensional time-marching aeroelastic analyses using an unstructured-grid Euler method

The modifications to a 3D implicit, upwind, unstructured-grid Euler code for aeroelastic analysis of complete aircraft configurations described involve the addition of structural equations of motion for their simultaneous time integration with the governing flow equations. A detailed description of the time-marching aeroelastic procedure and comparisons with experimental data are presented, to provide an assessment of capabilities. Flutter results are shown for an isolated 45-deg sweptback wing and an SST with a fuselage, clipped delta wing, and two identical rearward-mounted nacelles. Comparisons between computed and experimental flutter characteristics show good agreement.

Rausch, Russ D.↗
NASA NTRS1992

Determination and representation of electric charge distributions associated with adverse weather conditions

Algorithms are presented for determining the size and location of electric charges which model storm systems and lightning strikes. The analysis utilizes readings from a grid of ground level field mills and geometric constraints on parameters to arrive at a representative set of charges. This set is used to generate three dimensional graphical depictions of the set as well as contour maps of the ground level electrical environment over the grid. The composite, analytic and graphic package is demonstrated and evaluated using controlled input data and archived data from a storm system. The results demonstrate the packages utility as: an operational tool in appraising adverse weather conditions; a research tool in studies of topics such as storm structure, storm dynamics, and lightning; and a tool in designing and evaluating grid systems.

Rompala, John T.↗
NASA NTRS1994

Determination of ocean tides from the first year of TOPEX/POSEIDON altimeter measurements

An improved geocentric global ocean tide model has been determined using 1 year of TOPEX/POSEIDON altimeter measurements to provide corrections to the Cartwright and Ray (1991) model (CR91). The corrections were determined on a 3 deg x 3 deg grid using both the harmonic analysis method and the response method. The two approaches produce similar solutions. The effect on the tide solution of simultaneously adjusting radial orbit correction parameters using altimeter measurements was examined. Four semidiurnal (N(sub 2), M(sub 2), S(sub 2) and K(sub 2)), four diurnal (Q(sdub 1), O(sub 1), P(sub 1), and K(sub 1)), and three long-period (S(sub sa), M(sub m), and M(sub f)) constituents, along with the variations at the annual frequency, were included in the harmomnic analysis solution. The observed annual variations represents the first global measurement describing accurate seasonal changes of the ocean during an El Nino year. The corrections to the M(sub 2) constituent have an root mean square (RMS) of 3.6 cm and display a clear banding pattern with regional highs and lows reaching 8 cm. The improved tide model reduces the weighted altimeter crossover residual from 9.8 cm RMS, when the CR91 tide model is used, to 8.2 cm on RMS. Comparison of the improved model to pelagic tidal constants determined from 80 tide gauges gives RMS differences of 2.7 cm for M(sub 2) and 1.7 cm for K(sub 1). Comparable values when the CR91 model is used are 3.9 cm and 2.0 cm, respectively. Examination of TOPEX/POSEIDON sea level anomaly variations using the new tide model further confirms that the tide model has been improved.

Ma, X. C.↗