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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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↗
DOE OSTI2024

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↗
DOE OSTI2022

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↗
DOE OSTI2022

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↗
DOE OSTI2025

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↗
DOE OSTI2022

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↗
DOE OSTI2024

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↗
DOE OSTI2025

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]↗
DOE OSTI2026

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↗
DOE OSTI2022

Hyper-fidelity depletion coupled with discrete pebble motion in pebble bed reactors

Pebble bed reactors have raised new interest during the past decade due to their attractive characteristics. Therefore, accurate simulations must be performed to better understand these systems and ensure optimal and safe designs. Most current methods use lower fidelity approaches with representative unit-cells or macro-zones with uniform fluxes, which have accuracy and flexibility limitations. A novel hyper-fidelity method for pebble bed reactors depletion is presented and internally couples Serpent 2 and a pseudo-motion routine. Pseudo-motion is applied handling vertical shifts of compositions in a static pebble bed, random reinsertion of used pebbles, insertion of fresh pebbles and used pebbles discarding. Associated with individual depletion to correctly determine the flux spectrum and composition in each pebble, this hyper-fidelity approach paves the way towards more accurate depletion calculation in pebble bed reactors. Using this method, a demonstration is completed on a small-scale reactor. In this application, the core reaches equilibrium, and the following data is extracted: core-wise parameters evolution, pebble-wise spatial and statistical distribution. Discarded pebbles are analyzed, and relevant information is shown. This work proves the feasibility of hyper fidelity depletion with Serpent 2, and the range of use for this method: reactor design and analysis for equilibrium and slow transients, lower fidelity methods validation and feeding fuel performance, thermal-hydraulics, or waste management models. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

NEWTS Well Summary by Hydrologic Regions and Subbasins in the U.S.

Oil and gas well production data (2000-2022), including water production and injection, as well as summary information (e.g., well count by status, total vertical depth statistics, etc.) spatially summarized by Watershed Boundary Dataset's (WBD) region (Hydrologic Unit Code (HUC) 2) and subbasin (HUC8). Energy-related produced waters data were acquired and summarized in support of the development of the National Energy Technology Laboratory's NEWTS (National Energy Water Treatment and Speciation Database). Due to the proprietary nature of the wellbore data, this derived product has been spatially compiled by key areas to support research and stakeholder needs.

Energy Infrastructure↗
DOE OSTI2026

Anyonic Membranes and Pontryagin Statistics

Anyons, unique to two spatial dimensions, underlie extraordinary phenomena such as the fractional quantum Hall effect, but their generalization to higher dimensions has remained elusive. The topology of Eilenberg-MacLane spaces constrains the loop statistics to be only bosonic or fermionic in any dimension. In this work, we introduce the novel anyonic statistics for membrane excitations in four dimensions. Analogous to the $\mathbb{Z}_N$-particle exhibiting $\mathbb{Z}_{N\times \gcd(2,N)}$ anyonic statistics in two dimensions, we show that the $\mathbb{Z}_N$-membrane possesses $\mathbb{Z}_{N\times \gcd(3,N)}$ anyonic statistics in four dimensions. Given unitary volume operators that create membrane excitations on the boundary, we propose an explicit 56-step unitary sequence that detects the membrane statistics. We further analyze the boundary theory of $(5{+}1)$D 1-form $\mathbb{Z}_N$ symmetry-protected topological phases and demonstrate that their domain walls realize all possible anyonic membrane statistics. We then show that the $\mathbb{Z}_3$ subgroup persists in all higher dimensions. In addition to the standard fermionic $\mathbb{Z}_2$ membrane statistics arising from Stiefel-Whitney classes, membranes also exhibit $\mathbb{Z}_3$ statistics associated with Pontryagin classes. We explicitly verify that the 56-step process detects the nontrivial $\mathbb{Z}_3$ statistics in 5, 6, and 7 spatial dimensions. Furthermore, in 7 and higher dimensions, the statistics of membrane excitations stabilize to $\mathbb{Z}_{2} \times \mathbb{Z}_{3}$, with the $\mathbb{Z}_3$ sector consistently captured by this process.

Abstract algebra↗
DOE OSTI2025

Map-level baryonification: unified treatment of weak lensing two-point and higher-order statistics

