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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The concentration dependence of the galaxy–halo connection: Modeling assembly bias with abundance matching

Empirical methods for connecting galaxies to their dark matter halos have become essential for interpreting measurements of the spatial statistics of galaxies. In this work, we present a novel approach for parameterizing the degree of concentration dependence in the abundance matching method. Furthermore, this new parameterization provides a smooth interpolation between two commonly used matching proxies: the peak halo mass and the peak halo maximal circular velocity. This parameterization controls the amount of dependence of galaxy luminosity on halo concentration at a fixed halo mass. Effectively this interpolation scheme enables abundance matching models to have adjustable assembly bias in the resulting galaxy catalogs. With the new $400\,\mathrm{Mpc}\,{h}^{-1}$ DarkSky Simulation, whose larger volume provides lower sample variance, we further show that low-redshift two-point clustering and satellite fraction measurements from SDSS can already provide a joint constraint on this concentration dependence and the scatter within the abundance matching framework.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2010

Phenological Parameters Estimation Tool

The Phenological Parameters Estimation Tool (PPET) is a set of algorithms implemented in MATLAB that estimates key vegetative phenological parameters. For a given year, the PPET software package takes in temporally processed vegetation index data (3D spatio-temporal arrays) generated by the time series product tool (TSPT) and outputs spatial grids (2D arrays) of vegetation phenological parameters. As a precursor to PPET, the TSPT uses quality information for each pixel of each date to remove bad or suspect data, and then interpolates and digitally fills data voids in the time series to produce a continuous, smoothed vegetation index product. During processing, the TSPT displays NDVI (Normalized Difference Vegetation Index) time series plots and images from the temporally processed pixels. Both the TSPT and PPET currently use moderate resolution imaging spectroradiometer (MODIS) satellite multispectral data as a default, but each software package is modifiable and could be used with any high-temporal-rate remote sensing data collection system that is capable of producing vegetation indices. Raw MODIS data from the Aqua and Terra satellites is processed using the TSPT to generate a filtered time series data product. The PPET then uses the TSPT output to generate phenological parameters for desired locations. PPET output data tiles are mosaicked into a Conterminous United States (CONUS) data layer using ERDAS IMAGINE, or equivalent software package. Mosaics of the vegetation phenology data products are then reprojected to the desired map projection using ERDAS IMAGINE

McKellip, Rodney D.↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

Hurricane Directional Wave Spectrum Spatial Variation in the Open Ocean

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 deg 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 data presented were acquired on 24 August 1998 when hurricane Bonnie was east of the Bahamas and moving slowly to the north. Wave 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 and at times there were wave fields traveling at right angles to each other. The NOAA aircraft spent over five hours within 180 km of the hurricane Bonnie eye, and made five eye penetrations. A 2-minute animation of the directional wave spectrum spatial variation over this period will be shown.

Wright, C. W.↗
DOE OSTI2014

SU-E-J-08: A Hybrid Three Dimensional Registration Framework for Image-Guided Accurate Radiotherapy System ARTS-IGRT

Purpose: The purpose of this work was to develop a registration framework and method based on the software platform of ARTS-IGRT and implement in C++ based on ITK libraries to register CT images and CBCT images. ARTS-IGRT was a part of our self-developed accurate radiation planning system ARTS. Methods: Mutual information (MI) registration treated each voxel equally. Actually, different voxels even having same intensity should be treated differently in the registration procedure. According to their importance values calculated from self-information, a similarity measure was proposed which combined the spatial importance of a voxel with MI (S-MI). For lung registration, Firstly, a global alignment method was adopted to minimize the margin error and achieve the alignment of these two images on the whole. The result obtained at the low resolution level was then interpolated to become the initial conditions for the higher resolution computation. Secondly, a new similarity measurement S-MI was established to quantify how close the two input image volumes were to each other. Finally, Demons model was applied to compute the deformable map. Results: Registration tools were tested for head-neck and lung images and the average region was 128*128*49. The rigid registration took approximately 2 min and converged 10% faster than traditional MI algorithm, the accuracy reached 1mm for head-neck images. For lung images, the improved symmetric Demons registration process was completed in an average of 5 min using a 2.4GHz dual core CPU. Conclusion: A registration framework was developed to correct patient's setup according to register the planning CT volume data and the daily reconstructed 3D CBCT data. The experiments showed that the spatial MI algorithm can be adopted for head-neck images. The improved Demons deformable registration was more suitable to lung images, and rigid alignment should be applied before deformable registration to get more accurate result. Supported by National Natural Science Foundation of China (NO.81101132) and Natural Science Foundation of Anhui Province (NO.11040606Q55)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Real-Space Quantification of Exciton Localization in Acene Crystals Using Wannier Function Decomposition

