Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel simulation”

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,045 records · Page 58

Initial Kernel Timing Using a Simple PIM Performance Model

This presentation will describe some initial results of paper-and-pencil studies of 4 or 5 application kernels applied to a processor-in-memory (PIM) system roughly similar to the Cascade Lightweight Processor (LWP). The application kernels are: * Linked list traversal * Sun of leaf nodes on a tree * Bitonic sort * Vector sum * Gaussian elimination The intent of this work is to guide and validate work on the Cascade project in the areas of compilers, simulators, and languages. We will first discuss the generic PIM structure. Then, we will explain the concepts needed to program a parallel PIM system (locality, threads, parcels). Next, we will present a simple PIM performance model that will be used in the remainder of the presentation. For each kernel, we will then present a set of codes, including codes for a single PIM node, and codes for multiple PIM nodes that move data to threads and move threads to data. These codes are written at a fairly low level, between assembly and C, but much closer to C than to assembly. For each code, we will present some hand-drafted timing forecasts, based on the simple PIM performance model. Finally, we will conclude by discussing what we have learned from this work, including what programming styles seem to work best, from the point-of-view of both expressiveness and performance.

BRIEFING CHARTS↗

Implementation of a fully-balanced periodic tridiagonal solver on a parallel distributed memory architecture

While parallel computers offer significant computational performance, it is generally necessary to evaluate several programming strategies. Two programming strategies for a fairly common problem - a periodic tridiagonal solver - are developed and evaluated. Simple model calculations as well as timing results are presented to evaluate the various strategies. The particular tridiagonal solver evaluated is used in many computational fluid dynamic simulation codes. The feature that makes this algorithm unique is that these simulation codes usually require simultaneous solutions for multiple right-hand-sides (RHS) of the system of equations. Each RHS solutions is independent and thus can be computed in parallel. Thus a Gaussian elimination type algorithm can be used in a parallel computation and the more complicated approaches such as cyclic reduction are not required. The two strategies are a transpose strategy and a distributed solver strategy. For the transpose strategy, the data is moved so that a subset of all the RHS problems is solved on each of the several processors. This usually requires significant data movement between processor memories across a network. The second strategy attempts to have the algorithm allow the data across processor boundaries in a chained manner. This usually requires significantly less data movement. An approach to accomplish this second strategy in a near-perfect load-balanced manner is developed. In addition, an algorithm will be shown to directly transform a sequential Gaussian elimination type algorithm into the parallel chained, load-balanced algorithm.

Eidson, T. M.↗

Parallel processing methods for space based power systems

This report presents a method for doing load-flow analysis of a power system by using a decomposition approach. The power system for the Space Shuttle is used as a basis to build a model for the load-flow analysis. To test the decomposition method for doing load-flow analysis, simulations were performed on power systems of 16, 25, 34, 43, 52, 61, 70, and 79 nodes. Each of the power systems was divided into subsystems and simulated under steady-state conditions. The results from these tests have been found to be as accurate as tests performed using a standard serial simulator. The division of the power systems into different subsystems was done by assigning a processor to each area. There were 13 transputers available, therefore, up to 13 different subsystems could be simulated at the same time. This report has preliminary results for a load-flow analysis using a decomposition principal. The report shows that the decomposition algorithm for load-flow analysis is well suited for parallel processing and provides increases in the speed of execution.

Berry, F. C.↗

Propagation and absorption of electron-cyclotron maser radiation during solar flares

The propagation and absorption of the maser radiation during solar flares are examined through linear theory and electromagnetic particle simulations. It is shown using linear theory that strong absorption of the radiation should occur as it propagates toward the second harmonic layer, where the magnetic field is half as strong as in the emission region. Only radiation propagating nearly parallel to the magnetic field in a low-temperature plasma may be able to escape under certain limited conditions. Finite temperature effects can cause radiation propagating nearly perpendicular to the magnetic field to refract, causing enhanced absorption. Particle simulations are then used to evaluate the nonlinear response of the plasma as the maser radiation propagates through the absorption layer. It is shown that some of the maser radiation is able to escape through a process of absorption below the second harmonic of the local gyrofrequency and re-emission above it. The fraction able to escape is much higher than that predicted by linear theory, although the amount of escaping energy is only a small fraction of the incident energy.

Mckean, M. E.↗

Multiple curved descending approaches and the air traffic control problem

