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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Execution time supports for adaptive scientific algorithms on distributed memory machines

Optimizations are considered that are required for efficient execution of code segments that consists of loops over distributed data structures. The PARTI (Parallel Automated Runtime Toolkit at ICASE) execution time primitives are designed to carry out these optimizations and can be used to implement a wide range of scientific algorithms on distributed memory machines. These primitives allow the user to control array mappings in a way that gives an appearance of shared memory. Computations can be based on a global index set. Primitives are used to carry out gather and scatter operations on distributed arrays. Communications patterns are derived at runtime, and the appropriate send and receive messages are automatically generated.

Berryman, Harry↗
NASA NTRS1990

Execution time support for scientific programs on distributed memory machines

Optimizations are considered that are required for efficient execution of code segments that consists of loops over distributed data structures. The PARTI (Parallel Automated Runtime Toolkit at ICASE) execution time primitives are designed to carry out these optimizations and can be used to implement a wide range of scientific algorithms on distributed memory machines. These primitives allow the user to control array mappings in a way that gives an appearance of shared memory. Computations can be based on a global index set. Primitives are used to carry out gather and scatter operations on distributed arrays. Communications patterns are derived at runtime, and the appropriate send and receive messages are automatically generated.

Berryman, Harry↗
NASA NTRS1987

Mapping a battlefield simulation onto message-passing parallel architectures

Perhaps the most critical problem in distributed simulation is that of mapping: without an effective mapping of workload to processors the speedup potential of parallel processing cannot be realized. Mapping a simulation onto a message-passing architecture is especially difficult when the computational workload dynamically changes as a function of time and space; this is exactly the situation faced by battlefield simulations. This paper studies an approach where the simulated battlefield domain is first partitioned into many regions of equal size; typically there are more regions than processors. The regions are then assigned to processors; a processor is responsible for performing all simulation activity associated with the regions. The assignment algorithm is quite simple and attempts to balance load by exploiting locality of workload intensity. The performance of this technique is studied on a simple battlefield simulation implemented on the Flex/32 multiprocessor. Measurements show that the proposed method achieves reasonable processor efficiencies. Furthermore, the method shows promise for use in dynamic remapping of the simulation.

Nicol, David M.↗
NASA NTRS2002

Wave-Particle Interactions Associated with Nongyrotropic Distribution Functions: A Hybrid Simulation Study

Nongyrotropic plasma distribution functions can be formed in regions of space where guiding center motion breaks down as a result of strongly curved and weak ambient magnetic fields. Such are the conditions near the current sheet in the Earth's middle and distant magnetotail, where observations of nongyrotropic ion distributions have been made. Here a systematic parameter study of nongyrotropic proton distributions using electromagnetic hybrid simulations is made. We model the observed nongyrotropic distributions by removing a number of arc length segments from a cold ring distribution and find significant differences with the results of simulations that initially have a gyrotropic ring distribution. Model nongyrotropic distributions with initially small perpendicular thermalization produce growing fluctuations that diffuse the ions into a stable Maxwellian-like distribution within a few proton gyro periods. The growing waves produced by nongyrotropic distributions are similar to the electromagnetic proton cyclotron waves produced by a gyrotropic proton ring distribution in that they propagate parallel to the background magnetic field and occur at frequencies on the order of the proton gyrofrequency, The maximum energy of the fluctuating magnetic field increases as the initial proton distribution is made more nongyrotropic, that is, more highly bunched in perpendicular velocity space. This increase can be as much as twice the energy produced in the gyrotropic case.

Convery, P. D.↗
NASA NTRS2002

Observations of Closed Structures at the Magnetopause: A Case for Multiple Reconnections

