Search NASA⌕ Search

SEARCH · Search NASA

Results for “Regular Grid”

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 73 records · Page 4

Kriging as a Means of Improving WAAS Availability

The Wide Area Augmentation System (WAAS), an augmentation of the Global Positioning System (GPS), provides safe and reliable use of GPS signals for airline navigation over much of North America. Currently the largest source of positioning error in the system is signal delay caused by the ionosphere. To allow the user to take account of such error, WAAS computes and broadcasts ionospheric vertical delays at a set of regularly-spaced grid points. In addition, WAAS computes and broadcasts a safety-critical integrity bound at each ionospheric grid point (IGP) called the Grid Ionospheric Vertical Error (GIVE). GIVEs are constructed to be sufficiently large to protect the user against positioning error due to the presence of ionospheric irregularity. In the initial operating capability (IOC) of WAAS, the vertical delay estimate at each IGP is determined from a planar fit of neighboring slant delay measurements, projected to vertical using an obliquity factor specified by the standard thin-shell model of the ionosphere. In WAAS Follow-On (WFO) Release 3, however, the vertical delay will be estimated by an established, geo-statistical technique known as kriging. Compared to the planar fit model, the kriging model is found, in general, to match better the observed random structure of the vertical delay. This paper presents the kriging methodology that will be used to estimate the vertical delay and its uncertainty at each IGP, and it assesses the subsequent improvement in WAAS availability enabled by kriging.

GPS↗

Three dimensional imaging utilizing structured light

This paper describes a method of remote sensing 3-dimensional structure of the proximity utilizing a laser, a holographic grating, and a single regular CCD camera. Basically, the laser beam is split by a holographic grating to form a regular spaced grid of laser beams that are projected into the field of view of a CCD camera.

Chang, Johnny↗

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↗

Machine-Learning-Based Adaptive Thinning of CrIS Radiances to Improve Global Tropical Cyclone Analysis and Forecasts

This work is focused on optimizing the assimilation of hyperspectral infrared (IR) radiances from the Cross-track Infrared Sounder (CrIS) with the goal of improving the representation of tropical cyclones (TCs) in global analyses and forecasts. Current operational assimilation systems rely on subsampling IR radiances on a regular thinning grid. A new and improved adaptive methodology based on machine learning (ML) recognizes TCs from geostationary satellite imagery and is implemented in the Goddard Earth Observing System (GEOS) model and data assimilation framework. The ML methodology is extensively trained on existing TC data sets and creates for each TC a dynamic mask, based on the evolving shape and life cycle of that specific event. Once a TC mask is created, a switch is then activated in the data assimilation system to alter the thinning, ingesting more CrIS radiances within the moving mask, thus increasing the TC sampling. After the TC dissipates, the assimilation of CrIS radiances reverts to normal data density. Results of TC segmentation provided by a state-of-the-art generative machine learning model known as the Denoising Diffusion Probabilistic Model (DDPM) are compared to the previously used U-Net model. The new approach surpasses the performance of the previously developed one. The methodology is applied to both clear-sky and cloud-cleared radiances. Benefits from the latter methodology, particularly in improving the structure of TCs and the intensity forecasts, are presented.

Oreste Reale↗

Parallel architectures for iterative methods on adaptive, block structured grids

A parallel computer architecture well suited to the solution of partial differential equations in complicated geometries is proposed. Algorithms for partial differential equations contain a great deal of parallelism. But this parallelism can be difficult to exploit, particularly on complex problems. One approach to extraction of this parallelism is the use of special purpose architectures tuned to a given problem class. The architecture proposed here is tuned to boundary value problems on complex domains. An adaptive elliptic algorithm which maps effectively onto the proposed architecture is considered in detail. Two levels of parallelism are exploited by the proposed architecture. First, by making use of the freedom one has in grid generation, one can construct grids which are locally regular, permitting a one to one mapping of grids to systolic style processor arrays, at least over small regions. All local parallelism can be extracted by this approach. Second, though there may be a regular global structure to the grids constructed, there will be parallelism at this level. One approach to finding and exploiting this parallelism is to use an architecture having a number of processor clusters connected by a switching network. The use of such a network creates a highly flexible architecture which automatically configures to the problem being solved.

Gannon, D.↗

Wavenumber-frequency Spectra of Pressure Fluctuations on a Generic Space Vehicle Measured via Unsteady Pressure-Sensitive Paint

