Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial statistics”

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 325 records · Page 18

Single-Filament Imaging Mass Spectrometry Lipidomics in Arthrospira platensis

Elucidating intra-organismal biochemical and lipid organization in photosynthetic biological cell factories of filamentous cyanobacteria, i.e., Arthrospira platensis (Spirulina), is important to track physiological response mechanisms during growth. Little is known about the filaments' biochemical organization and cellular structure and no label-free imaging techniques exist that provide molecular mapping. We applied ultra-high resolution mass spectrometry (7T FT-ICR-MS) matrix-assisted laser desorption ionization (MALDI) imaging to immobilized Spirulina filaments to investigate the localization of lipids across distinct physiological regions. We optimized matrix selection and deposition methods with the goal of facilitating high spatial, and intra-filament, resolution using untargeted multivariate statistical spectral deconvolution across MS pixels. Our results demonstrate an improved two step application with an optimized procedure for intra-organismal lipid profiling to improve analyte sensitivity and achieve higher spatial resolution, whereby we evaluate three conventional matrices 2,5-dihydroxybenzoic acid (DHB), a 9:1 ratio of DHB:superDHB (sDHB), 1,5- diaminonaphthalene (DAN) and a 50:50 mix of DHB:sDHB and compare delineation and pixel-based elucidation of intra-filament lipidomics. We identified a total of 1,626 features that could be putatively assigned a lipid-like formula based on database query and 46 unique features, with associated lipid assignments that were significantly distinct in their intra-filament location. MALDI- imaging MS with untargeted statistical spectral deconvolution was used to visualize intra-filament lipidomics organization in Spirulina filaments. Improvements in matrix deposition, including sequential sublimation and pneumatic spraying, increased signal abundance at high spatial resolution and allowed for identification of distinct lipid composition regions. Here in this work outlines a methodology that may be used for micro-ecological untargeted molecular phenotyping.

09 BIOMASS FUELS↗

Scanning Backscatter Lidar Observations for Characterizing 4-D Cloud and Aerosol Fields to Improve Radiative Transfer Parameterizations

Clouds have a powerful influence on atmospheric radiative transfer and hence are crucial to understanding and interpreting the exchange of radiation between the Earth's surface, the atmosphere, and space. Because clouds are highly variable in space, time and physical makeup, it is important to be able to observe them in three dimensions (3-D) with sufficient resolution that the data can be used to generate and validate parameterizations of cloud fields at the resolution scale of global climate models (GCMs). Simulation of photon transport in three dimensionally inhomogeneous cloud fields show that spatial inhomogeneities tend to decrease cloud reflection and absorption and increase direct and diffuse transmission, Therefore it is an important task to characterize cloud spatial structures in three dimensions on the scale of GCM grid elements. In order to validate cloud parameterizations that represent the ensemble, or mean and variance of cloud properties within a GCM grid element, measurements of the parameters must be obtained on a much finer scale so that the statistics on those measurements are truly representative. High spatial sampling resolution is required, on the order of 1 km or less. Since the radiation fields respond almost instantaneously to changes in the cloud field, and clouds changes occur on scales of seconds and less when viewed on scales of approximately 100m, the temporal resolution of cloud properties should be measured and characterized on second time scales. GCM time steps are typically on the order of an hour, but in order to obtain sufficient statistical representations of cloud properties in the parameterizations that are used as model inputs, averaged values of cloud properties should be calculated on time scales on the order of 10-100 s. The Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE) provides exceptional temporal (100 ms) and spatial (30 m) resolution measurements of aerosol and cloud backscatter in three dimensions. HARLIE was used in a ground-based configuration in several recent field campaigns. Principal data products include aerosol backscatter profiles, boundary layer heights, entrainment zone thickness, cloud fraction as a function of altitude and horizontal wind vector profiles based on correlating the motions of clouds and aerosol structures across portions of the scan. Comparisons will be made between various cloud detecting instruments to develop a baseline performance metric.

Schwemmer, Geary K.↗

Interdisciplinary research on the application of ERTS-1 data to the regional land use planning process.

Although the degree to which ERTS-1 imagery can satisfy regional land use planning data needs is not yet known, it appears to offer means by which the data acquisition process can be immeasurably improved. This paper documents the initial experiences of an interdisciplinary group attempting to formulate ways of analyzing the effectiveness of ERTS-1 imagery as a base for environmental monitoring and the resolution of regional land allocation problems. Because of the need to describe and depict regional resource complexity in an interrelatable state, certain resources within the geographical regions have been inventoried and stored in a two-dimensional computer-based map form. Computer oriented processes were developed to provide for the economical storage, analysis and spatial display of natural and cultural data for regional land use planning purposes. Statistical programs have been developed that correlate interpreted data with stored data, both spatially and numerically.

