Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heterogeneous computing”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Dosimetric accuracy of a deterministic radiation transport based {sup 192}Ir brachytherapy treatment planning system. Part III. Comparison to Monte Carlo simulation in voxelized anatomical computational models

Purpose: To compare TG43-based and Acuros deterministic radiation transport-based calculations of the BrachyVision treatment planning system (TPS) with corresponding Monte Carlo (MC) simulation results in heterogeneous patient geometries, in order to validate Acuros and quantify the accuracy improvement it marks relative to TG43. Methods: Dosimetric comparisons in the form of isodose lines, percentage dose difference maps, and dose volume histogram results were performed for two voxelized mathematical models resembling an esophageal and a breast brachytherapy patient, as well as an actual breast brachytherapy patient model. The mathematical models were converted to digital imaging and communications in medicine (DICOM) image series for input to the TPS. The MCNP5 v.1.40 general-purpose simulation code input files for each model were prepared using information derived from the corresponding DICOM RT exports from the TPS. Results: Comparisons of MC and TG43 results in all models showed significant differences, as reported previously in the literature and expected from the inability of the TG43 based algorithm to account for heterogeneities and model specific scatter conditions. A close agreement was observed between MC and Acuros results in all models except for a limited number of points that lay in the penumbra of perfectly shaped structures in the esophageal model, or at distances very close to the catheters in all models. Conclusions: Acuros marks a significant dosimetry improvement relative to TG43. The assessment of the clinical significance of this accuracy improvement requires further work. Mathematical patient equivalent models and models prepared from actual patient CT series are useful complementary tools in the methodology outlined in this series of works for the benchmarking of any advanced dose calculation algorithm beyond TG43.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

A New LMP-Sensitivity-Based Heterogeneous Decomposition for Transmission and Distribution Coordinated Economic Dispatch

Transmission and distribution coordinated economic dispatch (TDCED) will be necessary to settle the problems caused by high penetrations of distributed generators in future grids. An efficient heterogeneous decomposition algorithm (HGD) was proposed to solve the TDCED, but it may diverge occasionally. To overcome this shortcoming, a new HGD (N-HGD) with enhanced convergency is proposed and studied in this paper. Compared with the previous HGD, the N-HGD similarly decomposes the TDCED into transmission ED and distribution ED subproblems, but it considers the sensitivity of the boundary locational marginal price to load injections to enhance convergency and avoid divergence. Two kinds of sensitivity computing methods, M-c and M-p, are studied. M-c is based on a centralized sensitivity equation, executed on the transmission side. M-p is based on a probing mechanism and can be executed distributedly on the distribution side. The reason for the N-HGD's improved convergency is theoretically analyzed, and the optimality is proven. Numerical tests show that N-HGD performs better than the original HGD, as well as a penalty-based modified HGD that needs a careful parametertuning process. Furthermore, compared with M-c, an M-p-based N-HGD is preferable because it has a much lower computational cost and no additional communicational cost.

distribution↗
DOE OSTI2018

Shock interactions with heterogeneous energetic materials

The complex physical phenomenon of shock wave interaction with material heterogeneities has significant importance and nevertheless remains little understood. In many materials, the observed macroscale response to shock loading is governed by characteristics of the microstructure. Yet the majority of computational studies aimed at predicting phenomena affected by these processes, such as initiation and propagation of detonation waves in explosives, or shock propagation in geological materials, employ continuum material and reactive burn model treatment. In an effort to highlight the grain-scale processes that underlie the observable effects in an energetic system, a grain-scale model for hexanitrostilbene (HNS) has been developed. Measured microstructures were used to produce synthetic computational representations of the pore structure, and a density functional theory molecular dynamics (DFT-MD) derived equation of state (EOS) was used for the fully dense HNS matrix. The explicit inclusion of microstructure along with a fully-dense EOS resulted in close agreement with historical shock compression experiments. More recent experiments on dynamic reaction threshold were also reproduced by inclusion of a global kinetics model. The complete model was shown to reproduce accurately the expected response of this heterogeneous material to shock loading. Mesoscale simulations were shown to provide clear insight into the nature of threshold behavior, and are a way to understand complex physical phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Convergence of micro-geochemistry and micro-geomechanics towards understanding proppant shale rock interaction: A Caney shale case study in southern Oklahoma, USA

