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,549 records · Page 86

Improving Memory Error Handling Using Linux

As supercomputers continue to get faster and more powerful in the future, they will also have more nodes. If nothing is done, then the amount of memory in supercomputer clusters will soon grow large enough that memory failures will be unmanageable to deal with by manually replacing memory DIMMs. "Improving Memory Error Handling Using Linux" is a process oriented method to solve this problem by using the Linux kernel to disable (offline) faulty memory pages containing bad addresses, preventing them from being used again by a process. The process of offlining memory pages simplifies error handling and results in reducing both hardware and manpower costs required to run Los Alamos National Laboratory (LANL) clusters. This process will be necessary for the future of supercomputing to allow the development of exascale computers. It will not be feasible without memory error handling to manually replace the number of DIMMs that will fail daily on a machine consisting of 32-128 petabytes of memory. Testing reveals the process of offlining memory pages works and is relatively simple to use. As more and more testing is conducted, the entire process will be automated within the high-performance computing (HPC) monitoring software, Zenoss, at LANL.

97 MATHEMATICS AND COMPUTING↗

Biological Electron Transfer and Catalysis Energy Frontier Research Center

The overarching objective of the BETCy EFRC project is to investigate the mechanisms and structural basis controlling electron bifurcation, electron-ion coupling, and redox catalysis in model enzymes. This detailed understanding of how biological systems control and exploit electron flow will be used to develop modular biochemical conversions for the production of hydrocarbon and hydrogen biofuels. This project will conduct an in-depth study of metalloenzyme structure-function relationships that enable retuning of catalytic sites, electron-transfer thermodynamics, and the evolution of tailor-designed chemical reactions involving two model enzymes: hydrogenase and nitrogenase. The center is organized into three major thrust areas: Thrust 1. Energy conservation and electron flow; Thrust 2. ATP-coupled electron transfer; and Thrust 3. Atomic level determinants of enzymatic redox properties and their relationship to catalytic bias. We will develop a collective knowledge in metalloenzymes by applying physical science and computational tools to characterize biochemical reactions catalyzed specifically by multi-subunit enzymes harboring iron-sulfur clusters and flavin cofactors, as models for redox reactions throughout biology. Understanding these mechanisms is central to overcoming the key thermodynamic barriers currently limiting production of reduced products and biofuels. The research plan of the BETCy EFRC has been designed to provide the basis to exploit unique biochemical mechanisms that have yet to be explored substantively in the context of bioenergy but have the potential for game-changing advancement. The successful outcome of this research will be a collection of fundamental principles that serve as a blueprint to enable the tailored re-design of biological systems and enzymes to control matter and energy at the level of electrons and molecules to provide the foundations to create new energy technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Multistage Stochastic Transmission Expansion Algorithm for Wide-Area Planning under Uncertainty

The overall objective for this project was to develop and demonstrate a set of methods for solving the transmission investment problem for a large network considering many possible scenarios of future conditions and multiple decision points when investments can be made. Project sub-objectives achieved this goal through a succession of extending the methods to apply to problems with increasing complexity or additional features, including the number of decision points, whether generation and transmission are co-optimized, and whether AC or DC power flow is used. A transmission model was developed for the Western Electric Coordinating Council (WECC) region, the high-voltage transmission system that serves the western third of the continental U.S. Using a dataset provided by WECC and by researchers from John Hopkins University, we have validated and demonstrated the model and used it to compare the new method for solving multi-stage stochastic transmission planning to several state-of-the-art techniques. The project has resulted in several key outcomes and achievements: The covariance-based method for choosing a small set of hours to represent short-term variability has superior performance in terms of accuracy to existing methods, including K-means clustering and Importance Sampling; The combined partitioning method for long-term uncertainty with the nested clustering approach for choosing representative hours for each long-term group has superior accuracy for equivalent computational effort compared with existing methods; Using the partitioning/clustering method combined with Sample Average Approximation provides both statistical bounds on the quality of the solution and at the same time, a complete investment plan for all contingencies in the full uncertainty set; no existing methods can provide both at the same time; The method is demonstrated to work well for choosing both transmission and generation investments; A variant on the method allows for both scenario selection and simultaneous correction for the error from the DC power flow approximation to provide a tractable method for AC power flow-based transmission planning under uncertainty; The method applied to the WECC case study demonstrates the additional value to the system operator and the consumer of identifying flexible investment options in the near-term decisions. In particular, the case study exhibits significant option value in postponing some transmission additions that appear useful but in some long-term system states create new congestion problems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Biological Electron Transfer and Catalysis EFRC

