Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial optimization”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Understanding the Design Space of Sparse/Dense Multiphase Dataflows for Mapping Graph Neural Networks on Spatial Accelerators

Graph Neural Networks (GNNs) have garnered a lot of recent interest because of their success in learning representations from graph-structured data across several critical applications in cloud and HPC. Owing to their unique compute and memory characteristics that come from an interplay between dense and sparse phases of computations, the emergence of reconfigurable dataflow (aka spatial) accelerators offers promise for acceleration by mapping optimized dataflows (i.e., computation order and parallelism) for both phases. The goal of this work is to characterize and understand the design-space of dataflow choices for running GNNs on spatial accelerators in order for the compilers to optimize the dataflow based on the workload. Specifically, we propose a taxonomy to describe all possible choices for mapping the dense and sparse phases of GNNs spatially and temporally over a spatial accelerator, capturing both the intra-phase dataflow and the inter-phase (pipelined) dataflow. Using this taxonomy, we do deep-dives into the cost and benefits of several dataflows and perform case studies on implications of hardware parameters for dataflows and value of flexibility to support pipelined execution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-fidelity holographic beam shaping with optimal transport and phase diversity

A phase-only spatial light modulator (SLM) provides a powerful way to shape laser beams into arbitrary intensity patterns but at the cost of a hard computational problem of determining an appropriate SLM phase. Here, we show that optimal transport methods can generate approximate solutions to this problem that serve as excellent initializations for iterative phase retrieval algorithms, yielding vortex-free solutions with superior accuracy and efficiency. Additionally, we show that analogous algorithms can be used to measure the intensity and phase of the input beam incident upon the SLM via phase diversity imaging. These techniques furnish flexible and convenient solutions to the computational challenges of beam shaping with an SLM.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Biorefinery siting and sizing to achieve the US Billion‐Ton Bioeconomy vision: A case study using a gasification–Fischer–Tropsch process

Achieving a secure, abundant, and affordable energy future requires a robust and adaptable energy strategy, with bioenergy playing a pivotal role. Biomass-based energy presents a promising pathway to use domestic resources while fostering economic opportunities in rural areas. Despite the potential to source more than 1 billion dry short tons of biomass annually in the US, significant infrastructure and economic barriers hinder full utilization for energy production. This study used the Biofuel Infrastructure, Logistics, and Transportation (BILT) model to assess biorefinery siting and scale and determine the number and size of facilities required to maximize use of the US biomass potential. A spatially agnostic approach first assessed the effects of facility capacity and transportation constraints on biomass use. Then, a spatially explicit analysis integrated county-level biomass availability from the US Department of Energy's 2023 Billion-Ton Report and technoeconomic assessments to evaluate different biorefinery deployment scenarios. The results indicate that an optimized mix of facility sizes is essential to leverage biomass resources fully across varying regional production densities to maximize use of the US biomass potential. Larger biorefineries or co-located smaller facilities significantly enhance biomass use while reducing costs through economies of scale. These findings underscore the importance of strategically balancing facility capacity and spatial distribution to optimize the bioenergy supply chain. In conclusion, this study provides critical insights for advancing the US bioenergy economy by aligning biorefinery deployment with biomass resource availability and economic viability.

BILT Model↗
DOE OSTI2021

Computational design of microarchitected porous electrodes for redox flow batteries

Porous electrodes are used as the core reactive component across electrochemical technologies. In flowing systems, controlling the fluid distribution, species transport, and reactive environment is critical to attaining high performance. However, conventional electrode materials like felts and papers provide few opportunities for precise engineering of the electrode and its microstructure. To address these limitations, architected electrodes composed of unit cells with spatially varying geometry determined via computational optimization are proposed. Resolved simulation is employed to develop a homogenized description of the constituent unit cells. These effective properties serve as inputs to a continuum model for the electrode when used in the negative half-cell of a vanadium redox flow battery. Porosity distributions minimizing power loss are then determined via computational design optimization to generate architected porosity electrodes. The architected electrodes are compared to bulk, uniform porosity electrodes and found to lead to increased power efficiency across operating flow rates and currents. The design methodology is further used to generate a scaled-up electrode with comparable power efficiency to the bench-scale systems. Finally, the variable porosity architecture and computational design methodology presented here thus offers a novel pathway for automatically generating spatially engineered electrode structures with improved power performance.

25 ENERGY STORAGE↗
DOE OSTI2025

Complexity Reduction Methods for Large-Scale Spatially Explicit Biofuels Network Design

