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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Geocellular model of St. Peter Sandstone for University of Illinois at Urbana-Champaign DDU Feasibility Study

The geocellular model of the St. Peter Sandstone was constructed for the University of Illinois at Urbana-Champaign DDU feasibility study. Starting with the initial area of review (18.0 km by 18.1 km [11.2 miles by 11.3 miles]) the boundaries of the model were trimmed down to 9.7 km by 9.7 km (6 miles by 6 miles) to ensure that the model enclosed a large enough volume so that the cones of depression of both the production and injection wells would not interact with each other, while at the same time minimizing the number of cells to model to reduce computational time. The grid-cell size was set to 61.0 m by 61.0 m (200 feet by 200 feet) for 160 nodes in the X and Y directions. The top surface of the St. Peter Sandstone was provided by geologists working on the project, and the average thickness of the formation was taken from the geologic prospectus they provided. An average thickness of 68.6 m (225 feet) was used for the St. Peter Sandstone, resulting in 45 layers for the model. Petrophysical data was taken from available rotary sidewall core data (Morrow et al., 2017). As geothermal properties (thermal conductivity, specific heat capacity) are closely related to mineralogy, specifically the percentage of quartz, available mineralogical data was assembled and used with published data of geothermal values to determine these properties (Waples and Waples, 2004; Robertson, 1988). The St. Peter Sandstone was divided into facies according to similar geothermal and petrophysical properties, and distributed according to available geophysical log data and prevailing interpretations of the depositional/diagenetic history (Will et al. 2014). Petrophysical and geothermal properties were distributed through geostatistical means according to the associated distributions for each lithofacies. The formation temperature was calculated, based on data from continuous temperature geophysical log from a deep well drilled into the Precambrian basement at the nearby Illinois Basin Decatur Project (IBDP) where CO2 is currently being sequestered (Schlumberger, 2012). Salinity values used in the model were taken from regional studies of brine chemistry in the St. Peter Sandstone, including for the IBDP (e.g., Panno et al. 2018). After being reviewed by the project's geologists, the model was then passed onto the geological engineers to begin simulations of the geothermal reservoir and wellbores.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Geocellular Model of Mt. Simon Sandstone for University of Illinois at Urbana-Champaign DDU feasibility study

The geocellular model of the Mt. Simon Sandstone was constructed for the University of Illinois at Urbana-Champaign DDU feasibility study. Starting with the initial area of review (18.0 km by 18.1 km [11.2 miles by 11.3 miles]) the boundaries of the model were trimmed down to 9.7 km by 9.7 km (6 miles by 6 miles) to ensure that the model enclosed a large enough volume so that the cones of depression of both the production and injection wells would not interact with each other, while at the same time minimizing the number of cells to model to reduce computational time. The grid-cell size was set to 61.0 m by 61.0 m (200 feet by 200 feet) for 160 nodes in the X and Y directions. Within the model, 67 layers are represented that are parameterized with their sediment/rock properties and petrophysical data. The top surface of the Mt. Simon Sandstone was provided by geologists working on the project, and the average thickness of the formation was taken from the geologic prospectus they provided. An average thickness of 762 m (2500 feet) was used for the Mt. Simon Sandstone, resulting in 60 layers for the model. Petrophysical data was taken from available rotary sidewall core data (Morrow et al., 2017). As geothermal properties (thermal conductivity, specific heat capacity) are closely related to mineralogy, specifically the percentage of quartz, available mineralogical data was assembled and used with published data of geothermal values to determine these properties (Waples and Waples, 2004; Robertson, 1988). The Mt. Simon Sandstone was divided into three separate units (lower, middle, upper) according to similar geothermal and petrophysical properties, and distributed according to available geophysical log data and prevailing interpretations of the depositional/diagenetic history (Freiburg et al. 2016). Petrophysical and geothermal properties were distributed through geostatistical means according to the associated distributions for each lithofacies. The formation temperature was calculated, based on data from continuous temperature geophysical log from a deep well drilled into the Precambrian basement at the nearby Illinois Basin Decatur Project (IBDP) where CO2 is currently being sequestered (Schlumberger, 2012). Salinity values used in the model were taken from regional studies of brine chemistry in the Mt. Simon Sandstone, including for the IBDP (e.g., Panno et al. 2018). After being reviewed by the project's geologists, the model was then passed onto the geological engineers to begin simulations of the geothermal reservoir and wellbores.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Combinatorial Algorithms in Scientific Computing

We provide the final report for this grant, detailing the publications, software produced, students trained who have joined the DOE workforce, and the impact our work has had on computational mathematics and related disciplines.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multi-Stage and Multi-Timescale Robust Co-Optimization Planning for Reliable and Sustainable Power Systems. Final Report

