Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel processing”

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,405 records · Page 78

I/O-Efficient Scientific Computation Using TPIE

In recent years, input/output (I/O)-efficient algorithms for a wide variety of problems have appeared in the literature. However, systems specifically designed to assist programmers in implementing such algorithms have remained scarce. TPIE is a system designed to support I/O-efficient paradigms for problems from a variety of domains, including computational geometry, graph algorithms, and scientific computation. The TPIE interface frees programmers from having to deal not only with explicit read and write calls, but also the complex memory management that must be performed for I/O-efficient computation. In this paper we discuss applications of TPIE to problems in scientific computation. We discuss algorithmic issues underlying the design and implementation of the relevant components of TPIE and present performance results of programs written to solve a series of benchmark problems using our current TPIE prototype. Some of the benchmarks we present are based on the NAS parallel benchmarks while others are of our own creation. We demonstrate that the central processing unit (CPU) overhead required to manage I/O is small and that even with just a single disk, the I/O overhead of I/O-efficient computation ranges from negligible to the same order of magnitude as CPU time. We conjecture that if we use a number of disks in parallel this overhead can be all but eliminated.

Vengroff, Darren Erik↗

Optimizing the hit finding algorithm for liquid argon TPC neutrino detectors using parallel architectures

Neutrinos are particles that interact rarely, so identifying them requires large detectors which produce lots of data. Processing this data with the computing power available is becoming even more difficult as the detectors increase in size to reach their physics goals. Liquid argon time projection chamber (LArTPC) neutrino experiments are expected to grow in the next decade to have 100 times more wires than in currently operating experiments, and modernization of LArTPC reconstruction code, including parallelization both at data- and instruction-level, will help to mitigate this challenge. The LArTPC hit finding algorithm is used across multiple experiments through a common software framework. In this paper we discuss a parallel implementation of this algorithm. Using a standalone setup we find speedup factors of two times from vectorization and 30–100 times from multi-threading on Intel architectures. The new version has been incorporated back into the framework so that it can be used by experiments. On a serial execution, the integrated version is about 10 times faster than the previous one and, once parallelization is enabled, more speedups comparable to the standalone program are achieved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Multi-Level Parallelization Concept for High-Fidelity Multi-Block Solvers

The integration of high-fidelity Computational Fluid Dynamics (CFD) analysis tools with the industrial design process benefits greatly from the robust implementations that are transportable across a wide range of computer architectures. In the present work, a hybrid domain-decomposition and parallelization concept was developed and implemented into the widely-used NASA multi-block Computational Fluid Dynamics (CFD) packages implemented in ENSAERO and OVERFLOW. The new parallel solver concept, PENS (Parallel Euler Navier-Stokes Solver), employs both fine and coarse granularity in data partitioning as well as data coalescing to obtain the desired load-balance characteristics on the available computer platforms. This multi-level parallelism implementation itself introduces no changes to the numerical results, hence the original fidelity of the packages are identically preserved. The present implementation uses the Message Passing Interface (MPI) library for interprocessor message passing and memory accessing. By choosing an appropriate combination of the available partitioning and coalescing capabilities only during the execution stage, the PENS solver becomes adaptable to different computer architectures from shared-memory to distributed-memory platforms with varying degrees of parallelism. The PENS implementation on the IBM SP2 distributed memory environment at the NASA Ames Research Center obtains 85 percent scalable parallel performance using fine-grain partitioning of single-block CFD domains using up to 128 wide computational nodes. Multi-block CFD simulations of complete aircraft simulations achieve 75 percent perfect load-balanced executions using data coalescing and the two levels of parallelism. SGI PowerChallenge, SGI Origin 2000, and a cluster of workstations are the other platforms where the robustness of the implementation is tested. The performance behavior on the other computer platforms with a variety of realistic problems will be included as this on-going study progresses.

Hatay, Ferhat F.↗

Density and Magnetic Field Asymmetric Kelvin‐Helmholtz Instability

Abstract The Kelvin‐Helmholtz (KH) instability can transport mass, momentum, magnetic flux, and energy between the magnetosheath and magnetosphere, which plays an important role in the solar‐wind‐magnetosphere coupling process for different planets. Meanwhile, strong density and magnetic field asymmetry are often present between the magnetosheath (MSH) and magnetosphere (MSP), which could affect the transport processes driven by the KH instability. Our magnetohydrodynamics simulation shows that the KH growth rate is insensitive to the density ratio between the MSP and the MSH in the compressible regime, which is different than the prediction from linear incompressible theory. When the interplanetary magnetic field (IMF) is parallel to the planet's magnetic field, the nonlinear KH instability can drive a double mid‐latitude reconnection (DMLR) process. The total double reconnected flux depends on the KH wavelength and the strength of the lower magnetic field. When the IMF is anti‐parallel to the planet's magnetic field, the nonlinear interaction between magnetic reconnection and the KH instability leads to fast reconnection (i.e., close to Petschek reconnection even without including kinetic physics). However, the peak value of the reconnection rate still follows the asymmetric reconnection scaling laws. We also demonstrate that the DMLR process driven by the KH instability mixes the plasma from different regions and consequently generates different types of velocity distribution functions. We show that the counter‐streaming beams can be simply generated via the change of the flux tube connection and do not require parallel electric fields.

