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 955 records · Page 53

The electron affinity of the uranium atom

The results of a combined experimental and computational study of the uranium atom are presented with the aim of determining its electron affinity. Experimentally, the electron affinity of uranium was measured via negative ion photoelectron spectroscopy of the uranium atomic anion, U – . Computationally, the electron affinities of both thorium and uranium were calculated by conducting relativistic coupled-cluster and multi-reference configuration interaction calculations. The experimentally determined value of the electron affinity of the uranium atom was determined to be 0.309 ± 0.025 eV. The computationally predicted electron affinity of uranium based on composite coupled cluster calculations and full four-component spin–orbit coupling was found to be 0.232 eV. Predominately due to a better convergence of the coupled cluster sequence for Th and Th – , the final calculated electron affinity of Th, 0.565 eV, was in much better agreement with the accurate experimental value of 0.608 eV. In both scenarios, the ground state of the anion corresponds to electron attachment to the 6d orbital.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Using Hybrid Model OpenSHMEM + CUDA to Implement the SHOC Benchmark Suite

This work describes the process of porting the Scalable HeterOgeneous Computing (SHOC) benchmark suite from the hybrid MPI + CUDA implementation to OpenSHMEM + CUDA. SHOC includes a wide variety of benchmark kernels used to measure accelerator performance in both single node and cluster configurations. The hybrid model implementation attempts to place all major computation on accelerator devices, and uses MPI to synchronize and aggregate results. In some cases, MPI Groups are used to gradually reduce the number of accelerators used for computation as the problem size drops. Porting this behavior to OpenSHMEM required implementing several synchronizing collective operations, and using SHMEM teams to replace MPI Group functionality. Benchmark results on a Cray XK7 system with one GPU per compute node show that SHMEM performance is equal to MPI performance in these hybrid tasks. These results and porting experience show that using OpenSHMEM for accelerator devices benefits from adding functionality for synchronization and teams, and would further benefit from adding support for communication within accelerator kernels.

97 MATHEMATICS AND COMPUTING↗

Diffusion Monte Carlo studies of MB-pol (H{sub 2}O){sub 2−6} and (D{sub 2}O){sub 2−6} clusters: Structures and binding energies

We employ the diffusion Monte Carlo (DMC) method in conjunction with the recently developed, ab initio-based MB-pol potential energy surface to characterize the ground states of small (H{sub 2}O){sub 2−6} clusters and their deuterated isotopomers. Observables, other than the ground state energies, are computed using the descendant weighting approach. Among those are various spatial correlation functions and relative isomer fractions. Interestingly, the ground states of all clusters considered in this study, except for the dimer, are delocalized over at least two conformations that differ by the orientation of one or more water monomers with the relative isomer populations being sensitive to the isotope substitution. Most remarkably, the ground state of the (H{sub 2}O){sub 6} hexamer is represented by four distinct cage structures, while that of (D{sub 2}O){sub 6} is dominated by the prism, i.e., the global minimum geometry, with a very small contribution from a prism-book geometry. In addition, for (H{sub 2}O){sub 6} and (D{sub 2}O){sub 6}, we performed DMC calculations to compute the ground states constrained to the cage and prism geometries. These calculations compared results for three different potentials, MB-pol, TTM3/F, and q-TIP4P/F.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On a simple way to calculate electronic resonances for polyatomic molecules

We propose a simple method for calculation of low-lying shape electronic resonances of polyatomic molecules. The method introduces a perturbation potential and requires only routine bound-state type calculations in the real domain of energies. Such a calculation is accessible by most of the free or commercial quantum chemistry software. The presented method is based on the analytical continuation in a coupling constant model, but unlike its previous variants, we experience a very stable and robust behavior for higher-order extrapolation functions. Moreover, the present approach is independent of the correlation treatment used in quantum many-electron computations and therefore we are able to apply Coupled Clusters (CCSD-T) level of the correlation model. We demonstrate these properties on determination of the resonance position and width of the {sup 2}Π{sub u} temporary negative ion state of diacetylene using CCSD-T level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Alkyltin clusters: the less symmetric Keggin isomers

