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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS1985

An experimental and analytical study of the aerodynamic interference effects between two Sears-Haack bodies at Mach 2.7

Aerodynamic interference effects were studied for two slender, streamlined bodies of revolution at Mach 2.7. A wind tunnel investigation produced force and moment data and measurements of pressure distributions on the bodies. As these bodies remained parallel with each other and with the freestream flow, their relative lateral and longitudinal spacing were varied. Results of theoretical methods were used in the analysis of results. The interference effects between the two bodies yielded less total drag than a single body of equal total volume and the same length.

Bantle, J. W.↗

NASA NTRS1986

Variation of the urban NO2 concentration during a ground inversion situation

In February 1985 the NO2 concentration over Goeteborg, a major Swedish city, was monitored during two inversion situations. The differential absorption technique was employed and the results show the build-up of the concentration. In this experiment the laser beam was scanned horizontally over Goeteborg. The results were presented as maps over the NO2 concentration distribution. For security reasons the beam was not parallel to the ground, which gives a height gradient in some measurements. During the measurement period the NO2 concentration variation build-up from the morning to the afternoon was monitored. Great concentration variations were seen. Later the mixing below the inversion layer, resulting in a more uniform, high concentration was seen. The breakdown of the inversion when stronger winds start was also seen. Comparisons with a traditional point-measuring chemiluminiscence instrument and an optical method DOAS (Differential Optical Absorption Spectroscopy), were made. The results are discussed. Strong spatial concentration variations are shown by the lidar measurements. For a point-measuring system, the instrument must be very carefully placed since its location will strongly influence the values it gives.

Sunesson, A.↗

NASA NTRS1990

Two-dimensional shape recognition using sparse distributed memory

Researchers propose a method for recognizing two-dimensional shapes (hand-drawn characters, for example) with an associative memory. The method consists of two stages: first, the image is preprocessed to extract tangents to the contour of the shape; second, the set of tangents is converted to a long bit string for recognition with sparse distributed memory (SDM). SDM provides a simple, massively parallel architecture for an associative memory. Long bit vectors (256 to 1000 bits, for example) serve as both data and addresses to the memory, and patterns are grouped or classified according to similarity in Hamming distance. At the moment, tangents are extracted in a simple manner by progressively blurring the image and then using a Canny-type edge detector (Canny, 1986) to find edges at each stage of blurring. This results in a grid of tangents. While the technique used for obtaining the tangents is at present rather ad hoc, researchers plan to adopt an existing framework for extracting edge orientation information over a variety of resolutions, such as suggested by Watson (1987, 1983), Marr and Hildreth (1980), or Canny (1986).

Kanerva, Pentti↗

NASA NTRS1992

Distribution in magnetotail of O(+) ions from cusp/cleft ionosphere - A possible substorm trigger

The transport of O(+) ions from the cusp/cleft ionosphere to the magnetotail during highly disturbed times was determined by computing the guiding-center trajectories of the ions to a distance of 6 R(E) from the ionosphere and the full-motion trajectories at later times. Case histories were tallied in six planes perpendicular to the X(GSM) axis, three planes perpendicular to the Y(GSM) axis, and in the center plane of the tail. At various times relative to the enhancement of the convection electric field, the following ion properties were constructed from the case histories: number density, mean energy, energy and pitch angle distributions of the flux, and ion pressure components parallel and perpendicular to the magnetic field. It was found that, after about 1.7 hours, the ion flux in the near-earth magnetotail increased dramatically and the spectrum hardened, much as observed during periods just preceding substorms. This increase is attributed to (1) the increase in the O(+) outflux from the ionosphere, (2) the increased energization of the ions by the convection electric field, and (3) ion trapping, which generally occurs because the ion magnetic moments generally increase after the ions first cross the geomagnetotail center plane.

Cladis, J. B.↗

NASA NTRS1992

Electron acceleration by Alfven waves in the magnetosphere