Astronomy & Astrophysics↗

Extended-Range Ultrarefractive 1D Photonic Crystal Prisms

A proposal has been made to exploit the special wavelength-dispersive characteristics of devices of the type described in One-Dimensional Photonic Crystal Superprisms (NPO-30232) NASA Tech Briefs, Vol. 29, No. 4 (April 2005), page 10a. A photonic crystal is an optical component that has a periodic structure comprising two dielectric materials with high dielectric contrast (e.g., a semiconductor and air), with geometrical feature sizes comparable to or smaller than light wavelengths of interest. Experimental superprisms have been realized as photonic crystals having three-dimensional (3D) structures comprising regions of amorphous Si alternating with regions of SiO2, fabricated in a complex process that included sputtering. A photonic crystal of the type to be exploited according to the present proposal is said to be one-dimensional (1D) because its contrasting dielectric materials would be stacked in parallel planar layers; in other words, there would be spatial periodicity in one dimension only. The processes of designing and fabricating 1D photonic crystal superprisms would be simpler and, hence, would cost less than do those for 3D photonic crystal superprisms. As in 3D structures, 1D photonic crystals may be used in applications such as wavelength-division multiplexing. In the extended-range configuration, it is also suitable for spectrometry applications. As an engineered structure or artificially engineered material, a photonic crystal can exhibit optical properties not commonly found in natural substances. Prior research had revealed several classes of photonic crystal structures for which the propagation of electromagnetic radiation is forbidden in certain frequency ranges, denoted photonic bandgaps. It had also been found that in narrow frequency bands just outside the photonic bandgaps, the angular wavelength dispersion of electromagnetic waves propagating in photonic crystal superprisms is much stronger than is the angular wavelength dispersion obtained by use of conventional prisms and diffraction gratings and is highly nonlinear.

Ting, David Z.↗

Improved silicon carbide for advanced heat engines

The development of silicon carbide materials of high strength was initiated and components of complex shape and high reliability were formed. The approach was to adapt a beta-SiC powder and binder system to the injection molding process and to develop procedures and process parameters capable of providing a sintered silicon carbide material with improved properties. The initial effort was to characterize the baseline precursor materials, develop mixing and injection molding procedures for fabricating test bars, and characterize the properties of the sintered materials. Parallel studies of various mixing, dewaxing, and sintering procedures were performed in order to distinguish process routes for improving material properties. A total of 276 modulus-of-rupture (MOR) bars of the baseline material was molded, and 122 bars were fully processed to a sinter density of approximately 95 percent. Fluid mixing techniques were developed which significantly reduced flaw size and improved the strength of the material. Initial MOR tests indicated that strength of the fluid-mixed material exceeds the baseline property by more than 33 percent. the baseline property by more than 33 percent.

Whalen, Thomas J.↗

Science of the Van Allen Probes Science Operations Centers

The Van Allen Probes mission operations materialized through a distributed model in which operational responsibility was divided between the Mission Operations Center (MOC) and separate instrument specific SOCs. The sole MOC handled all aspects of telemetering and receiving tasks as well as certain scientifically relevant ancillary tasks. Each instrument science team developed individual instrument specific SOCs proficient in unique capabilities in support of science data acquisition, data processing, instrument performance, and tools for the instrument team scientists. In parallel activities, project scientists took on the task of providing a significant modeling tool base usable by the instrument science teams and the larger scientific community. With a mission as complex as Van Allen Probes, scientific inquiry occurred due to constant and significant collaboration between the SOCs and in concert with the project science team. Planned cross-instrument coordinated observations resulted in critical discoveries during the seven-year mission. Instrument cross-calibration activities elucidated a more seamless set of data products. Specific topics include post-launch changes and enhancements to the SOCs, discussion of coordination activities between the SOCs, SOC specific analysis software, modeling software provided by the Van Allen Probes project, and a section on lessons learned. One of the most significant lessons learned was the importance of the original decision to implement individual team SOCs providing timely and well-documented instrument data for the NASA Van Allen Probes Mission scientists and the larger magnetospheric and radiation belt scientific community.

