Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster computer”

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 721 records · Page 40

Optimization of the Coupled Cluster Implementation in NWChem on Petascale Parallel Architectures

Coupled cluster singles and doubles (CCSD) algorithm has been optimized in NWChem software package. This modification alleviated the communication bottleneck and provided from 2- to 5-fold speedup in the CCSD iteration time depending on the problem size and available memory. Sustained 0.60 petaflop/sec performance on CCSD(T) calculation has been obtained on NCSA Blue Waters. This number included all stages of the calculation from initialization till termination, iterative computation of single and double excitations, and perturbative accounting for triple excitations. In the section of perturbative triples alone, the computation maintained 1.18 petaflop/sec performance level. CCSD computations have been performed on Guanine-Cytosine deoxydinucleotide monophosphate (GC-dDMP) to probe the conformational energy difference in DNA single strand in A- and B-conformations. The computation revealed significant discrepancy between CCSD and classical force fields in prediction of relative energy of A- and B-conformations of GC-dDMP.

calculations↗

Tensor-structured coupled cluster theory

We derive and implement a new way of solving coupled cluster equations with lower computational scaling. Our method is based on the decomposition of both amplitudes and two electron integrals, using a combination of tensor hypercontraction and canonical polyadic decomposition. While the original theory scales as O ( N 6 ) with respect to the number of basis functions, we demonstrate numerically that we achieve sub-millihartree difference from the original theory with O ( N 4 ) scaling. This is accomplished by solving directly for the factors that decompose the cluster operator. Here, the proposed scheme is quite general and can be easily extended to other many-body methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GROMACS 4.5: A high-throughput and highly parallel open source molecular simulation toolkit

In this study, molecular simulation has historically been a low-throughput technique, but faster computers and increasing amounts of genomic and structural data are changing this by enabling large-scale automated simulation of, for instance, many conformers or mutants of biomolecules with or without a range of ligands. At the same time, advances in performance and scaling now make it possible to model complex biomolecular interaction and function in a manner directly testable by experiment. These applications share a need for fast and efficient software that can be deployed on massive scale in clusters, web servers, distributed computing or cloud resources. As a result, we present a range of new simulation algorithms and features developed during the past 4 years, leading up to the GROMACS 4.5 software package. The software now automatically handles wide classes of biomolecules, such as proteins, nucleic acids and lipids, and comes with all commonly used force fields for these molecules built-in. GROMACS supports several implicit solvent models, as well as new free-energy algorithms, and the software now uses multithreading for efficient parallelization even on low-end systems, including windows-based workstations. Together with hand-tuned assembly kernels and state-of-the-art parallelization, this provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations.

97 MATHEMATICS AND COMPUTING↗

Ruprecht 147: A Paradigm of Dissolving Star Cluster

We employed recent Gaia/DR2 data to investigate the dynamical status of the nearby (300 pc), old (2.5 Gyr) open cluster Ruprecht 147. We found prominent leading and trailing tails of stars along the cluster orbit, which demonstrates that Ruprecht 147 is losing stars at fast pace. Star counts indicate the cluster has a core radius of 33.3 arcmin and a tidal radius of 137.5 arcmin. The cluster also possesses an extended corona, which cannot be reproduced by a simple King model. We computed the present-day cluster mass using its luminosity and mass function, and derived an estimate of 234 ± 52 M {sub ⊙}. We also estimated the cluster original mass using available recipes extracted from N-body simulations obtaining a mass at birth of 50,000 ± 6500 M {sub ⊙}. Therefore dynamical mass loss, mostly caused by tidal interaction with the Milky Way, reduced the cluster mass by about 99%. We then conclude that Ruprecht 147 is rapidly dissolving into the general Galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Gaussian mixture models as automated particle classifiers for fast neutron detectors

Pulse shape discrimination (PSD) is the task of classifying electronic pulse shapes for different particle types such as gamma rays and fast neutrons interacting in scintillators and read out by photo sensitive detectors. This field has been limited in its adoption of techniques found in the statistical learning community. Methods initially employed in the 1960s for analog electronic circuitry persist in the current PSD literature describing operations performed on digitized pulses, which are amenable to statistical rigor. Despite vast amounts of data collected at low energy levels, traditional PSD methods are unable to discriminate particles below a certain threshold. As such, in this work, Gaussian mixture models (GMMs) are used as a clustering technique for fast neutron detection in the absence of labeled data. GMMs yield improvements spanning the energy spectrum in a desirably efficient, unsupervised fashion. An extension, the Dirichlet Process GMM, provides further flexibility and classification improvement.