As a direct outcome of economic development coupled with an increase in population, global energy demand will continue to rise in the coming decades. Although renewable energy sources are increasingly investigated for optimal production, the immediate needs require focus on energy sources that are currently available and reliable, with a minimal environmental impact; the efficient exploration and production of unconventional hydrocarbon resources is bridging the energy needs and energy aspirations, during the current energy transition period. The main challenges are related to the accurate quantification of the critical rock properties that influence production, their heterogeneity and the multiscale driven physico-chemical nature of rock–fluid interactions. A key feature of shale reservoirs is their low permeability due to dominating nanoporosity of the clay-rich matrix. As a means of producing these reservoirs in a cost-effective manner, a prerequisite is creation of hydraulic fracture networks capable of the highest level of continued conductivity. Fracturing fluid chemical design, formation brine geochemical composition, and rock mineralogy all contribute to swelling-induced conductivity damage. The Caney Shale is an organic-rich, often calcareous mudrock. Many studies have examined the impact that clay has on different kinds of shale productivity but there is currently no data reported on the Caney Shale in relation to horizontal drilling; all reported data on the Caney Shale is on vertical wells which are shallow, compared to an emerging play that is at double the depth. Here, in this work we develop geochemical–geomechanical integration of rock properties at micro-and nanoscales that can provide insights into the potential proppant embedment and its mitigation. The novel methodology amalgamates the following: computed X-ray tomography, scanning electron microscopy, energy dispersive spectroscopy, micro-indentation, and Raman spectroscopy techniques. Here, our results show that due to the multiscale heterogeneity in the Caney Shale, these geochemical and structural properties translate into a variation in mechanical properties that will impact interaction between the proppant and the host shale rock.

03 NATURAL GAS↗
DOE OSTI2022

Programming approaches for scalability, performance, and portability of combustion physics codes

Here, this paper presents the process, strategy, and results associated with porting a typical combustion physics flow solver to current state-of-the-art and future massively-parallel computer architectures. Major focus is placed on the distinct algorithmic structure of these types of codes and how it can be integrated with modern programming paradigms for heterogeneous platforms (i.e., distributed many-core systems with accelerators). An end-to-end case study is presented that exemplifies the process in a generic manner, which then serves as a clear guide with respect to the strategy and best practices leading to a robust and adaptable framework that performs well, is durable over time, is portable, and requires minimal human-effort. This end is accomplished beginning with the use of a mature, validated, structured, multiblock code framework optimized for application of both Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS). This code has been ported to a variety of platforms over the past decade, including most recently the Oak Ridge Leadership Computing Facility’s “Summit” Platform. The experience gained on these multiple platforms provides general insights and thus the results presented are not specific to any one code or platform other than the overarching trend toward distributed many-core systems with accelerators in order to move toward exascale performance. The resultant performance and scalability of the ported code is demonstrated on a real-world application; a state-of-the-art rotating detonation rocket engine simulation that matches the complex geometry and boundary conditions imposed as part of a companion experimental campaign.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Impurity precipitation in atomized particles evidenced by nano x-ray diffraction computed tomography

Performances and physical properties of high technology materials are influenced or even determined by their initial microstructure and by the behavior of impurity phases. Characterizing these impurities and their relations with the surrounding matrix is therefore of primary importance but it unfortunately often requires a destructive approach, with the risk of misinterpreting the observations. The improvement we have done in high resolution X-ray diffraction computed tomography combined with the use of an X-ray nanoprobe allows non-destructive crystallographic description of materials with microscopic heterogeneous microstructure (with a grain size between 10 nm and 10 μm). In this study, the grain localization in a 2D slice of a 20 μm solidified atomized γU-Mo particle is shown and a minority U(C,O) phase (1 wt. %) with sub-micrometer sized grains was characterized inside. Evidence is presented showing that the onset of U(C,O) grain crystallization can be described by a precipitation mechanism since one single U-Mo grain has direct orientation relationship with more than one surrounding U(C,O) grains.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Multiphase Modeling of Geologic Carbon Sequestration in Saline Aquifers

