Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel clustering”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Large-scale Nanostructure Simulations from X-ray Scattering Data On Graphics Processor Clusters

X-ray scattering is a valuable tool for measuring the structural properties of materialsused in the design and fabrication of energy-relevant nanodevices (e.g., photovoltaic, energy storage, battery, fuel, and carbon capture andsequestration devices) that are key to the reduction of carbon emissions. Although today's ultra-fast X-ray scattering detectors can provide tremendousinformation on the structural properties of materials, a primary challenge remains in the analyses of the resulting data. We are developing novelhigh-performance computing algorithms, codes, and software tools for the analyses of X-ray scattering data. In this paper we describe two such HPCalgorithm advances. Firstly, we have implemented a flexible and highly efficient Grazing Incidence Small Angle Scattering (GISAXS) simulation code based on theDistorted Wave Born Approximation (DWBA) theory with C++/CUDA/MPI on a cluster of GPUs. Our code can compute the scattered light intensity from any givensample in all directions of space; thus allowing full construction of the GISAXS pattern. Preliminary tests on a single GPU show speedups over 125x compared tothe sequential code, and almost linear speedup when executing across a GPU cluster with 42 nodes, resulting in an additional 40x speedup compared to usingone GPU node. Secondly, for the structural fitting problems in inverse modeling, we have implemented a Reverse Monte Carlo simulation algorithm with C++/CUDAusing one GPU. Since there are large numbers of parameters for fitting in the in X-ray scattering simulation model, the earlier single CPU code required weeks ofruntime. Deploying the AccelerEyes Jacket/Matlab wrapper to use GPU gave around 100x speedup over the pure CPU code. Our further C++/CUDA optimization deliveredan additional 9x speedup.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Shergottite Lead Isotope Signature in Chassigny and the Nakhlites

The nakhlites/chassignites and the shergottites represent two differing suites of basaltic martian meteorites. The shergottites have ages less than or equal to 0.6 Ga and a large range of initial Sr-/Sr-86 and epsilon (Nd-143) ratios. Conversely, the nakhlites and chassignites cluster at 1.3-1.4 Ga and have a limited range of initial Sr-87/Sr-86 and epsilon (Nd-143). More importantly, the shergottites have epsilon (W-182) less than 1, whereas the nakhlites and chassignites have epsilon (W-182) approximately 3. This latter observation precludes the extraction of both meteorite groups from a single source region. However, recent Pb isotopic analyses indicate that there may have been interaction between shergottite and nakhlite/chassignite Pb reservoirs.Pb Analyses of Chassigny: Two different studies haveinvestigated 207Pb/204Pb vs. 206Pb/204Pb in Chassigny: (i)TIMS bulk-rock analyses of successive leaches and theirresidue [3]; and (ii) SIMS analysis of individual minerals[4]. The bulk-rock analyses fall along a regression of SIMSplagioclase analyses that define an errorchron that is olderthan the Solar System (4.61±0.1 Ga); i.e., these define amixing line between Chassigny’s principal Pb isotopic components(Fig. 1). Augites and olivines in Chassingy (notshown) also fall along or near the plagioclase regression [4].This agreement indicates that the whole-rock leachateslikely measure indigenous, martian Pb, not terrestrial contamination[5]. SIMS analyses of K-spars and sulfides definea separate, sub-parallel trend having higher 207Pb/206Pbvalues ([4]; Fig. 1). The good agreement between the bulkrockanalyses and the SIMS analyses of plagioclases alsoindicates that the Pb in the K-spars and sulfides cannot be amajor component of Chassigny.The depleted reservoir sampled by Chassigny plagioclaseis not the same as the solar system initial (PAT) andrequires a multi-stage origin. Here we show a two-stagemodel (Fig. 1) with a 238U/204Pb (μ) of 0.5 for 4.5-2.4 Gaand a μ of 7 for 2.4-1.4 Ga. This is not a unique model butdoes produce a Pb composition that falls on the plagioclaseregression at 1.4 Ga, the approximate igneous age of Chassigny [1]. It should be noted that low-μ single-stage modelsare not capable of producing sufficiently radiogenic 206Pb/204Pb at 1.4 Ga.Relation to Shergottites: The Chassigny K-spars and sulfides fall along a second mixing line defined by leachesand residues of depleted and intermediate shergottites [6]. This mixing line falls above the plagioclase regression.Therefore, we also interpret the radiogenic component of this mixing line to represent indigenous martian Pb. It ispossible that the depleted and intermediate shergottites and the Chassigny plagioclases sample radiogenic Pb from thethe same source, i.e., the mixing lines may intersect at high 206Pb/204Pb.Both K-spar and sulfide are late-stage phases. At the time of their crystallization, the Chassigny system appearsto have remained open to a depleted shergottite Pb reservoir. The depleted component of the shergottite mixing linecan be generated by a single-stage evolution from PAT (4.5 to 1.4 Ga) in a reservoir having a μ ~2. A similar modelfor the most depleted shergottites is also possible: μ = 1.5 for 4.5 to 0.3 Ga.Nakhlites: Nakhlite analyses plot between the shergottite and Chassigny plagioclase regressions [3]. So again,members of the nakhlite/chassignite suite show affinities to shergottite Pb.