47 OTHER INSTRUMENTATION↗

Revealing Deformation Mechanisms in Polymer-Grafted Thermoplastic Elastomers via In Situ Small-Angle X-ray Scattering

The tunable properties of thermoplastic elastomers (TPEs), through polymer chemistry manipulations, enable these technologically critical materials to be employed in a broad range of applications. The need to “dial-in” the mechanical properties and responses of TPEs generally requires the design and synthesis of new macromolecules. In these designs, TPEs with nonlinear macromolecular architectures outperform the mechanical properties of their linear copolymer counterparts, but the differences in the deformation mechanism providing enhanced performance are unknown. Here, in situ small-angle X-ray scattering (SAXS) measurements during uniaxial extension reveal distinct deformation mechanisms between a commercially available linear poly(styrene)–poly(butadiene)–poly(styrene) (SBS) triblock copolymer and the grafted SBS version containing grafted poly(styrene) (PS) chains from the poly(butadiene) (PBD) midblock. The neat SBS (φ SBS = 100%) sample deforms congruently with the macroscopic dimensions, with the domain spacing between spheres increasing and decreasing along and transverse to the stretch direction, respectively. At high extensions, end segment pullout from the PS-rich domains is detected, which is indicated by a disordering of SBS. Conversely, the PS-grafted SBS that is 30 vol % SBS and 70% styrene (φ SBS = 30%) exhibits a lamellar morphology, and in situ SAXS measurements reveal an unexpected deformation mechanism. During deformation, there are two simultaneous processes: significant lamellar domain rearrangement to preferentially orient the lamellae planes parallel to the stretch direction and crazing. The samples whiten at high strains as expected for crazing, which corresponds with the emergence of features in the 2D SAXS pattern during stretching consistent with fibril-like structures that bridge the voids in crazes. The significant domain rearrangement in the grafted copolymers is attributed to the new junctions formed across multiple PS domains by the grafting of a single chain. In conclusion, the in situ SAXS measurements provide insights into the enhanced mechanical properties of grafted copolymers that arise through improved physical cross-linking that leads to nanostructure domain reorientation for self-reinforcement and craze formation where fibrils help to strengthen the polymer.

36 MATERIALS SCIENCE↗

Generalized large optics fabrication multiplexing

High precision astronomical optics are manufactured through deterministic computer controlled optical surfacing processes, such as subaperture small tool polishing, magnetorheological finishing, bonnet tool polishing, and ion beam figuring. Due to the small tool size and the corresponding tool influence function, large optics fabrication is a highly time-consuming process. The framework of multiplexed figuring runs for the simultaneous use of two or more tools is presented. This multiplexing process increases the manufacturing efficiency and reduces the overall cost using parallelized subaperture tools.

36 MATERIALS SCIENCE↗

Multiphysics analysis system for heat pipe cooled micro-reactors employing PRAGMA-OpenFOAM-ANLHTP

A multiphysics analysis system for neutronics/thermo-mechanical/heat-pipe analysis of heat pipe cooled micro-reactors was developed using the PRAGMA code as the neutronics engine. PRAGMA, which used to be a GPU-based continuous-energy MC code for power reactor applications, now has an extended geometry package to handle geometries with unstructured meshes generated by the ANSYS Design-Modeler and Meshing. The NVIDIA ray tracing engine OptiX was exploited for efficient neutron transport on unstructured geometry. On the multiphysics side, the open-source CFD tool OpenFOAM and one-dimensional heat pipe analysis code ANLHTP were adopted. The manager-worker system based on the MPI dynamic process management (DPM) model enables efficient coupling of codes employing different parallelization schemes. With all the features, the multiphysics analysis of the one-sixth symmetrical MegaPower 2D core was performed and it demonstrated the benefit of the tight integration of the three-way coupling system and one-to-one geometry coupling strategy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Nonresonant particle acceleration in strong turbulence: Comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

ExTreeM: Scalable Augmented Merge Tree Computation via Extremum Graphs

Over the last decade merge trees have been proven to support a plethora of visualization and analysis tasks since they effectively abstract complex datasets. Here, this paper describes the ExTreeM-Algorithm: A scalable algorithm for the computation of merge trees via extremum graphs. The core idea of ExTreeM is to first derive the extremum graph G of an input scalar field f defined on a cell complex K, and subsequently compute the unaugmented merge tree of f on G instead of K; which are equivalent. Any merge tree algorithm can be carried out significantly faster on G, since K in general contains substantially more cells than G. To further speed up computation, ExTreeM includes a tailored procedure to derive merge trees of extremum graphs. The computation of the fully augmented merge tree, i.e., a merge tree domain segmentation of K, can then be performed in an optional post-processing step. All steps of ExTreeM consist of procedures with high parallel efficiency, and we provide a formal proof of its correctness. Our experiments, performed on publicly available datasets, report a speedup of up to one order of magnitude over the state-of-the-art algorithms included in the TTK and VTK-m software libraries, while also requiring significantly less memory and exhibiting excellent scaling behavior.