In this project, Clarkson University, in collaboration with Southern Methodist University and University of Pittsburgh, conducted the research to model, design, and implement a sophisticated generation and transmission co-optimization planning decision tool, called Multi-stage and Multi-timescale robust Co-Optimization Planning (MMCOP). The MMCOP decision tool intends to facilitate generation and transmission co-optimization planning of emerging power systems, while mitigating risks and uncertainties in both short-term operation dynamics and long-term policy and technology changes. Long-term power system planning aims at optimizing asset utilization by investing in a proper mix of various generation technologies and transmission lines to supply the future load growth. In particular, the Clean Power Plan (CPP), which is designed to combat climate change and reduce carbon emissions by setting a national limit on carbon pollution from power plants, may dramatically change the landscape of the power industry by further promoting clean energy and phasing out emissions-intensive generation technologies. In addition, novel non-wire alternatives (e.g., demand response (DR), distributed generation (DG), energy efficiency (EE), and smart grid technologies) and the computational complexity for large-scale systems significantly complicate the system planning procedure even further. However, existing conventional planning approaches neglect short-term variability and uncertainty of renewable energy, hourly chronological operation details, and physical nonlinear characteristics of the alternating current transmission network. As a result, existing conventional planning approaches may not work properly, and power systems reliability could be in jeopardy. In observing the limitations of existing conventional planning approaches and addressing new challenges of emerging power systems, the main scope of this project to develop co-optimization planning models within a multi-stage and multi-timescale framework. In particular, random contingencies, key uncertainty factors, and AC power flows are included to derive expansion plans while considering both long-term reliability and short-term flexibility.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

BACCHUS Analysis of Weak Lines in APOGEE Spectra (BAWLAS)

Elements with weak and blended spectral features in stellar spectra are challenging to measure and require specialized analysis methods to precisely measure their chemical abundances. In this work, we have created a catalog of approximately 120,000 giants with high signal-to-noise Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 17 (DR17) spectra, for which we explore weak and blended species to measure Na, P, S, V, Cu, Ce, and Nd abundances and 12 C/ 13 C isotopic ratios. We employ an updated version of the Brussels Automatic Code for Characterizing High-accuracy Spectra (BACCHUS) code to derive these abundances using the stellar parameters measured by APOGEE's DR17 Stellar Parameters and Chemical Abundances Pipeline, quality flagging to identify suspect spectral lines, and a prescription for upper limits. Combined, these allow us to provide our BACCHUS Analysis of Weak Lines in APOGEE Spectra catalog of precise chemical abundances for these weak and blended species, which agrees well with the literature and improves upon APOGEE abundances for these elements, some of which are unable to be measured with APOGEE's current, grid-based approach without computationally expensive expansions. This new catalog can be used alongside APOGEE and provides measurements for many scientific applications ranging from nuclear physics to Galactic chemical evolution and Milky Way population studies. To illustrate this we show some examples of uses for this catalog, such as showing that we observe stars with enhanced s-process abundances or that we can use the 12 C/ 13 C ratios to explore extra mixing along the red giant branch.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Systematic Framework for Tuning Open-Source Multifunctional IBR Models To Emulate OEM Black-Box Fault Dynamics

This paper presents a systematic framework to tune a generic IBR EMT model to match with an OEM provided balckbox inverter model based on the fault current responses. The key learnings and findings are summarized as follows: The tunable key parameters include inner control loops and current limiters to align the fault current magnitude, sequence content, and phase trajectories with the OEM models across diverse fault type and locations. The tuned model's fidelity is validated through comparative analysis with an OEM blackbox model, assessing both the fault current response and the responses of multiple relay elements. The results demonstrate the tuned generic model can trigger relay decision logic that is identical or near identical to that of the OEM model, thus generating very good match model for fault studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

A source-attractor approach to network detection of radiation sources

Radiation source detection using a network of detectors is an active field of research for homeland security and defense applications. We propose Source-attractor Radiation Detection (SRD) method to aggregate measurements from a network of detectors for radiation source detection. SRD method models a potential radiation source as a magnet -like attractor that pulls in pre-computed virtual points from the detector locations. A detection decision is made if a sufficient level of attraction, quantified by the increase in the clustering of the shifted virtual points, is observed. Compared with traditional methods, SRD has the following advantages: i) it does not require an accurate estimate of the source location from limited and noise-corrupted sensor readings, unlike the localizationbased methods, and ii) its virtual point shifting and clustering calculation involve simple arithmetic operations based on the number of detectors, avoiding the high computational complexity of grid-based likelihood estimation methods. We evaluate its detection performance using canonical datasets from Domestic Nuclear Detection Office s (DNDO) Intelligence Radiation Sensors Systems (IRSS) tests. SRD achieves both lower false alarm rate and false negative rate compared to three existing algorithms for network source detection.

Wu, Qishi↗
DOE OSTI2017

