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.

SourceResearch connectionFrom yearThrough year

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

At least 1,783 records · Page 99

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Advanced Photovoltaic Inverter Control Development and Validation in a Controller-Hardware-in-the-Loop Test Bed

Penetration of solar photovoltaic (PV) installations into the grid has increased in recent years. In the past, photovoltaic installations have not included grid-support functionality. Grid support functions including Volt-VAR (VVAR), Frequency-Watt (FW), Volt-Watt (VW), frequency/voltage-ride through, as well as other inverter functions such as anti-islanding are recommended by standards such as IEEE 1547 to allow PV inverters to contribute to acceptable operation and stability of distribution systems with high penetrations of PV and to ensure predictable response of the PV inverters. First, low-level and high-level grid support inverter control functions were developed and tested in simulation. The developed controller is able to integrate advanced controller functions in a coordinated approach, avoiding any conflict between the different control objectives. The algorithms were then programmed on an off-the-shelf embedded controller that included a dual-core CPU and FPGA. The details of the high-level grid support function and the lower level switching control will be presented. Second, this programmed controller is tested using a controller-hardware-in-the-loop (CHIL) test bed setup using an FPGA-based real-time simulator. The CHIL setup has the capability to run at a 500 nano-second time-step to accommodate the 20 kHz switching frequency of the developed controller. The advanced control function and CHIL test bed details provided here will aide future researchers to design, implement, and test advanced inverter functions for PV inverters.

14 SOLAR ENERGY↗

DOE OSTI2022

Distributed Energy Resource Cybersecurity Framework and Cyber Range Integration

Distributed energy resource (DER) systems feature complex, data-driven communications networks that require careful system coordination and constant vigilance to ensure that grid assets are secure. Because DERs are an important component of the decarbonization strategy, agencies need to secure energy data that could implicate issues of national security if compromised. To help federal energy managers assess, monitor, and manage cybersecurity while achieving decarbonization, the National Renewable Energy Laboratory's (NREL's) Distributed Energy Resource Cybersecurity Framework (DER-CF) offers a comprehensive, web-based assessment tool focusing on cyber governance or policies, technical management, and physical security. The DER-CF currently presents users with a series of pertinent cybersecurity questions that are used to generate a site-specific report and recommendations. This paper outlines a plan to integrate the DER-CF with another key asset-NREL's cyber range-to visualize cybersecurity resilience and compliance and to enhance the usability and accessibility of the DER-CF for federal facility energy managers and planners. This integration will result in a visualization environment to interpret and interact with compliance data. Its development will include regular conversations with stakeholders to assess the effectiveness of these efforts, refine the visualization capability, and ensure its value to our partners.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2015

OpenACC acceleration of an unstructured CFD solver based on a reconstructed discontinuous Galerkin method for compressible flows

Here, an OpenACC directive-based graphics processing unit (GPU) parallel scheme is presented for solving the compressible Navier–Stokes equations on 3D hybrid unstructured grids with a third-order reconstructed discontinuous Galerkin method. The developed scheme requires the minimum code intrusion and algorithm alteration for upgrading a legacy solver with the GPU computing capability at very little extra effort in programming, which leads to a unified and portable code development strategy. A face coloring algorithm is adopted to eliminate the memory contention because of the threading of internal and boundary face integrals. A number of flow problems are presented to verify the implementation of the developed scheme. Timing measurements were obtained by running the resulting GPU code on one Nvidia Tesla K20c GPU card (Nvidia Corporation, Santa Clara, CA, USA) and compared with those obtained by running the equivalent Message Passing Interface (MPI) parallel CPU code on a compute node (consisting of two AMD Opteron 6128 eight-core CPUs (Advanced Micro Devices, Inc., Sunnyvale, CA, USA)). Speedup factors of up to 24× and 1.6× for the GPU code were achieved with respect to one and 16 CPU cores, respectively. The numerical results indicate that this OpenACC-based parallel scheme is an effective and extensible approach to port unstructured high-order CFD solvers to GPU computing.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

DOE OSTI2021

Improving Prediction of Surface Solar Irradiance Variability by Integrating Observed Cloud Characteristics and Machine Learning

