Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial interpolation”

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 217 records · Page 12

Mean Winds and Tides over Urbana, Illinois During November 16-22, 1981

The Urbana Meteor Radar uses a high power pulse transmitter to achieve high echo rates, and an interferometer to spatially locate the individual echoes. The irregularly spaced observations of range, radial velocity, and direction cosines are reduced to a rectangular grid, evenly spaced in time and height, of estimates of the zonal and meridional wind. This is accomplished using the method of GROVES whereby a constant horizontal wind is fitted, using least squares, to all echoes occurring within a fixed time/height interval. This interpolation region is then stepped in time and height. the size of the interpolation region is 2 hours by 5 km, the time step is approximately 11 minutes, and the height step is 2 km. The filtering of high temporal and spatial frequencies is somewhat reduced by weighting the echoes according to their distance from the center of the interpolation box.

Tetenbaum, D.↗

Computation of three-dimensional inviscid flow over hypersonic missile configurations using the GIM code

A three-dimensional computational technique was used to obtain flowfield solutions to the Euler equations over selected hypersonic missile configurations. The General Interpolants Method (GIM) computer code was used with interpolation functions in an algebraic approach to generate a discrete computational grid for each configuration. The spatial marching version of the GIM code, which treats the parabolized Navier-Stokes (PNS) equations or the Euler equations with a shock capturing, 'MacCormack-like' scheme, was used to advance the solution hyperbolically over each configuration. The inviscid flowfield solutions over the two three-dimensional missile configurations, calculated using the GIM hyperbolic scheme, are presented here. The flow field over a wing/body configuration at zero degree angle of attack is presented. Flow over the fuselage of a tactical missile, termed the TAME 10, at both zero degree and 7.5 degree angles of attack is presented. In addition, an inviscid, two-dimensional analysis of an inlet configuration designed to mount on the TAME 10 is included. Contour maps of velocity and pressure are included for each configuration. Comparison of calculation and data show good agreement.

Xiques, K. E.↗

Estimating Long Term Surface Soil Moisture in the GCIP Area From Satellite Microwave Observations

Soil moisture is an important component of the water and energy balances of the Earth's surface. Furthermore, it has been identified as a parameter of significant potential for improving the accuracy of large-scale land surface-atmosphere interaction models. However, accurate estimates of surface soil moisture are often difficult to make, especially at large spatial scales. Soil moisture is a highly variable land surface parameter, and while point measurements are usually accurate, they are representative only of the immediate site which was sampled. Simple averaging of point values to obtain spatial means often leads to substantial errors. Since remotely sensed observations are already a spatially averaged or areally integrated value, they are ideally suited for measuring land surface parameters, and as such, are a logical input to regional or larger scale land process models. A nine-year database of surface soil moisture is being developed for the Central United States from satellite microwave observations. This region forms much of the GCIP study area, and contains most of the Mississippi, Rio Grande, and Red River drainages. Daytime and nighttime microwave brightness temperatures were observed at a frequency of 6.6 GHz, by the Scanning Multichannel Microwave Radiometer (SMMR), onboard the Nimbus 7 satellite. The life of the SMMR instrument spanned from Nov. 1978 to Aug. 1987. At 6.6 GHz, the instrument provided a spatial resolution of approximately 150 km, and an orbital frequency over any pixel-sized area of about 2 daytime and 2 nighttime passes per week. Ground measurements of surface soil moisture from various locations throughout the study area are used to calibrate the microwave observations. Because ground measurements are usually only single point values, and since the time of satellite coverage does not always coincide with the ground measurements, the soil moisture data were used to calibrate a regional water balance for the top 1, 5, and 10 cm surface layers in order to interpolate daily surface moisture values. Such a climate-based approach is often more appropriate for estimating large-area spatially averaged soil moisture because meteorological data are generally more spatially representative than isolated point measurements of soil moisture. Vegetation radiative transfer characteristics, such as the canopy transmissivity, were estimated from vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the 37 GHz Microwave Polarization Difference Index (MPDI). Passive microwave remote sensing presents the greatest potential for providing regular spatially representative estimates of surface soil moisture at global scales. Real time estimates should improve weather and climate modelling efforts, while the development of historical data sets will provide necessary information for simulation and validation of long-term climate and global change studies.

Owe, Manfred↗

