Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed parallelization”

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,135 records · Page 63

Whistler mode noise in Jupiter's inner magnetosphere

A study is made of the amplitude and spectral extent of whistler mode noise in the inner magnetosphere of Jupiter. It is found that the 'hat-shaped' pitch angle distributions of energetic electrons (21 and 31 MeV at L=3) are consistent with those predicted in the presence of a band-limited spectrum of whistler mode noise. The equatorial maximum linear growth rate of parallel propagating whistlers are consistent with those necessary to limit the energetic electron intensities by the whistler mode instability. It is noted that the wave phase speeds before wave reflection can occur at high latitudes and that wave growth is limited to a disk-like region centered around the magnetic equator. The frequency extent of the whistler mode noise spectrum may be estimated by the range of frequencies maximally unstable to equatorial linear growth. A value is found for the spectral density of the broadband whistler mode noise necessary to balance radial diffusion of energetic electrons above the critical range, and an expression is derived for the energetic electron system response to fluctuations about the limiting flux value.

Sentman, D. D.↗

Electromagnetic ion beam instabilities - Hot beams at interplanetary shocks

This paper considers Vlasov instabilities driven by a very hot ion beam streaming parallel to a magnetic field. The linear theory of electromagnetic instabilities driven by such a beam in a homogeneous, collisionless, nonrelativistic plasma is used. Numerical solutions of the full linear dispersion equation for bi-Maxwellian distribution functions are presented. If the thermal speed of the beam is greater than its drift speed, the two dominant modes are the right-hand and left-hand resonant ion beam instabilities. The parametric dependencies of the threshold drift speeds and growth rates for both modes are presented. In particular, the thresholds are shown to be sensitive functions of the beam temperature anisotropies. Arguments are presented that these two modes are driven by the suprathermal ion component at interplanetary shocks, and that the growth of these modes is sufficient to account for the amplitudes of MHD-like waves observed at such shocks.

Gary, S. P.↗

Simulating a small turboshaft engine in real-time multiprocessor simulator (RTMPS) environment

A Real-Time Multiprocessor Simulator (RTMPS) has been developed at NASA Lewis Research Center. The RTMPS uses parallel microprocessors to achieve computing speeds needed for real-time engine simulation. This report describes the use of the RTMPS system to simulate a small turboshaft engine. The process of programming the engine equations and distributing them over one, two, and four processors is discussed. Steady-state and transient results from the RTMPS simulation are compared with results from a main-frame-based simulation. Processor execution times and the associated execution time savings for the two and four processor cases are presented using actual data obtained from the RTMPS system. Included is a discussion of why the minimum achievable calculation time for the turboshaft engine model was attained using four processors. Finally, future enhancements to the RTMPS system are discussed including the development of a generalized partitioning algorithm to automatically distribute the system equations among the processors in optimum fashion.

Milner, E. J.↗

Software For Drawing Design Details Concurrently

Software system containing five computer-aided-design programs enables more than one designer to work on same part or assembly at same time. Reduces time necessary to produce design by implementing concept of parallel or concurrent detailing, in which all detail drawings documenting three-dimensional model of part or assembly produced simultaneously, rather than sequentially. Keeps various detail drawings consistent with each other and with overall design by distributing changes in each detail to all other affected details.

Crosby, Dewey C., III↗

CUSP Energetic Particles: Confinement, Acceleration and Implications