Jones, J. H.↗
NASA NTRS2015

A Data Exploration Tool for Large Sets of Spectra

We present an exploration tool for very large spectrum data sets such as the SDSS (Sloan Digital Sky Survey), LAMOST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope), and 4MOST (4-meter Multi-Object Spectroscopic Telescope) data sets. The tool works in two stages: the first uses batch processing and the second runs interactively. The latter employs the NASA hyperwall, a configuration of 128 workstation displays (8 by 16 array) controlled by a parallelized software suite running on NASA's Pleiades supercomputer. The stellar subset of the Sloan Digital Sky Survey, DR10, was chosen to show how the our tool may be used. In stage one, SDSS files for 569,740 stars are processed through our data pipeline. The pipeline fits each spectrum using an iterative continuum algorithm, distinguishing emission from absorption and handling molecular absorption bands correctly. It then measures 1659 discrete atomic and molecular spectral features that were carefully preselected based on their likelihood of being visible at some spectral type. The depths relative to the local continuum at each feature wavelength are determined for each spectrum: these depths, the local S/N (signal to noise ratio) level, and DR10-supplied variables such as magnitudes, colors, positions, and radial velocities are the basic measured quantities used on the hyperwall. In stage two, each hyperwall panel is used to display a 2-D scatter plot showing the depth of feature A vs the depth of feature B for all of the stars. A and B change from panel to panel. The relationships between the various (A,B) strengths and any distinctive clustering are immediately apparent when examining and inter-comparing the different panels on the hyperwall. The interactive software allows the user to select the stars in any interesting region of any 2-D plot on the hyperwall, immediately rendering the same stars on all the other 2-D plots in a unique color. The process may be repeated multiple times, each selection displaying a distinctive color on all the plots. At any time, the spectra of the selected stars may be examined in detail on a connected workstation display. We illustrate how our approach allows us to quickly isolate and examine such interesting stellar subsets as EMP (Extremely Metal‐Poor) stars, CV (Cataclymic Variable) stars and C (Carbon)-rich stars.

Data Exploration↗
DOE OSTI2019

iez: Resource Contention Aware Load Balancing for Large-Scale Parallel File Systems