The self-consistent electron kinetics of Alfven waves on the electron inertial scale is studied using a 2D hybrid-kinetic description. The ions follow a fluid description for Alfven waves at frequencies below the ion cyclotron frequency. The parallel electron dynamics is treated kinetically using particle-in-cell techniques. In this model, the electron plasma mode is eliminated, and only the physics of the Alfven waves is retained. At sufficiently large amplitudes, it is found that oblique Alfven waves break due to finite electron inertia in a cold plasma. The consequence of wave breaking is the formation of an electron beam which can be unstable to the beam-plasma instability. The electrons supporting the parallel current thermalize into a non-Maxwellian distribution with an energetic tail up to several keV, assuming a reasonable magnetospheric Alfven speed. In hot plasma simulations, electron trapping is the principal mechanism of electron acceleration. It is proposed that wave breaking or electron trapping of oblique Alfven waves at 1 R(E) can result in electron acceleration and may explain some observed auroral phenomena.

Hui, C.-H.↗

NASA NTRS1993

Technologies for Visualization in Computational Aerosciences

State-of-the-art research in computational aerosciences produces' complex, time-dependent datasets. Simulations can also be multidisciplinary in nature, coupling two or more physical disciplines such as fluid dynamics, structural dynamics, thermodynamics, and acoustics. Many diverse technologies are necessary for visualizing computational aerosciences simulations. This paper describes these technologies and how they contribute to building effective tools for use by domain scientists. These technologies include data management, distributed environments, advanced user interfaces, rapid prototyping environments, parallel computation, and methods to visualize the scalar and vector fields associated with computational aerosciences datasets.

Miceli, Kristina D.↗

NASA NTRS2002

Containerless Ripple Turbulence

One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Putterman, Seth↗

NASA NTRS2003

Computation of Engine Noise Propagation and Scattering Off an Aircraft

The paper presents a comparison of experimental noise data measured in flight on a two-engine business jet aircraft with Kulite microphones placed on the suction surface of the wing with computational results. Both a time-domain discontinuous Galerkin spectral method and a frequency-domain spectral element method are used to simulate the radiation of the dominant spinning mode from the engine and its reflection and scattering by the fuselage and the wing. Both methods are implemented in computer codes that use the distributed memory model to make use of large parallel architectures. The results show that trends of the noise field are well predicted by both methods.

Xu, J.↗

NASA NTRS1995

Slices: A Scalable Partitioner for Finite Element Meshes

A parallel partitioner for partitioning unstructured finite element meshes on distributed memory architectures is developed. The element based partitioner can handle mixtures of different element types. All algorithms adopted in the partitioner are scalable, including a communication template for unpredictable incoming messages, as shown in actual timing measurements.

partitioner finite element meshes parallel computi↗

NASA NTRS1995

Surface Displacement of the May 17, 1993 Eureka Valley, California, Earthquake

Satellite synthetic aperture radar (SAR) interferometry shows that the magnitude 6.1 Eureka Valley earthquake of 17 May 1993 produced an elongated subsidence basin oriented north-northwest, parallel to the trend defined by the aftershock distribution, whereas the source mechanism of the earthquake implies a north-northeast-striking normal fault. The ±3-millimeter accuracy of the radar-observed displacement map over short spatial scales allowed identification of the main surface rupture associated with the event. These observations suggest that the rupture began at depth and propagated diagonally upward and southward on a west-dipping, north-northeast fault plane, reactivating the largest escarpment in the Saline Range.

Gilles Peltzer↗

NASA NTRS2006

Fibonacci Grids

Recent years have seen a resurgence of interest in a variety of non-standard computational grids for global numerical prediction. The motivation has been to reduce problems associated with the converging meridians and the polar singularities of conventional regular latitude-longitude grids. A further impetus has come from the adoption of massively parallel computers, for which it is necessary to distribute work equitably across the processors; this is more practicable for some non-standard grids. Desirable attributes of a grid for high-order spatial finite differencing are: (i) geometrical regularity; (ii) a homogeneous and approximately isotropic spatial resolution; (iii) a low proportion of the grid points where the numerical procedures require special customization (such as near coordinate singularities or grid edges). One family of grid arrangements which, to our knowledge, has never before been applied to numerical weather prediction, but which appears to offer several technical advantages, are what we shall refer to as "Fibonacci grids". They can be thought of as mathematically ideal generalizations of the patterns occurring naturally in the spiral arrangements of seeds and fruit found in sunflower heads and pineapples (to give two of the many botanical examples). These grids possess virtually uniform and highly isotropic resolution, with an equal area for each grid point. There are only two compact singular regions on a sphere that require customized numerics. We demonstrate the practicality of these grids in shallow water simulations, and discuss the prospects for efficiently using these frameworks in three-dimensional semi-implicit and semi-Lagrangian weather prediction or climate models.

