Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse 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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗
NASA NTRS1990

Partitioning sparse matrices with eigenvectors of graphs

The problem of computing a small vertex separator in a graph arises in the context of computing a good ordering for the parallel factorization of sparse, symmetric matrices. An algebraic approach for computing vertex separators is considered in this paper. It is shown that lower bounds on separator sizes can be obtained in terms of the eigenvalues of the Laplacian matrix associated with a graph. The Laplacian eigenvectors of grid graphs can be computed from Kronecker products involving the eigenvectors of path graphs, and these eigenvectors can be used to compute good separators in grid graphs. A heuristic algorithm is designed to compute a vertex separator in a general graph by first computing an edge separator in the graph from an eigenvector of the Laplacian matrix, and then using a maximum matching in a subgraph to compute the vertex separator. Results on the quality of the separators computed by the spectral algorithm are presented, and these are compared with separators obtained from other algorithms for computing separators. Finally, the time required to compute the Laplacian eigenvector is reported, and the accuracy with which the eigenvector must be computed to obtain good separators is considered. The spectral algorithm has the advantage that it can be implemented on a medium-size multiprocessor in a straightforward manner.

Pothen, Alex↗
NASA NTRS2007

Terrain Portrayal for Synthetic Vision Systems Head-Down Displays Evaluation Results

A critical component of SVS displays is the appropriate presentation of terrain to the pilot. At the time of this study, the relationship between the complexity of the terrain presentation and resulting enhancements of pilot SA and pilot performance had been largely undefined. The terrain portrayal for SVS head-down displays (TP-HDD) simulation examined the effects of two primary elements of terrain portrayal on the primary flight display (PFD): variations of digital elevation model (DEM) resolution and terrain texturing. Variations in DEM resolution ranged from sparsely spaced (30 arc-sec) to very closely spaced data (1 arc-sec). Variations in texture involved three primary methods: constant color, elevation-based generic, and photo-realistic, along with a secondary depth cue enhancer in the form of a fishnet grid overlay.

Hughes, Monica F.↗
NASA NTRS1990

Application of the TRANAIR rectangular grid approach to the aerodynamic analysis of complex configurations

A numerical method is described which uses a rectangular grid to solve the nonlinear full potential equation about complex configurations. The grid is locally refined to resolve high velocity gradients arising from leading edge expansions or shock waves. The grid penetrates the boundary (described by networks of quadrilateral panels) and is generated automatically. Discrete operators are constructed using the finite element method. The system of nonlinear discrete equations is solved iteratively using a Krylov subspace method preconditioned by an exterior Poisson solver and a direct sparse solver. The primary emphasis is to provide design engineers with an aerodynamic analysis tool (the TRANAIR code) which is accurate, reliable, economical, and flexible to use. Computational results for many interesting configurations are presented.

Johnson, Forrester T.↗
NASA NTRS2016

Spatio-Temporal Variability of Groundwater Storage in India

Groundwater level measurements from 3907 monitoring wells, distributed within 22 major river basins of India, are assessed to characterize their spatial and temporal variability. Ground water storage (GWS) anomalies (relative to the long-term mean) exhibit strong seasonality, with annual maxima observed during the monsoon season and minima during pre-monsoon season. Spatial variability of GWS anomalies increases with the extent of measurements, following the power law relationship, i.e., log-(spatial variability) is linearly dependent on log-(spatial extent).In addition, the impact of well spacing on spatial variability and the power law relationship is investigated. We found that the mean GWS anomaly sampled at a 0.25 degree grid scale closes to unweighted average over all wells. The absolute error corresponding to each basin grows with increasing scale, i.e., from 0.25 degree to 1 degree. It was observed that small changes in extent could create very large changes in spatial variability at large grid scales. Spatial variability of GWS anomaly has been found to vary with climatic conditions. To our knowledge, this is the first study of the effects of well spacing on groundwater spatial variability. The results may be useful for interpreting large scale groundwater variations from unevenly spaced or sparse groundwater well observations or for siting and prioritizing wells in a network for groundwater management. The output of this study could be used to maintain a cost effective groundwater monitoring network in the study region and the approach can also be used in other parts of the globe.

Groundwater monitoring↗
NASA NTRS2003

Terrain Portrayal for Head-Down Displays Flight Test