The Keggin structure is prevalent in nature and synthesis, self-assembled from many metals across the periodic table as both isolated clusters and building blocks of condensed framework oxides. Here we present a one-step synthesis to obtain the sodium-centered butyltin Keggin ion in high yield and high purity, important for mechanistic nanolithography studies. Extensive solution characterization (small-angle X-ray scattering, 1 H, 13 C and 119 Sn nuclear magnetic resonance, electrospray mass spectrometry) also confirms solutions contain only the Na-centered dodecamers. We report three butyltin Keggin structures: the β-isomer (β-NaSn 12 ), the γ-isomer (γ-NaSn 12 ), and a γ-isomer capped with an additional butyltin (γ-NaSn 13 ). All Keggin ions presented here have the general formula [NaO 4 BuSn 12 (OCH 3 ) 12 (O,OH) 12 ] (Bu = butyl), and are of neutral charge. The lack of counterions (OH – ) facilitates mechanistic lithographic studies without inference from hydrolysis chemistry. The methanol reaction media enables solubility and ligates the cluster, both important to obtain high purity materials. Despite the monospecific nature of the NaSn 12 solutions, NMR reveals both isomer interconversion and ligand exchange. DFT computational comparisons of our three isolated structures, the capped β-isomer (β-NaSn 13 ), along with hypothetical α-isomers (α-NaSn 12 and α-NaSn 13 ), showed that the stability ranks β-NaSn 12 > γ-NaSn 12 > α-NaSn 12 , consistent with experimental observation. The uncapped isomers were computationally determined to be more stable than the respective capped analogues. These clusters provide a unique opportunity to investigate the lower-symmetry Keggin isomers, and to determine structural factors that control isomer selectivity as well as isomer labilization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Soret effect on hydrogen and helium retention in PFC tungsten under ELM-like conditions

Abstract In our previous work, we have demonstrated using nonequilibrium molecular-dynamics simulations that the fluxes of helium and self-interstitial atoms in the presence of a thermal gradient in tungsten are directed opposite to the heat flux, indicating that species transport is governed by a Soret effect, namely, thermal-gradient-driven diffusion, characterized by a negative heat of transport that drives species transport uphill, i.e. from the cooler to the hot regions of the tungsten sample. In this work, the findings of our thermal and species transport analysis have been implemented in our cluster-dynamics code, Xolotl, which has been used to compute temperature and species profiles over spatiotemporal scales representative of plasma-facing component (PFC) tungsten under typical reactor operating conditions, including extreme heat loads at the plasma-facing surface characteristic of plasma instabilities that induce edge localized modes (ELMs). We demonstrate that the steady-state species profiles, when properly accounting for the Soret effect, vary significantly from those where temperature-gradient-driven transport is not accounted for and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten. Although our cluster-dynamics simulations do not yet include self-clustering of helium or hydrogen blister formation, our simulation results show that the Soret effect substantially reduces helium and hydrogenic species retention inside PFC tungsten.

36 MATERIALS SCIENCE↗

Advancing Understanding of Fission Gas Behavior in Nuclear Fuel through Leadership Class Computing

This report describes the outcome of the pilot SciDAC project jointly funded by the DOE Office of Nuclear Energy and the DOE Office of Science, Advanced Scientific Computing Research to advance understanding of fission gas behavior in nuclear fuel. The computational approach involved utilizing both a bottom up and top down multiscale materials modeling to evaluate the processes of xenon fission gas diffusion, xenon clustering and precipitation and the mechanisms of fission gas resolution by energetic fission fragments in uranium dioxide nuclear fuel. The computer science and applied mathematics research activities were focused on optimizing the uranium – oxygen – xenon interatomic potential implementation in the LAMMPS molecular dynamics code for large-scale simulations of fission gas resolution, in addition to improvements in visualization of accelerated molecular dynamics simulations and performance optimization of the Xolotl continuum cluster dynamics with a new grouping scheme to reduce the computational cost of modeling very large fission gas bubbles containing more than 10,000 xenon atoms. The most impactful result from this project stems from a systematic series of atomistic molecular dynamics to evaluate fission gas resolution from xenon bubbles that, combined with the fission product yields and consideration of electron-phonon coupling, informed a first principles model that predicts the magnitude of fission gas resolution in good agreement with empirical estimates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Phonon localization in binary alloys with diagonal and off-diagonal disorder: A cluster Green's function approach

