Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel cluster”

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 235 records · Page 13

GUNDAM: a toolkit for fast spatial correlation functions in galaxy surveys

ABSTRACT We describe the capabilities of a new software package to calculate two-point correlation functions (2PCFs) of large galaxy samples. The code can efficiently estimate 3D/projected/angular 2PCFs with a variety of statistical estimators and bootstrap errors, and is intended to provide a complete framework (including calculation, storage, manipulation, and plotting) to perform this type of spatial analysis with large redshift surveys. Gundam implements a very fast skip list/linked list algorithm that efficiently counts galaxy pairs and avoids the computation of unnecessary distances. It is several orders of magnitude faster than a naive pair counter, and matches or even surpass other advanced algorithms. The implementation is also embarrassingly parallel, making full use of multicore processors or large computational clusters when available. The software is designed to be flexible, user friendly and easily extensible, integrating optimized, well-tested packages already available in the astronomy community. Out of the box, it already provides advanced features such as custom weighting schemes, fibre collision corrections and 2D correlations. Gundam will ultimately provide an efficient toolkit to analyse the large-scale structure ‘buried’ in upcoming extremely large data sets generated by future surveys.

Donoso, E.↗

The global star formation law by supernova feedback

ABSTRACT We address a simple model where the Kennicutt–Schmidt (KS) relation between the macroscopic densities of star formation rate (SFR, ρsfr) and gas (n) in galactic discs emerges from self-regulation of the SFR via supernova feedback. It arises from the physics of supernova bubbles, insensitive to the microscopic SFR recipe and not explicitly dependent on gravity. The key is that the filling factor of SFR-suppressed supernova bubbles self-regulates to a constant, f ∼ 0.5. Expressing the bubble fading radius and time in terms of n, the filling factor is $f\propto S\, n^{-s}$ with s ≃ 1.5, where S is the supernova rate density. A constant f thus refers to ρsfr ∝ n1.5, with a density-independent SFR efficiency per free-fall time ∼0.01. The self-regulation to f ∼ 0.5 and the convergence to a KS relation independent of the local SFR recipe are demonstrated in cosmological and isolated-galaxy simulations using different codes and recipes. In parallel, the spherical analysis of bubble evolution is generalized to clustered supernovae, analytically and via simulations, yielding s ≃ 1.5 ± 0.5. An analysis of photoionized bubbles about pre-supernova stars yields a range of KS slopes but the KS relation is dominated by the supernova bubbles. Superbubble blowouts may lead to an alternative self-regulation by outflows and recycling. While the model is oversimplified, its simplicity and validity in the simulations may argue that it captures the origin of the KS relation.

79 ASTRONOMY AND ASTROPHYSICS↗

CSRI Summer Proceedings 2020

The Computer Science Research Institute (CSRI) brings university faculty and students to Sandia for focused collaborative research on Department of Energy (DOE) computer and computational science problems. The institute provides an opportunity for university researchers to learn about problems in computer and computational science at DOE laboratories. Participants conduct leading-edge research, interact with scientists and engineers at the laboratories, and help transfer results of their research to programs at the labs. Some specific CSRI research interest areas are: scalable solvers, optimization, adaptivity and mesh refinement, graph-based, discrete, and combinatorial algorithms, uncertainty estimation, mesh generation, dynamic load-balancing, virus and other malicious-code defense, visualization, scalable cluster computers, data-intensive computing, environments for scalable computing, parallel input/output, advanced architectures, and theoretical computer science. The CSRI Summer Program is organized by CSRI and typically includes the organization of a weekly seminar series and the publication of a summer proceedings. In 2020, the CSRI summer program was executed completely virtually; all student interns worked from home, due to the COVID-19 pandemic.

97 MATHEMATICS AND COMPUTING↗

High order parallel numerical schemes for solving incompressible flows

