Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial Partitioning”

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 OSTI2025

Assessing Clouds in GFDL's AM4.0 With Different Microphysical Parameterizations Using the Satellite Simulator Package COSP

We evaluate cloud simulations using satellite simulators against multiple observational data sets. These simulators have been run within the Geophysical Fluid Dynamics Laboratory's Atmosphere Model version 4.0 (AM4.0), as well as an alternative configuration where a fully two‐moment Morrison‐Gettelman cloud microphysical parameterization with prognostic precipitation (MG2) is applied, denoted as AM4‐MG2. The modeled cloud spatial distributions, vertical profiles, phase partitioning, cloud‐to‐precipitation transitions, and radiative effects compare reasonably well with satellite observations. Model biases include the under‐prediction of total and low‐level clouds, especially optically thin/intermediate clouds with cloud optical depth of less than 23, but the over‐prediction of thick clouds, indicating “too few, too bright” biases. These biases counteract each other, and give rise to reasonable estimates of cloud radiative effects. The underestimate of low‐level clouds is associated with too early and too frequent drizzle/precipitation formation. The precipitation bias is improved in AM4‐MG2, where the autoconversion scheme initiates the precipitation more realistically. There also exist discrepancies between models and observations for midlevel and high‐level clouds. Additional biases include the underestimate of liquid cloud fraction and the overestimate of ice cloud fraction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Seismoacoustic Analysis of the Large Surface Explosion Coupling Experiment Using a Large- N Seismic Array

Seismoacoustic signals at local distance (<~10 km) are widely used as important constraints on source parameters for near-surface events, yet the seismoacoustic wave generation and energy partitioning are not fully understood. Spatially dense sensors could provide observations in high resolution to capture the full wavefield for better understanding wave propagation and improving source estimation. Recently, spatially dense observations of the local seismoacoustic wavefield produced by a pair of 1-ton surface explosions have been recorded using a large-N seismic array. This large-N array consists of 446 geophones and covers an area of approximately 2×2.5 km 2 . The two surface explosions occur at the same location but at different times with different atmospheric conditions. Both seismic and air–ground coupled acoustic waves from the two surface explosions are well observed. Analyses of signals recorded by the large-N seismic array show different acoustic wave speed and amplitude for the two explosions. A strong spatial variability in acoustic wave speed and amplitude for each explosion is also observed. Finally, the observations suggest the important role of local atmosphere state on wave propagation and source estimation and demonstrate how the use of a large-N capability can improve characterization of the propagation medium and source.

58 GEOSCIENCES↗
DOE OSTI2023

A seamless approach for evaluating climate models across spatial scales

In regions of the world where topography varies significantly with distance, most global climate models (GCMs) have spatial resolutions that are too coarse to accurately simulate key meteorological variables that are influenced by topography, such as clouds, precipitation, and surface temperatures. One approach to tackle this challenge is to run climate models of sufficiently high resolution in those topographically complex regions such as the North American Regionally Refined Model (NARRM) subset of the Department of Energy’s (DOE) Energy Exascale Earth System Model version 2 (E3SM v2). Although high-resolution simulations are expected to provide unprecedented details of atmospheric processes, running models at such high resolutions remains computationally expensive compared to lower-resolution models such as the E3SM Low Resolution (LR). Moreover, because regionally refined and high-resolution GCMs are relatively new, there are a limited number of observational datasets and frameworks available for evaluating climate models with regionally varying spatial resolutions. As such, we developed a new framework to quantify the added value of high spatial resolution in simulating precipitation over the contiguous United States (CONUS). To determine its viability, we applied the framework to two model simulations and an observational dataset. We first remapped all the data into Hierarchical Equal-Area Iso-Latitude Pixelization (HEALPix) pixels. HEALPix offers several mathematical properties that enable seamless evaluation of climate models across different spatial resolutions including its equal-area and partitioning properties. The remapped HEALPix-based data are used to show how the spatial variability of both observed and simulated precipitation changes with resolution increases. This study provides valuable insights into the requirements for achieving accurate simulations of precipitation patterns over the CONUS. It highlights the importance of allocating sufficient computational resources to run climate models at higher temporal and spatial resolutions to capture spatial patterns effectively. Furthermore, the study demonstrates the effectiveness of the HEALPix framework in evaluating precipitation simulations across different spatial resolutions. This framework offers a viable approach for comparing observed and simulated data when dealing with datasets of varying spatial resolutions. By employing this framework, researchers can extend its usage to other climate variables, datasets, and disciplines that require comparing datasets with different spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Mixing of the Connecticut River Plume During Ambient Flood Tides: Spatial Heterogeneity and Contributions of Bottom–Generated and Interfacial Mixing

