Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel, scaling”

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 685 records · Page 38

IN13B-1660: Analytics and Visualization Pipelines for Big Data on the NASA Earth Exchange (NEX) and OpenNEX

We are developing capabilities for an integrated petabyte-scale Earth science collaborative analysis and visualization environment. The ultimate goal is to deploy this environment within the NASA Earth Exchange (NEX) and OpenNEX in order to enhance existing science data production pipelines in both high-performance computing (HPC) and cloud environments. Bridging of HPC and cloud is a fairly new concept under active research and this system significantly enhances the ability of the scientific community to accelerate analysis and visualization of Earth science data from NASA missions, model outputs and other sources. We have developed a web-based system that seamlessly interfaces with both high-performance computing (HPC) and cloud environments, providing tools that enable science teams to develop and deploy large-scale analysis, visualization and QA pipelines of both the production process and the data products, and enable sharing results with the community. Our project is developed in several stages each addressing separate challenge - workflow integration, parallel execution in either cloud or HPC environments and big-data analytics or visualization. This work benefits a number of existing and upcoming projects supported by NEX, such as the Web Enabled Landsat Data (WELD), where we are developing a new QA pipeline for the 25PB system.

visualization↗

DSMC analysis in a heterogeneous parallel computing environment

A methodology for implementing parallel DSMC codes in a heterogeneous computing environment is described. The methodology involves the use of a common message-passing software library together with recently developed software that handles the actual interprocessor communications in a standard manner across a variety of computing platforms. Benchmark tests using a simple DSMC model problem were performed on an Intel iPSC/860, a Cray-YMP and a group of Sun workstations. The approach was found to give speedups that scaled linearly with problem size on all the computing platforms tested. This methodology was then incorporated into a production-type DSMC code to allow the simulation of problems that would not otherwise have been practical. The application of this production code to simulations of hypersonic shear flows and shock-lip interactions under near-continuum conditions is described. Synchronous and asynchronous models for implementing parallelism into DSMC simulations are also described and both models are shown to produce the same steady-state result.

Wilmoth, R. G.↗

Magnetic clouds in the solar wind

Two interplanetary magnetic clouds, characterized by anomalous magnetic field directions and unusually high magnetic field strengths with a scale of the order of 0.25 AU, are identified and described. As the clouds moved past a spacecraft located in the solar wind near Earth, the magnetic field direction changed by rotating approximately 180 deg nearly parallel to a plane which was essentially perpendicular to the ecliptic. The configuration of the magnetic field in the clouds might be that of a tightly wound cylindrical helix or a series of closed circular loops. One of the magnetic clouds was in a cold stream preceded by a shock, and it caused both a geomagnetic storm and a depression in the galactic cosmic ray intensity. No stream, geomagnetic storm, or large cosmic ray decrease was associated with the other magnetic cloud.

Burlaga, L. F.↗

The suitability of the ILLIAC IV architecture for image processing

The major architectural features of the ILLIAC IV large scale, array processor are summarized along with their applicability to image processing. Several image processing algorithms are considered, including multispectral classification, texture feature extraction, two-dimensional Fourier transform, and synthetic aperture radar processing. The basic parallelism of the ILLIAC IV (64 processing elements acting in lock-step) is usually fully utilized by the image processing applications. The major architectural aspect of the system with respect to image processing is the relatively small local scratch-pad memory and the long latency time to access the main storage device. The major precision used for the image processing applications is the 32-bit floating point, given a choice of 8-bit integers and 64-bit floating point.

Stevenson, D. K.↗

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

Hubble Space Telescope imaging of Eta Carinae

