Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid stability”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗
DOE OSTI2019

Global Phase and Magnitude Synchronization of Coupled Oscillators With Application to the Control of Grid-Forming Power Inverters

In this paper, we explore a new approach to synchronization of coupled oscillators. In contrast to the celebrated Kuramoto model, we do not work in polar coordinates and do not consider oscillations of fixed magnitude. We propose a synchronizing feedback based on relative state information and local measurements that induces consensus-like dynamics. We show that, under a mild stability condition, the combination of the synchronizing feedback with a decentralized magnitude control law renders the oscillators' almost globally asymptotically stable with respect to set points for the phase shift, frequency, and magnitude. We apply these result to rigorously solve an open problem in control of inverter-based ac power systems. In this context, the proposed control strategy can be implemented using purely local information, induces a grid-forming behavior, and ensures that a network of ac power inverters is almost globally asymptotically stable with respect to a prespecified solution of the ac power-flow equations. Moreover, we show that the controller exhibits a droop-like behavior around the standard operating point, thus, making it backward compatible with the existing power system operation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Cycle-to-Cycle Variations in Multi-Cycle Engine RANS Simulations

Reynolds-averaged Navier-Stokes (RANS) modeling is expected to deliver an ensemble-averaged result for the majority of turbulent flows. This could lead to the conclusion that multi-cycle internal combustion engine (ICE) simulations performed using RANS must exhibit a converging numerical solution after a certain number of consecutive cycles. However, for some engine configurations unsteady RANS simulations are not guaranteed to deliver an ensemble-averaged result. In this paper it is shown that, when using RANS modeling to simulate multiple engine cycles, the cycle-to-cycle variations (CCV) generated from different initial conditions at each cycle are not damped out even after a large number of cycles. A single-cylinder GDI research engine is simulated using RANS modeling and the numerical results for 20 consecutive engine cycles are evaluated for two specific operating conditions. One condition is characterized by stoichiometric operation and stable combustion (COVIMEP < 2%) and the other features dilute combustion resulting in misfire events and much higher COVIMEP. In both cases, multi-cycle RANS simulation results show cyclic variability. An in-depth analysis of the most significant physical and chemical quantities highlights that CCV are caused primarily by the variability of the in-cylinder flow. For stoichiometric combustion this does not greatly affect flame propagation, and the resulting fluctuations of the pressure traces are narrow. However, in the event of dilute combustion, the variability of the flow gains more importance, therefore delivering large cyclic fluctuations of in-cylinder pressure. Even though unsteady RANS simulations are not expected to predict as much variability in engine flows as LES, in particular for stable cases, they can still be used to capture typical combustion stability features in an internal combustion engine. This paper shows that the occurrence of CCV when using multi-cycle RANS modeling is not a numerical artifact or the effect of changes in the computational grid from cycle to cycle. Cold-flow analysis reveals that the variability of the simulated flow is intrinsic for the unsteady RANS approach and can only be damped out by increasing numerical viscosity, either by increasing the cell size or by adding up-winding. This is proved in a simplified non-ICE case such as the simulation of the vortex shredding from a cylinder in cross flow. Also, the effect of numerical viscosity on CCV is shown for multi-cycle engine simulations by using a simplified combustion model. Finally, it is shown that multi-cycle RANS can successfully predict the effect of changing specific engine parameters, such as the properties of the ignition system, on cyclic variability and combustion stability.

Scarcelli, Riccardo↗
DOE OSTI2017

Stability and scalability of the CMS Global Pool: Pushing HTCondor and glideinWMS to new limits

The CMS Global Pool, based on HTCondor and glideinWMS, is the main computing resource provisioning system for all CMS workflows, including analysis, Monte Carlo production, and detector data reprocessing activities. The total resources at Tier-1 and Tier-2 grid sites pledged to CMS exceed 100,000 CPU cores, while another 50,000 to 100,000 CPU cores are available opportunistically, pushing the needs of the Global Pool to higher scales each year. These resources are becoming more diverse in their accessibility and configuration over time. Furthermore, the challenge of stably running at higher and higher scales while introducing new modes of operation such as multi-core pilots, as well as the chaotic nature of physics analysis workflows, places huge strains on the submission infrastructure. This paper details some of the most important challenges to scalability and stability that the CMS Global Pool has faced since the beginning of the LHC Run II and how they were overcome.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A Review on Protection Challenges on Transmission Lines Connected to Inverter Based Resources