Abstract Geologic carbon sequestration ( GCS ) is being considered as a climate change mitigation option in many future energy scenarios. Mathematical modeling is routinely used to predict subsurface CO 2 and resident brine migration for the design of injection operations, to demonstrate the permanence of CO 2 storage, and to show that other subsurface resources will not be degraded. Many processes impact the migration of CO 2 and brine, including multiphase flow dynamics, geochemistry, and geomechanics, along with the spatial distribution of parameters such as porosity and permeability. In this article, we review a set of multiphase modeling approaches with different levels of conceptual complexity that have been used to model GCS . Model complexity ranges from coupled multiprocess models to simplified vertical equilibrium (VE) models and macroscopic invasion percolation models. The goal of this article is to give a framework of conceptual model complexity, and to show the types of modeling approaches that have been used to address specific GCS questions. Application of the modeling approaches is shown using five ongoing or proposed CO 2 injection sites. For the selected sites, the majority of GCS models follow a simplified multiphase approach, especially for questions related to injection and local‐scale heterogeneity. Coupled multiprocess models are only applied in one case where geomechanics have a strong impact on the flow. Owing to their computational efficiency, VE models tend to be applied at large scales. A macroscopic invasion percolation approach was used to predict the CO 2 migration at one site to examine details of CO 2 migration under the caprock.

Bandilla, Karl W.↗
DOE OSTI2017

HotSpots: Failure Cascades on Heterogeneous Critical Infrastructure Networks

Critical Infrastructure Systems such as transportation, water and power grid systems are vital to our national security, economy, and public safety. Recent events, like the 2012 hurricane Sandy, show how the interdependencies among different CI networks lead to catastrophic failures among the whole system. Hence, analyzing these CI networks, and modeling failure cascades on them becomes a very important problem. However, traditional models either do not take multiple CIs or the dynamics of the system into account, or model it simplistically. In this paper, we study this problem using a heterogeneous network viewpoint. We first construct heterogeneous CI networks with multiple components using national-level datasets. Then we study novel failure maximization problems on these networks, to compute critical nodes in such systems. We then provide HotSpots, a scalable and effective algorithm for these problems, based on careful transformations. Finally, we conduct extensive experiments on real CIS data from multiple US states, and show that our method HotSpots outperforms non-trivial baselines, gives meaningful results and that our approach gives immediate benefits in providing situational-awareness during large-scale failures.

Chen, Liangzhe↗
DOE OSTI2022

Characterization and Modeling of Deformation Induced Damage in Titanium Alloys

Over the lifetime of their DOE support, and specifically during the past funding cycle, the PIs have primarily focused on ambient temperature studies of the mechanisms of heterogeneous deformation and shear transfer at GBs. The PIs have developed and applied a broad array of experimental approaches complemented by advanced computational methods and simulations to enable new understanding of GB shear transfer, and a new approach for modeling this behavior at the mechanistic level has been established. The PIs have examined how slip behavior differs on the surface and interior of polycrystals using a variety of approaches. Some of the associated data has been used to generate 3-D models of experimentally-characterized microstructures that both support interpretation of deformation mechanisms as well as reveal where modeling needs to be improved to more closely agree with what is experimentally observed. The PIs progress in constitutive modeling of crystal plasticity has introduced explicit GB resistance to shear that depends on the mutual geometry of active slip systems. This past work places the PIs in a good position to take on the challenges associated with mechanistically understanding the interplay between GBS and slip interactions.

42 ENGINEERING↗
DOE OSTI2023

Opportunities for enhancing MLCommons efforts while leveraging insights from educational MLCommons earthquake benchmarks efforts

