Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel codes”

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 829 records · Page 46

Extending HPF for advanced data parallel applications

The stated goal of High Performance Fortran (HPF) was to 'address the problems of writing data parallel programs where the distribution of data affects performance'. After examining the current version of the language we are led to the conclusion that HPF has not fully achieved this goal. While the basic distribution functions offered by the language - regular block, cyclic, and block cyclic distributions - can support regular numerical algorithms, advanced applications such as particle-in-cell codes or unstructured mesh solvers cannot be expressed adequately. We believe that this is a major weakness of HPF, significantly reducing its chances of becoming accepted in the numeric community. The paper discusses the data distribution and alignment issues in detail, points out some flaws in the basic language, and outlines possible future paths of development. Furthermore, we briefly deal with the issue of task parallelism and its integration with the data parallel paradigm of HPF.

Chapman, Barbara↗

A general and computationally fast formulation for radiative transfer with scattering

A general formulation of monocromatic radiative transfer with scattering has been developed for plane-parallel geometry. The inhomogeneous and nonisothermal medium absorbs, emits, and anisotropically scatters radiation. Surfaces can emit and scatter radiation in any specified manner. The solution procedure uses the fact that phase incoherent scattering is linear in radiative sources. Certain basic scattering functions are then defined and calculated by an adding computer code using matrix algebra. These scattering functions are weighted by the temperature field and summed (superimposed) to obtain the solution for any specific problem. Numerical results for exiting intensities and one-sided heat fluxes from general media bound by one arbitrary surface are presented. These parametric studies demonstrate the effects of scattering particles and surfaces on radiative transfer from inhomogeneous and nonisothermal media. Application of the formulation to radiative equilibrium is also discussed. The conclusion is that all problems in plane-parallel radiative transfer with scattering can be solved by a common and computationally fast algorithm based on this formulation.

Cogley, A. C.↗

NavP: Structured and Multithreaded Distributed Parallel Programming

We present Navigational Programming (NavP) -- a distributed parallel programming methodology based on the principles of migrating computations and multithreading. The four major steps of NavP are: (1) Distribute the data using the data communication pattern in a given algorithm; (2) Insert navigational commands for the computation to migrate and follow large-sized distributed data; (3) Cut the sequential migrating thread and construct a mobile pipeline; and (4) Loop back for refinement. NavP is significantly different from the current prevailing Message Passing (MP) approach. The advantages of NavP include: (1) NavP is structured distributed programming and it does not change the code structure of an original algorithm. This is in sharp contrast to MP as MP implementations in general do not resemble the original sequential code; (2) NavP implementations are always competitive with the best MPI implementations in terms of performance. Approaches such as DSM or HPF have failed to deliver satisfying performance as of today in contrast, even if they are relatively easy to use compared to MP; (3) NavP provides incremental parallelization, which is beyond the reach of MP; and (4) NavP is a unifying approach that allows us to exploit both fine- (multithreading on shared memory) and coarse- (pipelined tasks on distributed memory) grained parallelism. This is in contrast to the currently popular hybrid use of MP+OpenMP, which is known to be complex to use. We present experimental results that demonstrate the effectiveness of NavP.

navigational programming (NavP)↗

Modular VLSI Reed-Solomon Decoder

Proposed Reed-Solomon (RS) decoder assembled from very-large-scale integrated-circuit (VLSI) building blocks. Decoder exploits recursive forms in RS decoding algorithms. RS codes capable of correcting random or burst errors in telemetry and other data-communication signals. Because of small size and low power consumption, advantageous to employ several such decoders in parallel-processing scheme to increase decoding speed.

Liu, K. Y.↗

The Navier-Stokes computer

The Navier-Stokes computer (NSC) has been developed for solving problems in fluid mechanics involving complex flow simulations that require more speed and capacity than provided by current and proposed Class VI supercomputers. The machine is a parallel processing supercomputer with several new architectural elements which can be programmed to address a wide range of problems meeting the following criteria: (1) the problem is numerically intensive, and (2) the code makes use of long vectors. A simulation of two-dimensional nonsteady viscous flows is presented to illustrate the architecture, programming, and some of the capabilities of the NSC.