New high spatial resolution observations of the material around Eta Carinae, obtained with the Hubble Space Telescope Wide Field/Planetary Camera, are presented. The star Eta Carinae is one of the most massive and luminous stars in the Galaxy, and has been episodically expelling significant quantities of gas over the last few centuries. The morphology of the bright central nebulosity (the homunculus) indicates that it is a thin shell with very well defined edges, and is clumpy on 0.2 arcsec (about 10 to the 16th cm) scales. An extension to the northeast of the star (NN/NS using Walborn's 1976 nomenclature) appears to be a stellar jet and its associated bow shock. The bow shock is notable for an intriguing series of parallel linear features across its face. The S ridge and the W arc appear to be part of a 'cap' of emission located to the SW and behind the star. Together, the NE jet and the SW cap suggest that the symmetry axis for the system runs NE-SW rather than SE-NW, as previously supposed. Overall, the data indicate that the material around the star may represent an oblate shell with polar blowouts, rather than a bipolar flow.

Hester, J. J.↗

From Oxygen Generation to Metals Production: In Situ Resource Utilization by Molten Oxide Electrolysis

For the exploration of other bodies in the solar system, electrochemical processing is arguably the most versatile technology for conversion of local resources into usable commodities: by electrolysis one can, in principle, produce (1) breathable oxygen, (2) silicon for the fabrication of solar cells, (3) various reactive metals for use as electrodes in advanced storage batteries, and (4) structural metals such as steel and aluminum. Even so, to date there has been no sustained effort to develop such processes, in part due to the inadequacy of the database. The objective here is to identify chemistries capable of sustaining molten oxide electrolysis in the cited applications and to examine the behavior of laboratory-scale cells designed to generate oxygen and to produce metal. The basic research includes the study of the underlying high-temperature physical chemistry of oxide melts representative of lunar regolith and of Martian soil. To move beyond empirical approaches to process development, the thermodynamic and transport properties of oxide melts are being studied to help set the limits of composition and temperature for the processing trials conducted in laboratory-scale electrolysis cells. The goal of this investigation is to deliver a working prototype cell that can use lunar regolith and Martian soil to produce breathable oxygen along with metal by-product. Additionally, the process can be generalized to permit adaptation to accommodate different feedstock chemistries, such as those that will be encountered on other bodies in the solar system. The expected results of this research include: (1) the identification of appropriate electrolyte chemistries; (2) the selection of candidate anode and cathode materials compatible with electrolytes named above; and (3) performance data from a laboratory-scale cell producing oxygen and metal. On the strength of these results it should be possible to assess the technical viability of molten oxide electrolysis for in situ resource utilization on the Moon and Mars. In parallel, there may be commercial applications here on earth, such as new green technologies for metals extraction and for treatment of hazardous waste, e.g., fixing heavy metals.

Khetpal, Deepak↗

IRAS 16293-2422: Evidence for Infall onto a Counter-Rotating Protostellar Accretion Disk

We report high spatial resolution VLA observations of the low-mass star-forming region IRAS 16293-2422 using four molecular probes: ethyl cyanide (CH3CH2CN)) methyl formate (CH3OCHO), formic acid (HCOOH), and the ground vibrational state of silicon monoxide (SiO). Ethyl cyanide emission has a spatial scale of approx. 20" and encompasses binary cores A and B as determined by continuum emission peaks. Surrounded by formic acid emission, methyl formate emission has a spatial scale of approx. 6" and is confined to core B. SiO emission shows two velocity components with spatial scales less than 2" that map approx. 2" northeast of the A and B symmetry axis. The redshifted SiO is approx. 2" northwest of blueshifted SiO along a position angle of approx. 135deg which is approximately parallel to the A and B symmetry axis. We interpret the spatial position offset in red and blueshifted SiO emission as due to rotation of a protostellar accretion disk and we derive approx. 1.4 Solar Mass, interior to the SiO emission. In the same vicinity, Mundy et al. (1986) also concluded rotation of a nearly edge-on disk from OVRO observations of much stronger and ubiquitous CO-13 emission but the direction of rotation is opposite to the SiO emission findings. Taken together, SiO and CO-13 data suggest evidence for a counter-rotating disk. Moreover, archival BIMA array CO-12C data show an inverse P Cygni profile with the strongest absorption in close proximity to the SiO emission, indicating unambiguous material infall toward the counter-rotating protostellar disk at a new source location within the IRAS 16293-2422 complex. The details of these observations and our interpretations are discussed.

