Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

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

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

The Use of Embedded Electric Grid Storage for Resilience, Operational Flexibility, and Cyber-Security

In recent years, bulk energy storage has been applied to electric power systems as an auxiliary device for the support of grid reliability via grid services. This approach is useful but only extracts value from storage on a marginal basis because grid services involve only a tiny fraction of the power flowing in a grid. While storage is flexible enough to perform many services, assessing its value on a stacked-services basis obscures its real value: storage is a buffer for electric energy flows. All of its capabilities stem from this one fundamental property. Limiting storage to use as an ancillary services device leaves this primary value untapped. The real value of storage is as a means to provide a key characteristic missing from power grids: the ability to absorb stresses with little or no loss of performance – the essence of resilience. Storage applied systematically throughout the grid can provide the missing “shock absorber” springiness that the grid is missing. To provide this value, storage must be incorporated into the grid as core infrastructure and must be deeply integrated into grid operations. Doing so will provide far-reaching benefits to users of electricity at all levels, including vastly increased system resilience, expanded system operational flexibility, support for critical lifeline functions during critical events, and even improved cyber security.

25 ENERGY STORAGE↗
DOE OSTI2021

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

AC and DC Hybrid Distribution Grids with Solar Integration: Architecture, Stabilization and Cost Assessment

In this project, a holistic analysis of architecture, stabilization, and cost/efficiency analysis in hybrid AC and DC distribution grids are conducted. Particularly, versatile and cost-efficient multiport converters are developed to not only integrate solar sources and other DERs in DC grids, but also facilitate the interactive operation of DC sub-grids and conventional AC distribution grids. Meanwhile, a universal and extended impedance-based stabilization approach with a decentralized and adaptive virtual impedance loop is developed in hybrid AC and DC distribution grids, which comprehensively covers the active stabilization of DC sections, AC sections, and interface inverters interlinking both AC and DC sections. Furthermore, to quantitatively evaluate the cost and efficiency of hybrid AC and DC distribution grids and quantify the improvement of hybrid AC and DC grids over conventional pure AC grids, the project team develops an alpha-version tool to monitor and calculate the efficiency and cost of the DC, AC, or hybrid AC and DC grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Impedance-Based Stability Analysis in Grid Interconnection Impact Study Owing to the Increased Adoption of Converter-Interfaced Generators

This study investigates the emerging harmonic stability concerns to be addressed by grid planners in generation interconnection studies, owing to the increased adoption of renewable energy resources connected to the grid via power electronic converters. The wideband and high-frequency electromagnetic transient (EMT) characteristics of these converter-interfaced generators (CIGs) and their interaction with the grid impedance are not accurately captured in the typical dynamic studies conducted by grid planners. This paper thus identifies the desired components to be studied and subsequently develops a practical process for integrating a new CIG into a grid with the existing CIGs. The steps of this process are as follows: the impedance equation of a CIG using its control dynamics and an interface filter to the grid, for example, an LCL filter (inductor-capacitor-inductor type), is developed; an equivalent impedance model including the existing CIGs nearby and the grid observed from the point of common coupling are derived; the system stability for credible operating scenarios is assessed. Detailed EMT simulations validate the accuracy of the impedance models and stability assessment for various connection scenarios. Here, by complementing the conventional EMT simulation studies, the proposed analytical approach enables grid planners to identify critical design parameters for seamlessly integrating a new CIG and ensuring the reliability of the grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

SymPix: A spherical grid for efficient sampling of rotationally invariant operators

We present SymPix, a special-purpose spherical grid optimized for efficiently sampling rotationally invariant linear operators. This grid is conceptually similar to the Gauss–Legendre (GL) grid, aligning sample points with iso-latitude rings located on Legendre polynomial zeros. Unlike the GL grid, however, the number of grid points per ring varies as a function of latitude, avoiding expensive oversampling near the poles and ensuring nearly equal sky area per grid point. The ratio between the number of grid points in two neighboring rings is required to be a low-order rational number (3, 2, 1, 4/3, 5/4, or 6/5) to maintain a high degree of symmetries. Our main motivation for this grid is to solve linear systems using multi-grid methods, and to construct efficient preconditioners through pixel-space sampling of the linear operator in question. As a benchmark and representative example, we compute a preconditioner for a linear system that involves the operator D-hat + B-hat {sup T}N{sup −1} B-hat , where B-hat and D-hat may be described as both local and rotationally invariant operators, and N is diagonal in the pixel domain. For a bandwidth limit of ℓ{sub max} = 3000, we find that our new SymPix implementation yields average speed-ups of 360 and 23 for B-hat {sup T}N{sup −1} B-hat and D-hat , respectively, compared with the previous state-of-the-art implementation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

An investigation of the sub-grid variability of trace gases and aerosols for global climate modeling

