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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

New modelling capabilities in IDT

This work concerns the enhanced modelling capabilities of the discrete ordinates transport solver IDT. The novelties introduced allow for modelling unstructured geometries composed by a collection of X/Y segments and circles, and the use of reciprocity and conservation relations reduce the memory imprint as well as the computational cost of the method. IDT decomposes geometries in modular Cartesian patterns, which are the so-called Heterogeneous Cartesian Cells (HCCs), containing a chunk of the original unstructured geometries. Each HCC can be then discretized by superimposing a XY grid to refine locally the HCC. Unlike the most popular MOC, IDT performs the spatial sweeping by directional collision probabilities instead of trajectories. The sources and interface angular fluxes are expanded up to linear order. The accuracy of ray-tracing, the memory imprints together with the novel mesh refinement capabilities have been verified. A first set of preliminary results on PWR lattice problems will be presented in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗
DOE OSTI2017

Patient-specific modeling of individual sickle cell behavior under transient hypoxia

Sickle cell disease (SCD) is a highly complex genetic blood disorder in which red blood cells (RBC) exhibit heterogeneous morphology changes and decreased deformability. We employ a kinetic model for cell morphological sickling that invokes parameters derived from patient-specific data. This model is used to investigate the dynamics of individual sickle cells in a capillary-like microenvironment in order to address various mechanisms associated with SCD. We show that all RBCs, both hypoxia-unaffected and hypoxia-affected ones, regularly pass through microgates under oxygenated state. However, the hypoxia-affected cells undergo sickling which significantly alters cell dynamics. In particular, the dense and rigid sickle RBCs are obstructed thereby clogging blood flow while the less dense and deformable ones are capable of circumnavigating dead (trapped) cells ahead of them by choosing a serpentine path. Informed by recent experiments involving microfluidics that provide in vitro quantitative information on cell dynamics under transient hypoxia conditions, we have performed detailed computational simulations of alterations to cell behavior in response to morphological changes and membrane stiffening. Our model reveals that SCD exhibits substantial heterogeneity even within a particular density-fractionated subpopulation. These findings provide unique insights into how individual sickle cells move through capillaries under transient hypoxic conditions, and offer novel possibilities for designing effective therapeutic interventions for SCD.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

On proton CT reconstruction using MVCT-converted virtual proton projections

Purpose: To describe a novel methodology of converting megavoltage x-ray projections into virtual proton projections that are otherwise missing due to the proton range limit. These converted virtual proton projections can be used in the reconstruction of proton computed tomography (pCT). Methods: Relations exist between proton projections and multispectral megavoltage x-ray projections for human tissue. Based on these relations, these tissues can be categorized into: (a) adipose tissue; (b) nonadipose soft tissues; and (c) bone. These three tissue categories can be visibly identified on a regular megavoltage x-ray computed tomography (MVCT) image. With an MVCT image and its projection data available, the x-ray projections through heterogeneous anatomy can be converted to the corresponding proton projections using predetermined calibration curves for individual materials, aided by a coarse segmentation on the x-ray CT image. To show the feasibility of this approach, mathematical simulations were carried out. The converted proton projections, plotted on a proton sinogram, were compared to the simulated ground truth. Proton stopping power images were reconstructed using either the virtual proton projections only or a blend of physically available proton projections and virtual proton projections that make up for those missing due to the range limit. These images were compared to a reference image reconstructed from theoretically calculated proton projections. Results: The converted virtual projections had an uncertainty of {+-}0.8% compared to the calculated ground truth. Proton stopping power images reconstructed using a blend of converted virtual projections (48%) and physically available projections (52%) had an uncertainty of {+-}0.86% compared with that reconstructed from theoretically calculated projections. Reconstruction solely from converted virtual proton projections had an uncertainty of {+-}1.1% compared with that reconstructed from theoretical projections. If these images are used for treatment planning, the average proton range uncertainty is estimated to be less than 1.5% for an imaging dose in the milligray range. Conclusions: The proposed method can be used to convert x-ray projections into virtual proton projections. The converted proton projections can be blended with existing proton projections or can be used solely for pCT reconstruction, addressing the range limit problem of pCT using current therapeutic proton machines.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A Stochastic Intracellular Model of Anthrax Infection With Spore Germination Heterogeneity

