Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid analysis”

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 523 records · Page 29

External-Compression Supersonic Inlet Design Code

A computer code named SUPIN has been developed to perform aerodynamic design and analysis of external-compression, supersonic inlets. The baseline set of inlets include axisymmetric pitot, two-dimensional single-duct, axisymmetric outward-turning, and two-dimensional bifurcated-duct inlets. The aerodynamic methods are based on low-fidelity analytical and numerical procedures. The geometric methods are based on planar geometry elements. SUPIN has three modes of operation: 1) generate the inlet geometry from a explicit set of geometry information, 2) size and design the inlet geometry and analyze the aerodynamic performance, and 3) compute the aerodynamic performance of a specified inlet geometry. The aerodynamic performance quantities includes inlet flow rates, total pressure recovery, and drag. The geometry output from SUPIN includes inlet dimensions, cross-sectional areas, coordinates of planar profiles, and surface grids suitable for input to grid generators for analysis by computational fluid dynamics (CFD) methods. The input data file for SUPIN and the output file from SUPIN are text (ASCII) files. The surface grid files are output as formatted Plot3D or stereolithography (STL) files. SUPIN executes in batch mode and is available as a Microsoft Windows executable and Fortran95 source code with a makefile for Linux.

Slater, John W.↗

Ultra-Large-Scale Power System Control and Coordination Architecture: A Strategic Framework for Integrating Advanced Grid Functionality

Electric power grids are being stressed by integration of intermittent renewable resources and significant adoption of distributed energy resources. The complexity of the grid is growing rapidly as we attempt to support technical, business, and societal goals for which power grids were not originally designed. Today, we largely take stability of the grid for granted. Stability could significantly degrade or collapse because of new dynamics introduced to the grid, and because the emerging extreme complexity makes traditional control analysis intractable, so that grid behavior is more unpredictable than in the past. This paper describes emerging issues in grid control and provides reasons why the present path of grid control evolution is problematic and presents an ultra-large scale framework for grid control that can solve today’s problems and those expected over the next 30 years. Failure to address these issues in the short term will result in rapidly escalating system deployment and maintenance costs, potential stranded assets related to replacement of the “ad hoc” systems, along with substantial operational risks that are unacceptable under current utility and regulatory practices.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Error analysis of finite element results on plates with nonuniform grids

Discrete Dirac delta functions in two dimensions and the governing plate differential equations are used to produce continuous approximations from discrete finite-element data on nonuniform grids. The resulting continuous approximation when compared to the discrete data provides an error analysis. Moreover, the continuous approximation can be differentiated to compute continuous stresses needed for component failure prediction. Stresses can then be calculated accurately at boundaries and at interfaces where finite-element analysis alone may not provide reliable results. The continuous solution along with the differential equations also provides a residual error which is used in the error analysis. This is a unique feature of this approach. Two application examples of the procedure are provided. One example involves the linear response of a clamped square plate under uniform transverse load. Results from finite-element analysis and the continuous solution from the present analysis are compared to a closed-form elasticity solution. The continuous solution is used in a Zienkiewicz-Zhu error estimator to assess the performance of the present analysis. The continuous solution is observed to dramatically correct the finite-element results. In a second example, results for the nonlinear response of a stiffened panel are compared to the finite-element solution.

Sistla, Rajaram↗

Gridded thermionic gun and integral superconducting ballistic bunch compression cavity

