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

A comparison of different continuum approaches in modeling mixed-type dislocations in Al

Mixed-type dislocations are common in metals and play a vital role in their plastic deformation. Key characteristics of mixed-type dislocations cannot simply be extrapolated from those of dislocations with pure edge or pure screw characters. Yet, mixed-type dislocations traditionally received disproportionately less attention in the modeling and simulation community. Here, we explore core structures of mixed-type dislocations in Al using three continuum approaches, namely, the phase-field dislocation dynamics (PFDD) method, the atomistic phase-field microelasticity (APFM) method, and the concurrent atomistic-continuum (CAC) method. Results are benchmarked against molecular statics. We advance the PFDD and APFM methods in several aspects such that they can better describe the dislocation core structure. In particular, in these two approaches, the gradient energy coefficients for mixed-type dislocations are determined based on those for pure-type ones using a trigonometric interpolation scheme, which is shown to provide better prediction than a linear interpolation scheme. The dependence of the in-slip-plane spatial numerical resolution in PFDD and CAC is also quantified.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Antarctic sea ice concentrations inferred from Nimbus 5 ESMR and Landsat imagery

Landsat Multispectral Scanner near-infrared band images (MSS-7) are used to calculate Antarctic ice concentration values and these results are quantitatively compared to those derived from Nimbus 5 electrically scanning microwave radiometer (ESMR) brightness temperature data. The set of images includes regions near the ice edge and near the continental boundary during late winter or spring. The percentage of ice cover is derived from the microwave data, using an algorithm that incorporates the brightness temperature emissivity and interpolated climatological physical temperatures and that is not dependent on spatial resolution. Cloud-free Landsat images are used to provide an alternative determination of large-scale ice concentrations, but this method is limited in its ability to resolve individual ice flows, especially near the ice edge. A proportional classification procedure is applied to Landsat data based on observed reflection and the results are correlated with similarly-derived ESMR ice concentration values. The corresponding ice concentration values are found to agree to + or - 15%, with the largest uncertainties occurring in regions of low ice concentration.

Comiso, J. C.↗
NASA NTRS2000

BOREAS AFM-5 Level-2 Upper Air Network Standard Pressure Level Data

The BOREAS AFM-5 team collected and processed data from the numerous radiosonde flights during the project. The goals of the AFM-05 team were to provide large-scale definition of the atmosphere by supplementing the existing AES aerological network, both temporally and spatially. This data set includes basic upper-air parameters interpolated at 0.5 kiloPascal increments of atmospheric pressure from data collected from the network of upper-air stations during the 1993, 1994, and 1996 field campaigns over the entire study region. The data are contained in tabular ASCII files. The data files are available on a CD-ROM (see document number 20010000884) or from the Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center (DAAC).

Barr, Alan↗
NASA NTRS2010

Sources of Uncertainty in Predicting Land Surface Fluxes Using Diverse Data and Models

In the domain of predicting land surface fluxes, models are used to bring data from large observation networks and satellite remote sensing together to make predictions about present and future states of the Earth. Characterizing the uncertainty about such predictions is a complex process and one that is not yet fully understood. Uncertainty exists about initialization, measurement and interpolation of input variables; model parameters; model structure; and mixed spatial and temporal supports. Multiple models or structures often exist to describe the same processes. Uncertainty about structure is currently addressed by running an ensemble of different models and examining the distribution of model outputs. To illustrate structural uncertainty, a multi-model ensemble experiment we have been conducting using the Terrestrial Observation and Prediction System (TOPS) will be discussed. TOPS uses public versions of process-based ecosystem models that use satellite-derived inputs along with surface climate data and land surface characterization to produce predictions of ecosystem fluxes including gross and net primary production and net ecosystem exchange. Using the TOPS framework, we have explored the uncertainty arising from the application of models with different assumptions, structures, parameters, and variable definitions. With a small number of models, this only begins to capture the range of possible spatial fields of ecosystem fluxes. Few attempts have been made to systematically address the components of uncertainty in such a framework. We discuss the characterization of uncertainty for this approach including both quantifiable and poorly known aspects.