Here, we introduce the Wannier function decomposition of excitons (WFDX) method to quantify exciton localization in solids within the ab initio Bethe-Salpeter equation framework. By decomposing each Bloch exciton wave function into products of single-particle electron and hole maximally localized Wannier functions, this real-space approach provides well-defined orbital- and spatial-resolved measures of both Frenkel and charge-transfer excitons at low computational cost. We apply WFDX to excitons in acene crystals, quantifying how the number of rings, the exciton spin state, and the center-of-mass momentum affect spatial localization. Additionally, we show how this real-space representation reflects structural nonsymmorphic symmetries that are hidden in standard reciprocal-space descriptions. We demonstrate how the WFDX framework can be used to efficiently interpolate exciton expansion coefficients in reciprocal-space and outline how it may facilitate evaluation of observables involving position operators, highlighting its potential as a general tool for both analyzing and computing excitonic properties in solids.

Tao, Zui [University of California, Berkeley, CA (↗
DOE OSTI2025

High-resolution national mapping of natural gas composition substantially updates methane leakage impacts

Methane is emitted from oil and gas operations alongside heavier hydrocarbons and non-hydrocarbon gases, shaping emissions management decision-making, including air quality impacts. Yet, most assessments assume fixed gas composition, overlooking significant spatial and temporal variations. Here, we generate a high-resolution, data-driven map of natural gas composition across the United States, reconstructing methane, heavier hydrocarbons, and non-hydrocarbon species using spatio-temporal interpolation and oil-and-gas production patterns. Our approach is able to reduce composition prediction errors by 39% in terms of Mean Absolute Error (MAE) compared to standard techniques and reveals that methane loss rates have been underestimated by more than 50% in some regions. Beyond methane, we uncover substantial variability in co-emitted gases, exposing blind spots in current emissions inventories and emissions management frameworks. Our work enables more accurate emissions assessments, guides targeted measurement strategies, and informs emissions management decision-making. It also provides a general framework for prediction in environmental applications that integrate sparse measurements with auxiliary variables.

03 NATURAL GAS↗
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.↗
DOE OSTI2014

AN ALGORITHM FOR RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT TRANSFER EQUATION

We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics in the regime that light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity-dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the ATHENA code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

In situ compressive damage of cement paste characterized by lab source X-ray computer tomography

This paper aims at illustrating the potential of lab source X-ray CT for studying the damage behavior of cement based materials through in situ load experiments. This approach permits quantifying the microstructure prior and during loading. The load damage is separated from the specimen deformation using an image interpolation method. A quantitative relationship between external load and internal specimen damage is analyzed using the statistical information of gray scale values of the CT data. Local damage degrees are defined on 3D subset, and the 3D spatial distribution of damage information is clarified in this research. - Highlights: • On line damage is characterized by lab source X-ray CT. • Loading damage is separated with the specimen deformation. • Local damage is analyzed using gray scale values of the CT data. • 3D spatial distribution of the local damage information is clarified.

36 MATERIALS SCIENCE↗
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.↗
DOE OSTI2014

CHROMOSPHERIC AND CORONAL OBSERVATIONS OF SOLAR FLARES WITH THE HELIOSEISMIC AND MAGNETIC IMAGER

We report observations of white-light ejecta in the low corona, for two X-class flares on 2013 May 13, using data from the Helioseismic and Magnetic Imager (HMI) of the Solar Dynamics Observatory. At least two distinct kinds of sources appeared (chromospheric and coronal), in the early and later phases of flare development, in addition to the white-light footpoint sources commonly observed in the lower atmosphere. The gradual emissions have a clear identification with the classical loop-prominence system, but are brighter than expected and possibly seen here in the continuum rather than line emission. We find the HMI flux exceeds the radio/X-ray interpolation of the bremsstrahlung produced in the flare soft X-ray sources by at least one order of magnitude. This implies the participation of cooler sources that can produce free-bound continua and possibly line emission detectable by HMI. One of the early sources dynamically resembles {sup c}oronal rain{sup ,} appearing at a maximum apparent height and moving toward the photosphere at an apparent constant projected speed of 134 ± 8 km s{sup –1}. Not much literature exists on the detection of optical continuum sources above the limb of the Sun by non-coronagraphic instruments and these observations have potential implications for our basic understanding of flare development, since visible observations can in principle provide high spatial and temporal resolution.

79 ASTRONOMY AND ASTROPHYSICS↗