A terminal area air traffic control simulation was designed to study ways of accommodating increased air traffic density. The concepts that were investigated assumed the availability of the microwave landing system and data link and included: (1) multiple curved descending final approaches; (2) parallel runways certified for independent and simultaneous operation under IFR conditions; (3) closer spacing between successive aircraft; and (4) a distributed management system between the air and ground. Three groups each consisting of three pilots and two air traffic controllers flew a combined total of 350 approaches. Piloted simulators were supplied with computer generated traffic situation displays and flight instruments. The controllers were supplied with a terminal area map and digital status information. Pilots and controllers also reported that the distributed management procedure was somewhat more safe and orderly than the centralized management procedure. Flying precision increased as the amount of turn required to intersect the outer mark decreased. Pilots reported that they preferred the alternative of multiple curved descending approaches with wider spacing between aircraft to closer spacing on single, straight in finals while controllers preferred the latter option. Both pilots and controllers felt that parallel runways are an acceptable way to accommodate increased traffic density safely and expeditiously.

Hart, S. G.↗

Simulation of Sweep-Jet Flow Control, Single Jet and Full Vertical Tail

This work is a simulation technology demonstrator, of sweep jet flow control used to suppress boundary layer separation and increase the maximum achievable load coefficients. A sweep jet is a discrete Coanda jet that oscillates in the plane parallel to an aerodynamic surface. It injects mass and momentum in the approximate streamwise direction. It also generates turbulent eddies at the oscillation frequency, which are typically large relative to the scales of boundary layer turbulence, and which augment mixing across the boundary layer to attack flow separation. Simulations of a fluidic oscillator, the sweep jet emerging from a nozzle downstream of the oscillator, and an array of sweep jets which suppresses boundary layer separation are performed. Simulation results are compared to data from a dedicated validation experiment of a single oscillator and its sweep jet, and from a wind tunnel test of a full-scale Boeing 757 vertical tail augmented with an array of sweep jets. A critical step in the work is the development of realistic time-dependent sweep jet inflow boundary conditions, derived from the results of the single-oscillator simulations, which create the sweep jets in the full-tail simulations. Simulations were performed using the computational fluid dynamics (CFD) solver Overow, with high-order spatial discretization and a range of turbulence modeling. Good results were obtained for all flows simulated, when suitable turbulence modeling was used.

flow control↗

Modal analysis for Liapunov stability of rotating elastic bodies

This study consisted of four parallel efforts: (1) modal analyses of elastic continua for Liapunov stability analysis of flexible spacecraft; (2) development of general purpose simulation equations for arbitrary spacecraft; (3) evaluation of alternative mathematical models for elastic components of spacecraft; and (4) examination of the influence of vehicle flexibility on spacecraft attitude control system performance. A complete record is given of achievements under tasks (1) and (3), in the form of technical appendices, and a summary description of progress under tasks two and four.

Colin, A. D.↗

Determination of the geometrical factor and of the spatial distributions of particles for a cosmic ray telescope using a Monte Carlo approach

A program is presented which uses a Monte Carlo technique to compute the geometrical factor of a detector, and the spatial and angular distributions of simulated particles at designated layers within the detector. As the program is written, each layer is assumed to be rectangular, parallel to every other layer with its edges parallel to the edges of every other layer, and centered about a common axis. This program is written in Fortran 4 for the IBM 360 computers and is currently being run under the Fortran 4 (h) compiler. Same features such as the subroutine ERRSET, RANDU, and REMTIM and the type declarations REAL *8 and LOGICAL *1 may need to be changed for use with a different computer system. The use of input data cards, and a sample data set are presented. Lists of the symbol definition, card images of the FORTRAN object deck and subroutines are included along with a flow diagram.

Hemmer, M.↗

Research on spacecraft electrical power conversion

The steady state characteristics and starting behavior of some widely used self-oscillating magnetically coupled square wave inverters were studied and the development of LC-tuned square wave inverters is reported. An analysis on high amplitude voltage spikes which occur in dc-to-square-wave parallel converters shows the importance of various circuit parameters for inverter design and for the suppression of spikes. A computerized simulation of an inductor energy storage dc-to-dc converter with closed loop regulators and of a preregulating current step-up converter are detailed. Work continued on the computer aided design of two-winding energy storage dc-to-dc converters.

Wilson, T. G.↗

Structure of scintillations in Neptune's occultation shadow

An exceptionally high-quality data set from a Neptune occultation is used here to derive a number of new results about the statistical properties of the fluctuations of the intensity distribution in various parts of Neptune's occultation shadow. An approximate numerical ray-tracing model which successfully accounts for many of the qualitative aspects of the observed intensity fluctuation distribution is introduced. Strong refractive scintillation is simulated by including the effects of 'turbulence' with projected atmospheric properties allowed to vary in both the direction perpendicular and parallel to the limb, and an explicit two-dimensional picture of a typical intensity distribution throughout an occulting planet's shadow is presented. The results confirm the existence of highly anisotropic turbulence.

Hubbard, W. B.↗

Convex-profile inversion of asteroid lightcurves - Theory and applications