A 5-year, 1-minute resolution observational dataset of clouds and solar radiation was produced that includes two metrics of the variability in surface solar irradiance due to cloud type and fractional sky cover. Multiple regression models were trained to fit observations of surface solar irradiance variability from those two cloud property predictors. We found that ensemble tree-based methods, Random Forest and Gradient Boosting Machine, have the least overfitting issues and showed the best performance with an R2 of 0.42. While the observational data trained in this study was only from one site, the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in Oklahoma, initial comparisons of the seasonality of the statistics suggest that these results are relatively weather regime independent; the generality of such a finding across sites will be tested in future work. The observational data and developed machine learning model are being used to create a numerical weather prediction model parameterization to enable day-ahead solar variability prediction in a computationally efficient way. This is a first step towards creating a new paradigm of predicting day-ahead variability with the potential to provide a new tool to improve grid operation, planning, and resilience.

Riihimaki, Laura↗

DOE OSTI2025

Representing Fine‐Scale Topographic Effects on Surface Radiation Balance in Hyper‐Resolution Land Surface Models

Land surface models are increasingly used to simulate land surface processes at hyper-spatial resolutions (e.g., ∼1 km). As model resolution increases, grid-scale topographic effects on surface radiation fluxes and their interactions between adjacent grids become more pronounced. However, current land surface models routinely neglect the fine-scale topographic effects on surface radiation balance. This study developed physically-based and computationally-efficient parameterizations (fineTOP) that explicitly resolve fine-scale topographic effects on downward shortwave and longwave radiation as well as land surface radiative properties. The newly developed parameterizations were implemented and tested in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Multi-decadal km-resolution ELM simulations over the California Sierra Nevada show that fine-scale topography significantly impacts the surface energy balance and snow processes across seasons. Slope determines the magnitude of topographic effects, while aspect controls their sign. For slopes larger than 30°, topography-induced change in annual surface temperature can be as large as 3.3 K. Regionally, the mean value and standard deviation of topography-induced changes in annual surface temperature are −0.22 ± 0.38 K and +0.25 ± 0.37 K over north-facing and south-facing slopes, respectively. Topography-induced changes in surface radiative properties account for 3.5% ± 13.8% of total topographic effects on annual net radiation. With fineTOP, ELM captures the aspect-dependence of snow cover fraction, snow water equivalent, and land surface temperature found in MODIS satellite observations and a snow reanalysis data set, while the default ELM fails to capture this phenomenon. The enhanced capability to represent fine-scale topographic effects on surface radiation balance can be used to advance understanding of the role of fine-scale topography in land surface processes and land-atmosphere interactions over mountainous regions.

