SEARCH · Engineering Papers

Results for “large 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 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Speed to Power: Solutions for Accelerating Large Load Connections

Rapid growth in demand from data centers and other large loads is creating a range of new challenges for electricity planners, investors, system operators, and regulators, leading to bottlenecks that have slowed connection of large loads to the electric grid. In response, innovative solutions for accelerating large load connections are beginning to emerge across the U.S. Drawing on an extensive document and literature review, this report identifies more than 40 potential solutions for accelerating large load connections, organized into five functional areas: load forecasting, interconnection, resource planning and procurement, markets and operations, and cost allocation and ratemaking. The five functional areas provide a framework for organizing challenges and solutions to large load connection bottlenecks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS2010

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Inviscid Fluxes

Cell-centered and node-centered approaches have been compared for unstructured finite-volume discretization of inviscid fluxes. The grids range from regular grids to irregular grids, including mixed-element grids and grids with random perturbations of nodes. Accuracy, complexity, and convergence rates of defect-correction iterations are studied for eight nominally second-order accurate schemes: two node-centered schemes with weighted and unweighted least-squares (LSQ) methods for gradient reconstruction and six cell-centered schemes two node-averaging with and without clipping and four schemes that employ different stencils for LSQ gradient reconstruction. The cell-centered nearest-neighbor (CC-NN) scheme has the lowest complexity; a version of the scheme that involves smart augmentation of the LSQ stencil (CC-SA) has only marginal complexity increase. All other schemes have larger complexity; complexity of node-centered (NC) schemes are somewhat lower than complexity of cell-centered node-averaging (CC-NA) and full-augmentation (CC-FA) schemes. On highly anisotropic grids typical of those encountered in grid adaptation, discretization errors of five of the six cell-centered schemes converge with second order on all tested grids; the CC-NA scheme with clipping degrades solution accuracy to first order. The NC schemes converge with second order on regular and/or triangular grids and with first order on perturbed quadrilaterals and mixed-element grids. All schemes may produce large relative errors in gradient reconstruction on grids with perturbed nodes. Defect-correction iterations for schemes employing weighted least-square gradient reconstruction diverge on perturbed stretched grids. Overall, the CC-NN and CC-SA schemes offer the best options of the lowest complexity and secondorder discretization errors. On anisotropic grids over a curved body typical of turbulent flow simulations, the discretization errors converge with second order and are small for the CC-NN, CC-SA, and CC-FA schemes on all grids and for NC schemes on triangular grids; the discretization errors of the CC-NA scheme without clipping do not converge on irregular grids. Accurate gradient reconstruction can be achieved by introducing a local approximate mapping; without approximate mapping, only the NC scheme with weighted LSQ method provides accurate gradients. Defect correction iterations for the CC-NA scheme without clipping diverge; for the NC scheme with weighted LSQ method, the iterations either diverge or converge very slowly. The best option in curved geometries is the CC-SA scheme that offers low complexity, second-order discretization errors, and fast convergence.

Diskin, Boris↗

NASA NTRS1989

A numerically exact full wave packet approach to molecule-surface scattering

A numerically exact spectral method for solving the time-dependent Schroedinger equation in spherical coordinates is described. The angular dependence of the wave function is represented on a two-dimensional grid of evenly spaced points. The fast Fourier transform algorithm is used to transform between the angle space representation of the wave function and its conjugate representation in momentum space. The time propagation of the wave function is evaluated using an expansion of the time evolution operator as a series of Chebyshev polynomials. Calculations performed for a model system representing H2 scattering from a rectangular corrugated surface yield transition probabilities that are in excellent agreement with those obtained using the close-coupling wave packet (CCWP) method. However, the new method is found to require substantially more computation time than the CCWP method because of the large number of grid points needed to represent the angular dependence of the wave function and the variation in the number of terms required in the Chebyshev representation of the time evolution operator.

Mowrey, R. C.↗