A theory for asteroid light curve interpretation is derived from the reformulation, extension, and calibration of Ostro and Connely's (1984) convex-profile inversion, and then used to constrain the shapes of selected asteroids. The ideal conditions for the estimation of an asteroid's main cross-section from a light curve are (1) that the viewing-illumination geometry be equatorial, (2) that the scattering be uniform and geometric, (3) that the asteroid surface contours parallel to the equatorial plane be convex, (4) and that the solar angle not equal zero. Results are presented from simulations calibrating the kind, severity, and predictability of systematic errors.

Ostro, Steven J.↗

Monte-Carlo modeling of polar wind photoelectron distributions with anomalous heat flux

In situ measurements above the polar cap by the DE satellites show asymmetric field-aligned electron velocity distributions in the photoelectron energy range 5-60 eV. The formation of these anisotropic distributions is described by a collisional kinetic equation that includes the effects of Coulomb collisions (via a Fokker-Planck operator), those of a parallel electric field, and the magnetic mirror effects of the geomagnetic field. This equation is solved using a Monte Carlo simulation method. The simulation illustrates the formation of the portion of the observed electron distributions originating in the ionosphere and makes it possible to evaluate the significance of the suprathermal electrons at higher altitude. Both observations and calculations indicate that these electrons possess an energy flux which, by its effect on the polar-wind electric field, may influence the dynamics of the polar-wind outflow.

Yasseen, F.↗

Time-accurate Navier-Stokes computations of self-excited two-dimensional unsteady cavity flows

Time-accurate Navier Stokes computations were performed to study self-excited oscillatory shear layer behavior over a two-dimensional rectangular cavity and to gain insight into the vortex-acoustic coupling that sustains these oscillations. Results were obtained solving the full, unsteady compressible Navier-Stokes equations for Reynolds-averaged turbulent air flow past a cavity with L/d = 2. Conditions duplicated test conditions in experiments by Rossiter at Mach numbers varied from 0.6 to 1.2. The computational scheme was the USA 2D series code developed by Rockwell's Science Center with an algebraic turbulence model for the shear layer and recirculating region. Two physical conditions were simulated: an isolated cavity situated (1) in an infinite plane wall and (2) in an infinite channel with an opposing parallel wall. Computed results for (1) showed excellent agreement with Rossiter's measured amplitudes and frequencies. A tuned acoustic duct resonance effect was observed in (2), demonstrating a controlled change in the acoustic feedback mechanism and resultant amplification in computed oscillations.

Dougherty, N. Sam↗

ROM/RMF control of the flexible, articulated-truss space crane

The large size and extreme flexibility of an articulated space crane entail great difficulties in precision positioning. Control strategies for the flexible multibody system are here developed using a Timoshenko-beam finite-element model (FEM) which captures the system's low frequency characteristics. From the FEM, a reduced-order model (ROM) controller is designed which rotates the fixed shape crane about the base alpha joint to actuate positioning control. Eigenvalue analysis of the closed-loop system shows the destabilizing effect of controller-structure interaction (CSI). Making use of a parallel bank of residual mode filters (RMFs) compensated for CSI, restoring system stability and performance. The linear control law is simulated in closed-loop with the beam model. The ROM-based control law with RMF compensation is an effective, low-order controller.

Reisenauer, Brian T.↗

Job Management Requirements for NAS Parallel Systems and Clusters

A job management system is a critical component of a production supercomputing environment, permitting oversubscribed resources to be shared fairly and efficiently. Job management systems that were originally designed for traditional vector supercomputers are not appropriate for the distributed-memory parallel supercomputers that are becoming increasingly important in the high performance computing industry. Newer job management systems offer new functionality but do not solve fundamental problems. We address some of the main issues in resource allocation and job scheduling we have encountered on two parallel computers - a 160-node IBM SP2 and a cluster of 20 high performance workstations located at the Numerical Aerodynamic Simulation facility. We describe the requirements for resource allocation and job management that are necessary to provide a production supercomputing environment on these machines, prioritizing according to difficulty and importance, and advocating a return to fundamental issues.

Saphir, William↗

Multiple Satellite Observations of High-Latitude Ionospheric Outflows

We will present reported observations of up- and down-flows of topside ionospheric thermal plasmas from multiple near-simultaneous tracks through the high-latitude topside ionosphere. From several Southern polar passes, it has been possible to construct plots of field-aligned flows of O+ observed by the Thermal Ion Dynamics Experiment(TIDE) on the POLAR spacecraft near 5000 km altitude together with vertical ion flow observations from one or more DMSP spacecraft near 800 km altitude. We will also involve simulations from our Dynamic Fluid-Kinetic(DyFK) modeling of polar plasma transport and display resulting altitude profiles of ion parallel velocities and densities, and overlay those profiles with the "conjunction" measurements by DMSP(800 km) and POLAR(approx. 5000 km altitude). We also will present simultaneous observations of POLAR auroral UVI images with field-aligned flows.

Horwitz, J.↗

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.↗