The size and complexity of energy system optimization models have increased significantly in recent years, driven by the availability of high-resolution spatial data. We present complexity reduction and solution methods that enable us to efficiently represent high-resolution spatial data in the network design of large-scale energy systems. We aim to reduce the size and enhance the computational efficiency of network design models without sacrificing solution accuracy. Specifically, we first present how to aggregate highly granular data into larger resolutions without averaging out their specific properties through a composite-curve-based approach and then develop a method to linearly represent these curves. Second, we utilize a general clustering method to determine groups of geographically proximate biomass fields and establish a single transportation arc for all of them, reducing the number of transportation-related variables while maintaining an accurate representation of the system. Finally, we introduce a two-step algorithm that decomposes large-scale network design problems into two smaller, more manageable subproblems. We demonstrate the application of our methods using a case study of switchgrass-to-biofuels network design in the eight states of the U.S. Midwest, using realistic and highly explicit spatial data.

09 BIOMASS FUELS↗
DOE OSTI2020

Laser-induced digital oxidation for copper-based flexible photodetectors

Copper oxide compounds (Cu x O) with bandgaps of 1.3–2.1 eV (CuO) and 2.1–2.6 eV (Cu 2 O) have been investigated as promising p-type semiconducting materials. Cu x O is generally obtained by deposition or thermal oxidation, but those methods are not optimal for flexible substrates. Furthermore, additional patterning steps are required to fabricate devices. In this work, we present an easy, controllable method to fabricate a metal-semiconductor-metal (MSM) photodetector using laser-induced oxidation of a thin Cu film. After laser irradiation, the Cu film is heated under ambient conditions, and this leads to a thermal oxidation reaction, in which Cu oxide is monolithically formed in the Cu film and a Cu-Cu x O-Cu MSM structure is produced. Since the laser offers localized heating, an arbitrary CuxO pattern can be written in the Cu film by spatially controlled heating. In addition, by optimizing the heating time, the laser-induced oxidation can be successfully performed even on a flexible substrate. To study the laser-induced oxidation, we examined the correlation between laser parameters and the oxidation pattern and analyzed the composition using scanning electron microscopy, Raman spectroscopy, and X-ray photoemission spectroscopy. Furthermore, we measured the transient photoresponse and employed scanning photocurrent microscopy to investigate the mechanism of carrier transport behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Alchemy: A Model-Based Approach for 2D to 3D Autonomous Nuclear System Design

Engineering design of nuclear power plant (NPP) piping and equipment systems frequently bypasses crucial 2D system planning, instead moving straight to 3D modeling. This often leads to designs that exceed building envelope constraints, forcing expensive and time-consuming redesigns. When 2D modeling is employed, it typically involves labor-intensive manual workflows that convert 2D drawings into 3D models, resulting in inefficiencies and errors across design iterations. These workflows further suffer from poor software interoperability and dependence on proprietary software ecosystems, thereby contributing to schedule delays and cost overruns. This paper presents Alchemy, an autonomous framework that transforms 2D system definitions into Industry Foundation Classes (IFC)-compliant 3D building information models (BIMs) for expediting nuclear facility design at the conceptual preliminary phase. Using a model-based approach, the framework treats the 2D system diagram as the central reference model employed to automatically generate all subsequent outputs, ensuring consistency between the system definition and the resulting physical design. A web-based interface enables engineers to define hierarchical system topologies including associated equipment, geometric properties, and connectivity requirements. A two-phase equipment layout optimization algorithm automatically computes collision-free spatial configurations within predefined building envelopes. An artificial intelligence (AI)-assisted pipe routing module then generates orthogonal, collision-free routing paths, allowing the user to select either an A* search-based method or an Ant Colony Optimization (ACO)-based method. All outputs are authored natively in IFC format, relying on open-source technologies and standardized formats in order to ensure extensibility and eliminate proprietary software dependencies. The proposed framework is validated on two representative pressurized-water reactor (PWR)-based case studies, for which it autonomously generates IFC-compliant 3D models in minutes, drastically reducing workflows that typically require hours of manual effort. The generated model demonstrates topologically correct equipment placement, physically plausible spatial relationships, and collision-free pipe routing consistent with known PWR loop configurations. This work represents a foundational step toward digital engineering for nuclear facility preliminary design, with future ongoing development targeting design code compliance and expanded system complexity.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Sequence of polyhedral relaxations for nonlinear univariate functions

Here, given a nonlinear, univariate, bounded, and differentiable function f(x), this article develops a sequence of Mixed Integer Linear Programming (MILP) and Linear Programming (LP) relaxations that converge to the graph of f(x) and its convex hull, respectively. Theoretical convergence of the sequence of relaxations to the graph of the function and its convex hull is established. For nonlinear non-convex optimization problems, the relaxations presented in this article can be used to construct tight MILP and LP relaxations. These MILP and the LP relaxations can also be used with MILP-based and spatial branch-and-bound based global optimization algorithms, respectively.

