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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Linear Static Structural and Vibration Analysis on High-Performance Computers

Parallel computers offer the opportunity to significantly reduce the computation time necessary to analyze large-scale aerospace structures. This paper presents algorithms developed for and implemented on a massively-parallel computers hereafter referred to as Scalable High Performance Computers (SHPC) for the most computationally intensive tasks involved in structural analysis, namely, generation and assembly of system matrices, solution of systems of equations and calculation of the eigenvalues and eigenvectors. Results on SHPC are presented for large-scale structural problems (i.e. Models of high speed civil transport). The goal of this research is to develop new efficient technique which extend structural analysis to SHPC and make large-scale structural analyses tractable.

Baddourah, Majdi↗
NASA NTRS1994

The streaming-trapped ion interface in the equatorial inner magnetosphere

Spacecraft measurements of core ions on L=4-7 field-lines typically show trapped ion distributions near the magnetic equator, and frequently indicate field-aligned ion streams at higher latitudes. The nature of the transition between them may indicate both the microphysics of hot-cold plasma interactions and overall consequences for core plasma evolution. We have undertaken a statistical analysis and characterization of this interface and its relation to the equatorial region of the inner magnetosphere. In this analysis, we have characterized such features as the equatorial ion flux anisotropy, the penetration of field-aligned ionospheric streams into the equatorial region, the scale of the transition into trapped ion populations, and the transition latitude. We found that most transition latitudes occur within 13 deg of the equator. The typical values of equatorial ion anisotropies are consistent with bi-Maxwellian temperature ratios of T(sub perpendicular)/T(sub parallel) in the range of 3-5. The latitudinal scales for the edges of the trapped ion populations display a rather strong peak in the 2-3 deg range. We also found that there is a trend for the penetration ratio, the anisotropy half width, and the transition scale length to decrease with a higher equatorial ion anisotropy. We may interpret these features in terms of Liouville mapping of equatorially trapped ions and the reflection of the incoming ionospheric ion streams from the equatorial potential peaks associated with such trapped ions.

Lin, J.↗
NASA NTRS1994

Scalability and Portability of Two Parallel Implementations of ADI

Two domain decompositions for the implementation of the NAS Scalar Penta-diagonal Parallel Benchmark on MIMD systems are investigated, namely transposition and multi-partitioning. Hardware platforms considered are the Intel iPSC/860 and Paragon XP/S-15, and clusters of SGI workstations on ethernet, communicating through PVM. It is found that the multi-partitioning strategy offers the kind of coarse granularity that allows scaling up to hundreds of processors on a massively parallel machine. Moreover, efficiency is retained when the code is ported verbatim (save message passing syntax) to a PVM environment on a modest size cluster of workstations.

Phung, Thanh↗
NASA NTRS1993

Using parallel banded linear system solvers in generalized eigenvalue problems

Subspace iteration is a reliable and cost effective method for solving positive definite banded symmetric generalized eigenproblems, especially in the case of large scale problems. This paper discusses an algorithm that makes use of two parallel banded solvers in subspace iteration. A shift is introduced to decompose the banded linear systems into relatively independent subsystems and to accelerate the iterations. With this shift, an eigenproblem is mapped efficiently into the memories of a multiprocessor and a high speed-up is obtained for parallel implementations. An optimal shift is a shift that balances total computation and communication costs. Under certain conditions, we show how to estimate an optimal shift analytically using the decay rate for the inverse of a banded matrix, and how to improve this estimate. Computational results on iPSC/2 and iPSC/860 multiprocessors are presented.

Zhang, Hong↗
NASA NTRS1994

Using parallel banded linear system solvers in generalized eigenvalue problems

Subspace iteration is a reliable and cost effective method for solving positive definite banded symmetric generalized eigenproblems, especially in the case of large scale problems. This paper discusses an algorithm that makes use of two parallel banded solvers in subspace iteration. A shift is introduced to decompose the banded linear systems into relatively independent subsystems and to accelerate the iterations. With this shift, an eigenproblem is mapped efficiently into the memories of a multiprocessor and a high speedup is obtained for parallel implementations. An optimal shift is a shift that balances total computation and communication costs. Under certain conditions, we show how to estimate an optimal shift analytically using the decay rate for the inverse of a banded matrix, and how to improve this estimate. Computational results on iPSC/2 and iPSC/860 multiprocessors are presented.

