Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

BonVoision: Leveraging Spatial Data Smoothness For Recovery From Memory Soft Errors

The increasing soft error rates in memory systems raise an emerging concern for modern computing systems. As a result, detectable but uncorrectable errors (DUEs) become potentially more frequent and affect HPC applications. Today, upon encountering a DUE, HPC applications crash, incurring significant performance, storage, and energy overheads. In this paper, we propose a technique to continue application execution past a DUE through the repair of the corrupted memory data by leveraging spatial data smoothness. We present BonVoision, a run-time system that intercepts DUE events, analyzes the binary to identify data elements in the structural neighborhood of the event, and fixes the corrupted data elements by interpolating from the values in their neighborhood. Our evaluation demonstrates that BonVoision incurs negligible overhead and outperforms other recovery strategies by a factor of 2×, on average. We demonstrate that BonVoision also improves the efficiency of existing checkpointing/restart schemes by approximately increasing the optimal checkpoint interval by 23%.

Fang, Bo↗
NASA NTRS1987

Geometric processing of digital images of the planets

New procedures and software have been developed for geometric transformation of images to support digital cartography of the planets. The procedures involve the correction of spacecraft camera orientation of each image with the use of ground control and the transformation of each image to a Sinusoidal Equal-Area map projection with an algorithm which allows the number of transformation calculations to vary as the distortion varies within the image. When the distortion is low in an area of an image, few transformation computations are required, and most pixels can be interpolated. When distortion is extreme, the location of each pixel is computed. Mosaics are made of these images and stored as digital databases. Completed Sinusoidal databases may be used for digital analysis and registration with other spatial data. They may also be reproduced as published image maps by digitally transforming them to appropriate map projections.

Edwards, Kathleen↗
NASA NTRS2000

Hydrological Balance in Tropical Cyclones with Scatterometer and TRMM Data

Precipitation over oceans can be estimated from the radar and the microwave radiometer of the Tropical Rain Measuring Mission (TRMM). It can also be estimated from the divergence of the vertically integrated water vapor transport, through the conservation principle, assuming evaporation is relatively small. In tropical cyclones, the divergence of vertically integrated water vapor is highly dependent on the vertical transport and, therefore, on the wind divergence. Spaceborne scatterometers provide surface wind velocity and, therefore, surface wind divergence at spatial resolutions that are much higher than products of numerical weather prediction (NWP). In this study, ocean surface winds derived from the observations of space-based scatterometers and surface precipitation measured by TRMM were objectively interpolated to the same time and location during the passage of a tropical cyclone. Surface precipitation distribution was derived from wind and humidity profiles provided by NWP. When the surface level winds of NWP were replaced by the scatterometer winds, the surface precipitation patterns computed with the conservation method were found to be significantly changed and the new patterns are much closer in agreement with the patterns observed by TRMM.

Liu, W. Timothy↗
NASA NTRS2000

Hurricane Directional Wave Spectrum Spatial Variation in the Open Ocean and at Landfall

The sea surface directional wave spectrum was measured for the first time in all quadrants of a hurricane in open water using the NASA airborne scanning radar altimeter (SRA) carried aboard one of the NOAA WP-3D hurricane hunter aircraft at 1.5 km height. The SRA measures the energetic portion of the directional wave spectrum by generating a topographic map of the sea surface. At 8 Hz, the SRA sweeps a radar beam of 1 E half-power width (two-way) across the aircraft ground track over a swath equal to 0.8 of the aircraft height, simultaneously measuring the backscattered power at its 36 GHz (8.3 mm) operating frequency and the range to the sea surface at 64 positions. These slant ranges are multiplied by the cosine of the incidence angles to determine the vertical distances from the aircraft to the sea surface. Subtracting these distances from the aircraft height produces the sea surface elevation map. The sea surface topography is interpolated to a uniform grid, transformed by a two-dimensional FFT, and Doppler corrected. The open-ocean data were acquired on 24 August 1998 when hurricane Bonnie was east of the Bahamas and moving slowly to the north. Individual waves with heights up to 18 m were observed and the spatial variation of the wave field was dramatic. The dominant waves generally propagated at significant angles to the downwind direction. At some positions there were three different wave fields of comparable energy crossing each other. The NOAA aircraft spent over five hours within 180 km of the hurricane Bonnie eye, and made five eye penetrations. A 3-minute animation of the directional wave spectrum spatial variation over this period will be shown as well as summary plots of the wave field spatial variation. On 26 August 1998, the NOAA aircraft flew at 2.2 km height when hurricane Bonnie was making landfall near Wilmington, NC, documenting the directional wave spectrum in the region between Charleston, SC and Cape Hatteras, NC. The aircraft ground track included both segments along the shoreline and Pamlico Sound as well as far offshore. An animation of the directional wave spectrum spatial variation at landfall will be presented and contrasted with the spatial variation when Bonnie was in the open ocean on 24 August 1998.

