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 37 records · Page 2

Morphostatistical characterization of the spatial galaxy distribution through Gibbs point processes

ABSTRACT This paper proposes a morphostatistical characterization of the galaxy distribution through spatial statistical modelling based on inhomogeneous Gibbs point processes. The galaxy distribution is supposed to exhibit two components. The first one is related to the major geometrical features exhibited by the observed galaxy field, here, its corresponding filamentary pattern. The second one is related to the interactions exhibited by the galaxies. Gibbs point processes are statistical models able to integrate these two aspects in a probability density, controlled by some parameters. Several such models are fitted to real observational data via the ABC shadow algorithm. This algorithm provides simultaneous parameter estimation and posterior-based inference, hence allowing the derivation of the statistical significance of the obtained results.

79 ASTRONOMY AND ASTROPHYSICS↗

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗

Data Summarization and Inference at Scale

This is the final report for the DOE ASCR grant SC-0022260, Data Summarization and Inference at Scale, PI: Alex Pothen, Purdue University. The goal of the project was to solve data-intensive and compute-intensive problems in the physical sciences, engineering, information science, data science, etc. by designing and implementing new algorithms that could work with a subset of the data. The four subgoals were: (a) The solution of problems where the data is too large to be stored in the memory of a computer. In this streaming model of computation, the data arrives as a stream of elements to the computer, each element is processed as it arrives, and a decision is made to discard the data or to store it; only a small subset of the data proportional to the size of the output solution is stored, and when all the data has been streamed, a solution to the problem is computed from the stored subset. (b) The use of machine learning methods to compute solutions to data-intensive problems. The use of GPUs is critical to obtain high performance on machine learning tasks, but their memory sizes are smaller relative to that of CPUs. For large-scale problems, the data is sampled many times, and small samples are used with repetition, for robustness, to compute solutions to inference tasks. This sampling reduces the memory required to solve the problem, but attention is needed to avoid slow convergence to the solutions, and reduced accuracy of inference. We propose submodular optimization, Large Language Models, and physics-informed neural networks to enable GPU computations here. (c) Modeling and visualization of high-dimensional data using interpretable features. Clinical proteomic data sets from immunology for the detection of cancer and other diseases are temporal and high-dimensional, and algorithms for visualizing these data sets using clinically interpretable features are lacking. We propose methods that compute distances based on the optimal transportation problem and graph edit distances to address this problem. We also propose the use of optimal transport-based distances, spatial statistics, and network structure to classify image data sets, We apply these algorithms to electron micrographs of the peripheral nervous system in the digestive tract. (d) The design of data-intensive algorithms on emerging architectures, specifically, noisy, intermediate-scale quantum (NISQ) devices. Quantum computers offer the possibility of exploring large solution spaces due to the principle of superposition, but current quantum computers are limited by few qubits, short coherence times due to noise, poor interconections among the qubits, etc. We propose the use of the divide and conquer paradigm to solve large-scale problems, wherein collections of small subproblems are solved on the quantum devices, and the solutions to the subproblems are integrated into a solution for the original problem on a classical computer.

97 MATHEMATICS AND COMPUTING↗

Data from: Understanding the biogeochemical and spatial drivers of methane and carbon dioxide fluxes in a large temperate reservoir

This dataset contains spatially resolved measurements of CO₂ and CH₄ fluxes and associated environmental variables collected across 200 sites in Douglas Reservoir (Tennessee, USA) between July 29-August 2, 2024. Measurements include diffusive fluxes of CO₂ and CH₄, CH₄ ebullition, and biogeochemical and spatial variables such as dissolved oxygen, temperature, conductivity, chlorophyll-a, pH, water depth, and distance from the dam. Sampling was conducted using a spatially balanced design to capture longitudinal and depth-related gradients throughout the reservoir. The dataset is structured to support analyses of spatial variability, flux pathway comparisons, and modeling approaches (e.g., spatial statistics and machine learning) aimed at understanding controls on reservoir CO₂ and CH₄ fluxes and improving upscaling to whole-reservoir and regional estimates.

Neeper, Jamie [ORNL] (ORCID:0009000842342101)↗

Assessing Current and Future Infrastructure Hazards