The cusp energetic particle (CEP) event is a new magnetospheric phenomenon. The events were detected in the dayside cusp for hours, in which the measured helium ions had energies up to 8 MeV. All of these events were associated with a dramatic decrease and large fluctuations in the local magnetic field strength. During January 1999 - December 1999 covered by this report, I have studied the CEP events by analyzing the POLAR, GEOTAIL, and WIND particle and magnetic field data measured during the geomagnetic quiet periods in 1996 and one geomagnetic storm period in 1998. The simultaneous observations indicated that the ion fluxes in the CEP events were higher than that in both the upstream and the downstream from the bow shock. The pitch angle distribution of the helium ions in the CEP events was found to peak around 90 deg. It was found that the mirror parameter, defined as the ratio of the square root of the integration of the parallel turbulent power spectral component over the ultra-low frequency (ULF) ranges to the mean field in the cusp, is correlated with the intensity of the cusp MeV helium flux, which is a measure of the influence of mirroring interactions and an indication of local effect. It was also found that the turbulent power of the local magnetic field in the ultra-low frequency (ULF) ranges is correlated with the intensity of the cusp energetic helium ions. Such ULF ranges correspond to periods of about 0.33-500 seconds that cover the gyroperiods, the bounce periods, and the drift periods of the tens keV to MeV charged particles when they are temporarily confined in the high-altitude dayside cusp. These observations represent a discovery that the high-altitude dayside cusp is a new acceleration and dynamic trapping region of the magnetosphere. The cusp geometry is connected via gradient and curvature drift of these energized ions to the equatorial plasma sheet as close as the geostationary orbit at local midnight. It implies that the dayside cusp is potentially an important source of magnetospheric particles. The discovery of the CEP events has been recognized as one of the most significant results from POLAR. I was invited to give a talk at 1999 IUGG meeting to interpret the CEP events. This discovery has also been written into the web-based Space Physics Text Book (http://www.oulu.fi/- spaceweb/textbook/cusp.html).

Chen, Jiasheng↗

Efficacy of Code Optimization on Cache-based Processors

The current common wisdom in the U.S. is that the powerful, cost-effective supercomputers of tomorrow will be based on commodity (RISC) micro-processors with cache memories. Already, most distributed systems in the world use such hardware as building blocks. This shift away from vector supercomputers and towards cache-based systems has brought about a change in programming paradigm, even when ignoring issues of parallelism. Vector machines require inner-loop independence and regular, non-pathological memory strides (usually this means: non-power-of-two strides) to allow efficient vectorization of array operations. Cache-based systems require spatial and temporal locality of data, so that data once read from main memory and stored in high-speed cache memory is used optimally before being written back to main memory. This means that the most cache-friendly array operations are those that feature zero or unit stride, so that each unit of data read from main memory (a cache line) contains information for the next iteration in the loop. Moreover, loops ought to be 'fat', meaning that as many operations as possible are performed on cache data-provided instruction caches do not overflow and enough registers are available. If unit stride is not possible, for example because of some data dependency, then care must be taken to avoid pathological strides, just ads on vector computers. For cache-based systems the issues are more complex, due to the effects of associativity and of non-unit block (cache line) size. But there is more to the story. Most modern micro-processors are superscalar, which means that they can issue several (arithmetic) instructions per clock cycle, provided that there are enough independent instructions in the loop body. This is another argument for providing fat loop bodies. With these restrictions, it appears fairly straightforward to produce code that will run efficiently on any cache-based system. It can be argued that although some of the important computational algorithms employed at NASA Ames require different programming styles on vector machines and cache-based machines, respectively, neither architecture class appeared to be favored by particular algorithms in principle. Practice tells us that the situation is more complicated. This report presents observations and some analysis of performance tuning for cache-based systems. We point out several counterintuitive results that serve as a cautionary reminder that memory accesses are not the only factors that determine performance, and that within the class of cache-based systems, significant differences exist.

VanderWijngaart, Rob F.↗

Radiative Heat Transfer Capability Implemented in OpenNCC for Conjugate Heat Transfer Applications

Thermal efficiency of gas turbine engine increases as the temperature and pressure at the combustor increases. Consequently, the materials used inside a combustor must survive an increasingly challenging environment. For this reason, accurate assessment of heat transfer is crucial for combustor design. While all three modes of heat transfer are present inside a combustor, the focus of this paper is the thermal radiation. Radiative heat transfer in a gas turbine combustors are particularly interesting from three reasons. Firstly, the radiative heat loss from the combustion region may affect the emission performance. Secondly, the cooling air will protect the liner from convection but not necessary from radiation. Finally, it is less frequently incorporated in CFD analysis than other forms of heat transfer. In this work, radiative heat transfer using discrete ordinate method has been incorporated in OpenNCC (a publicly releasable version of the National Combustion Code) developed at NASA Glenn Research Center. Aside from massively parallel computation capability using MPI and the ability to utilize unstructured mesh, the current implementation includes two types of spectral models, namely, the weighted some of gray gas model and the full spectrum correlated k-distribution model. After presenting the theory and the strategy of implementation, results of validation cases for gray gas and spectral models will be presented. While the implementation of the radiation solver is intended for gas turbine application, the radiation solver can run independently from the convection/combustion solver and the same theory can be applied to other application.

OpenNCC↗

Radiative Heat Transfer Capability Implemented in OpenNCC for Conjugate Heat Transfer Applications

Thermal efficiency of gas turbine engine increases as the temperature and pressure at the combustor increases. Consequently, the materials used inside a combustor must survive an increasingly challenging environment. For this reason, accurate assessment of heat transfer is crucial for combustor design. While all three modes of heat transfer are present inside a combustor, the focus of this paper is the thermal radiation. Radiative heat transfer in a gas turbine combustors are particularly interesting from three reasons. Firstly, the radiative heat loss from the combustion region may affect the emission performance. Secondly, the cooling air will protect the liner from convection but not necessary from radiation. Finally, it is less frequently incorporated in CFD analysis than other forms of heat transfer. In this work, radiative heat transfer using discrete ordinate method has been incorporated in OpenNCC (a publicly releasable version of the National Combustion Code) developed at NASA Glenn Research Center. Aside from massively parallel computation capability using MPI and the ability to utilize unstructured mesh, the current implementation includes two types of spectral models, namely, the weighted some of gray gas model and the full spectrum correlated k-distribution model. After presenting the theory and the strategy of implementation, results of validation cases for gray gas and spectral models will be presented. While the implementation of the radiation solver is intended for gas turbine application, the radiation solver can run independently from the convection/combustion solver and the same theory can be applied to other application.

OpenNCC↗

Online assessment of a distributed processor

ORT (Operational Readiness Test) software allows one engineer to test readiness of 64 minicomputers and their peripherals from single console. Software makes roll call of computers and peripherals via common data buffer to check readiness of system in morning "wake up" or at other important times. Subsystems are tested in parallel to save time. "Watchdog" terminates test of any system that does not respond in time, so one failed system does not halt test sequence. Entire rollcall is complete in about 15 minutes. Software is designed for Space Shuttle prelaunch checkout, but approach should interest users of similar equipment.

Ehrlich, L. F.↗

Physiological asymmetry in etiolated pea epicotyls: relation to patterns of auxin distribution and phototropic behavior

Etiolated pea seedlings require transformation of Pr phytochrome to Pfr before they display optimal phototropic response to unilateral blue light. This study investigates the possible role of auxin transport in explaining these phenomena. Labeled [2-14C]IAA applied to the intact terminal buds of dark-grown and red light-treated pea seedlings was measured 210 min later on the shaded and illuminated sides of the epicotyl as a function of direction and duration of irradiation with blue light. Totally darkened epicotyls show an asymmetry in distribution of radioactivity in the upper growth zone of the epicotyl, in favor of the side under the concave part of the apical hook. Red light, which greatly potentiates curvature toward subsequent unilateral blue light, lowers this asymmetry. Blue light directed to the epicotyl of red-pretreated plants in a plane parallel to the hook and from the side bearing the convex portion of the hook induces positive phototropic curvature as well as a surplus of radioactivity on the illuminated side of the upper epicotyl and on the shaded side of the lower growth zone of the epicotyl. Light directed to the side bearing the concave part of the hook also causes an accumulation of counts in the upper part of the lighted side but produces neither curvature of the epicotyl nor accumulation of counts in the lower shaded side. Because of this built-in physiological asymmetry in the growth zone just below the apical hook, it is difficult to explain the effects of red and blue light on curvature in terms of patterns of auxin distribution alone.

NASA Program Space Biology↗

Unsteady supersonic aerodynamic theory by the method of potential gradient

A generalized solution of the hyperbolic wave equation has been derived. The method used has been further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals have been evaluated in closed form, while numerical integration methods are suggested for treating more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient because line integrations along the other two sides need only be considered. Furthermore, all the integrands vanish on the Mach cone and the need for determining the hyperbolic curves of intersection of the cone with the lifting surface is avoided. Generalized aerodynamic coefficients for three AGARD planforms have been calculated and compared with the available results.

Jones, W. P.↗

A Comparison of the Chordwise Pressure Distribution and Spanwise Distribution of Loading at Subsonic Speeds on Two Triangular Wings of Aspect Ratio 2 Having NACA 0005 and 0008 Sections

Measurements of the distribution of surface pressures were made on two wing-body combinations employing a triangular wing having an aspect ratio ratio of 2. The first of these combinations had a wing with the NACA 0008-63 section parallel to the air stream while the wing of the second had the NACA 0005-63 section in the stream direction. The measurements were made for Mach numbers from 0.11 to 0.95 at a constant Reynolds number of 3.0 million. Data were also obtained for Reynolds numbers up to 15.0 million at the lower Mach number. A summary of data obtained from previously published force tests of these same models is also included for completeness of the presentation of the effect of wing thickness.

Smith, Donald W↗

NAS Experiences of Porting CM Fortran Codes to HPF on IBM SP2 and SGI Power Challenge

Current Connection Machine (CM) Fortran codes developed for the CM-2 and the CM-5 represent an important class of parallel applications. Several users have employed CM Fortran codes in production mode on the CM-2 and the CM-5 for the last five to six years, constituting a heavy investment in terms of cost and time. With Thinking Machines Corporation's decision to withdraw from the hardware business and with the decommissioning of many CM-2 and CM-5 machines, the best way to protect the substantial investment in CM Fortran codes is to port the codes to High Performance Fortran (HPF) on highly parallel systems. HPF is very similar to CM Fortran and thus represents a natural transition. Conversion issues involved in porting CM Fortran codes on the CM-5 to HPF are presented. In particular, the differences between data distribution directives and the CM Fortran Utility Routines Library, as well as the equivalent functionality in the HPF Library are discussed. Several CM Fortran codes (Cannon algorithm for matrix-matrix multiplication, Linear solver Ax=b, 1-D convolution for 2-D datasets, Laplace's Equation solver, and Direct Simulation Monte Carlo (DSMC) codes have been ported to Subset HPF on the IBM SP2 and the SGI Power Challenge. Speedup ratios versus number of processors for the Linear solver and DSMC code are presented.

Saini, Subhash↗

Orchestrator Telemetry Processing Pipeline

Orchestrator is a software application infrastructure for telemetry monitoring, logging, processing, and distribution. The architecture has been applied to support operations of a variety of planetary rovers. Built in Java with the Eclipse Rich Client Platform, Orchestrator can run on most commonly used operating systems. The pipeline supports configurable parallel processing that can significantly reduce the time needed to process a large volume of data products. Processors in the pipeline implement a simple Java interface and declare their required input from upstream processors. Orchestrator is programmatically constructed by specifying a list of Java processor classes that are initiated at runtime to form the pipeline. Input dependencies are checked at runtime. Fault tolerance can be configured to attempt continuation of processing in the event of an error or failed input dependency if possible, or to abort further processing when an error is detected. This innovation also provides support for Java Message Service broadcasts of telemetry objects to clients and provides a file system and relational database logging of telemetry. Orchestrator supports remote monitoring and control of the pipeline using browser-based JMX controls and provides several integration paths for pre-compiled legacy data processors. At the time of this reporting, the Orchestrator architecture has been used by four NASA customers to build telemetry pipelines to support field operations. Example applications include high-volume stereo image capture and processing, simultaneous data monitoring and logging from multiple vehicles. Example telemetry processors used in field test operations support include vehicle position, attitude, articulation, GPS location, power, and stereo images.

Powell, Mark↗

Analysis of composite ablators using massively parallel computation

In this work, the feasibility of using massively parallel computation to study the response of ablative materials is investigated. Explicit and implicit finite difference methods are used on a massively parallel computer, the Thinking Machines CM-5. The governing equations are a set of nonlinear partial differential equations. The governing equations are developed for three sample problems: (1) transpiration cooling, (2) ablative composite plate, and (3) restrained thermal growth testing. The transpiration cooling problem is solved using a solution scheme based solely on the explicit finite difference method. The results are compared with available analytical steady-state through-thickness temperature and pressure distributions and good agreement between the numerical and analytical solutions is found. It is also found that a solution scheme based on the explicit finite difference method has the following advantages: incorporates complex physics easily, results in a simple algorithm, and is easily parallelizable. However, a solution scheme of this kind needs very small time steps to maintain stability. A solution scheme based on the implicit finite difference method has the advantage that it does not require very small times steps to maintain stability. However, this kind of solution scheme has the disadvantages that complex physics cannot be easily incorporated into the algorithm and that the solution scheme is difficult to parallelize. A hybrid solution scheme is then developed to combine the strengths of the explicit and implicit finite difference methods and minimize their weaknesses. This is achieved by identifying the critical time scale associated with the governing equations and applying the appropriate finite difference method according to this critical time scale. The hybrid solution scheme is then applied to the ablative composite plate and restrained thermal growth problems. The gas storage term is included in the explicit pressure calculation of both problems. Results from ablative composite plate problems are compared with previous numerical results which did not include the gas storage term. It is found that the through-thickness temperature distribution is not affected much by the gas storage term. However, the through-thickness pressure and stress distributions, and the extent of chemical reactions are different from the previous numerical results. Two types of chemical reaction models are used in the restrained thermal growth testing problem: (1) pressure-independent Arrhenius type rate equations and (2) pressure-dependent Arrhenius type rate equations. The numerical results are compared to experimental results and the pressure-dependent model is able to capture the trend better than the pressure-independent one. Finally, a performance study is done on the hybrid algorithm using the ablative composite plate problem. It is found that there is a good speedup of performance on the CM-5. For 32 CPU's, the speedup of performance is 20. The efficiency of the algorithm is found to be a function of the size and execution time of a given problem and the effective parallelization of the algorithm. It also seems that there is an optimum number of CPU's to use for a given problem.

Shia, David↗

Ion sheath effects on RF plasma probes - Experimental results in laboratory plasmas

Laboratory work dealing with the frequency characteristic of the plasma impedance of spherical and cylindrical electrode systems is reported. The influence of the ion sheath on various features of the impedance characteristic is emphasized. Those features are the series and parallel resonance as well as additional resonances due to the excitation of electroacoustic and cyclotron harmonic waves. The dependence of the series and parallel resonances on dc biasing leads to a method of determining the ion sheath capacity for a cylindrical electrode system. The obtained values agree fairly well with those obtained from a theoretical model for the density and potential distribution in the sheath of a cylindrical sensor aligned with a supersonic plasma flow. The amplitude of resonances due to excitation of longitudinal plasma waves (electroacoustic and cyclotron harmonic) is reduced or even vanishes for sufficiently negative dc bias. Positive bias first leads to an increased amplitude up to a certain dc bias value above which, however, the amplitude decreases again due to electron absorption at the sensor surface.

Kist, R.↗

Infrared radiation emerging from smoke produced by brush fires

The IR radiative transport properties of brush fire smoke clouds, computed for a model with finite horizontal dimensions as well as the more common plane-parallel model, are presented. The finite model is a three-dimensional version of the two-stream approximation applied to cubic clouds of steam, carbon, and silicates. Assumptions are made with regard to the shape and size distributions of the smoke particles. It is shown that 11.5-micron radiometry can detect fires beneath smoke clouds if the path integrated mass density of the smoke is less than or equal to 3 g/sq m.

Weinman, J. A.↗

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