Remijan, Anthony J.↗

A balanced submatrix merging algorithm for multiprocessor architectures

In this article, a parallel algorithm which applies Givens rotations to selectively annihilate k(k + 1)/2 nonzero elements from two k x n(k not more than n) upper trapezoidal submatrices is described. The new algorithm is suitable for implementation on either a pair of directly connected local-memory processors or two clusters of multiple tightly-coupled processors. Analyses show that in both cases the proposed algorithms achieve optimal speed-up by balancing the work load distribution and masking interprocessor or intercluster communication by computation if k is much small than n. In the context of solving large scale least squares problems, this submatrix merging step is repetitively needed during the entire computation and, furthermore, there are usually many pairs of such submatrices to be merged with each submatrix stored in the memory of a processor or a cluster of processors. The proposed algorithm can be applied to each pair of submatrices concurrently, and thus parallelizes an important step in solving the least squares problems.

Chu, Eleanor↗

Optoelectronic Tool Adds Scale Marks to Photographic Images

A simple, easy-to-use optoelectronic tool projects scale marks that become incorporated into photographic images (including film and electronic images). The sizes of objects depicted in the images can readily be measured by reference to the scale marks. The role played by the scale marks projected by this tool is the same as that of the scale marks on a ruler placed in a scene for the purpose of establishing a length scale. However, this tool offers the advantage that it can put scale marks quickly and safely in any visible location, including a location in which placement of a ruler would be difficult, unsafe, or time-consuming. The tool (see Figure 1) includes an aluminum housing, within which are mounted four laser diodes that operate at a wavelength of 670 nm. The laser diodes are spaced 1 in. (2.54 cm) apart along a baseline. The laser diodes are mounted with setscrews, which are used to adjust their beams to make them all parallel to each other and perpendicular to the baseline. During the adjustment process, the effect of the adjustments is observed by measuring the positions of the laser-beam spots on a target 80 ft (approx.24 m) away. Once the adjustments have been completed, the laser beams define three 1-in. (2.54-cm) intervals and the location of each beam is defined to within 1/16 in. (approx.1.6 mm) at any target distance out to about 80 ft (approx.24 m). The distance between the laser-beam spots as seen in an image is strictly defined only along an axis parallel to the baseline and perpendicular to the laser beam (also perpendicular to the line of sight of the camera, assuming that the camera-to-target distance is much greater than the distance between the tool and the camera lens). If a flat target surface illuminated by the laser beams is tilted with respect to the aforesaid axis, then the distance along the target surface between scale marks is proportional to the secant of the tilt angle. If one knows the tilt angle, one can correct for it. Even if one does not know the tilt angle precisely, it may not matter: For example, at a tilt of 10 , the secant is approximately 1.0154, so that the tilt error is only about 1.54 percent, which is negligibly small for a typical application in which only approximate measurements are needed.

Stevenson, Charlie↗

Numerical simulations of double layers and auroral electric fields

Recent one-dimensional and two-dimensional numerical simulations of double layers (DLs) in the electric fields of the auroral plasma are reviewed, with reference to observational data. It is found that two-dimensional DLs driven by current sheets of finite thickness have different characteristics, depending on whether the layer thickness is less than or much greater than the ion gyroradius: When thickness is less than ion gyroradius, V-shaped DLs form with nearly equal parallel and perpendicular potential drops; when layer thickness is much greater than ion gyroradius the major parallel potential drop occurs outside the current sheet and the perpendicular electric fields are localized at the edges of the current sheet. It is shown that some features of the simulated fields, such as the amplitudes and scale lengths, are qualitatively similar to those observed in space.

Singh, N.↗

Computational Issues in Damping Identification for Large Scale Problems