We report the development and application of a method for carrying out computational investigations of the effects of mass and force-constant (FC) disorder on phonon spectra. The method is based on the recently developed typical medium dynamical cluster approach, which is a Green's function approach. Excellent quantitative agreement with previous exact diagonalization results establishes the veracity of the method. Application of the method to a model system of binary mass and an FC-disordered system leads to several findings. A narrow resonance, significantly below the van Hove singularity, that has been termed the boson peak, is seen to emerge for low soft particle concentrations. We show, using the typical phonon spectrum, that the states constituting the boson peak cross over from being completely localized to being extended as a function of increasing soft particle concentration. In general, an interplay of mass and FC disorder is found to be cooperative in nature, enhancing phonon localization over all frequencies. However, for a certain range of frequencies, and depending on material parameters, FC disorder can delocalize the states that were localized by mass disorder, and vice versa. Modeling vacancies as weakly bonded sites with vanishing mass, we find that vacancies, even at very low concentrations, are extremely effective in localizing phonons. Thus, inducing vacancies is proposed as a promising route for efficient thermoelectrics. Lastly, we use model parameters corresponding to the alloy system, Ni 1- x Pt x , and we show that mass disorder alone is insufficient to explain the pseudogap in the phonon spectrum; the concomitant presence of FC disorder is necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

THE SWIFT X-RAY TELESCOPE CLUSTER SURVEY. III. CLUSTER CATALOG FROM 2005-2012 ARCHIVAL DATA

We present the Swift X-ray Cluster Survey (SWXCS) catalog obtained using archival data from the X-ray telescope (XRT) on board the Swift satellite acquired from 2005 February to 2012 November, extending the first release of the SWXCS. The catalog provides positions, soft fluxes, and, when possible, optical counterparts for a flux-limited sample of X-ray group and cluster candidates. We consider the fields with Galactic latitude |b| > 20° to avoid high H I column densities. We discard all of the observations targeted at groups or clusters of galaxies, as well as particular extragalactic fields not suitable to search for faint extended sources. We finally select ∼3000 useful fields covering a total solid angle of ∼400 deg{sup 2}. We identify extended source candidates in the soft-band (0.5-2 keV) images of these fields using the software EXSdetect, which is specifically calibrated for the XRT data. Extensive simulations are used to evaluate contamination and completeness as a function of the source signal, allowing us to minimize the number of spurious detections and to robustly assess the selection function. Our catalog includes 263 candidate galaxy clusters and groups down to a flux limit of 7 × 10{sup –15} erg cm{sup –2} s{sup –1} in the soft band, and the logN-logS is in very good agreement with previous deep X-ray surveys. The final list of sources is cross-correlated with published optical, X-ray, and Sunyaev-Zeldovich catalogs of clusters. We find that 137 sources have been previously identified as clusters in the literature in independent surveys, while 126 are new detections. Currently, we have collected redshift information for 158 sources (60% of the entire sample). Once the optical follow-up and the X-ray spectral analysis of the sources are complete, the SWXCS will provide a large and well-defined catalog of groups and clusters of galaxies to perform statistical studies of cluster properties and tests of cosmological models.

79 ASTRONOMY AND ASTROPHYSICS↗

Revealing EDL-driven reduction mechanisms in binary, ternary, and quaternary fluorinated electrolytes via an integrated MD–DFT–ML framework

Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.

DFT-MD-ML workflow↗

Unitary Selective Coupled-Cluster Method

Simulating molecules using the Variational Quantum Eigensolver method is one of the promising applications for NISQ-era quantum computers. Designing an efficient ansatz to represent the electronic wave function is crucial in such simulations. Standard unitary coupled-cluster with singles and doubles (UCCSD) ansatz tends to have a large number of insignificant terms that do not lower the energy of the system. In this work, we present a unitary selective coupled-cluster method, a way to construct a unitary coupled-cluster ansatz iteratively using a selection procedure with excitations up to fourth order. This approach uses the electronic Hamiltonian matrix elements and the amplitudes for excitations already present in the ansatz to find the important excitations of higher order and to add them to the ansatz. The important feature of the method is that it systematically reduces the energy error with increasing ansatz size for a set of test molecules. {The main advantage of the proposed method is that the effort to increase the ansatz does not require any additional measurements on a quantum computer.}

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comprehensive Benchmark Suite for Simulation of Particle Laden Flows Using the Discrete Element Method with Performance Profiles from the Multiphase Flow with Interface eXchanges (MFiX) Code

Five benchmark problems are developed and simulated with the computational fluid dynamics and discrete element model code MFiX. The benchmark problems span dilute and dense regimes, consider statistically homogeneous and inhomogeneous (both clusters and bubbles) particle concentrations and a range of particle and fluid dynamic computational loads. Several variations of the benchmark problems are also discussed to extend the computational phase space to cover granular (particles only), bidisperse and heat transfer cases. A weak scaling analysis is performed for each benchmark problem and, in most cases, the scalability of the code appears reasonable up to approx. 103 cores. Profiling of the benchmark problems indicate that the most substantial computational time is being spent on particle-particle force calculations, drag force calculations and interpolating between discrete particle and continuum fields. Hardware performance analysis was also carried out showing significant Level 2 cache miss ratios and a rather low degree of vectorization. These results are intended to serve as a baseline for future developments to the code as well as a preliminary indicator of where to best focus performance optimizations.

