SEARCH · Engineering Papers

Results for “Structured Grids”

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.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Average Incremental Cost Pricing for the AC Unit Commitment Problem [Rev. 1]

Unit Commitment (UC) problems that consider the Alternating Current (AC) model of the transmission network have long been considered intractable to solve at scale by the power system community. Recently, the Grid-Optimization (GO) Competition held by the Advanced Research Project Agency-Energy (ARPA-E) has facilitated the development of the first algorithms to solve large-scale ACUC problems. This new capability opens a path towards the explicit consideration of the AC transmission network model in UC problems used to clear day-ahead electricity markets. This calls for the analysis of electricity market structures that accommodate both the continuous non-linearity of the AC transmission network and the discrete non-linearity of the UC problem simultaneously. This paper serves as an initial effort to do so by proposing an Average Incremental Cost (AIC) pricing structure that is designed around the ACUC problem. In particular, an AIC one-pass pricing problem is proposed that represents a continuously constrained variant of the ACUC problem and allows for the computation of Locational Incremental Prices (LIPs) for both real and reactive power as the local optimal Lagrange multipliers of the power balance constraints. To avoid degeneracy, the pricing problem includes a small parameter ϵ > 0. Under certain assumptions market participants are shown to realize profit that converges to a non-negative value as ϵ approaches zero, practically ensuring profitability for small values of ϵ. We additionally provide many simple and important examples that provide intuition and insights into the proposed prices. Examples illustrate the basic concept of AIC pricing, the derived profitability results, the existence of multiple LIPs, the importance of including reactive power in the dispatch and pricing problems, the need for reactive power prices, and the improved incentives exhibited by LIPs as compared to traditional Locational Marginal Prices (LMPs). We additionally indicate many directions for future work including analysis of larger test cases.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2023

An Open-Source Parallel EMT Simulation Framework: Preprint

As the integration level of inverter-based resources (IBR) increases, ensuring the reliable operation of the bulk power systems requires the use of electromagnetic transient (EMT) simulation tools to identify and mitigate system-wide stability risks. Conducting EMT studies for large-scale, IBR-rich grids, however, is challenging due to the inherent computational bottleneck caused by the underlying high-fidelity models and required small time steps. This paper introduces ParaEMT: an open-source, generic EMT simulation framework designed to accelerate simulations by leveraging advanced parallel computational technologies, such as high-performance computers. This paper presents a comprehensive exposition of ParaEMT, covering its modeling library, simulation strategy, framework structure, operational procedures, and auxiliary features, alongside its extensible parallel computational architecture. Notably, ParaEMT is a publicly accessible and modularized framework written in Python, thereby facilitating future development and the integration of new models and algorithms. The accuracy and efficiency of ParaEMT are demonstrated by rigorous validations via multiple case studies.

electromagnetic transient simulation↗

DOE OSTI2024

An Accelerated Clip Algorithm for Unstructured Meshes: A Batch-Driven Approach

The clip technique is a popular method for visualizing complex structures and phenomena within 3D unstructured meshes. Meshes can be clipped by specifying a scalar isovalue to produce an output unstructured mesh with its external surface as the isovalue. Similar to isocontouring, the clipping process relies on scalar data associated with the mesh points, including scalar data generated by implicit functions such as planes, boxes, and spheres, which facilitates the visualization of results interior to the grid. In this paper, we introduce a novel batch-driven parallel algorithm based on a sequential clip algorithm designed for high-quality results in partial volume extraction. Our algorithm comprises five passes, each progressively processing data to generate the resulting clipped unstructured mesh. The novelty lies in the use of fixed-size batches of points and cells, which enable rapid workload trimming and parallel processing, leading to a significantly improved memory footprint and run-time performance compared to the original version. On a 32-core CPU, the proposed batch-driven parallel algorithm demonstrates a run-time speed-up of up to 32.6x and a memory footprint reduction of up to 4.37x compared to the existing sequential algorithm. The software is currently available under an open-source license in the VTK visualization system.

Tsalikis, Spiros↗

DOE OSTI2022

Advanced Hydropower and PSH Capacity Expansion Modeling (Final Report on HydroWIRES D1 Improvements to Capacity Expansion Modeling)