MLCommons is an effort to develop and improve the artificial intelligence (AI) ecosystem through benchmarks, public data sets, and research. It consists of members from start-ups, leading companies, academics, and non-profits from around the world. The goal is to make machine learning better for everyone. In order to increase participation by others, educational institutions provide valuable opportunities for engagement. In this article, we identify numerous insights obtained from different viewpoints as part of efforts to utilize high-performance computing (HPC) big data systems in existing education while developing and conducting science benchmarks for earthquake prediction. As this activity was conducted across multiple educational efforts, we project if and how it is possible to make such efforts available on a wider scale. This includes the integration of sophisticated benchmarks into courses and research activities at universities, exposing the students and researchers to topics that are otherwise typically not sufficiently covered in current course curricula as we witnessed from our practical experience across multiple organizations. As such, we have outlined the many lessons we learned throughout these efforts, culminating in the need for benchmark carpentry for scientists using advanced computational resources. The article also presents the analysis of an earthquake prediction code benchmark while focusing on the accuracy of the results and not only on the runtime; notedly, this benchmark was created as a result of our lessons learned. Energy traces were produced throughout these benchmarks, which are vital to analyzing the power expenditure within HPC environments. Additionally, one of the insights is that in the short time of the project with limited student availability, the activity was only possible by utilizing a benchmark runtime pipeline while developing and using software to generate jobs from the permutation of hyperparameters automatically. It integrates a templated job management framework for executing tasks and experiments based on hyperparameters while leveraging hybrid compute resources available at different institutions. The software is part of a collection called cloudmesh with its newly developed components, cloudmesh-ee (experiment executor) and cloudmesh-cc (compute coordinator).

58 GEOSCIENCES↗
DOE OSTI2024

Leveraging regional mesh refinement to simulate future climate projections for California using the Simplified Convection-Permitting E3SM Atmosphere Model Version 0

Abstract. The spatial heterogeneity related to complex topography in California demands high-resolution (< 5 km) modeling, but global convection-permitting climate models are computationally too expensive to run multi-decadal simulations. We developed a 3.25 km California climate modeling framework by leveraging regional mesh refinement (CARRM) using the U.S. Department of Energy (DOE)'s global Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) version 0. Four 5-year time periods (2015–2020, 2029–2034, 2044–2049, and 2094–2099) were simulated by nudging CARRM outside California to 1° coupled simulation of E3SMv1 under the Shared Socioeconomic Pathways (SSP)5-8.5 future scenario. The 3.25 km grid spacing adds considerable value to the prediction of the California climate changes, including more realistic high temperatures in the Central Valley and much improved spatial distributions of precipitation and snowpack in the Sierra Nevada and coastal stratocumulus. Under the SSP5-8.5 scenario, CARRM simulation predicts widespread warming of 6–10 °C over most of California, a 38 % increase in statewide average 30 d winter–spring precipitation, a near-complete loss of the alpine snowpack, and a sharp reduction in shortwave cloud radiative forcing associated with marine stratocumulus by the end of the 21st century. We note a climatological wet precipitation bias for the CARRM and discuss possible reasons. We conclude that SCREAM RRM is a technically feasible and scientifically valid tool for climate simulations in regions of interest, providing an excellent bridge to global convection-permitting simulations.

58 GEOSCIENCES↗
DOE OSTI2012

Validation of the U.S. NRC NGNP evaluation model with the HTTR

The High Temperature Test Reactor (HTTR) was modeled with TRITON/PARCS. Traditional light water reactor (LWR) homogenization methods rely on the short mean free paths of neutrons in LWR. In gas-cooled, graphite-moderated reactors like the HTTR neutrons have much longer mean free paths and penetrate further into neighboring assemblies than in LWRs. Because of this, conventional lattice calculations with a single assembly may not be valid. In addition to difficulties caused by the longer mean free paths, the HTTR presents unique axial and radial heterogeneities that require additional modifications to the single assembly homogenization method. To handle these challenges, the homogenization domain is decreased while the computational domain is increased. Instead of homogenizing a single hexagonal fuel assembly, the assembly is split into six triangles on the radial plane and five blocks axially in order to account for the placement of burnable poisons. Furthermore, the radial domain is increased beyond a single fuel assembly to account for spectrum effects from neighboring fuel, reflector, and control rod assemblies. A series of five two-dimensional cases, each closer to the full core, were calculated to evaluate the effectiveness of the homogenization method and cross-sections. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

