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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A parallel evolutionary multiple-try metropolis Markov chain Monte Carlo algorithm for sampling spatial partitions

We develop an Evolutionary Markov Chain Monte Carlo (EMCMC) algorithm for sampling spatial partitions that lie within a large, complex, and constrained spatial state space. Our algorithm combines the advantages of evolutionary algorithms (EAs) as optimization heuristics for state space traversal and the theoretical convergence properties of Markov Chain Monte Carlo algorithms for sampling from unknown distributions. Local optimality information that is identified via a directed search by our optimization heuristic is used to adaptively update a Markov chain in a promising direction within the framework of a Multiple-Try Metropolis Markov Chain model that incorporates a generalized Metropolis-Hastings ratio. We further expand the reach of our EMCMC algorithm by harnessing the computational power afforded by massively parallel computing architecture through the integration of a parallel EA framework that guides Markov chains running in parallel.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Scaling neural simulations in STACS

Abstract As modern neuroscience tools acquire more details about the brain, the need to move towards biological-scale neural simulations continues to grow. However, effective simulations at scale remain a challenge. Beyond just the tooling required to enable parallel execution, there is also the unique structure of the synaptic interconnectivity, which is globally sparse but has relatively high connection density and non-local interactions per neuron. There are also various practicalities to consider in high performance computing applications, such as the need for serializing neural networks to support potentially long-running simulations that require checkpoint-restart. Although acceleration on neuromorphic hardware is also a possibility, development in this space can be difficult as hardware support tends to vary between platforms and software support for larger scale models also tends to be limited. In this paper, we focus our attention on Simulation Tool for Asynchronous Cortical Streams (STACS), a spiking neural network simulator that leverages the Charm++ parallel programming framework, with the goal of supporting biological-scale simulations as well as interoperability between platforms. Central to these goals is the implementation of scalable data structures suitable for efficiently distributing a network across parallel partitions. Here, we discuss a straightforward extension of a parallel data format with a history of use in graph partitioners, which also serves as a portable intermediate representation for different neuromorphic backends. We perform scaling studies on the Summit supercomputer, examining the capabilities of STACS in terms of network build and storage, partitioning, and execution. We highlight how a suitably partitioned, spatially dependent synaptic structure introduces a communication workload well-suited to the multicast communication supported by Charm++. We evaluate the strong and weak scaling behavior for networks on the order of millions of neurons and billions of synapses, and show that STACS achieves competitive levels of parallel efficiency.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Substrate Partitioning into Protein Macromolecular Frameworks for Enhanced Catalytic Turnover

Spatial partitioning of chemical processes is an important attribute of many biological systems, the effect of which is reflected in the high efficiency of enzymes found within otherwise chaotic cellular environments. Barriers, often provided through the formation of compartments or phase segregation, gate the access of macromolecules and small molecules within the cell and provide an added level of metabolic control. Taking inspiration from nature, we have designed virus-like particles (VLPs) as nanoreactor compartments that sequester enzyme catalysts and have used these as building blocks to construct 3D protein macromolecular framework (PMF) materials, which are structurally characterized using small-angle X-ray scattering (SAXS). Additionally, the highly charged PMFs form a separate phase in suspension, and by tuning the ionic strength, we show positively charged molecules preferentially partition into the PMF, while negatively charged molecules are excluded. This molecular partitioning was exploited to tune the catalytic activity of enzymes enclosed within the individual particles in the PMF, the results of which showed that positively charged substrates had turnover rates that were 8500× faster than their negatively charged counterparts. Moreover, the catalytic PMF led to cooperative behavior resulting in charge dependent trends opposite to those observed with individual P22 nanoreactor particles.

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

Investigating Thaw and Plant Productivity Constraints on Old Soil Carbon Respiration From Permafrost