The use of parallel computers for numerically solving flow fields has gained much importance in recent years. This paper introduces a new high order numerical scheme for computational fluid dynamics (CFD) specifically designed for parallel computational environments. A distributed MIMD system gives the flexibility of treating different elements of the governing equations with totally different numerical schemes in different regions of the flow field. The parallel decomposition of the governing operator to be solved is the primary parallel split. The primary parallel split was studied using a hypercube like architecture having clusters of shared memory processors at each node. The approach is demonstrated using examples of simple steady state incompressible flows. Future studies should investigate the secondary split because, depending on the numerical scheme that each of the processors applies and the nature of the flow in the specific subdomain, it may be possible for a processor to seek better, or higher order, schemes for its particular subcase.

Lin, Avi↗

Ensuring correct rollback recovery in distributed shared memory systems

Distributed shared memory (DSM) implemented on a cluster of workstations is an increasingly attractive platform for executing parallel scientific applications. Checkpointing and rollback techniques can be used in such a system to allow the computation to progress in spite of the temporary failure of one or more processing nodes. This paper presents the design of an independent checkpointing method for DSM that takes advantage of DSM's specific properties to reduce error-free and rollback overhead. The scheme reduces the dependencies that need to be considered for correct rollback to those resulting from transfers of pages. Furthermore, in-transit messages can be recovered without the use of logging. We extend the scheme to a DSM implementation using lazy release consistency, where the frequency of dependencies is further reduced.

Janssens, Bob↗

Parallel Implementation of the Discontinuous Galerkin Method

This paper describes a parallel implementation of the discontinuous Galerkin method. Discontinuous Galerkin is a spatially compact method that retains its accuracy and robustness on non-smooth unstructured grids and is well suited for time dependent simulations. Several parallelization approaches are studied and evaluated. The most natural and symmetric of the approaches has been implemented in all object-oriented code used to simulate aeroacoustic scattering. The parallel implementation is MPI-based and has been tested on various parallel platforms such as the SGI Origin, IBM SP2, and clusters of SGI and Sun workstations. The scalability results presented for the SGI Origin show slightly superlinear speedup on a fixed-size problem due to cache effects.

Baggag, Abdalkader↗

The Hydrogen Abstraction from A Diamond(111) Surface in A Uniform Electric Field

Bond breaking in a strong electric field is shown to arise from a crossing of the ionic and covalent asymptotes. The specific example of hydrogen abstraction from a diamond(111) surface is studied using a cluster model. The addition of nearby atoms in both the parallel and perpendicular direction to the electric field are found to have an effect. It is also shown that the barrier is not only related to the position of the ionic and covalent asymptotes.

Ricca, Alessandra↗

A Parallel Particle Swarm Optimization Algorithm Accelerated by Asynchronous Evaluations

A parallel Particle Swarm Optimization (PSO) algorithm is presented. Particle swarm optimization is a fairly recent addition to the family of non-gradient based, probabilistic search algorithms that is based on a simplified social model and is closely tied to swarming theory. Although PSO algorithms present several attractive properties to the designer, they are plagued by high computational cost as measured by elapsed time. One approach to reduce the elapsed time is to make use of coarse-grained parallelization to evaluate the design points. Previous parallel PSO algorithms were mostly implemented in a synchronous manner, where all design points within a design iteration are evaluated before the next iteration is started. This approach leads to poor parallel speedup in cases where a heterogeneous parallel environment is used and/or where the analysis time depends on the design point being analyzed. This paper introduces an asynchronous parallel PSO algorithm that greatly improves the parallel e ciency. The asynchronous algorithm is benchmarked on a cluster assembled of Apple Macintosh G5 desktop computers, using the multi-disciplinary optimization of a typical transport aircraft wing as an example.

Venter, Gerhard↗