We present a stochastic mathematical model of the intracellular infection dynamics of Bacillus anthracis in macrophages. Following inhalation of B. anthracis spores, these are ingested by alveolar phagocytes. Ingested spores then begin to germinate and divide intracellularly. This can lead to the eventual death of the host cell and the extracellular release of bacterial progeny. Some macrophages successfully eliminate the intracellular bacteria and will recover. Here, a stochastic birth-and-death process with catastrophe is proposed, which includes the mechanism of spore germination and maturation of B. anthracis . The resulting model is used to explore the potential for heterogeneity in the spore germination rate, with the consideration of two extreme cases for the rate distribution: continuous Gaussian and discrete Bernoulli. We make use of approximate Bayesian computation to calibrate our model using experimental measurements from in vitro infection of murine peritoneal macrophages with spores of the Sterne 34F2 strain of B. anthracis . The calibrated stochastic model allows us to compute the probability of rupture, mean time to rupture, and rupture size distribution, of a macrophage that has been infected with one spore. We also obtain the mean spore and bacterial loads over time for a population of cells, each assumed to be initially infected with a single spore. Our results support the existence of significant heterogeneity in the germination rate, with a subset of spores expected to germinate much later than the majority. Furthermore, in agreement with experimental evidence, our results suggest that most of the spores taken up by macrophages are likely to be eliminated by the host cell, but a few germinated spores may survive phagocytosis and lead to the death of the infected cell. Finally, we discuss how this stochastic modelling approach, together with dose-response data, allows us to quantify and predict individual infection risk following exposure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

An Efficient Algorithm for Mapping Imaging Data to 3D Unstructured Grids in Computational Biomechanics

Geometries for organ scale and multiscale simulations of organ function are now routinely derived from imaging data. However, medical images may also contain spatially heterogeneous information other than geometry that are relevant to such simulations either as initial conditions or in the form of model parameters. In this manuscript, we present an algorithm for the efficient and robust mapping of such data to imaging based unstructured polyhedral grids in parallel. We then illustrate the application of our mapping algorithm to three different mapping problems: 1) the mapping of MRI diffusion tensor data to an unstuctured ventricular grid; 2) the mapping of serial cyro-section histology data to an unstructured mouse brain grid; and 3) the mapping of CT-derived volumetric strain data to an unstructured multiscale lung grid. Execution times and parallel performance are reported for each case.

Einstein, Daniel R.↗
DOE OSTI2014

Harnessing the Heterogeneity of T Cell Differentiation Fate to Fine-Tune Generation of Effector and Memory T Cells

Recent studies show that naïve T cells bearing identical T cell receptors experience heterogeneous differentiation and clonal expansion processes. The factors controlling this outcome are not well characterized, and their contributions to immune cell dynamics are similarly poorly understood. In this study, we develop a computational model to elaborate mechanisms occurring within and between two important physiological compartments, lymph nodes and blood, to determine how immune cell dynamics are controlled. Our multiorgan (multi-compartment) model integrates cellular and tissue level events and allows us to examine the heterogeneous differentiation of individual precursor cognate naïveT cells to generate both effector and memory T lymphocytes. Using this model, we simulate a hypothetical immune response and reproduce both primary and recall responses to infection. Increased numbers of antigen-bearing dendritic cells (DCs) are predicted to raise production of both effector and memoryT cells, and distinct “sweet spots” of peptide-MHC levels on those DCs exist that favor CD4+ or CD8+ T cell differentiation toward either effector or memory cell phenotypes. This has important implications for vaccine development and immunotherapy.

cognate↗
DOE OSTI2019

Impacts of Spatial Heterogeneity and Temporal Non-Stationarity on Intensity-Duration-Frequency Estimates—A Case Study in a Mountainous California-Nevada Watershed