TABULATED EQUIVALENT SDR FLAMELET (TESF) MODEFL

The code consists of an implementation of a novel tabulated combustion model for non-premixed flames in CFD solvers. This novel technique/model is used to implement an unsteady flamelet tabulation without using progress variables for non-premixed flames. It also has the capability to include history effects which is unique within tabulated flamelet models. The flamelet table generation code can be run in parallel to generate tables with large chemistry mechanisms in relatively short wall clock times. The combustion model/code reads these tables. This framework can be coupled with any CFD solver with RANS as well as LES turbulence models. This framework enables CFD solvers to run large chemistry mechanisms with large number of grids at relatively lower computational costs. Currently it has been coupled with the Converge CFD code and validated against available experimental data. This model can be used to simulate non-premixed combustion in a variety of applications like reciprocating engines, gas turbines and industrial burners operating over a wide range of fuels.

KUNDU, PRITHWISH↗
DOE OSTI2017

Transmission Line Ampacity Improvements of AltaLink Wind Plant Overhead Tie-Lines Using Weather-Based Dynamic Line Rating

Abstract—Overhead transmission lines (TLs) are conventionally given seasonal ratings based on conservative environmental assumptions. Such an approach often results in underutilization of the line ampacity as the worst conditions prevail only for a short period over a year/season. We presents dynamic line rating (DLR) as an enabling smart grid technology that adaptively computes ratings of TLs based on local weather conditions to utilize additional headroom of existing lines. In particular, general line ampacity state solver utilizes measured weather data for computing the real-time thermal rating of the TLs. The performance of the presented method is demonstrated from a field study of DLR technology implementation on four TL segments at AltaLink, Canada. The performance is evaluated and quantified by comparing the existing static and proposed dynamic line ratings, and the potential benefits of DLR for enhanced transmission assets utilization. For the given line segments, the proposed DLR results in real-time ratings above the seasonal static ratings for most of the time; up to 95.1% of the time, with a mean increase of 72% over static rating.

17 WIND ENERGY↗
DOE OSTI2020

Conditional Karhunen-Loève expansion for uncertainty quantification and active learning in partial differential equation models

We use a conditional Karhunen-Lo` eve (KL) model to quantify and reduce uncertainty in a stochastic partial differential equation (SPDE) problem with partially-known space-dependent coefficient, Y (x). We assume that a small number of Y (x) measurements are available and model Y (x) with a KL expansion. We achieve reduction in uncertainty by conditioning the KL expansion coefficients on measurements. We consider two approaches for conditioning the KL expansion: In Approach 1, we condition the KL model first and then truncate it. In Approach 2, we first truncate the KL expansion and then condition it. We employ the conditional KL expansion together with Monte Carlo and sparse grid collocation methods to compute the moments of the solution of the SPDE problem. Uncertainty of the problem is further reduced by adaptively selecting additional observation locations using two active learning methods. Method 1 minimizes the variance of the PDE coefficient, while Method 2 minimizes the variance of the solution of the PDE. We demonstrate that conditioning leads to dimension reduction of the KL representation of Y (x). For a linear diffusion SPDE with uncertain log-normal coefficient, we show that Approach 1 provides a more accurate approximation of the conditional log-normal coefficient and solution of the SPDE than Approach 2 for the same number of random dimensions in a conditional KL expansion. Furthermore, Approach 2 provides a good estimate for the number of terms of the truncated KL expansion of the conditional field of Approach 1. Finally, we demonstrate that active learning based on Method 2 is more efficient for uncertainty reduction in the SPDE’s states (i.e., it leads to a larger reduction of the variance) than active learning using Method 2.

Conditioned Karhunen-Lo` eve expanion, machine lea↗
DOE OSTI2017

TABULATED EQUIVALENT SDR FLAMELET (TESF)

The code consists of an implementation of a novel tabulated combustion model for non-premixed flames in CFD solvers. This novel technique/model is used to implement an unsteady flamelet tabulation without using progress variables for non-premixed flames. It also has the capability to include history effects which is unique within tabulated flamelet models. The flamelet table generation code can be run in parallel to generate tables with large chemistry mechanisms in relatively short wall clock times. The combustion model/code reads these tables. This framework can be coupled with any CFD solver with RANS as well as LES turbulence models. This framework enables CFD solvers to run large chemistry mechanisms with large number of grids at relatively lower computational costs. Currently it has been coupled with the Converge CFD code and validated against available experimental data. This model can be used to simulate non-premixed combustion in a variety of applications like reciprocating engines, gas turbines and industrial burners operating over a wide range of fuels.

Kundu, Prithwish↗
DOE OSTI2020

Validation of blade-resolved computational fluid dynamics for a MW-scale turbine rotor in atmospheric flow