The Synthetic Vision Systems General Aviation (SVS-GA) element of NASA's Aviation Safety Program is developing technology to eliminate low visibility induced General Aviation (GA) accidents through the application of synthetic vision techniques. SVS displays present computer generated 3-dimensional imagery of the surrounding terrain to greatly enhance pilot's situation awareness (SA), reducing or eliminating Controlled Flight into Terrain (CFIT), as well as Low-Visibility Loss of Control (LVLOC) accidents. In addition to substantial safety benefits, SVS displays have many potential operational benefits that can lead to flight in instrument meteorological conditions (IMC) resembling those conducted in visual meteorological conditions (VMC). Potential benefits could include lower landing minimums, more approach options, reduced training time, etc. SVS conducted research will develop display concepts providing the pilot with an unobstructed view of the outside terrain, regardless of weather conditions and time of day. A critical component of SVS displays is the appropriate presentation of terrain to the pilot. The relationship between the realism of the terrain presentation and resulting enhancements of pilot SA and pilot performance has been largely undefined. Comprised of coordinated simulation and flight test efforts, the terrain portrayal for head-down displays (TP-HDD) test series examined the effects of two primary elements of terrain portrayal: variations of digital elevation model (DEM) resolution and terrain texturing. Variations in DEM resolution ranged from sparsely spaced (30 arc-sec/2,953ft) to very closely spaced data (1 arc-sec/98 ft). Variations in texture involved three primary methods: constant color, elevation-based generic, and photo-realistic, along with a secondary depth cue enhancer in the form of a fishnet grid overlay. The TP-HDD test series was designed to provide comprehensive data to enable design trades to optimize all SVS applications, as well as develop requirements and recommendations to facilitate the implementation and certification of SVS displays. The TP-HDD flight experiment utilized the NASA LaRC Cessna 206 Stationaire and evaluated eight terrain portrayal concepts in an effort to confirm and extend results from the previously conducted TP-HDD simulation experiment. A total of 15 evaluation pilots, of various qualifications, accumulated over 75 hours of dedicated research flight time at Newport News (PHF) and Roanoke (ROA), VA, airports from August through October, 2002. This report will present results from the portion of testing conducted at Roanoke, VA.

Hughes, Monica F.↗
NASA NTRS2017

Development and Assessment of the SMAP Enhanced Passive Soil Moisture Product

Launched in January 2015, the National Aeronautics and Space Administration (NASA) Soil Moisture Active Passive (SMAP) observatory was designed to provide frequent global mapping of high-resolution soil moisture and freeze-thaw state every two to three days using a radar and a radiometer operating at L-band frequencies. Despite a hardware mishap that rendered the radar inoperable shortly after launch, the radiometer continues to operate nominally, returning more than two years of science data that have helped to improve existing hydrological applications and foster new ones.Beginning in late 2016 the SMAP project launched a suite of new data products with the objective of recovering some high-resolution observation capability loss resulting from the radar malfunction. Among these new data products are the SMAP Enhanced Passive Soil Moisture Product that was released in December 2016, followed by the SMAPSentinel-1 Active-Passive Soil Moisture Product in April 2017.This article covers the development and assessment of the SMAP Level 2 Enhanced Passive Soil Moisture Product (L2_SM_P_E). The product distinguishes itself from the current SMAP Level 2 Passive Soil Moisture Product (L2_SM_P) in that the soil moisture retrieval is posted on a 9 km grid instead of a 36 km grid. This is made possible by first applying the Backus-Gilbert optimal interpolation technique to the antenna temperature (TA) data in the original SMAP Level 1B Brightness Temperature Product to take advantage of the overlapped radiometer footprints on orbit. The resulting interpolated TA data then go through various correctioncalibration procedures to become the SMAP Level 1C Enhanced Brightness Temperature Product (L1C_TB_E). The L1C_TB_E product, posted on a 9 km grid, is then used as the primary input to the current operational SMAP baseline soil moisture retrieval algorithm to produce L2_SM_P_E as the final output. Images of the new product reveal enhanced visual features that are not apparent in the standard product. Based on in situ data from core validation sites and sparse networks representing different seasons and biomes all over the world, comparisons between L2_SM_P_E and in situ data were performed for the duration of April 1, 2015 October 30, 2016. It was found that the performance of the enhanced 9 km L2_SM_P_E is equivalent to that of the standard 36 km L2_SM_P, attaining a retrieval uncertainty below 0.040 m(exp 3)/m(exp 3) unbiased root-mean-square error (ubRMSE) and a correlation coefficient above 0.800. This assessment also affirmed that the Single Channel Algorithm using the V-polarized TB channel (SCA-V) delivered the best retrieval performance among the various algorithms implemented for L2_SM_P_E, a result similar to a previous assessment for L2_SM_P.