The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z = 4 to 8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Hubble Frontier Field deep parallel observations near the Abell 2744 and MACS J0416.1- 2403 clusters. The combination of these surveys provides an effective volume of 0.6-1.2 ×10(exp 6) Mpc(exp 3) over this epoch, allowing us to perform a robust search for bright (M(sub UV) less than −21) and faint (M(sub UV) = −18) galaxies. We select galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5 less than z less than 8.5, with more than 1000 galaxies at z of approximately 6 - 8. We measure both a stepwise luminosity function for galaxies in our redshift samples, as well as a Schechter function, using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our UV luminosity functions agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z of greater than or equal to 6. Our bestfit value of the characteristic magnitude M* is consistent with −21 at z of greater than or equal to 5, different than that inferred based on previous trends at lower redshift. At z = 8, a single power-law provides an equally good fit to the UV luminosity function, while at z = 6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare our luminosity functions to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift, though a decreasing impact of feedback may also be possible. We measure the evolution of the cosmic star-formation rate (SFR) density by integrating our observed luminosity functions to M(sub UV) = −17, correcting for dust attenuation, and find that the SFR density declines proportionally to (1 + z)((exp −4.3)(+/-)(0.5)) at z greater than 4, consistent with observations at z greater than or equal to 9. Our observed luminosity functions are consistent with a reionization history that starts at redshift of approximately greater than 10, completes at z greater than 6, and reaches a midpoint (x(sub HII) = 0.5) at 6.7 less than z less than 9.4. Finally, using a constant cumulative number density selection and an empirically derived rising star-formation history, our observations predict that the abundance of bright z = 9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z similar to 10 galaxies.

Galaxy↗

Motions in Prominence Barbs Observed on the Solar Limb

We analyze and discuss an example of prominence barbs observed on the limb on 2016 January 7 by the Hinode/Solar Optical Telescope in Ca II and H alpha, the Interface Region Imaging Spectrograph, with slit jaw images and Mg II spectral data, and the Solar Dynamics Observatory's Atmospheric Imaging Assembly. In the recent literature there has been a debate concerning whether these features, sometimes referred to as "tornadoes," are rotating. Our data analysis provides no evidence for systematic rotation in the barbs. We do find line-of-sight motions in the barbs that vary with location and time. We also discuss observations of features moving along the barbs. These moving features are elongated parallel to the solar limb and tend to come in clusters of features moving along the same or similar paths in the plane of the sky during a period of 10 minutes to an hour, moving toward or away from the limb. The motion may have a component along the line of sight as well. The spectral data indicate that the features are Doppler shifted. We discuss possible explanations for these features.

Kucera, T. A.↗

NETRA: A parallel architecture for integrated vision systems. 1: Architecture and organization

Computer vision is regarded as one of the most complex and computationally intensive problems. An integrated vision system (IVS) is considered to be a system that uses vision algorithms from all levels of processing for a high level application (such as object recognition). A model of computation is presented for parallel processing for an IVS. Using the model, desired features and capabilities of a parallel architecture suitable for IVSs are derived. Then a multiprocessor architecture (called NETRA) is presented. This architecture is highly flexible without the use of complex interconnection schemes. The topology of NETRA is recursively defined and hence is easily scalable from small to large systems. Homogeneity of NETRA permits fault tolerance and graceful degradation under faults. It is a recursively defined tree-type hierarchical architecture where each of the leaf nodes consists of a cluster of processors connected with a programmable crossbar with selective broadcast capability to provide for desired flexibility. A qualitative evaluation of NETRA is presented. Then general schemes are described to map parallel algorithms onto NETRA. Algorithms are classified according to their communication requirements for parallel processing. An extensive analysis of inter-cluster communication strategies in NETRA is presented, and parameters affecting performance of parallel algorithms when mapped on NETRA are discussed. Finally, a methodology to evaluate performance of algorithms on NETRA is described.

Choudhary, Alok N.↗

Discovery of Extended Tidal Tails around the Globular Cluster Palomar 13