NASA NTRS1998

MOLA Topography of Impact Basins in the Northern Hemisphere of Mars

Coverage of the northern hemisphere of Mars by the Mars Orbiter Laser Altimeter (MOLA) during the aerobraking hiatus and the two Science Phasing Operation periods provides improved definition and characterization of large impact basins. Gridded MOLA data show the Utopia Basin has a pronounced bowl-like structure, as opposed to the interior rises suggested by the earlier USGS DEM. The elevation structure is concentric about the basin center as mapped by McGill. In particular, the proposed inner ring closely follows the -4 km contour over much of the southern, western and northwestern sides. Higher topography along portions of the dichotomy boundary aligns with the basin's outer ring. High topography in the polar region also occurs where the outer ring should lie, raising the possibility that perhaps some of the polar topography is due to basin structure as well as ice. Two MOLA passes near Phison Rupes provide evidence for a large "stealth" hole where Viking imagery show little evidence of any major structure. The 2 km deep, 600 km wide depression at 31OW, 3ON is as large as the Cassini impact basin 1000 km to the SW. While Cassini is easily recognized in image data, the "MOLA Hole" is not. If this depression is a deeply eroded and buried impact basin (as perhaps suggested by a decrease in the crater density and somewhat smoother terrain than in adjacent areas), it is not clear why it has managed to maintain its great depth. In Tempe at the dichotomy boundary a 300 km wide impact basin is revealed by pronounced bowl-like topography centered at 87W, 47N, even though only about 1/3 of the basin rim structure is obvious. The basin lies on a sloping boundary zone, with the more buried N rim up to 2 km below the rugged S rim. A similar N-S asymmetry in basin ring structure occurs for the much larger Isidis Basin, where the S rim rises 6 km but the subdued N rim rises barely 2 km above the floor. There is essentially no topographic expression of the main ring in the NE quadrant of Isidis where, if it exists, it lies below Hesperian-age plains.

Frey, Herbert↗

NASA NTRS2000

Ordering Unstructured Meshes for Sparse Matrix Computations on Leading Parallel Systems

The ability of computers to solve hitherto intractable problems and simulate complex processes using mathematical models makes them an indispensable part of modern science and engineering. Computer simulations of large-scale realistic applications usually require solving a set of non-linear partial differential equations (PDES) over a finite region. For example, one thrust area in the DOE Grand Challenge projects is to design future accelerators such as the SpaHation Neutron Source (SNS). Our colleagues at SLAC need to model complex RFQ cavities with large aspect ratios. Unstructured grids are currently used to resolve the small features in a large computational domain; dynamic mesh adaptation will be added in the future for additional efficiency. The PDEs for electromagnetics are discretized by the FEM method, which leads to a generalized eigenvalue problem Kx = AMx, where K and M are the stiffness and mass matrices, and are very sparse. In a typical cavity model, the number of degrees of freedom is about one million. For such large eigenproblems, direct solution techniques quickly reach the memory limits. Instead, the most widely-used methods are Krylov subspace methods, such as Lanczos or Jacobi-Davidson. In all the Krylov-based algorithms, sparse matrix-vector multiplication (SPMV) must be performed repeatedly. Therefore, the efficiency of SPMV usually determines the eigensolver speed. SPMV is also one of the most heavily used kernels in large-scale numerical simulations.

Oliker, Leonid↗

NASA NTRS1994

Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation

In this paper finite-difference solutions of the Helmholtz equation in an open domain are considered. By using a second-order central difference scheme and the Bayliss-Turkel radiation boundary condition, reasonably accurate solutions can be obtained when the number of grid points per acoustic wavelength used is large. However, when a smaller number of grid points per wavelength is used excessive reflections occur which tend to overwhelm the computed solutions. Excessive reflections are due to the incompability between the governing finite difference equation and the Bayliss-Turkel radiation boundary condition. The Bayliss-Turkel radiation boundary condition was developed from the asymptotic solution of the partial differential equation. To obtain compatibility, the radiation boundary condition should be constructed from the asymptotic solution of the finite difference equation instead. Examples are provided using the improved radiation boundary condition based on the asymptotic solution of the governing finite difference equation. The computed results are free of reflections even when only five grid points per wavelength are used. The improved radiation boundary condition has also been tested for problems with complex acoustic sources and sources embedded in a uniform mean flow. The present method of developing a radiation boundary condition is also applicable to higher order finite difference schemes. In all these cases no reflected waves could be detected. The use of finite difference approximation inevita bly introduces anisotropy into the governing field equation. The effect of anisotropy is to distort the directional distribution of the amplitude and phase of the computed solution. It can be quite large when the number of grid points per wavelength used in the computation is small. A way to correct this effect is proposed. The correction factor developed from the asymptotic solutions is source independent and, hence, can be determined once and for all. The effectiveness of the correction factor in providing improvements to the computed solution is demonstrated in this paper.

Tam, Christopher K. W.↗

DOE OSTI2021

Economic Analysis of an Electric Thermal Energy Storage System Using Solid Particles for Grid Electricity Storage

As renewable power generation becomes the mainstream new-built energy source, energy storage will become an indispensable need to complement the uncertainty of renewable resources to firm the power supply. When phasing out fossil-fuel power plants to meet the carbon neutral utility target in the midcentury around the world, large capacity of energy storage will be needed to provide reliable grid power. The renewable power integration with storage can support future carbon-free utility and has several significant impacts including increasing the value of renewable generation to the grid, improving the peak-load response, and balancing the electricity supply and demand. Long-duration energy storage (10–100 hours duration) can potentially complement the reduction of fossil-fuel baseload generation that otherwise would risk grid security when a large portion of grid power comes from variable renewable sources. Current energy storage methods based on pumped storage hydropower or batteries have many limitations. Thermal energy storage (TES) has unique advantages in scale and siting flexibility to provide grid-scale storage capacity. A particle-based TES system has promising cost and performance for the future growing energy storage needs. This paper introduces the system and components required for the particle TES to be technically and economically competitive. A technoeconomic analysis based on preliminary component designs and performance shows that the particle TES integrated with an efficient air-Brayton combined cycle power system can provide power for several days by low-cost, high-performance storage cycles. It addresses grid storage needs by enabling large-scale grid integration of intermittent renewables like wind and solar, thereby increasing their grid value. The design specifications and cost estimations of major components in a commercial scale system are presented in this paper. The cost model provides insights for further development and cost comparison with competing technologies.

27 ARPA - Advanced Research Projects Agency-Energy↗

DOE OSTI2021

Economic Analysis of a Novel Thermal Energy Storage System Using Solid Particles for Grid Electricity Storage: Preprint

As renewable power generation becomes the mainstream new-built energy source, energy storage will become an indispensable need to complement the uncertainty of renewable resources to firm the power supply. When phasing out fossil-fuel power plants to meet the carbon neutral utility target in the midcentury around the world, large capacity of energy storage will be needed to provide reliable grid power. The integration of renewable power and storage to support future carbon-free utility has several significant and positive impacts including expanding the renewable generation into the grid, improving the peak-load response, and balancing the electricity supply and demand. Long-duration energy storage (10–100 hours duration) can potentially complement the reduction of fossil-fuel baseload generation that otherwise would risk grid security when a large portion of grid power comes from variable renewable sources. Current energy storage methods based on pumped storage hydropower or batteries have many limitations. Thermal energy storage (TES) has unique advantages in scale and siting flexibility to provide grid-scale storage capacity. A particle-based TES system has promising cost and performance for the future growing energy storage needs. This paper introduces the system and components required for the particle TES to be technically and economically competitive. A technoeconomic analysis based on preliminary component designs and performance shows that the particle TES integrated with an efficient air-Brayton combined cycle power system can provide power for several days by low-cost, high-performance storage cycles. It addresses grid storage needs by enabling large-scale grid integration of intermittent renewables like wind and solar, thereby increasing their grid value. The design specifications and cost estimations of major components in a commercial scale system are presented in this paper. The cost model provides insights for further development and cost comparison with competing technologies.