We further analyze a case of Interball LLBL crossing on the dusk flank of geomagnetosphere under southward magnetosheath magnetic field, previously categorized as an interval of highly structured LLBL. These conditions of highly structured LLBL include reconnection signatures. Observed ion velocity distributions with LLBL are quite variable. D-shaped distributions that are associated with the open reconnected flux tube are observed at the boundaries of LLBL transients and sometimes within the LLBL transients. In most cases the ion velocity distributions consist of two magnetosheath-type components with different velocities parallel to the magnetic field, or of three components one of which has nearly zero Vpar. The shapes of ion velocity distributions and their evolution with decreasing number density in LLBL indicate that most of LLBL is located on closed magnetic field lines. These observations strongly favor multiple reconnections between magnetosheath and magnetosphereric flux tubes, creating long spiral flux tube islands at the magnetopause. We report evidence for the simultaneous occurrence of magnetic reconnection at multiple points across the magnetopause, as has been proposed and found to occur in magnetopause simulations. The evidence is in the form of highly structured distributions of ions in velocity parallel to the local magnetic field direction, within the magnetopause and low latitude boundary layer region, from the Interball-Tall spacecraft. We interpret these distributions as a natural consequence of the formation of spiral magnetic flux tubes consisting of a mixture of alternating segments originating from the magnetosheath or interplanetary plasma and from the low latitude boundary layer or magnetospheric plasma. We further analyze a case of Interball LLBL crossing on the dusk flank of geomagnetosphere under southward magnetosheath magnetic field, previously categorized as an interval of highly structured LLBL. These conditions of highly structured LLBL include reconnection signatures. Observed ion velocity distributions with LLBL are quite variable. D-shaped distributions that are associated with the open reconnected flux tube are observed at the boundaries of LLBL transients and sometimes within the LLBL transients. In most cases the ion velocity distributions consist of two magnetosheath-type components with different velocities parallel to the magnetic field, or of three components one of which has nearly zero Vpar. The shapes of ion velocity distributions and their evolution with decreasing number density in LLBL indicate that most of LLBL is located on closed magnetic field lines. These observations strongly favor multiple reconnection between magnetosheath and magnetospheric flux tubes, creating long spiral flux tube islands at the magnetopause. We report evidence for the simultaneous occurrence of magnetic reconnection at multiple points across the magnetopause, as has been proposed and found to occur in magnetopause simulations. The evidence is in the form of highly structured distributions of ions in velocity parallel to the local magnetic field direction, within the magnetopause and low latitude boundary layer region, from the Interball-Tail spacecraft. We interpret these distributions as a natural consequence of the formation of spiral. We further analyze a case of Interball LLBL crossing on the dusk flank of geomagnetosphere under southward magnetosheath magnetic field, previously categorized as an interval of highly structured LLBL. These conditions of highly structured LLBL include reconnection signatures. Observed ion velocity distributions with LLBL are quite variable. D-shaped distributions that are associated with the open reconnected flux tube are observed at the boundaries of LLBL transients and sometimes within the LLBL transients. In most cases the ion velocity distributions consist of two magnetosheath-type components with different velocities parallel to the magnetic field, or of three components one of which has nearly zero Vpar.

Vaisberg, O. L.↗
NASA NTRS1987

Molecular dynamics of the water liquid-vapor interface

The results of molecular dynamics calculations on the equilibrium interface between liquid water and its vapor at 325 K are presented. For the TIP4P model of water intermolecular pair potentials, the average surface dipole density points from the vapor to the liquid. The most common orientations of water molecules have the C2 nu molecular axis roughly parallel to the interface. The distributions are quite broad and therefore compatible with the intermolecular correlations characteristic of bulk liquid water. All near-neighbor pairs in the outermost interfacial layers are hydrogen bonded according to the common definition adopted here. The orientational preferences of water molecules near a free surface differ from those near rigidly planar walls which can be interpreted in terms of patterns found in hexagonal ice 1. The mean electric field in the interfacial region is parallel to the mean polarization which indicates that attention cannot be limited to dipolar charge distributions in macroscopic descriptions of the electrical properties of this interface. The value of the surface tension obtained is 132 +/- 46 dyn/cm, significantly different from the value for experimental water of 68 dyn/cm at 325 K.

NASA Center ARC↗
NASA NTRS2000

Towards Flange-to-Flange Turbopump Simulations for Liquid Rocket Engines

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbopump geometry through numerical simulation will be of significant value toward design. This will help to improve safety of future space missions. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the MPI and MLP versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology. INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbopump simulations, moving boundary capability and an efficient time-accurate integration methods are build in the flow solver. To handle the geometric complexity and moving boundary problems, overset grid scheme is incorporated with the solver that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. The performance of the two time integration schemes for time-accurate computations is investigated. For an unsteady flow which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive. The current geometry for the LOX boost turbopump has various rotating and stationary components, such as inducer, stators, kicker, hydrolic turbine, where the flow is extremely unsteady. Figure 1 shows the geometry and computed surface pressure of the inducer. The inducer and the hydrolic turbine rotate in different rotational speed.