We use photometry from the DECam Legacy Survey to detect candidate tidal tails extending ~5° on either side of the Palomar 13 globular cluster. The tails are aligned with the proper motion of Palomar 13 and are consistent with its old, metal-poor stellar population. We identify three RR Lyrae stars (RRLs) that are plausibly associated with the tails, in addition to four previously known in the cluster. From these RRLs, we find that the mean distance to the cluster and tails is 23.6 ± 0.2 kpc and estimate the total (initial) luminosity of the cluster to be , consistent with previous claims that its initial luminosity was higher than its current luminosity. Combined with previously determined proper motion and radial velocity measurements of the cluster, we find that Palomar 13 is on a highly eccentric orbit (e ~ 0.8) with a pericenter of ~9 kpc and an apocenter of ~69 kpc, and a recent pericentric passage of the cluster ~75 Myr ago. We note a prominent linear structure in the interstellar dust map that runs parallel to the candidate tidal features, but conclude that reddening due to dust is unlikely to account for the structure that we observe. If confirmed, the Palomar 13 stellar stream would be one of very few streams with a known progenitor system, making it uniquely powerful for studying the disruption of globular clusters, the formation of the stellar halo, and the distribution of matter within our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Chandra Follow-up of the SDSS DR8 Redmapper Catalog Using the MATCha Pipeline

In order to place constraints on cosmology through optical surveys of galaxy clusters, one must first understand the properties of those clusters. To this end, we introduce the Mass Analysis Tool for Chandra (MATCha), a pipeline that uses a parallellized algorithm to analyze archival Chandra data. MATCha simultaneously calculates X-ray temperatures and luminosities and performs centering measurements for hundreds of potential galaxy clusters using archival X-ray exposures. We run MATCha on the redMaPPer SDSS DR8 cluster catalog and use MATCha's output X-ray temperatures and luminosities to analyze the galaxy cluster temperature–richness, luminosity–richness, luminosity–temperature, and temperature–luminosity scaling relations. We detect 447 clusters and determine 246 r 2500 temperatures across all redshifts. Within 0.1 < z < 0.35, we find that r 2500 T X scales with optical richness (λ) as $\mathrm{ln}\left(\tfrac{{k}_{B}{T}_{{\rm{X}}}}{1.0\,\mathrm{keV}}\right)=(0.52\pm 0.05)\mathrm{ln}\left(\tfrac{\lambda }{70}\right)+(1.85\pm 0.03)$ with an intrinsic scatter of $0.27\pm 0.02$ ($1\sigma $). We investigate the distribution of offsets between the X-ray center and redMaPPer center within 0.1 < z < 0.35, finding that 68%.3 ± 6.5% of clusters are well-centered. However, we find a broad tail of large offsets in this distribution, and we explore some of the causes of redMaPPer miscentering.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring the Hubble constant and our Virgo-infall velocity independently

A sample of spiral galaxies with B(sub T) less than 14.5 located in two local volumes, one in the direction of, but behind, the Virgo Cluster (behind-Virgo volume (BV)) and the other in the opposite direction (anti-Virgo volume (AV)), were used via a Tully-Fisher (TF) relation to derive the following two parameters: H(sub AB), the mean Hubble ratio between AV and BV, and delta v(sub parallel), the peculiar velocity of the Local Group in the direction of the Virgo Cluster (VC) with respect to a uniformly expanding reference system defined by our AV and BV sub-samples. The two sampled volumes, separated by a velocity interval of 5600 km/s, form an antipodal pair. This particular geometry not only allows us to derive the two parameters independently but also reduces the dynamical effect of the Local Supercluster on H(sub AB) without increasing the Malmquist bias. By limiting our sample to spiral galaxies having large velocity widths W(sub R), we effectively reduce the TF scatter and Malmquist bias in our sample. The TF zero point and dispersion were then determined by further correcting for the small residual Malmquist bias. An additional sample of fainter galaxies was used to test for a non-Gaussian tail to the TF disperison. We found no evidence for such a tail and formally give an upper limit of about 18% for the fractional contribution of an unseen tail. The average intrinsic TF dispersion for the dominant Gaussian component is sigma(sub TF)(sup 0) approximately 0.33 mag for W(sub R) approximately equal to or greater than 180 km/s. Our numerical results are delta v(sub parallel) approximately equals 414 +/- 82 km/s and H(sub AB) approximately equals (84.0 +/- 2.4)(1 + epsilon) km/s Mpc, where (1 + epsilon) accounts for any systematic error between the calibrators and the sample galaxies. Various dynamical models were tested to explore the effect on H(sub AB) of the uncertainties in the local velocity field. Constrained by our observed delta v(sub parallel) as well as other observational quantities, we found that the rms deviation from unity of H(sub AB)/H(sub 0) (where H(sub 0) is the Hubble constant for each model) is 5%, making H(sub AB) a good indicator for H(sub 0). Taking this variation as an additional error, our formal estimate for the Hubble constant is H(sub 0) approximately equals (84 +/- 5)(1 + epsilon) km/s Mpc.