Parallel I/O performance is crucial to sustaining scientific applications on large-scale High-Performance Computing (HPC) systems. However, I/O load imbalance in the underlying distributed and shared storage systems can significantly reduce overall application performance. There are two conflicting challenges to mitigate this load imbalance: (i) optimizing systemwide data placement to maximize the bandwidth advantages of distributed storage servers, i.e., allocating I/O resources efficiently across applications and job runs; and (ii) optimizing client-centric data movement to minimize I/O load request latency between clients and servers, i.e., allocating I/O resources efficiently in service to a single application and job run. Moreover, existing approaches that require application changes limit wide-spread adoption in commercial or proprietary deployments. We propose iez, an “end-to-end control plane” where clients transparently and adaptively write to a set of selected I/O servers to achieve balanced data placement. Our control plane leverages realtime load information for distributed storage server global data placement while our design model leverages trace-based optimization techniques to minimize I/O load request latency between clients and servers. We evaluate our proposed system on an experimental cluster for two common use cases: synthetic I/O benchmark IOR for large sequential writes and a scientific application I/O kernel, HACC-I/O. Results show read and write performance improvements of up to 34% and 32%, respectively, compared to the state of the art.

Wadhwa, Bharti↗
DOE OSTI2015

Lattice cluster theory for dense, thin polymer films

While the application of the lattice cluster theory (LCT) to study the miscibility of polymer blends has greatly expanded our understanding of the monomer scale molecular details influencing miscibility, the corresponding theory for inhomogeneous systems has not yet emerged because of considerable technical difficulties and much greater complexity. Here, we present a general formulation enabling the extension of the LCT to describe the thermodynamic properties of dense, thin polymer films using a high dimension, high temperature expansion. Whereas the leading order of the LCT for bulk polymer systems is essentially simple Flory-Huggins theory, the highly non-trivial leading order inhomogeneous LCT (ILCT) for a film with L layers already involves the numerical solution of 3(L − 1) coupled, highly nonlinear equations for the various density profiles in the film. The new theory incorporates the essential “transport” constraints of Helfand and focuses on the strict imposition of excluded volume constraints, appropriate to dense polymer systems, rather than the maintenance of chain connectivity as appropriate for lower densities and as implemented in self-consistent theories of polymer adsorption at interfaces. The ILCT is illustrated by presenting examples of the computed profiles of the density, the parallel and perpendicular bonds, and the chain ends for free standing and supported films as a function of average film density, chain length, temperature, interaction with support, and chain stiffness. The results generally agree with expected general trends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

MPIRUN: A Portable Loader for Multidisciplinary and Multi-Zonal Applications

Multidisciplinary and multi-zonal applications are an important class of applications in the area of Computational Aerosciences. In these codes, two or more distinct parallel programs or copies of a single program are utilized to model a single problem. To support such applications, it is common to use a programming model where a program is divided into several single program multiple data stream (SPMD) applications, each of which solves the equations for a single physical discipline or grid zone. These SPMD applications are then bound together to form a single multidisciplinary or multi-zonal program in which the constituent parts communicate via point-to-point message passing routines. One method for implementing the message passing portion of these codes is with the new Message Passing Interface (MPI) standard. Unfortunately, this standard only specifies the message passing portion of an application, but does not specify any portable mechanisms for loading an application. MPIRUN was developed to provide a portable means for loading MPI programs, and was specifically targeted at multidisciplinary and multi-zonal applications. Programs using MPIRUN for loading and MPI for message passing are then portable between all machines supported by MPIRUN. MPIRUN is currently implemented for the Intel iPSC/860, TMC CM5, IBM SP-1 and SP-2, Intel Paragon, and workstation clusters. Further, MPIRUN is designed to be simple enough to port easily to any system supporting MPI.

Fineberg, Samuel A.↗
NASA NTRS2006

Forecasting of Storm-Surge Floods Using ADCIRC and Optimized DEMs