27 ARPA - Advanced Research Projects Agency-Energy↗

DOE OSTI2026

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗

DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗

DOE OSTI2022

Sodium trithiocarbonate cathode for high-performance sodium–sulfur batteries

The high abundance and low cost of sodium and sulfur make room-temperature sodium-sulfur (RT Na-S) batteries an attractive technology compared to the current lithium-ion batteries for large-scale grid-storage applications. However, the commercialization of RT Na-S batteries is impeded by the slow kinetics of Na-S chemistry, severe sodium polysulfide shuttling, and uncontrollable growth of dendritic Na. Herein, sodium trithiocarbonate (Na 2 CS 3 ) is applied as a cathode material to facilitate concurrent improvement in both electrodes, leading to a high-rate performance with extended cycle life. The conductive characteristic of the carbon-sulfur resonance bond provides fast ion and electron transfer throughout the cathode, resulting in superior electrochemical reactivity. At the cathode, the presence of Na2CS3 forms an oligomer-structured layer to suppress the dissolution and shuttling of active materials. Meanwhile, when small portions of Na 2 CS 3 intermediates migrate to the anode, a stable solid electrolyte interphase (SEI) layer with uniform Na-ion flux is formed, enabling improved Na stripping and plating performance. A series of electrochemical and materials characterization, accompanied by density functional theory calculations demonstrate that Na 2 CS 3 is a promising candidate to realize high-rate performance long cycle life RT Na-S batteries.

25 ENERGY STORAGE↗

DOE OSTI2024

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗

DOE OSTI2025

A Participation Factor-Based Approach for Defining the EMT Model Boundary for Power System Simulations with Inverter-Based Resources

The increasing penetration of inverter-based resources (IBRs) introduces new challenges to power system simulations, particularly with the emergence of fast electromagnetic transient (EMT) dynamics and sub-synchronous oscillations (SSO) that require time-consuming EMT simulations. To reduce the time cost for simulating a large-scale power grid with IBRs, this paper proposes a novel participation factor-based approach for defining a critical zone for detailed EMT modeling and simulations, which includes the IBRs, synchronous generators, and the network components participating significantly in simulated contingencies. Both model-based and response-based methods are introduced for the estimation of participation factors (PFs). The case study on the 240-bus Western Electricity Coordinating Council (WECC) system demonstrates that the EMT zone determined by the proposed approach can effectively capture power system dynamics involving IBRs.

EMT simulation↗

DOE OSTI2024

Modeling and Automation Framework for High IBR Integration in Large-Scale Power Systems

The increasing prevalence of power electronics- interfaced renewable generation sources is leading to a gradual replacement of traditional thermal generation-based synchronous machines. In this context, the modeling of a large-scale power grid that incorporates a significant number of inverter-based resources is crucial for understanding the dynamics and effects of these resources on the power system. This study investigates the positive sequence model of grid-following and grid-forming inverters. Additionally, this work explores the integration of distributed energy resources using population as an indicator of their relative geographic locations. To address challenge to integrate these inverter based resources into a realistic grid of the US Western interconnection, automation scripts are developed to streamline the process of replacing conventional generators with grid-following and grid-forming inverters, as well as allocating distributed energy resources. Different penetration levels of these inverters are considered, and their frequency regulation support following a disturbance is compared through dynamic simulations.

Lyu, Xue↗

DOE OSTI2024

Benchmark High-Fidelity EMT Models for Power Grid with PV Plants