Development of a Shuttle Infrared Telescope Facility /SIRTF/ fine guidance sensor

Fine guidance technology development for the Shuttle Infrared Telescope Facility (SIRTF) centers upon the use of a single multiple-star-tracking sensor to provide the position information necessary to produce three-axis attitude control signals for precision payload pointing. The effort described in this paper is concerned with the development of a fine guidance sensor that employs a high-density charge-coupled imaging device for producing position information signals by using star fields. Multiple star position information produces three-axis position error signals that are used to update inertial reference gyros. The sensor employs advanced position interpolation algorithms to enhance field-of-view resolution and to correct for optical aberrations inherent in spatially chopped star images resulting from the telescope's movable secondary mirror. Operation of the sensor is under the control of a high-performance microcomputer that provides both autonomy and flexibility in a guidance application.

Salomon, P. M.↗

Kriging in the Shadows: Geostatistical Interpolation for Remote Sensing

It is often useful to estimate obscured or missing remotely sensed data. Traditional interpolation methods, such as nearest-neighbor or bilinear resampling, do not take full advantage of the spatial information in the image. An alternative method, a geostatistical technique known as indicator kriging, is described and demonstrated using a Landsat Thematic Mapper image in southern Chiapas, Mexico. The image was first classified into pasture and nonpasture land cover. For each pixel that was obscured by cloud or cloud shadow, the probability that it was pasture was assigned by the algorithm. An exponential omnidirectional variogram model was used to characterize the spatial continuity of the image for use in the kriging algorithm. Assuming a cutoff probability level of 50%, the error was shown to be 17% with no obvious spatial bias but with some tendency to categorize nonpasture as pasture (overestimation). While this is a promising result, the method's practical application in other missing data problems for remotely sensed images will depend on the amount and spatial pattern of the unobscured pixels and missing pixels and the success of the spatial continuity model used.

Rossi, Richard E.↗

Investigating Uncertainty in Predicting Carbon Dynamics in North American Biomes: Putting Support-Effect Bias in Perspective

A fundamental strategy in NASA's Earth Observing System's (EOS) monitoring of vegetation and its contribution to the global carbon cycle is to rely on deterministic, process-based ecosystem models to make predictions of carbon flux over large regions. These models are parameterized (that is, the input variables are derived) using remotely sensed images such as those from the Moderate Resolution Imaging Spectroradiometer (MODIS), ground measurements and interpolated maps. Since early applications of these models, investigators have noted that results depend partly on the spatial support of the input variables. In general, the larger the support of the input data, the greater the chance that the effects of important components of the ecosystem will be averaged out. A review of previous work shows that using large supports can cause either positive or negative bias in carbon flux predictions. To put the magnitude and direction of these biases in perspective, we must quantify the range of uncertainty on our best measurements of carbon-related variables made on equivalent areas. In other words, support-effect bias should be placed in the context of prediction uncertainty from other sources. If the range of uncertainty at the smallest support is less than the support-effect bias, more research emphasis should probably be placed on support sizes that are intermediate between those of field measurements and MODIS. If the uncertainty range at the smallest support is larger than the support-effect bias, the accuracy of MODIS-based predictions will be difficult to quantify and more emphasis should be placed on field-scale characterization and sampling. This talk will describe methods to address these issues using a field measurement campaign in North America and "upscaling" using geostatistical estimation and simulation.

Dungan, Jennifer L.↗

Fusion of Surface Ceilometer Data and Satellite Cloud Retrievals in 2D Mesh Interpolating Model with Clustering

For accurate cloud ceiling information, a data fusion approach is proposed that utilizes satellite data to extend surface station information to much wider areas. Cloud base height (CBH) retrieved from satellite observations provides for much larger spatial coverage and higher resolution. The direct comparison of GOES-16 CBH with surface station ceiling yields a local bias that has to be corrected for in the initial GOES-16 cloud base information. This sparsely sampled bias correction presents an irregular 2D mesh of control points, which is then interpolated by constructing a continuous smooth field using polyharmonic splines. The influence of remote stations is restricted by grouping the control points into clusters depending on an effective distance. This cluster-based approach allows for constructing separate spline surfaces corresponding to physically different clouds. The obtained continuous bias correction function is then applied to the entire GOES-16 pixel level CBH except for areas far away from surface stations in data sparse regions such as offshore. The described method is currently being tested using daytime-only observations over the central and eastern United States. Overall, this approach has potential to provide more accurate, high spatial resolution cloud ceiling information for the aviation community.