Hao, Dalei [Pacific Northwest National Laboratory ↗

DOE OSTI2020

Analysis and visualization of energy densities. II. Insights from linear-response time-dependent density functional theory calculations

Inspired by the analysis of Kohn–Sham energy densities by Nakai and coworkers, we extended the energy density analysis to linear-response time-dependent density functional theory (LR-TDDFT) calculations. Using ethylene–tetrafluoroethylene and oxyluciferin–water complexes as examples, distinctive distribution patterns were demonstrated for the excitation energy densities of local excitations (within a molecular fragment) and charge-transfer excitations (between molecular fragments). It also provided a simple way to compute the effective energy of both hot carriers (particle and hole) from charge-transfer excitations via an integration of the excitation energy density over the donor and acceptor grid points.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2014

Aspherical nucleosynthesis in a core-collapse supernova with 25 M {sub ☉} standard progenitor

The problem of nucleosynthesis was studied within an aspherical supernova model. The explosive burning was computed in a star of 25 M {sub ☉} initial mass on its final stage of evolution. The chemical composition of a presupernova was taken from realistic evolutionary computations. A piecewise parabolic method on a local stencil was applied to simulate the hydrodynamics of the explosion. The gravity was recomputed by a Poisson solver on a fine grid as the explosion developed. A detailed yield of chemical elements was performed as a post-processing step using the tracer particles method. The produced nuclei formed a layer-like structure enclosing large fragments of nickel and iron-group isotopes that were pushed away from the central region by an explosion along the polar direction. The light nuclei were preferentially moving along the equatorial plane forming a torus-like structure.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

IGMS: An Integrated ISO-to-Appliance Scale Grid Modeling System

This paper describes the Integrated Grid Modeling System (IGMS), a novel electric power system modeling platform for integrated transmission-distribution analysis that co-simulates off-the-shelf tools on high performance computing (HPC) platforms to offer unprecedented resolution from ISO markets down to appliances and other end uses. Specifically, the system simultaneously models hundreds or thousands of distribution systems in co-simulation with detailed Independent System Operator (ISO) markets and AGC-level reserve deployment. IGMS uses a new MPI-based hierarchical co-simulation framework to connect existing sub-domain models. Our initial efforts integrate opensource tools for wholesale markets (FESTIV), bulk AC power flow (MATPOWER), and full-featured distribution systems including physics-based end-use and distributed generation models (many instances of GridLAB-D[TM]). The modular IGMS framework enables tool substitution and additions for multi-domain analyses. This paper describes the IGMS tool, characterizes its performance, and demonstrates the impacts of the coupled simulations for analyzing high-penetration solar PV and price responsive load scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2017

Algebraic Multigrid Benchmark

AMG is a parallel algebraic multigrid solver for linear systems arising from problems on unstructured grids. It has been derived directly from the BoomerAMG solver in the hypre library, a large linear solvers library that is being developed in the Center for Applied Scientific Computing (CASC) at LLNL and is very similar to the AMG2013 benchmark with additional optimizations. The driver provided in the benchmark can build various test problems. The default problem is a Laplace type problem with a 27-point stencil, which can be scaled up and is designed to solve a very large problem. A second problem simulates a time dependent problem, in which successively various smnllcr systems are solved.

Yang, Ulrike↗

DOE OSTI2013

Verification & Validation of High-Order Short-Characteristics-Based Deterministic Transport Methodology on Unstructured Grids

The research team has developed a practical, high-order, discrete-ordinates, short characteristics neutron transport code for three-dimensional configurations represented on unstructured tetrahedral grids that can be used for realistic reactor physics applications at both the assembly and core levels. This project will perform a comprehensive verification and validation of this new computational tool against both a continuous-energy Monte Carlo simulation (e.g. MCNP) and experimentally measured data, an essential prerequisite for its deployment in reactor core modeling. Verification is divided into three phases. The team will first conduct spatial mesh and expansion order refinement studies to monitor convergence of the numerical solution to reference solutions. This is quantified by convergence rates that are based on integral error norms computed from the cell-by-cell difference between the code’s numerical solution and its reference counterpart. The latter is either analytic or very fine- mesh numerical solutions from independent computational tools. For the second phase, the team will create a suite of code-independent benchmark configurations to enable testing the theoretical order of accuracy of any particular discretization of the discrete ordinates approximation of the transport equation. For each tested case (i.e. mesh and spatial approximation order), researchers will execute the code and compare the resulting numerical solution to the exact solution on a per cell basis to determine the distribution of the numerical error. The final activity comprises a comparison to continuous-energy Monte Carlo solutions for zero-power critical configuration measurements at Idaho National Laboratory’s Advanced Test Reactor (ATR). Results of this comparison will allow the investigators to distinguish between modeling errors and the above-listed discretization errors introduced by the deterministic method, and to separate the sources of uncertainty.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2017

Techbelt Energy Innovation Center

This project consisted of three main components 1) The primary goal of the project was to renovate and upgrade an existing commercial building to the highest possible environmentally sustainable level for the purpose of creating an energy incubator. This initiative was part of the Infrastructure Technologies Program, through which a sustainable energy demonstration facility was to be created and used as a research and community outreach base for sustainable energy product and process incubation; 2) In addition, fundamental energy related research on wind energy was performed; a shrouded wind turbine on the Youngstown State University campus was commissioned; and educational initiatives were implemented; and 3) The project also included an education and outreach component to inform and educate the public in sustainable energy production and career opportunities. Youngstown State University and the Tech Belt Energy Innovation Center (TBEIC) renovated a 37,000 square foot urban building which is now being used as a research and development hub for the region’s energy technology innovation industry. The building houses basic research facilities and business development in an incubator format. In addition, the TBEIC performs community outreach and education initiatives in advanced and sustainable energy. The building is linked to a back warehouse which will eventually be used as a build-out for energy laboratory facilities. The projects research component investigated shrouded wind turbines, and specifically the “Windcube” which was renamed the “Wind Sphere” during the course of the project. There was a specific focus on the development in the theory of shrouded wind turbines. The goal of this work was to increase the potential efficiency of wind turbines by improving the lift and drag characteristics. The work included computational modeling, scale models and full-sized design and construction of a test turbine. The full-sized turbine was built on the YSU campus as a grid-tie system that supplies the YSU research facility. Electrical power meters and weather monitors were installed to record the power generated and aid in continued study. In addition, an education/outreach component to help elicit creative engineering and design from amongst area students, faculty, entrepreneurs, and small business in the energy related fields was performed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2018

Evaluating the Feasibility of MCNP/FDTD EMP Calculations at Altitudes above 20km