Walsh, E. J.↗
NASA NTRS2000

Hurricane Directional Wave Spectrum Spatial Variation in the Open Ocean and at Landfall

The sea surface directional wave spectrum was measured for the first time in all quadrants of a hurricane in open water using the NASA airborne scanning radar altimeter (SRA) carried aboard one of the NOAA WP-3D hurricane hunter aircraft at 1.5 km height. The SRA measures the energetic portion of the directional wave spectrum by generating a topographic map of the sea surface. At 8 Hz, the SRA sweeps a radar beam of 1' half-power width (two-way) across the aircraft ground track over a swath equal to 0.8 of the aircraft height, simultaneously measuring the backscattered power at its 36 GHz (8.3 mm) operating frequency and the range to the sea surface at 64 positions. These slant ranges are multiplied by the cosine of the off-nadir angles to determine the vertical distances from the aircraft to the sea surface. Subtracting these distances from the aircraft height produces the sea surface elevation map. The sea surface topography is interpolated to a uniform grid, transformed by a two dimensional FFT, and Doppler corrected. The open-ocean data were acquired on 24 August 1998 when hurricane Bonnie was east of the Bahamas and moving toward 330 deg at about 5 m/s. Individual waves up to 18 m height were observed and the spatial variation of the wave field was dramatic. The dominant waves generally propagated at significant angles to the downwind direction. At some positions there were three different wave fields of comparable energy crossing each other. The NOAA aircraft spent over five hours within 180 km of the eye, and made five eye penetrations. On 26 August 1998, the NOAA aircraft flew at 2.2 km height when hurricane Bonnie was making landfall near Wilmington, NC, documenting the directional wave spectrum in the region between Charleston, SC and Cape Hatteras, NC. The aircraft flight lines included segments near and along the shoreline as well as far offshore. Animations of the directional wave spectrum spatial variation along the aircraft tracks on the two flights will be presented using a 100: 1 time compression.

Walsh, Edward J.↗
NASA NTRS1983

A multigrid method for the Euler equations

A multigrid algorithm has been developed for the numerical solution of the steady two-dimensional Euler equations. Flux vector splitting and one-sided differencing are employed to define the spatial discretization. Newton's method is used to solve the nonlinear equations, and a multigrid solver is used on each linear problem. The relaxation scheme for the linear problems is symmetric Gauss-Seidel. Standard restriction and interpolation operators are employed. Local mode analysis is used to predict the convergence rate of the multigrid process on the linear problems. Computed results for transonic flows over airfoils are presented.

Jespersen, D. C.↗
DOE OSTI2017

The Curvature-Augmented Closest Point method with vesicle inextensibility application

Here, the Closest Point method, initially developed by Ruuth and Merriman, allows for the numerical solution of surface partial differential equations without the need for a parameterization of the surface itself. Surface quantities are embedded into the surrounding domain by assigning each value at a given spatial location to the corresponding value at the closest point on the surface. This embedding allows for surface derivatives to be replaced by their Cartesian counterparts (e.g. ∇ s =∇). This equivalence is only valid on the surface, and thus, interpolation is used to enforce what is known as the side condition away from the surface. To improve upon the method, this work derives an operator embedding that incorporates curvature information, making it valid in a neighborhood of the surface. With this, direct enforcement of the side condition is no longer needed. Comparisons in R 2 and R 3 show that the resulting Curvature-Augmented Closest Point method has better accuracy and requires less memory, through increased matrix sparsity, than the Closest Point method, while maintaining similar matrix condition numbers. To demonstrate the utility of the method in a physical application, simulations of inextensible, bi-lipid vesicles evolving toward equilibrium shapes are also included.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1984

