Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallel applications”

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 685 records · Page 38

Element-topology-independent preconditioners for parallel finite element computations

A family of preconditioners for the solution of finite element equations are presented, which are element-topology independent and thus can be applicable to element order-free parallel computations. A key feature of the present preconditioners is the repeated use of element connectivity matrices and their left and right inverses. The properties and performance of the present preconditioners are demonstrated via beam and two-dimensional finite element matrices for implicit time integration computations.

Park, K. C.↗

LiNbO{sub 3}: A photovoltaic substrate for massive parallel manipulation and patterning of nano-objects

The application of evanescent photovoltaic (PV) fields, generated by visible illumination of Fe:LiNbO{sub 3} substrates, for parallel massive trapping and manipulation of micro- and nano-objects is critically reviewed. The technique has been often referred to as photovoltaic or photorefractive tweezers. The main advantage of the new method is that the involved electrophoretic and/or dielectrophoretic forces do not require any electrodes and large scale manipulation of nano-objects can be easily achieved using the patterning capabilities of light. The paper describes the experimental techniques for particle trapping and the main reported experimental results obtained with a variety of micro- and nano-particles (dielectric and conductive) and different illumination configurations (single beam, holographic geometry, and spatial light modulator projection). The report also pays attention to the physical basis of the method, namely, the coupling of the evanescent photorefractive fields to the dielectric response of the nano-particles. The role of a number of physical parameters such as the contrast and spatial periodicities of the illumination pattern or the particle deposition method is discussed. Moreover, the main properties of the obtained particle patterns in relation to potential applications are summarized, and first demonstrations reviewed. Finally, the PV method is discussed in comparison to other patterning strategies, such as those based on the pyroelectric response and the electric fields associated to domain poling of ferroelectric materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Automated problem scheduling and reduction of synchronization delay effects

It is anticipated that in order to make effective use of many future high performance architectures, programs will have to exhibit at least a medium grained parallelism. A framework is presented for partitioning very sparse triangular systems of linear equations that is designed to produce favorable preformance results in a wide variety of parallel architectures. Efficient methods for solving these systems are of interest because: (1) they provide a useful model problem for use in exploring heuristics for the aggregation, mapping and scheduling of relatively fine grained computations whose data dependencies are specified by directed acrylic graphs, and (2) because such efficient methods can find direct application in the development of parallel algorithms for scientific computation. Simple expressions are derived that describe how to schedule computational work with varying degrees of granularity. The Encore Multimax was used as a hardware simulator to investigate the performance effects of using the partitioning techniques presented in shared memory architectures with varying relative synchronization costs.

Saltz, Joel H.↗

Machine and Application Aware Partitioning for Adaptive Mesh Refinement Applications. In: HPDC '17 Proceedings of the 26th International Symposium on High-Performance Parallel and Distributed Computing

Load balancing and partitioning are critical when it comes to parallel computations. Popular partitioning strategies based on space filling curves focus on equally dividing work. The partitions produced are independent of the architecture or the application. Given the ever-increasing relative cost of data movement and increasing heterogeneity of our architectures, it is no longer sufficient to only consider an equal partitioning of work. Minimizing communication costs are equally if not more important. Our hypothesis is that an unequal partitioning that minimizes communication costs significantly can scale and perform better than conventional equal-work partitioning schemes. This tradeoff is dependent on the architecture as well as the application. We validate our hypothesis in the context of a finite-element computation utilizing adaptive mesh-refinement. Our central contribution is a new partitioning scheme that minimizes the overall runtime of subsequent computations by performing architecture and application-aware non-uniform work assignment in order to decrease time to solution, primarily by minimizing data-movement. We evaluate our algorithm by comparing it against standard space-filling curve based partitioning algorithms and observing time-to-solution as well as energy-to-solution for solving Finite Element computations on adaptively refined meshes. We demonstrate excellent scalability of our new partition algorithm up to $262,144$ cores on ORNL's Titan and demonstrate that the proposed partitioning scheme reduces overall energy as well as time-to-solution for application codes by up to 22.0%.