97 MATHEMATICS AND COMPUTING↗

Charliecloud

Charliecloud is a set of scripts to let users run a virtual cluster of virtual machines (VMs) on a desktop or supercomputer. Key functions include: 1. Creating (typically by installing an operating system from vendor media) and updating VM images; 2. Running a single VM; 3. Running multiple VMs in a virtual cluster. The virtual machines can talk to one another over the network and (in some cases) the outside world. This is accomplished by calling external programs such as QEMU and the Virtual Distributed Ethernet (VDE) suite. The goal is to let users have a virtual cluster containing nodes where they have privileged access, while isolating that privilege within the virtual cluster so it cannot affect the physical compute resources. Host configuration enforces security; this is not included in Charliecloud, though security guidelines are included in its documentation and Charliecloud is designed to facilitate such configuration. Charliecloud manages passing information from host computers into and out of the virtual machines, such as parameters of the virtual cluster, input data specified by the user, output data from virtual compute jobs, VM console display, and network connections (e.g., SSH or X11). Parameters for the virtual cluster (number of VMs, RAM and disk per VM, etc.) are specified by the user or gathered from the environment (e.g., SLURM environment variables). Example job scripts are included. These include computation examples (such as a "hello world" MPI job) as well as performance tests. They also include a security test script to verify that the virtual cluster is appropriately sandboxed. Tests include: 1. Pinging hosts inside and outside the virtual cluster to explore connectivity; 2. Port scans (again inside and outside) to see what services are available; 3. Sniffing tests to see what traffic is visible to running VMs; 4. IP address spoofing to test network functionality in this case; 5. File access tests to make sure host access permissions are enforced. This test script is not a comprehensive scanner and does not test for specific vulnerabilities. Importantly, no information about physical hosts or network topology is included in this script (or any of Charliecloud); while part of a sensible test, such information is specified by the user when the test is run. That is, one cannot learn anything about the LANL network or computing infrastructure by examining Charliecloud code.

Priedhorsky, Reid↗

Charliecloud

Charliecloud is a set of scripts to let users run a virtual cluster of virtual machines (VMs) on a desktop or supercomputer. Key functions include: 1. Creating (typically by installing an operating system from vendor media) and updating VM images; 2. Running a single VM; 3. Running multiple VMs in a virtual cluster. The virtual machines can talk to one another over the network and (in some cases) the outside world. This is accomplished by calling external programs such as QEMU and the Virtual Distributed Ethernet (VDE) suite. The goal is to let users have a virtual cluster containing nodes where they have privileged access, while isolating that privilege within the virtual cluster so it cannot affect the physical compute resources. Host configuration enforces security; this is not included in Charliecloud, though security guidelines are included in its documentation and Charliecloud is designed to facilitate such configuration. Charliecloud manages passing information from host computers into and out of the virtual machines, such as parameters of the virtual cluster, input data specified by the user, output data from virtual compute jobs, VM console display, and network connections (e.g., SSH or X11). Parameters for the virtual cluster (number of VMs, RAM and disk per VM, etc.) are specified by the user or gathered from the environment (e.g., SLURM environment variables). Example job scripts are included. These include computation examples (such as a "hello world" MPI job) as well as performance tests. They also include a security test script to verify that the virtual cluster is appropriately sandboxed. Tests include: 1. Pinging hosts inside and outside the virtual cluster to explore connectivity; 2. Port scans (again inside and outside) to see what services are available; 3. Sniffing tests to see what traffic is visible to running VMs; 4. IP address spoofing to test network functionality in this case; 5. File access tests to make sure host access permissions are enforced. This test script is not a comprehensive scanner and does not test for specific vulnerabilities. Importantly, no information about physical hosts or network topology is included in this script (or any of Charliecloud); while part of a sensible test, such information is specified by the user when the test is run. That is, one cannot learn anything about the LANL network or computing infrastructure by examining Charliecloud code.