Dungan, Jennifer L.↗
DOE OSTI2024

Estimation and Visualization of Isosurface Uncertainty from Linear and High-Order Interpolation Methods

Isosurface visualization is fundamental for exploring and analyzing 3D volumetric data. Marching cubes (MC) algorithms with linear interpolation are commonly used for isosurface extraction and visualization. Although linear interpolation is easy to implement, it has limitations when the underlying data is complex and high-order, which is the case for most real-world data. Linear interpolation can output vertices at the wrong location. Its inability to deal with sharp features and features smaller than grid cells can lead to an incorrect isosurface with holes and broken pieces. Despite these limitations, isosurface visualizations typically do not include insight into the spatial location and the magnitude of these errors. We utilize high-order interpolation methods with MC algorithms and interactive visualization to highlight these uncertainties. Our visualization tool helps identify the regions of high interpolation errors. It also allows users to query local areas for details and compare the differences between isosurfaces from different interpolation methods. In addition, we employ high-order methods to identify and reconstruct possible features that linear methods cannot detect. We showcase how our visualization tool helps explore and understand the extracted isosurface errors through synthetic and real-world data.

Ouermi, Timbwaoga↗
DOE OSTI2018

Spatial associations between social groups and ozone air pollution exposure in the Beijing urban area

Highlights: • Greater percentages of persons living alone and migrants were associated with higher O{sub 3} exposure. • Social groups divided by Hukou status had distinct spatial distributions in the Beijing urban area. • Migrants bore a double burden of social deprivation and O{sub 3} exposure due to city development planning and the Hukou system. • Areas with greater percentages of seniors had lower O{sub 3} exposure but higher NO{sub 2} exposure. Few studies have linked social factors to air pollution exposure in China. Unlike the race or minority concepts in western countries, the Hukou system (residential registration system) is a fundamental reason for the existence of social deprivation in China. To assess the differences in ozone (O{sub 3}) exposure among social groups, especially groups divided by Hukou status, we assigned estimates of O{sub 3} exposure to the latest census data of the Beijing urban area using a kriging interpolation model. We developed simultaneous autoregressive (SAR) models that account for spatial autocorrelation to identify the associations between O{sub 3} exposure and social factors. Principal component regression was used to control the multicollinearity bias as well as explore the spatial structure of the social data. The census tracts (CTs) with higher proportions of persons living alone and migrants with non-local Hukou were characterized by greater exposure to ambient O{sub 3}. The areas with greater proportions of seniors had lower O{sub 3} exposure. The spatial distribution patterns were similar among variables including migrants, agricultural population and household separation (population status with separation between Hukou and actual residences), which fit the demographic characteristics of the majority of migrants. Migrants bore a double burden of social deprivation and O{sub 3} pollution exposure due to city development planning and the Hukou system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

An Immersed Boundary method with divergence-free velocity interpolation and force spreading

The Immersed Boundary (IB) method is a mathematical framework for constructing robust numerical methods to study fluid–structure interaction in problems involving an elastic structure immersed in a viscous fluid. The IB formulation uses an Eulerian representation of the fluid and a Lagrangian representation of the structure. The Lagrangian and Eulerian frames are coupled by integral transforms with delta function kernels. The discretized IB equations use approximations to these transforms with regularized delta function kernels to interpolate the fluid velocity to the structure, and to spread structural forces to the fluid. It is well-known that the conventional IB method can suffer from poor volume conservation since the interpolated Lagrangian velocity field is not generally divergence-free, and so this can cause spurious volume changes. In practice, the lack of volume conservation is especially pronounced for cases where there are large pressure differences across thin structural boundaries. The aim of this paper is to greatly reduce the volume error of the IB method by introducing velocity-interpolation and force-spreading schemes with the properties that the interpolated velocity field in which the structure moves is at least $\mathscr{C}$ 1 and satisfies a continuous divergence-free condition, and that the force-spreading operator is the adjoint of the velocity-interpolation operator. We confirm through numerical experiments in two and three spatial dimensions that this new IB method is able to achieve substantial improvement in volume conservation compared to other existing IB methods, at the expense of a modest increase in the computational cost. Further, the new method provides smoother Lagrangian forces (tractions) than traditional IB methods. Finally, the method presented here is restricted to periodic computational domains. Its generalization to non-periodic domains is important future work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Variable Grid Method