DISTRIBUTED MEMORY MULTIPROCES↗
NASA NTRS1991

Mesoscale and large-scale variability of the Antarctic circumpolar current

This investigation of the physical oceanography of the Southern Ocean will carry out two parallel efforts during the years preceding the launch of TOPEX/POSEIDON. First, the Geosat data will be used to develop a preliminary descriptive picture of the mesoscale and large-scale, low-frequency surface circulation of the Southern Ocean. Some of this analysis of Geosat data has already begun. For example, as a measure of the geographical distribution of mesoscale variability, a color-coded map of the standard deviation of sea level from two years of Geosat data is shown. Efforts are presently under way to investigate the seasonal and year-to-year variability of this mesoscale energy. The data are also being used to generate low-pass filtered fields of sea level from which the temporal evolution of large-scale variability in the Southern Ocean may be investigated. The second parallel effort is the development and test of modeling and data assimilation techniques that will later be applied to TOPEX/POSEIDON data during the postlaunch phase. One objective of the modeling and data assimilation is to investigate the relation between mesoscale sea level variations and eddy flux in the Southern Ocean. Uncertainties in present estimates of the various components of meridional oceanic heat transport are large. The evidence presented indicates that the very energetic mesoscale variability in the ACC apparently accounts for much of the estimated 0.45x10(exp 15) watts of poleward heat transport across the ACC required to balance the heat budget. Eddy variability is strongly coherent vertically in the ACC, at least in the vicinity of Drake Passage where nearly all of the historical in situ data have been collected.

Chelton, Dudley B.↗
NASA NTRS1981

Planetary magnetospheres

A concise overview is presented of our understanding of planetary magnetospheres (and in particular, of that of the Earth), as of the end of 1981. Emphasis is placed on processes of astrophysical interest, e.g., on particle acceleration, collision-free shocks, particle motion, parallel electric fields, magnetic merging, substorms, and large scale plasma flows. The general morphology and topology of the Earth's magnetosphere are discussed, and important results are given about the magnetospheres of Jupiter, Saturn and Mercury, including those derived from the Voyager 1 and 2 missions and those related to Jupiter's satellite Io. About 160 references are cited, including many reviews from which additional details can be obtained.

Stern, D. P.↗
NASA NTRS1981

Observations of the magnetic field and plasma flow in Jupiter's magnetosheath

A comprehensive description is given of the Jovian MS magnetic fields, and explanations of these phenomena are proposed. While emphasizing Voyager 1 and 2 magnetic field observations and their relations to the plasma observations, it is also shown that the same phenomena are present in the Pioneer 10 magnetic field data. An unusually high occurrence of nearly north or south fields is observed in the outbound MS, especially in the vicinity of the MP. It is noted that the outbound MS fields and their variations tend to occur in a plane parallel to the local MP, according to large scale MP models.

Lepping, R. P.↗
NASA NTRS1984

A numerical study of the vertical transport of momentum in a tropical rainband

The vertical transport of horizontal momentum in a convective tropical rainband is studied using a two-dimensional cloud ensemble model. Twelve simulations are made under the same large-scale conditions. The vertical transports of v momentum (parallel to the rainband) are essentially the same in all of the simulations, even though the structure of the clouds is different in each of the runs. The magnitude of the v-momentum transport by clouds is fairly large. It takes only half of a day to smooth out the tropical low-level easterly jet parallel to the rainband if no other processes are operating. The vertical transports of u momentum (perpendicular to the rainband) are quite different in all of the simulations. This difference can be explained by the dissimilarities in the distributions of horizontal momentum associated with various cloud configurations. The simulated vertical transports of horizontal momentum are compared with those computed with the Schneider and Lindzen scheme. The results suggest that their scheme is basically correct and usable if some improvements are made.

Soong, S.-T.↗
NASA NTRS1986

Recent observations of the formation of filaments