Swinbank, Richard↗

NASA NTRS2016

Finite Gyroradius Effects in the Electron Outflow of Asymmetric Magnetic Reconnection

We present observations of asymmetric magnetic reconnection showing evidence of electron demagnetization in the electron outflow. The observations were made at the magnetopause by the four Magnetospheric Multiscale (MMS) spacecraft, separated by approximately 15 km. The reconnecting current sheet has negligible guide field, and all four spacecraft likely pass close to the electron diffusion region just south of the X line. In the electron outflow near the X line, all four spacecraft observe highly structured electron distributions in a region comparable to a few electron gyroradii. The distributions consist of a core with T(sub parallel) greater than T(sub perpendicular) and a nongyrotropic crescent perpendicular to the magnetic field. The crescents are associated with finite gyroradius effects of partly demagnetized electrons. These observations clearly demonstrate the manifestation of finite gyroradius effects in an electron-scale reconnection current sheet.

Norgren, C.↗

DOE OSTI2020

Scalable Deep-Learning-Accelerated Topology Optimization for Additively Manufactured Materials

Topology optimization (TO) is a popular and powerful computational approach for designing novel structures, materials, and devices. Two computational challenges have limited the applicability of TO to a variety of industrial applications. First, a TO problem often involves a large number of design variables to guarantee sufficient expressive power. Second, many TO problems require a large number of expensive physical model simulations, and those simulations cannot be parallelized. To address these issues, we propose a general scalable deep-learning (DL) based TO framework, referred to as SDL-TO, which utilizes parallel schemes in high performance computing (HPC) to accelerate the TO process for designing additively manufactured (AM) materials. Unlike the existing studies of DL for TO, our framework accelerates TO by learning the iterative history data and simultaneously training on the mapping between the given design and its gradient. The surrogate gradient is learned by utilizing parallel computing on multiple CPUs incorporated with a distributed DL training on multiple GPUs. The learned TO gradient enables a fast online update scheme instead of an expensive update based on the physical simulator or solver. Using a local sampling strategy, we achieve to reduce the intrinsic high dimensionality of the design space and improve the training accuracy and the scalability of the SDL-TO framework. The method is demonstrated by benchmark examples and AM materials design for heat conduction. The proposed SDL-TO framework shows competitive performance compared to the baseline methods but significantly reduces the computational cost by a speed up of around 8.6x over the standard TO implementation.

Bi, Sirui↗

NASA NTRS2002

Software Engineering Support of the Third Round of Scientific Grand Challenge Investigations: Earth System Modeling Software Framework Survey