Electron-Ion colliders such as the Medium energy Electron Ion Collider (MEIC) being developed by JLAB require high current electrons with low energy spread for electron cooling of the collider ring. Accelerator techniques for improving bunch charge, average current, emittance, and energy spread are required for Energy Recovery Linacs (ERLs) and Circulator Rings (CR) for next generation colliders for nuclear physics experiments. Example candidates include thermionic-cathode electron guns with RF accelerating structures. Thermionic cathodes are known to produce high currents and have excellent lifetime. The success of the IR and THz Free-Electron Laser (FEL) designed and installed by Advanced Energy Systems at the Fritz Haber Institute (FHI) of the Max Planck Society in Berlin [1,2] demonstrates that gridded thermionic cathodes and rf systems be considered for next generation collider technology. In Phase 1 Advanced Energy Systems (AES) developed and analyzed a design concept using a superconducting cavity pair and gridded thermionic cathode. Analysis included Beam Dynamics and thermal analysis to show that a design of this type is feasible. The latest design goals for the MEIC electron cooler were for electron bunches of 420 pC at a frequency of 952.6 MHz with a magnetic field on the cathode of 2kG. This field magnetizes the beam imparting angular momentum that provides for helical motion of the electrons in the cooling solenoid. The helical motion increases the interaction time and improves the cooling efficiency. A coil positioned around the cathode providing 2kG field was developed. Beam dynamics simulations were run to develop the particle dynamics near the cathode and grid. Lloyd Young added capability to Tstep to include space charge effects between two plates and include image charge effects from the grid. He also added new pepper-pot geometry capability to account for honeycomb grids. These additions were used to develop the beam dynamics for this gun. The general design is a modified ballistic compression cavity pair with two independently powered cells [3]. The first is a cathode cell that includes the thermionic cathode and grid to provide for beam bunching. The second is a full cell with independent phasing and field levels designed to minimize energy spread. The primary goal for Phase II is to manufacture a superconducting gun with a thermionic cathode and imbedded coil. The system developed here is applicable to many high current electron accelerators. The analysis and design constraints imposed by the magnetized cathode make the cathode system developed here more complicated and limited than one without the magnetized beam constraints. High power ERLs would benefit by a gun with the capabilities shown here, 400 mA or more of current. ERLs hold great promise for electron cooling experiments, advanced light sources and Free Electron Lasers. This high current electron injector is a technological advance that will place the requirements for an ERL capable of providing quality bunches needed for cooling within the MEIC circulator ring within reach. This injector would have application to future ERLs around the world.

43 PARTICLE ACCELERATORS↗

Satellite-based Investigation of Power-Line Vegetation Encroachment in the US (SILVANUS)

Rapid wide-area assessments of vegetation encroachment on transmission and distribution line rights-of-way (ROW) is a highly desirable capability for understanding risks to the power grid during severe weather and wildfire events. Conventional assessments are time-consuming and expensive due to the need for in-situ inspections and the use of aerial assets. Performing conventional assessments on a wide area would require immense resources and time that might not be available within the horizon of an expected adverse event. Developing a capability to accurately assess vegetation encroachment into ROWs will enable faster analysis of potential grid vulnerabilities in NAERM. This project sought to develop a prototype capability for rapid ROW vegetation encroachment assessments by using Puerto Rico as a test case. Puerto Rico is a heavily forested island territory frequently impacted by tropical cyclones that threaten the electric grid by downing trees across transmission and distribution lines. Multispectral satellite imagery (MSI) enable very high resolution (i.e., 0.5 - 2 meter) assessment of vegetation conditions at scale and with revisit times appropriate for regular monitoring (e.g., weekly to quarterly, depending on cloud cover) of the entire island. Synthetic aperture radar (SAR) data from satellites was also investigated as solution to the cloud-cover issue as they are active sensors that emit and receive a microwave signal rather than relying on solar illumination, and are therefore unaffected by cloud cover and can collect data during day or night. Finally, MSI-derived digital surface models (DSMs) map the height of objects relative to sea level, and were assessed for their ability to estimate the height of forest canopies relative to coincident transmission lines. The results of the mapping investigation are described in this report.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Advanced Klystrons for High Efficiency Accelerator Systems (Final Report)