This report investigates the use of MCNP/FDTD for 2D EMP calculations at altitudes above 20km. The objective is to determine the feasibility of performing EMP calculations for high altitude nuclear detonations using the MCNP/FDTD codes. The work in this report indicates that the issue with running such EMP calculations lies more in MCNP than in FDTD. In general, when performing EMP calculations with MCNP/FDTD, the majority of the computation time is in getting the results from MCNP. Typically, for EMP calculations with spatial domains of about 5 to 10km in both the x and z directions, the MCNP calculations take weeks to complete and the FDTD portion of the computation only takes one to two days. The problems with running MCNP on large problems are largely due to memory issues. One issue is that MCNP requires all tally grids of a simulation to reside in the memory of a single node (128GB on snow) and this, in turn, requires multiple MCNP simulations to complete large spatial/time domain problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Challenging the Curse of Dimensionality in Multidimensional Numerical Integration by Using a Low-Rank Tensor-Train Format

Numerical integration is a basic step in the implementation of more complex numerical algorithms suitable, for example, to solve ordinary and partial differential equations. The straightforward extension of a one-dimensional integration rule to a multidimensional grid by the tensor product of the spatial directions is deemed to be practically infeasible beyond a relatively small number of dimensions, e.g., three or four. In fact, the computational burden in terms of storage and floating point operations scales exponentially with the number of dimensions. This phenomenon is known as the curse of dimensionality and motivated the development of alternative methods such as the Monte Carlo method. The tensor product approach can be very effective for high-dimensional numerical integration if we can resort to an accurate low-rank tensor-train representation of the integrand function. In this work, we discuss this approach and present numerical evidence showing that it is very competitive with the Monte Carlo method in terms of accuracy and computational costs up to several hundredths of dimensions if the integrand function is regular enough and a sufficiently accurate low-rank approximation is available.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

DMTN-221: Periodicity Analysis in Alert Production

The baselined timeseries features to be computed in Alert Production include a Lomb-Scargle periodogram ran on two classes of variable systems: RR Lyrae and Eclipsing Binaries. Based on a simulated LSST-like cadence light curves taken from the Extended LSST Astronomical Time-series Classification Challenge (ELAsTiCC) we perform an end-to-end test to characterize the periodicity recovery on the Alert Production multi-band light curves. In both variable classes, we found that a single-band Lomb-Scargle implementation yields to a low fraction of recovered periods, with a significant preference on the simple periodic phenomena such as RR Lyrae. We also investigated the results from a multi-band Lomb-Scargle and found an increased fraction of recovered periodicities above 15% for the eclipsing binaries, and over 80$\%$ for the RR Lyrae stars. Our findings suggest that a multi-band Lomb-Scargle should be implemented for searching periodic phenomena through AP. We also asses the computational and scientific performance of several configurations on simulated alert data and find that our current configuration scales linearly with the number of detections while assuming an heuristic frequency grid.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Algebraic Multigrid Benchmark

Yang, Ulrike↗

DOE OSTI2018

A Heuristic Approach to the Post-disturbance and Stochastic Pre-disturbance Microgrid Formation Problem

Microgrid formation is a potential solution in post-disaster electric grid recovery efforts. Recent works propose distribution level microgrid formation models using mixed-integer linear programming techniques. However, these models can only be solved for small and medium size power systems due to their computational intractability. In this article we introduce a heuristic approach that allows to approximately solve the post-disturbance microgrid formation problem for medium to large, more realistic, instances. Furthermore, the proposed approach allows to approximately solve the pre-disturbance microgrid formation problem (a stochastic version of the post-disturbance problem), in which the aim is to allocate extra generation capacity to the network to immunize it, as best as possible, against an uncertain cascading failure. Our results are illustrated by solving versions of the problem on ten test cases with up to 3012 buses.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

An Adaptive Virtual Oscillator Control Structure for Grid-Forming Inverters

The electrical grid is facing unprecedented challenges due to the increasing penetration of inverter-based resources. Grid forming inverters (GFMIs) are a promising technology to address these challenges. Recently, the virtual oscillator based GFMI control is attracting more attention due to its superior performance over other control strategies. In this article, an adaptive control strategy is proposed to provide flexible operation and transition between grid-connected and islanded modes. The controller adapts the virtual oscillator's parameter values depending on the operation mode. It also provides a feedback signal using a measured frequency to account for any steady-state errors and to allow a seamless transition from grid-connected to islanded mode. Finally, to show the feasibility of the proposed controller, this article discusses the simulation results from the implementation of the controller on a single inverter system and on a group of inverters on a large practical utility feeder, the IEEE 13 node feeder, using the DIgSILENT simulation environment.

adaptive control↗