Fernando, Milinda↗

Numerically stable, scalable formulas for parallel and online computation of higher-order multivariate central moments with arbitrary weights

Formulas for incremental or parallel computation of second order central moments have long been known, and recent extensions of these formulas to univariate and multivariate moments of arbitrary order have been developed. Such formulas are of key importance in scenarios where incremental results are required and in parallel and distributed systems where communication costs are high. We survey these recent results, and improve them with arbitrary-order, numerically stable one-pass formulas which we further extend with weighted and compound variants. We also develop a generalized correction factor for standard two-pass algorithms that enables the maintenance of accuracy over nearly the full representable range of the input, avoiding the need for extended-precision arithmetic. We then empirically examine algorithm correctness for pairwise update formulas up to order four as well as condition number and relative error bounds for eight different central moment formulas, each up to degree six, to address the trade-offs between numerical accuracy and speed of the various algorithms. Finally, we demonstrate the use of the most elaborate among the above mentioned formulas, with the utilization of the compound moments for a practical large-scale scientific application.

97 MATHEMATICS AND COMPUTING↗

An O ( N ) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles

Large classes of materials systems in physics and engineering are governed by magnetic and electrostatic interactions. Continuum or mesoscale descriptions of such systems can be cast in terms of integral equations, whose direct computational evaluation requires O ( N 2 ) operations, where N is the number of unknowns. Such a scaling, which arises from the many-body nature of the relevant Green's function, has precluded wide-spread adoption of integral methods for solution of large-scale scientific and engineering problems. In this work, a parallel computational approach is presented that relies on using scalable open source libraries and utilizes a kernel-independent Fast Multipole Method (FMM) to evaluate the integrals in O ( N ) operations, with O ( N ) memory cost, thereby substantially improving the scalability and efficiency of computational integral methods. We demonstrate the accuracy, efficiency, and scalability of our approach in the context of two examples. In the first, we solve a boundary value problem for a ferroelectric/ferromagnetic volume in free space. In the second, we solve an electrostatic problem involving polarizable dielectric bodies in an unbounded dielectric medium. Lastly, the results from these test cases show that our proposed parallel approach, which is built on a kernel-independent FMM, can enable highly efficient and accurate simulations and allow for considerable flexibility in a broad range of applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Relaxation Matrix for Symmetric Tops with Inversion Symmetry: Line Coupling and Line Mixing Effects on NH3 Lines in the V4 Band

Line shape parameters including the half-widths and the off-diagonal elements of the relaxation matrix have been calculated for self-broadened NH3 lines in the perpendicular v4 band. As in the pure rotational and the parallel v1 bands, the small inversion splitting in this band causes a complete failure of the isolated line approximation. As a result, one has to use formalisms not relying on this approximation. However, due to differences between parallel and perpendicular bands of NH3, the applicability of the formalism used in our previous studies of the v1 band and other parallel bands must be carefully verified. We have found that, as long as potential models only contain components with K1 equals K2 equals 0, whose matrix elements require the selection rule delta k equals 0, the formalism is applicable for the v4 band with some minor adjustments. Based on both theoretical considerations and results from numerical calculations, the non-diagonality of the relaxation matrices in all the PP, RP, PQ, RQ, PR, and RR branches is discussed. Theoretically calculated self-broadened half-widths are compared with measurements and the values listed in HITRAN 2012. With respect to line coupling effects, we have compared our calculated intra-doublet off-diagonal elements of the relaxation matrix with reliable measurements carried out in the PP branch where the spectral environment is favorable. The agreement is rather good since our results do well reproduce the observed k and j dependences of these elements, thus validating our formalism.

Absorption coefficient↗

Understanding the Machine Learning Needs of ECP Applications