42 ENGINEERING↗
DOE OSTI2024

Visual Analytics of Performance of Quantum Computing Systems and Circuit Optimization

Driven by potential exponential speedups in business, security, and scientific scenarios, interest in quantum computing is surging. This interest feeds the development of quantum computing hardware, but several challenges arise in optimizing application performance for hardware metrics (e.g., qubit coherence and gate fidelity). In this work, we describe a visual analytics approach for analyzing the performance properties of quantum devices and quantum circuit optimization. Our approach allows users to explore spatial and temporal patterns in quantum device performance data and it computes similarities and variances in key performance metrics. Detailed analysis of the error properties characterizing individual qubits is also supported. We also describe a method for visualizing the optimization of quantum circuits. The resulting visualization tool allows researchers to design more efficient quantum algorithms and applications by increasing the interpretability of quantum computations.

Chae, Junghoon↗
DOE OSTI2024

Detection and Signal Processing for Near-Field Nanoscale Fourier Transform Infrared Spectroscopy

Researchers from a broad spectrum of scientific and engineering disciplines are increasingly using scattering-type near-field infrared spectroscopic techniques to characterize materials non-destructively with nanoscale spatial resolution. However, a sub-optimal understanding of a technique's implementation can complicate data interpretation and act as a barrier to entering the field. Here the key detection and processing steps involved in producing scattering-type near-field nanoscale Fourier transform infrared spectra (nano-FTIR) are outlined. The self-contained mathematical and experimental work derives and explains: i) how normalized complex-valued nano-FTIR spectra are generated, ii) why the real and imaginary components of spectra qualitatively relate to dispersion and absorption respectively, iii) a new and generally valid equation for spectra which can be used as a springboard for additional modeling of the scattering processes, and iv) an algebraic expression that can be used to extract an approximation to the sample's local extinction coefficient from nano-FTIR. The algebraic model for weak oscillators is validated with nano-FTIR and attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectra on samples of polystyrene and Kapton and further provides a pedagogical pathway to cementing some of the technique's key qualitative attributes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

4D Imaging of ZnO-Coated Nanoporous Al 2 O 3 Aerogels by Chemically Sensitive Ptychographic Tomography: Implications for Designer Catalysts

The 3D chemical structure (4D spectromicroscopy) of nanoporous Al 2 O 3 aerogels coated with ZnO by atomic layer deposition (ALD) was studied by multienergy scanning transmission X-ray microscopy. These materials are representative of a class of designer catalysts in which the nanoporous support is prepared separately from the active catalyst material, which is subsequently introduced by ALD, thereby allowing independent optimization of the morphology, chemistry, and spatial distribution of the support and catalyst. The samples studied were prepared by Ga ion and Xe plasma focused ion beam (FIB) milling as well as drop casting from water suspension. Zn L and Al K edge spectra of six samples with three different ZnO loadings were measured to investigate how loading and different sample preparation methods affect the 3D distribution of the ZnO and Al 2 O 3 . Scanning transmission X-ray microscopy (STXM) and ptychographic imaging at two energies each at the Zn L 3 and Al K edge were measured. The ptychography data were analyzed by using the SHARP reconstruction code to generate quantitative 2D chemical maps of the Al 2 O 3 and the ZnO. The STXM and ptychography maps were then measured at a sequence of tilt angles, covering up to 160° of rotation. The 3D structure of the ZnO and Al 2 O 3 was derived from the tilt series data by tomographic reconstruction using a compressed sensing algorithm. A twodimensional spatial resolution (half-period) of 6 nm, measured by Fourier ring correlation, and a 3D spatial resolution (half-period) of 9 nm, measured by Fourier shell correlation, were achieved when using the COSMIC beamline at the Advanced Light Source (ALS). The results show that for all of the ZnO loadings studied there is nonuniform coverage of the ZnO on the Al 2 O 3 aerogel framework. In addition, we found that both FIB methods create sample artifacts, although the distortion was less with Xe plasma than Ga ion FIB.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Middle to Late Holocene Sea Surface Temperature and Productivity Changes in the Northeast Pacific