One fundamental property and limitation of grid based models is their inability to identify spatial details smaller than the grid cell size. While decades of work have gone into developing sub-grid treatments for clouds and land surface processes in climate models, the quantitative understanding of sub-grid processes and variability for aerosols and their precursors is much poorer. In this study, WRF-Chem is used to simulate the trace gases and aerosols over central Mexico during the 2006 MILAGRO field campaign, with multiple spatial resolutions and emission/terrain scenarios. Our analysis focuses on quantifying the sub-grid variability (SGV) of trace gases and aerosols within a typical global climate model grid cell, i.e. 75x75 km2. Our results suggest that a simulation with 3-km horizontal grid spacing adequately reproduces the overall transport and mixing of trace gases and aerosols downwind of Mexico City, while 75-km horizontal grid spacing is insufficient to represent local emission and terrain-induced flows along the mountain ridge, subsequently affecting the transport and mixing of plumes from nearby sources. Therefore, the coarse model grid cell average may not correctly represent aerosol properties measured over polluted areas. Probability density functions (PDFs) for trace gases and aerosols show that secondary trace gases and aerosols, such as O3, sulfate, ammonium, and nitrate, are more likely to have a relatively uniform probability distribution (i.e. smaller SGV) over a narrow range of concentration values. Mostly inert and long-lived trace gases and aerosols, such as CO and BC, are more likely to have broad and skewed distributions (i.e. larger SGV) over polluted regions. Over remote areas, all trace gases and aerosols are more uniformly distributed compared to polluted areas. Both CO and O3 SGV vertical profiles are nearly constant within the PBL during daytime, indicating that trace gases are very efficiently transported and mixed vertically by turbulence. But, simulated horizontal variability indicates that trace gases and aerosols are not well mixed horizontally in the PBL. During nighttime the SGV for trace gases is maximum at the surface, and quickly decreases with height. Unlike the trace gases, the SGV of BC and secondary aerosols reaches a maximum at the PBL top during the day. The SGV decreases with distance away from the polluted urban area, has a more rapid decrease for long-lived trace gases and aerosols than for secondary ones, and is greater during daytime than nighttime. The SGV of trace gases and aerosols is generally larger than for meteorological quantities. Emissions can account for up to 50% of the SGV over urban areas such as Mexico City during daytime for less-reactive trace gases and aerosols, such as CO and BC. The impact of emission spatial variability on SGV decays with altitude in the PBL and is insignificant in the free troposphere. The emission variability affects SGV more significantly during daytime (rather than nighttime) and over urban (rather than rural or remote) areas. The terrain, through its impact on meteorological fields such as wind and the PBL structure, affects dispersion and transport of trace gases and aerosols and their SGV.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A parameterization of sub-grid topographical effects on solar radiation in the E3SM Land Model (version 1.0): implementation and evaluation over the Tibetan Plateau

Abstract. Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth system models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes that the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to compare with the ELM simulations. The results show that topography has non-negligible effects on surface energy budget, snow cover, snow depth, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20, 10, and 5 W m−2, respectively. The differences in land surface albedo, snow cover fraction, snow depth, and surface temperature between TOP and PP exceed 0.1, 0.1, 10 cm, and 1 K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are larger with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover, and surface temperature, especially in the high-elevation and snow-covered regions over the TP. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

Hao, Dalei (ORCID:0000000334979774)↗
DOE OSTI2021

A Parameterization of Sub-grid Topographical Effects on Solar Radiation in the E3SM Land Model (Version 1.0): Implementation and Evaluation Over the Tibetan Plateau

Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth System Models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to evaluate the performance of ELM. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover and surface temperature especially in the high-elevation and snow-cover regions over the TP. Topography has non-negligible effects on surface energy budget, snow cover, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20 W/m 2 , 10 W/m 2 , and 5 W/m 2 , respectively. The differences in land surface albedo, snow cover fraction, and surface temperature between TOP and PP exceed 0.1, 20%, and 1K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations, and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are more significant with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

58 GEOSCIENCES↗
NASA NTRS1980

Grid generation for general three-dimensional configurations

A suitable grid system for complex three dimensional configurations such as wing/body/nacelle shape for the solution of nonlinear transonic flow problems was constructed. Two approaches were explored based on Thompson's body fitted coordinate concept. The most general approach was to divide the computational domain into multiple rectangular blocks where the configuration itself was also represented by a set of blocks whose structure follows the natural lines of the configuration. The block structured grid system was adaptable to complex configurations and gives good grid quality near physical corners. However, it introduced algorithm issues for the flow solution concerning the treatment of nonanalytic grid block boundaries and nonstandard grid cells. These issues were explored in relation to the grid generation. A more limited approach treats a wing/body configuration with only a single rectangular block in computational space. The issues involving nonstandard cells were avoided, but other limitations on grid resolution appear. Both a linear and a nonlinear system of grid generation equations were developed including methods of grid control.