Isotopic radiocarbon (Δ 14 C) signatures of ecosystem respiration (Reco) can identify old soil carbon (C) loss and serve as an early indicator of permafrost destabilization in a warming climate. Warming also stimulates plant productivity causing plant respiration to dominate Reco Δ 14 C signatures and potentially obscuring old soil C loss. Here, we investigate how a wide spatio-temporal gradient of permafrost thaw and plant productivity affects Reco Δ 14 C patterns and isotopic partitioning. Spatial gradients came from a warming experiment with doubling thaw depth and variable biomass, and a vegetation removal manipulation to eliminate plant contributions. We sampled in August and September to capture transitions from high to low plant productivity, decreased surface soil temperature, and relatively small seasonal thaw extensions. We found that surface processes dominate spatial variation in old soil C loss and a process-based partitioning approach was crucial for constraining old soil C loss. Resampling the same plots in different times of the year revealed that old soil C losses tripled with cooling surface temperature, and the largest old soil C losses were detected when the organic-to-mineral soil horizons thawed (~50–60 cm). We suggest that the measured increase in old soil respiration over the season and when the organic-to-mineral horizon thawed, may be explained by mobilization of nitrogen that stimulates microbial decomposition at depth. Our results suggest that soil C in the organic to mineral horizon may be an important source of soil C loss as the entire Arctic region warms and could lead to nonlinearities in projected permafrost climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The holographic entropy cone from marginal independence

The holographic entropy cone characterizes the relations between entanglement entropies for a spatial partitioning of the boundary spacetime of a holographic CFT in any state describing a classical bulk geometry. We argue that the holographic entropy cone, for an arbitrary number of parties, can be reconstructed from more fundamental data determined solely by subadditivity of quantum entropy. We formulate certain conjectures about graph models of holographic entanglement, for which we provide strong evidence, and rigorously prove that they all imply that such a reconstruction is possible. Our conjectures (except only for the weakest) further imply that the necessary data is remarkably simple. In essence, all one needs to know to reconstruct the holographic entropy cone, is a certain subset of the extreme rays of this simpler “subadditivity cone”, namely those which can be realized in holography. This recasting of the bewildering entanglement structure of geometric states into primal building blocks paves the way to distilling the essence of holography for the emergence of a classical bulk spacetime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗
DOE OSTI2024

Elucidating mechanisms of Photosystem I assembly and repair (Final Technical Report)

Oxygenic photosynthesis generates the carbohydrates and oxygen that support life on earth. Central to this process are photosystems I and II (PSI and PSII), large protein-pigment complexes that are among the most complex macromolecular structures in nature. Impressive advances have been made in understanding the structure and photochemical mechanisms of PSI and PSII. However, less is known about the mechanisms underlying their biogenesis. The overarching goal of this project was to elucidate mechanisms of photosystem assembly. Several proteins that are not found in mature PSI or PSII are essential for their assembly, but the assembly pathways and the role of each assembly factor are just starting to emerge. During the course of this project, we advanced understanding of this problem in several ways. First, we addressed the spatial organization of photosystem assembly by localizing many known assembly factors to specific subdomains of thylakoid membranes in plant chloroplasts. This work revealed that the "curvature" domain - the highly curved periphery of stacked granal membrane disks - is a center of photosystem assembly. However, whereas some PSI assembly factors localized exclusively to the curvature domain, another was found exclusively in granal margins and unstacked membranes. This spatial partitioning suggests an assembly line along which different assembly intermediates are physically segregated. Most assembly factors localized similarly in dark-adapted and illuminated plants, but one PSII assembly factor shifted from the curvature domain to unstacked membranes upon illumination, suggesting an unanticipated role in the repair of light-induced PSII damage. Whereas most PSII subunits were found primarily in grana stacks (as expected), one PSII subunit, PsbP, was prominent also in the curvature fraction, suggesting that PsbP has functions beyond its role in mature PSII, possibly in biogenesis or repair.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1989

Multitasking a three-dimensional Navier-Stokes algorithm on the Cray-2