In order to support the codesign needs of ECP applications in current and future hardware in the area of machine learning, the ExaLearn team at Sandia studied the different machine learning use cases in three different ECP applications. This report is a summary of the needs of the three applications. The Sandia ExaLearn team will develop a proxy application representative of ECP application needs, specifically the ExaSky and EXAALT ECP projects. The proxy application will allow us to demonstrate performance portable kernels within machine learning codes. Furthermore, current training scalability of machine learning networks in these applications is negatively affected by large batch sizes. Training throughput of the network will increase as batch size increases, but network accuracy and generalization worsens. The proxy application will contain hybrid model- and data-parallelism to improve training efficiency while maintaining network accuracy. The proxy application will also target optimizing 3D convolutional layers, specific to scientific machine learning, which have not been as thoroughly explored by industry.

97 MATHEMATICS AND COMPUTING↗

On the adequacy of lightweight thread approaches for high-level parallel programming models

High-level parallel programming models (PMs) are becoming crucial in order to extract the computational power of current on-node multi-threaded parallelism. The most popular PMs, such as OpenMP or OmpSs, are directive-based: the complexity of the hardware is hidden by the underlying runtime system, improving coding productivity. The implementations of OpenMP usually rely on POSIX threads (pthreads), offering excellent performance for coarse-grained parallelism and a perfect match with the current hardware. OmpSs is a task oriented PM based on an ad hoc runtime solution called Nanos++; it is the precursor of the tasking parallelism in the OpenMP tasking specification. A recent trend in runtimes and applications points to leveraging massive on-node parallelism in conjunction with fine-grained and dynamic scheduling paradigms. In this paper we analyze the behavior of the OpenMP and OmpSs PMs on top of the recently emerged Generic Lightweight Threads (GLT) API. GLT exposses a common API for lightweight thread (LWT) libraries that offers the possibility of running the same application over different native LWT solutions. We describe the design details of those high-level PMs implemented on top of GLT and analyze different scenarios in order to assess where the use of LWTs may benefit application performance. Our work reveals those scenarios where LWTs overperform pthread-based solutions and compares the performance between an ad hoc solution and a generic implementation.

97 MATHEMATICS AND COMPUTING↗

Development of a Navier-Stokes algorithm for parallel-processing supercomputers

An explicit flow solver, applicable to the hierarchy of model equations ranging from Euler to full Navier-Stokes, is combined with several techniques designed to reduce computational expense. The computational domain consists of local grid refinements embedded in a global coarse mesh, where the locations of these refinements are defined by the physics of the flow. Flow characteristics are also used to determine which set of model equations is appropriate for solution in each region, thereby reducing not only the number of grid points at which the solution must be obtained, but also the computational effort required to get that solution. Acceleration to steady-state is achieved by applying multigrid on each of the subgrids, regardless of the particular model equations being solved. Since each of these components is explicit, advantage can readily be taken of the vector- and parallel-processing capabilities of machines such as the Cray X-MP and Cray-2.

Swisshelm, Julie M.↗

Development of an aerodynamics algorithm for parallel-processing supercomputers

An explicit flow solver, applicable to the hierarchy of model equations ranging from Euler to full Navier-Stokes, is combined with several techniques designed to reduce computational expense. The computational domain consists of local grid refinements embedded in a global coarse mesh, where the locations of these refinements are defined by the physics of the flow. Flow characteristics are also used to determine which set of model equations is appropriate for solution in each region, thereby reducing not only the number of grid points at which the solution must be obtained, but also the computational effort required to get that solution. Acceleration to steady-state is achieved by applying multigrid on each of the subgrids, regardless of the particular model equations being solved. Since each of these components is explicit, advantage can readily be taken of the vector- and parallel-processing capabilities of machines such as the Cray X-MP and Cray-2.

Swisshelm, Julie M.↗

A Tradeoff Analysis of Series / Parallel Three-Phase Converter Topologies for Wireless Extreme Chargers