Enhanced↗
NASA NTRS1985

Thickness of Ridged Plains Materials in Hesperia Planum, Mars

The thickness distribution of plains-forming materials in Hesperia Planum is determined by the diameter of partially buried craters. As with all Martian thickness studies, the distribution of measured thicknesses is sparse due to a low number of suitable partially buried craters. The resulting isopach of plains materials has a low level of confidence but is sufficient for generalized thickness trends. The mean of the thickness estimates is 360 + or - 120 m. The isopach of plains-forming materials exhibits an uneven thickness distribution. Average thickness, determined by grid sampling over the region, is 216 + or - 153 m. Local lenses exceed 500 m. A fifth order trend surface provides a partial match (coefficient of correlation = 0.657) to the isopach.

Dehon, R. A.↗
NASA NTRS2012

Multi-Angle Implementation of Atmospheric Correction for MODIS (MAIAC). Part 3: Atmospheric Correction

This paper describes the atmospheric correction (AC) component of the Multi-Angle Implementation of Atmospheric Correction algorithm (MAIAC) which introduces a new way to compute parameters of the Ross-Thick Li-Sparse (RTLS) Bi-directional reflectance distribution function (BRDF), spectral surface albedo and bidirectional reflectance factors (BRF) from satellite measurements obtained by the Moderate Resolution Imaging Spectroradiometer (MODIS). MAIAC uses a time series and spatial analysis for cloud detection, aerosol retrievals and atmospheric correction. It implements a moving window of up to 16 days of MODIS data gridded to 1 km resolution in a selected projection. The RTLS parameters are computed directly by fitting the cloud-free MODIS top of atmosphere (TOA) reflectance data stored in the processing queue. The RTLS retrieval is applied when the land surface is stable or changes slowly. In case of rapid or large magnitude change (as for instance caused by disturbance), MAIAC follows the MODIS operational BRDF/albedo algorithm and uses a scaling approach where the BRDF shape is assumed stable but its magnitude is adjusted based on the latest single measurement. To assess the stability of the surface, MAIAC features a change detection algorithm which analyzes relative change of reflectance in the Red and NIR bands during the accumulation period. To adjust for the reflectance variability with the sun-observer geometry and allow comparison among different days (view geometries), the BRFs are normalized to the fixed view geometry using the RTLS model. An empirical analysis of MODIS data suggests that the RTLS inversion remains robust when the relative change of geometry-normalized reflectance stays below 15%. This first of two papers introduces the algorithm, a second, companion paper illustrates its potential by analyzing MODIS data over a tropical rainforest and assessing errors and uncertainties of MAIAC compared to conventional MODIS products.

Lyapustin, A.↗
NASA NTRS1990

Potential, velocity, and density fields from sparse and noisy redshift-distance samples - Method

A method for recovering the three-dimensional potential, velocity, and density fields from large-scale redshift-distance samples is described. Galaxies are taken as tracers of the velocity field, not of the mass. The density field and the initial conditions are calculated using an iterative procedure that applies the no-vorticity assumption at an initial time and uses the Zel'dovich approximation to relate initial and final positions of particles on a grid. The method is tested using a cosmological N-body simulation 'observed' at the positions of real galaxies in a redshift-distance sample, taking into account their distance measurement errors. Malmquist bias and other systematic and statistical errors are extensively explored using both analytical techniques and Monte Carlo simulations.

Dekel, Avishai↗
NASA NTRS1987

Solving very large, sparse linear systems on mesh-connected parallel computers

The implementation of Pan and Reif's Parallel Nested Dissection (PND) algorithm on mesh connected parallel computers is described. This is the first known algorithm that allows very large, sparse linear systems of equations to be solved efficiently in polylog time using a small number of processors. How the processor bound of PND can be matched to the number of processors available on a given parallel computer by slowing down the algorithm by constant factors is described. Also, for the important class of problems where G(A) is a grid graph, a unique memory mapping that reduces the inter-processor communication requirements of PND to those that can be executed on mesh connected parallel machines is detailed. A description of an implementation on the Goodyear Massively Parallel Processor (MPP), located at Goddard is given. Also, a detailed discussion of data mappings and performance issues is given.