Increasing the accuracy of storm-surge flood forecasts is essential for improving preparedness for hurricanes and other severe storms and, in particular, for optimizing evacuation scenarios. An interactive database, developed by WorldWinds, Inc., contains atlases of storm-surge flood levels for the Louisiana/Mississippi gulf coast region. These atlases were developed to improve forecasting of flooding along the coastline and estuaries and in adjacent inland areas. Storm-surge heights depend on a complex interaction of several factors, including: storm size, central minimum pressure, forward speed of motion, bottom topography near the point of landfall, astronomical tides, and, most importantly, maximum wind speed. The information in the atlases was generated in over 100 computational simulations, partly by use of a parallel-processing version of the ADvanced CIRCulation (ADCIRC) model. ADCIRC is a nonlinear computational model of hydrodynamics, developed by the U.S. Army Corps of Engineers and the US Navy, as a family of two- and three-dimensional finite-element-based codes. It affords a capability for simulating tidal circulation and storm-surge propagation over very large computational domains, while simultaneously providing high-resolution output in areas of complex shoreline and bathymetry. The ADCIRC finite-element grid for this project covered the Gulf of Mexico and contiguous basins, extending into the deep Atlantic Ocean with progressively higher resolution approaching the study area. The advantage of using ADCIRC over other storm-surge models, such as SLOSH, is that input conditions can include all or part of wind stress, tides, wave stress, and river discharge, which serve to make the model output more accurate. To keep the computational load manageable, this work was conducted using only the wind stress, calculated by using historical data from Hurricane Camille, as the input condition for the model. Hurricane storm-surge simulations were performed on an eight-node Linux computer cluster. Each node contained dual 2-GHz processors, 2GB of memory, and a 40GB hard drive. The digital elevation model (DEM) for this region was specified using a combination of Navy data (over water), NOAA data (for the coastline), and optimized Interferometric Synthetic Aperture Radar data (over land). This high-resolution topographical data of the Mississippi coastal region provided the ADCIRC model with improved input with which to calculate improved storm-surge forecasts.

Valenti, Elizabeth↗
DOE OSTI2014

Final Report for Project DE-FC02-06ER25755 [Pmodels2]

In this report, we describe the research accomplished by the OSU team under the Pmodels2 project. The team has worked on various angles: designing high performance MPI implementations on modern networking technologies (Mellanox InfiniBand (including the new ConnectX2 architecture and Quad Data Rate), QLogic InfiniPath, the emerging 10GigE/iWARP and RDMA over Converged Enhanced Ethernet (RoCE) and Obsidian IB-WAN), studying MPI scalability issues for multi-thousand node clusters using XRC transport, scalable job start-up, dynamic process management support, efficient one-sided communication, protocol offloading and designing scalable collective communication libraries for emerging multi-core architectures. New designs conforming to the Argonne’s Nemesis interface have also been carried out. All of these above solutions have been integrated into the open-source MVAPICH/MVAPICH2 software. This software is currently being used by more than 2,100 organizations worldwide (in 71 countries). As of January ’14, more than 200,000 downloads have taken place from the OSU Web site. In addition, many InfiniBand vendors, server vendors, system integrators and Linux distributors have been incorporating MVAPICH/MVAPICH2 into their software stacks and distributing it. Several InfiniBand systems using MVAPICH/MVAPICH2 have obtained positions in the TOP500 ranking of supercomputers in the world. The latest November ’13 ranking include the following systems: 7th ranked Stampede system at TACC with 462,462 cores; 11th ranked Tsubame 2.5 system at Tokyo Institute of Technology with 74,358 cores; 16th ranked Pleiades system at NASA with 81,920 cores; Work on PGAS models has proceeded on multiple directions. The Scioto framework, which supports task-parallelism in one-sided and global-view parallel programming, has been extended to allow multi-processor tasks that are executed by processor groups. A quantum Monte Carlo application is being ported onto the extended Scioto framework. A public release of Global Trees (GT) has been made, along with the Global Chunks (GC) framework on which GT is built. The Global Chunks (GC) layer is also being used as the basis for the development of a higher level Global Graphs (GG) layer. The Global Graphs (GG) system will provide a global address space view of distributed graph data structures on distributed memory systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