Changes in extreme precipitation events may require revisions of civil engineering standards to prevent water infrastructures from performing below the designated guidelines. Climate change may invalidate the intensity-duration-frequency (IDF) computation that is based on the assumption of data stationarity. Efforts in evaluating non-stationarity in the annual maxima series are inadequate, mostly due to the lack of long data records and convenient methods for detecting trends in the higher moments. In this study, using downscaled high resolution climate simulations of the historical and future periods under different carbon emission scenarios, we tested two solutions to obtain reliable IDFs under non-stationarity: (1) identify quasi-stationary time windows from the time series of interest to compute the IDF curves using data for the corresponding time windows; (2) introduce a parameter representing the trend in the means of the extreme value distributions. Focusing on a mountainous site, the Walker Watershed, the spatial heterogeneity and variability of IDFs or extremes are evaluated, particularly in terms of the terrain and elevation impacts. We compared observations-based IDFs that use the stationarity assumption with the two approaches that consider non-stationarity. The IDFs directly estimated based on the traditional stationarity assumption may underestimate the 100-year 24-h events by 10% to 60% towards the end of the century at most grids, resulting in significant under-designing of the engineering infrastructure at the study site. Strong spatial heterogeneity and variability in the IDF estimates suggest a preference for using high resolution simulation data for the reliable estimation of exceedance probability over data from sparsely distributed weather stations. Discrepancies among the three IDFs analyses due to non-stationarity are comparable to the spatial variability of the IDFs, underscoring a need to use an ensemble of non-stationary approaches to achieve unbiased and comprehensive IDF estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

SU-F-19A-01: APBI Brachytherapy Treatment Planning: The Impact of Heterogeneous Dose Calculations

Purpose: To quantify the difference in dose to PTV-EVAL and OARs (skin and rib) as calculated by (TG43) and heterogeneous calculations (CCC). Methods: 25 patient plans (5 Contura and 20 SAVI) were selected for analysis. Clinical dose distributions were computed with a commercially available treatment planning algorithm (TG43-D-(w,w)) and then recomputed with a pre-clinical collapsed cone convolution algorithm (CCCD-( m,m)). PTV-EVAL coverage (V90%, V95%), and rib and skin maximum dose were compared via percent difference. Differences in dose to normal tissue (V150cc, V200cc of PTV-EVAL) were also compared. Changes in coverage and maximum dose to organs at risk are reported in percent change, (100*(TG43 − CCC) / TG43)), and changes in maximum dose to normal tissue are absolute change in cc (TG43 − CCC). Results: Mean differences in V90, V95, V150, and V200 for the SAVI cases were −0.2%, −0.4%, −0.03cc, and −0.14cc, respectively, with maximum differences of −0.78%, −1.7%, 1.28cc, and 1.01cc, respectively. Mean differences in the 0.1cc dose to the rib and skin were −1.4% and −0.22%, respectively, with maximum differences of −4.5% and 16%, respectively. Mean differences in V90, V95, V150, and V200 for the Contura cases were −1.2%, −2.1%, −1.8cc, and −0.59cc, respectively, with maximum differences of −2.0%, −3.16%, −2.9cc, and −0.76cc, respectively. Mean differences in the 0.1cc dose to the rib and skin were −2.6% and −3.9%, respectively, with maximum differences of −3.2% and −5.7%, respectively. Conclusion: The effects of translating clinical knowledge based on D-(w,w) to plans reported in D-(m,m) are minimal (2% or less) on average, but vary based on the type and placement of the device, source, and heterogeneity information.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-406: Verification of Total Body Irradiation Commissioned MU Lookup Table Accuracy Using Treatment Planning System for Wide Range of Patient Sizes

Purpose: To verify the accuracy of total body irradiation (TBI) measurement commissioning data using the treatment planning system (TPS) for a wide range of patient separations. Methods: Our institution conducts TBI treatments with an 18MV photon beam at 380cm extended SSD using an AP/PA technique. Currently, the monitor units (MU) per field for patient treatments are determined using a lookup table generated from TMR measurements in a water phantom (75 × 41 × 30.5 cm3). The dose prescribed to an umbilicus midline point at spine level is determined based on patient separation, dose/ field and dose rate/MU. One-dimensional heterogeneous dose calculations from Pinnacle TPS were validated with thermoluminescent dosimeters (TLD) placed in an average adult anthropomorphic phantom and also in-vivo on four patients with large separations. Subsequently, twelve patients with various separations (17–47cm) were retrospectively analyzed. Computed tomography (CT) scans were acquired in the left and right decubitus positions from vertex to knee. A treatment plan for each patient was generated. The ratio of the lookup table MU to the heterogeneous TPS MU was compared. Results: TLD Measurements in the anthropomorphic phantom and large TBI patients agreed with Pinnacle calculated dose within 2.8% and 2%, respectively. The heterogeneous calculation compared to the lookup table agreed within 8.1% (ratio range: 1.014–1.081). A trend of reduced accuracy was observed when patient separation increases. Conclusion: The TPS dose calculation accuracy was confirmed by TLD measurements, showing that Pinnacle can model the extended SSD dose without commissioning a special beam model for the extended SSD geometry. The difference between the lookup table and TPS calculation potentially comes from lack of scatter during commissioning when compared to extreme patient sizes. The observed trend suggests the need for development of a correction factor between the lookup table and TPS dose calculations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-131: Dosimetric Impact and Evaluation of Different Heterogenity Algorithm in Volumetric Modulated Arc Therapy Plan for Stereotactic Ablative Radiotherapy Lung Treatment with the Flattening Filter Free Beam