Opsahl, Torstein↗
NASA NTRS2013

The Utility of the Real-Time NASA Land Information System Data for Drought Monitoring Applications

Measurements of soil moisture are a crucial component for the proper monitoring of drought conditions. The large spatial variability of soil moisture complicates the problem. Unfortunately, in situ soil moisture observing networks typically consist of sparse point observations, and conventional numerical model analyses of soil moisture used to diagnose drought are of coarse spatial resolution. Decision support systems such as the U.S. Drought Monitor contain drought impact resolution on sub-county scales, which may not be supported by the existing soil moisture networks or analyses. The NASA Land Information System, which is run with 3 km grid spacing over the eastern United States, has demonstrated utility for monitoring soil moisture. Some of the more useful output fields from the Land Information System are volumetric soil moisture in the 0-10 cm and 40-100 cm layers, column-integrated relative soil moisture, and the real-time green vegetation fraction derived from MODIS (Moderate Resolution Imaging Spectroradiometer) swath data that are run within the Land Information System in place of the monthly climatological vegetation fraction. While these and other variables have primarily been used in local weather models and other operational forecasting applications at National Weather Service offices, the use of the Land Information System for drought monitoring has demonstrated utility for feedback to the Drought Monitor. Output from the Land Information System is currently being used at NWS Huntsville to assess soil moisture, and to provide input to the Drought Monitor. Since feedback to the Drought Monitor takes place on a weekly basis, weekly difference plots of column-integrated relative soil moisture are being produced by the NASA Short-term Prediction Research and Transition Center and analyzed to facilitate the process. In addition to the Drought Monitor, these data are used to assess drought conditions for monthly feedback to the Alabama Drought Monitoring and Impact Group and the Tennessee Drought Task Force, which are comprised of federal, state, and local agencies and other water resources professionals.

White, Kristopher D.↗
NASA NTRS1994

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗
NASA NTRS1981

Sensitivity of wind analyses to barotropic energetics during summer

Three analysis schemes - the extended correction method of the University of Hawaii at Manoa (UHM), the National Meteorological Center (NMC) Hough analysis, and Sadler's subjective stream analysis (SUB) - were used to obtain grid-point winds using data for a five-day period in late summer, and barotropic energetics computations obtained using the three schemes for six limited regions in the Pacific and Indian Oceans were compared. Few differences were observed between UHM and NMC stream-function analyses, but significant differences were found in the wind fields and energetics obtained using the subjective method and the two objective methods in data-sparse regions. Major features in the local barotropic criteria field were found to be less sensitive to the analysis scheme used.

Sumi, A.↗
NASA NTRS1993

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗
NASA NTRS2007

Using High Frequency Passive Microwave, A-train, and TRMM Data to Evaluate Hydrometer Structure in the NASA GEOS-5 Data Assimilation System

Validating water vapor and prognostic condensate in global models remains a challenging research task. Model parameterizations are still subject to a large number of tunable parameters; furthermore, accurate and representative in situ observations are very sparse, and satellite observations historically have significant quantitative uncertainties. Progress on improving cloud / hydrometeor fields in models stands to benefit greatly from the growing inventory ofA-Train data sets. ill the present study we are using a variety of complementary satellite retrievals of hydrometeors to examine condensate produced by the emerging NASA Modem Era Retrospective Analysis for Research and Applications, MERRA, and its associated atmospheric general circulation model GEOS5. Cloud and precipitation are generated by both grid-scale prognostic equations and by the Relaxed Arakawa-Schubert (RAS) diagnostic convective parameterization. The high frequency channels (89 to 183.3 GHz) from AMSU-B and MRS on NOAA polar orbiting satellites are being used to evaluate the climatology and variability of precipitating ice from tropical convective anvils. Vertical hydrometeor structure from the Tropical Rainfall Measuring Mission (TRMM) and CloudSat radars are used to develop statistics on vertical hydrometeor structure in order to better interpret the extensive high frequency passive microwave climatology. Cloud liquid and ice water path data retrieved from the Moderate Resolution Imaging Spectroradiometer, MODIS, are used to investigate relationships between upper level cloudiness and tropical deep convective anvils. Together these data are used to evaluate cloud / ice water path, gross aspects of vertical hydrometeor structure, and the relationship between cloud extent and surface precipitation that the MERRA reanalysis must capture.