In this study, simulation results of two different computational fluid dynamics codes, Nalu-Wind and EllipSys3D, are presented for a wind turbine rotor in complex yawed and sheared inflow. The results are compared to measurements from the DanAero experiments, to validate computed pressures and azimuthal trends. Despite different code methodologies and grid setups, the codes agree well in computed pressures and integrated forces along the blade for all blade azimuthal positions, however with some discrepancy in the very yawed case. Additionally, both codes capture well the azimuthal trends and force levels seen in measurements. Investigation into discrepancies shows that expanding grids before the rotor, lead to smearing of the wind profiles, which is likely the main cause of the differences in the results between the codes. Additionally, omission of the ground constraint cause discrepancies in relative velocity seen by the passing blade, due to an over speeding beneath the rotor.

17 WIND ENERGY↗
DOE OSTI2023

A Transductive Graph Neural Network learning for Grid Resilience Analysis

Power grids are critical infrastructures that require robust resilience analysis to ensure reliable and uninterrupted electricity supply. Traditional simulation-based methods for grid resilience analysis suffer from computational complexity and limited ability to capture the full spectrum of potential disruptions. This paper presents a novel approach to enhance grid resilience by leveraging transductive graph neural network (GNN) learning to identify critical nodes and links. By leveraging the graph structure and system features, GNNs effectively learn resilience metrics and accurately identify critical nodes based on actual grid operational behavior. The efficacy of the proposed approach is demonstrated through case studies on node criticality scoring and critical node/line identification in cascading outage scenarios. The results highlight the advantages of learning-based methods over traditional simulation-based approaches and their potential to revolutionize grid resilience analysis. The contributions of this paper include a graph-based scalable approach for fast cascading analysis, an inductive formulation for training GNN models, and a transfer learning-based approach to scale the model to largescale power systems.

grid resilience, graph neural networks, transducti↗
DOE OSTI2017

Standard grids for high-precision integration of modern density functionals: SG-2 and SG-3

Density-functional approximations developed in the past decade necessitate the use of quadrature grids that are far more dense than those required to integrate older generations of functionals. This category of difficult-to-integrate functionals includes meta-generalized gradient approximations, which depend on orbital gradients and/or the Laplacian of the density, as well as functionals based on B97 and the popular “Minnesota” class of functionals, each of which contain complicated and/or oscillatory expressions for the exchange inhomogeneity factor. Following a strategy introduced previously by Gill and co-workers to develop the relatively sparse “SG-0” and “SG-1” standard quadrature grids, we introduce two higher-quality grids that we designate SG-2 and SG-3, obtained by systematically “pruning” medium- and high-quality atom-centered grids. The pruning procedure affords computational speedups approaching a factor of two for hybrid functionals applied to systems of ~100 atoms, without significant loss of accuracy. The grid dependence of several popular density functionals is characterized for various properties.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Final Technical Report: Integrated Distribution-Transmission Analysis for Very High Penetration Solar PV

Transmission and distribution simulations have historically been conducted separately, echoing their division in grid operations and planning while avoiding inherent computational challenges. Today, however, rapid growth in distributed energy resources (DERs)--including distributed generation from solar photovoltaics (DGPV)--requires understanding the unprecedented interactions between distribution and transmission. To capture these interactions, especially for high-penetration DGPV scenarios, this research project developed a first-of-its-kind, high performance computer (HPC) based, integrated transmission-distribution tool, the Integrated Grid Modeling System (IGMS). The tool was then used in initial explorations of system-wide operational interactions of high-penetration DGPV.

14 SOLAR ENERGY↗
DOE OSTI2023

Machine learning from RANS and LES to inform coarse grid simulations

Nuclear system thermal hydraulic analysis has historically relied on computationally inexpensive 1D codes. However, such tools are unable to capture multiscale multidimensional effects in large nuclear reactor enclosures. On the other hand, simulations with higher fidelity can be too expensive for such purposes. One of the ways to reduce computational cost is to perform simulations on a coarse grid, which, unfortunately, introduces large discretization errors. In this paper, two high-to-low data-driven approaches are investigated: (1) a coarse grid turbulence model to predict eddy viscosity and (2) correction of errors in coarse grid velocity fields. The approaches aim to reduce grid- and turbulence model-induced errors in coarse grid Reynolds-averaged Navier–Stokes (RANS) simulations. Two sources of high-fidelity data, RANS and large eddy simulations (LES), are explored. To extract the eddy viscosity from the LES data, an inverse optimization problem is solved. However, the LES eddy viscosity is shown to be comparable to the RANS eddy viscosity in terms of error reduction. Therefore, the directly available RANS eddy viscosity was used to develop a coarse grid data-driven turbulence model. Additionally, error correction in velocity is used to reduce the remaining uncertainties and bring the results closer to reality. In conclusion, the performance of the frameworks is demonstrated for a scaled upper plenum of a gas-cooled reactor facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