Multigrid techniques for the numerical solution of the diffusion equation

An accurate numerical solution of diffusion problems containing large local gradients can be obtained with a significant reduction in computational time by using a multigrid computational scheme. The spatial domain is covered with sets of uniform square grids of different sizes. The finer grid patterns overlap the coarse grid patterns. The finite-difference expressions for each grid pattern are solved independently by iterative techniques. Two interpolation methods were used to establish the values of the potential function on the fine grid boundaries with information obtained from the coarse grid solution. The accuracy and computational requirements for solving a test problem by a simple multigrid and a multilevel-multigrid method were compared. The multilevel-multigrid method combined with a Taylor series interpolation scheme was found to be best.

Phillips, R. E.↗
NASA NTRS1995

Numerical solution of the full potential equation using a chimera grid approach

A numerical scheme utilizing a chimera zonal grid approach for solving the full potential equation in two spatial dimensions is described. Within each grid zone a fully-implicit approximate factorization scheme is used to advance the solution one interaction. This is followed by the explicit advance of all common zonal grid boundaries using a bilinear interpolation of the velocity potential. The presentation is highlighted with numerical results simulating the flow about a two-dimensional, nonlifting, circular cylinder. For this problem, the flow domain is divided into two parts: an inner portion covered by a polar grid and an outer portion covered by a Cartesian grid. Both incompressible and compressible (transonic) flow solutions are included. Comparisons made with an analytic solution as well as single grid results indicate that the chimera zonal grid approach is a viable technique for solving the full potential equation.

Holst, Terry L.↗
NASA NTRS2022

Tropospheric Ozone Retrieval By A Combination of TROPOMI/S5P Measurements With BASCOE Assimilated Data

We present a new tropospheric ozone dataset based on TROPOspheric Monitoring Instrument (TROPOMI)/Sentinel-5 Precursor (S5P) total ozone measurements combined with stratospheric ozone data from the Belgian Assimilation System for Chemical ObsErvations (BASCOE) constrained by assimilating ozone observations from the Microwave Limb Sounder (MLS). The BASCOE stratospheric data are interpolated to the S5P observations and subtracted from the TROPOMI total ozone data. The difference is equal to the tropospheric ozone residual column from the surface up to the tropopause. The tropospheric ozone columns are retrieved at the full spatial resolution of the TROPOMI sensor (5.5×3.5 km 2 ) with daily global coverage. Compared to the Ozone Mapping and Profiler Suite Modern-Era Retrospective analysis for Research and Applications 2 (OMPS-MERRA-2) data, a global mean positive bias of 3.3 DU is found for the analysed period April 2018 to June 2020. A small negative bias of about −0.91 DU is observed in the tropics relative to the operational TROPOMI tropical tropospheric data based on the convective cloud differential (CCD) algorithm throughout the same period. The new tropospheric ozone data (S5P-BASCOE) are compared to a set of globally distributed ozonesonde data integrated up to the tropopause level. We found 2254 comparisons with cloud-free TROPOMI observations within 25 km of the stations. In the global mean, S5P-BASCOE deviates by 2.6 DU from the integrated ozonesondes. Depending on the latitude the S5P-BASCOE deviate from the sondes and between −4.8 and 7.9 DU, indicating a good agreement. However, some exceptional larger positive deviations up to 12 DU are found, especially in the northern polar regions (north of 70∘). The monthly mean tropospheric column and time series for selected areas showed the expected spatial and temporal pattern, such as the wave one structure in the tropics or the seasonal cycle, including a summer maximum, in the mid-latitudes.

Satellite Retrieval↗
DOE OSTI2023

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part I - Analysis with a Small Sample Size