Robertson, Franklin↗
NASA NTRS1997

Global Analysis of Empirical Relationships Between Annual Climate and Seasonality of NDVI

This study describes the use of satellite data to calibrate a new climate-vegetation greenness function for global change studies. We examined statistical relationships between annual climate indexes (temperature, precipitation, and surface radiation) and seasonal attributes of the AVHRR Normalized Difference Vegetation Index (NDVI) time series for the mid-1980s in order to refine our empirical understanding of intraannual patterns and global abiotic controls on natural vegetation dynamics. Multiple linear regression results using global l(sup o) gridded data sets suggest that three climate indexes: growing degree days, annual precipitation total, and an annual moisture index together can account to 70-80 percent of the variation in the NDVI seasonal extremes (maximum and minimum values) for the calibration year 1984. Inclusion of the same climate index values from the previous year explained no significant additional portion of the global scale variation in NDVI seasonal extremes. The monthly timing of NDVI extremes was closely associated with seasonal patterns in maximum and minimum temperature and rainfall, with lag times of 1 to 2 months. We separated well-drained areas from l(sup o) grid cells mapped as greater than 25 percent inundated coverage for estimation of both the magnitude and timing of seasonal NDVI maximum values. Predicted monthly NDVI, derived from our climate-based regression equations and Fourier smoothing algorithms, shows good agreement with observed NDVI at a series of ecosystem test locations from around the globe. Regions in which NDVI seasonal extremes were not accurately predicted are mainly high latitude ecosystems and other remote locations where climate station data are sparse.

Potter, C. S.↗
NASA NTRS2013

A Global Ozone Climatology from Ozone Soundings via Trajectory Mapping: A Stratospheric Perspective

This study explores a domain-filling trajectory approach to generate a global ozone climatology from sparse ozonesonde data. Global ozone soundings of 51,898 profiles at 116 stations over 44 years (1965-2008) are used, from which forward and backward trajectories are performed for 4 days, driven by a set of meteorological reanalysis data. Ozone mixing ratios of each sounding from the surface to 26 km altitude are assigned to the entire path along the trajectory. The resulting global ozone climatology is archived monthly for five decades from the 1960s to the 2000s with grids of 5 degree 5 degree 1 km (latitude, longitude, and altitude). It is also archived yearly from 1965 to 2008. This climatology is validated at 20 ozonesonde stations by comparing the actual ozone sounding profile with that found through the trajectories, using the ozone soundings at all the stations except one being tested. The two sets of profiles are in good agreement, both individually with correlation coefficients between 0.975 and 0.998 and root mean square (RMS) differences of 87 to 482 ppbv, and overall with a correlation coefficient of 0.991 and an RMS of 224 ppbv. The ozone climatology is also compared with two sets of satellite data, from the Satellite Aerosol and Gas Experiment (SAGE) and the Optical Spectrography and InfraRed Imager System (OSIRIS). Overall, the ozone climatology compares well with SAGE and OSIRIS data by both seasonal and zonal means. The mean difference is generally under 20 above 15 km. The comparison is better in the northern hemisphere, where there are more ozonesonde stations, than in the southern hemisphere; it is also better in the middle and high latitudes than in the tropics, where assimilated winds are imperfect in some regions. This ozone climatology can capture known features in the stratosphere, as well as seasonal and decadal variations of these features. Furthermore, it provides a wealth of detail about longitudinal variations in the stratosphere such as the spring ozone maximum over the Canadian Arctic. It also covers higher latitudes than current satellite data. The climatology shows clearly the depletion of ozone from the 1970s to the mid 1990s and ozone recovery in the 2000s. When this climatology is used as the upper boundary condition in an Environment Canada operational chemical forecast model, the forecast is improved in the vicinity of the upper tropospherelower stratosphere region. As this ozone climatology is neither dependent on a priori data or photochemical modeling, it provides independent information and insight that can supplement satellite data and model simulations and enhance our understanding of stratospheric ozone.

SAGE↗