Lee, K. D.↗
NASA NTRS1986

Algebraic grid generation about wing-fuselage bodies

An algebraic procedure for the generation of boundary-fitted grids about wing-fuselage configurations is presented. A wing-fuselage configuration is specified by cross sections and mathematically represented by Coons' patches. A configuration is divided into sections so that several grid blocks that either adjoin each other or partially overlap each other can be generated. Each grid has six exterior surfaces that map into a computational cube. Grids are first determined on the six boundary surfaces and then in the interior. Grid curves that are on the surface of the configuration are derived from the intersection of planes with the Coons' patch definition. Single-valued functions relating approximate arc lengths along the grid curves to a computational coordinate define the distribution of grid points. The two-boundary technique and transfinite interpolation are used to determine the boundary surface grids that are not on the configuration, and transfinite interpolation with linear blending functions is used to determine the interior grid.

Smith, R. E.↗
NASA NTRS1989

A Simple Algebraic Grid Adaptation Scheme with Applications to Two- and Three-dimensional Flow Problems

An algebraic adaptive grid scheme based on the concept of arc equidistribution is presented. The scheme locally adjusts the grid density based on gradients of selected flow variables from either finite difference or finite volume calculations. A user-prescribed grid stretching can be specified such that control of the grid spacing can be maintained in areas of known flowfield behavior. For example, the grid can be clustered near a wall for boundary layer resolution and made coarse near the outer boundary of an external flow. A grid smoothing technique is incorporated into the adaptive grid routine, which is found to be more robust and efficient than the weight function filtering technique employed by other researchers. Since the present algebraic scheme requires no iteration or solution of differential equations, the computer time needed for grid adaptation is trivial, making the scheme useful for three-dimensional flow problems. Applications to two- and three-dimensional flow problems show that a considerable improvement in flowfield resolution can be achieved by using the proposed adaptive grid scheme. Although the scheme was developed with steady flow in mind, it is a good candidate for unsteady flow computations because of its efficiency.

Hsu, Andrew T.↗
NASA NTRS1990

An overlapped grid method for multigrid, finite volume/difference flow solvers: MaGGiE

The objective is to develop a domain decomposition method via overlapping/embedding the component grids, which is to be used by upwind, multi-grid, finite volume solution algorithms. A computer code, given the name MaGGiE (Multi-Geometry Grid Embedder) is developed to meet this objective. MaGGiE takes independently generated component grids as input, and automatically constructs the composite mesh and interpolation data, which can be used by the finite volume solution methods with or without multigrid convergence acceleration. Six demonstrative examples showing various aspects of the overlap technique are presented and discussed. These cases are used for developing the procedure for overlapping grids of different topologies, and to evaluate the grid connection and interpolation data for finite volume calculations on a composite mesh. Time fluxes are transferred between mesh interfaces using a trilinear interpolation procedure. Conservation losses are minimal at the interfaces using this method. The multi-grid solution algorithm, using the coaser grid connections, improves the convergence time history as compared to the solution on composite mesh without multi-gridding.

Baysal, Oktay↗
NASA NTRS1990

A solution adaptive grid procedure for an upwind parabolized flow solver

A solution adaptive grid procedure based on an error equi-distribution scheme is developed and applied to a Parabolized Navier-Stokes solver. An improved method for selecting weighting functions is introduced which involves normalizing a combination of flowfield gradients and curvature of a number of dependent variables and then selecting the largest at each point. The scheme re-distributes grid points line-by-line, with grid point motion controlled by forces analogous to tensional and torsional spring forces with the spring constants set equal to the weighting functions. Torsional terms are functions of the grid point positions along neighboring grid lines and provide grid smoothness and stability. A grid-fitting scheme is introduced for external flows in which the number of grid points in the freestream are reduced to a minimum. Results for several problems are presented to demonstrate the improvements obtainable with the solution adaptive grid procedure.

Harvey, Albert D.↗
NASA NTRS1992

TCGRID: A three dimensional C-grid generator for turbomachinery

A fast 3-D grid code for turbomachinery was developed. The code, TCGRID (Turbomachinery C-GRID), can generate either C or H type grids for fairly arbitrary axial or radical turbomachinery geometries. The code also has limited blocked grid capability and can generate an axisymmetric H type grid upstream of the blade row or an O type grid within the tip clearance region. Hub and tip geometries are input as a simple list of pairs. All geometric data is handled using parametric splines so that geometries that turn 90 degrees can be handled without difficulty. Blade input is in standard MERIDL or Lewis compressor design code format. TCGRID adds leading and trailing edge circles to MERIDL geometries and intersects the blade with the hub and tip if necessary using a novel intersection algorithm. The procedure used to generate the grid is given. Output is in PLOT3D format, which can also be read by the RVC3D (Rotor Viscous Code 3-D) Navier-Stokes code for turbomachinery. Intermediate 2-D or 3-D grids useful for debug and other purposes can also be output using a convenient output flag. A grid generated figure is given.