Priedhorsky, Reid↗

Uniformly Ordered Binary Decision Algorithm for Benchmark Experiment Correlations in Whisper Validation

When performing a validation exercise for determining the upper subcritical limit of a nuclear criticality safety application, an analyst should select and perform a statistical analysis on a population of benchmark experiments that are neutronically similar to the application. The size of this population should be sufficiently large such that the statistical analysis has a high degree of confidence that the bias plus bias uncertainty (calculational margin) has been accurately quantified. A complication arises because many benchmark experiments share common components, leading to correlations in their measured effective multiplication factors. Correlations between benchmark experiments within the population reduces its predictive power. This motivates the need for methods that consider benchmark experiment correlations and ensure adequate statistical significance of results. The Whisper code is a statistical analysis pack- age that incorporates nuclear data sensitivity coefficients from MCNP to assess benchmark experiment similarity and then performs an extreme-value analysis to estimate the bias plus bias uncertainty. The original methodology in Whisper does not consider the effect of benchmark experiment correlations when making this estimation, and this summary proposes the uniformly ordered binary decision algorithm to address this shortcoming. The original methodology in Whisper computes similarity coefficients ck for an application compared to all benchmark experiments in its library and develops weighting factors for a selected population proportional to the ck values. The methodology can be interpreted as statistically emulating a validation exercise for a particular application where the weighting factors may be viewed as the likelihood that an analyst would include a particular benchmark experiment within the population. The effective sample size of the population is the expected or mean number of benchmark experiments in the population. The uniformly ordered binary decision algorithm identifies clusters of correlated benchmark experiments within the population and then computes adjusted weighting factors based on the magnitude of the correlation coefficients within the cluster to compute a reduced effective sample size accounting for the lower information content because of correlations. Benchmark experiments within the cluster are ordered randomly with equal probability and probabilistic decisions are made as to whether a benchmark. experiment within the cluster should treated as redundant with a previous one; if two redundant benchmark experiments are included, then the conservative worst case bias plus bias uncertainty is used and the pair is counted as a single benchmark experiment in the population. Results are provided for HEU solutions in a research version of the Whisper software using benchmark experiment correlations provided by DICE, the Database for the International Criticality Safety Benchmark Evaluation Project (ICSBEP). These show that there can be a significant increase in the bias plus bias uncertainty because the effective sample size is reduced, and therefore the algorithm, needing to meet sample size requirements, expands the benchmark experiment population by accepting less similar benchmark experiments that would have otherwise not been included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A computational study of Al n - and Al n -1 1Pt - clusters: the effects of doping and a uniform tuning gauge for single-atom nanocatalysts

Results of density functional theory calculations on Al n - and Al n -1 1Pt - , n = 2-8, clusters are presented and analyzed. The analysis includes different structural forms of the clusters characterized in terms of binding energy, spin and symmetry, and a comparative evaluation of various properties of the two systems viewed as connected through a single-Pt substitutional doping and examined in terms of their respective most stable structures. Furthermore, the Al n -1 1Pt - clusters are then used as a paradigmatic (model) case of single-atom nanocatalysts, with Pt as the catalytic center and Al n-1 as its support, to implement a uniform descriptor for gauging the tuning effects of all the parameters ( "knobs ") of a nanocatalyst that include the identity of the active center and the material and size of its support.

97 MATHEMATICS AND COMPUTING↗

The antimony-group 11 chemical bond: Dissociation energies of the diatomic molecules CuSb, AgSb, and AuSb