Khlopenkov, Konstantin↗

Quantifying Libya-4 Surface Reflectance Heterogeneity With WorldView-1, 2 and EO-1 Hyperion

The land surface imaging (LSI) virtual constellation approach promotes the concept of increasing Earth observations from multiple but disparate satellites. We evaluated this through spectral and spatial domains, by comparing surface reflectance from 30-m Hyperion and 2-m resolution WorldView-2 (WV-2) data in the Libya-4 pseudoinvariant calibration site. We convolved and resampled Hyperion to WV-2 bands using both cubic convolution and nearest neighbor (NN) interpolation. Additionally, WV-2 and WV-1 same-date imagery were processed as a cross-track stereo pair to generate a digital terrain model to evaluate the effects from large (>70 m) linear dunes. Agreement was moderate to low on dune peaks between WV-2 and Hyperion (R2 <; 0.4) but higher in areas of lower elevation and slope (R2 > 0.6). Our results provide a satellite sensor intercomparison protocol for an LSI virtual constellation at high spatial resolution, which should start with geolocation of pixels, followed by NN interpolation to avoid tall dunes that enhance surface reflectance differences across this internationally utilized site.

Hyperion↗

Spatial characteristics of observed precipitation fields: A catalog of summer storms in Arizona, Volume 2

The parameters of the conceptual model are evaluated from the analysis of eight years of summer rainstorm data from the dense raingage network in the Walnut Gulch catchment near Tucson, Arizona. The occurrence of measurable rain at any one of the 93 gages during a noon to noon day defined a storm. The total rainfall at each of the gages during a storm day constituted the data set for a single storm. The data are interpolated onto a fine grid and analyzed to obtain: an isohyetal plot at 2 mm intervals, the first three moments of point storm depth, the spatial correlation function, the spatial variance function, and the spatial distribution of the total storm depth. The description of the data analysis and the computer programs necessary to read the associated data tapes are presented.

Fennessey, N. M.↗

ArborX 2.0

ArborX library tackles a problem of efficiently finding geometric objects that are close in space. Variations of this problem, such as finding the nearest neighbors of a point, or finding all objects within a certain distance, are inherent components of applications in many fields. The data may be large so that solving the problem efficiently may require significant computational resources, such as multiple processors or accelerators such as general purpose GPUs. ArborX' main advantage in its ability to solve large problems efficiently utilizing a combination of distributed and on-node parallelism. ArborX can be run efficiently on a wide variety of hardware, including GPUs from different vendors, which distinguishes it from other available libraries which typically choose only few of these. The other advantage is that it supports both types of user problems: spatial problems (useful for intersections and finding objects within certain distance), and nearest neighbor problems. ArborX also supports flexible interface in its interaction with a user. Particularly, it allows a user to call user's own function on a positive match, a functionality not rarely available in other libraries. ArborX implements construction and traversal algorithms using efficient tree structures, such as bounding volume hierarchy (BVH). At its core, ArborX uses linear BVH for its low construction cost and sufficient quality. ArborX implements both spatial and nearest-neighbor traversal algorithms. ArborX also provides several clustering algorithms (minimum spanning tree, DBSCAN, HDBSCAN*), interpolation using minimum least squares and ray tracing. ArborX is written using C++, and is parallelized using the message passing interface (MPI) for the distributed communication, and the Kokkos library for on-node parallelism. This approach allows ArborX to be run on a wide variety of hardware, from common laptops and desktops to supercomputers while using the same codebase.