Hydropower and pumped-storage hydropower (PSH) have played a key role in providing flexible, low-carbon electricity to the U.S. electricity system for over a century. As variable generation (VG) deployment increases the demand for flexible, dispatchable generation, it is important to use all available methods for understanding how hydropower and PSH interfaces with VG in a future low-carbon grid. Capacity expansion models (CEMs) of electricity systems are often used to study future electricity scenarios, but these tools often have difficulty representing site- and technology-details of hydropower and PSH due to limited spatial, temporal, or process resolution. This report demonstrates a set of model advancements to improve hydropower and PSH representations in CEMs and other models that consider hydropower and PSH's role in electricity systems. Model advancements include new data integration to define closed-loop PSH resource availability and cost, along with site-level parameterizations of existing PSH capacity and energy storage specifications. Shifting energy across seasons for both hydropower and PSH is explored to demonstrate the potential value of long-duration storage. New representations of hydropower upgrades offer opportunities to increase dispatchability, add pumps, or independently add capacity or energy depending on what is most valuable. New model structures and data are demonstrated individually and in combination to observe their impact on model results and provide initial insights into what is most important for analysts, electricity system planners, and hydropower decision-makers to consider when assessing future roles of hydropower. Improving flexibility of existing hydropower assets has the potential to reduce CO 2 emissions through complementing VG technologies and/or improve electricity system economics by reducing the need to deploy higher-cost systems. Large-scale energy storage is also shown to provide opportunities for balancing seasonal differences in VG, particularly solar photovoltaics (PV). Initial analysis with a new closed-loop PSH resource and cost dataset also demonstrate the potential for new PSH deployment to offer new grid flexibility opportunities. Future analysis and modeling of hydropower's role in the U.S. and other electricity systems could include methodological improvements based on those discussed here along with sensitivity analysis to better represent current and future potential hydropower and PSH flexibility. Improvements to geospatial and site-level data can further improve the understanding of how hydropower and PSH can be upgraded and operated in response to future electricity system needs. This work lays the groundwork to use advanced planning models to explore the range of roles hydropower and PSH can play in the grid of the future.

13 HYDRO ENERGY↗

DOE OSTI2024

Inverter Model Validation and Calibration Using Phasor Measurement Unit Data

As the penetration of inverter-based renewable energy resources increases in the power grid, especially at the distribution and microgrid levels, the need to accurately represent them in planning studies increases as well. However, due to the lack of well-established standard procedures, and vendor reluctance towards the detailed sharing of proprietary models, automated dynamic model validation and parameter calibration tools for inverter based resources (IBRs) remain scarce. This work presents a model validation and parameter calibration platform for representing IBRs with generic phasor-domain models. Phasor measurements of power system events are used for continuous validation using the data playback method, and model parameters are re-calibrated if a significant mismatch between measurements and model response is observed. Unique features of the proposed platform include- (a) an iterative Bayesian optimization approach towards parameter calibration to address a possible mismatch between the structures of generic models implemented in simulation softwares and actual commercial inverters, (b) error metrics designed to account for a possible mismatch between the time resolution of simulation and measurements, and (c) analysis of the measurement-simulation mismatch to provide guidance to engineering personnel regarding model shortcomings. The performance of the platform has been illustrated using both simulated data and field measurements to validate/calibrate inverter models in GridLAB-D.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Oscillatory Spreading and Inertia in Power Grids

The increase in variable renewable generators (VRGs) in power systems has altered the dynamics from a historical experience. VRGs introduce new sources of power oscillations, and the stabilizing response provided by synchronous generators (SGs, e.g., natural gas, coal, etc.), which help avoid some power fluctuations, will lessen as VRGs replace SGs. These changes have led to the need for new methods and metrics to quickly assess the likely oscillatory behavior for a particular network without performing computationally expensive simulations. This work studies the impact of a critical dynamical parameter - the inertia value - on the rest of a power system's oscillatory response to representative VRG perturbations. We use a known localization metric in a novel way to quantify the number of nodes responding to a perturbation and the magnitude of those responses. This metric allows us to relate the spread and severity of a system's power oscillations with inertia. We find that as inertia increases, the system response to node perturbations transitions from localized (only a few close nodes respond) to delocalized (many nodes across the network respond). We introduce a heuristic computed from the network Laplacian to relate this oscillatory transition to the network structure. We show that our heuristic accurately describes the spread of oscillations for a realistic power-system test case. Using a heuristic to determine the likely oscillatory behavior of a system given a set of parameters has wide applicability in power systems, and it could decrease the computational workload of planning and operation.

dynamical systems↗

DOE OSTI2024

Sustainable layered cathode with suppressed phase transition for long-life sodium-ion batteries