NETL’s Variable Grid Method (VGM) is a novel approach that leverages GIS capabilities to simultaneously visualize and quantify spatial data trends and underlying data uncertainty. The flexible VGM approach utilizes a range of spatial datasets and uncertainty quantifications, which can be calculated using data related to sample density, sample variance, interpolation error, multiple simulations, etc., to create an integrated visualization of data and uncertainty. The intuitive manner of the VGM helps communicate the relationship between uncertainty and spatial data to effectively guide research, support advanced computation analyses, and helps inform research, management and policy decisions. The VGM approach is being developed into a standalone tool designed to integrate into a variety of research workflows, providing users the capabilities to utilize the VGM for their datasets, analysis, and models in support of their decision making needs. The VGM approach was developed as part of NETL’s Offshore portfolio. It was further matured and utilized in projects associated with both NETL’s Unconventional Resources and Carbon Storage portfolios, and is now being further developed by Varigrid, LLC.

Carbon Storage↗
DOE OSTI2018

Improving the spatial dynamic range of digital inline particle holography

Digital in-line holography (DIH) has been proven to provide three-dimensional droplet position, size, and velocity distributions with a single-camera. This data is crucial for understanding multi-phase flows. In this work, we examine the performance of this diagnostic in the limit of very small particles, on the order of a pixel in diameter and smaller, and propose a post-processing method to improve them: Lanczos interpolation. The Lanczos interpolation kernel is the digital implementation of the Whittaker sinc filter, and strikes a compromise between maintaining the spatial frequency ceiling of the original digital image and computational cost of the interpolation. Without Lanczos interpolation, or super-sampling, the ultimate detectable particle size floor is on the order of 4 pixel widths. We show in this work that this limit can be reduced by 50% or more with super-sampling, depending upon the desired diameter accuracy, and examine the effect of super-sampling on the resulting accuracy of the extracted size and position of spherical particles. Extending this resolution limit increases the overall detection efficiency of the diagnostic. Since this increases the spatial dynamic range of the diagnostic, it can also allow a larger field-of-view to be captured with the same particle size floor.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Validation of finite-element models using full-field experimental data: Levelling finite-element analysis data through a digital image correlation engine

Full-field data from digital image correlation (DIC) provide rich information for finite-element analysis (FEA) validation. However, there are several inherent inconsistencies between FEA and DIC data that must be rectified before meaningful, quantitative comparisons can be made, including strain formulations, coordinate systems, data locations, strain calculation algorithms, spatial resolutions and data filtering. As such, in this paper, we investigate two full-field validation approaches: (1) the direct interpolation approach, which addresses the first three inconsistencies by interpolating the quantity of interest from one mesh to the other, and (2) the proposed DIC-levelling approach, which addresses all six inconsistencies simultaneously by processing the FEA data through a stereo-DIC simulator to ‘level’ the FEA data to the DIC data in a regularisation sense. Synthetic ‘experimental’ DIC data were generated based on a reference FEA of an exemplar test specimen. The direct interpolation approach was applied, and significant strain errorswere computed, even though therewas no model form error, because the filtering effect of theDIC enginewas neglected. In contrast, the levelling approach provided accurate validation results, with no strain error when no model form error was present. Next, model form error was purposefully introduced via a mismatch of boundary conditions. With the direct interpolation approach, the mismatch in boundary conditions was completely obfuscated, while with the levelling approach, it was clearly observed. Finally, the ‘experimental’ DIC datawere purposefully misaligned slightly fromthe FEA data. Both validation techniques suffered from the misalignment, thus motivating continued efforts to develop a robust alignment process. In conclusion, direct interpolation is insufficient, and the proposed levelling approach is required to ensure that the FEA and the DIC data have the same spatial resolution and data filtering. Only after the FEA data have been ‘levelled’ to the DIC data can meaningful, quantitative error maps be computed.