The focus of the Biological Electron Transfer and Catalysis (BETCy) EFRC research is elucidating mechanisms of conversion of electrochemical potential into chemical bond energy and is organized into two integrated Thrusts: Electron Bifurcation Catalysis and Electron Transfer Kinetics. We are developing a collective knowledge of metalloenzymes as models for redox reactions by applying physical science and computational tools to characterize biochemical reactions catalyzed by multi- subunit enzymes harboring arrays of iron-sulfur clusters and flavin cofactors. Understanding these mechanisms is central to overcoming the thermodynamic barriers that currently limit production of reduced products and fuels. Biological systems have elegant strategies for converting electrochemical potential energy into chemical bond energy (e.g., C-H, H-H, and N-H) stored in reduced compounds that can serve as advanced biofuels. One significant limitation for the production of highly reduced compounds is that their production in natural and industrial processes relies on low oxidation-reduction potential chemicals as feedstocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THE FORMATION OF SECONDARY STELLAR GENERATIONS IN MASSIVE YOUNG STAR CLUSTERS FROM RAPIDLY COOLING SHOCKED STELLAR WINDS

We study a model of rapidly cooling shocked stellar winds in young massive clusters and estimate the circumstances under which secondary star formation, out of the reinserted winds from a first stellar generation (1G), is possible. We have used two implementations of the model: a highly idealized, computationally inexpensive, spherically symmetric semi-analytic model, and a complex, three-dimensional radiation-hydrodynamic, simulation; they are in a good mutual agreement. The results confirm our previous findings that, in a cluster with 1G mass 10{sup 7} M {sub ⊙} and half-mass–radius 2.38 pc, the shocked stellar winds become thermally unstable, collapse into dense gaseous structures that partially accumulate inside the cluster, self-shield against ionizing stellar radiation, and form the second generation (2G) of stars. We have used the semi-analytic model to explore a subset of the parameter space covering a wide range of the observationally poorly constrained parameters: the heating efficiency, η {sub he}, and the mass loading, η {sub ml}. The results show that the fraction of the 1G stellar winds accumulating inside the cluster can be larger than 50% if η {sub he} ≲ 10%, which is suggested by the observations. Furthermore, for low η {sub he}, the model provides a self-consistent mechanism predicting 2G stars forming only in the central zones of the cluster. Finally, we have calculated the accumulated warm gas emission in the H30 α recombination line, analyzed its velocity profile, and estimated its intensity for super star clusters in interacting galaxies NGC4038/9 (Antennae) showing that the warm gas should be detectable with ALMA.

79 ASTRONOMY AND ASTROPHYSICS↗

A New CIERRA Gridded Product Climatology for LIS/OTD

The GLM-CIERRA processing generates two products: A level-2 cluster feature dataset that merges flashes that should have been a single flash following the clustering algorithm flash definition (for GLM, this is due to the LCFA thresholds); A set of level-3 gridded products including Flash Extent Density, Convective Probability, etc. The CIERRA Level-3 products are generated on a 0.1 degree grid. This is suboptimal given the spatial variations in GLM pixel size AND the size of gridpoints in the output grid. As a stepping-stone towards developing improved GLM-CIERRA grids, I am testing pixel matching techniques on the LIS/OTD data. The goal is to maintain accuracy while vastly improving computational efficiency. The CIERRA reclustering codes are also being applied to LIS / OTD to merge flashes split by the “first fit” clustering technique. The end result will be standardized grids at the nominal resolution of each instrument (5 km for LIS, 10 km for OTD) that take into account the pixel geometry for each event that comprises a given reclustered flash.

54 ENVIRONMENTAL SCIENCES↗

The role of silicon, vacancies, and strain in carbon distribution for low temperature bainite

Here, we investigated the phenomenon of carbon supersaturation and carbon clustering in bainitic ferrite with atom probe tomography (APT) and ab-initio density functional theory (DFT) calculations. The experimental results show a homogeneous distribution of silicon in the microstructure, which contains both ferrite and retained austenite. This distribution is mimicked well by the computational approach. In addition, an accumulation of C in certain regions of the bainitic ferrite with C concentrations up to 13 at % is observed. Based on the DFT results, these clusters are explained as strained, tetragonal regions in the ferritic bainite, in which the solution enthalpy of C can reach large, negative values. It seems that Si itself only has a minor influence on this phenomenon.

36 MATERIALS SCIENCE↗

Trend Analysis in Biomedical Texts via Vector Space Model Synonym Extraction

Accurate synonym extraction is a highly desired tool in computational semantics and natural language processing (NLP). Intrinsic to many current language models is the Vector Space Model (VSM) whereby linguistic properties are modeled by real vectors of small length compared to the document size. Here, several vector models are discussed. Utilizing the open source Word2Vec tool, temporal trends in biomedical abstracts from MEDLINE have been analyzed via a unique transformation matrix method. Additionally, accuracy of the word vector method and clustering of temporal trends are discussed.

59 BASIC BIOLOGICAL SCIENCES↗

Performance of the fusion code GYRO on four generations of Cray computers

