Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid size”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A correction procedure for self-induced velocity of a finite-sized particle in two-way coupled Euler–Lagrange simulations

The importance of incorporating a correction to undo the self-induced perturbation velocity of a particle, when its size becomes comparable to the Eulerian grid, in a two-way coupled Euler–Lagrange (EL) simulation is now well appreciated. The present work improves upon the prior correction procedures in a few important ways. First, the past correction procedures have been scalar-based with the assumption that the quasi-steady force is the source of self-induced velocity perturbation. Here we generalize to a vector correction procedure and thereby the directions of feedback force and relative velocity can be different. This allows the correction procedure to be used even in the presence of added-mass, history, and lift forces. Second, the effect of a nearby wall has been systematically included in the correction procedure. The correction procedure depends on fundamental Oseen solutions of streamwise and transverse regularized feedback forces. We present a Fourier transform-based analytical approach to obtaining these regularized Oseen solutions. We also present a step-by-step numerical procedure for obtaining the Oseen solutions in any EL code. With the analytical or numerical Oseen functions, the correction procedure can be easily implemented in any EL code. Iterations are required in solving the implicit correction equations and it is demonstrated that the correction procedure converges rapidly within three or four iterations. In conclusion, a simple empirical approach is also presented to account for unsteady effects in the correction procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Reducing the computational burden of a microgrid energy management system

As renewable technology advances and decreases in cost, microgrids are becoming an appealing means of distributed generation both for isolated communities and integrated with existing electrical grid systems. Due to their small size, however, microgrids may have financial limitations which preclude them from using commercial software to optimize control of their assets. Open-source optimization solvers are a viable alternative, but increase computation time. This work expands on a rolling horizon optimization framework for economic dispatch within an existing residential microgrid located in Hoover, Alabama. The microgrid has an open-source solver requirement and a need for quick solution time on a rolling horizon as opposed to a day-ahead commitment. We present a method of reducing integer variables by relaxation which completes two goals: reduction in computation time for real-time operations, and reduction in daily operational cost for the microgrid. Seasonal data for load and photovoltaic (PV) power was also collected from the microgrid to facilitate simulation testing. Computation time was successfully reduced using multiple variations of the relaxation method, while obtaining solution quality with operational cost similar to or better than the original model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evolution of the SLATE linear algebra library

SLATE (Software for Linear Algebra Targeting Exascale) is a distributed, dense linear algebra library targeting both CPU-only and GPU-accelerated systems, developed over the course of the Exascale Computing Project (ECP). While it began with several documents setting out its initial design, significant design changes occurred throughout its development. In some cases, these were anticipated: an early version used a simple consistency flag that was later replaced with a full-featured consistency protocol. In other cases, performance limitations and software and hardware changes prompted a redesign. Sequential communication tasks were parallelized; host-to-host MPI calls were replaced with GPU device-to-device MPI calls; more advanced algorithms such as Communication Avoiding LU and the Random Butterfly Transform (RBT) were introduced. Early choices that turned out to be cumbersome, error prone, or inflexible have been replaced with simpler, more intuitive, or more flexible designs. Applications have been a driving force, prompting a lighter weight queue class, nonuniform tile sizes, and more flexible MPI process grids. Of paramount importance has been building a portable library that works across several different GPU architectures – AMD, Intel, and NVIDIA – while keeping a clean and maintainable codebase. Here we explore the evolving design choices and their effects, both in terms of performance and software sustainability.

Gates, Mark↗
DOE OSTI2020

Microgrid Design Toolkit (MDT) User Guide. Software v1.3.

The Microgrid Design Toolkit (MDT) supports decision analysis for new ("greenfield") microgrid designs as well as microgrids with existing infrastructure. The current version of MDT includes two main capabilities. The first capability, the Microgrid Sizing Capability (MSC), is used to determine the size and composition of a new, grid connected microgrid in the early stages of the design process. MSC is focused on developing a microgrid that is economically viable when connected to the grid. The second capability is focused on designing a microgrid for operation in islanded mode. This second capability relies on two models: the Technology Management Optimization (TMO) model and Performance Reliability Model (PRM).

42 ENGINEERING↗
DOE OSTI2024

Predicting Missing Regions in Charged Particle Tracks Using a Sparse 3D Convolutional Neural Network