Clapp, J. L.↗

Incorporation of spatial information in Bayesian image reconstruction - The maximum residual likelihood criterion

We have developed a new figure of merit, a 'maximum-residual-likelihood' (MRL) statistic, for the goodness of fit for Bayesian image restoration which explicitly incorporates spatial information. The MRL constraint provides a natural means of incorporating the prior knowledge that the residuals contain no spatial structure through the autocorrelation function of the residuals. We demonstrate that this statistic follows a Chi-square distribution and that forcing this statistic to have its most probable value leads to a restored image whose residuals are consistent with the noise model. Our numerical experiments suggest that image restoration using the MRL statistic alone is numerically robust and produces results which are independent of the initial guess for the restored image. However, we caution that using the MRL statistic without an image prior can result in overresolution in low SNR portions of the image.

Pina, R. K.↗

Time-resolved particle image velocimetry

Time-resolved particle image velocimetry (TR-PIV) has become widespread in fluid dynamics. Essentially a velocity field movie, the dynamic content provides temporal as well as spatial information, in contrast to conventional PIV offering only statistical ensembles of flow quantities. From these time series arise further analyses such as accelerometry, space-time correlations, frequency spectra of turbulence including spatial variability, and derivation of pressure fields and forces. The historical development of TR-PIV is chronicled, culminating in an assessment of the current state of technology in high-repetition-rate lasers and high-speed cameras. Commercialization of pulse-burst lasers has expanded TR-PIV into more flows, including the compressible regime, and has achieved MHz rates. Particle response times and peak locking during image interrogation require attention but generally are not impediments to success. Accuracy considerations are discussed, including the risks of noise and aliasing in spectral content. Oversampled TR-PIV measurements allow use of multi-frame image interrogation methods, which improve the precision of the correlation and raise the velocity dynamic range of PIV. In combination with volumetric methods and data assimilation, a full four-dimensional description of a flow is not only achievable but becoming standardized. Lastly, a survey of exemplary applications is followed by a few predictions concerning the future of TR-PIV.

47 OTHER INSTRUMENTATION↗

Is the spatial distribution of brain lesions associated with closed-head injury predictive of subsequent development of attention-deficit/hyperactivity disorder? Analysis with brain-image database

PURPOSE: To determine whether there is an association between the spatial distribution of lesions detected at magnetic resonance (MR) imaging of the brain in children after closed-head injury and the development of secondary attention-deficit/hyperactivity disorder (ADHD). MATERIALS AND METHODS: Data obtained from 76 children without prior history of ADHD were analyzed. MR images were obtained 3 months after closed-head injury. After manual delineation of lesions, images were registered to the Talairach coordinate system. For each subject, registered images and secondary ADHD status were integrated into a brain-image database, which contains depiction (visualization) and statistical analysis software. Using this database, we assessed visually the spatial distributions of lesions and performed statistical analysis of image and clinical variables. RESULTS: Of the 76 children, 15 developed secondary ADHD. Depiction of the data suggested that children who developed secondary ADHD had more lesions in the right putamen than children who did not develop secondary ADHD; this impression was confirmed statistically. After Bonferroni correction, we could not demonstrate significant differences between secondary ADHD status and lesion burdens for the right caudate nucleus or the right globus pallidus. CONCLUSION: Closed-head injury-induced lesions in the right putamen in children are associated with subsequent development of secondary ADHD. Depiction software is useful in guiding statistical analysis of image data.

NASA Discipline Neuroscience↗

Large-eddy simulation of helical- and straight-bladed vertical-axis wind turbines in boundary layer turbulence

Here, turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines (VAWTs) in boundary layer turbulence are numerically studied using the large-eddy simulation (LES) method combined with the actuator line model. Based on the LES data, systematic statistical analyses are performed to explore the effects of blade geometry on the characteristics of the turbine wake. The time-averaged velocity fields show that the helical-bladed VAWT generates a mean vertical velocity along the center of the turbine wake, which causes a vertical inclination of the turbine wake and alters the vertical gradient of the mean streamwise velocity. Consequently, the intensities of the turbulent fluctuations and Reynolds shear stresses are also affected by the helical-shaped blades when compared with those in the straight-bladed VAWT case. The LES results also show that reversing the twist direction of the helical-bladed VAWT causes the spatial patterns of the turbulent wake flow statistics to be reversed in the vertical direction. Moreover, the mass and kinetic energy transports in the turbine wakes are directly visualized using the transport tube method, and the comparison between the helical- and straight-bladed VAWT cases show significant differences in the downstream evolution of the transport tubes.

17 WIND ENERGY↗

Fast Grain Mapping with Sub-Nanometer Resolution Using 4D-STEM with Grain Classification by Principal Component Analysis and Non-Negative Matrix Factorization