The intermetallic molecules CuSb, AgSb, and AuSb were identified in the effusive molecular beam produced at high temperature under equilibrium conditions in a double-cell-like Knudsen source. Several gaseous equilibria involving these species were studied by mass spectrometry as a function of temperature in the overall range 1349–1822 K, and the strength of the chemical bond formed between antimony and the group 11 metals was for the first time measured deriving the following thermochemical dissociation energies (D{sub 0}{sup ∘}, kJ/mol): 186.7 ± 5.1 (CuSb), 156.3 ± 4.9 (AgSb), 241.3 ± 5.8 (AuSb). The three species were also investigated computationally at the coupled cluster level with single, double, and noniterative quasiperturbative triple excitations (CCSD(T)). The spectroscopic parameters were calculated from the potential energy curves and the dissociation energies were evaluated at the Complete Basis Set limit, resulting in an overall good agreement with experimental values. An approximate evaluation of the spin-orbit effect was also performed. CCSD(T) calculations were further extended to the corresponding group 11 arsenide species which are here studied for the first time and the following dissociation energies (D{sub 0}{sup ∘}, kJ/mol): 190 ± 10 (CuAs), 151 ± 10 (AgAs), 240 ± 15 (AuAs) are proposed. Taking advantage of the new experimental and computational information here presented, the bond energy trends along group 11 and 4th and 5th periods of the periodic table were analyzed and the bond energies of the diatomic species CuBi and AuBi, yet experimentally unobserved, were predicted on an empirical basis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formulation of Multicomponent Lattice Gas Model Cluster Expansions Parameterized on Ab Initio Data: An Introduction to the Ab Initio Mean-Field Augmented Lattice Gas Modeling Code

A formalism is presented for the construction of multicomponent lattice gas models parameterized with ab initio (typically density functional theory) data. The leave-multiple-out and the leave-one-out cross validation score are showcased, and practical algorithms are developed and implemented in a new code called the Ab initio Mean-field Augmented Lattice Gas Modeling (AMALGM) code. Here, the functionality of these algorithms is demonstrated with a fully worked out example using the O/Fe(100) system. AMALGM and the formalism on which it is based are envisioned as a surface-oriented lattice gas alternative to other cluster expansion codes that are typically geared toward bulk systems whose lateral interactions between components are parameterized using an Ising model. While the formalism is created within the context of surfaces, it is equally applicable to bulk materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parameters associated with design effect of child anthropometry indicators in small-scale field surveys

Background: Cluster surveys provide rapid but representative estimates of key nutrition indicators in humanitarian crises. For these surveys, an accurate estimate of the design effect is critical to calculate a sample size that achieves adequate precision with the minimum number of sampling units. This paper describes the variability in design effect for three key nutrition indicators measured in small-scale surveys and models the association of design effect with parameters hypothesized to explain this variability. Methods: 380 small-scale surveys from 28 countries conducted between 2006 and 2013 were analyzed. We calculated prevalence and design effect of wasting, underweight, and stunting for each survey as well as standard deviations of the underlying continuous Z-score distribution. Mean cluster size, survey location and year were recorded. To describe design effects, median and interquartile ranges were examined. Generalized linear regression models were run to identify potential predictors of design effect. Results: Median design effect was under 2.00 for all three indicators; for wasting, the median was 1.35, the lowest among the indicators. Multivariable linear regression models suggest significant, positive associations of design effect and mean cluster size for all three indicators, and with prevalence of wasting and underweight, but not stunting. Standard deviation was positively associated with design effect for wasting but negatively associated for stunting. Survey region was significant in all three models. Conclusions: This study supports the current field survey guidance recommending the use of 1.5 as a benchmark for design effect of wasting, but suggests this value may not be large enough for surveys with a primary objective of measuring stunting or underweight. The strong relationship between design effect and region in the models underscores the continued need to consider country- and locality-specific estimates when designing surveys. These models also provide empirical evidence of a positive relationship between design effect and both mean cluster size and prevalence, and introduces standard deviation of the underlying continuous variable (Z-scores) as a previously unexplored factor significantly associated with design effect. The magnitude and directionality of this association differed by indicator, underscoring the need for further investigation into the relationship between standard deviation and design effect.

60 APPLIED LIFE SCIENCES↗

Observational and Numerical Diagnostics of Galaxy Cluster Outer Regions

Aims. We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r(sub 200) and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. Results. As opposed to several recent results, we observe a steepening of the density profiles beyond approximately 0.3r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or gas clumping are in better agreement with the observed gas distribution. We note a systematic difference between cool-core and non-cool core clusters beyond approximately 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond approximately r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. Conclusions. The general trend of steepening density around the virial radius indicates that the shallow density profiles found in several recent works were probably obtained along particular directions (e.g., filaments) and are not representative of the typical behavior of clusters. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Therefore, a detailed treatment of gas cooling, star formation, clumping, and AGN feedback is required to construct realistic models of cluster outer regions.