The 2x2 Demonstrator is a prototype detector for the Deep Underground Neutrino Experiment (DUNE)'s Near Detector. Both the 2x2 Demonstrator and the Near Detector itself will have inactive regions wherein there is no sensitivity to charge deposition and light signals that arise from charged particle interactions with liquid argon. In the 2x2, these inactive regions are positioned in-between the active detector modules, which introduces the challenge of inferring what charge signals ought to look like in these regions. This study explores the use of a Sparse 3D Convolutional Neural Network (ConvNet) to infer missing regions in charged particle tracks. Hits corresponding to energy depositions are voxelized into a three-dimensional (3D) grid for each track. Inactive regions within the tracks are replaced with a dense, rectangular 3D grid of voxels, ensuring consistent step sizes in X, Y, and Z directions. Voxels in these dense regions are initialized with an energy value of -1, indicating nonphysical energy or charge. The model is trained to predict which voxels should activate as part of the track and which should not, with the goal of eventually inferring the missing charge or energy values in these voxels. Results indicate that the model accurately predicts track voxels within ±1 unit in X, Y, or Z directions and effectively identifies non-track voxels, despite some overprediction. The approach shows promise in prediction of missing track regions with some accuracy.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

CITRC Cable Info

Working on a DOE project with university sub-contractors to build power grid simulations. We need the cable size and spacing to build accurate models. Other related documents on the physical set up and electrical components have been released, INL/MIS-20-59900 and INL/MIS-20-60003.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Predicting Missing Regions in Charged Particle Tracks Using a Sparse 3D Convolutional Neural Network

The 2x2 Demonstrator is a prototype detector for the Deep Underground Neutrino Experiment (DUNE)'s Near Detector. Both the 2x2 Demonstrator and the Near Detector itself will have inactive regions wherein there is no sensitivity to charge deposition and light signals that arise from charged particle interactions with liquid argon. In the 2x2, these inactive regions are positioned in-between the active detector modules, which introduces the challenge of inferring what charge signals ought to look like in these regions. This study explores the use of a Sparse 3D Convolutional Neural Network (ConvNet) to infer missing regions in charged particle tracks. Hits corresponding to energy depositions are voxelized into a three-dimensional (3D) grid for each track. Inactive regions within the tracks are replaced with a dense, rectangular 3D grid of voxels, ensuring consistent step sizes in X, Y, and Z directions. Voxels in these dense regions are initialized with an energy value of -1, indicating nonphysical energy or charge. The model is trained to predict which voxels should activate as part of the track and which should not, with the goal of eventually inferring the missing charge or energy values in these voxels. Results indicate that the model accurately predicts track voxels within 1 unit in X, Y, or Z directions and effectively identifies non-track voxels, despite some overprediction. The approach shows promise in prediction of missing track regions with some accuracy.

Utaegbulam, Hilary↗
DOE OSTI2024

Webinar: Modular Nuclear Power: The Microreactor Revolution

Many times smaller than traditional nuclear reactors, microreactors can be made in a factory and shipped wherever needed. This innovative approach to energy production has huge potential to bring electricity to remote areas, support regions affected by natural disasters, and provide right-sized, scalable power to our energy grid. Learn what they are, how they work, and their benefits in this webinar, featuring Piyush Sabharwall, distinguished staff nuclear research scientist and technical area lead for the DOE Office of Nuclear Energy's Microreactor R&D Program at Idaho National Laboratory.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

AGILE GETRS

AGILE GETRS provides an optimized implementation of DETRS that is faster than vendor optimized libraries for small (< 4000) problem sizes. These small problems sizes are of interest to researchers in the Grid Research Integration and Deployment Center (GRID-C) researching real-time power electronics simulations. In addition to our implementation of DETRS, this code provides a test using Google Benchmark which generates matrices of various sizes, fills them with random numbers, runs GETRS, verifies the solution is correct and reports average runtime and other performance metrics.