Two examples of the formation of small filaments in H alpha are described and illustrated. In both cases, the formation is seen to be the spontaneous appearance of strands of absorbing mass that evolve from no previous structure. The initial development of the filaments appears to consist of the accumulation of these absorptive strands along approximately parallel paths in a channel between large-scale, opposite polarity magnetic fields on either side of the filaments. The strands exhibit continuous changes in shape and degree of absorption which can be due to successive condensations resulting in new strands, mass motions within the strands, and outflow of the mass from the strands. For at least several hours before the formation of both filaments, small-scale fragments of opposite polarity, line-of-sight magnetic flux adjacent to or immediately below the filaments, and at the ends of the filaments, were cancelling. This type of magnetic flux disappearance continued during the development of the filaments and is commonly observed in association with established filaments. Cancellation is interpreted as an important evolutionary change in the magnetic field that can lead to configurations suitable for the formation of filaments.

Martin, Sara F.↗
NASA NTRS1987

Plasma simulation using the massively parallel processor

Two dimensional electrostatic simulation codes using the particle-in-cell model are developed on the Massively Parallel Processor (MPP). The conventional plasma simulation procedure that computes electric fields at particle positions by means of a gridded system is found inefficient on the MPP. The MPP simulation code is thus based on the gridless system in which particles are assigned to processing elements and electric fields are computed directly via Discrete Fourier Transform. Currently, the gridless model on the MPP in two dimensions is about nine times slower that the gridded system on the CRAY X-MP without considering I/O time. However, the gridless system on the MPP can be improved by incorporating a faster I/O between the staging memory and Array Unit and a more efficient procedure for taking floating point sums over processing elements. The initial results suggest that the parallel processors have the potential for performing large scale plasma simulations.

Lin, C. S.↗
NASA NTRS1991

Structure of medium Mach number quasi-parallel shocks - Upstream and downstream waves

The transition from steady low-Mach-number to unsteady high-Mach-number quasi-parallel shocks was investigated by performing large-scale 1D hybrid code simulations at increasing Mach numbers. It was found that only at very low Mach number shocks the steepening is limited by upstream phase-standing whistlers, as predicted by the classical theory (Tidman and Northrop, 1968). In the intermediate region of Mach numbers between 1.5 and 3.5, a very diverse behavior is observed. Backstreaming ions generate fast magnetosonic waves which dominate the upstream, with wavelengths longer than phase-standing whistlers. At increasing Mach numbers, the phase and group velocities of the dominant waves are reduced until they point back toward the shock; when there is sufficient energy flux in these waves, they lead to unsteady shock behavior and eventually to shock reformation.

Krauss-Varban, D.↗
NASA NTRS1991

Theory and simulation of cometary shocks

In this paper, recent kinetic simulation studies of cometary bow shocks are reviewed. Cometary shocks are formed due to solar wind mass loading by water group cometary ions. This process is kinetic in nature and varies as a function of the angle between the solar wind flow velocity and the interplanetary magnetic field, as do the properties of cometary shocks. For perpendicular and parallel orientations, quasi-steady shocks with dissipation scales on the order of proton inertial length form. At oblique orientations, no steady shocks are formed; instead, the transition from supersonic to subsonic flow takes place through multiple shocklets (steepened magnetosonic waves) which are generated by the pickup ions via the resonant electromagnetic ion beam instability. This new, time dependent model of cometary bow shocks is further investigated using a large scale kinetic simulation and is compared to the observations at Comet Giacobini-Zinner.

Omidi, N.↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗
NASA NTRS1999

Large-Scale Distributed Computational Fluid Dynamics on the Information Power Grid Using Globus

This paper describes an experiment in which a large-scale scientific application development for tightly-coupled parallel machines is adapted to the distributed execution environment of the Information Power Grid (IPG). A brief overview of the IPG and a description of the computational fluid dynamics (CFD) algorithm are given. The Globus metacomputing toolkit is used as the enabling device for the geographically-distributed computation. Modifications related to latency hiding and Load balancing were required for an efficient implementation of the CFD application in the IPG environment. Performance results on a pair of SGI Origin 2000 machines indicate that real scientific applications can be effectively implemented on the IPG; however, a significant amount of continued effort is required to make such an environment useful and accessible to scientists and engineers.

Barnard, Stephen↗
NASA NTRS2004

Particle Acceleration, Magnetic Field Generation, and Emission in Relativistic Pair Jets

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating into an ambient plasma. We find that the growth times depend on the Lorenz factors of jets. The jets with larger Lorenz factors grow slower. Simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The small scale magnetic field structure generated by the Weibel instability is appropriate to the generation of "jitter" radiation from deflected electrons (positrons) as opposed to synchrotron radiation. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks.

Nishikawa, K. I.↗