Nosenchuck, D. M.↗

Automating the multiprocessing environment

An approach to automate the programming and operation of tree-structured networks of multiprocessor systems is discussed. A conceptual, knowledge-based operating environment is presented, and requirements for two major technology elements are identified as follows: (1) An intelligent information translator is proposed for implementating information transfer between dissimilar hardware and software, thereby enabling independent and modular development of future systems and promoting a language-independence of codes and information; (2) A resident system activity manager, which recognizes the systems capabilities and monitors the status of all systems within the environment, is proposed for integrating dissimilar systems into effective parallel processing resources to optimally meet user needs. Finally, key computational capabilities which must be provided before the environment can be realized are identified.

Arpasi, Dale J.↗

An experimental and computational study of rotor-vortex interactions

An experimental and computational study has been performed on a close rotor-blade/vortex interaction. Surface pressure data was obtained from a rotor operating close to the tip-vortex from an upstream wing in a wind tunnel. Data was obtained for a wide range of blade-vortex proximities, orientations, and blade-tip Mach numbers (up to the transonic regime). A numerical model of these interactions was constructed using the unsteady, three-dimensional, full-potential rotor code called FPR. The model employed an undistorted full-field representation of the measured vortex. This simple model gave excellent comparisons with the data for a wide range of conditions, including parallel head-on interactions. Computational studies have also been performed on the manner of vortex representation and the influence of vortex-core size.

Caradonna, Francis X.↗

Concurrent Finite-Element Analysis On Hypercube Computers

Improved approach to execution of finite-element codes on hypercube and similar concurrent data processors increases efficiency of computation for many different types of problems. Based on flexible general model of computation on, and communication among, large-node parallel processors. Hybrid combining direct methods within subdomains and preconditioned-conjugate-gradient (PCG) iteration on remaining boundary system to obtain method of solution both robust and efficient. Preserves general structure and function of conventional sequential finite-element-method software. Generalized to both distributed- and shared-memory multicomputers, eliminating degree of machine specificity restricting general usefulness.

Lyzenga, G. A.↗

Root damage analysis of aircraft engine blade subject to ice impact

The blade root response due to ice impact on an engine blade is simulated using the NASA in-house code BLASIM. The ice piece is modeled as an equivalent spherical object impacting on the leading edge of the blade and has the velocity opposite to that of the aircraft with direction parallel to the engine axis. The effect of ice impact is considered to be an impulse load on the blade with its amplitude computed based on the momentum transfer principle. The blade response due to the impact is carried out by modal superposition using the first three modes. The maximum dynamic stresses at the blade root are computed at the quarter cycle of the first natural frequency. A combined stress failure function based on modified distortion energy is used to study the spanwise bending damage response at the blade root. That damage function reaches maximum value for very low ice speeds and increases steeply with increases in engine speed.

Reddy, E. S.↗

Progress in Parallel Schur Complement Preconditioning for Computational Fluid Dynamics

We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a non-overlapping Schur complement procedure for arbitrary triangulated domains. The ultimate goal of this research is to develop scalable preconditioning algorithms for fluid flow discretizations on parallel computing architectures. In our implementation of the Schur complement preconditioning technique, the triangulation is first partitioned into a number of subdomains using the METIS multi-level k-way partitioning code. This partitioning induces a natural 2X2 partitioning of the p.d.e. discretization matrix. By considering various inverse approximations of the 2X2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 and the SGI Power Challenge array using MPI message passing protocol. A number of example CFD calculations will be presented to illustrate and assess various Schur complement approximations.

Barth, Timothy J.↗