Hahn, Steven [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Designing Yield Grids for Airborne Nuclear Explosions to Maintain Decrements in Rayleigh Waveform Amplitudes

This document details a method for selecting grid points that sample the yields of airborne, low yield nuclear explosions that constrains the relative size of seismic waveforms sourced by explosions that are parameterized by this grid. In detail, we impose this constraint so that the vertical component of ground displacement that is sourced by explosions with a fixed height of burst (HoB) at grid point k change by a fixed factor of a when compared to the waveform amplitudes sourced by explosions at grid point k + 1, as observed at a fixed standoff distance. This gridding thereby preserves amplitude ratios along the grid so that they are agnostic to k. We assert that this method is useful in tests that apply algorithms against synthetic or augmented datasets for which the performance of the algorithm must be measured against source size, and waveform amplitudes must change uniformly (e.g., Berg et al, 2025).

58 GEOSCIENCES↗
DOE OSTI2021

Analytical Expression for DC Link Capacitor Current in a Cascaded H-Bridge Multi-Level Active Front-End Converter

Medium-voltage grid-tied systems often use a cascaded H-bridge multi-level active front-end. In this converter, dc link bus capacitors play an important role in stabilizing the converter and enabling both active and reactive power injections. The present work provides analytical expressions for the capacitor current, which are essential for optimizing system design (especially capacitor size vs. lifetime). Then, the expression is incorporated into the grid connected bidirectional power system model. Consequently, this work contributes to the guiding principles to choose accurate dc link capacitor ratings against grid-side power delivery requirements. The analytical results have been validated with detailed simulations and hardware results.

cascaded H-bridge, carrier-based PWM, dc-link capa↗
DOE OSTI2024

Parallel computing for power system climate resiliency: Solving a large-scale stochastic capacity expansion problem with mpi-sppy

Here we propose a nodal stochastic generation and transmission expansion planning model that incorporates the output from high-resolution global climate models through load and generation availability scenarios. We implement our model in Pyomo and perform computational studies on a realistically-sized test case of the California electric grid in a high performance computing environment. We propose model reformulations and algorithm tuning to efficiently solve this large problem using a variant of the Progressive Hedging Algorithm. We utilize the parallelization capabilities and overall versatility of mpi-sppy, exploiting its hub-and-spoke architecture to concurrently obtain inner and outer bounds on an optimal expansion plan. Initial results show that instances with 360 representative days on a system with over 8,000 buses can be solved to within 5% of optimality in under 4 h of wall clock time, a first step towards solving a large-scale power system expansion planning problem across a wide range of climate-informed operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

A cancellation problem in hybrid particle-in-cell schemes due to finite particle size

The quasi-neutral hybrid particle-in-cell algorithm with kinetic ions and fluid electrons is a popular model to study multi-scale problems in laboratory, space, and astrophysical plasmas. Here, it is shown that the different spatial discretizations of ions as finite-spatial-size particles and electrons as a grid-based fluid can lead to significant numerical wave dispersion errors in the long wavelength limit (kd i «1, where k is the wavenumber and di is the ion skin-depth). The problem occurs when high-order particle-grid interpolations, or grid-based smoothing, spreads the electric field experienced by the ions across multiple spatial cells and leads to inexact cancellation of electric field terms in the total (ion + electron) momentum equation. Practical requirements on the mesh spacing Δx/d i are suggested to bound these errors from above. The accuracy impact of not respecting these resolution constraints is shown for a non-linear shock problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

LANL Institutional Computing Highlight: Local Time-Stepping in Global to Regional Ocean Modeling (Annual Report) [Slides]

Multi-resolution global ocean simulations are possible with the Model for Prediction Across Scales-Ocean (MPAS-Ocean), although the model is currently constrained to a uniform time step. The time step is determined by the size of the smallest cell in the grid, according to the CFL condition. To overcome this issue, a local time stepping (LTS) scheme has been designed for MPAS-Ocean. We have made progress on implementing the LTS scheme, as well as quantitative results and performance tests.

58 GEOSCIENCES↗
DOE OSTI2022

Deep Electrification Analysis: The Role of the U.S. Power Grid for Sustainable Transportation

This project attempts to quantify the size of electric generation for the entire nation to transition from a fossil fuel based transportation sector to a zero GHG emission-based energy source. The scope of this analysis is limited to decarbonizing the transportation sector, leaving the remaining sectors, such as power (for those that are still fossil based), industry, and building sectors, for later phases of study. The study year for this analysis is 2050, with expected escalation in transportation services and naturally occurring evolutions in the electric power sector and the entire economy. This analysis uses the projections of the Energy Information Administration’s (EIA’s) Annual Energy Outlook (AEO 2020) Reference Case for study year 2050 [EIA/AEO2020] as a base-case. The transportation sector is disaggregated by the following modes and classes: (1) on-road (divided into light-duty, medium-duty, heavy-duty vehicles), (2) aviation, (3) maritime, and (4) rail. The decarbonization case was based on only 2 pathways: (1) electrification of on-road transportation except for 30% of heavy-duty vehicles, and (2) power-to-liquid for the remaining transportation modes. The study estimated for 11 US regions what the additional wind and storage capacities requirements are to replace the fossil-based fuels with renewable wind capacity. Considered were the utilization of the existing idle capacity particularly during the load valley at night and any additional new generation capacity in EIA projections for the reference case. To balance the additional wind capacity required significant energy storage capabilities which were estimated in terms of power capacity (GW) and energy capacity (GWh). The paper further characterizes the energy requirements by a relation of power capacity to duration, allowing the analyst to gain insights into what the best technology portfolio might be to meet the new balancing or flexibility needs.

24 POWER TRANSMISSION AND DISTRIBUTION↗