A three-dimensional computational aerodynamics algorithm has been multitasked for efficient parallel execution on the Cray-2. It provides a means for examining the multitasking performance of a complete CFD application code. An embedded zonal multigrid scheme is used to solve the Reynolds-averaged Navier-Stokes equations for an internal flow model problem. The explicit nature of each component of the method allows a spatial partitioning of the computational domain to achieve a well-balanced task load for MIMD computers with vector-processing capability. Experiments have been conducted with both two- and three-dimensional multitasked cases. The best speedup attained by an individual task group was 3.54 on four processors of the Cray-2, while the entire solver yielded a speedup of 2.67 on four processors for the three-dimensional case. The multiprocessing efficiency of various types of computational tasks is examined, performance on two Cray-2s with different memory access speeds is compared, and extrapolation to larger problems is discussed.

Swisshelm, Julie M.↗
NASA NTRS2010

Geology of the Terra Cimmeria-Utopia Planitia Highland Lowland Transitional Zone: Final Technical Approach and Scientific Results

The southern Utopia highland-lowland transitional zone extends from northern Terra Cimmeria to southern Utopia Planitia and contains broad, bench-like platforms with depressions, pitted cones, tholi, and lobate flows. The locally occurring geologic units and landforms contrast other transitional regions and record a spatially partitioned geologic history. We systematically delineated and described the geologic units and landforms of the southern Utopia-Cimmeria highland-lowland transitional zone for the production of a 1:1,000,000-scale geologic map (MTMs 10237, 15237, 20237, 10242, 15242, 20242, 10247, 15247, and 20247). Herein, we present technical and scientific results of this mapping project.

Skinner, J. A., Jr.↗
NASA NTRS2015

Hybrid Glacier Inventory, Gravimetry and Altimetry (HIGA) Mass Balance Product for Greenland and the Canadian Arctic

We present a novel inversion algorithm that generates a mass balance field that is simultaneously consistent with independent observations of glacier inventory derived from optical imagery, cryosphere-attributed mass trends derived from satellite gravimetry, and ice surface elevation trends derived from airborne and satellite altimetry. We use this algorithm to assess mass balance across Greenland and the Canadian Arctic over the Sep-2003 to Oct- 2009 period at 26 kilometers resolution. We evaluate local algorithm-inferred mass balance against forty in situ point observations. This evaluation yields a root mean squared error (RMSE) of 0.15 mWE/a ( 0.15 meters (water equivalent) per annum), and highlights a paucity of in situ observations from regions of high dynamic mass loss and peripheral glaciers. We assess mass losses of 212 plus or minus 67 Gigatons per annum to the Greenland ice sheet proper, 38 Gigatons per annum to peripheral glaciers in Greenland, and 42 Gigatons per annum to glaciers in the Canadian Arctic. These magnitudes of mass loss are dependent on the gravimetry-derived spherical harmonic mass trend we invert. We spatially partition the transient glacier continuity equation by differencing algorithm inferred mass balance from modeled surface mass balance, in order to solve the horizontal divergence of ice flux as a residual. This residual ice dynamic field infers flux divergence (or submergent flow) in the ice sheet accumulation area and at tidewater margins, and flux convergence (or emergent flow) in land-terminating ablation areas, which is consistent with continuum mechanics theory.

Altimetry↗
DOE OSTI2021

AutoAtlas: Neural Network for 3D Unsupervised Partitioning and Representation Learning

Here we present a novel neural network architecture called AutoAtlas for fully unsupervised partitioning and representation learning of 3D brain Magnetic Resonance Imaging (MRI) volumes. AutoAtlas consists of two neural network components: one neural network to perform multi-label partitioning based on local texture in the volume, and a second neural network to compress the information contained within each partition. We train both of these components simultaneously by optimizing a loss function that is designed to promote accurate reconstruction of each partition, while encouraging spatially smooth and contiguous partitioning, and discouraging relatively small partitions. We show that the partitions adapt to the subject specific structural variations of brain tissue while consistently appearing at similar spatial locations across subjects. AutoAtlas also produces very low dimensional features that represent local texture of each partition. We demonstrate prediction of metadata associated with each subject using the derived feature representations and compare the results to prediction using features derived from FreeSurfer anatomical parcellation. Since our features are intrinsically linked to distinct partitions, we can then map values of interest, such as partition-specific feature importance scores onto the brain for visualization.