High-throughput grain mapping with sub-nanometer spatial resolution is demonstrated using scanning nanobeam electron diffraction (also known as 4D scanning transmission electron microscopy, or 4D-STEM) combined with high-speed direct-electron detection. An electron probe size down to 0.5 nm in diameter is used and the sample investigated is a gold–palladium nanoparticle catalyst. Computational analysis of the 4D-STEM data sets is performed using a disk registration algorithm to identify the diffraction peaks followed by feature learning to map the individual grains. Two unsupervised feature learning techniques are compared: principal component analysis (PCA) and non-negative matrix factorization (NNMF). The characteristics of the PCA versus NNMF output are compared and the potential of the 4D-STEM approach for statistical analysis of grain orientations at high spatial resolution is discussed.

47 OTHER INSTRUMENTATION↗

Dual spectra and mixed energy cascade of turbulence in the wavelet representation

The wavelet-transformed Navier-Stokes equations are used to define quantities such as the transfer of kinetic energy and the flux of kinetic energy through scale r at position x. Direct numerical simulations of turbulent shear flow reveal that although their mean spatial values agree with their traditional counterparts in Fourier space, their spatial variability at every scale is very large, exhibiting non-Gaussian statistics. The local flux of energy involving scales smaller than some r also exhibits large spatial intermittency, and it is negative quite often, indicative of local inverse cascades.

Meneveau, Charles↗

Asymptotic modal analysis and statistical energy analysis

Asymptotic Modal Analysis (AMA) is a method which is used to model linear dynamical systems with many participating modes. The AMA method was originally developed to show the relationship between statistical energy analysis (SEA) and classical modal analysis (CMA). In the limit of a large number of modes of a vibrating system, the classical modal analysis result can be shown to be equivalent to the statistical energy analysis result. As the CMA result evolves into the SEA result, a number of systematic assumptions are made. Most of these assumptions are based upon the supposition that the number of modes approaches infinity. It is for this reason that the term 'asymptotic' is used. AMA is the asymptotic result of taking the limit of CMA as the number of modes approaches infinity. AMA refers to any of the intermediate results between CMA and SEA, as well as the SEA result which is derived from CMA. The main advantage of the AMA method is that individual modal characteristics are not required in the model or computations. By contrast, CMA requires that each modal parameter be evaluated at each frequency. In the latter, contributions from each mode are computed and the final answer is obtained by summing over all the modes in the particular band of interest. AMA evaluates modal parameters only at their center frequency and does not sum the individual contributions from each mode in order to obtain a final result. The method is similar to SEA in this respect. However, SEA is only capable of obtaining spatial averages or means, as it is a statistical method. Since AMA is systematically derived from CMA, it can obtain local spatial information as well.

Dowell, Earl H.↗

Additional studies of forest classification accuracy as influenced by multispectral scanner spatial resolution

First, an analysis of forest feature signatures was used to help explain the large variation in classification accuracy that can occur among individual forest features for any one case of spatial resolution and the inconsistent changes in classification accuracy that were demonstrated among features as spatial resolution was degraded. Second, the classification rejection threshold was varied in an effort to reduce the large proportion of unclassified resolution elements that previously appeared in the processing of coarse resolution data when a constant rejection threshold was used for all cases of spatial resolution. For the signature analysis, two-channel ellipse plots showing the feature signature distributions for several cases of spatial resolution indicated that the capability of signatures to correctly identify their respective features is dependent on the amount of statistical overlap among signatures. Reductions in signature variance that occur in data of degraded spatial resolution may not necessarily decrease the amount of statistical overlap among signatures having large variance and small mean separations. Features classified by such signatures may thus continue to have similar amounts of misclassified elements in coarser resolution data, and thus, not necessarily improve in classification accuracy.

Sadowski, F. E.↗

Vibration Response Models of a Stiffened Aluminum Plate Excited by a Shaker

Numerical models of structural-acoustic interactions are of interest to aircraft designers and the space program. This paper describes a comparison between two energy finite element codes, a statistical energy analysis code, a structural finite element code, and the experimentally measured response of a stiffened aluminum plate excited by a shaker. Different methods for modeling the stiffeners and the power input from the shaker are discussed. The results show that the energy codes (energy finite element and statistical energy analysis) accurately predicted the measured mean square velocity of the plate. In addition, predictions from an energy finite element code had the best spatial correlation with measured velocities. However, predictions from a considerably simpler, single subsystem, statistical energy analysis model also correlated well with the spatial velocity distribution. The results highlight a need for further work to understand the relationship between modeling assumptions and the prediction results.

Cabell, Randolph H.↗

Characterizing multiscale variability of zero intermittency in spatial rainfall