In recent times electromagnetic transient (EMT) modeling tools have been identified as one of the most important requirements in replicating, analyzing, and investigating the dynamics of the power grid with photovoltaic (PV) plants. However, there are no benchmark models for power grid with PVs to investigate emerging challenges with higher penetration of PVs (like trips and momentary cessations during faults from a region far away). To this end, in this paper, synthetic benchmark high-fidelity EMT dynamic models of power grid with large-scale PV plants are presented. The models are developed in PSCAD and PSCAD/Fortran. Simulation results for different use cases (events) and scenarios are presented.

Marthi, Phani Ratna Vanamali↗

DOE OSTI2022

PowerNet: Multi-agent Deep Reinforcement Learning for Scalable Powergrid Control

This paper develops an efficient multi-agent deep reinforcement learning algorithm for cooperative controls in powergrids. Specifically, we consider the decentralized inverter-based secondary voltage control problem in distributed generators (DGs), which is first formulated as a cooperative multi-agent reinforcement learning (MARL) problem. We then propose a novel on-policy MARL algorithm, PowerNet, in which each agent (DG) learns a control policy based on (sub-)global reward but local states and encoded communication messages from its neighbors. Motivated by the fact that a local control from one agent has limited impact on agents distant from it, we exploit a novel spatial discount factor to reduce the effect from remote agents, to expedite the training process and improve scalability. Furthermore, a differentiable, learning-based communication protocol is employed to foster the collaborations among neighboring agents. In addition, to mitigate the effects of system uncertainty and random noise introduced during on-policy learning, we utilize an action smoothing factor to stabilize the policy execution. To facilitate training and evaluation, we develop PGSim, an efficient, high-fidelity powergrid simulation platform. Here, experimental results in two microgrid setups show that the developed PowerNet outperforms the conventional model-based control method, as well as several state-of-the-art MARL algorithms. The decentralized learning scheme and high sample efficiency also make it viable to large-scale power grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Enhancing Active Distribution Systems Resilience by Fully Distributed Self-Healing Strategy

Distributed restoration can exploit smart grid technologies to enhance the resilience of active distribution networks toward a self-healing smart grid. However, the large number of decision variables, especially the binary ones for reconfiguration, bring challenges to developing scalable distributed distribution service restoration (DDSR) strategies. This paper proposes a fully distributed solution procedure based on the alternating direction method of multipliers (ADMM) for mixed-integer programming problems and applies to develop the DDSR framework. The method consists of relax-drive-polish phases, 1) relaxing binary variables, and applying the convex ADMM as a warm start; 2) driving the solutions toward Boolean values through a proximal operator; 3) fixing the obtained binding binary variables and solving the rest of the problem to polish results and achieve a high-quality suboptimal solution. Then, an autonomous clustering strategy and consensus ADMM are integrated with the proposed method to realize the fully distributed cluster-based framework of DDSR. This framework can first determine DER scheduling and switch status for reconfiguration to energize the out-of-service areas from local faults, and then provide the load restoration solution in a distributed manner for total blackouts in large-scale distribution networks. Furthermore, the effectiveness and scalability of the proposed DDSR framework are demonstrated through testing on the IEEE 123-node, IEEE 8500-node, and synthetic 100k-node test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

HEMP Transformer Defense Through Power Electronics

High altitude electromagnetic pulses and geo-magnetic disturbances have the potential to severely impact the electric power grid by damaging large power transformers and causing severe power quality issues. This impact comes as a result of a quasi-static bias induced on transmission lines by geomagnetically induced currents which saturate magnetic components in the electric power system. This paper introduces the concept of utilizing a h-bridge inverter on the neutral of a LPT to inject a DC bias equivalent voltage onto the neutral side of the transformer windings. This biasing floats the transformer windings, eliminating the effect of the DC current and keeping the transformer from saturating. Schematic diagrams will be presented, along with simulation model data using Typhoon and PLECS, and finally test results from a benchtop hardware test.

24 POWER TRANSMISSION AND DISTRIBUTION↗