Project Objective Execute intelligent analytics via an advanced analytical framework, to assess the current state of offshore infrastructure, evaluate infrastructure life, and identify technologies to reduce infrastructure hazards, costs, and extend infrastructure life. Approach &amp; Results thus Far• Build comprehensive dataset• Perform data-driven analytics to evaluate infrastructure integrity<p> 1. Remaining lifespan</p><p> 2. Likelihood of future risk</p><p>• Apply data-driven advanced spatial, statistical, and Machine Learning (ML) models to quantify existing infrastructure integrity</p><p>• Release data and models through a smart, online platform hosted by Energy Data eXchange (EDX)</p>

Romeo, Lucy F.↗

Statistical analysis of spatial distribution of pores in metal additive manufacturing

We report mechanical properties of additively manufactured parts are sensitive to the presence of pores form during the manufacturing process. The impact of pores on the mechanical performance has been investigated extensively with respect to different parameters, such as pore volume fraction, shape, and size. However, statistical investigations focusing on the relationships between the spatial distribution of pores and process parameters; and consequently, the performance of the manufactured parts are scattered and limited. Also, these sparse investigations usually suffer from an ambiguous definition of terminologies. For instance, the required criteria to consider a point pattern as complete spatial randomness (CSR) are generally not clarified. Moreover, no numerical formalism is yet developed to show how much the observed spatial results are statistically significant. To address these shortcomings, the statistical definition of CSR in a point pattern and the procedure to quantitatively determine the deviation of a pattern from CSR were discussed. The explained statistical approach was used to investigate the effect of scanning speed parameter on the spatial distribution of gas pores in laser powder bed fusion manufactured stainless steel parts. Furthermore, and to highlight the impact of the spatial distribution of pores on mechanical properties, fatigue performances of parts with clustered and randomly distributed pores were simulated by finite element analysis. It is shown that by reducing the scanning speed, the spatial distribution of gas pores deviated more from CSR, and correspondingly fatigue performance deteriorated.

36 MATERIALS SCIENCE↗

A sensitivity analysis to predict the neutronics behavior of samples irradiated in the VTR rabbit system

We report a low-order neutronics model is developed to carry out hundreds of simulations efficiently and investigate the neutronics behavior of samples being irradiated in a test reactor setting under different geometrical constraints. The low-order model allowed for simulations that yield the expected neutronics behavior of any irradiated sample in any environment and allows for the calculation of highly accurate spatially averaged statistics and idealized spatial distributions in the neutron flux. Several benchmarks are performed to evaluate the performance and limitations of the low-order model revealing many important findings. The low-order model predicted the LHGR in the EBR-II driver fuel to within 2.34% by only simulating the fuel rod by itself, which served as a validation for the model. Sensitivity studies investigated 3% enriched UO 2 and U-10Zr being irradiated in the Versatile Test Reactor rabbit system. The analyses investigated a range of combinations of 15 radii and 5 heights for each sample in the rabbit system. Similar data sets are also provided for irradiations in the Advanced Test Reactor’s B-10 irradiation position, which is a thermal neutron spectrum environment. Generalized fits and fit coefficients are obtained for sample heating, reaction rate densities, and local multiplication rate characteristics, allowing the predictions of the neutronics behavior of the samples based on their geometrical constraints. The analyses and fits laid the groundwork for developing a user-end Multiphysics analysis framework to assist and accelerate irradiation experiment design and optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Spatially-resolved soot evolution and statistics in high-pressure diesel spray flames using two-color pyrometry

Previous studies of diesel spray flames have focused on details of the sooting behavior mainly during the quasi-steady period, but few have considered the spatially-resolved transient evolution in combination with injection-to-injection variations. In this study, a 500-injection data set is utilized to investigate the temporal evolution of a spray flame during auto-ignition, the premixed burn phase, and the quasi-steady period. Spatially-resolved OH* chemiluminescence data provide ignition delay times and reaction zone locations. Two-color pyrom-etry with a vastly improved optical system is used to quantify spatially-resolved soot evolution and its statistical variations. The ambient thermodynamic con-ditions are slightly below those in modern diesel engines, resulting in longer lift-off lengths and lower overall soot production. Spatially, soot formation in the lift-off region is small, but appreciable soot forms in the jet core and jet head regions, while it oxidizes gradually on the jet periphery. Total soot mass profiles indicate that regions with larger local soot mass take longer to form, but are the first to oxidize. Probability distributions of soot mass in localized regions indicate that a few injections with high soot mass bias the average soot mass data towards higher values. Altogether, results show good agreement with previous studies employing different diagnostic techniques while providing statistical details of transient and localized soot behavior in high pressure diesel spray flames in support of the related modeling efforts.