Purpose: The present study aimed to investigate the dosimetric impacts of the anisotropic analytic algorithm (AAA) and the Acuros XB (AXB) plan for lung stereotactic ablative radiation therapy using flattening filter-free (FFF) beam. We retrospectively analyzed 10 patients. Methods: We retrospectively analyzed 10 patients. The dosimetric parameters for the target and organs at risk (OARs) from the treatment plans calculated with these dose calculation algorithms were compared. The technical parameters, such as the computation times and the total monitor units (MUs), were also evaluated. Results: A comparison of DVHs from AXB and AAA showed that the AXB plan produced a high maximum PTV dose by average 4.40% with a statistical significance but slightly lower mean PTV dose by average 5.20% compared to the AAA plans. The maximum dose to the lung was slightly higher in the AXB compared to the AAA. For both algorithms, the values of V5, V10 and V20 for ipsilateral lung were higher in the AXB plan more than those of AAA. However, these parameters for contralateral lung were comparable. The differences of maximum dose for the spinal cord and heart were also small. The computation time of AXB was found fast with the relative difference of 13.7% than those of AAA. The average of monitor units (MUs) for all patients was higher in AXB plans than in the AAA plans. These results indicated that the difference between AXB and AAA are large in heterogeneous region with low density. Conclusion: The AXB provided the advantages such as the accuracy of calculations and the reduction of the computation time in lung stereotactic ablative radiotherapy (SABR) with using FFF beam, especially for VMAT planning. In dose calculation with the media of different density, therefore, the careful attention should be taken regarding the impacts of different heterogeneity correction algorithms. The authors report no conflicts of interest.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

The Heterogeneous Integration of Electronic Components

Heterogeneous integration (HI) of electronics components is broadly recognized as a powerful and crucial enabler for the continued growth of computing and communication. From 2010 onwards, the value of HI is increasingly visible in the advanced packaging used in artificial intelligence, high-performance computing, smartphones and communications product implementations. In this Perspective, we argue that HI is crucial to semiconductors and more broadly to the continued evolution of computing and communications. We use leading-edge advanced packaging examples to represent the value, advancements and opportunities for HI. To succeed, it is critical to develop comprehensive HI roadmaps that inform collaborations across the design, manufacturing and reliability spectrum between systems architects, packaging and semiconductor technologists to common goals. Although this article does not provide a full roadmap, we instead detail additional parameters for artificial intelligence, smartphone and other cellular communication devices, and their constituent building blocks including interconnects, power electronics, photonics, thermal management, reliability, modelling and co-design, to foster greater collaboration opportunities among academia, research laboratories and industry.

42 ENGINEERING↗
DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗
DOE OSTI2019

Heterogeneous damage in Li-ion batteries: Experimental analysis and theoretical modeling

We assess the heterogeneous electrochemistry and mechanics in a composite electrode of commercial batteries using synchrotron X-ray tomography analysis and microstructure-resolved computational modeling. We visualize the morphological defects at multi-scales ranging from the macroscopic composite, particle ensembles, to individual single particles. Particle fracture and interfacial debonding are identified in a large set of tomographic data of active particles. Mechanical failure in the regime near the separator is more severe than toward the current collector. The active particles close to the separator experience deeper charge and discharge over cycles and thus are more mechanically loaded. The difference in the Li activity originates from the polarization of the electrolyte potential and the non-uniform distribution of the activation energy for the charge transfer reaction. We model the kinetics of intergranular fracture and interfacial degradation to confirm that the various Li activities are the major cause of the heterogeneous damage. The interfacial failure may reconstruct the conductive network and redistribute the electrochemical activities that render a dynamic nature of electrochemistry and mechanics evolving over time in the composite electrode. We further quantify the influence of the mechanical damage on the metrics of battery performance. We simulate the electrochemical impedance profile to build a relationship between the interfacial debonding and the impedance of electron transport and surface charge transfer. Lastly, the mechanical failure disrupts the conduction path of electrons and results in significant polarization and capacity loss in batteries.