The artificial retina processor for track reconstruction at the LHC crossing rate

We present results of an R&D study for a specialized processor capable of precisely reconstructing, in pixel detectors, hundreds of charged-particle tracks from high-energy collisions at 40 MHz rate. We apply a highly parallel pattern-recognition algorithm, inspired by studies of the processing of visual images by the brain as it happens in nature, and describe in detail an efficient hardware implementation in high-speed, high-bandwidth FPGA devices. This is the first detailed demonstration of reconstruction of offline-quality tracks at 40 MHz and makes the device suitable for processing Large Hadron Collider events at the full crossing frequency.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

A REVISED PARALLEL-SEQUENCE MORPHOLOGICAL CLASSIFICATION OF GALAXIES: STRUCTURE AND FORMATION OF S0 AND SPHEROIDAL GALAXIES

We update van den Bergh's parallel-sequence galaxy classification in which S0 galaxies form a sequence S0a-S0b-S0c that parallels the sequence Sa-Sb-Sc of spiral galaxies. The ratio B/T of bulge-to-total light defines the position of a galaxy in this tuning-fork diagram. Our classification makes one major improvement. We extend the S0a-S0b-S0c sequence to spheroidal ('Sph') galaxies that are positioned in parallel to irregular galaxies in a similarly extended Sa-Sb-Sc-Im sequence. This provides a natural 'home' for spheroidals, which previously were omitted from galaxy classification schemes or inappropriately combined with ellipticals. To motivate our juxtaposition of Sph and Im galaxies, we present photometry and bulge-disk decompositions of four rare, late-type S0s that bridge the gap between the more common S0b and Sph galaxies. NGC 4762 is an edge-on SB0bc galaxy with a very small classical-bulge-to-total ratio of B/T = 0.13 {+-} 0.02. NGC 4452 is an edge-on SB0 galaxy with an even tinier pseudobulge-to-total ratio of PB/T = 0.017 {+-} 0.004. It is therefore an SB0c. VCC 2048, whose published classification is S0, contains an edge-on disk, but its 'bulge' plots in the structural parameter sequence of spheroidals. It is therefore a disky Sph. And NGC 4638 is similarly a 'missing link' between S0s and Sphs-it has a tiny bulge and an edge-on disk embedded in an Sph halo. In the Appendix, we present photometry and bulge-disk decompositions of all Hubble Space Telescope Advanced Camera for Surveys Virgo Cluster Survey S0s that do not have published decompositions. We use these data to update the structural parameter correlations of Sph, S+Im, and E galaxies. We show that Sph galaxies of increasing luminosity form a continuous sequence with the disks (but not bulges) of S0c-S0b-S0a galaxies. Remarkably, the Sph-S0-disk sequence is almost identical to that of Im galaxies and spiral galaxy disks. We review published observations for galaxy transformation processes, particularly ram-pressure stripping of cold gas. We suggest that Sph galaxies are transformed, 'red and dead' Scd-Im galaxies in the same way that many S0 galaxies are transformed, red and dead Sa-Sc spiral galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Discovery of a Ferromagnetic Nickel Chalcogenide Nanocluster Ni 3 S 3 H(PEt 3 ) 5