This program explored tailoring of RF pulses used to drive accelerator cavities. Simulations indicated that properly shaping the pulse risetime to match accelerator cavity characteristics reduced reflected power and increased total efficiency. Tailoring the pulse requires a high power, gridded, klystron to shape the risetime while also controlling the beam current. The Phase I program generated a preliminary design of a gridded electron gun for a klystron producing 5-10 MW of RF power. This required design of a segmented cathode using Controlled Porosity Reservoir cathodes to limit power deposition on the grid. The program was successful in computationally designing a gun producing a high quality electron beam with grid control. Additional analysis of pulse tailoring indicated that technique would only be useful for cavity drive pulses that were less than approximately 2-3 times the risetime. Otherwise, the efficiency gained during the risetime of the pulse became insignificant when considering the efficiency over the entire pulse. Consequently, it was determined that a Phase II program would not provide sufficient return to justify the cost. Never the less, other applications for a high power gridded gun are currently being pursued. This klystron, for example, would facilitate development inverse Comptom x-ray sources by providing a high repetition rate (10 -100 kHz) RF source.

43 PARTICLE ACCELERATORS↗

Phasor Measurement Units Optimal Placement and Performance Limits for Fault Localization

Here, the performance limits of faults localization are investigated using synchrophasor data. The focus is on a non-trivial operating regime where the number of Phasor Measurement Unit (PMU) sensors available is insufficient to have full observability of the grid state. Proposed analysis uses the Kullback Leibler (KL) divergence between the distributions corresponding to different fault location hypotheses associated with the observation model. This analysis shows that the most likely locations are concentrated in clusters of buses more tightly connected to the actual fault site akin to graph communities. Consequently, a PMU placement strategy is derived that achieves a near-optimal resolution for localizing faults for a given number of sensors. The problem is also analyzed from the perspective of sampling a graph signal, and how the placement of the PMUs i.e. the spatial sampling pattern and the topological characteristic of the grid affect the ability to successfully localize faults. To highlight the superior performance of presented fault localization and placement algorithms, the proposed strategy is applied to a modified IEEE 34, IEEE-123 bus test cases and to data from a real distribution grid. Additionally, the detection of cyber-physical attacks is also examined where PMU data and relevant Supervisory Control and Data Acquisition (SCADA) network traffic information are compared to determine if a network breach has affected the integrity of the system information and/or operations.

42 ENGINEERING↗

Large-Scale Data Challenges in Future Power Grids

This paper describes technical challenges in supporting large-scale real-time data analysis for future power grid systems and discusses various design options to address these challenges. Even though the existing U.S. power grid has served the nation remarkably well over the last 120 years, big changes are in the horizon. The widespread deployment of renewable generation, smart grid controls, energy storage, plug-in hybrids, and new conducting materials will require fundamental changes in the operational concepts and principal components. The whole system becomes highly dynamic and needs constant adjustments based on real time data. Even though millions of sensors such as phase measurement units (PMUs) and smart meters are being widely deployed, a data layer that can support this amount of data in real time is needed. Unlike the data fabric in cloud services, the data layer for smart grids must address some unique challenges. This layer must be scalable to support millions of sensors and a large number of diverse applications and still provide real time guarantees. Moreover, the system needs to be highly reliable and highly secure because the power grid is a critical piece of infrastructure. No existing systems can satisfy all the requirements at the same time. We examine various design options. In particular, we explore the special characteristics of power grid data to meet both scalability and quality of service requirements. Our initial prototype can improve performance by orders of magnitude over existing general-purpose systems. The prototype was demonstrated with several use cases from PNNL’s FPGI and was shown to be able to integrate huge amount of data from a large number of sensors and a diverse set of applications.

Smart Grid, Real Time, Commercial Products, Big Da↗

Updates to Implicit Edge-Based Gradient Methods

In this paper, we report updates to the implicit edge-based gradient methods originally introduced in [H. Nishikawa, AIAA Paper 2020-3048, 2020]. First, we clarify the relationship between gradient accuracy and truncation error and show that the quadratic method involves a free parameter. Then, we provide a complete description and a simplified matrix form of the implicit gradient systems including a consistent boundary treatment, and derive a set of parameters for achieving fourth-order gradient accuracy on regular tetrahedral grids. A stability analysis is performed for a relaxation scheme used to solve the implicit gradient systems, and the result serves as a guide for choosing parameters. Numerical results are shown for accuracy verification and also for realistic inviscid flow problems in three dimensions, including flows with shock waves.