In this paper, extreme fast charging (XFC) technology is studied considering the charge rates of 300 kW for wireless power transfer (WPT) applications. Tradeoff analysis of series and parallel connection of three-phase WPT system are presented by comparisons of voltage and current stresses on power electronics active / passive components. In addition, star (Y) / delta (Δ) connection configurations for three-phase wireless power transfer coupling coils are analyzed with series and LCC resonant compensation circuits. The system series and parallel connection controllability is also reviewed considering voltage and current balance techniques with output control. In a conclusion of evaluation analysis, it is revealed that each component of 300 kW wireless charging network must be designed for high fast charging system and the overall system operation need to be strategically planned for high power charging and infrastructure deployment.

Asa, Erdem↗

Putting the 'super' in supercomputers

Computers used for numerical simulations of physical phenomena, e.g., flowfields, meteorology, structural analysis, etc., replace physical experiments that are too expensive or impossible to perform. The problems considered continually become increasingly more complex and thus demand faster processing times to do all necessary computations. The effects components technologies have on computer speed are leveling off, leaving new architectures and programming as the only currently viable means to upgrade speed. Parallel computations, either in the form of array processors, assembly line processing or multiprocessors are being explored using existing microprocessor technologies. Slower hardware configurations can also be made equivalent to faster supercomputers by economic programming. The availability of rudimentary parallel architecture supercomputers for general industrial use is increasing. Scientific applications continue to drive the development of more sophisticated parallel machines.

Schulbach, C.↗

Recent developments and applications of the SEC vidicon for astronomy

The engineering development of the SEC vidicon as an astronomical sensor has continued in parallel with its operational use. One scientific application was the six hour exposure of the quasar PHL-957 at the Coude spectrograph of the 200-inch Hale telescope. The developmental effort includes both the enhancement of the basic attributes that make the SEC an appropriate sensor, namely, high quantum efficiency, low threshold, and long integration; as well as work to broaden its scientific usefulness, such as the development of a MgF2 photocathode window for vacuum ultraviolet sensitivity, and a permanent magnet focus design for thermal compatibility with proposed large space telescopes. Additional details on the characteristics of the SEC tube are discussed, as well as plans to make a larger and higher resolution version.

Zucchino, P.↗

Gratings and Random Reflectors for Near-Infrared PIN Diodes

Crossed diffraction gratings and random reflectors have been proposed as means to increase the quantum efficiencies of InGaAs/InP positive/intrinsic/ negative (PIN) diodes designed to operate as near-infrared photodetectors. The proposal is meant especially to apply to focal-plane imaging arrays of such photodetectors to be used for near-infrared imaging. A further increase in quantum efficiency near the short-wavelength limit of the near-infrared spectrum of such a photodetector array could be effected by removing the InP substrate of the array. The use of crossed diffraction gratings and random reflectors as optical devices for increasing the quantum efficiencies of quantum-well infrared photodetectors (QWIPs) was discussed in several prior NASA Tech Briefs articles. While the optical effects of crossed gratings and random reflectors as applied to PIN photodiodes would be similar to those of crossed gratings and random reflectors as applied to QWIPs, the physical mechanisms by which these optical effects would enhance efficiency differ between the PIN-photodiode and QWIP cases: In a QWIP, the multiple-quantum-well layers are typically oriented parallel to the focal plane and therefore perpendicular or nearly perpendicular to the direction of incidence of infrared light. By virtue of the applicable quantum selection rules, light polarized parallel to the focal plane (as normally incident light is) cannot excite charge carriers and, hence, cannot be detected. A pair of crossed gratings or a random reflector scatters normally or nearly normally incident light so that a significant portion of it attains a component of polarization normal to the focal plane and, hence, can excite charge carriers. A pair of crossed gratings or a random reflector on a PIN photodiode would also scatter light into directions away from the perpendicular to the focal plane. However, in this case, the reason for redirecting light away from the perpendicular is to increase the length of the optical path through the detector to increase the probability of absorption of photons and thereby increase the resulting excitation of charge carriers. A pair of crossed gratings or a random reflector according to the proposal would be fabricated as an integral part of photodetector structure on the face opposite the focal plane (see figure). In the presence of crossed gratings, light would make four passes through the device before departing. In the presence of a random reflector, a significant portion of the light would make more than four passes: After each bounce, light would be scattered at a different random angle, and would have a chance to escape only when it was reflected, relative to the normal, at an angle less than the critical angle for total internal reflection. Given the indices of refraction of the photodiode materials, this angle would be about 17 . This amounts to a very narrow cone for escape of trapped light.