GYRO is a code used for the direct numerical simulation of plasma microturbulence. It has been ported to a variety of modern MPP platforms including several modern commodity clusters, IBM SPs, and Cray XC, XT, and XE series machines. We briefly describe the mathematical structure of the equations, the data layout, and the redistribution scheme. Also, while the performance and scaling of GYRO on many of these systems has been shown before, here we show the comparative performance and scaling on four generations of Cray supercomputers including the newest addition - the Cray XC30. The more recently added hybrid OpenMP/MPI imple- mentation also shows a great deal of promise on custom HPC systems that utilize fast CPUs and proprietary interconnects. Four machines of varying sizes were used in the experiment, all of which are located at the National Institute for Computational Sciences at the University of Tennessee at Knoxville and Oak Ridge National Laboratory. The advantages, limitations, and performance of using each system are discussed.

Fahey, Mark R↗

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗

Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework: Theory, implementation, and validation against multi-reference methods

In this work, we report an implementation of non-adiabatic coupling (NAC) forces within the equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) framework via the summed-state approach. Using illustrative examples, we compare NAC forces computed with EOM-CCSD and multi-reference (MR) wave functions (for selected cases, we also consider configuration interaction singles). In addition to the magnitude of the NAC vectors, we analyze their direction, which is important for the calculations of the rate of non-adiabatic transitions. Our benchmark set comprises three doublet radical-cations (hexatriene, cyclohexadiene, and uracil), neutral uracil, and sodium-doped ammonia clusters. When the characters of the states agree among different methods, we observe good agreement between the respective NAC vectors, both in the Franck-Condon region and away. In the cases of large discrepancies between the methods, the disagreement can be attributed to the difference in the states’ character, which, in some cases, is very sensitive to electron correlation, both within single-reference and multi-reference frameworks. Finally, the numeric results confirm that the accuracy of NAC vectors depends critically on the quality of the underlying wave functions. Within their domain of applicability, EOM-CC methods provide a viable alternative to MR approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using anion photoelectron spectroscopy of cluster models to gain insights into mechanisms of catalyst-mediated H 2 production from water

We report that metal oxide cluster models of catalyst materials offer a powerful platform for probing the molecular-scale features and interactions that govern catalysis. This perspective gives an overview of studies implementing the combination of anion photoelectron (PE) spectroscopy and density functional theory calculations toward exploring cluster models of metal oxides and metal-oxide supported Pt that catalytically drive the hydrogen evolution reaction (HER) or the water–gas shift reaction. The utility in the combination of these experimental and computational techniques lies in our ability to unambiguously determine electronic and molecular structures, which can then connect to results of reactivity studies. In particular, we focus on the activity of oxygen vacancies modeled by suboxide clusters, the critical mechanistic step of forming proximal metal hydride and hydroxide groups as a prerequisite for H2 production, and the structural features that lead to trapped dihydroxide groups. The pronounced asymmetric oxidation found in heterometallic group 6 oxides and near-neighbor group 5/group 6 results in higher activity toward water, while group 7/group 6 oxides form very specific stoichiometries that suggest facile regeneration. Studies on the trans-periodic combination of cerium oxide and platinum as a model for ceria supported Pt atoms and nanoparticles reveal striking negative charge accumulation by Pt, which, combined with the ionic conductivity of ceria, suggests a mechanism for the exceptionally high activity of this system towards the water–gas shift reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the accuracy of density functional theory and wave function methods for calculating vertical ionization energies

The best practice in computational methods for determining vertical ionization energies (VIEs) is assessed, via reference to experimentally determined VIEs that are corroborated by highly accurate coupled-cluster calculations. These reference values are used to benchmark the performance of density functional theory (DFT) and wave function methods: Hartree-Fock theory, second-order Møller-Plesset perturbation theory, and Electron Propagator Theory (EPT). The core test set consists of 147 small molecules. An extended set of six larger molecules, from benzene to hexacene, is also considered to investigate the dependence of the results on molecule size. The closest agreement with experiment is found for ionization energies obtained from total energy difference calculations. In particular, DFT calculations using exchange-correlation functionals with either a large amount of exact exchange or long-range correction perform best. The results from these functionals are also the least sensitive to an increase in molecule size. In general, ionization energies calculated directly from the orbital energies of the neutral species are less accurate and more sensitive to an increase in molecule size. For the single-calculation approach, the EPT calculations are in closest agreement for both sets of molecules. For the orbital energies from DFT functionals, only those with long-range correction give quantitative agreement with dramatic failing for all other functionals considered. The results offer a practical hierarchy of approximations for the calculation of vertical ionization energies. In addition, the experimental and computational reference values can be used as a standardized set of benchmarks, against which other approximate methods can be compared.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using first principles calculations to interpret XANES experiments: extracting the size-dependence of the ( p , T ) phase diagram of sub-nanometer Cu clusters in an O 2 environment