Precision cosmology benefits from extracting maximal information from cosmic structures, motivating the use of higher-order statistics (HOS) at small spatial scales. However, predicting how baryonic processes modify matter statistics at these scales has been challenging. The baryonic correction model (BCM) addresses this by modifying dark-matter-only simulations to mimic baryonic effects, providing a flexible, simulation-based framework for predicting both two-point and HOS. We show that a 3-parameter version of the BCM can jointly fit weak lensing maps' two-point statistics, wavelet phase harmonics coefficients, scattering coefficients, and the third and fourth moments to within 2% accuracy across all scales ℓ < 2000 and tomographic bins for a DES-Y3-like redshift distribution ( z ≲ 2), using the FLAMINGO simulations. These results demonstrate the viability of BCM-assisted, simulation-based weak lensing inference of two-point and HOS, paving the way for robust cosmological constraints that fully exploit non-Gaussian information on small spatial scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Correlation-aware binning for small-angle neutron scattering via Gaussian-process inference

Binning in small-angle neutron scattering (SANS) is typically performed empirically, with fixed parameters chosen for convenience rather than statistical optimality. Such practices often fail to balance statistical precision and spatial resolution, leading to inconsistencies across instruments and datasets. Here we establish a correlation-aware framework that determines the optimal bin width from first principles by extending the classical Freedman–Diaconis (FD) rule to account for inter-bin correlations with a Gaussian process. In this formulation, the scattering intensity is treated as a smooth stochastic field whose statistical coherence is described by a covariance matrix. Analytical expressions of errors derived from this model yield closed-form criteria that separate the total deviation into contributions from counting noise, aliasing distortion and curvature-dependent correlation effects. Expressed in reduced variables, the resulting dimensionless error surface reveals a continuous transition from the uncorrelated FD regime to the correlation-dominated limit, providing a unified description of noise suppression and resolution control. Because the formulation depends only on the profile characteristics of scattering intensity I(Q), specifically its average intensity and first- and second-order derivatives, it applies generally to any SANS measurement regardless of sample, instrument or geometry. Experimental validation using small- and ultra-small-angle neutron scattering data confirms the predicted scaling behavior, demonstrating that correlation-aware inference systematically reduces mean-squared error and enables information-efficient reproducible data reduction across materials and instruments.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗
DOE OSTI2022

The complex structural and chemical nature of monolithic U-10Mo fuel and Zr barrier layer

Nanoscale microstructural characterization by advanced transmission electron microscopy techniques on a U-10Mo/Zr barrier layer monolithic fuel plate was performed to evaluate the microstructural evaluation after high burn-up. Gas bubble superlattice evolution, grain restructuring, and evolution of the Zr interaction layer is investigated through detailed electron microscopy characterization. The use of automated crystallographic orientation mapping to irradiated U-10Mo fuel highlights that the restructured ultra-fine grains are separated by high angle grain boundaries at a burn up of 4.42 × 10 21 fissions/cm 3 . Additionally, advanced chemical analysis and multi-variable statistical analysis shows spatial clustering of solid fission product precipitates. Finally, characterization of a newly observed porous nanocrystalline Zr region in the barrier layer is studied. Finally, this work provides insights into the grain subdivision and restructuring process while using advanced microscopy techniques to analyze fission products in neutron-irradiated U-10Mo fuel.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data-Driven Insights into the Structural Essence of Plasticity in High-Entropy Alloys

The heterogeneous mechanical response of a crystalline alloy with multiple principal elements was investigated using molecular dynamics simulations. The local configuration of the alloy in its quiescent state was characterized by the variables derived from the gyration tensor and the atomic electronegativity. A multivariate analysis identified the geometric and chemical factors that influenced the atomic packing variations. Further, upon straining, the non-affine displacement exhibited spatial heterogeneity. A statistical correlation was established between the local yield events and the specific features of the local configuration. Our findings, validated by the performance metrics analysis, provided a structural criterion for the instability mechanisms in high-entropy alloys (HEAs) and enhanced the understanding of their plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Light scattering by V 4 O 7 film across the metal–insulator transition

The experimental study of the angle-resolved hemispherical light scattering by V 4 O 7 film within a broad temperature range across metal–insulator transition reveals complex structural reorganization of the film deposited on the c-cut sapphire crystal. The bidirectional scattering distribution function and the surface autocorrelation function were obtained from scattering data to visualize statistics of the spatially resolved contributions of optical inhomogeneities in normal and lateral directions to the surface. The measurements reveal an anisotropic surface roughness distribution due to the twinned domain structure, with significant anisotropy changes across the phase transition. The V 4 O 7 film deposited on sapphire leads to a polydomain structure, minimizing elastic strain energy with distinct multiscale distributions of surface domains. Near T c , the material shows the lowest roughness but the highest lateral disorder of the surface.

36 MATERIALS SCIENCE↗