Search NASA⌕ Search

SEARCH · Search NASA

Results for “Global Arrays”

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 37 records · Page 2

Mars Meteorolgical Network

Exploring and ultimately establishing a permanent presence on the surface of Mars will necessitate an understanding the weather conditions and the ability to forecast its dynamic behavior. The meteorology of Mars will need to be developed. This abstract puts forth a concept for a Mars Meteorological Network that will be used to investigate the Mars atmosphere behavior, explore the surface environment, and prepare for operational activities. It is proposed that the long term and the dynamic nature of the lower atmosphere and surface of Mars be observed with a distributed global array of simple automated surface nodes. The data would be ingested into the Mars Global Reference Atmospheric Model (Mars-GRAM) and other research tools for analyses to gain a better understanding of the atmospheric conditions on Mars.

Justh, H. L.↗

The Aquarius Mission: Sea Surface Salinity from Space

Aquarius is a new satellite mission concept to study the impact of the global water cycle on the ocean, including the response of the ocean to buoyancy forcing and the subsequent feedback of the ocean on the climate. The measurement objective of Aquarius is sea surface salinity, which reflects the concentration of freshwater at the ocean surface. Salinity affects the dielectric constant of sea water and, consequently, the radiometric emission of the sea surface to space. Rudimentary space observations with an L-band radiometer were first made from Skylab in the mid-70s and numerous aircraft missions of increasing quality and improved technology have been conducted since then. Technology is now available to carry out a global mission, which includes both an accurate L band (1.413 Ghz) radiometer and radar system in space and a global array of in situ observations for calibration and validation, in order to address key NASA Earth Science Enterprise questions about the global cycling of water and the response of the ocean circulation to climate change. The key scientific objectives of Aquarius examine the cycling of water at the ocean's surface, the response of the ocean circulation to buoyancy forcing, and the impact of buoyancy forcing on the ocean's thermal feedback to the climate. Global surface salinity will also improve our ability to model the surface solubility chemistry needed to estimate the air-sea exchange of CO2. In order to meet these science objectives, the NASA Salinity Sea Ice Working Group over the past three years has concluded that the mission measurement goals should be better than 0.2 practical salinity units (psu) accuracy, 100 km resolution, and weekly to revisits. The Aquarius mission proposes to meet these measurement requirements through a real aperture dual-polarized L band radiometer and radar system. This system can achieve the less than 0.1 K radiometric temperature measurement accuracy that is required. A 3 m antenna at approx. 600km altitude in a sun-synchronous orbit and 300 km swath can provide the desired 100 km resolution global coverage every week. Within this decade, it may be possible to combine satellite sea surface salinity measurements with ongoing satellite observations of temperature, surface height, air-sea fluxes; vertical profiles of temperature and salinity from the Argo program; and modern ocean/atmosphere modeling and data assimilation tools, in order to finally address the complex influence of buoyancy on the ocean circulation and climate.

Koblinsky, Chester↗

Versioned distributed arrays for resilience in scientific applications: Global view resilience

Exascale studies project reliability challenges for future high-performance computing (HPC) systems. We propose the Global View Resilience (GVR) system, a library that enables applications to add resilience in a portable, application-controlled fashion using versioned distributed arrays. We describe GVR’s interfaces to distributed arrays, versioning, and cross-layer error recovery. Using several large applications (OpenMC, the preconditioned conjugate gradient solver PCG, ddcMD, and Chombo), we evaluate the programmer effort to add resilience. The required changes are small (<2% LOC), localized, and machine-independent, requiring no software architecture changes. We also measure the overhead of adding GVR versioning and show that generally overheads <2% are achieved. We conclude that GVR’s interfaces and implementation are flexible and portable and create a gentle-slope path to tolerate growing error rates in future systems.

97 MATHEMATICS AND COMPUTING↗

Runtime Techniques to Enable a Highly-Scalable Global Address Space Model for Petascale Computing

Over the past decade, the trajectory to the petascale has been built on increased complexity and scale of the underlying parallel architectures. Meanwhile, software developers have struggled to provide tools that maintain the productivity of computational science teams using these new systems. In this regard, Global Address Space (GAS) programming models provide a straightforward and easy to use addressing model, which can lead to improved productivity. However, the scalability of GAS depends directly on the design and implementation of the runtime system on the target petascale distributed-memory architecture. In this paper, we describe the design, implementation, and optimization of the Aggregate Remote Memory Copy Interface (ARMCI) runtime library on the Cray XT5 2.3 PetaFLOPs computer at Oak Ridge National Laboratory. We optimized our implementation with the flow intimation technique that we have introduced in this paper. Our optimized ARMCI implementation improves scalability of both the Global Arrays programming model and a real-world chemistry application—NWChem—from small jobs up through 180,000 cores.