97 MATHEMATICS AND COMPUTING↗

Processing unit performance projection using dynamic hardware behaviors

Methods for estimating accelerator performance for dynamic hardware behaviors are disclosed. Computer program code to be executed on a first processing unit is received, and an execution of the computer code on the first processing unit is monitored to determine a plurality of performance characteristics. A plurality of dynamic hardware behaviors is determined by applying a clustering algorithm to the performance characteristics, and an equivalent accelerator portion of computer code to be executed on a second processing unit is generated by translating a set of instructions in a first portion of computer code corresponding to a first one of the plurality of dynamic hardware behaviors to an equivalent set of instructions to be executed on the second processing unit. An estimated measure of performance for executing the equivalent accelerator portion on the second processing unit is determined for the first one of the plurality of dynamic hardware behaviors.

97 MATHEMATICS AND COMPUTING↗

Challenge Online Time Series Clustering For Demand Response A Theory to Break the ‘Curse of Dimensionality'

The advent of smart meters and advanced communication infrastructures catalyzes numerous smart grid applications such as dynamic demand response, and paves the way to solve challenging research problems in sustainable energy consumption. The space of solution possibilities are restricted primarily by the huge amount of generated data requiring considerable computational resources and efficient algorithms. To overcome this Big Data challenge, data clustering techniques have been proposed. Current approaches however do not scale in the face of the “increasing dimensionality” problem where a cluster point is represented by the entire customer consumption time series. To overcome this aspect we first rethink the way cluster points are created and designed, and then design an efficient online clustering technique for demand response (DR) in order to analyze high volume, high dimensional energy consumption time series data at scale, and on the fly. Our online algorithm is randomized in nature, and provides optimal performance guarantees in a computationally efficient manner. Unlike prior work we (i) study the consumption properties of the whole population simultaneously rather than developing individual models for each customer separately, claiming it to be a ‘killer’ approach that breaks the “curse of dimensionality” in online time series clustering, and (ii) provide tight performance guarantees in theory to validate our approach. Our insights are driven by the field of sociology, where collective behavior often emerges as the result of individual patterns and lifestyles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Three dimensional cluster analysis for atom probe tomography using Ripley’s K-function and machine learning

The size and structure of spatial molecular and atomic clustering can significantly impact material properties and is therefore important to accurately quantify. Ripley’s K-function (K(r)), a measure of spatial correlation, can be used to perform such quantification when the material system of interest can be represented as a marked point pattern. This work demonstrates how machine learning models based on K (r)-derived metrics can accurately estimate cluster size and intra-cluster density in simulated three dimensional (3D) point patterns containing spherical clusters of varying size; over 90% of model estimates for cluster size and intra-cluster density fall within 11% and 18% error of the true values, respectively. These K (r)-based size and density estimates are then applied to an experimental APT reconstruction to characterize MgZn clusters in a 7000 series aluminum alloy. Here we find that the estimates are more accurate, consistent, and robust to user interaction than estimates from the popular maximum separation algorithm. Using K (r) and machine learning to measure clustering is an accurate and repeatable way to quantify this important material attribute.

36 MATERIALS SCIENCE↗

Structure, stability, and dissociation of small ionic silicon oxide clusters [SiO{sub n}{sup +}(n = 3, 4)]: Insight from density functional and topological exploration

The structures, energies, isomerization, and decomposition pathways of small ionic silicon oxide clusters, SiO{sub n}{sup +} (n = 3, 4), on doublet and quartet energy surfaces are investigated by density functional theory. New structural isomers of these ionic clusters have been obtained with this systematic study. The energy ordering of the isomeric cluster ions on doublet spin surface is found to follow the same general trend as that of the neutral ones, while it differs on the quartet surface. Our computational results reveal the energetically most preferred decomposition pathways of the ionic clusters on both spin surfaces. To comprehend the reaction mechanism, bonding evolution theory has also been employed using atoms in molecules formalism. The possible reasons behind the structural deformation of some isomers on quartet surface have also been addressed. Our results are expected to provide important insight into the decomposition mechanism and relative stability of the SiO{sub n}{sup +} clusters on both the energy surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