42 ENGINEERING↗
NASA NTRS2004

The Influence of Tungsten on the Chemical Composition of a Temporally Evolving Nanostructure of a Model Ni-Al-Cr Superalloy

The influence of W on the temporal evolution of gamma' precipitation toward equilibrium in a model Ni-Al-Cr alloy is investigated by three-dimensional atom-probe (3DAP) microscopy and transmission electron microscopy (TEM). We report on the alloys Ni-10 Al-8.5 Cr (at.%) and Ni-10 Al-8.5 Cr-2 W (at.%), which were aged isothermally in the gamma + gamma' two-phase field at 1073 K, for times ranging from 0.25 to 264 h. Spheroidal-shaped gamma' precipitates, 5-15 nm diameter, form during quenching from above the solvus temperature in both alloys at a high number density (approx. 10(exp 23/cu m). As gamma' precipitates grow with aging at 1073 K, a transition from spheriodal- to cuboidal-shaped precipitates is observed in both alloys. The elemental partitioning and spatially resolved concentration profiles across the gamma' precipitates are obtained as a function of aging time from three-dimensional atom-by-atom reconstructions. Proximity histogram concentration profiles of the quaternary alloy demonstrate that W concentration gradients exist in gamma' precipitates in the as-quenched and 0.25-h aging states, which disappear after 1 h of aging. The diffusion coefficient of W in gamma' is estimated to be 6.2 x 10(exp -20) sq m/s at 1073 K. The W addition decreases the coarsening rate constant, and leads to stronger partitioning of Al to gamma' and Cr to gamma.

Sudbrack, Chantal K.↗
DOE OSTI2022

Variability of Snow and Rainfall Partitioning Into Evapotranspiration and Summer Runoff Across Nine Mountainous Catchments

Abstract Understanding the partitioning of snow and rain contributing to either catchment streamflow or evapotranspiration (ET) is of critical relevance for water management in response to climate change. To investigate this partitioning, we use endmember splitting and mixing analyses based on stable isotope ( 18 O) data from nine headwater catchments in the East River, Colorado. Our results show that one third of the snow partitions to ET and 13% of the snowmelt sustains summer streamflow. Only 8% of the rainfall contributes to the summer streamflow, because most of the rain (67%) partitions to ET. The spatial variability of precipitation partitioning is mainly driven by aspect and tree cover across the sub‐catchments. Catchments with higher tree cover have a higher share of snow becoming ET, resulting in less snow in summer streamflow. Summer streamflow did not contain more rain with higher rainfall sums, but more rain was taken up in ET.

58 GEOSCIENCES↗
DOE OSTI2020

The Case for Strong Scaling in Deep Learning: Training Large 3D CNNs with Hybrid Parallelism

Here, we present scalable hybrid-parallel algorithms for training large-scale 3D convolutional neural networks. Deep learning-based emerging scientific workflows often require model training with large, high-dimensional samples, which can make training much more costly and even infeasible due to excessive memory usage. We solve these challenges by extensively applying hybrid parallelism throughout the end-to-end training pipeline, including both computations and I/O. Our hybrid-parallel algorithm extends the standard data parallelism with spatial parallelism, which partitions a single sample in the spatial domain, realizing strong scaling beyond the mini-batch dimension with a larger aggregated memory capacity. We evaluate our proposed training algorithms with two challenging 3D CNNs, CosmoFlow and 3D U-Net. Our comprehensive performance studies show that good weak and strong scaling can be achieved for both networks using up to 2K GPUs. More importantly, we enable training of CosmoFlow with much larger samples than previously possible, realizing an order-of-magnitude improvement in prediction accuracy.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

Sources of Variability in the Stratospheric Column of Nitrogen Dioxide