09 BIOMASS FUELS↗

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗

Classification of Intensity Distributions of Transmission Eigenchannels of Disordered Nanophotonic Structures Using Machine Learning

Light-matter interaction optimization in complex nanophotonic structures is a critical step towards the tailored performance of photonic devices. The increasing complexity of such systems requires new optimization strategies beyond intuitive methods. For example, in disordered photonic structures, the spatial distribution of energy densities has large random fluctuations due to the interference of multiply scattered electromagnetic waves, even though the statistically averaged spatial profiles of the transmission eigenchannels are universal. Classification of these eigenchannels for a single configuration based on visualization of intensity distributions is difficult. However, successful classification could provide vital information about disordered nanophotonic structures. Emerging methods in machine learning have enabled new investigations into optimized photonic structures. In this work, we combine intensity distributions of the transmission eigenchannels and the transmitted speckle-like intensity patterns to classify the eigenchannels of a single configuration of disordered photonic structures using machine learning techniques. Specifically, we leverage supervised learning methods, such as decision trees and fully connected neural networks, to achieve classification of these transmission eigenchannels based on their intensity distributions with an accuracy greater than 99%, even with a dataset including photonic devices of various disorder strengths. Simultaneous classification of the transmission eigenchannels and the relative disorder strength of the nanophotonic structure is also possible. Our results open new directions for machine learning assisted speckle-based metrology and demonstrate a novel approach to classifying nanophotonic structures based on their electromagnetic field distributions. These insights can be of paramount importance for optimizing light-matter interactions at the nanoscale.

79 ASTRONOMY AND ASTROPHYSICS↗

Closures and Simulation for Thermal Radiation Transport in Stochastic Media with Nonlinear Temperature Dependence

Because of the practical challenge of rendering very complex realistic spatial structures for numerical work, it is common practice to resort to characterizing such media as stochastic mixtures of materials, ideally parametrized with low order statistics such as the mean, variance, and correlation functions of the now random material properties. This enables realizations of the medium to be repeatedly generated and radiation transport computations to, in principle, be performed for a large ensemble of these realizations to obtain a statistically well-characterized radiation field. Statistical post-processing yields desired quantities such as conditional and unconditional mean radiation flux and probability distributions of transmitted radiation. However, such computations prove expensive for all but the simplest stochastic geometries and are most suited for benchmarking approximate models. The most common approximations lead to homogenized media so that transport computations are required only on a single medium realization but by construct provide only limited statistical information on the radiation field. Almost all approximate approaches to this problem attempt to develop equations for low order moments of the radiation intensity (mean, second moment, correlation function) but inevitably encounter a closure problem: the equation for any statistical moment will contain terms depending on unknown higher-order moments. Thus, the challenge shifts to one of developing closure relations that relate the unknown moments to the lower order moments. Under very special conditions, an exact closure can be derived but in general closures are heuristically stated constitutive relations. Also, closure approaches depend on whether the mixing statistics are spatially and/or temporally continuous as in fluctuating turbulent fields, or discontinuous as in randomly mixed solid chunks of material. Thus, unconditional averaging is generally applied in the former case but conditional averaging is more appropriate when the mixing is discontinuous. In this work, the emphasis is on binary statistical mixtures of immiscible fluids as as such the quantities of interest are averages (flux, temperature) conditioned on the material type.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Geospatial analysis of preterm and small-for-gestational age births in Washington D.C.