Prokopenko, Andrey [Oak Ridge National Laboratory ↗

UFNet: Joint U-Net and Fully Connected Neural Network to Bias Correct Precipitation Predictions from

Paper information. Shuang Yu, Indrasis Chakraborty, Gemma J. Anderson, Donald D. Lucas, Yannic Lops, and Daniel Galea. UFNet: Joint U-Net and fully connected neural network to bias correct precipitation predictions from climate models. Artificial Intelligence for the Earth Systems, 2024. Overview. This work develops the UFNet methodology to correct E3SM historical precipitation projection bias. The UFNet deep learning framework consists of a two-part architecture: a U-Net convolutional network to capture the spatiotemporal distribution of precipitation and a fully connected network to capture the distribution of higher-order statistics. The joint network, termed UFNet, can simultaneously improve the spatial structure of the modeled precipitation and capture the distribution of extreme precipitation values. Below we provide guidance for applying UFNet to correct the Energy Exascale Earth System Model (E3SM; Golaz et al. 2019) daily precipitation projection over the contiguous United States (CONUS). Getting started 1. Obtain the historical climate simulation and observation data. The E3SM historical simulation data are available through https://aims2.llnl.gov/search/cmip6/. The CPC unified gauge-based analysis of daily precipitation can be found through https://psl.noaa.gov/data/gridded/data.cpc.globalprecip.html. The ECMWF atmospheric reanalysis of the 20th century (ERA-20C) data are available through https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era-20c. The spatial resolution of E3SM and observed datasets are both regridded to a common 1° resolution grid using conservative interpolation. The regridded E3SM, CPC and ERA-20C with 1° resolution can be found throught ./data/. 2. Train the fully connected network (DNN) Python train_dnn.py 3. Train the UFNet Python train_ufnet.py 4. Evaluation and compared with the baseline Python evaluation.py

Lucas, DonaldD↗

Appalachian Basin Play Fairway Analysis Thermal Risk Factor and Quality Analyses

*This submission revises the analysis and products for Thermal Quality Analysis for the northern half of the Appalachian Basin (https://gdr.openei.org/submissions/638)* This submission is one of five major parts of a Low Temperature Geothermal Play Fairway Analysis. Phase 1 of the project identified potential Geothermal Play Fairways within the Appalachian basin of Pennsylvania, West Virginia and New York. This submission includes a subset of the necessary shapefiles, rasters, datasets, code, and references to code repositories that were used to create the thermal resource and risk factor maps as part of the project. This subset is those contents that were improved upon during calendar year 2016. Figures are provided as examples of some shapefiles and rasters. See also: Final Report: Low Temperature Geothermal Play Fairway Analysis for the Appalachian Basin (https://gdr.openei.org/submissions/899). The 2015 data submission should be visited to obtain: 1) the regional standardized 1 square km grid used in the project as points (cell centers), polygons, and as a raster, 2) the raw well data for the state well temperature databases, 3) the COSUNA section shapefile and formation thermal conductivities by state as *.xlsx files, 4) the sediment thickness map and 30 m Digital Elevation Model for the Appalachian Basin as GeoTIFF raster files, 5) the BHT correction sections shapefile and drilling fluid databases as *.csv files, 6) the unbuffered interpolation regions as shapefiles, 7) several 50 km buffered interpolation regions as shapefiles, 8) several gridded interpolation regions as raster files, 9) an R script for organizing the thermal data and running the local spatial outlier analysis, 10) shapefiles and rasters for the prediction, uncertainty, and cross validation of the temperature at 1.5 km, 2.5 km, and 3.5 km depth, 11) shapefiles and rasters for the prediction, uncertainty, and cross validation depth to 100 degrees C, 12) an ArcGIS toolbox for thermal risk factor models, 13) an ArcGIS model for extracting results specific to each county of interest, 14) thermal resource cross section plots, 15) the geothermal Play Fairways.

15 GEOTHERMAL ENERGY↗

Solar physics applications of computer graphics and image processing

Computer graphics devices coupled with computers and carefully developed software provide new opportunities to achieve insight into the geometry and time evolution of scalar, vector, and tensor fields and to extract more information quickly and cheaply from the same image data. Two or more different fields which overlay in space can be calculated from the data (and the physics), then displayed from any perspective, and compared visually. The maximum regions of one field can be compared with the gradients of another. Time changing fields can also be compared. Images can be added, subtracted, transformed, noise filtered, frequency filtered, contrast enhanced, color coded, enlarged, compressed, parameterized, and histogrammed, in whole or section by section. Today it is possible to process multiple digital images to reveal spatial and temporal correlations and cross correlations. Data from different observatories taken at different times can be processed, interpolated, and transformed to a common coordinate system.

Altschuler, M. D.↗

Accuracy of diagnostic heat and moisture budgets using SESAME-79 field data as revealed by observing system simulation experiments

Observing system simulation experiments are used to investigate the accuracy of diagnostic heat and moisture budgets which employ the AVE-SESAME 1979 data. The time-including, four-dimensional data set of a mesoscale model is used to simulate rawinsonde observations from AVE-SESAME 1979. The 5 C/day (heat budget) and 2 g/kg per day (moisture budget) magnitudes of error obtained indicate difficulties in the diagnosis of the heating rate in weak convective systems. The influences exerted by observational frequency, objective analysis, observational density, vertical interpolation, and observational errors on the budget are also studied, and it is found that the temporal and spatial resolution of the SESAME regional network is marginal for diagnosing convective effects on a horizontal time scale of 550 x 550 km.

Kuo, Y.-H.↗

Finite-volume application of high-order ENO schemes to two-dimensional boundary-value problems

Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.

Casper, Jay↗

bifacialVF - Bifacial PV View Factor model [SWR 17-58]

bifacialvf is a self-contained view factor (or configuration factor) model which replicates a 5-row PV system of infinite extent perpendicular to the module rows. Single-axis tracking is supported, and hourly output files based on TMY inputs are saved. Spatial nonuniformity is reported, with multiple rear-facing irradiances collected on the back of each module row. Bilinear interpolation code add-on to bifacialvf (description below) to pre-generate IV arrays and bifacial coefficients, and to examine the energy production with back side irradiance mismatch for either a portrait or landscape module. Included are IV curves and bifacial info for a Yingli (standard) module. Based on the publication: "A Practical Irradiance Model for Bifacial PV Modules" B. Marion, S. MacAlpine, C. Deline, A. Asgharzadeh, F. Toor, D. Riley, J. Stein, C. Hansen 2017 IEEE Photovoltaic Specialists Conference, Washington DC, 2017 URL: https://www.nrel.gov/docs/fy17osti/67847.pdf

Deline, Chris↗

bifacialVF - Bifacial PV View Factor model 2022 [SWR 17-58]

bifacialvf is a self-contained view factor (or configuration factor) model which replicates a 5-row PV system of infinite extent perpendicular to the module rows. Single-axis tracking is supported, and hourly output files based on TMY inputs are saved. Spatial nonuniformity is reported, with multiple rear-facing irradiances collected on the back of each module row. Bilinear interpolation code add-on to bifacialvf (description below) to pre-generate IV arrays and bifacial coefficients, and to examine the energy production with back side irradiance mismatch for either a portrait or landscape module. Included are IV curves and bifacial info for a Yingli (standard) module. Based on the publication: "A Practical Irradiance Model for Bifacial PV Modules" B. Marion, S. MacAlpine, C. Deline, A. Asgharzadeh, F. Toor, D. Riley, J. Stein, C. Hansen 2017 IEEE Photovoltaic Specialists Conference, Washington DC, 2017 URL: https://www.nrel.gov/docs/fy17osti/67847.pdf

Deline, Chris↗

Deep Learning-based Non-Stationary Bias Correction (NSBC)

This work develops the NSBC (non-stationary bias correction) methodology to correct temperature projection bias from E3SM. The NSBC deep learning framework consists of a three-part architecture: an auto-encoder for compressing the spatial information, an LSTM for predicting annual temperature mean, and a U-Net for capturing the residual bias in temperature. The non-stationary bias correction (NSBC) framework can correct the non-stationarity of the biases of the climate models, which significantly improves the accuracy of future temperature prediction and improves the overestimation of extreme high temperatures that many existing bias correction methods suffer from. Getting started 1. Obtain the historical climate simulation and observation data. The E3SM simulation data are available through https://aims2.llnl.gov/search/cmip6/. The pseudo observations, the Geophysical Fluid Dynamics Laboratory (GFDL)-ESM4 model (Krasting et al., 2018) are available through https://aims2.llnl.gov/search/cmip6/. The spatial resolution of E3SM and pseudo observation datasets are both regridded to a common 1° resolution grid using conservative interpolation. The regridded E3SM and pseudo observation with 1° resolution can be found throught ./data/. 2. Train the Auto-encoder model. Python 0-autoencoder.py 3. Train the LSTM Python 1-LSTM.py 4. Generate the annual mean temperature based on trained LSTM Python 2-generate_annual_mean_LSTM.py 5. Train the U-Net. Python 3-unet.py 6. Evaluation and compared with the baseline Python 4_evaluation.py Is there a deadline approaching that requires the release of yo

Lucas, Donald↗