The Connecticut River plume is influenced by energetic ambient tides in the Long Island Sound receiving waters. The objectives of this modeling study are (a) characterizing the spatial heterogeneity of turbulent buoyancy fluxes, (b) partitioning turbulent buoyancy fluxes into bottom-generated and interfacial shear contributions, and (c) quantifying contributions to plume-integrated mixing within the tidal plume. The plume formed during ambient flood tides under low river discharge, spring tides, and no winds is analyzed. Turbulent buoyancy fluxes (B) and depth-integrated B through the plume (Bd) are characterized by pronounced spatial heterogeneity. Strong mixing (Bd ~ 10 –5 -10 –4 m 3 /s 3 ) occurs near the mouth, in the nearfield plume turning region, over shoals, and nearshore shallow areas. Low to moderate mixing (Bd ~ 10 –8 -10 –6 m 3 /s 3 ) occupies half the plume. Buoyancy fluxes are first partitioned based on the depth of the shear stress minimum between plume-generated and bottom-generated shear maxima. Four other tested partitioning methods are based on open channel flow and stratified shear flow parameterizations. Interfacial and bottom-generated shear contribute to different areas of intense and moderate mixing. All methods indicate a significant plume mixing role for bottom-generated mixing, but interfacial mixing is a bigger contributor. Plume-integrated total and interfacial mixing peak at max ambient flood and the timing of peak bottom-generated mixing varies among partitioning methods. Two-thirds of the mixing occurs in concentrated intense mixing areas. A parameter space with the ambient tidal Froude number and plume thickness to depth ratio as axes indicates many tidally modulated plumes are moderately to dominantly influenced by bottom-generated tidal mixing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2026

The biogeography of soil and airborne fungi in the Southwestern USA in relation to climate and vegetation

To assess how fungal dispersal might respond to climate change, we examined how climate and geography influence the regional distribution of fungi in soil and air. Specifically, we hypothesized that neighboring fungal communities should be more similar than distant communities (i.e. spatially autocorrelated) and that fungal dispersal should be more limited in soil than in air. We collected soil and air samples from 60 sites across five states in the Southwestern USA. Then, we sequenced the ITS2 region to identify fungal taxa in each sample. Next, we used distance-based redundancy analysis to partition variation in fungal community composition between climate variables and spatial structure. Fungi were indeed spatially autocorrelated. Moreover, precipitation, maximum vapor pressure deficit, and soil moisture were significantly related to fungal community composition in soils. In comparison, only precipitation was significantly related to community composition in the air. After accounting for climate, the strength of spatial autocorrelation did not differ significantly in soilborne versus airborne fungi. Dispersal limitation was evident in soilborne fungi at short distances (<100 km) and was not observed at any distance in airborne fungi. Altogether, climate may influence which fungal taxa are present in soil and air, and fungi could feasibly wind disperse over regional scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effects of exogenous citric acid on the concentration and spatial distribution of Ni, Zn, Co, Cr, Mn and Fe in leaves of Noccaea caerulescens grown on a serpentine soil

The aim of this study was to show the potential of citric acid in increasing the concentration of Ni, Zn, Co, Cr, Mn and Fe in leaves of the hyperaccumulator Noccaea caerulescens. Synchrotron x-ray fluorescence (μ-XRF) images were collected to assess the distribution of metals in leaves. Applying citric acid (20 mmol kg -1 ) to soil increased in 14-, 10-, 7-, 2- and 1.4- fold the concentration of Mn, Fe, Co, Ni, and Cr, respectively, compared to the control. The μ-XRF imaging revealed that Ni and Zn were not spatially correlated across the leaf. We observed a clear partitioning of Zn between veins and surrounding leaf cells while Ni was more evenly distributed between veins and leaf blade. The accumulation of metals in citric acid treated plants did not change the Ni and Zn distribution pattern in leaves but altered the Mn distribution. It seems that Mn reached toxic concentrations in leaves and we hypothesize that a mechanism driven by transpiration through the xylem was used to excrete the metal. Our results show that citric acid can enhance metal accumulation by N. caerulescens and have impact for soil remediation by either decreasing the time for clean up or increasing the access to non-labile pools of metals in soil.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Probing Nanoscale Chemical Environments of Zinc in Diatoms