Background: This study is based on the recognition that adverse pregnancy outcomes significantly affect maternal and infant health, leading to increased morbidity and mortality. These outcomes are shaped by a complex interplay of individual-level factors—like maternal age and education—and community-level influences, including socio-economic status and access to healthcare. Understanding these determinants is crucial for developing effective public health strategies, especially for marginalized populations, by identifying high-risk areas and informing targeted interventions that address both individual and structural barriers. Methods: We utilized geospatial analysis to explore the association between individual- and community-level factors and adverse pregnancy outcomes, specifically preterm birth (PTB) and small-for-gestational-age (SGA) birthweight in Washington, D.C. We used Empirical Bayes smoothing methods to calculate rates of adverse birth outcomes from 2010 to 2018 at the U.S. Census tract–level. Spatial scan statistics were used to investigate if adverse birth outcomes clustered in specific areas. ANOVA tests were conducted for individual- and community-level factors within identified clusters. Results: Spatial analysis identified significant high-risk clusters for PTB and SGA infants primarily in southeastern Washington, D.C., particularly in Wards 7 and 8. Individuals residing within these clusters experienced a 47% increased risk of PTB (RR = 1.467) and a 56% increased risk of SGA (RR = 1.560) compared to those outside clusters. Space–time analysis revealed temporal variation, with PTB clusters persisting from 2011 to 2014 and SGA clusters extending through 2017. Compared to low-risk clusters, high-risk clusters had younger birthing individuals (mean age ~26.5 vs. ~33 years), lower maternal college degree attainment (~20% vs. ~80%), higher rates of late or no prenatal care (~16% vs. 11%), and increased prevalence of smoking and hypertension (all P < 0.001). Community-level indicators showed lower median household incomes ($\$40,000$ vs. ~$\$105,000$), greater poverty (~16% vs. ~7% below $\$10,000$/year), higher public assistance use (~32% vs. ~5%), and reduced healthcare access (greater distances to emergency and specialty care) in high-risk areas (all P < 0.001). Neighborhood deprivation indices were significantly elevated, commutes were longer, and population density was lower in these clusters. These findings highlight that adverse birth outcomes cluster in neighborhoods with pronounced socioeconomic and health disparities. Conclusion: High-risk birth clusters highlight intertwined factors: individual, socio-economic, and geographic. Addressing these requires comprehensive interventions focusing on social and structural determinants of health.

Birth outcomes↗

Boron-Based Neutron Scintillator Screens for Neutron Imaging

In digital neutron imaging, the neutron scintillator screen is a limiting factor of spatial resolution and neutron capture efficiency and must be improved to enhance the capabilities of digital neutron imaging systems. Commonly used neutron scintillators are based on 6 LiF and gadolinium oxysulfide neutron converters. This work explores boron-based neutron scintillators because 10 B has a neutron absorption cross-section four times greater than 6Li, less energetic daughter products than Gd and 6 Li, and lower γ-ray sensitivity than Gd. These factors all suggest that, although borated neutron scintillators may not produce as much light as 6 Li-based screens, they may offer improved neutron statistics and spatial resolution. This work conducts a parametric study to determine the effects of various boron neutron converters, scintillator and converter particle sizes, converter-to-scintillator mix ratio, substrate materials, and sensor construction on image quality. The best performing boron-based scintillator screens demonstrated an improvement in neutron detection efficiency when compared with a common 6 LiF/ZnS scintillator, with a 125% increase in thermal neutron detection efficiency and 67% increase in epithermal neutron detection efficiency. The spatial resolution of high-resolution borated scintillators was measured, and the neutron tomography of a test object was successfully performed using some of the boron-based screens that exhibited the highest spatial resolution. For some applications, boron-based scintillators can be utilized to increase the performance of a digital neutron imaging system by reducing acquisition times and improving neutron statistics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗

Shrubs Strongly Influence Snow Properties in Two Subarctic Watersheds