In this paper the authors study how zero intermittency in spatial rainfall, as described by the fraction of area covered by rainfall, changes with spatial scale of rainfall measurement or representation. A statistical measure of intermittency that describes the size distribution of 'voids' (nonrainy areas imbedded inside rainy areas) as a function of scale is also introduced. Morphological algorithms are proposed for reconstructing rainfall intermittency at fine scales given the intermittency at coarser scales. These algorithms are envisioned to be useful in hydroclimatological studies where the rainfall spatial variability at the subgrid scale needs to be reconstructed from the results of synoptic- or mesoscale meteorological numerical models. The developed methodologies are demsonstrated and tested using data from a severe springtime midlatitude squall line and a mild midlatitude winter storm monitored by a meteorological radar in Norman, Oklahoma.

Kumar, Praveen↗

Severe Dirac Mass Gap Suppression in Sb 2 Te 3 -Based Quantum Anomalous Hall Materials

The quantum anomalous Hall (QAH) effect appears in ferromagnetic topological insulators (FMTIs) when a Dirac mass gap opens in the spectrum of the topological surface states (SSs). Unaccountably, although the mean mass gap can exceed 28 meV (or ~320 K), the QAH effect is frequently only detectable at temperatures below 1 K. Using atomic-resolution Landau level spectroscopic imaging, we compare the electronic structure of the archetypal FMTI Cr 0.08 (Bi 0.1 Sb 0.9 ) 1.92 Te 3 to that of its nonmagnetic parent (Bi 0.1 Sb 0.9 ) 2 Te 3 , to explore the cause. In (Bi 0.1 Sb 0.9 ) 2 Te 3 , we find spatially random variations of the Dirac energy. Statistically equivalent Dirac energy variations are detected in Cr 0.08 (Bi 0.1 Sb 0.9 ) 1.92 Te 3 with concurrent but uncorrelated Dirac mass gap disorder. Additionally, these two classes of SS electronic disorder conspire to drastically suppress the minimum mass gap to below 100 μeV for nanoscale regions separated by <1 μm. This fundamentally limits the fully quantized anomalous Hall effect in Sb 2 Te 3 -based FMTI materials to very low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computational budget optimization for Bayesian parameter estimation in heavy-ion collisions

Abstract Bayesian parameter estimation provides a systematic approach to compare heavy-ion collision models with measurements, leading to constraints on the properties of nuclear matter with proper accounting of experimental and theoretical uncertainties. Aside from statistical and systematic model uncertainties, interpolation uncertainties can also play a role in Bayesian inference, if the model’s predictions can only be calculated at a limited set of model parameters. This uncertainty originates from using an emulator to interpolate the model’s prediction across a continuous space of parameters. In this work, we study the trade-offs between the emulator (interpolation) and statistical uncertainties. We perform the analysis using spatial eccentricities from the T R ENTo model of initial conditions for nuclear collisions. Given a fixed computational budget, we study the optimal compromise between the number of parameter samples and the number of collisions simulated per parameter sample. For the observables and parameters used in the present study, we find that the best constraints are achieved when the number of parameter samples is slightly smaller than the number of collisions simulated per parameter sample.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Data & Code from Phoenix CPPP Phase 2 Analysis

This data and code package supports the analysis presented in “Beyond Surface Cooling: Comprehensive Field Assessment of Reflective Pavement Thermal Performance in Phoenix, Arizona” and provides fully reproducible workflows for evaluating the thermal performance of cool pavement treatments in a hot urban environment. The dataset integrates multi-modal field measurements collected across residential and nonresidential settings, including mobile air temperature traverses, stationary air temperature monitoring, residential mean radiant temperature (MRT) measurements, subsurface temperature profiles, and controlled testbed observations. The data package contains raw and processed datasets in comma-separated value (CSV) format, accompanying metadata files describing site characteristics and measurement protocols, and R scripts (.R files) used for data cleaning, time synchronization, spatial and temporal matching, quality control filtering, statistical comparison, and figure generation. All analyses were conducted using R (version ≥ 4.2.0) with commonly available packages (e.g., tidyverse, lubridate, data.table, ggplot2). No proprietary software is required to reproduce results. Field campaigns were designed to quantify the effects of high-reflectance pavement coatings on surface temperature, near-surface air temperature, subsurface heat propagation, and radiative heat exposure. Temporal alignment procedures include standardized timestamp conversion and nearest-neighbor matching of high-frequency sensor measurements to stop-based metadata within defined tolerance windows to ensure comparability across instruments. The workflows generate summary statistics, treatment–control contrasts, depth-dependent thermal gradients, and time-series visualizations used in the associated publication. By integrating mobile, stationary, radiative, and subsurface measurements within a unified and transparent processing framework, this package enables comprehensive evaluation of cool pavement performance across multiple thermal exposure pathways and supports reuse in future urban heat mitigation and climate resilience studies.

AIR TEMPERATURE↗

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