Kiris, Cetin↗
NASA NTRS1993

A parallel dynamic load balancing algorithm for 3-D adaptive unstructured grids

Adaptive local grid refinement and coarsening results in unequal distribution of workload among the processors of a parallel system. A novel method for balancing the load in cases of dynamically changing tetrahedral grids is developed. The approach employs local exchange of cells among processors in order to redistribute the load equally. An important part of the load balancing algorithm is the method employed by a processor to determine which cells within its subdomain are to be exchanged. Two such methods are presented and compared. The strategy for load balancing is based on the Divide-and-Conquer approach which leads to an efficient parallel algorithm. This method is implemented on a distributed-memory MIMD system.

Vidwans, A.↗
NASA NTRS1993

Electron velocity space hole modes

In space plasmas the local electron distribution is formed by ballistic motion from different source regions. Since electrons with positive and negative parallel (v-parallel) velocities have different sources and experience different interactions with the large-scale electromagnetic fields, the electron distribution function near v-parallel = 0 should possess considerable structure, including the possibility of velocity space holes, i.e., regions where the phase space density is zero, and sharp discontinuities. A particular hole distribution, which may occur upstream of slow shocks in the distant geomagnetic tail, is shown to support a nonstandard plasma wave mode with frequencies between the ion and electron plasma frequencies and parallel phase speed above the ion, but well below the electron, thermal speed. The hole mode can be destabilized by a positive slope in the low parallel velocity electron distribution, and the growth rate is strongly enhanced if the electron parallel temperature exceeds the perpendicular temperature; the mode can also be excited by a hot, fast ion beam. The hole modes exist only for a narrow range of propagation angles about the magnetic field direction.

Coroniti, Ferdinand V.↗
NASA NTRS1990

Balancing Loads Among Parallel Data Processors

Heuristic algorithm minimizes amount of memory used by multiprocessor system. Distributes load of many identical, short computations among multiple parallel digital data processors, each of which has its own (local) memory. Each processor operates on distinct and independent set of data in larger shared memory. As integral part of load-balancing scheme, total amount of space used in shared memory minimized. Possible applications include artificial neural networks or image processors for which "pipeline" and vector methods of load balancing inappropriate.

Baffes, Paul Thomas↗
NASA NTRS1995

Parallel processing for scientific computations

The scope of this project dealt with the investigation of the requirements to support distributed computing of scientific computations over a cluster of cooperative workstations. Various experiments on computations for the solution of simultaneous linear equations were performed in the early phase of the project to gain experience in the general nature and requirements of scientific applications. A specification of a distributed integrated computing environment, DICE, based on a distributed shared memory communication paradigm has been developed and evaluated. The distributed shared memory model facilitates porting existing parallel algorithms that have been designed for shared memory multiprocessor systems to the new environment. The potential of this new environment is to provide supercomputing capability through the utilization of the aggregate power of workstations cooperating in a cluster interconnected via a local area network. Workstations, generally, do not have the computing power to tackle complex scientific applications, making them primarily useful for visualization, data reduction, and filtering as far as complex scientific applications are concerned. There is a tremendous amount of computing power that is left unused in a network of workstations. Very often a workstation is simply sitting idle on a desk. A set of tools can be developed to take advantage of this potential computing power to create a platform suitable for large scientific computations. The integration of several workstations into a logical cluster of distributed, cooperative, computing stations presents an alternative to shared memory multiprocessor systems. In this project we designed and evaluated such a system.

Alkhatib, Hasan S.↗
NASA NTRS1991

Pitch angle scattering of cometary ions into monospherical and bispherical distributions

Low frequency magnetic fluctuations generated by solar wind/cometary ion distributions cause pitch-angle scattering of cometary ions. The nature of this pitch angle scattering is examined by means of one-dimensional hybrid simulations, in which newborn ions are created at a constant rate, and for various injection angles, alpha, relative to the magnetic field. Pitch angle scattering in the quasi-perpendicular regime (alpha above 60 deg) results in a bispherical velocity distribution while in the quasi-parallel regime (alpha between 0 and 60 deg), pitch angle scattering can result in either a monospherical or bispherical distribution, depending on the injection angle and initial beam velocity.

Miller, Ronald H.↗
DOE OSTI2024