Understanding changes in snow distribution in permafrost ecosystems is fundamental to predicting their response to future climate change. The expansion of tall shrubs into tundra ecosystems can trap snow and insulate permafrost ecosystems during the winter, but the overall insulation effect is dependent upon many ecosystem properties. To study shrub–snow–ground interactions, small temperature sensors were deployed at two research sites on the Seward Peninsula of Alaska, USA, during the 2019–2020 winter. Snow temperatures were used to extrapolate multiple metrics, including freezing n-factors, the snow insulation effect, snow cover duration, and the length of the snowmelt period. Statistical and spatial analysis showed that shrub patches were a dominant control on all snow metrics. Within shrub patches, average ground temperatures were 2.1°C warmer, snow persisted 50 days longer, snow insulation was double, and a longer, later spring snowmelt period occurred compared to nonshrubby areas. Site-level differences contributed relatively little to variation in snow metrics, indicating that shrub presence is an overarching driver of snow–ground interactions at the locations examined. Shrub expansion, which is anticipated under climate change, will strongly impact future permafrost distribution and Arctic energy, water, and carbon cycles through snow–shrub–ground feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Critical impact of experimentally-driven strut level anisotropic material models in advanced stress analysis of additively manufactured lattice structures

The rapid acceleration in materials discovery may overshadow the importance of thoroughly understanding the mechanical performance of newly developed materials in demanding environments. The recent interest in combining parametric studies with machine learning techniques to explore how changes in specific processing parameters or model inputs affect the overall behavior of a material system can only be truly beneficial if the governing constitutive relations describing material behavior are accurately established. In this study, we demonstrate the critical impact of accurately representing strut-level anisotropic material behavior in advanced stress analysis of additively manufactured lattice structures (AMLS). We introduce a systematic experimental and modeling approach for developing strut-level anisotropic elastoplastic material models that account for the influence of microstructural features such as porosity, texture, and surface roughness on the development of local anisotropic mechanical properties, which vary with strut orientation relative to the build direction (BD). As a result the presented material model captures and relates the statistics of spatially varying struts’ microstructural features to the local stress distribution. Our findings suggest that incorporating strut-level anisotropic material behavior into unit cell analysis significantly influences the load distribution and evolution of local stresses within the structure. Therefore, accounting for this anisotropy is critical for developing an understanding of unit cell behavior and performance, including subsequent topology/component design optimization based on this analysis.

Sahoo, Subhadip [University of Arizona]↗

Three dimensional cluster analysis for atom probe tomography using Ripley’s K-function and machine learning

The size and structure of spatial molecular and atomic clustering can significantly impact material properties and is therefore important to accurately quantify. Ripley’s K-function (K(r)), a measure of spatial correlation, can be used to perform such quantification when the material system of interest can be represented as a marked point pattern. This work demonstrates how machine learning models based on K (r)-derived metrics can accurately estimate cluster size and intra-cluster density in simulated three dimensional (3D) point patterns containing spherical clusters of varying size; over 90% of model estimates for cluster size and intra-cluster density fall within 11% and 18% error of the true values, respectively. These K (r)-based size and density estimates are then applied to an experimental APT reconstruction to characterize MgZn clusters in a 7000 series aluminum alloy. Here we find that the estimates are more accurate, consistent, and robust to user interaction than estimates from the popular maximum separation algorithm. Using K (r) and machine learning to measure clustering is an accurate and repeatable way to quantify this important material attribute.

36 MATERIALS SCIENCE↗

Multidimensional Distributional Neural Network Output Demonstrated in Super‐Resolution of Surface Wind Speed

Accurate quantification of uncertainty in neural network predictions remains a central challenge for scientific applications involving high-dimensional, correlated data. While existing methods capture either aleatoric or epistemic uncertainty, few offer closed-form, multidimensional distributions that preserve spatial correlation while remaining computationally tractable. In this work, we present a framework for training neural networks with a multidimensional Gaussian loss, generating a closed-form predictive distribution over outputs informed by non-identically distributed training data. Our approach captures aleatoric uncertainty by iteratively estimating the means and covariance matrices, and is demonstrated on a super-resolution example out-of-training-sample. We leverage a Fourier representation of the covariance matrix to stabilize network training and preserve spatial correlation. We introduce a novel regularization strategy—referred to as information sharing—that interpolates between image-specific and global covariance estimates, enabling convergence of the super-resolution downscaling network trained on image-specific distributional loss functions. This framework allows for efficient sampling, explicit correlation modeling, and extensions to more complex distribution families all without disrupting prediction performance. We demonstrate the method on a surface wind speed downscaling task and discuss its broader applicability to uncertainty-aware prediction in scientific models.

17 WIND ENERGY↗