Time histories of pressure fluctuations on a generic, hammerhead space vehicle model were measured using unsteady Pressure-Sensitive Paint (uPSP). The test was conducted in the 11-foot transonic wind tunnel of NASA Ames Research Center over a Mach number range of 0.6 M 1.2, and angles of attack of -4 4. The model was coated with a porous binder and PtTFPP-based porous polymer paint. An elaborate system of four high-speed cameras, and forty LED lamps was used for image acquisition. Various steps for image registration, reduction of shot noise, photogrammetry procedure to map images from the four cameras on a grid for the model, and finally a calibration procedure to convert the measured fluctuations in light intensity to fluctuating pressure, are discussed in the paper. The calibration process using a set of unsteady pressure sensors mounted on the model, was found to overcome some of the inherent problems of the fast response paint, such as rapid photo-degradation, non-linearity in pressure response, and significant temperature sensitivity. Comparison of spectra of pressure fluctuations between UPSP and pressure sensors demonstrated the ability of the paint to faithfully follow fluctuations up to 10 kHz, the maximum attempted. It was also found that the camera bit-depth and the illumination level limited the lowest measurable levels of pressure fluctuations to around 140dB. The large data set exposed various critical transonic flow physics not seen before, such as a coupling of the shock motion on the Payload Fairing (PF) with the separated flow region on the upper stage of the launch vehicle, and upstream convection of pressure fluctuation on PF at certain Mach numbers. The data also confirmed the expectation of a general lowering of the coefficient of pressure fluctuation with Mach number. The availability of the data set on a dense, regularly-spaced, surface grid allowed for the calculation of wavenumber-frequency (k-) spectra via straightforward applications of Fourier transform. The k- spectra were compared for the separated flow regions on the Second Stage, and the shock-boundary layer interactions on PF. The former showed self-similarity with Mach number while the latter was distinctly different, and confirmed the upstream propagation of pressure fluctuations. The k- spectra were dominated by the convected fluctuations; the acoustic domain was not discernable. These data, valuable for the vibro-acoustics analysis of aerospace vehicles, are believed to be the first obtained for the transonic flight regime, and pave the path for application on production models of aerospace vehicles.

buffet↗

Data traffic reduction schemes for sparse Cholesky factorizations

Load distribution schemes are presented which minimize the total data traffic in the Cholesky factorization of dense and sparse, symmetric, positive definite matrices on multiprocessor systems with local and shared memory. The total data traffic in factoring an n x n sparse, symmetric, positive definite matrix representing an n-vertex regular 2-D grid graph using n (sup alpha), alpha is equal to or less than 1, processors are shown to be O(n(sup 1 + alpha/2)). It is O(n(sup 3/2)), when n (sup alpha), alpha is equal to or greater than 1, processors are used. Under the conditions of uniform load distribution, these results are shown to be asymptotically optimal. The schemes allow efficient use of up to O(n) processors before the total data traffic reaches the maximum value of O(n(sup 3/2)). The partitioning employed within the scheme, allows a better utilization of the data accessed from shared memory than those of previously published methods.

Naik, Vijay K.↗

Impact of Assimilating AIRS Cloud Cleared Radiances on the Analysis and Forecast of Polar Lows

Previous work by this team (Reale et al. 2018) has found that the current assimilation of AIRS (Atmospheric InfraRed Sounder) radiances on a regularly spaced thinning grid is suboptimal, probably because of horizontal error correlation over meteorologically inactive areas. Moreover, cloud-cleared radiances appear to be a better product than clear-sky radiances, but need to be assimilated at a much lower density globally, because of the higher information content. Specifically: 1. Assimilation of AIRS cloud-cleared radiances at a density of about one quarter of the clear-sky radiances improves global forecast skill; 2. An adaptive thinning strategy assimilating cloud-cleared radiances at reduced density globally except around tropical cyclones (TCs), leads to substantial improvement in the structure and intensity forecast of TCs without damaging global skill.

Polar Low↗

Comparing Level-2 and Level-3 Satellite Ocean Color Retrieval Validation Methodologies

Many ocean color data applications leverage global spatially composited Level-3(L3) satellite data, because of their regular Earth-grid frame of reference. However, ocean color satellite retrieval performance is routinely evaluated on Level-2 (L2) data at the native satellite swath resolution and geometries. This study assesses how accurately binned and gridded L3 data represent L2 satellite data products via satellite-to-in situ match-up activities.L2 and L3 satellite data retrievals of the photosynthetic pigment chlorophyll-a are compared with a common in situ data set, revealing similar L2 and L3 satellite-to-in situ performance for both MODIS-Aqua and VIIRS-SNPP. This agreement implies that L2 validation results are generally applicable to L3 data. However, uncertainties are introduced during the generation of L3 data from L2 data. L3 data comparisons introduce a wider temporal window between the time of in situ measurement and the time of the satellite observation, which can unintentionally reflect on the quality of the satellite retrieval or algorithm performance. The choice of L3 map projection may introduce additional uncertainty by spatially distorting the true location of the satellite retrievals. Each manipulation of satellite data beyond theinstrument's native spatio temporal reference (L2) reduces the applicability of L2 validation results to higher data processing levels.

Ocean Color↗

Fibonacci Grids