Computer Science↗

An Uncertainty Quantification Framework for Studying the Effect of Spatial Heterogeneity in Reservoir Permeability on CO2 Sequestration

In this paper, we introduce an uncertainty quantification (UQ) software framework for carbon sequestration, focused on the effect of spatial heterogeneity of reservoir properties on CO2 migration. We use a sequential Gaussian method (SGSIM) to generate realizations of permeability fields with various spatial statistical attributes. To deal with the computational difficulties, we integrate the following ideas/approaches. First, we use three different sampling approaches (probabilistic collocation, quasi-Monte Carlo, and adaptive sampling) to reduce the number of forward calculations while trying to explore the parameter space and quantify the input uncertainty. Second, we use eSTOMP as the forward modeling simulator. eSTOMP is implemented with the Global Arrays toolkit that is based on one-sided inter-processor communication and supports a shared memory programming style on distributed memory platforms, providing a highly-scalable performance. Third, we built an adaptive system infrastructure to select the best possible data transfer mechanisms, to optimally allocate system resources to improve performance and to integrate software packages and data for composing carbon sequestration simulation, computation, analysis, estimation and visualization. We demonstrate the framework with a given CO2 injection scenario in heterogeneous sandstone reservoirs.

Carbon sequestration, uncertainty quantification, ↗

Development of a coupled thermo-hydro-mechanical model in discontinuous media for carbon sequestration

Geomechanical alteration of porous media is generally ignored for most shallow subsurface applications, whereas CO 2 injection, migration, and trapping in deep saline aquifers will be controlled by coupled multifluid flow, energy transfer, and geomechanical processes. The accurate assessment of the risks associated with potential leakage of injected CO 2 and the design of effective injection systems requires that we represent these coupled processes within numerical simulators. The objectives of this study were to develop a coupled thermal-hydro-mechanical model into a single software, and to examine the coupling of thermal, hydrological, and geomechanical processes for simulation of CO 2 injection into the subsurface for carbon sequestration. A numerical model is developed to couple nonisothermal multiphase hydrological and geomechanical processes for prediction of multiple interconnected processes for carbon sequestration in deep saline aquifers. The geomechanics model was based on Rigid Body-Spring Model (RBSM), one of the discrete methods to model discontinuous rock system. Poisson’s effect that was often ignored by RBSM was considered in the model. The simulation of large-scale and long-term coupled processes in carbon capture and storage projects requires large memory and computational performance. Global Array Toolkit was used to build the model to permit the high performance simulations of the coupled processes. The model was used to simulate a case study with several scenarios to demonstrate the impacts of considering coupled processes and Poisson’s effect for the prediction of CO 2 sequestration.

54 ENVIRONMENTAL SCIENCES↗

On the Suitability of MPI as a PGAS Runtime

Partitioned Global Address Space (PGAS) models are emerging as a popular alternative to MPI models for designing scalable applications. At the same time, MPI remains a ubiquitous communication subsystem due to its standardization, high performance, and availability on leading platforms. In this paper, we explore the suitability of using MPI as a scalable PGAS communication subsystem. We focus on the Remote Memory Access (RMA) communication in PGAS models which typically includes {\em get, put,} and {\em atomic memory operations}. We perform an in-depth exploration of design alternatives based on MPI. These alternatives include using a semantically-matching interface such as MPI-RMA, as well as not-so-intuitive interfaces such as MPI two-sided with a combination of multi-threading and dynamic process management. With an in-depth exploration of these alternatives and their shortcomings, we propose a novel design which is facilitated by the data-centric view in PGAS models. This design leverages a combination of highly tuned MPI two-sided semantics and an automatic, user-transparent split of MPI communicators to provide asynchronous progress. We implement the asynchronous progress ranks approach and other approaches within the Communication Runtime for Exascale which is a communication subsystem for Global Arrays. Our performance evaluation spans pure communication benchmarks, graph community detection and sparse matrix-vector multiplication kernels, and a computational chemistry application. The utility of our proposed PR-based approach is demonstrated by a 2.17x speed-up on 1008 processors over the other MPI-based designs.

Computational modeling↗

Barrier elision for production parallel programs