From RNNs to Foundation Models: An Empirical Study on Commercial Building Energy Consumption

Accurate short-term energy consumption forecasting for commercial buildings is crucial for smart grid operations. While smart meters and deep learning models enable forecasting using past data from multiple buildings, data heterogeneity from diverse buildings can reduce model performance. The impact of increasing dataset heterogeneity in time series forecasting, while keeping size and model constant, is understudied. We tackle this issue using the ComStock dataset, which provides synthetic energy consumption data for U.S. commercial buildings. Two curated subsets, identical in size and region but differing in building type diversity, are used to assess the performance of various time series forecasting models, including finetuned open-source foundation models (FMs). The results show that dataset heterogeneity and model architecture have a greater impact on post-training forecasting performance than the parameter count. Moreover, despite the higher computational cost, finetuned FMs demonstrate competitive performance compared to base models trained from scratch.

commercial buildings↗
DOE OSTI2017

Electron Transfer Pathways Facilitating U(VI) Reduction by Fe(II) on Al- vs Fe-Oxides

This study continues mechanistic development of heterogeneous electron transfer (ET) pathways at mineral surfaces in aquatic environments that enable the reduction U(VI) by surface-associated Fe(II). Using computational molecular simulation within the framework of Marcus Theory, our findings highlight the importance of the configurations and interaction of the electron donor and acceptor species with the substrate, with respect to influencing its electronic structure and thereby the ability of semiconducting minerals to facilitate ET. U(VI) reduction by surface-associated Fe(II) (adsorbed or structurally incorporated into the lattice) on an insulating, corundum (001) surface (α-Al2O3) occurs when proximal inner-sphere (IS) surface complexes are formed, such that ET occurs through a combination of direct exchange (i.e., Fe d- and U f-orbitals overlap through space) and superexchange via intervening surface oxygen atoms. U(VI) reduction by coadsorbed Fe(II) on the isostructural semiconducting hematite (α-Fe2O3) basal surface requires either their direct electronic interaction (e.g., IS complexation) or mediation of this interaction indirectly through the surface via an intrasurface pathway. Conceptually possible longer-range ET by charge-hopping through surface Fe atoms was investigated to determine whether this indirect pathway is competitive with direct ET. The calculations show that energy barriers are large for this conduction-based pathway; interfacial ET into the hematite surface is endothermic (+80.1 kJ/mol) and comprises the rate-limiting step (10–6 s–1). The presence of the IS adsorbates appears to weaken the electronic coupling between underlying Fe ions within the surface, resulting in slower intra-surface ET (10–5 s–1) than expected in the bulk basal plane. Our findings lay out first insights into donor-acceptor communication via a charge-hopping pathway through the surface for heterogeneous reduction of U(VI) by Fe(II) and help provide a basis for experimental interrogation of this important process at mineral-water interfaces.

Taylor, S. D.↗
DOE OSTI2014

Inference of tumor evolution during chemotherapy by computational modeling and in situ analysis of genetic and phenotypic cellular diversity

Cancer therapy exerts a strong selection pressure that shapes tumor evolution, yet our knowledge of how tumors change during treatment is limited. Here, we report the analysis of cellular heterogeneity for genetic and phenotypic features and their spatial distribution in breast tumors pre- and post-neoadjuvant chemotherapy. We found that intratumor genetic diversity was tumor-subtype specific, and it did not change during treatment in tumors with partial or no response. However, lower pretreatment genetic diversity was significantly associated with pathologic complete response. In contrast, phenotypic diversity was different between pre- and post-treatment samples. We also observed significant changes in the spatial distribution of cells with distinct genetic and phenotypic features. We used these experimental data to develop a stochastic computational model to infer tumor growth patterns and evolutionary dynamics. Our results highlight the importance of integrated analysis of genotypes and phenotypes of single cells in intact tissues to predict tumor evolution.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

An asynchronous traversal engine for graph-based rich metadata management