Atomically precise ligated nanoclusters (NC) are promising cluster-based materials with novel molecular architectures and tunable magnetic properties. Herein, the synthesis and characterization of a nickel sulfide NC Ni 3 S 3 H(PEt 3 ) 5 (PEt 3 = triethylphosphine) with distinct magnetic properties are reported. Magnetization measurements reveal its magnetic moment of 1.5 µ B in the solid phase, consistent with the existence of one unpaired electron predicted by density functional theory (DFT) calculations. Additionally, experimental measurements indicate the presence of ferromagnetic ordering within each Ni 3 S 3 H(PEt 3 ) 5 NC and strong coercivity at temperatures below 20 K. Ion mobility-mass spectrometry is employed in conjunction with DFT calculations and collision cross-section simulations to investigate the structure of the isolated Ni 3 S 3 H(PEt 3 ) 5 . Theoretical studies show that [Ni 3 S 3 H(PEt 3 ) 5 ] + has a planar Ni 3 S 3 core where three Ni atoms are arranged in a triangle with three bridging S atoms residing in the same plane. This structure is preserved in both solution and solid phases, which is confirmed by spectroscopic studies of Ni 3 S 3 H(PEt 3 ) 5 . Additionally, DFT calculations indicate that all spins at the Ni sites are aligned parallel, confirming the presence of ferromagnetic coupling. Overall, this study provides key insights into the structure and magnetic properties of Ni 3 S 3 H(PEt 3 ) 5 , which will facilitate the design of new NC-based magnetic materials.

Nickel sulfide nanocluster↗
DOE OSTI2012

Combining Partial Redundancy and Checkpointing for HPC

Today's largest High Performance Computing (HPC) systems exceed one Petaflops (10^15 floating point operations per second) and exascale systems are projected within seven years. But reliability is becoming one of the major challenges faced by exascale computing. With billion-core parallelism, the mean time to failure is projected to be in the range of minutes or hours instead of days. Failures are becoming the norm rather than the exception during execution of HPC applications. Current fault tolerance techniques in HPC focus on reactive ways to mitigate faults, namely via checkpoint and restart (C/R). Apart from storage overheads, C/R-based fault recovery comes at an additional cost in terms of application performance because normal execution is disrupted when checkpoints are taken. Studies have shown that applications running at a large scale spend more than 50\% of their total time saving checkpoints, restarting and redoing lost work. Redundancy is another fault tolerance technique, which employs redundant processes performing the same task. If a process fails, a replica of it can take over its execution. Thus, redundant copies can decrease the overall failure rate. The downside of redundancy is that extra resources are required and there is an additional overhead on communication and synchronization. This work contributes a model and analyzes the benefit of C/R in coordination with redundancy at different degrees to minimize the total wallclock time and resources utilization of HPC applications. We further conduct experiments with an implementation of redundancy within the MPI layer on a cluster. Our experimental results confirm the benefit of dual and triple redundancy - but not for partial redundancy - and show a close fit to the model. At 80,000 processes, dual redundancy requires twice the number of processing resources for an application but allows two jobs of 128 hours wallclock time to finish within the time of just one job without redundancy. For narrow ranges of processor counts, partial redundancy results in the lowest time. Once the count exceeds 770, 000, triple redundancy has the lowest overall cost. Thus, redundancy allows one to trade-off additional resource requirements against wallclock time, which provides a tuning knob for users to adapt to resource availabilities.

Elliott, James↗
DOE OSTI2020

Bridging the Transuranics with Uranium(IV) Sulfate Aqueous Species and Solid Phases

Isolating isomorphic compounds of tetravalent actinides (i.e., Th IV , U IV , Np IV , and Pu IV ) improve our understanding of the bonding behavior across the series, in addition to their relationship with tetravalent transition metals (Zr and Hf) and lanthanides (Ce). Similarities between these tetravalent metals are particularly illuminated in their hydrolysis and condensation behavior in aqueous systems, leading to polynuclear clusters typified by the hexamer [M IV 6 O 4 (OH) 4 ] 12+ building block. Prior studies have shown the predominance and coexistence of smaller species for Th IV (monomers, dimers, and hexamers) and larger species for U IV , Np IV , and Pu IV (including 38-mers and 70-mers). We show here that aqueous uranium(IV) sulfate also displays behavior similar to that of Th IV (and Zr IV ) in its isolated solidphase and solution speciation. Two single-crystal X-ray structures are described: a dihydroxide-bridged dimer (U 2 ) formulated as U 2 (OH) 2 (SO 4 ) 3 (H 2 O) 4 and a monomer-linked hexamer framework (U-U 6 ) as (U(H 2 O) 3.5 ) 2 U 6 O 4 (OH) 4 (SO 4 ) 10 (H 2 O) 9 . These structures are similar to those previously described for ThIV. Moreover, cocrystallization of monomer and dimer and of dimer and monomer–hexamer phases for both Th IV (prior) and U IV (current) indicates the coexistence of these species in solution. Because it was not possible to effectively study the sulfate-rich solutions via X-ray scattering from which U 2 and U-U 6 crystallized, we provide a parallel solution speciation study in low sulfate conditions, as a function of the pH. Raman spectroscopy, UV–vis spectroscopy, and small-angle X-ray scattering of these show decreasing sulfate binding, increased hydrolysis, increased species size, and increased complexity, with increasing pH. Furthermore, this study describes a bridge across the first half the actinide series, highlighting U IV similarities to Th IV , in addition to the previously known similarities to the transuranic elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