Sodium-ion batteries are among the most promising alternatives to lithium-based technologies for grid and other energy storage applications due to their cost benefits and sustainable resource supply. For the cathode—the component that largely determines the energy density of a sodium-ion battery cell—one major category of materials is P2-type layered oxides. Unfortunately, at high state-of-charge, such materials tend to undergo a phase transition with a very large volume change and consequent structural degradation during long-term cycling. Here we address this issue by introducing vacancies into the transition metal layer of P2-Na$_{0.7}$Fe$_{0.1}$Mn$_{0.75\square}$$_{0.15}$O2 (‘$_{\square}$’ represents a vacancy). The transition metal vacancy serves to suppress migration of neighbouring Na ions and therefore maintain structural and thermal stability in Na-depleted states. Moreover, the specific Na-O-$_{\square}$ configuration triggers a reversible anionic redox reaction and boosts the energy density. As a result, the cathode design here enables pouch cells with energy densities of 170 Wh kg -1 and 120 Wh kg -1 that can operate for over 600 and 1,000 cycles, respectively. In conclusion, our work not only suggests a feasible strategy for cathode design but also confirms the possibility of developing a battery chemistry that features a reduced need for critical raw materials.

25 ENERGY STORAGE↗

DOE OSTI2024

Time-resolved cryogenic electron tomography for the study of transient cellular processes

Cryogenic electron tomography (cryo-ET) is the highest resolution imaging technique applicable to the life sciences, enabling sub-nanometer visualization of specimens preserved in their near native states. The rapid plunge freezing process used to prepare samples lends itself to time-resolved studies, which researchers have pursued for in vitro samples for decades. Here, we focus on developing a freezing apparatus for time-resolved studies in situ. The device mixes cellular samples with solution-phase stimulants before spraying them directly onto an electron microscopy grid that is transiting into cryogenic liquid ethane. By varying the flow rates of cell and stimulant solutions within the device, we can control the reaction time from tens of milliseconds to over a second prior to freezing. In a proof-of-principle demonstration, the freezing method is applied to a model bacterium, Caulobacter crescentus, mixed with an acidic buffer. Through cryo-ET we resolved structural changes throughout the cell, including surface-layer protein dissolution, outer membrane deformation, and cytosolic rearrangement, all within 1.5 seconds of reaction time. This new approach, Time-Resolved cryo-ET (TR-cryo-ET), enhances the capabilities of cryo-ET by incorporating a sub-second temporal axis and enables the visualization of induced structural changes at the molecular, organelle, or cellular level.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

SMC2021 Where to go in atomic world

'Graphene_CrSi.npy' contains a scanning transmission electron microscopy (STEM) movie from graphene monolayer. The movie is a sequence of atom-resolved images from the same sample region that undergoes chemical and structural transformations due to interaction with electron beam (which is used to perform imaging). 'topo_defects.npy' contains coordinates of some of the objects of interest (topological defects in graphene). It is a dictionary, where keys are frame numbers and values are xy coordinates. These objects are usually localized in relatively small areas of the image and we are interested in identifying them without having to scan an entire grid (which leads to fast degradation of the sample).

36 MATERIALS SCIENCE↗

DOE OSTI2024

Machine Learning-Based Identification of the Interface Regions for Coupling Local and Nonlocal Models

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. However, the coupling process can be challenging, requiring expertise to identify the interface between local and nonlocal regions. Here, this study introduces a machine learning-based approach to automatically detect the regions in which the local and nonlocal models should be used. The method uses loading functions evaluated at grid points to decide the model selection at those points. Training of the networks is based on datasets provided by classes of loading functions for which reference coupling configurations are computed using accurate coupled solutions, where accuracy is measured in terms of the relative error between the solution to the coupling approach and the solution to the nonlocal model. We study two approaches that vary in data structure. The first, the full-domain input data approach, uses the entire load vector and outputs a complete label vector, performing a global classification. The second, a window-based approach, processes loads into windows and addresses the problem as a node-wise classification where each window's central point is classified individually. The classification problems are solved via deep learning algorithms based on convolutional neural networks. The performance of these approaches is studied on one-dimensional numerical examples using F1-scores and accuracy metrics. Notably, the windowing approach achieves an accuracy of 0.96 and an F1-score of 0.97, highlighting its potential to automate coupling processes effectively and enhance computational efficiency in material science applications.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A Consistent Representation of Cloud Overlap and Cloud Subgrid Vertical Heterogeneity

Many global climate models underestimate the cloud cover and overestimate the cloud albedo, especially for low-level clouds. We determine how a correct representation of the vertical structure of clouds can fix part of this bias. We use the 1D McICA framework and focus on low-level clouds. Using Large Eddy Simulations results as reference, we propose a method based on exponential-random overlap that represents the cloud overlap between layers and the subgrid cloud properties over several vertical scales, with a single value of the overlap parameter. Starting from a coarse vertical grid, representative of atmospheric models, this algorithm is used to generate the vertical profile of the cloud fraction with a finer vertical resolution, or to generate it on the coarse grid but with subgrid heterogeneity and cloud overlap that ensures a correct cloud cover. Doing so we find decorrelation lengths are dependent on the vertical resolution, except if the vertical subgrid heterogeneity and interlayer overlap are taken into account coherently. We confirm that the frequently used maximum-random overlap leads to a significant error by underestimating the low-level cloud cover with a relative error of about 50%, that can lead to an error of SW cloud albedo as big as 70%. Not taking into account the subgrid vertical heterogeneity of clouds can cause a relative error of 20% in brightness, assuming the cloud cover is correct.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{Lagrangian LES}, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗

DOE OSTI2020

A Review of Sandia Energy Storage Research Capabilities and Opportunities (2020 to 2030)

Large-scale integration of energy storage on the electric grid will be essential to enabling greater penetration of intermittent renewable energy sources, modernizing the grid for increased flexibility security, reliability, and resilience, and enabling cleaner forms of transportation. The purpose of this report is to summarize Sandia's research and capabilities in energy storage and to provide a preliminary roadmap for future efforts in this area that can address the ongoing program needs of DOE and the nation. Mission and vision statements are first presented followed by an overview of the organizational structure at Sandia that provides support and activities in energy storage. Then, a summary of Sandia's energy storage capabilities is presented by technology, including battery storage and materials, power conversion and electronics, subsurface-based energy storage, thermal/thermochemical energy storage, hydrogen storage, data analytics/systems optimization/controls, safety of energy storage systems, and testing/demonstrations/model validation. A summary of identified gaps and needs is also presented for each technology and capability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

WRF-simulated low-level jets over Iowa: characterization and sensitivity studies

Abstract. Output from 6 months of high-resolution simulations with the Weather Research and Forecasting (WRF) model are analyzed to characterize local low-level jets (LLJs) over Iowa for winter and spring in the contemporary climate. Low-level jets affect rotor plane aerodynamic loading, turbine structural loading and turbine performance, and thus accurate characterization and identification are pertinent. Analyses using a detection algorithm wherein the wind speed above and below the jet maximum must be below 80 % of the jet wind speed within a vertical window of approximately 20–530 m a.g.l. (above ground level) indicate the presence of an LLJ in at least one of the 14 700 4 km×4 km grid cells over Iowa on 98 % of nights. Nocturnal LLJs are most frequently associated with stable stratification and low turbulent kinetic energy (TKE) and hence are more frequent during the winter months. The spatiotemporal mean LLJ maximum (jet core) wind speed is 9.55 m s−1, and the mean height is 182 m. Locations of high LLJ frequency and duration across the state are seasonally varying, with a mean duration of 3.5 h. The highest frequency occurs in the topographically complex northwest of the state in winter and in the flatter northeast of the state in spring. Sensitivity of LLJ characteristics to the (i) LLJ definition and (ii) vertical resolution at which the WRF output is sampled is examined. LLJ definitions commonly used in the literature are considered in the first sensitivity analysis. These sensitivity analyses indicate that LLJ characteristics are highly variable with definition. Use of different definitions identifies both different frequencies of LLJs and different LLJ events. Further, when the model output is down-sampled to lower vertical resolution, the mean jet core wind speed height decreases, but spatial distributions of regions of high frequency and duration are conserved. Implementation of a polynomial interpolation to extrapolate down-sampled output to full-resolution results in reduced sensitivity of LLJ characteristics to down-sampling.

17 WIND ENERGY↗

DOE OSTI2021

Failure prediction and estimation of failure parameters

Machine-learning methods and apparatus are disclosed to determine frictional state or other parameters in an earthquake zone or other failing medium, using acoustic emission, seismic waves, or other detectable indicators of microscopic processes. Predictions of future failures are demonstrated in different regimes. A classifier is trained using time series of acoustic emission data along with historic data of frictional state or failure events. In disclosed examples, random forests and gradient boost trees are used, and grid-search or EGO procedures are used for hyperparameter tuning. Once trained, the classifier can be applied to testing or live data in order to assess a frictional state, assess seismic hazard, or make predictions regarding a future failure event. The technology has been developed in a double direct shear apparatus, but can be widely applied to seismic faults, other terrestrial failures, or failures in man-made structures. Variations are disclosed.

Johnson, Paul Allan↗

DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗

DOE OSTI2026

SPAROW: Stochastic Programming and Related Optimization Workflows

SAND2026-16703O SPAROW: Stochastic Programming and Related Optimization Workflows is a Python library tool that facilitates the development and solution of stochastic programming problems. It provides a user-friendly class structure for defining stochastic programs through scenario-based representations of uncertainties. SPAROW incorporates multiple optimization strategies, including integer programming with all scenarios, progressive hedging, Benders decomposition, and Snoglode, a novel technique developed by Carnegie Mellon University. It also features interfaces to external solvers and functions that are commonly used in analysis workflows, making it applicable to a wide range of scientific and engineering design challenges, particularly in power grid planning. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Hart, William [Sandia National Lab. (SNL-NM), Albu↗