Weighted Least-Squares↗

Region-Based Stability Analysis of Resilient Distribution Systems with Hybrid Grid-Forming and Grid-Following Inverters

In this paper, distribution system stability considering mixed grid-forming (GFM) and grid-following (GFL) inverters is discussed. A holistic small-signal model of the distribution feeder accommodating both GFM and GFL inverters is derived. In contrast to conventional point-by-point stability analysis, a region-based stability assessment approach is proposed and implemented to evaluate small-signal system stability with both GFM and GFL controls. Furthermore, to enhance the computational efficiency, an Artificial Intelligence (AI) assisted Kernel Ridge Regression (KRR) approach is developed to obtain the stability region boundary. Therefore, with the derived stability region, the impact of critical parameters on system stability is intuitively shown, and more importantly, the relationship between the maximum penetration level and the corresponding feasible range of the selected control parameters can be quantified. The proposed approach has been validated using MATLAB/Simulink.

grid-following↗

Sample IEEE123 Bus system for OEDI SI

Time series load and PV data from an IEEE123 bus system. An example electrical system, named the OEDI SI feeder, is used to test the workflow in a co-simulation. The system used is the IEEE123 test system, which is a well studied test system (see link below to IEEE PES Test Feeder), but some modifications were made to it to add some solar power modules and measurements on the system. The aim of this project is to create an easy-to-use platform where various types of analytics can be performed on a wide range of electrical grid datasets. The aim is to establish an open-source library of algorithms that universities, national labs and other developers can contribute to which can be used on both open-source and proprietary grid data to improve the analysis of electrical distribution systems for the grid modeling community. OEDI Systems Integration (SI) is a grid algorithms and data analytics API created to standardize how data is sent between different modules that are run as part of a co-simulation. The readme file included in the S3 bucket provides information about the directory structure and how to use the algorithms. The sensors.json file is used to define the measurement locations.

123 bus↗

Numerical prediction of the Mid-Atlantic states cyclone of 18-19 February 1979

A series of forecast experiments was conducted to assess the accuracy of the GLAS model, and to determine the importance of large scale dynamical processes and diabatic heating to the cyclogenesis. The GLAS model correctly predicted intense coastal cyclogenesis and heavy precipitation. Repeated without surface heat and moisture fluxes, the model failed to predict any cyclone development. An extended range forecast, a forecast from the NMC analysis interpolated to the GLAS grid, and a forecast from the GLAS analysis with the surface moisture flux excluded predicted weak coastal low development. Diabatic heating resulting from oceanic fluxes significantly contributed to the generation of low level cyclonic vorticity and the intensification and slow rate of movement of an upper level ridge over the western Atlantic. As an upper level short wave trough approached this ridge, diabatic heating associated with the release of latent heat intensified, and the gradient of vorticity, vorticity advection and upper level divergence in advance of the trough were greatly increased, providing strong large scale forcing for the surface cyclogenesis.

Atlas, R.↗

Wind-US Users Guide Version 3.0

Wind-US is a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Currently, the code supports the solution of the Euler and Navier-Stokes equations of fluid mechanics, along with supporting equation sets governing turbulent and chemically reacting flows. Wind-US is a product of the NPARC Alliance, a partnership between the NASA Glenn Research Center (GRC) and the Arnold Engineering Development Complex (AEDC) dedicated to the establishment of a national, applications-oriented flow simulation capability. The Boeing Company has also been closely associated with the Alliance since its inception, and represents the interests of the NPARC User's Association. The "Wind-US User's Guide" describes the operation and use of Wind-US, including: a basic tutorial; the physical and numerical models that are used; the boundary conditions; monitoring convergence; the files that are read and/or written; parallel execution; and a complete list of input keywords and test options. For current information about Wind-US and the NPARC Alliance, please see the Wind-US home page at http://www.grc.nasa.gov/WWW/winddocs/ and the NPARC Alliance home page at http://www.grc.nasa.gov/WWW/wind/. This manual describes the operation and use of Wind-US, a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Wind-US represents a merger of the capabilities of four CFD codes - NASTD (a structured grid flow solver developed at McDonnell Douglas, now part of Boeing), NPARC (the original NPARC Alliance structured grid flow solver), NXAIR (an AEDC structured grid code used primarily for store separation analysis), and ICAT (an unstructured grid flow solver developed at the Rockwell Science Center and Boeing).