MemXCT: Design, Optimization, Scaling, and Reproducibility of X-Ray Tomography Imaging

Here, this work extends our previous research entitled "MemXCT: Memory-centric X-ray CT Reconstruction with Massive Parallelization" that was originally published at SC19 conference (Hidayetoglu et al., 2019) with reproducibility of the computational imaging performance. X-ray computed tomography (XCT) is regularly used at synchrotron light sources to study the internal morphology of materials at high resolution. However, experimental constraints, such as radiation sensitivity, can result in noisy or undersampled measurements. Further, depending on the resolution, sample size and data acquisition rates, the resulting noisy dataset can be in the order of terabytes. Advanced iterative reconstruction techniques can produce high-quality images from noisy measurements, but their computational requirements have made their use an exception rather than the rule. We propose a novel memory-centric approach that avoids redundant computations at the expense of additional memory complexity. We develop a memory-centric iterative reconstruction system, MemXCT, that uses an optimized SpMV implementation with two-level pseudo-Hilbert ordering and multi-stage input buffering. We evaluate MemXCT on various supercomputer architectures involving KNL and GPU. MemXCT can reconstruct a large (11Kx11K) mouse brain tomogram in 10 seconds using 4096 KNL nodes (256K cores). The results presented in our original article at the SC19 were based on large-scale supercomputing resources. The MemXCT application was selected for the Student Cluster Competition (SCC) Reproducibility Challenge and evaluated on a variety of cloud computing resources by universities around the world in the SC20 conference. We summarize the results of the top-ranked SCC Reproducibility Challenge teams and identify the most pertinent measures for ensuring the reproducibility of our experiments in this article.

47 OTHER INSTRUMENTATION↗
NASA NTRS1996

Non-Linear Cosmological Power Spectra in Real and Redshift Space

We present an expression for the non-linear evolution of the cosmological power spectrum based on Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self-gravitating fields. Some exact solutions are found for power-law initial spectra. We have extended this analysis into red-shift space and found a solution for the non-linear, anisotropic redshift-space power spectrum in the limit of plane-parallel redshift distortions. The quadrupole-to-monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, beta. The point of zero-crossing of the quadrupole, kappa(sub o), is found to obey a simple scaling relation and we calculate this scale in the Zel'dovich approximation. This model is found to be in good agreement with a series of N-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole-to-monopole ratio found in the merged QDOT plus 1.2-Jy-IRAS redshift survey. Using a likelihood technique we have estimated that the distortion parameter is constrained to be beta greater than 0.5 at the 95 percent level. Our results are fairly insensitive to the local primordial spectral slope, but the likelihood analysis suggests n = -2 un the translinear regime. The zero-crossing scale of the quadrupole is k(sub 0) = 0.5 +/- 0.1 h Mpc(exp -1) and from this we infer that the amplitude of clustering is sigma(sub 8) = 0.7 +/- 0.05. We suggest that the success of this model is due to non-linear redshift-space effects arising from infall on to caustic and is not dominated by virialized cluster cores. The latter should start to dominate on scales below the zero-crossing of the quadrupole, where our model breaks down.