One of the most significant challenges in large-scale climate modeling, as well as in high-performance computing in other scientific fields, is that of effectively integrating many software models from multiple contributors. A software framework facilitates the integration task, both in the development and runtime stages of the simulation. Effective software frameworks reduce the programming burden for the investigators, freeing them to focus more on the science and less on the parallel communication implementation. while maintaining high performance across numerous supercomputer and workstation architectures. This document surveys numerous software frameworks for potential use in Earth science modeling. Several frameworks are evaluated in depth, including Parallel Object-Oriented Methods and Applications (POOMA), Cactus (from (he relativistic physics community), Overture, Goddard Earth Modeling System (GEMS), the National Center for Atmospheric Research Flux Coupler, and UCLA/UCB Distributed Data Broker (DDB). Frameworks evaluated in less detail include ROOT, Parallel Application Workspace (PAWS), and Advanced Large-Scale Integrated Computational Environment (ALICE). A host of other frameworks and related tools are referenced in this context. The frameworks are evaluated individually and also compared with each other.

Talbot, Bryan↗

NASA NTRS1930

Calculation of Pressure Distribution on Airship Hulls

These calculations were based on the shape of the ZR III, with the following simplifications: cars, fins, and rudders removed; all cross sections replaced by equivalent circular cross sections. Under these assumptions the pressure distribution was calculated for the following cases: symmetrical case, or flow parallel to the axis; unsymmetrical case, or flow at an angle to the axis. In both cases the simple potential flow first forms the basis for the determination of the pressure distribution.

Von Karman, Theodor↗

NASA NTRS1994

A Navier-Strokes Chimera Code on the Connection Machine CM-5: Design and Performance

We have implemented a three-dimensional compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the 'chimera' approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. A parallel machine like the CM-5 is well-suited for finite-difference methods on structured grids. The regular pattern of connections of a structured mesh maps well onto the architecture of the machine. So the first design choice, finite differences on a structured mesh, is natural. We use centered differences in space, with added artificial dissipation terms. When numerically solving the Navier-Stokes equations, there are liable to be some mesh cells near a solid body that are small in at least one direction. This mesh cell geometry can impose a very severe CFL (Courant-Friedrichs-Lewy) condition on the time step for explicit time-stepping methods. Thus, though explicit time-stepping is well-suited to the architecture of the machine, we have adopted implicit time-stepping. We have further taken the approximate factorization approach. This creates the need to solve large banded linear systems and creates the first possible barrier to an efficient algorithm. To overcome this first possible barrier we have considered two options. The first is just to solve the banded linear systems with data spread over the whole machine, using whatever fast method is available. This option is adequate for solving scalar tridiagonal systems, but for scalar pentadiagonal or block tridiagonal systems it is somewhat slower than desired. The second option is to 'transpose' the flow and geometry variables as part of the time-stepping process: Start with x-lines of data in-processor. Form explicit terms in x, then transpose so y-lines of data are in-processor. Form explicit terms in y, then transpose so z-lines are in processor. Form explicit terms in z, then solve linear systems in the z-direction. Transpose to the y-direction, then solve linear systems in the y-direction. Finally transpose to the x direction and solve linear systems in the x-direction. This strategy avoids inter-processor communication when differencing and solving linear systems, but requires a large amount of communication when doing the transposes. The transpose method is more efficient than the non-transpose strategy when dealing with scalar pentadiagonal or block tridiagonal systems. For handling geometrically complex problems the chimera strategy was adopted. For multiple zone cases we compute on each zone sequentially (using the whole parallel machine), then send the chimera interpolation data to a distributed data structure (array) laid out over the whole machine. This information transfer implies an irregular communication pattern, and is the second possible barrier to an efficient algorithm. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran. We make use of the Connection Machine Scientific Software Library (CMSSL) for the linear solver and array transpose operations.

Jespersen, Dennis C.↗

NASA NTRS1990

Parallel simulated annealing algorithms for cell placement on hypercube multiprocessors

Two parallel algorithms for standard cell placement using simulated annealing are developed to run on distributed-memory message-passing hypercube multiprocessors. The cells can be mapped in a two-dimensional area of a chip onto processors in an n-dimensional hypercube in two ways, such that both small and large cell exchange and displacement moves can be applied. The computation of the cost function in parallel among all the processors in the hypercube is described, along with a distributed data structure that needs to be stored in the hypercube to support the parallel cost evaluation. A novel tree broadcasting strategy is used extensively for updating cell locations in the parallel environment. A dynamic parallel annealing schedule estimates the errors due to interacting parallel moves and adapts the rate of synchronization automatically. Two novel approaches in controlling error in parallel algorithms are described: heuristic cell coloring and adaptive sequence control.

Banerjee, Prithviraj↗

NASA NTRS2013

THEMIS Observations of the Magnetopause Electron Diffusion Region: Large Amplitude Waves and Heated Electrons

We present the first observations of large amplitude waves in a well-defined electron diffusion region based on the criteria described by Scudder et al at the subsolar magnetopause using data from one Time History of Events and Macroscale Interactions during Substorms (THEMIS) satellite. These waves identified as whistler mode waves, electrostatic solitary waves, lower hybrid waves, and electrostatic electron cyclotron waves, are observed in the same 12 s waveform capture and in association with signatures of active magnetic reconnection. The large amplitude waves in the electron diffusion region are coincident with abrupt increases in electron parallel temperature suggesting strong wave heating. The whistler mode waves, which are at the electron scale and which enable us to probe electron dynamics in the diffusion region were analyzed in detail. The energetic electrons (approx. 30 keV) within the electron diffusion region have anisotropic distributions with T(sub e(right angle))/T(sub e(parallel)) > 1 that may provide the free energy for the whistler mode waves. The energetic anisotropic electrons may be produced during the reconnection process. The whistler mode waves propagate away from the center of the "X-line" along magnetic field lines, suggesting that the electron diffusion region is a possible source region of the whistler mode waves.

magnetopause and boundary layers↗