Two damping identification methods are tested for efficiency in large-scale applications. One is an iterative routine, and the other a least squares method. Numerical simulations have been performed on multiple degree-of-freedom models to test the effectiveness of the algorithm and the usefulness of parallel computation for the problems. High Performance Fortran is used to parallelize the algorithm. Tests were performed using the IBM-SP2 at NASA Ames Research Center. The least squares method tested incurs high communication costs, which reduces the benefit of high performance computing. This method's memory requirement grows at a very rapid rate meaning that larger problems can quickly exceed available computer memory. The iterative method's memory requirement grows at a much slower pace and is able to handle problems with 500+ degrees of freedom on a single processor. This method benefits from parallelization, and significant speedup can he seen for problems of 100+ degrees-of-freedom.

Pilkey, Deborah L.↗

Oxygen-rich Mira variables: Near-infrared luminosity calibrations. Populations and period-luminosity relations

Hipparcos astrometric and kinematical data of oxygen-rich Mira variables are used to calibrate absolute near-infrared magnitudes and kinematic parameters. Three distinct classes of stars with different kinematics and scale heights were identified. The two most significant groups present characteristics close to those usually assigned to extended/thick disk-halo populations and old disk populations, respectively, and thus they may differ by their metallicity abundance. Two parallel period-luminosity relations are found, one for each population. The shift between these relations is interpreted as the consequence of the effects of metallicity abundance on the luminosity.

Alvarez, R.↗

Aerodynamic Characteristics of a Large-Scale Unswept Wing-Body-Tail Configuration with Blowing Applied Over the Flap and Wind Leading Edge

An investigation has been conducted in the Langley full-scale tunnel to determine the effects of a blowing boundary-layer-control lift-augmentation system on the aerodynamic characteristics of a large-scale model of a fighter-type airplane. The wing was unswept at the 70-percent- chord station, had an aspect ratio of 2.86, a taper ratio of 0.40, and 4-percent-thick biconvex airfoil sections parallel to the plane of symmetry. The tests were conducted over a range of angles of attack from approximately -4 deg to 23 deg for a Reynolds number of approximately 5.2 x 10(exp 6) which corresponds to a Mach number of 0.08. Blowing rates were normally restricted to values just sufficient to control air-flow separation. The results of this investigation showed that wing leading-edge blowing in combination with large values of wing leading-edge-flap deflection was a very effective leading-edge flow-control device for wings having highly loaded trailing-edge flaps. With leading-edge blowing there was no hysteresis of the lift, drag, and pitching-moment characteristics upon recovery from stall. End plates were found to improve the lift and drag characteristics of the test configuration in the moderate angle-of-attack range, and blockage to one-quarter of the blowing-slot area was not detrimental to the aerodynamic characteristics. Blowing boundary-layer control resulted in a considerably reduced landing speed and reduced landing and take-off distances. The ailerons were very effective lateral-control devices when used with blowing flaps.

McLemore, H. Clyde↗

Characterization of Clouds and the Anisotropy of Emitted and Reflected Radiances for the Purpose of Obtaining the Radiative Heating of the Atmosphere

The goal of the work supported through this grant was to assess the validity of the assumptions underlying the CERES Strategy for determining radiative fluxes. Specifically, the work focused on the determination of scene type and the use of anisotropic factors to derive radiative fluxes from observed broadband radiances. The work revealed a dependence of the anisotropy of reflected and emitted broadband radiances on the spatial resolution of the observations that had been overlooked in the formulation of the CERES strategy. This dependence on spatial resolution coupled with errors in scene identification led to view zenith angle dependent biases in the ERBE derived radiative fluxes. Scene identification will be greatly improved in CERES thereby alleviating somewhat the biases arising from the dependence of the anisotropy of the radiances on spatial resolution. Attention was then focused on the validity of plane-parallel radiative transfer theory which is relied on to characterize the scene types viewed by the CERES scanner. Again, viewing geometry dependent biases were found even for single-layered, overcast cloud systems. Such systems are taken to be the closest examples of plane-parallel clouds. At least some of the departures from plane-parallel behavior were evidently due to relatively small bumps on the tops of extensive stratus layers. The bumps cannot be resolved in the imagery that will be used to characterize the scenes viewed by the CERES scanner. As part of this investigation, the ice sheets of Greenland and Antarctica were shown to provide radiometrically stable targets for determining the visible and near infrared calibrations of radiometers. These targets were used to calibrate the reflected sunlight at visible wavelengths used in this study. Finally, the limitations of plane-parallel theory notwithstanding, the common practice of ignoring fractional cloud cover within the fields of view of imaging radiometers was shown to lead to biases in the retrieved cloud properties. The development of retrievals for pixel-scale cloud cover fraction is an attempt to reduce such bases. Work on these retrievals continues.