The increased penetration of inverter-based resources (IBRs) connected to transmission grids is challenging how traditional electric power systems operate. Such systems were originally designed based on the dominant presence of synchronous generators. Some main challenges include the increased concerns about system stability due to the reduced system inertia caused by IBRs and the reduced protection reliability caused by the non-deterministic fault current response of IBRs. The protection relays’ algorithms have been designed and validated based on well-known fault current characteristics of synchronous generations, while the short-circuit responses of IBRs heavily depend on their vendor-specific control schemes, which are mainly set to protect inverter semiconductor switches and DC components. The magnitude of short current that IBRs contribute is typically between 1.1–1.5 times the rated current that may behave unstably and incoherently with the voltages. The control scheme of the IBRs also controls the active (P) and reactive (Q) support during system voltage sag and swell conditions. These characteristics challenge not only protective relay schemes but also existing short-circuit software tools to correctly model IBRs for fault studies and protection coordination analyses. Utilities, university researchers, research institutes, and relay manufacturers have performed studies to investigate different aspects of high IBR penetration. Almost all studies support 1) creating international and national standards to regulate the IBRs, 2) revising engineering tools to address shortcomings, and 3) adjusting the protective relay settings to resolve protection problems. This paper summarizes the studies and research conducted in recent years to quantify the impact of IBR penetration on transmission line protection and to identify potential solutions to improve the dependability and reliability of grid protection systems.

14 SOLAR ENERGY↗
DOE OSTI2017

Alfvén eigenmode stability and critical gradient energetic particle transport using the Trapped-Gyro-Landau-Fluid model

The Trapped-Gyro-Landau-Fluid (TGLF) transport model is a physically realistic and comprehensive theory based on a local quasilinear transport model fitted to linear and nonlinear GYRO gyrokinetic simulations [Staebler et al., Phys. Plasmas 14, 55909 (2007)]. This work presents the first use of the TGLF model to treat low-n Alfvén eigenmode (AE) stability and energetic particle (EP) transport. TGLF accurately recovers the local GYRO toroidicity-induced AE (TAE) and energetic particle mode (EPM) linear growth and frequency rates for a fusion alpha case. With a very high grid resolution, TGLF can quickly find the critical EP pressure gradient profile for stiff EP transport based on an AE linear threshold given the background thermal plasma profiles in DIII-D. The TGLF critical gradient profile using the recipe γ AE =0, that is the linear AE growth rate without additional driving rates from the background plasma gradients, matches the more expensive linear GYRO results with a single worst toroidal mode number n. TGLF can easily find the minimum critical gradient profile with testing multiple ns. From a database of runs using a newly developed TGLFEP code, a rough but insightful parametric “power law” scaling for critical EP beta is demonstrated. An important toroidal stabilization condition on the EP pressure gradient p EP /L$EP\atop{p}$ drive is isolated: R/L$EP\atop{p}$ > C R ~ 3, where L$EP\atop{p}$ is the EP pressure gradient length and R is the tokamak major radius. This paper also demonstrates that relaxation of the fixed slowing down EP profile shape approximation often used to find the critical EP density profile has little effect on the resulting EP transport. The single EP species critical gradient model is generalized to handle two EP species.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evaluation of Optimal Net Load Management in Microgrids Using Hardware-in-the-Loop Simulation

This paper presents the performance evaluation of a net load management (NLM) engine that balances load and generation in an isolated community to power a critical facility after a grid interruption event (e.g., the loss of a large generation unit). This NLM engine is particularly important for microgrid systems because it provides a high-speed, cost-optimal control solution to coordinate grid-forming inverters and to dispatch grid-following inverters and deferrable loads in microgrid systems to enhance grid resilience and reliability. The NLM algorithm cost-optimally dispatches the grid-following inverters and deferrable loads based on the demanded power and load priorities, and the grid-forming inverters use droop control to form system voltages and share active and reactive power. A controller-hardware-in-the-loop platform is developed to evaluate the control performance of the NLM algorithm with two sequential contingency events of lost generation units. The experimental results indicate that the NLM engine can maintain system stability, achieve the targeted system voltage and frequency, and balance load and generation to serve the critical facility with improved system resilience and reliability.

grid-following inverter↗
NASA NTRS2006

Using Flow Regime Lightning and Sounding Climatologies to Initialize Gridded Lightning Threat Forecasts for East Central Florida