Computer simulations, especially of complex phenomena, can be expensive, requiring high-performance computing resources. Often, to understand a phenomenon, multiple simulations are run, each with a different set of simulation input parameters. These data are then used to create an interpolant, or surrogate, relating the simulation outputs to the corresponding inputs. When the inputs and outputs are scalars, a simple machine learning model can suffice. However, when the simulation outputs are vector valued, available at locations in two or three spatial dimensions, often with a temporal component, creating a surrogate is more challenging. In this report, we use a two-dimensional problem of a jet interacting with high explosives to understand how we can build high-quality surrogates. The characteristics of our data set are unique - the vector-valued outputs from each simulation are available at over two million spatial locations; each simulation is run for a relatively small number of time steps; the size of the computational domain varies with each simulation; and resource constraints limit the number of simulations we can run. We show how we analyze these extremely large data-sets, set the parameters for the algorithms used in the analysis, and use simple ways to improve the accuracy of the spatio-temporal surrogates without substantially increasing the number of simulations required.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

Spatiotemporal Interpolation of Elevation Changes Derived from Satellite Altimetry for Jakobshavn Isbrae, Greenland

Estimation of ice sheet mass balance from satellite altimetry requires interpolation of point-scale elevation change (dHdt) data over the area of interest. The largest dHdt values occur over narrow, fast-flowing outlet glaciers, where data coverage of current satellite altimetry is poorest. In those areas, straightforward interpolation of data is unlikely to reflect the true patterns of dHdt. Here, four interpolation methods are compared and evaluated over Jakobshavn Isbr, an outlet glacier for which widespread airborne validation data are available from NASAs Airborne Topographic Mapper (ATM). The four methods are ordinary kriging (OK), kriging with external drift (KED), where the spatial pattern of surface velocity is used as a proxy for that of dHdt, and their spatiotemporal equivalents (ST-OK and ST-KED).

Jakobshavn Isbrae↗
NASA NTRS2018

Super Resolution of Remote Sensing Images Using Edge-Directed Radial Basis Functions

Edge-Directed Radial Basis Functions (EDRBF) are used to compute super resolution(SR) image from a given set of low resolution (LR) images differing in subpixel shifts. The algorithm is tested on remote sensing images and compared for accuracy with other well-known algorithms such as Iterative Back Projection (IBP), Maximum Likelihood (ML) algorithm, interpolation of scattered points using Nearest Neighbor (NN) and Inversed Distance Weighted (IDW) interpolation, and Radial Basis Functin(RBF) . The accuracy of SR depends on various factors besides the algorithm (i) number of subpixel shifted LR images (ii) accuracy with which the LR shifts are estimated by registration algorithms (iii) and the targeted spatial resolution of SR. In our studies, the accuracy of EDRBF is compared with other algorithms keeping these factors constant. The algorithm has two steps: i) registration of low resolution images and (ii) estimating the pixels in High Resolution (HR) grid using EDRBF. Experiments are conducted by simulating LR images from a input HR image with different sub-pixel shifts. The reconstructed SR image is compared with input HR image to measure the accuracy of the algorithm using sum of squared errors (SSE). The algorithm has outperformed all of the algorithms mentioned above. The algorithm is robust and is not overly sensitive to the registration inaccuracies.

Mareboyana, Manohar M.↗
DOE OSTI2026

Systematic analysis of melt pool dynamics in laser processing of mixed powder feedstocks

Functionally graded materials (FGMs) fabricated via additive manufacturing of blended powders offer the potential to spatially tailor properties for new technologies, such as fusion first-wall systems, turbine blades, and spacecraft. However, processing these materials is difficult due to the multiplicity of processing parameters to optimize, all of which must be changed as substrate material, powder feedstock compositions, and melt pool dynamics evolve. Here, this work systematically evaluates the qualitative and quantitative effects of these variables on the melt pool size, shape, composition, and particle distribution in an exemplar Ti-Ta system, and connects the experimental results to Marangoni flow behavior and phenomena observed in other systems. Increasing laser power linearly increases melt pool size and layer thickness, driving engineering considerations such as part/geometrical tolerances. Decreasing laser velocity changes the melt pool shape from lenticular to convex and reduces chemical homogeneity due to extreme thermal and compositional gradients between the melt pool center and boundaries. Thermophysical property differences between the powder feedstock and substrate material, as well as the directionality of the gradient, affect dilution and melt pool dynamics, which in turn affect the melt pool boundary characteristics, shape, and uniformity. Mixed powder feedstocks of intermediate compositions do not behave according to linear interpolations between single-material endpoints, instead building taller and wider melt pools. As such, it is recommended to quantify process maps for at least one intermediate composition in the FGM or multi-material system of interest to ensure optimized processing parameters, predictable melt pool sizes and shapes, and compositional and spatial precision.

