Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster 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,531 records · Page 85

Advances in ArborX to support exascale applications

ArborX is a performance portable geometric search library developed as part of the Exascale Computing Project (ECP). In this paper, we explore a collaboration between ArborX and a cosmological simulation code HACC. Large cosmological simulations on exascale platforms encounter a bottleneck due to the in-situ analysis requirements of halo finding, a problem of identifying dense clusters of dark matter (halos). This problem is solved by using a density-based DBSCAN clustering algorithm. With each MPI rank handling hundreds of millions of particles, it is imperative for the DBSCAN implementation to be efficient. In addition, the requirement to support exascale supercomputers from different vendors necessitates performance portability of the algorithm. We describe how this challenge problem guided ArborX development, and enhanced the performance and the scope of the library. We explore the improvements in the basic algorithms for the underlying search index to improve the performance, and describe several implementations of DBSCAN in ArborX. Further, we report the history of the changes in ArborX and their effect on the time to solve a representative benchmark problem, as well as demonstrate the real world impact on production end-to-end cosmology simulations.

97 MATHEMATICS AND COMPUTING↗

The apparent reversal of the Law of Mass Action in concentrated multicomponent aqueous solutions

Although it is well known that aqueous electrolyte solutions behave non-ideally, few are so non-ideal that adding an electrolyte to a saturated solution with a common ion enhances rather than depresses solubility. Nonetheless, this apparent reversal of the Law of Mass Action (LMA) has been observed in simple nitrate solutions at high concentrations. A leading hypothesis is that ion clusters are formed and stabilized by ions having different charge densities. We report the present study examines this concept, in part by reviewing relevant data from multicomponent aqueous solutions containing sodium nitrate (NaNO 3 ), sodium nitrite (NaNO 2 ), sodium hydroxide (NaOH), and sodium aluminate (NaAl(OH) 4 ) - the major constituents in alkaline nuclear waste. Here, NaOH and NaAl(OH) 4 did not enhance the solubility of NaNO 3 or NaNO 2 , whereas NaNO 2 and NaNO 3 enhanced rather than depressed the solubility of each other despite each having the sodium cation (Na + ) in common. Solutions evaluated in this study have more than 20 molal total Na + concentration, and most have less than one mole of water per mole of ion. Thus, this reversal of the LMA occurs in solutions where there is not enough water to fully hydrate the ions, pointing to the importance of ion cluster formation. Within the composition range of the solutions analyzed here, this reversal of the LMA occurs regardless of NaAl(OH) 4 and NaOH concentrations. Elevated temperatures also result in the reversal of the LMA in the subsystems NaNO 2 -NaNO 3 -H 2 O and NaOH-NaNO 3 -H 2 O, consistent with spectroscopic and computational studies showing enhanced interactions at higher temperatures. Although speciation in these highly concentrated electrolyte solutions is not well understood, the stabilization of ion clusters by: (i) optimizing charge density around the ions through different combinations of mixed cations and anions in solution; and (ii) high temperatures where ions experience prolonged contact, offers important clues for future research into how these species control solubility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ORCA: Outlier detection and Robust Clustering for Attributed graphs

Here, a framework is proposed to simultaneously cluster objects and detect anomalies in attributed graph data. Our objective function along with the carefully constructed constraints promotes interpretability of both the clustering and anomaly detection components, as well as scalability of our method. In addition, we developed an algorithm called Outlier detection and Robust Clustering for Attributed graphs (ORCA) within this framework. ORCA is fast and convergent under mild conditions, produces high quality clustering results, and discovers anomalies that can be mapped back naturally to the features of the input data. The efficacy and efficiency of ORCA is demonstrated on real world datasets against multiple state-of-the-art techniques.

97 MATHEMATICS AND COMPUTING↗

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗

Simulating X-ray Absorption Spectra with Linear-Response Density Cumulant Theory