Large scientific code bases are often composed of several layers of runtime libraries, implemented in multiple programming languages. In such situation, programmers often choose conservative synchronization patterns leading to suboptimal performance. In this paper, we present context-sensitive dynamic optimizations that elide barriers redundant during the program execution. In our technique, we perform data race detection alongside the program to identify redundant barriers in their calling contexts; after an initial learning, we start eliding all future instances of barriers occurring in the same calling context. We present an automatic on-the-fly optimization and a multi-pass guided optimization. We apply our techniques to NWChem - a 6 million line computational chemistry code written in C/C++/Fortran that uses several runtime libraries such as Global Arrays, ComEx, DMAPP, and MPI. Our technique elides a surprisingly high fraction of barriers (as many as 63%) in production runs. This redundancy elimination translates to application speedups as high as 14% on 2048 cores. Our techniques also provided valuable insight about the application behavior, later used by NWChem developers. Overall, we demonstrate the value of holistic context-sensitive analyses that consider the domain science in conjunction with the associated runtime software stack.

97 MATHEMATICS AND COMPUTING↗

Scalable fault tolerance runtime technology for petascale computers (Final Project Report)

With the massive number of components comprising the forthcoming petascale computer systems, hardware failures will be routinely encountered during execution of large-scale applications. Due to the multidisciplinary, multiresolution, and multiscale nature of scientific problems that drive the demand for high end systems, applications place increasingly differing demands on the system resources: disk, network, memory, and CPU. In addition to MPI, future applications are expected to use advanced programming models such as those developed under the DARPA HPCS program as well as existing global address space programming models such as Global Arrays, UPC, and Co-Array Fortran. While there has been a considerable amount of work in fault tolerant MPI with a number of strategies and extensions for fault tolerance proposed, virtually none of advanced models proposed for emerging petascale systems is currently fault aware. To achieve fault tolerance, development of underlying runtime and OS technologies able to scale to petascale level is needed. This project has evaluated range of runtime techniques for fault tolerance for advanced programming models.

97 MATHEMATICS AND COMPUTING↗

Determination of Venus winds by ground-based radio tracking of the VEGA balloons

A global array of 20 radio observatories was used to measure the three-dimensional position and velocity of the two meteorological balloons that were injected into the equatorial region of the Venus atmosphere near Venus midnight by the VEGA spacecraft on 11 and 15 June 1985. Initial analysis of only radial velocities indicates that each balloon was blown westward about 11,500 kilometers (8,000 kilometers on the night side) by zonal winds with a mean speed of about 70 meters per second. Excursions of the data from a model of constant zonal velocity were generally less than 3 meters per second; however, a much larger variation was evident near the end of the flight of the second balloon. Consistent systematic trends in the residuals for both balloons indicate the possibility of a solar-fixed atmospheric feature. Rapid variations in balloon velocity were often detected within a single transmission (330 seconds); however, they may represent not only atmospheric motions but also self-induced aerodynamic motions of the balloon.

R A Preston↗

The Vega balloon experiment - Initial results from the global radio tracking

A unique global array of 20 radio telescopes provided 24-h telemetry acquisition of meteorological data from the Vega balloons and differential VLBI measurements of their trajectories. Initial Doppler-tracking analysis indicates mean zonal wind velocities of 69 + or - 1 and 66 + or - 1 m/sec at the Vega 1 and Vega 2 float heights, and discloses an anomaly in the Vega 2 trajectory above the mountains in Aphrodite Terra.

Preston, R. A.↗

Initial estimates of Venus winds from ground-based radio tracking of the Vega balloons

A global array of 20 radio observatories was used to measure the three-dimensional position and velocity of the two meteorological balloons that were injected into the equatorial region of the Venus atmosphere by the Vega spacecraft. Initial analysis of only radial velocities indicates that each balloon was blown westward about 11,500 kilometers (8000 kilometers on the night side) by zonal winds with a mean speed of about 70 meters per second. Excursions of the data from a model of constant zonal velocity are generally less than 3-meters per second; however, a much larger variation is evident near the end of the flight of the second balloon. Consistent systematic trends in the residual for both balloons indicate the possibility of a solar-fixed atmospheric feature.

Preston, Robert A.↗

Unconditionally stable concurrent procedures for transient finite-element analysis

A family of algorithms was outlined which would appear to be particularly well-suited for implementation in a parallel environment. This is due to the fact that for any partition of the mesh each subdomain in the partition can be processed over a time step simultaneously and independently of the rest. The method eliminates the need for assembling and factorizing large global arrays while retaining the unconditional stability properties of the algorithms used at the local level. To critically appraise the proposed methodology, two limiting cases were considered: element-by-element mesh partitions, and coarse mesh partitions. It was concluded that while the proposed methodology can be useful in sequential machines, it would appear to be promising as it bears on computation. It should also be emphasized that extensions of the method to nonlinear problems are possible.

Ortiz, Michael↗

The parsec-scale jets in 3C 273 and 3C 345