Each morning, the forecasters at the National Weather Service in Melbourne, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (http://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in East Central Florida, especially during the warm season months of May September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Currently, the forecasters create each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent is to improve consistency between forecasters while allowing them to focus on the mesoscale detail of the forecast, ultimately benefiting the end-users of the product. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) in which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities, or number of strikes per specified area. The soundings used to determine the flow regimes were taken at Miami (MIA), Tampa (TBW), and Jacksonville (JAX), FL, and the lightning data for the strike densities came from the National Lightning Detection Network (NLDN). The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at FSU and NWS TAE provided this data and supporting software for the work performed by the AMU.

Lambert, Winifred↗
NASA NTRS2007

Using Flow Regime Lightning and Sounding Climatologies to Initialize Gridded Lightning Threat Forecasts for East Central Florida

Each morning, the forecasters at the National Weather Service in Melbourne, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (httl://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in East Central Florida, especially during the warm season months of May September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Until recently, the forecasters created each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent was to improve consistency between forecasters while allowing them to focus on the mesoscale detail of the forecast. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) in which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities, or number of strikes per specified area. The soundings used to determine the flow regimes were taken at Miami (MIA), Tampa (TBW), and Jacksonville (JAX), FL, and the lightning data for the strike densities came from the National Lightning Detection Network (NLDN). The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at FSU and NWS TAE provided this data and supporting software for the work performed by the AMU.

Lambert, Winifred↗
DOE OSTI2018

On a Future for Smart Inverters with Integrated System Functions

The modern electric power grid is gradually evolving to include high penetration power electronics in a distributed architecture which will pose pressing new challenges on power quality, system control, and operation. The state-of-the-art smart inverters for solar photovoltaics may not be sufficient to address some of the stability and control challenges such as harmonic resonances due to the interactions of multiple power electronics converters. This paper discusses a future of smart inverters with integrated system functions and its associated challenges, and identifies smart inverter research trends and features in the context of flexible control architectures. Future power electronics research will shift from hardware solutions to more control solutions with integrated system functions under distributed power grid scenarios with small-scale power electronics interfaced with distributed generations and loads.

Xue, Yaosuo↗
NASA NTRS1980

Rectangular beam (5 X 40 cm multipole ion source)

A 5 x 40 cm rectangular-beam ion source was designed and fabricated. A multipole field configuration was used to facilitate design of the modular rectangular chamber, while a three-grid ion optics system was used for increased ion current densities. For the multipole chamber, a magnetic integral of 0.000056 Tesla-m was used to contain the primary electrons. This integral value was reduced from the initial design value, with the reduction found necessary for discharge stability. The final value of magnetic integral resulted in discharge losses at typical operating conditions which ranged from 600 to 1000 eV/ion, in good agreement with the design value of 800 eV/ion. The beam current density at the ion optics was limited to about 3.2 mA/sq cm at 500 eV and to about 3.5 mA/sq cm at 1000 ev. The effects of nonuniform ion current, dimension tolerance, and grid thermal warping were considered. The use of multiple rectangular-beam ion sources to process wider areas than would be possible with a single source (approx. 40 cm) was also studied. Beam profiles were surveyed at a variety of operating conditions and the results of various amounts of beam overlap calculated.

Haynes, C. M.↗
DOE OSTI2023

Understanding the Role of Zinc Hydroxide Sulfate and its Analogues in Mildly Acidic Aqueous Zinc Batteries: A Review

Abstract Mildly acidic aqueous zinc batteries (AZBs) have attracted tremendous attention for grid storage applications with the expectation to tackle the issues of Li‐ion batteries on high cost and poor safety. However, the performance, particularly energy density and cycle stability of AZBs are still unsatisfactory when compared with LIBs. To help the development of AZBs, a lot of effort have been made to understand the battery reaction mechanisms and precedent microscopic and spectroscopic analyses have shown flake‐like large particles of zinc hydroxide sulfate (ZHS) and its analogues formed on the surfaces of cathodes and anodes in sulfate and other electrolyte systems during cycling. However, because of the complexity of the thermodynamics and kinetics of aqueous reactions to understand different battery conditions, controversies still exist. This article will review the roles of ZHS discussed in recent representative references aiming to shine light on the fundamental mechanisms of AZBs and pave ways to further improve the battery performance.

25 ENERGY STORAGE↗
NASA NTRS1981

Two-dimensional integrating matrices on rectangular grids

The use of integrating matrices in solving differential equations associated with rotating beam configurations is examined. In vibration problems, by expressing the equations of motion of the beam in matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the spatial dependence is removed from the governing partial differential equations and the resulting ordinary differential equations can be cast into standard eigenvalue form. Integrating matrices are derived based on two dimensional rectangular grids with arbitrary grid spacings allowed in one direction. The derivation of higher dimensional integrating matrices is the initial step in the generalization of the integrating matrix methodology to vibration and stability problems involving plates and shells.

Lakin, W. D.↗
NASA NTRS1983

Numerical Simulation for Droplet Combustion Using Lagrangian Hydrodynamics

A predictive model of spray combustion must incorporate models for the wide variety of physical environments in a practical combustor. In regions where droplets are closely spaced, combustion resembles a diffusion flame; where they are well separated, an envelope or wake flame results. The relative velocity field between the fuel droplets and oxidizer in influences boundary layer development about the droplet, recirculating flow patterns, and droplet shape and stability. A model must encompass these interacting temporal and spatial effects as well as complicated combustor boundaries. The objective of the current work is to develop the triangular gridding method for describing the individual and collective properties of vaporizing and burning fuel droplets.

Fritts, M. J.↗
NASA NTRS1988

External store carriage loads prediction at transonic speeds

A computational method for prediction of external store carriage loads at transonic speeds is described. The geometric flexibility required for treatment of isolated and underwing, pylon-mounted stores is achieved by computing solutions on a five-level embedded grid arrangement. A completely automated grid generation procedure facilitates applications. Store modeling capability consists of bodies of revolution with multiple fore and aft fins. A body-conforming grid improves the accuracy of the computed store body flow field. A nonlinear finite difference relaxation scheme, developed specifically for modified transonic small disturbance flow equations, enhances numerical stability and accuracy. As a result, more accurate treatment of low aspect ratio, highly swept and tapered wing planforms is possible. A limited supersonic freestream capability is also provided. Pressure, load distribution, force and moment correlations show good agreement for several test cases.

Rosen, Bruce S.↗
NASA NTRS1989

Method to predict external store carriage loads at transonic speeds

A computational method for prediction of external store carriage loads at transonic speeds is described. The geometric flexibility required for treatment of isolated and underwing, pylon mounted stores is achieved by computing solutions on a five level embedded grid arrangement. A completely automated grid generation procedure facilitates applications. Store modeling capability consists of bodies of revolution with multiple fore and aft fins. A body conforming grid improves the accuracy of the computed store body flow field. A nonlinear finite difference relaxation scheme, developed specifically for modified transonic small disturbance flow equations, enhances numerical stability and accuracy. As a result, more accurate treatment of low aspect ratio, highly swept and tapered wing planforms is possible. A limited supersonic freestream capability is also provided. Pressure, load distribution, force and moment correlation show good agreement for several test cases.

Rosen, Bruce S.↗
DOE OSTI2024

Identifying Potential Sub-Synchronous Oscillations Using Impedance Scan Approach: Preprint

This paper presents an impedance scan study on the real-world power system in Australia that has observed 17 to 20 Hz intermittent subsynchronous oscillations. Through the impedance scan of each of the IBRs, both individually and collectively, potential resonance modes are identified. The impedance scans were carried out using electromagnetic transient PSCAD models of the network. The network comprises of site-specific, black-boxed models of IBRs supplied by the generators. The impedance scan approach was divided in three major steps: (1) the first step is to identify IBRs where impedance analysis needs to be performed based on the magnitude of oscillations observed at their points of interconnection (POIs); (2) the second step focuses on performing impedance scans at selected IBRs in single-machine infinite-bus (SMIB) configuration to identify internal resonance modes of an IBR and to evaluate if any of these modes become unstable under certain grid conditions; (3) the third step performs impedance scans at a few selected IBRs while connected to the wider network model to obtain the impedance response of both the IBR and the grid - this step evaluates control interactions among IBRs. The impedance scan study found that a few IBRs have an underdamped resonance mode at around 17 Hz, which becomes unstable under a certain operating condition. Another interesting finding was that certain IBRs increase the effective grid impedance seen by another IBR in proximity for a particular operating condition. Under this operating condition, the resulting subsynchronous oscillation mode is more pronounced due to a combination of control interactions among IBRs through the transmission network under study and a resonance mode inside the IBRs.

control interactions↗