Gunapala, Sarath↗

Metal Side Reflectors for Trapping Light in QWIPs

Focal-plane arrays of quantum-well infrared photodetectors (QWIPs) equipped with both light-coupling diffraction gratings and metal side reflectors have been proposed, and prototypes are expected to be fabricated soon. The purpose served by the metal side reflectors is to increase quantum efficiency by helping to trap light in the photosensitive material of each pixel. The reasons for using diffraction gratings were discussed in several prior NASA Tech Briefs articles. To recapitulate: In an array of QWIPs, the quantum-well layers are typically oriented parallel to the focal plane and therefore perpendicular or nearly perpendicular to the direction of incidence of infrared light. By virtue of the applicable quantum selection rules, light polarized parallel to the focal plane (as normally incident light is) cannot excite charge carriers and, hence, cannot be detected. Diffraction gratings scatter normally or nearly normally incident light into directions more nearly parallel to the focal plane, so that a significant portion of the light attains a component of polarization normal to the focal plane and, hence, can excite charge carriers. Unfortunately, light scattered in directions parallel or nearly parallel to the focal plane can escape sideways from the QWIP of a given pixel, as illustrated in Figure 1. The escaped light has made only a single pass through the interior photosensitive volume of the QWIP. The quantum efficiency of the QWIP would be increased by trapping light so that it makes multiple passes through the photosensitive volume. As shown in Figure 2, the sides of the QWIP of each pixel would be coated with gold to reflect escaping light back into the interior.

Gunapala, Sarath↗

Ion streaming instabilities with application to collisionless shock wave structure

The electromagnetic dispersion relation for two counterstreaming ion beams of arbitrary relative strength flowing parallel to a dc magnetic field is derived. The beams flow through a stationary electron background and the dispersion relation in the fluid approximation is unaffected by the electron thermal pressure. The dispersion relation is solved with a zero net current condition applied and the regions of instability in the k-U space (U is the relative velocity between the two ion beams) are presented. The parameters are then chosen to be applicable for parallel shocks. It was found that unstable waves with zero group velocity in the shock frame can exist near the leading edge of the shock for upstream Alfven Mach numbers greater than 5.5. It is suggested that this mechanism could generate sufficient turbulence within the shock layer to scatter the incoming ions and create the required dissipation for intermediate strength shocks.

Golden, K. I.↗

Axom

The Axom library consists of n collection of software components that provide core computer science infrastructure capabilities that can be shared by diverse high performance computing (HPC) applications. The current set of capabilities that Axom provides includes: Customizable support for error/warning and diagnostic message reporting and coordination among components of integrated applications (e.g., physics packages, libraries, etc.); Scalable parallel aggregation and filtering of message logging; In memory datastore for hierarchical, mesh-aware simulation data. The datastore supports data description, allocation, deallocation. parallel l/0, etc; Mesh data model abstractions that enable the development of mesh-based computational algorithms that work with many diverse, application· specific mesh types; and Computational geometry primitives for developing high performance parallel algorithms based on spatial geometry and meshes, e.g., spatial acceleration data structures, sets, maps, relations, etc.

Hornung, RichardD.↗