Lu, Nanyao Y.↗

Parallel computational fluid dynamics - Implementations and results

The present volume on parallel CFD discusses implementations on parallel machines, numerical algorithms for parallel CFD, and performance evaluation and computer science issues. Attention is given to a parallel algorithm for compressible flows through rotor-stator combinations, a massively parallel Euler solver for unstructured grids, a fast scheme to analyze 3D disk airflow on a parallel computer, and a block implicit multigrid solution of the Euler equations. Topics addressed include a 3D ADI algorithm on distributed memory multiprocessors, clustered element-by-element computations for fluid flow, hypercube FFT and the Fourier pseudospectral method, and an investigation of parallel iterative algorithms for CFD. Also discussed are fluid dynamics using interface methods on parallel processors, sorting for particle flow simulation on the connection machine, a large grain mapping method, and efforts toward a Teraflops capability for CFD.

Simon, Horst D.↗

Impact of Magnetic Draping, Convection, and Field Line Tying on Magnetopause Reconnection Under Northward IMF

We simulate a northward IMF cusp reconnection event at the magnetopause using the OpenGGCM resistive MHD code. The ACE input data, solar wind parameters, and dipole tilt belong to a 2002 reconnection event observed by IMAGE and Cluster. Based on a fully three-dimensional skeleton separators, nulls, and parallel electric fields, we show magnetic draping, convection, ionospheric field line tying play a role in producing a series of locally reconnecting nulls with flux ropes. The flux ropes in the cusp along the global separator line of symmetry. In 2D projection, the flux ropes the appearance of a tearing mode with a series of 'x's' and 'o's' but bearing a kind of 'guide field' that exists only within the magnetopause. The reconnecting field lines in the string of ropes involve IMF and both open and closed Earth magnetic field lines. The observed magnetic geometry reproduces the findings of a superposed epoch impact parameter study derived from the Cluster magnetometer data for the same event. The observed geometry has repercussions for spacecraft observations of cusp reconnection and for the imposed boundary conditions reconnection simulations.

Wendel, Deirdre E.↗

Measuring the homogeneity of the universe using polarization drift

We propose a method to probe the homogeneity of a general universe, without assuming symmetry. We show that isotropy can be tested at remote locations on the past lightcone by comparing the line-of-sight and transverse expansion rates, using the time dependence of the polarization of Cosmic Microwave Background photons that have been inverse-Compton scattered by the hot gas in massive clusters of galaxies. This probes a combination of remote transverse and parallel components of the expansion rate of the metric, and we may use radial baryon acoustic oscillations or cosmic clocks to measure the parallel expansion rate. Thus we can test remote isotropy, which is a key requirement of a homogeneous universe. We provide explicit formulas that connect observables and properties of the metric.

79 ASTRONOMY AND ASTROPHYSICS↗

DAAP (Data Analytics Application Profiling)

The Data Analytics Application Profiling (DAAP) utility is a software program comprised of a library and scripts to provide basic instrumenting of parallel scientific software applications while they are running on an HPC cluster, and transmit this data over a secure connection as individual records to aggregators. When used in conjunction with additional open source software components (Telegraf, FluentD, and RabbitMQ), the data collected by DAAP can be transmitted off-cluster to an analytics system of the user’s choice for further analysis.

Shereda, Charles↗