We present a new approach for simulating X-ray absorption spectra based on linear-response density cumulant theory (LR-DCT) [Copan, A. V.; Sokolov, A. Yu. J. Chem. Theory Comput. 2018, 14, 4097–4108]. Our new method combines the LR-ODC-12 formulation of LR-DCT with core- valence separation approximation (CVS) that allows us to efficiently access high-energy core-excited states. We describe our computer implementation of the CVS-approximated LRODC-12 method (CVS-ODC-12) and benchmark its performance by comparing simulated X-ray absorption spectra to those obtained from experiment for several small molecules. Our results demonstrate that the CVS-ODC-12 method shows good agreement with experiment for relative spacings between transitions and their intensities, but the excitation energies are systematically overestimated. When compared to results from excited-state coupled cluster methods with single and double excitations, the CVS-ODC-12 method shows a similar performance for intensities and peak separations, while coupled cluster spectra are less shifted, relative to experiment. An important advantage of CVS-ODC-12 is that its excitation energies are computed by diagonalizing a Hermitian matrix, which enables efficient computation of transition intensities.

Peng, Ruojing↗

CaMn 3 IV O 4 Cubane Models of the Oxygen‐Evolving Complex: Spin Ground States S <9/2 and the Effect of Oxo Protonation

Abstract We report the single crystal XRD and MicroED structure, magnetic susceptibility, and EPR data of a series of CaMn 3 IV O 4 and YMn 3 IV O 4 complexes as structural and spectroscopic models of the cuboidal subunit of the oxygen‐evolving complex (OEC). The effect of changes in heterometal identity, cluster geometry, and bridging oxo protonation on the spin‐state structure was investigated. In contrast to previous computational models, we show that the spin ground state of CaMn 3 IV O 4 complexes and variants with protonated oxo moieties need not be S =9/2. Desymmetrization of the pseudo ‐ C 3 ‐symmetric Ca(Y)Mn 3 IV O 4 core leads to a lower S =5/2 spin ground state. The magnitude of the magnetic exchange coupling is attenuated upon oxo protonation, and an S =3/2 spin ground state is observed in CaMn 3 IV O 3 (OH). Our studies complement the observation that the interconversion between the low‐spin and high‐spin forms of the S 2 state is pH‐dependent, suggesting that the (de)protonation of bridging or terminal oxygen atoms in the OEC may be connected to spin‐state changes.

Lee, Heui Beom↗

CaMn 3 IV O 4 Cubane Models of the Oxygen‐Evolving Complex: Spin Ground States S <9/2 and the Effect of Oxo Protonation

Abstract We report the single crystal XRD and MicroED structure, magnetic susceptibility, and EPR data of a series of CaMn 3 IV O 4 and YMn 3 IV O 4 complexes as structural and spectroscopic models of the cuboidal subunit of the oxygen‐evolving complex (OEC). The effect of changes in heterometal identity, cluster geometry, and bridging oxo protonation on the spin‐state structure was investigated. In contrast to previous computational models, we show that the spin ground state of CaMn 3 IV O 4 complexes and variants with protonated oxo moieties need not be S =9/2. Desymmetrization of the pseudo ‐ C 3 ‐symmetric Ca(Y)Mn 3 IV O 4 core leads to a lower S =5/2 spin ground state. The magnitude of the magnetic exchange coupling is attenuated upon oxo protonation, and an S =3/2 spin ground state is observed in CaMn 3 IV O 3 (OH). Our studies complement the observation that the interconversion between the low‐spin and high‐spin forms of the S 2 state is pH‐dependent, suggesting that the (de)protonation of bridging or terminal oxygen atoms in the OEC may be connected to spin‐state changes.

Lee, Heui Beom↗

Role of Interactions and Correlations on Collective Dynamics of Molecular Motors Along Parallel Filaments

Cytoskeletal motors known as motor proteins are molecules that drive cellular transport along several parallel cytoskeletal filaments and support many biological processes. Experimental evidences suggest that they interact with the nearest molecules of their filament while performing any mechanical work. These interactions modify the microscopic level properties of motor proteins. In this work, a new version of two-channel totally asymmetric simple exclusion process, that incorporates the intra-channel interactions in a thermodynamically consistent way, is proposed. As the existing approaches for multi-channel systems deviate from analyzing the combined effect of inter and intra-channel interactions, a new approach known as modified vertical cluster mean field is developed. The approach along with Monte Carlo simulations successfully encounters some correlations and computes the complex dynamic properties of the system. Role of symmetry of interactions and inter-channel coupling is observed on the phase diagrams, maximal particle current and its corresponding optimal interaction strength. Surprisingly, for all values of coupling rate and most of the interaction splittings, the optimal interaction strength corresponding to maximal current belongs to the case of weak repulsive interactions. Moreover, for weak interaction splittings and with an increase in the coupling rate, the optimal interaction strength tends towards the known experimental results. The effect of coupling as well as interaction energy is also measured for correlations. They are found to be short-range and weaker for repulsive and weak attractive interactions while they are long-range and stronger for large attractions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward Discord : Code for Simulating Continuous Spin Systems