Coakley, James A., Jr.↗

Large-scale electric fields resulting from magnetic reconnection in the corona

The method of Forbes and Priest (2-D model) is applied to the large two-ribbon flare of July 29, 1973, for which both detailed H observations and magnetic data are available. For this flare the ribbons were long, nearly straight, and parallel to each other, and the 2-D model for the coronal field geometry may be adequate. The temporal profile E(t) is calculated and indicates that reconnection sets in at the beginning of the decay phase. From this time the electric field grows rapidly to a maximum value of about 2 V/cm within just a few minutes. Thereafter it decreases monotonically with time.

Kopp, R. A.↗

Preliminary large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids

Large-eddy simulation (LES) has matured to the point where application to complex flows is desirable. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential. In the past two years, important building blocks have been put into place making possible a careful study of LES on unstructured grids. The first building block was an efficient mesh generator which allowed the placement of points according to smooth variation of physical length scales. This variation of length scales is in all three directions independently, which allows a large reduction in points when compared to structured-grid methods, which can only vary length scales in one direction at a time. The second building block was the development of a dynamic model appropriate for unstructured grids. The principle obstacle was the development of an unstructured-grid filtering operator. In the past year, some of the new filters developed by Jansen have been implemented into a highly parallelized finite element code based on the Galerkin/least-squares finite element method. We have chosen the NACA 4412 airfoil at maximum lift as the first simulation for a variety of reasons. First, it is a problem of significant interest since it would be the first LES of an aircraft component. Second, this flow has been the subject of three experimental studies. The third reason for considering this flow is the variety of flow features which provide an important test of the dynamic model. Only the dynamic model can be expected to perform satisfactorily in this variety of situations: from the laminar regions where it must not modify the flow at all to the turbulent boundary layers and wake where it must represent a wide variety of subgrid-scale structures. The flow configuration we have chosen is that of Wadcock (1987) at Reynolds number based on chord Re(sub c) = u(sub infinity)c/v = 1.64 x 10(exp 6), Mach number M = 0.2, and 12 deg angle of attack.

Jansen, Kenneth↗

The force on the flex: Global parallelism and portability

A parallel programming methodology, called the force, supports the construction of programs to be executed in parallel by an unspecified, but potentially large, number of processes. The methodology was originally developed on a pipelined, shared memory multiprocessor, the Denelcor HEP, and embodies the primitive operations of the force in a set of macros which expand into multiprocessor Fortran code. A small set of primitives is sufficient to write large parallel programs, and the system has been used to produce 10,000 line programs in computational fluid dynamics. The level of complexity of the force primitives is intermediate. It is high enough to mask detailed architectural differences between multiprocessors but low enough to give the user control over performance. The system is being ported to a medium scale multiprocessor, the Flex/32, which is a 20 processor system with a mixture of shared and local memory. Memory organization and the type of processor synchronization supported by the hardware on the two machines lead to some differences in efficient implementations of the force primitives, but the user interface remains the same. An initial implementation was done by retargeting the macros to Flexible Computer Corporation's ConCurrent C language. Subsequently, the macros were caused to directly produce the system calls which form the basis for ConCurrent C. The implementation of the Fortran based system is in step with Flexible Computer Corporations's implementation of a Fortran system in the parallel environment.

Jordan, H. F.↗