Variations of the sea surface temperature (SST) and primary productivity in the northeast Pacific have far‐reaching implications. In addition to influencing the regional and global temperature and hydroclimate, these conditions also control marine ecosystems and their services, which subsequently impact regional economies. Yet, our understanding of the variability and controls of northeast Pacific SST and productivity on timescales exceeding observational records remains limited. Here, we use marine sediment records from seven locations, spanning 25.2°N–59.6°N, in the northeast Pacific to characterize the millennial‐scale variability of SST and productivity from 9,000 to 1,000 years BP. We explore the dynamics of their spatiotemporal evolution and compare these data with transient climate model outputs to identify potential drivers. Through a heat budget analysis and optimal fingerprinting analysis, we characterize the spatial pattern of forcings. We find that SST varied spatially in the northeast Pacific, with higher latitudes exhibiting greater magnitude changes than lower latitudes, which differs from previous work suggesting regional synchronicity and coherence during the Holocene. Our analysis did not find evidence for coherent variability of primary producer community nor carbon export, highlighting the difficulty of identifying the complex interactions between environmental conditions, producers, and carbon export. Model‐proxy disagreement demonstrates the need for higher resolution model frameworks, but shows nonetheless that observed variability in the proxy records can be explained by a combination of greenhouse gas and orbital forcing. Here we suggest that the complex SST variations and marine ecosystem responses to forced changes are important factors that can drive disagreements in model projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗
DOE OSTI2022

NRAP Recommended Practices for Containment Assurance and Leakage Risk Quantification

This report presents recommendations for quantitatively assessing and managing leakage risk and effective containment assurance at geologic carbon storage (GCS) sites, amidst uncertainty in site characteristics and performance. These recommended practices are a product of the U.S. Department of Energy’s (DOE) Office of Fossil Energy and Carbon Management’s National Risk Assessment Partnership (NRAP)—a multi-year collaborative research effort of Lawrence Berkeley National Laboratory, Lawrence Livermore National Laboratory, Los Alamos National Laboratory, National Energy Technology Laboratory, and Pacific Northwest National Laboratory. These recommended practices are organized as a set of conceptual workflows that correspond with the stages of the GCS project lifecycle. They consider aspects of site characterization, quantification of potential leakage, assessing the spatial extent of potential environmental impacts, optimizing monitoring design for detection of potential leakage, evaluation of concordance between monitored and forecasted reservoir response, assessing conformance with expected/accepted performance, and determination of the timeline over which risk of impactful leakage can be expected to abate. These conceptual workflows are intended to aid stakeholders in making risk management decisions at GCS sites and to support development of sound, risk-based justifications for characterizing the subsurface, defining area of potential leakage impact, designing site monitoring plans, selecting leakage risk management/mitigation alternatives, and estimating site closure timelines.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2024

MemFriend: Understanding Memory Performance with Spatial-Temporal Affinity

In HPC applications, memory access behavior is one of the main factors affecting performance. Improving an application’s memory access behavior involves optimizing data layout and/or restructuring code, and requires studying spatial-temporal data locality. Existing data locality analyses focus on single-location metrics and are restricted to evaluating temporal locality. We introduce spatial-temporal affinity metrics that quantify temporal access proximity, forward access correlation, and nearby access correlation between pairs of memory locations. We describe methods for distinguishing between potential vs. realized affinity and for reasoning about affinity at multiple resolutions (3D, 2D, 1D). Finally, we construct spatial-temporal affinity signatures that classify memory behavior and that be used to reason about changes in software (data relayout, code refactoring) or hardware (caching, prefetching). We describe methods for signature visualization, interpretation, and quantitative comparison of signatures. We evaluate our methodology using applications with variants that contrast data structures, data layouts and algorithms. We show that spatial-temporal affinity analysis provides novel insights and enables predictive reasoning about application performance when contrasted with reuse distance analysis.

Suriyakumar, Yasodhadevi↗
DOE OSTI2024

Floating Wind Farm Layout Optimization Considering Moorings and Seabed Variations

This paper presents a method for optimizing the layout of floating wind farms that accounts for realistic seabed variations and the consequent adjustments to the mooring systems required for different turbine positions. The mooring lines of floating wind farms create large spatial constraints that are depth-dependent, since mooring designs must adapt to variations in seabed conditions over the array area. We develop a layout optimization methodology that addresses this, adjusting mooring system designs based on the local seabed characteristics as the layout changes and using steady-state models for the wake effects and mooring lines. The approach includes design algorithms that adjust the anchor positions and line length to achieve the desired mooring line profile for different water depths, and a layout optimization framework that implements spatial constraints between the turbines, mooring lines, and lease area boundaries. Demonstrating the method on several cases shows the effect of the seabed and spatial-constraint factors, as well as their interactions, on the optimal array layout. This demonstration paves the way for scaling up the method, using more powerful optimization algorithms to handle larger farm sizes and situations with more intensely varied seabed conditions.

17 WIND ENERGY↗