Navier-Stokes equation↗

Numerical simulation of steady and unsteady viscous flow in turbomachinery using pressure based algorithm

The objective of this research is to simulate steady and unsteady viscous flows, including rotor/stator interaction and tip clearance effects in turbomachinery. The numerical formulation for steady flow developed here includes an efficient grid generation scheme, particularly suited to computational grids for the analysis of turbulent turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and convergency in the solution. The numerical model is used to investigate the structure of tip clearance flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to-blade variation of all three velocity components, pitch and yaw angles, losses and blade static pressures in the tip clearance region. The simulation also includes evaluation of such quantities of leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial dissipation, that the tip clearance flow field can be captured accurately. The above numerical formulation was modified to incorporate time accurate solutions. An inner loop iteration scheme is used at each time step to account for the non-linear effects. The computation of unsteady flow through a flat plate cascade subjected to a transverse gust reveals that the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the computation upstream of the stator, and the upstream rotor wake is specified from the experimental data. The results show that the stator potential effects have appreciable influence on the upstream rotor wake. The predicted unsteady wake profiles are compared with the available experimental data and the agreement is good. The numerical results are interpreted to draw conclusions on the unsteady wake transport mechanism in the blade passage.

Lakshminarayana, B.↗

Analysis of Distributed Energy Storage as a Core Grid Infrastructure via Production Cost Modeling

Energy storage plays a pivotal role in enabling power system operation with more flexibility and resilience. Unlike current practice that considers energy storages as attached ancillary devices, this paper focuses on storages as a core infrastructure by looking at their spatial distribution in the system. A sensitivity analysis based on production cost modeling is conducted to demonstrate the benefits of distributed energy storages. First, the modeling of energy storages in production cost modeling is presented. Second, potential optimal locations of distributed energy storages in a power system are discussed. Finally, multiple scenarios with various numbers and locations of additional distributed energy storages in the WECC 2030 model are created. The production cost modeling results of these scenarios show that distributed energy storages have higher utilization compared to the centralized ES units and therefore provide significantly more benefits in terms of reduction in generation cost, emission cost, and volatility of location marginal prices. A saturation effect is observed suggesting the selection of optimal locations will further improve the benefits.

Nguyen, Quan H.↗

AMPS: An Augmented Matrix Formulation for Principal Submatrix Updates with Application to Power Grids

We present an augmented matrix approach to update the solution to a linear system of equations when the coefficient matrix is modified by a few elements within a principal submatrix. This problem arises in the dynamic security analysis of a power grid, where operators need to perform $N-x$ contingency analysis, i.e., determine the state of the system when up to $x$ links from $N$ fail. Our algorithms augment the coefficient matrix to account for the changes in it, and then compute the solution to the augmented system without refactoring the modified matrix. We provide two algorithms, a direct method, and a hybrid direct-iterative method for solving the augmented system. We also exploit the sparsity of the matrices and vectors to accelerate the overall computation. Our algorithms are compared on three power grids with PARDISO, a parallel direct solver, and CHOLMOD, a direct solver with the ability to modify the Cholesky factors of the coefficient matrix. We show that our augmented algorithms outperform PARDISO (by two orders of magnitude), and CHOLMOD (by a factor of up to 5). Further, our algorithms scale better than CHOLMOD as the number of elements updated increases. The solutions are computed with high accuracy. Our algorithms are capable of computing $N-x$ contingency analysis on a $778K$ bus grid, updating a solution with $x=20$ elements in $1.6 \times 10^{-2}$ seconds on an Intel Xeon processor.

Yeung, Yu-Hong↗