Eckert, D.↗

Beyond the plane-parallel approximation for redshift surveys

Redshift -space distortions privilege the location of the observer in cosmological redshift surveys, breaking the translational symmetry of the underlying theory. This violation of statistical homogeneity has consequences for the modelling of clustering observables, leading to what are frequently called 'wide-angle effects'. We study these effects analytically, computing their signature in the clustering of the multipoles in configuration and Fourier space. We take into account both physical wide-angle contributions as well as the terms generated by the galaxy selection function. Similar considerations also affect the way power spectrum estimators are constructed. We quantify in an analytical way the biases that enter and clarify the relation between what we measure and the underlying theoretical modelling. The presence of an angular window function is also discussed. Motivated by this analysis, we present new estimators for the three dimensional Cartesian power spectrum and bispectrum multipoles written in terms of spherical Fourier-Bessel coefficients. In conclusion, we show how the latter have several interesting properties, allowing in particular a clear separation between angular and radial modes.

79 ASTRONOMY AND ASTROPHYSICS↗

The Gas Distribution in the Outer Regions of Galaxy Clusters

Aims. We present our analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We have exploited the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius, We stacked the density profiles to detect a signal beyond T200 and measured the typical density and scatter in cluster outskirts. We also computed the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compared our average density and scatter profiles with the results of numerical simulations. Results. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond approximately r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict density profiles that are too steep, whereas runs including additional physics and/ or treating gas clumping agree better with the observed gas distribution. We report high-confidence detection of a systematic difference between cool-core and non cool-core clusters beyond approximately 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond approximately r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only small differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the ENZO simulations. Conclusions. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Although their general behavior agrees more closely with the observations, simulations including cooling and star formation convert a large amount of gas into stars, which results in a low gas fraction with respect to the observations. Consequently, a detailed treatment of gas cooling, star formation, AGN feedback, and consideration of gas clumping is required to construct realistic models of the outer regions of clusters.

Eckert, D.↗

Computational strategies for three-dimensional flow simulations on distributed computer systems

This research effort is directed towards an examination of issues involved in porting large computational fluid dynamics codes in use within the industry to a distributed computing environment. This effort addresses strategies for implementing the distributed computing in a device independent fashion and load balancing. A flow solver called TEAM presently in use at Lockheed Aeronautical Systems Company was acquired to start this effort. The following tasks were completed: (1) The TEAM code was ported to a number of distributed computing platforms including a cluster of HP workstations located in the School of Aerospace Engineering at Georgia Tech; a cluster of DEC Alpha Workstations in the Graphics visualization lab located at Georgia Tech; a cluster of SGI workstations located at NASA Ames Research Center; and an IBM SP-2 system located at NASA ARC. (2) A number of communication strategies were implemented. Specifically, the manager-worker strategy and the worker-worker strategy were tested. (3) A variety of load balancing strategies were investigated. Specifically, the static load balancing, task queue balancing and the Crutchfield algorithm were coded and evaluated. (4) The classical explicit Runge-Kutta scheme in the TEAM solver was replaced with an LU implicit scheme. And (5) the implicit TEAM-PVM solver was extensively validated through studies of unsteady transonic flow over an F-5 wing, undergoing combined bending and torsional motion. These investigations are documented in extensive detail in the dissertation, 'Computational Strategies for Three-Dimensional Flow Simulations on Distributed Computing Systems', enclosed as an appendix.

Sankar, Lakshmi N.↗

The Gas Distribution in Galaxy Cluster Outer Regions

Aims. We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r200 and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. Results. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond approximately r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or treating gas clumping are in better agreement with the observed gas distribution. We report for the first time the high-confidence detection of a systematic difference between cool-core and non-cool core clusters beyond 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. Conclusions. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Although their general behavior is in better agreement with the observations, simulations including cooling and star formation convert a large amount of gas into stars, which results in a low gas fraction with respect to the observations. Consequently, a detailed treatment of gas cooling, star formation, AGN feedback, and taking into account gas clumping is required to construct realistic models of cluster outer regions.

Eckert, D.↗