Distributed Augmentation, Hypersweeps, and Branch Decomposition of Contour Trees for Scientific Exploration

Contour trees describe the topology of level sets in scalar fields and are widely used in topological data analysis and visualization. A main challenge of utilizing contour trees for large-scale scientific data is their computation at scale using highperformance computing. To address this challenge, recent work has introduced distributed hierarchical contour trees for distributed computation and storage of contour trees. However, effective use of these distributed structures in analysis and visualization requires subsequent computation of geometric properties and branch decomposition to support contour extraction and exploration. In this work, we introduce distributed algorithms for augmentation, hypersweeps, and branch decomposition that enable parallel computation of geometric properties, and support the use of distributed contour trees as query structures for scientific exploration. Finally, we evaluate the parallel performance of these algorithms and apply them to identify and extract important contours for scientific visualization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Parallel Variable Population Multi-Objective Optimizer (pvpmoo) v1.0

This is a parallel variable population multi-objective optimizer with an adaptive unified differential evolution algorithm or a genetic algorithm. It can also be used for single objective optimization. Some features of this code include: 1) The population size varies from generation to generation to save the total # of objective function evaluations. 2) The population is uniformly distributed to a number of parallel processors for simultaneous objective function evaluation. 3) The objective function evaluation can be attained from an external simulation program with control variables in its input file and objectives calculated from its output files. 4) The optimizer includes an adaptive unified differential evolution algorithm and a real value genetic algorithm. The parameters in the unified differential evolution algorithm can be chosen to attain any mutation schemes in the published literature.

Qiang, Ji↗
NASA NTRS1991

A Parallel Rendering Algorithm for MIMD Architectures

Applications such as animation and scientific visualization demand high performance rendering of complex three dimensional scenes. To deliver the necessary rendering rates, highly parallel hardware architectures are required. The challenge is then to design algorithms and software which effectively use the hardware parallelism. A rendering algorithm targeted to distributed memory MIMD architectures is described. For maximum performance, the algorithm exploits both object-level and pixel-level parallelism. The behavior of the algorithm is examined both analytically and experimentally. Its performance for large numbers of processors is found to be limited primarily by communication overheads. An experimental implementation for the Intel iPSC/860 shows increasing performance from 1 to 128 processors across a wide range of scene complexities. It is shown that minimal modifications to the algorithm will adapt it for use on shared memory architectures as well.

Crockett, Thomas W.↗
NASA NTRS2005

In Vitro Experimental Model to Investigate the Biological Effects across the Bragg Curve of High-LET Radiation

The space environment consists of a varying field of radiation particles including high energy ions, with a spacecrafts shielding material providing the only major protection to astronauts from harmful exposure. Unlike lowLET gamma or Xrays, the presence of shielding does not always reduce the radiation risks for energetic charged particle exposure since the dose delivered by the charged particle increases sharply as the particle approaches the end of its range, a position known as the Bragg peak and the correlating spatial dose distribution identified as the Bragg curve. The Bragg curve does not necessarily represent the biological damage along the particle traversal since biological effects are influenced by the track structure of both primary and secondary particles. Therefore, the biological Bragg curve is dependent of the energy and the type of the primary particle, and may vary for different biological endpoints. Here we describe a unique irradiation geometry and experimental system to measure the biological response across the Bragg curve in one consistent biological sample. Polyethylene shielding was used to achieve a Bragg curve distribution with the beam geometry parallel to a monolayer of fibroblast cells. We present data that highlights the differential formation of DNA double strand breaks (DSBs) and chromosomal deletions across the Bragg curve in human fibroblasts irradiated with 600 MeV/nucleon iron ion beams. Qualitative analyses of gammaH2AX fluorescence, a known marker of DSBs, indicated potentially increased clustering of DNA damage before the Bragg peak, enhanced homogenous distribution at the peak, and provided visual evidence of high linear energy transfer (LET) particle traversal of cells beyond the Bragg peak in agreement with one-dimensional transport approximations. A biological response curve generated for micronuclei induction across the Bragg curve for 600 MeV/n Fe ions did not reveal an increase in the yield of micronuclei at the Bragg peak location. Assessment of such biological parameters employing the described in vitro experimental system may provide improved platforms to measure a number of biological consequences of shielding materials across the Bragg curve for high charge and energy (HZE) ions.

Desai, N.↗