42 ENGINEERING↗
DOE OSTI2020

Benzene homologues contaminants in a former herbicide factory site: distribution, attenuation, risk, and remediation implication

Benzene homologues often used as organic raw materials or as detergents in chemical industry are prone to accidental release into the environment which can cause serious long-term soil pollutions. In a large former herbicide factory site, we investigated 43 locations for benzene homologues contaminations in soil, soil gas, and groundwater; and studied the hydrogeological conditions. An inverse distance weighted interpolation method was employed to determine the pollutants three-dimensional spatial distribution in the soils. Results showed that benzene homologues residues were mainly originated from the herbicide production workshop; and that the pollution had horizontally expanded at the deeper soil layer. Contaminants had already migrated 15m downward from ground surface. Contaminant phase distribution study showed that NAPL was the primary phase (>99%) for the pollutants accumulated in the unsaturated zone; while it had not migrated to groundwater. The primary mechanism for contaminant transport and attenuation included dissolution of “occluded” NAPL into pore water and pollutant volatilization into soil pore space. Risk assessment revealed that the pollutants brought unacceptable high carcinogenic and non-carcinogenic risks to public health. In order to convert this former chemical processing factory site into a residential area, a remediation to the polluted production workshop sites is urgently required.

Yang, Shuo↗
DOE OSTI2022

Nanoindentation mapping defects filtration for heterogeneous materials using generative adversarial networks

Advanced composite materials with multiple phases and heterogeneous microstructure necessitate spatial mapping characterization of elastic modulus to develop constitutive relations and overall mechanical response. Such modulus mapping can be obtained using the nanoindentation technique, where the indenter tip raster over the selected microstructure region. Typically, a surface preparation procedure is done in the specimens to ensure proper contact between the indenter tip and sample surface. However, a near-perfect surface finish is unachievable in heterogeneous materials, primarily with ceramic reinforcements, due to the differential material removal rate during polishing. Thus, the nanoindenter records localized erroneous measurements due to differences in surface roughness and corresponding force response. This study establishes a novel deep learning-based strategy to rectify incorrect experimental spatial measurements acquire during nanoindentation modulus mapping. Here, the integrated bicubic interpolation and generative adversarial networks (GANs) model was trained using 14 ceramic and 18 metallic data sets, each comprising 65,536 measurements. The developed algorithm was validated against experimental measurements on four unknown specimens. The standard deviation in measured elastic modulus reduces by ~50% in ceramics and ~72% in metallic samples. This computational framework proposes a novel approach to reducing uncertainty in materials’ properties using state-of-the-art computer vision techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Scalable Volume Visualization for Big Scientific Data Modeled by Functional Approximation

Considering the challenges posed by the space and time complexities in handling extensive scientific volumetric data, various data representations have been developed for the analysis of large-scale scientific data. Multivariate functional approximation (MFA) is an innovative data model designed to tackle substantial challenges in scientific data analysis. It computes values and derivatives with high-order accuracy throughout the spatial domain, mitigating artifacts associated with zero- or first-order interpolation. However, the slow query time through MFA makes it less suitable for interactively visualizing a large MFA model. In this work, we develop the first scalable interactive volume visualization pipeline, MFA-DVV, for the MFA model encoded from large-scale datasets. Our method achieves low input latency through distributed architecture, and its performance can be further enhanced by utilizing a compressed MFA model while still maintaining a high-quality rendering result for scientific datasets. We conduct comprehensive experiments to show that MFA-DVV can decrease the input latency and achieve superior visualization results for big scientific data compared with existing approaches.

big scientific dataset↗
NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗
NASA NTRS1992

Numerical simulation of aerothermal loads in hypersonic engine inlets due to shock impingement

The effect of shock impingement on an axial corner simulating the inlet of a hypersonic vehicle engine is modeled using a finite-difference procedure. A three-dimensional dynamic grid adaptation procedure is utilized to move the grids to regions with strong flow gradients. The adaptation procedure uses a grid relocation stencil that is valid at both the interior and boundary points of the finite-difference grid. A linear combination of spatial derivatives of specific flow variables, calculated with finite-element interpolation functions, are used as adaptation measures. This computational procedure is used to study laminar and turbulent Mach 6 flows in the axial corner. The description of flow physics and qualitative measures of heat transfer distributions on cowl and strut surfaces obtained from the analysis are compared with experimental observations. Conclusions are drawn regarding the capability of the numerical scheme for enhanced modeling of high-speed compressible flows.