Here, we have used ab initio density functional theory together with ab initio atomistic thermodynamics, and in situ X-ray absorption near edge spectroscopy (XANES) experiments, to study the oxidation of sub-nanometer clusters of Cu n O x supported on a hydroxylated amorphous alumina substrate in an O 2 -rich environment. We obtain ( p , T ) phase diagrams; these differ notably for the nanoclusters compared to the bulk. Both theory and experiment suggest that in the presence of oxygen, the cluster will oxidize from its elemental state to the oxidized state as temperature decreases. We obtain a clear trend for the transition Cu n → Cu n O n/2 : we see that smaller the cluster, greater is the tendency toward oxidation. However, we do not see a monotonic size-dependent trend for the transition Cu n O n/2 → Cu n O n . We suggest that theoretically computed Bader charges constitute a simple yet quantitative way to align experimental measures of XANES edges with theoretical calculations, so as to yield oxidation states for nanoclusters. Our results have important implications for the use of small clusters in fields such as nanocatalysis and nanomedicine.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Predicting phase behavior of grain boundaries with evolutionary search and machine learning

The study of grain boundary phase transitions is an emerging field until recently dominated by experiments. The major bottleneck in the exploration of this phenomenon with atomistic modeling has been the lack of a robust computational tool that can predict interface structure. Here we develop a computational tool based on evolutionary algorithms that performs efficient grand-canonical grain boundary structure search and we design a clustering analysis that automatically identifies different grain boundary phases. Its application to a model system of symmetric tilt boundaries in Cu uncovers an unexpected rich polymorphism in the grain boundary structures. We find new ground and metastable states by exploring structures with different atomic densities. Our results demonstrate that the grain boundaries within the entire misorientation range have multiple phases and exhibit structural transitions, suggesting that phase behavior of interfaces is likely a general phenomenon.

36 MATERIALS SCIENCE↗

Fluorescence Color by Data-Driven Design of Genomic Silver Clusters

DNA nucleobase sequence controls the size of DNA-stabilized silver clusters, leading to their well-known yet little understood sequence-tuned colors. The enormous space of possible DNA sequences for templating clusters has challenged the understanding of how sequence selects cluster properties and has limited the design of applications that employ these clusters. Here, we investigate the genomic role of DNA sequence for fluorescent silver clusters using a data-driven approach. Employing rapid parallel silver cluster synthesis and fluorimetry, we determine the fluorescence spectra of silver cluster products stabilized by 1432 distinct DNA oligomers. By applying pattern recognition algorithms to this large experimental data set, we discover certain DNA base patterns, or “motifs,” that correlate to silver clusters with similar fluorescence spectra. These motifs are employed in machine learning classifiers to predictively design DNA template sequences for specific fluorescence color bands. Our method improves selectivity of templates by 330% for silver clusters with peak emission wavelengths beyond 660 nm. The discovered base motifs also provide physical insights into how DNA sequence controls silver cluster size and color. Finally, this predictive design approach for color of DNA-stabilized silver clusters exhibits the potential of machine learning and data mining to increase the precision and efficiency of nanomaterials design, even for a soft-matter-inorganic hybrid system characterized by an extremely large parameter space.

59 BASIC BIOLOGICAL SCIENCES↗

Computational stochastic model of ions implantation

Implantation flux ions into crystal leads to phase transition /PT/ 1-st kind. Damaging lattice is associated with processes clustering vacancies and gaseous bubbles as well their brownian motion. System of stochastic differential equations /SDEs/ Ito for evolution stochastic dynamical variables corresponds to the superposition Wiener processes. The kinetic equations in partial derivatives /KE/, Kolmogorov-Feller and Einstein-Smolukhovskii, were formulated for nucleation into lattice of weakly soluble gases. According theory, coefficients of stochastic and kinetic equations uniquely related. Radiation stimulated phase transition are characterized by kinetic distribution functions /DFs/ of implanted clusters versus their sizes and depth of gas penetration into lattice. Macroscopic parameters of kinetics such as the porosity and stress calculated in thin layers metal/dielectric due to Xe{sup ++} irradiation are attracted as example. Predictions of porosity, important for validation accumulation stresses in surfaces, can be applied at restoring of objects the cultural heritage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

libEnsemble: A complete Python toolkit for dynamic ensembles of calculations

Almost all science and engineering applications eventually stop scaling: their runtime no longer decreases as available computational resources increase. Therefore, many applications will struggle to efficiently use emerging extreme-scale high-performance, parallel, and distributed systems. libEnsemble is a complete Python toolkit and workflow system for intelligently driving ensembles of experiments or simulations at massive scales. It enables and encourages multidisciplinary design, decision, and inference studies portably running on laptops, clusters, and supercomputers.

97 MATHEMATICS AND COMPUTING↗