Instruments such as the Global Ozone Monitoring Experiment (GOME, on the European Remote Sensing Satellite (ERS-2), launched 1995), the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY, on ENVISAT, to be launched July 2001) and the Ozone Monitoring Instrument (OMI on EOS Aura, to be launched 2003) make measurements of the total column of NO2. There is interest in separating the stratospheric and tropospheric contributions to the column, as the tropospheric column provides a measure of pollution. We are using a 3D chemistry and transport model driven by winds from the Goddard Space Flight Center Data Assimilation System to examine variability in the stratospheric NO2 column. Model results for NOx = NO + NO2 + 2N2O5 will be shown to compare well with sunset observations from the Halogen Occultation Experiment on the Upper Atmosphere Research Satellite, and to exhibit similar temporal and spatial dependence. Partitioning between NO, NO2, and N2O5 is also shown to compare well with observations. This good agreement supports the use of simulated fields in the stratosphere to derive the tropospheric column from the total column. Preliminary comparisons of the tropospheric column with model simulations for the troposphere will also be shown.

Douglass, Anne↗
DOE OSTI2025

Spatially Aware Linear Transformer (SAL-T) for Particle Jet Tagging

Transformers are very effective in capturing both global and local correlations within high-energy particle collisions, but they present deployment challenges in high-data-throughput environments, such as the CERN LHC. The quadratic complexity of transformer models demands substantial resources and increases latency during inference. In order to address these issues, we introduce the Spatially Aware Linear Transformer (SAL-T), a physics-inspired enhancement of the linformer architecture that maintains linear attention. Our method incorporates spatially aware partitioning of particles based on kinematic features, thereby computing attention between regions of physical significance. Additionally, we employ convolutional layers to capture local correlations, informed by insights from jet physics. In addition to outperforming the standard linformer in jet classification tasks, SAL-T also achieves classification results comparable to full-attention transformers, while using considerably fewer resources with lower latency during inference. Experiments on a generic point cloud classification dataset (ModelNet10) further confirm this trend. Our code is available at https://github.com/aaronw5/SAL-T4HEP.

Wang, Aaron [Illinois U., Chicago] (ORCID:00000003↗
NASA NTRS1995

Intelligent automated surface grid generation

The goal of our research is to produce a flexible, general grid generator for automated use by other programs, such as numerical optimizers. The current trend in the gridding field is toward interactive gridding. Interactive gridding more readily taps into the spatial reasoning abilities of the human user through the use of a graphical interface with a mouse. However, a sometimes fruitful approach to generating new designs is to apply an optimizer with shape modification operators to improve an initial design. In order for this approach to be useful, the optimizer must be able to automatically grid and evaluate the candidate designs. This paper describes and intelligent gridder that is capable of analyzing the topology of the spatial domain and predicting approximate physical behaviors based on the geometry of the spatial domain to automatically generate grids for computational fluid dynamics simulators. Typically gridding programs are given a partitioning of the spatial domain to assist the gridder. Our gridder is capable of performing this partitioning. This enables the gridder to automatically grid spatial domains of wide range of configurations.

Yao, Ke-Thia↗
NASA NTRS2022

Quantifying Spatial Drought Propagation Potential in North America Using Complex Network Theory

Droughts have a dominant three-dimensional (3-D) spatiotemporal structure typically spanning hundreds of kilometers and often lasting for months to years. Here, we introduced a novel framework to explore the 3-D structure of the evolution of droughts based on network theory concepts. The proposed framework is applied to identify critical source regions responsible for large-scale drought onsets during 1901–2014 for the North American continent using the Standardized Precipitation Evaporation Index (SPEI). We built a spatial network connecting the drought onset timings for the North American continent. Using a spatially weighted network partitioning algorithm, the whole continent is then classified into regional spatial drought networks (RSN), where droughts are more likely to propagate within these regional systems. Finally, a customized network metric was applied to identify locations (source regions) where the drought onsets further propagate to other areas within the regional spatial network. Our results indicated that the West coast, Texas coastal region, and Southeastern Arkansas as major source regions through which atmospheric drought propagates to Western, South Central, and Eastern North America. The formation of drought source regions are due to presence of high pressure ridges and anomalous wind patterns. Furthermore, our results indicate that the drought propagation from these source regions may be due to inadequate moisture transport. The proposed framework can help to develop an early warning detection system for droughts and other spatially extensive extreme events such as heatwaves and floods.

Goutam Konapala↗