A new computational tool to simulate classical spin systems with frustrated crystal structures is presented. Complementary single- and cluster-spin flip algorithms are implemented to calculate the diffuse scattering patterns, spin-pair correlations, and thermodynamic quantities. Test cases of geometrically frustrated kagome, pyrochlore, and cubic systems are detailed. Two recent scientific cases are also shown here. This new method, together with recent developments of the rmc-discord package (https://github.com/zjmorgan/rmc-discord), represent integrated and strategic step in a complete forward and reverse Monte Carlo framework discord.

36 MATERIALS SCIENCE↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Enthalpies of formation for Criegee intermediates: A correlation energy convergence study

Criegee intermediates, formed from the ozonolysis of alkenes, are known to have a role in atmospheric chemistry, including the modulation of the oxidizing capacity of the troposphere. Although studies have been conducted since their discovery, the synthesis of these species in the laboratory has ushered in a new wave of investigations of these structures, both theoretically and experimentally. In some of these theoretical studies, high-order corrections for correlation energy are included to account for the mid multi-reference character found in these systems. Many of these studies include a focus on kinetics; therefore, the calculated energies should be accurate (<1 kcal/mol in error). In this research, we compute the enthalpies of formation for a small set of Criegee intermediates, including higher-order coupled cluster corrections for correlation energy up to coupled cluster with perturbative quintuple excitations. The enthalpies of formation for formaldehyde oxide, anti-acetaldehyde oxide, syn-acetaldehyde oxide, and acetone oxide are presented at 0 K as 26.5, 15.6, 12.2, and 0.1 kcal mol−1, respectively. Additionally, we do not recommend the coupled cluster with perturbative quadruple excitations [CCSDT(Q)] energy correction, as it is approximately twice as large as that of the coupled cluster with full quadruple excitations (CCSDTQ). Half of the CCSDT(Q) energy correction may be included as a reliable, cost-effective estimation of CCSDTQ energies for Criegee intermediates.

Chemistry↗

Helium segregation on surfaces of plasma-exposed tungsten

Here we report a hierarchical multi-scale modeling study of implanted helium segregation on surfaces of tungsten, considered as a plasma facing component in nuclear fusion reactors. We employ a hierarchy of atomic-scale simulations based on a reliable interatomic interaction potential, including molecular-statics simulations to understand the origin of helium surface segregation, targeted molecular-dynamics (MD) simulations of near-surface cluster reactions, and large-scale MD simulations of implanted helium evolution in plasma-exposed tungsten. We find that small, mobile He-n (1 <= n <= 7) clusters in the near-surface region are attracted to the surface due to an elastic interaction force that provides the thermodynamic driving force for surface segregation. Elastic interaction force induces drift fluxes of these mobile Hen clusters, which increase substantially as the migrating clusters approach the surface, facilitating helium segregation on the surface. Moreover, the clusters' drift toward the surface enables cluster reactions, most importantly trap mutation, in the near-surface region at rates much higher than in the bulk material. Moreover, these near-surface cluster dynamics have significant effects on the surface morphology, near-surface defect structures, and the amount of helium retained in the material upon plasma exposure. We integrate the findings of such atomic-scale simulations into a properly parameterized and validated spatially dependent, continuum-scale reaction-diffusion cluster dynamics model, capable of predicting implanted helium evolution, surface segregation, and its near-surface effects in tungsten. This cluster-dynamics model sets the stage for development of fully atomistically informed coarse-grained models for computationally efficient simulation predictions of helium surface segregation, as well as helium retention and surface morphological evolution, toward optimal design of plasma facing components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of a pair-based, joint-likelihood association approach for regional infrasound event identification

SUMMARY A Bayesian framework for the association of infrasonic detections is presented and evaluated for analysis at regional propagation scales. A pair-based, joint-likelihood association approach is developed that identifies events by computing the probability that individual detection pairs are attributable to a hypothetical common source and applying hierarchical clustering to identify events from the pair-based analysis. The framework is based on a Bayesian formulation introduced for infrasonic source localization and utilizes the propagation models developed for that application with modifications to improve the numerical efficiency of the analysis. Clustering analysis is completed using hierarchical analysis via weighted linkage for a non-Euclidean distance matrix defined by the negative log-joint-likelihood values. The method is evaluated using regional synthetic data with propagation distances of hundreds of kilometres in order to study the sensitivity of the method to uncertainties and errors in backazimuth and time of arrival. The method is found to be robust and stable for typical uncertainties, able to effectively distinguish noise detections within the data set from those in events, and can be made numerically efficient due to its ease of parallelization.

58 GEOSCIENCES↗

Resilience-Oriented DG Siting and Sizing Considering Stochastic Scenario Reduction

In this paper, a fuel-based distributed generator (DG) allocation strategy is proposed to enhance the distribution system resilience against extreme weather. The long-term planning problem is formulated as a two-stage stochastic mixed-integer programming (SMIP). The first stage is to make decisions of DG siting and sizing under the given budget constraint. In the second stage, a post-extreme-event-restoration (PEER) is employed to minimize the operating cost in an uncertain fault scenario. In particular, this study proposes a method to select the most representative scenarios for the SMIP. First, a Monte Carlo Simulation (MCS) is introduced to generate sufficient scenarios considering random fault locations and load profiles. Then, the number of scenarios is reduced by the K-means clustering algorithm. The advantage of scenario reduction is to make a trade-off between accuracy and computational efficiency. Finally, the SMIP is solved by the progressive hedging algorithm. Here, the case studies of the IEEE 33-bus and 123-bus test systems demonstrate the effectiveness of the proposed algorithm in reducing the expected energy not served (EENS), which is a critical criterion of resilience.

42 ENGINEERING↗

Scalar Field Comparison with Topological Descriptors: Properties and Applications for Scientific Visualization

In topological data analysis and visualization, topological descriptors such as persistence diagrams, merge trees, contour trees, Reeb graphs, and Morse–Smale complexes play an essential role in capturing the shape of scalar field data. Herein we present a state–of–the–art report on scalar field comparison using topological descriptors. We provide a taxonomy of existing approaches based on visualization tasks associated with three categories of data: single fields, time–varying fields, and ensembles. These tasks include symmetry detection, periodicity detection, key event/feature detection, feature tracking, clustering, and structure statistics. Our main contributions include the formulation of a set of desirable mathematical and computational properties of comparative measures, and the classification of visualization tasks and applications that are enabled by these measures.

97 MATHEMATICS AND COMPUTING↗

SU-E-J-107: Supervised Learning Model of Aligned Collagen for Human Breast Carcinoma Prognosis

Purpose: Our goal is to develop and apply a set of optical and computational tools to enable large-scale investigations of the interaction between collagen and tumor cells. Methods: We have built a novel imaging system for automating the capture of whole-slide second harmonic generation (SHG) images of collagen in registry with bright field (BF) images of hematoxylin and eosin stained tissue. To analyze our images, we have integrated a suite of supervised learning tools that semi-automatically model and score collagen interactions with tumor cells via a variety of metrics, a method we call Electronic Tumor Associated Collagen Signatures (eTACS). This group of tools first segments regions of epithelial cells and collagen fibers from BF and SHG images respectively. We then associate fibers with groups of epithelial cells and finally compute features based on the angle of interaction and density of the collagen surrounding the epithelial cell clusters. These features are then processed with a support vector machine to separate cancer patients into high and low risk groups. Results: We validated our model by showing that eTACS produces classifications that have statistically significant correlation with manual classifications. In addition, our system generated classification scores that accurately predicted breast cancer patient survival in a cohort of 196 patients. Feature rank analysis revealed that TACS positive fibers are more well aligned with each other, generally lower density, and terminate within or near groups of epithelial cells. Conclusion: We are working to apply our model to predict survival in larger cohorts of breast cancer patients with a diversity of breast cancer types, predict response to treatments such as COX2 inhibitors, and to study collagen architecture changes in other cancer types. In the future, our system may be used to provide metastatic potential information to cancer patients to augment existing clinical assays.

60 APPLIED LIFE SCIENCES↗