Particle acceleration magnetic field generation, and emission in Relativistic pair jets

Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) are responsible for particle acceleration in relativistic pair jets. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic pair jet propagating through a pair plasma. Simulations show that the Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. Simulation results show that this instability generates and amplifies highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The "jitter' I radiation from deflected electrons can have different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. The growth rate of the Weibel instability and the resulting particle acceleration depend on the magnetic field strength and orientation, and on the initial particle distribution function. In this presentation we explore some of the dependencies of the Weibel instability and resulting particle acceleration on the magnetic field strength and orientation, and the particle distribution function.

Nishikawa, K.-I.↗

Impact of Airfoils on Aerodynamic Optimization of Heavy Lift Rotorcraft

Rotor airfoils were developed for two large tiltrotor designs, the Large Civil Tilt Rotor (LCTR) and the Military Heavy Tilt Rotor (MHTR). The LCTR was the most promising of several rotorcraft concepts produced by the NASA Heavy Lift Rotorcraft Systems Investigation. It was designed to carry 120 passengers for 1200 nm, with performance of 350 knots cruise at 30,000 ft altitude. A parallel design, the MHTR, had a notional mission of 40,000 Ib payload, 500 nm range, and 300 knots cruise at 4000 ft, 95 F. Both aircraft were sized by the RC code developed by the U. S. Army Aeroflightdynamics Directorate (AFDD). The rotors were then optimized using the CAMRAD II comprehensive analysis code. Rotor airfoils were designed for each aircraft, and their effects on performance analyzed by CAMRAD II. Airfoil design criteria are discussed for each rotor. Twist and taper optimization are presented in detail for each rotor, with discussions of performance improvements provided by the new airfoils, compared to current technology airfoils. Effects of stall delay and blade flexibility on performance are also included.

Acree, Cecil W., Jr.↗

Receiver Would Control Phasing of a Phased-Array Antenna

In a proposed digital signal-processing technique, a radio receiver would control the phasing of a phased-array antenna to aim the peaks of the antenna radiation pattern toward desired signal sources while aiming the nulls of the pattern toward interfering signal sources. The technique was conceived for use in a Global Positioning System (GPS) receiver, for which the desired signal sources would be GPS satellites and typical interference sources would be terrestrial objects that cause multipath propagation. The technique could also be used to optimize reception in spread-spectrum cellular-telephone and military communication systems. During reception of radio signals in a conventional phased-array antenna system, received signals at their original carrier frequencies are phase-shifted, then combined by analog circuitry. The combination signal is then subjected to down-conversion and demodulation. In a system according to the proposed technique (see figure), the signal received by each antenna would be subjected to down-conversion, spread-spectrum demodulation, and correlation; this processing would be performed separately from, and simultaneously with, similar processing of signals received by the other antenna elements. Following analog down-conversion to baseband, the signals would be digitized, and all subsequent processing would be digital. In the digital process, residual carriers would be removed and each signal would be correlated with a locally generated model pseudorandum-noise code, all following normal GPS procedure. As part of this procedure, accumulated values would be added in software and the resulting signals would be phase-shifted in software by the amounts necessary to synthesize the desired antenna directional gain pattern of peaks and nulls. The principal advantage of this technique over the conventional radio-frequency-combining technique is that the parallel digital baseband processing of the signals from the various antenna elements would be a relatively inexpensive and flexible means for exploiting the inherent multiple-peak/multiple-null aiming capability of a phased-array antenna. In the original intended GPS application, the peaks and nulls could be directed independently for each GPS signal being tracked by the GPS receiver. The technique could also be applied to other code-division multiple-access communication systems.

Dunn, Charles E.↗

Divergence study of a high-aspect ratio, forward-swept wing