Chima, Rodrick V.↗
NASA NTRS1992

Cartesian based grid generation/adaptive mesh refinement

Grid adaptation has recently received attention in the computational fluid dynamics (CFD) community as a means to capture the salient features of a flowfield by either moving grid points of a structured or by adding cells in an unstructured manner. An approach based on a background cartesian mesh is investigated from which the geometry is 'cut' out of the mesh. Once the mesh is obtained, a solution on this coarse grid is found, that indicates which cells need to be refined. This process of refining/solving continues until the flow is grid refined in terms of a user specified global parameter (such as drag coefficient etc.). The advantages of this approach are twofold: the generation of the base grid is independent of the topology of the bodies or surfaces around/through which the flow is to be computed, and the resulting grid (in uncut regions) is highly isotropic, so that the truncation error is low. The flow solver (which, along with the grid generation is still under development) uses a completely unstructured data base, and is a finite volume, upwinding scheme. Current and future work will address generating Navier-Stokes suitable grids by using locally aligned and normal face/cell refining. The attached plot shows a simple grid about two turbine blades.

Coirier, William J.↗
NASA NTRS1992

Interactive solution-adaptive grid generation procedure

TURBO-AD is a two dimensional interactive solution adaptive grid generation program. The code uniquely combines a grid adaptation technique that uses parametric mapping with control sources and an algebraic grid generation that uses control points into a single software package. The grid adaptation is achieved by first adapting the control points to a numerical solution in the parametric domain using the control sources obtained from flow properties. Then a new grid is generated from the adapted control net using the control point formulation. The new adapted grid in the parametric domain is then mapped back to the physical domain. This solution adaptive grid generation process is efficient because the number of control points is much less than the number of grid points and the grid generation from the adapted control net is an efficient algebraic process.

Choo, Yung K.↗
NASA NTRS1993

A solution-adaptive hybrid-grid method for the unsteady analysis of turbomachinery

A solution-adaptive method for the time-accurate analysis of two-dimensional flows in turbomachinery is described. The method employs a hybrid structured-unstructured zonal grid topology in conjunction with appropriate modeling equations and solution techniques in each zone. The viscous flow region in the immediate vicinity of the airfoils is resolved on structured O-type grids while the rest of the domain is discretized using an unstructured mesh of triangular cells. Implicit, third-order accurate, upwind solutions of the Navier-Stokes equations are obtained in the inner regions. In the outer regions, the Euler equations are solved using an explicit upwind scheme that incorporates a second-order reconstruction procedure. An efficient and robust grid adaptation strategy, including both grid refinement and coarsening capabilities, is developed for the unstructured grid regions. Grid adaptation is also employed to facilitate information transfer at the interfaces between unstructured grids in relative motion. Results for grid adaptation to various features pertinent to turbomachinery flows are presented. Good comparisons between the present results and experimental measurements and earlier structured-grid results are obtained.

Mathur, Sanjay R.↗
NASA NTRS1993

Unstructured grid large eddy simulation of wall bounded turbulent flows

Historically, large eddy simulations (LES) have been restricted to simple geometries where spectral or finite difference methods have dominated due to their efficient use of structured grids. Structured grids, however, not only have difficulty representing complex domains and adapting to complicated flow features, but also are rather inefficient for simulating flows at high Reynolds numbers. The lack of efficiency stems from the need to resolve the viscous sub layer which requires very fine resolution in all three directions near the wall. Structured grids make use of a stretching to reduce the normal grid spacing but must carry the fine resolution in the streamwise and spanwise directions throughout the domain. The unnecessarily fine grid for much of the domain leads to disturbingly high grid estimates. Chapman (1979), and later Moin & Jimenez (1993), pointed out that, in order to advance the technology to airfoils at flight Reynolds numbers, structured grids must be abandoned in lieu of what are known as nested or unstructured grids. We illustrate the ability of an unstructured mesh to refine only the near wall region. Note the large number of points near the wall (where the fine vortical features need better resolution) and the coarseness in all directions away from the wall (where the scales are much larger). The important difference between this approach and the usual structured grid stretching is that the number of elements used to discretize the spanwise and streamwise features of the flow is reduced in each successive layer coming off the wall. This is due to the fact that the elements not only grow in the normal direction, but in the other directions as well. This greatly reduces the total number of points or elements required for a given Reynolds number flow.

Jansen, Kenneth↗