Recent years have seen a resurgence of interest in a variety of non-standard computational grids for global numerical prediction. The motivation has been to reduce problems associated with the converging meridians and the polar singularities of conventional regular latitude-longitude grids. A further impetus has come from the adoption of massively parallel computers, for which it is necessary to distribute work equitably across the processors; this is more practicable for some non-standard grids. Desirable attributes of a grid for high-order spatial finite differencing are: (i) geometrical regularity; (ii) a homogeneous and approximately isotropic spatial resolution; (iii) a low proportion of the grid points where the numerical procedures require special customization (such as near coordinate singularities or grid edges). One family of grid arrangements which, to our knowledge, has never before been applied to numerical weather prediction, but which appears to offer several technical advantages, are what we shall refer to as "Fibonacci grids". They can be thought of as mathematically ideal generalizations of the patterns occurring naturally in the spiral arrangements of seeds and fruit found in sunflower heads and pineapples (to give two of the many botanical examples). These grids possess virtually uniform and highly isotropic resolution, with an equal area for each grid point. There are only two compact singular regions on a sphere that require customized numerics. We demonstrate the practicality of these grids in shallow water simulations, and discuss the prospects for efficiently using these frameworks in three-dimensional semi-implicit and semi-Lagrangian weather prediction or climate models.

Swinbank, Richard↗

Automation of assertion testing - Grid and adaptive techniques

Assertions can be used to automate the process of testing software. Two methods for automating the generation of input test data are described in this paper. One method selects the input values of variables at regular intervals in a 'grid'. The other, adaptive testing, uses assertion violations as a measure of errors detected and generates new test cases based on test results. The important features of assertion testing are that: it can be used throughout the entire testing cycle; it provides automatic notification of error conditions; and it can be used with automatic input generation techniques which eliminate the subjectivity in choosing test data.

Andrews, D. M.↗

Development and testing of a simple assimilation technique to derive average wind fields from simulated scatterometer data

A simple algorithm is developed and tested to derive a regularly spaced wind field in a limited area from simulated multiorbit scatterometer data. The data are generated by sampling a time-varying known wind field, the 1000-mb FGGE data, from a simulated scatterometer. A simple assimilation technique is used to derive a regularly spaced (100-km grid) wind field representation of two-day averages from the simulated data. Several test cases are considered, and it is noted that the assimilation technique might be applicable to large-scale ocean or atmospheric models.

Legler, D. M.↗

Communication requirements of sparse Cholesky factorization with nested dissection ordering

Load distribution schemes for minimizing the communication requirements of the Cholesky factorization of dense and sparse, symmetric, positive definite matrices on multiprocessor systems are presented. The total data traffic in factoring an n x n sparse symmetric positive definite matrix representing an n-vertex regular two-dimensional grid graph using n exp alpha, alpha not greater than 1, processors are shown to be O(n exp 1 + alpha/2). It is O(n), when n exp alpha, alpha not smaller than 1, processors are used. Under the conditions of uniform load distribution, these results are shown to be asymptotically optimal.

Naik, Vijay K.↗

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a by-pass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by an ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled 3-component problem were developed in 1994. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 has been developed. It is planned to use the steady-state global solution provided by ENG10 as input to a localized three-dimensional FSI analysis for engine regions where aeroelastic effects may be important.

Felippa, C. A.↗

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

Applications are described of high-performance computing methods to the numerical simulation of complete jet engines. The methodology focuses on the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field elements. New partitioned analysis procedures to treat this coupled three-component problem were developed. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. The NASA-sponsored ENG10 program was used for the global steady state analysis of the whole engine. This program uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 was developed as well as the capability for the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames.

Felippa, C. A.↗

Direct Volume Rendering of Curvilinear Volumes

Direct volume rendering can visualize sampled 3D scalar data as a continuous medium, or extract features. However, it is generally slow. Furthermore, most algorithms for direct volume rendering have assumed rectilinear gridded data. This paper discusses methods for using direct volume rendering when the original volume is curvilinear, i.e. is divided into six-sided cells which are not necessarily equilateral hexahedra. One approach is to ray-cast such volumes directly. An alternative approach is to interpolate the sample volumes to a rectilinear grid, and use this regular volume for rendering. Advantages and disadvantages of the two approaches in terms of speed and image quality are explored.

Vaziri, Arsi↗

Extension of a Kolmogorov Atmospheric Turbulence Model for Time-Based Simulation Implementation

The development of any super/hypersonic aircraft requires the interaction of a wide variety of technical disciplines to maximize vehicle performance. For flight and engine control system design and development on this class of vehicle, realistic mathematical simulation models of atmospheric turbulence, including winds and the varying thermodynamic properties of the atmosphere, are needed. A model which has been tentatively selected by a government/industry group of flight and engine/inlet controls representatives working on the High Speed Civil Transport is one based on the Kolmogorov spectrum function. This report compares the Dryden and Kolmogorov turbulence forms, and describes enhancements that add functionality to the selected Kolmogorov model. These added features are: an altitude variation of the eddy dissipation rate based on Dryden data, the mapping of the eddy dissipation rate database onto a regular latitude and longitude grid, a method to account for flight at large vehicle attitude angles, and a procedure for transitioning smoothly across turbulence segments.

McMinn, John D.↗