Dissimilar↗
DOE OSTI2023

Fitting Matérn smoothness parameters using automatic differentiation

The Mat$\acute{e}$rn covariance function is ubiquitous in the application of Gaussian processes to spatial statistics and beyond. Perhaps the most important reason for this is that the smoothness parameter $\nu$ gives complete control over the mean-square differentiability of the process, which has significant implications for the behavior of estimated quantities such as interpolants and forecasts. Unfortunately, derivatives of the Mat$\acute{e}$rn covariance function with respect to $\nu$ require derivatives of the modified second-kind Bessel function $K$ $\nu$ with respect to $\nu$. While closed form expressions of these derivatives do exist, they are prohibitively difficult and expensive to compute. For this reason, many software packages require fixing $\nu$ as opposed to estimating it, and all existing software packages that attempt to offer the functionality of estimating $\nu$ use finite difference estimates for $\partial$ $\nu$ $K$ $\nu$ . In this work, we introduce a new implementation of $K$$\nu$ that has been designed to provide derivatives via automatic differentiation (AD), and whose resulting derivatives are significantly faster and more accurate than those computed using finite differences. Here, we provide comprehensive testing for both speed and accuracy and show that our AD solution can be used to build accurate Hessian matrices for second-order maximum likelihood estimation in settings where Hessians built with finite difference approximations completely fail.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An a priori evaluation of a principal component and artificial neural network based combustion model in diesel engine conditions

A principal component analysis (PCA) and artificial neural network (ANN) based chemistry tabulation approach is presented. ANNs are used to map the thermochemical state onto a low-dimensional manifold consisting of five control variables that have been identified using PCA. Three canonical configurations are considered to train the PCA-ANN model: a series of homogeneous reactors, a nonpremixed flamelet, and a two-dimensional lifted flame. The performance of the model in predicting the thermochemical manifold of a spatially-developing turbulent jet flame in diesel engine thermochemical conditions is a priori evaluated using direct numerical simulation (DNS) data. The PCA-ANN approach is compared with a conventional tabulation approach (tabulation using ad hoc defined control variables and linear interpolation). The PCA-ANN model provides higher accuracy and requires several orders of magnitude less memory. Here, these observations indicate that the PCA-ANN model is superior for chemistry tabulation, especially for modelling complex chemistries that present multiple combustion modes as observed in diesel combustion. The performance of the PCA-ANN model is then compared to the optimal estimator, i.e. the conditional mean from the DNS. The results indicate that the PCA-ANN model gives high prediction accuracy, comparable to the optimal estimator, especially for major species and the thermophysical properties. Higher errors are observed for the minor species and reaction rate predictions when compared to the optimal estimator. It is shown that the prediction of minor species and reaction rates can be improved by using training data that exhibits a variation of parameters as observed in the turbulent flame. The output of the ANN is analysed to assess mass conservation. It is observed that the ANN incurs a mean absolute error of 0.05% in mass conservation. Furthermore, it is demonstrated that this error can be reduced by modifying the cost function of the ANN to penalise for deviation from mass conservation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Distributed Fiber Optics Measurements of Rock Deformation and Failure in Triaxial Tests

The implementation of distributed strain measurement methods in triaxial and uniaxial tests have demonstrated the development of strain localization, even at early stages. This implies that single point measurement methods are location dependent. The use of distributed methods is required not only to improve the interpreted constitutive parameters obtained from triaxial tests, but also to understand the implications of strain localization in the failure process. Here, this work uses optical fibers in triaxial tests. The developed distributed measurement method was implemented on granite, gabbro, and sandstone samples and tested under different confining pressures, reaching 200 MPa on the granite sample. Using a temporal resolution of 0.25 Hz and 5 mm of spatial resolution, the strain evolution at over 300 locations at the sample surface was measured during testing. When compared to point sensing methods, the use of optical fiber greatly increases the number of measurements at the surface of the sample, and their interpolation provides the entire deformation of the sample surface. In all the tests performed, strain localizations were revealed before failure. A three-dimensional interpretation of a test, combining an optical scan of the sample and the distributed measurements, show good correlation between the fractures and the strain localization.

58 GEOSCIENCES↗