Diatoms impact the biogeochemical cycling of zinc (Zn) due to elevated cellular Zn levels over other phytoplankton, and their strong biological uptake of Zn in the Southern Ocean sets the global distribution of Zn in ocean water columns. Past studies have revealed the abundance and spatial distribution of Zn in individual diatoms, while others have shown diverse Zn chemical species in whole-cell aggregates. However, intracellular and intercellular variations in Zn chemical forms remained unknown. For the first time, we applied a synchrotron X-ray nanoprobe to cultured diatoms (Phaeodactylum tricornutum and Chaetoceros muelleri) and resolved spatial heterogeneities of Zn chemistry. The result shows the partitioning of Zn between intracellular contents (soft parts) and frustules (mainly in the intracellular contents for P. tricornutum and similar Zn concentrations in the intracellular contents and frustules for C. muelleri). We found multiple Zn chemical species unevenly distributed within individual cells, where Zn–phosphoryl complexes were the most abundant followed by cysteine, histidine, biogenic silica, and carboxyl complexes. Cellular Zn species varied between individual organisms and were influenced by Zn availability. To our knowledge, this work presents the first X-ray measurements of trace metal speciation at sub-100 nm resolution in biological samples in their natural state. In conclusion, this research improves our ability to examine the biogeochemistry of Zn on the nanoscale and can help us understand the role of Zn in plankton growth and the Zn cycle in ocean waters.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Unexpected Oceanic Peak in Energy Input to the Atmosphere and Its Consequences for Monsoon Rainfall

Abstract Monsoons have historically been understood to be caused by the low thermal inertia of land, allowing more energy from summer insolation to be transferred to the overlying atmosphere than over adjacent ocean. Here, we show that during boreal summer, the global maximum net energy input (NEI) to the atmosphere unexpectedly lies over the Indian Ocean, not over land. Observed radiative fluxes suggest that cloud‐radiative effects (CRE) almost double the NEI over ocean, shifting the NEI peak from land to ocean. Global climate model experiments with both land and interactive sea surface temperatures confirm that CRE create the oceanic NEI maximum. Interactions between CRE, NEI, circulation, and land‐sea contrast in surface heat capacity shift precipitation from Southeast to South Asia. CRE thus alter the global partitioning of precipitation between land and ocean and the spatial structure of Earth's strongest monsoon, in ways that can be understood through the NEI.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Uranium and Nickel Partitioning in a Contaminated Riparian Wetland

Uranium (U) and nickel (Ni) released 50 years ago have been immobilized in the Tims Branch wetlands located on the Savannah River Site in the United States. Sediments were collected from seven locations to identify the factors responsible for this attenuation. Ni and U contents in the solids were significantly correlated, suggesting that depositional as opposed to chemical processes contributed to their spatial distribution. Based on sequential extractions, 63 ± 16% of the U was partitioned into the organic fraction, whereas Ni was distributed between several sediment fractions. An inverse pH-organic matter (OM) correlation and positive correlations of OM with total U and organic-bound U/Ni suggest that increased OM preservation and binding to the mineral surfaces were likely responsible for Ni- and especially U-sediment retention (Tims Branch pH = 4.84 ± 0.68). EXAFS analysis indicated the predominance of U(VI) coordinated with clay minerals (~65%), together with ~35% coordinated to either OM (in areas with elevated OM levels) or iron oxides. The desorption-Kd coefficients of U (3972 ± 1370 L/kg) and Ni (30 ± 8 L/kg) indicate that dissolved Ni poses a greater long-term risk than dissolved U for migrating downstream. This study suggests that a delicate balance of geochemical properties controls whether wetlands behave as sinks or sources of contaminants.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2026

Many-Body Benchmark of Electronic Charge and Spin Densities for Li 1–x NiO 2

Accurate benchmarks are particularly important for highly correlated oxides as mean-field approximations often fail to describe the subtle balance of charge transfer and magnetism in these materials with an accuracy comparable to experimental needs. Here we present accurate diffusion Monte Carlo (DMC) results of the electronic charge and spin densities for the tunable highly correlated oxide Li 1–x NiO 2 for x = 0, 1/2, and 1. To enable quantitative comparisons, we introduce a robust density-partitioning scheme, extending Voronoi analysis to assign atomic charges from spatially noisy DMC densities. We then benchmark common approximations used in density functional theory (DFT). Comparison against DMC shows that r 2 SCAN delivers the most balanced performance across charge, spin, and radial density descriptors, nearly reproducing DMC results for LiNiO 2 and apical Ni sites in Li 0.5 NiO 2 . Hybrid functionals (PBE0, SCAN0) perform unexpectedly poorly, and PBE + U + V yields inconsistent trends between charge and spin densities. Therefore, the r 2 SCAN functional minimizes errors relative to DMC while capturing the variable valence of the Ni ion and also retaining the computational efficiency of DFT for large-scale simulations of the tunable structural and electronic phases of Li1−xNiO2. Our study highlights the importance of accurate benchmarking of the fundamental quantities involved in DFT to select appropriate DFT approximations in order to advance the predictive modeling of charge-transfer-driven phenomena in correlated electron systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Poisson tensor completion non-parametric differential entropy estimator