We present recent centimeter-wavelength 'global-array' VLBI images of the quasars 3C 273 and 3C 345 with dynamic ranges in excess of 3000:1. They trace the jet emission on scales from r about 5 parsecs out to about 200 parsecs. The maps of 3C273 at 18 cm wavelength show a well-collimated one-sided jet with a wavy ridge line; these wiggles exist on size scales ranging from about 1 pc to over 10 kpc. We show that the well-known superluminal flow extends to r about 120 pc. Images of 3C 345 at 6 cm wavelength from data taken in 1989 and 1990 show surprising features not seen in lower dynamic-range maps of this otherwise well-studied quasar: the inner part of the jet shows edge-brightening, which is important for modeling of jet confinement. The jet fades out very abruptly at r about 40 pc, then reappears at about 70 pc; beyond 70 pc, the resumed jet flares and is more diffuse than an extrapolation of the inner jet would predict. This morphology is reminiscent of M 87, and is suggestive of a shock wave.

Unwin, Stephen C.↗

The second Caltech-Jodrell Bank VLBI survey. 1: Observations of 91 of 193 sources

We define the sample for the second Caltech-Jodrell Bank very long base interferometry (VLBI) survey. This is a sample of 193 flat- or gigahertz-peaked-spectrum sources selected at 4850 MHz. This paper presents images of 91 sources with a resolution of approximately 1 mas, obtained using VLBI observations at 4992 MHz with a global array. The remaining images and the integrated radio spectra will be presented in a forthcoming paper by Henstock et al.

Taylor, G. B.↗

Commissioning Results on the JWST Testbed Telescope

The one-meter 18 segment JWST Testbed Telescope (TBT) has been developed at Ball Aerospace to facilitate commissioning operations for the JWST Observatory. Eight different commissioning activities were tested on the TBT: telescope focus sweep, segment ID and Search, image array, global alignment, image stacking, coarse phasing, fine phasing, and multi-field phasing. This paper describes recent commissioning results from experiments performed on the TBT.

Dean, Bruce H.↗

Radioastron Discovery of A Mini-Cocoon Around the Restarted Parsec-Scale Jet in 3C 84

We present RadioAstron space-based very long baseline interferometry (VLBI) observations of the nearby radio galaxy 3C 84 (NGC 1275) at the centre of the Perseus cluster. The observations were carried out during a perigee passage of the Spektr-R spacecraft on September 21–22, 2013 and involved a global array of 24 ground radio telescopes observing at 5 GHz and 22 GHz, together with the Space Radio Telescope (SRT). Furthermore, the Very Long Baseline Array (VLBA) and the phased Very Large Array (VLA) observed the source quasi-simultaneously at 15 GHz and 43 GHz. Fringes between the ground array and the SRT were detected on baseline lengths up to 8.1 times the Earth’s diameter, providing unprecedented resolution for 3C 84 at these wavelengths. We note that the corresponding fringe spacing is 125 μas at 5 GHz and 27 μas at 22 GHz. Our space-VLBI images reveal a previously unseen sub-structure inside the compact ∼1 pc long jet that was ejected about ten years earlier. In the 5 GHz image, we detected, for the first time, low-intensity emission from a cocoon-like structure around the restarted jet. Our results suggest that the increased power of the young jet is inflating a bubble of hot plasma as it carves its way through the ambient medium of the central region of the galaxy. Here, we estimate the minimum energy stored in the mini-cocoon, along with its pressure, volume, expansion speed, and the ratio of heavy particles to relativistic electrons, as well as the density of the ambient medium. About half of the energy delivered by the jet is dumped into the mini-cocoon and the quasi-spherical shape of the bubble suggests that this energy may be transferred to a significantly larger volume of the interstellar medium than what would be accomplished by the well-collimated jet on its own. The pressure of the hot mini-cocoon also provides a natural explanation for the almost cylindrical jet profile seen in the 22 GHz RadioAstron image.

galaxies: jets↗

Global-View Coefficients: A Data Management Solution for Parallel Quantum Monte Carlo Applications

Quantum Monte Carlo (QMC) applications perform simulation with respect to an initial state of the quantum mechanical system, which is often captured by using a cubic B-spline basis. This representation is stored as a read-only table of coefficients, and accesses to the table are generated at random as part of the Monte Carlo simulation. Current QMC applications, such as QWalk and QMCPACK, replicate this table at every process or node, which limits scalability because increasing the number of processors does not enable larger systems to be run. We present a partitioned global address space (PGAS) approach to transparently managing this data using Global Arrays in a manner that allows the memory of multiple nodes to be aggregated. We develop an automated data management system that significantly reduces communication overheads, enabling new capabilities for QMC codes. Experimental results with QWalk and QMCPACK demonstrate the effectiveness of the data management system.

Niu, Qingpeng↗