25 ENERGY STORAGE↗
DOE OSTI2012

Three-Dimensional Bayesian Geostatistical Aquifer Characterization at the Hanford 300 Area using Tracer Test Data

Tracer testing under natural or forced gradient flow holds the potential to provide useful information for characterizing subsurface properties, through monitoring, modeling and interpretation of the tracer plume migration in an aquifer. Non-reactive tracer experiments were conducted at the Hanford 300 Area, along with constant-rate injection tests and electromagnetic borehole flowmeter (EBF) profiling. A Bayesian data assimilation technique, the method of anchored distributions (MAD) [Rubin et al., 2010], was applied to assimilate the experimental tracer test data with the other types of data and to infer the three-dimensional heterogeneous structure of the hydraulic conductivity in the saturated zone of the Hanford formation. In this study, the Bayesian prior information on the underlying random hydraulic conductivity field was obtained from previous field characterization efforts using the constant-rate injection tests and the EBF data. The posterior distribution of the conductivity field was obtained by further conditioning the field on the temporal moments of tracer breakthrough curves at various observation wells. MAD was implemented with the massively-parallel three-dimensional flow and transport code PFLOTRAN to cope with the highly transient flow boundary conditions at the site and to meet the computational demands of MAD. A synthetic study proved that the proposed method could effectively invert tracer test data to capture the essential spatial heterogeneity of the three-dimensional hydraulic conductivity field. Application of MAD to actual field data shows that the hydrogeological model, when conditioned on the tracer test data, can reproduce the tracer transport behavior better than the field characterized without the tracer test data. This study successfully demonstrates that MAD can sequentially assimilate multi-scale multi-type field data through a consistent Bayesian framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Probabilistic Forecast of Single–Phase Flow in Porous Media With Uncertain Properties

Uncertainty about geologic makeup and properties of the subsurface renders inadequate a unique quantitative prediction of flow and transport. Instead, multiple alternative scenarios have to be explored within the probabilistic framework, typically by means of Monte Carlo simulations (MCS). These can be computationally expensive, and often prohibitively so, especially when the goal is to compute the tails of a distribution, i.e., probabilities of rare events, which are necessary for risk assessment and decision making under uncertainty. We deploy the method of distributions to derive a deterministic equation for the cumulative distribution function (CDF) of hydraulic head in an aquifer with uncertain (random) hydraulic conductivity. The CDF equation relies on a self-consistent closure approximation, which ensures that the resulting CDF of hydraulic head has the same mean and variance as those computed with either statistical moment equation (the approach used herein) or MCS. We conduct a series of numerical experiments dealing with steady-state two-dimensional flow driven by either a natural hydraulic head gradient or a pumping well. Furthermore, these experiments reveal that the CDF method remains accurate and robust for highly heterogeneous formations with the variance of log conductivity as large as five. For the same accuracy, it is also up to four orders of magnitude faster than MCS in computing hydraulic head with a required degree of confidence (probability).

42 ENGINEERING↗
DOE OSTI2020

Scalable Heterogeneous Execution of a Coupled-Cluster Model with Perturbative Triples

The CCSD(T) coupled-cluster model with perturbative triples is considered a gold standard for computational modeling of the correlated behavior of electrons in molecular systems. A fundamental constraint is the relatively small global-memory capacity in GPUs compared to the main-memory capacity on host nodes, necessitating relatively smaller tile sizes for high-dimensional tensor contractions in NWChem's GPU-accelerated implementation of the CCSD(T) method. A coordinated redesign is described to address this limitation and associated data movement overheads, including a novel fused GPU kernel for a set of tensor contractions, along with inter-node communication optimization and data caching. The new implementation of GPU-accelerated CCSD(T) improves overall performance by 3.4x. Finally, we discuss the trade-offs in using this fused algorithm on current and future supercomputing platforms.

Kim, Jinsung↗