An experimental wind-tunnel study to determine the divergence characteristics of a high-aspect ratio, forward-swept wing has been conducted in the NASA Langley Research Center (LaRC) Transonic Dynamics Tunnel (TDT). The rectangular wing used for this study had a panel aspect ratio of 9.16 (lambda = 0 deg.) and the sweep angle could be set at lambda = 0 deg., -15 deg., -30 deg., -45 deg., or -60 deg. A rectangular wing tip shape was tested at each of these sweep angles. In addition, a tip shape parallel to the freestream flow was tested for a wing sweep angle of lambda = -45 deg. The root of the wing was cantilever mounted to the wall of the wind tunnel. Divergence conditions were measured at M = 0.4 for each sweep angle and tip configuration tested. Subcritical response techniques were used to extrapolate to the divergence conditions during the wind-tunnel test. The primary objective of this test was to obtain data which could be used to verify for this configuration the divergence prediction capability of an aeroelastic analysis code. Subsonic lifting surface theory (kernel function) aerodynamics are utilized by this particular code. The analytical predictions of divergence were found to be significantly conservative at all forward sweep angles. At lambda = -45 deg., the analysis was 14 percent conservative. The effect of the two tip shapes on the divergence dynamic pressure was predicted accurately by the analysis. The divergence condition for the tip shape parallel to the flow occurred at a dynamic pressure 14 percent higher than the divergence condition with a rectangular tip shape.

Cole, S. R.↗

Analysis of supersonic plug nozzle flowfield and heat transfer

A number of problems pertaining to the flowfield in a plug nozzle, designed as a supersonic thruster nozzle, with provision for cooling the plug with a coolant stream admitted parallel to the plug wall surface, were studied. First, an analysis was performed of the inviscid, nonturbulent, gas dynamic interaction between the primary hot stream and the secondary coolant stream. A numerical prediction code for establishing the resulting flowfield with a dividing surface between the two streams, for various combinations of stagnation and static properties of the two streams, was utilized for illustrating the nature of interactions. Secondly, skin friction coefficient, heat transfer coefficient and heat flux to the plug wall were analyzed under smooth flow conditions (without shocks or separation) for various coolant flow conditions. A numerical code was suitably modified and utilized for the determination of heat transfer parameters in a number of cases for which data are available. Thirdly, an analysis was initiated for modeling turbulence processes in transonic shock-boundary layer interaction without the appearance of flow separation.

Murthy, S. N. B.↗

Fault-tolerant multichannel demultiplexer subsystems

Fault tolerance in future processing and switching communication satellites is addressed by showing new methods for detecting hardware failures in the first major subsystem, the multichannel demultiplexer. An efficient method for demultiplexing frequency slotted channels uses multirate filter banks which contain fast Fourier transform processing. All numerical processing is performed at a lower rate commensurate with the small bandwidth of each bandbase channel. The integrity of the demultiplexing operations is protected by using real number convolutional codes to compute comparable parity values which detect errors at the data sample level. High rate, systematic convolutional codes produce parity values at a much reduced rate, and protection is achieved by generating parity values in two ways and comparing them. Parity values corresponding to each output channel are generated in parallel by a subsystem, operating even slower and in parallel with the demultiplexer that is virtually identical to the original structure. These parity calculations may be time shared with the same processing resources because they are so similar.

Redinbo, Robert↗

Loop heating by D.C. electric current and electromagnetic wave emissions simulated by 3-D EM particle zone

We have shown that a current-carrying plasma loop can be heated by magnetic pinch driven by the pressure imbalance between inside and outside the loop, using a 3-dimensional electromagnetic (EM) particle code. Both electrons and ions in the loop can be heated in the direction perpendicular to the ambient magnetic field, therefore the perpendicular temperature can be increased about 10 times compared with the parallel temperature. This temperature anisotropy produced by the magnetic pinch heating can induce a plasma instability, by which high-frequency electromagnetic waves can be excited. The plasma current which is enhanced by the magnetic pinch can also excite a kinetic kink instability, which can heat ions perpendicular to the magnetic field. The heating mechanism of ions as well as the electromagnetic emission could be important for an understanding of the coronal loop heating and the electromagnetic wave emissions from active coronal regions.

Sakai, J. I.↗