Rich metadata in high-performance computing (HPC) systems contains extended information about users, jobs, data files, and their relationships. Property graphs are a promising data model to represent heterogeneous rich metadata flexibly. Specifically, a property graph can use vertices to represent different entities and edges to record the relationships between vertices with unique annotations. The high-volume HPC use case, with millions of entities and relationships, naturally requires an out-of-core distributed property graph database, which must support live updates (to ingest production information in real time), low-latency point queries (for frequent metadata operations such as permission checking), and large-scale traversals (for provenance data mining). Among these needs, large-scale property graph traversals are particularly challenging for distributed graph storage systems. Most existing graph systems implement a "level synchronous" breadth-first search algorithm that relies on global synchronization in each traversal step. This performs well in many problem domains; but a rich metadata management system is characterized by imbalanced graphs, long traversal lengths, and concurrent workloads, each of which has the potential to introduce or exacerbate stragglers (i.e., abnormally slow steps or servers in a graph traversal) that lead to low overall throughput for synchronous traversal algorithms. Previous research indicated that the straggler problem can be mitigated by using asynchronous traversal algorithms, and many graph-processing frameworks have successfully demonstrated this approach. Such systems require the graph to be loaded into a separate batch-processing framework instead of being iteratively accessed, however. In this work, we investigate a general asynchronous graph traversal engine that can operate atop a rich metadata graph in its native format. We outline a traversal-aware query language and key optimizations (traversal-affiliate caching and execution merging) necessary for efficient performance. We further explore the effect of different graph partitioning strategies on the traversal performance for both synchronous and asynchronous traversal engines. Our experiments show that the asynchronous graph traversal engine is more efficient than its synchronous counterpart in the case of HPC rich metadata processing, where more servers are involved and larger traversals are needed. Furthermore, the asynchronous traversal engine is more adaptive to different graph partitioning strategies.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2019

FFT, FMM, and multigrid on the road to exascale: Performance challenges and opportunities

FFT, FMM, and multigrid methods are widely used fast and highly scalable solvers for elliptic PDEs. However, emerging large-scale computing systems are introducing challenges in comparison to current petascale computers. Recent efforts (Dongarra et al. 2011) have identified several constraints in the design of exascale software that include massive concurrency, resilience management, exploiting the high performance of heterogeneous systems, energy efficiency, and utilizing the deeper and more complex memory hierarchy expected at exascale. In this paper, we perform a model-based comparison of the FFT, FMM, and multigrid methods in the context of these projected constraints. Finally, we use performance models to offer predictions about the expected performance on upcoming exascale system configurations based on current technology trends.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

SharP Data Constructs: Data Constructs to Enable Data-Centric Computing

Extreme-scale applications (i.e., Big-Compute) are becoming increasingly data-intensive, i.e., producing and consuming increasingly large amounts of data. The HPC systems traditionally used for these applications are now used for Big-Data applications such as data analytics, social network analysis, machine learning, and genomics. As a consequence of these trends, the system architecture should be flexible and data-centric. This can already be witnessed in the pre-exascale systems with TBs of on-node hierarchical and heterogeneous memories, PBs of system memory, low-latency, high-throughput networks, and many threaded cores. As such, the pre-exascale systems suit the needs of both Big-Compute and Big-Data applications. Though the system architecture is flexible enough to support both Big-Compute and Big-Data, we argue there is a software gap. Particularly, we need data-centric abstractions to leverage the full potential of the system, i.e., there is a need for native support for data resilience, the ability to express data locality and affinity, mechanisms to reduce data movement, the ability to share data, and abstractions to express User's data usage and data access patterns. In this paper, we (i) show the need for taking a holistic approach towards data-centric abstractions, (ii) show how these approaches were realized in the SHARed data-structure centric Programming abstraction (SharP) library, a data-structure centric programming abstraction, and (iii) apply these approaches to a variety of applications that demonstrate its usefulness. Particularly, we apply these approaches to QMCPack and the Graph500 benchmark and demonstrate the advantages of this approach on extreme-scale systems.

Aderholdt, William Ferrol↗