97 MATHEMATICS AND COMPUTING↗

How to form a wormhole

Abstract We provide a simple but very useful description of the process of wormhole formation. We place two massive objects in two parallel universes (modeled by two branes). Gravitational attraction between the objects competes with the resistance coming from the brane tension. For sufficiently strong attraction, the branes are deformed, objects touch and a wormhole is formed. Our calculations show that more massive and compact objects are more likely to fulfill the conditions for wormhole formation. This implies that we should be looking for wormholes either in the background of black holes and compact stars, or massive microscopic relics. Our formation mechanism applies equally well for a wormhole connecting two objects in the same universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Magnetopause and boundary layer

A brief overview is given of our present knowledge, observational and theoretical, of the structure of the magnetopause and the adjoining plasma boundary layer. Particular attention is given to the relationship between these electromagnetic and plasma structures on the front lobe of the magnetosphere and the magnetic field reconnection process. Items discussed include: magnetopause thickness; behavior of magnetic field components parallel and perpendicular to the magnetopause; particle energization; structure of the boundary layer from reconnection theory.

Sonnerup, B. U. O.↗

Beyond the supercomputer

A NASA-directed development of massively parallel processor (MPP) computers is outlined, noting intended applications for data processing for near term earth resource and environment mapping, radar, and television transmissions. The MPP is designed to perform 100 billion operations/sec to obtain satisfactory image processing, while separate processing units correct distortions, register images, calculate correlation functions, and classify multispectral characteristics. Arrays of 1s and 0s will be manipulated in analog-to-digital conversions generating separate planes corresponding to powers of binaries. Data wires are replaced by fiber-optic tubes or thousands of wires, and single logic gates are replaced by thousands of logic gates and every memory element by thousands of memory elements. Features of the interconnections and the images control processor units are detailed, along with implementation of sliders for program flexibility.

Schaefer, D. H.↗

Real-time data compressor for Eos-class missions

A conceptual design for a real-time VLSI compressor capable of processing rate up to one gigabit per second is presented. This scheme is capable of providing a three-to-one distortion-free data reduction factor to both the High Resolution Imaging Spectrometer and processed SAR imaging data. The design uses a VLSI parallel/piplined architecture capable of processing at a real time rate. The design consists of a parallel array of VLSI compressor modules. Each module is built on a single customized VLSI chip using existing state-of-the-art semiconductor technology.

Lee, Jun-Ji↗

Excitation of whistlers and waves with mixed polarization by newborn cometary ions

The present study has been motivated by the ICE wave measurements. It is found that the newborn cometary ions, particularly the protons, can excite whistlers and waves with frequencies much higher than the proton gyrofrequency but with mixed electrostatic and electromagnetic polarization. For the case of oblique propagation the newborn ions are treated as if they are unmagnetized. This is justified not only because the wave frequencies are high but also because the growth rates are large. On the other hand in the case of parallel propagation the growth rate is much smaller, and the excitation process seems to be unimportant.

Wu, C. S.↗

A tesselated probabilistic representation for spatial robot perception and navigation

The ability to recover robust spatial descriptions from sensory information and to efficiently utilize these descriptions in appropriate planning and problem-solving activities are crucial requirements for the development of more powerful robotic systems. Traditional approaches to sensor interpretation, with their emphasis on geometric models, are of limited use for autonomous mobile robots operating in and exploring unknown and unstructured environments. Here, researchers present a new approach to robot perception that addresses such scenarios using a probabilistic tesselated representation of spatial information called the Occupancy Grid. The Occupancy Grid is a multi-dimensional random field that maintains stochastic estimates of the occupancy state of each cell in the grid. The cell estimates are obtained by interpreting incoming range readings using probabilistic models that capture the uncertainty in the spatial information provided by the sensor. A Bayesian estimation procedure allows the incremental updating of the map using readings taken from several sensors over multiple points of view. An overview of the Occupancy Grid framework is given, and its application to a number of problems in mobile robot mapping and navigation are illustrated. It is argued that a number of robotic problem-solving activities can be performed directly on the Occupancy Grid representation. Some parallels are drawn between operations on Occupancy Grids and related image processing operations.

Elfes, Alberto↗