Ramakrishnan, R.↗
NASA NTRS1993

Dynamic-overlapped-grid simulation of aerodynamically determined relative motion

Currently, there is a need to develop a means of analyzing and studying unsteady flowfields which involve multiple component configurations with at least one of the components in relative motion with respect to the others. Two of the important phenomena that such analyses can help to understand are the unsteady aerodynamic interference and the boundary-induced component of the flowfield. With this motivation, a computational method is developed which couples the governing equations of the unsteady flowfield and the rigid-body dynamics in six degrees-of-freedom. These equations are solved on composite meshes of overlapped subdomain grids which can move with respect to each other. Initially, several measures that reduce the numerical error are studied and compared with the exact solution of a moving normal shock in a tube. It is concluded that a second-order accurate method, for spatial and temporal discretizations as well as for the moving subdomain interpolations, is needed as a minimum measure. Furthermore, the CFL numbers should be restricted to unity. Then, the method is used to simulate the flowfield history and predict the aerodynamically determined trajectory of a store dropped from its initial position under a wing.

Yen, Guan-Wei↗
NASA NTRS1993

Video Animation of Ocean Topography From TOPEX/POSEIDON

Three video loops showing various aspects of the dynamic ocean topography obtained from the TOPEX/POSEIDON radar altimetry data will be presented. The first shows the temporal change of the global ocean topography during the first year of the mission. The time-averaged mean is removed to reveal the temporal variabilities. Temporal interpolation is performed to create daily maps for the animation. A spatial smoothing is also performed to retain only the large-sale features. Gyre-scale seasonal changes are the main features. The second shows the temporal evolution of the Gulf Stream. The high resolution gravimetric geoid of Rapp is used to obtain the absolute ocean topography. Simulated drifters are used to visualize the flow pattern of the current. Meanders and rings of the current are the main features. The third is an animation of the global ocean topography on a spherical earth. The JGM-2 geoid is used to obtain the ocean topography...

Fu, Lee-Lueng↗
DOE OSTI2024

Efficient continuous Energy-Multigroup hybrid depletion scheme using the Shift Monte Carlo code. Part I: Energy condensation sensitivity analysis

Monte Carlo (MC) codes coupled to depletion solvers are increasingly used to provide high fidelity fuel cycle modeling capabilities. Here, these coupled depletion-MC tools produce accurate results in general but can experience nonphysical spatial oscillations when time steps are large or when a system’s dominance ratio approaches unity. Two substepping techniques have been developed previously to remedy and dampen these spatial oscillations without needing to reduce step sizes. The first approach relied on higher-order techniques to account for spectral changes within steps (extrapolation and interpolation techniques). The second approach used the first order perturbation (FOP) theory to account for the change in the one-group spatial flux distribution within steps. This paper develops a hybrid depletion methodology which, in a way, combines how the flux is handled in both substepping techniques. Specifically, the multigroup (MG) MC Shift code is used to update the flux distribution within steps rather than a one-group FOP solver. A fully reflected pincell is investigated, which is not spatially dependent in the MG representation. Thus, the analysis in this paper is an initial demonstration of hybrid depletion. An upcoming companion paper will focus on how the hybrid depletion dampens spatial oscillations. The hybrid depletion approach is verified to be consistent with previous constant extrapolation depletion (CED) methods. This paper finds that the hybrid CED exhibits some error in the eigenvalue and one group constants within macro steps. To address this discrepancy, a simple interpolation scheme (CELI) is investigated. This work found that CELI sufficiently addresses the discrepancy in spectrum for macro steps up to 100 days. Overall, this work demonstrates that the hybrid depletion method can significantly reduce the number of high fidelity MC executions in a MC-coupled depletion with an acceptable eigenvalue error.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