Taylor, A. N.↗
NASA NTRS2012

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗
NASA NTRS2016

DeepSAT: A Deep Learning Approach to Tree-Cover Delineation in 1-m NAIP Imagery for the Continental United States

High resolution tree cover classification maps are needed to increase the accuracy of current land ecosystem and climate model outputs. Limited studies are in place that demonstrates the state-of-the-art in deriving very high resolution (VHR) tree cover products. In addition, most methods heavily rely on commercial softwares that are difficult to scale given the region of study (e.g. continents to globe). Complexities in present approaches relate to (a) scalability of the algorithm, (b) large image data processing (compute and memory intensive), (c) computational cost, (d) massively parallel architecture, and (e) machine learning automation. In addition, VHR satellite datasets are of the order of terabytes and features extracted from these datasets are of the order of petabytes. In our present study, we have acquired the National Agriculture Imagery Program (NAIP) dataset for the Continental United States at a spatial resolution of 1-m. This data comes as image tiles (a total of quarter million image scenes with ~60 million pixels) and has a total size of ~65 terabytes for a single acquisition. Features extracted from the entire dataset would amount to ~8-10 petabytes. In our proposed approach, we have implemented a novel semi-automated machine learning algorithm rooted on the principles of "deep learning" to delineate the percentage of tree cover. Using the NASA Earth Exchange (NEX) initiative, we have developed an end-to-end architecture by integrating a segmentation module based on Statistical Region Merging, a classification algorithm using Deep Belief Network and a structured prediction algorithm using Conditional Random Fields to integrate the results from the segmentation and classification modules to create per-pixel class labels. The training process is scaled up using the power of GPUs and the prediction is scaled to quarter million NAIP tiles spanning the whole of Continental United States using the NEX HPC supercomputing cluster. An initial pilot over the state of California spanning a total of 11,095 NAIP tiles covering a total geographical area of 163,696 sq. miles has produced true positive rates of around 88 percent for fragmented forests and 74 percent for urban tree cover areas, with false positive rates lower than 2 percent for both landscapes.

Imagery↗
DOE OSTI2016

Structural basis for promiscuous PAM recognition in type I–E Cascade from E. coli

Clustered regularly interspaced short palindromic repeats (CRISPRs) and the cas (CRISPR-associated) operon form an RNA-based adaptive immune system against foreign genetic elements in prokaryotes. Type I accounts for 95% of CRISPR systems, and has been used to control gene expression and cell fate. During CRISPR RNA (crRNA)-guided interference, Cascade (CRISPR-associated complex for antiviral defence) facilitates the crRNA-guided invasion of double-stranded DNA for complementary base-pairing with the target DNA strand while displacing the non-target strand, forming an R-loop. Cas3, which has nuclease and helicase activities, is subsequently recruited to degrade two DNA strands. A protospacer adjacent motif (PAM) sequence flanking target DNA is crucial for self versus foreign discrimination. Here we present the 2.45 Å crystal structure of Escherichia coli Cascade bound to a foreign double-stranded DNA target. The 5'-ATG PAM is recognized in duplex form, from the minor groove side, by three structural features in the Cascade Cse1 subunit. The promiscuity inherent to minor groove DNA recognition rationalizes the observation that a single Cascade complex can respond to several distinct PAM sequences. Optimal PAM recognition coincides with wedge insertion, initiating directional target DNA strand unwinding to allow segmented base-pairing with crRNA. The non-target strand is guided along a parallel path 25 Å apart, and the R-loop structure is further stabilized by locking this strand behind the Cse2 dimer. Furthermore, these observations provide the structural basis for understanding the PAM-dependent directional R-loop formation process.

59 BASIC BIOLOGICAL SCIENCES↗