We introduce the Poisson tensor completion (PTC) estimator, a non-parametric differential entropy estimator. The PTC estimator leverages inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram. Our crucial observation is that the histogram bins are an instance of a space partitioning of counts and thus can be identified with a spatial Poisson process. The Poisson tensor decomposition leads to a completion of the intensity measure over all bins—including those containing few to no samples—and leads to our proposed PTC differential entropy estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values and so can be safely used directly in entropy estimation. Our estimator is the first tensor-based estimator that exploits the underlying spatial Poisson process related to the histogram explicitly when estimating the probability density with low-rank tensor decompositions for the purpose of tensor completion. Furthermore, we demonstrate that our PTC estimator is a substantial improvement over standard histogram-based estimators for sub-Gaussian probability distributions because of the concentration of norm phenomenon.

42 ENGINEERING↗
DOE OSTI2026

The Poisson tensor completion parametric estimator

We introduce the Poisson tensor completion (PTC) estimator that exploits inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram for samples of a multivariate distribution. Our crucial observation is that the histogram bins are an instance of a space partitioning of counts and thus can be identified with a spatial non-homogeneous Poisson process. The Poisson tensor decomposition leads to a completion of the mean measure over all bins—including those containing few to no samples—and leads to our proposed estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values obviating the need for additional constraints to ensure non-negativity. Furthermore, we demonstrate that our PTC estimator is a substantial improvement over standard histogram-based estimators for sub-Gaussian probability distributions because of the concentration of norm phenomenon.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A decreasing carbon allocation to belowground autotrophic respiration in global forest ecosystems

Belowground autotrophic respiration (RAsoil) depends on carbohydrates from photosynthesis flowing to roots and rhizospheres, and is one of the most important but least understood components in forest carbon cycling. Carbon allocation plays an important role in forest carbon cycling and reflects forest adaptation to changing environmental conditions. However, carbon allocation to RAsoil has not been fully examined at the global scale. To fill this knowledge gap, first, the spatio-temporal patterns of RAsoil from 1981 to 2017 were predicted by a Random Forest (RF) algorithm using the most updated Global Soil Respiration Database (v5) with global environmental variables; second, carbon allocation from photosynthesis to RAsoil (CAB), was calculated as the ratio of RAsoil to gross primary production; and its temporal and spatial patterns were assessed in global forest ecosystems. . Globally, mean RAsoil from forests was 8.9 ± 0.08 Pg C yr-1 (mean ± standard deviation) from 1981 to 2017 with strong spatial variabilities. Temporally, RAsoil increased at a rate of 0.0059 Pg C yr-2, paralleling broader soil respiration changes and indicating increasing carbon respired by roots. Mean CAB was 0.243 ± 0.016 and decreased over time. The temporal trend of CAB varied greatly in space, reflecting uneven responses of CAB to environmental changes. This study is the first attempt to predict global CAB and analyze its temporal and spatial patterns. With the linkage of carbon use efficiency, the developed CAB offers an completely independent approach to quantify global aboveground autotropic respiration spatially and temporally, which could provide crucial insights into carbon flux partition and global carbon cycling under climate change.

Tang, Xiaolu↗
DOE OSTI2022

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2025

Advancing Organized Convection Representation in the Unified Model: Implementing and Enhancing Multiscale Coherent Structure Parameterization

To address the effect of stratiform latent heating on meso- to large-scale circulations, an enhanced implementation of the Multiscale Coherent Structure Parameterization (MCSP) is developed for the Met Office Unified Model. MCSP represents the top-heavy stratiform latent heating from under-resolved organized convection in general circulation models. We couple the MCSP with a mass-flux convection scheme (CoMorph-A) to improve storm lifecycle continuity. The improved MCSP trigger is specifically designed for mixed-phase deep convective cloud, combined with a background vertical wind shear, both known to be crucial for stratiform development. We also test a cloud top temperature dependent convective-stratiform heating partitioning, in contrast to the earlier fixed partitioning. Assessments from ensemble weather forecasts and decadal simulations demonstrate that MCSP directly reduces cloud deepening and precipitation areas by moderating mesoscale circulations. Indirectly, it amends tropical precipitation biases, notably correcting dry and wet biases over India and the Indian Ocean, respectively. Remarkably, the scheme outperforms a climate model ensemble by improving seasonal precipitation cycle predictions in these regions. The scheme also improves Madden-Julian Oscillation (MJO) spectra, achieving better alignment with observational and reanalysis data by intensifying the simulated MJO over the Indian Ocean during phases 4 to 5. However, the scheme increases precipitation overestimation over the Western Pacific. Shifting from fixed to temperature-dependent convective-stratiform partitioning reduces the Pacific precipitation overestimation and further improves the seasonal cycle in India. Spatially correlated biases highlight the necessity for advances beyond deterministic approaches to align MCSP with environmental conditions.

54 ENVIRONMENTAL SCIENCES↗