Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel data structures”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

GridOPTICS/GridPACK

GridPACK is a software framework consisting of a set of modules designed to simplify the development of programs that model the power grid and run on parallel, high performance computing platforms. It also contains several fully developed applications, including powerflow, dynamic simulation, state estimation, Kalman filter analysis (dynamic state estimation), contingency analysis and real time path rating. These applications can be used either standalone or as components in more complicated workflows that combine several different types of application together. The framework modules are available as a combination of libraries and software templates and consist of components for setting up and distributing power grid networks, support for modeling the behavior of individual buses and branches in the network, converting the network models to the corresponding algebraic equations, and parallel routines for manipulating and solving large algebraic systems. The framework also contains a module for distributing tasks evenly amongst computing resources, even if individual tasks vary widely in their execution times. Additional modules support input and output, basic statistical analysis of contingency based calculations, distributed data structures, as well as basic profiling and error management.

Palmer, Bruce↗
NASA NTRS1986

Systolic VLSI for Kalman filters

A novel two-dimensional parallel computing method for real-time Kalman filtering is presented. The mathematical formulation of a Kalman filter algorithm is rearranged to be the type of Faddeev algorithm for generalizing signal processing. The data flow mapping from the Faddeev algorithm to a two-dimensional concurrent computing structure is developed. The architecture of the resulting processor cells is regular, simple, expandable, and therefore naturally suitable for VLSI chip implementation. The computing methodology and the two-dimensional systolic arrays are useful for Kalman filter applications as well as other matrix/vector based algebraic computations.

Yeh, H.-G.↗
NASA NTRS1989

Parallel processing environment for multi-flexible body dynamics

The implementation of a dynamics solution algorithm with inherent parallelism which is applicable to the dynamics of large flexible space structures is described. The algorithm is unique in that parts of the solution can be computed simultaneously by working with different branches of its tree topology. The algorithm exhibits close to 0(n) type behavior. The data flow within the solution algorithm is discussed along with results from its implementation in a multiprocessing environment. A model of the United States Space Station is used as an example. The results show that, with fast multiple scalar processors, an efficient algorithm, and symbolically generated equations of motion, real-time performance can be achieved with present-day hardware technology, even with complex dynamical models.

Venugopal, Ravi↗
NASA NTRS2006

Computational Analysis of a Chevron Nozzle Uniquely Tailored for Propulsion Airframe Aeroacoustics

A computational flow field and predicted jet noise source analysis is presented for asymmetrical fan chevrons on a modern separate flow nozzle at take off conditions. The propulsion airframe aeroacoustic asymmetric fan nozzle is designed with an azimuthally varying chevron pattern with longer chevrons close to the pylon. A baseline round nozzle without chevrons and a reference nozzle with azimuthally uniform chevrons are also studied. The intent of the asymmetric fan chevron nozzle was to improve the noise reduction potential by creating a favorable propulsion airframe aeroacoustic interaction effect between the pylon and chevron nozzle. This favorable interaction and improved noise reduction was observed in model scale tests and flight test data and has been reported in other studies. The goal of this study was to identify the fundamental flow and noise source mechanisms. The flow simulation uses the asymptotically steady, compressible Reynolds averaged Navier-Stokes equations on a structured grid. Flow computations are performed using the parallel, multi-block, structured grid code PAB3D. Local noise sources were mapped and integrated computationally using the Jet3D code based upon the Lighthill Acoustic Analogy with anisotropic Reynolds stress modeling. In this study, trends of noise reduction were correctly predicted. Jet3D was also utilized to produce noise source maps that were then correlated to local flow features. The flow studies show that asymmetry of the longer fan chevrons near the pylon work to reduce the strength of the secondary flow induced by the pylon itself, such that the asymmetric merging of the fan and core shear layers is significantly delayed. The effect is to reduce the peak turbulence kinetic energy and shift it downstream, reducing overall noise production. This combined flow and noise prediction approach has yielded considerable understanding of the physics of a fan chevron nozzle designed to include propulsion airframe aeroacoustic interaction effects.

Massey, Steven J.↗
DOE OSTI2017

Distributed-Memory Breadth-First Search on Massive Graphs. IN: Parallel Graph Algorithms

This chapter studies the problem of traversing large graphs using the breadth-first search order on distributed-memory supercomputers. We consider both the traditional level-synchronous top-down algorithm as well as the recently discovered direction optimizing algorithm. We analyze the performance and scalability trade-offs in using different local data structures such as CSR and DCSC, enabling in-node multithreading, and graph decompositions such as 1D and 2D decomposition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Crystal structure of vismodegib, C 19 H 14 Cl 2 N 2 O 3 S

The crystal structure of vismodegib has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Vismodegib crystallizes in space group P2 1 /a (#14) with a = 16.92070(20), b = 10.20235(4), c = 12.16161(10) Å, β = 108.6802(3)°, V = 1988.873(9) Å 3 , and Z = 4. The crystal structure consists of corrugated layers of molecules parallel to the bc-plane. There is only one classical hydrogen bond in the structure, between the amide nitrogen atom and the N atom of the pyridine ring. Pairs of these hydrogen bonds link the molecules into dimers, with a graph set R2,2(14) > a > a. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal structure of oxibendazole, C 12 H 15 N 3 O 3

The crystal structure of oxibendazole has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Oxibendazole crystallizes in space group C2/c (#15) with a = 23.18673(22), b = 5.35136(5), c = 19.88932(13) Å, β = 97.0876(9)°, V = 2449.018(17) Å 3 , and Z = 8. The structure consists of hydrogen-bonded layers of planar molecules parallel to the bc-plane. Strong N–H∙∙∙N hydrogen bonds link the molecules into dimers, with a graph set R2,2(8). N–H∙∙∙O hydrogen bonds further link these dimers into layers parallel to the bc-plane. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal structure of elvitegravir Form II, C 23 H 23 ClFNO 5

The crystal structure of elvitegravir Form II has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Elvitegravir Form II crystallizes in space group P2 1 (#4) with a = 11.54842(7), b = 14.04367(5), c = 13.33333(8) Å, β = 90.0330(6)°, V = 2162.427(14) Å 3 , and Z = 4. The crystal structure consists of alternating layers of parallel molecules perpendicular to the b-axis. The mean planes of the oxoquinoline ring systems in molecules 1 and 2 are 1(22)-1 and -1(22)1. Between the stacks are layers of the halogenated phenyl rings. These exhibit herringbone stacking. In each molecule, the carboxylic acid group forms a strong intramolecular O–H∙∙∙O hydrogen bond to the nearby carbonyl group. The hydroxyl group of each molecule forms a strong hydrogen bond to the carbonyl group of the carboxylic acid of the other molecule. These O–H∙∙∙O hydrogen bonds link the molecules into dimers, with a graph set R2,2(18) > a > c. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
NASA NTRS1976

Phenomenology of the earth's bow shock system - A summary description of experimental results

Observational data on the earth's bow shock system are classified and characterized. Foreshock components, midshock components, and aftershock components are discussed separately. Schematic representations of the field and plasma particle parameters are elaborated, with attention given to quasi-perpendicular geometry and quasi-parallel geometry. Magnetic pulsation structure is delineated. Schematic profiles of field, particle, and wave behavior through a representative quasi-perpendicular shock crossing are displayed.

Greenstadt, E. W.↗
NASA NTRS1995

Exploratory flow visualization investigation of mast-mounted sights in presence of a rotor

A flow visualization investigation with a laser light sheet system was conducted on a 27-percent-scale AH-64 attack helicopter model fitted with two mast-mounted sights in the langley 14- by 22-foot subsonic tunnel. The investigation was conducted to identify aerodynamic phenomena that may have contributed to adverse vibration encountered during full-scale flight of the AH-64D apache/longbow helicopter with an asymmetric mast-mounted sight. Symmetric and asymmetric mast-mounted sights oriented at several skew angles were tested at simulated forward and rearward flight speeds of 30 and 45 knots. A laser light sheet system was used to visualize the flow in planes parallel to and perpendicular to the free-stream flow. Analysis of these flow visualization data identified frequencies of flow patterns in the wake shed from the sight, the streamline angle at the sight, and the location where the shed wake crossed the rotor plane. Differences in wake structure were observed between the sight configurations and various skew angles. Analysis of lateral light sheet plane data implied significant vortex structure in the wake of the asymmetric mast-mounted sight in the configuration that produced maximum in-flight vibration. The data showed no significant vortex structure in the wake of the asymmetric and symmetric configurations that produced no increase in in-flight adverse vibration.

Ghee, Terence A.↗
NASA NTRS1994

Comparison between Liouville's theorem and observed latitudinal distributions of trapped ions in the plasmapause region

The presence of anisotropic plasma distributions, trapped at the Earth's magnetic equator, has consequences for the electric field structure which must exist in equilibrium along the magnetic field line. Data from SCATHA and Dynamics Explorer 1 indicate that the core ion distributions at the magnetic equator can be well described as bi-Maxwellian distributions, with a perpendicular temperature an order of magnitude larger than the parallel temperature. A collisionless model is developed for the variation in plasma parameters, following the forms developed by Whipple (1977). If the core electron anisotropy is low, the resulting electric field of approximately 0.1 microV/m is pointed away from the equator. Under these conditions the self-consistent electric field will not overcome the effects of magnetic trapping. The resulting potential distribution results in a local maximum in total plasma density at the equator. Only when the electron distribution is primarily field-aligned can there be a density minimum at the equator. Comparisons are made between this model and the observed variations in DE 1 plasma parameters with latitude.

Olsen, R. C.↗
DOE OSTI2019

Experimental setup combining in situ hard X-ray photoelectron spectroscopy and real-time surface X-ray diffraction for characterizing atomic and electronic structure evolution during complex oxide heterostructure growth

We describe the next-generation system for in situ characterization of a complex oxide thin film and heterostructure growth by pulsed laser deposition (PLD) using synchrotron hard X-rays. The system consists of a PLD chamber mounted on a diffractometer allowing both real-time surface X-ray diffraction (SXRD) and in situ hard X-ray photoelectron spectroscopy (HAXPES). HAXPES is performed in the incident X-ray energy range from 4 to 12 keV using a Scienta EW4000 electron energy analyzer mounted on the PLD chamber fixed parallel with the surface normal. In addition to the standard application mode of HAXPES for disentangling surface from bulk properties, the increased penetration depth of high energy photoelectrons is used for investigation of the electronic structure changes through thin films grown deliberately as variable thickness capping layers. Such heterostructures represent model systems for investigating a variety of critical thickness and dead layer phenomena observed at complex oxide interfaces. In this new mode of operation, in situ HAXPES is used to determine the electronic structure associated with unique structural features identified by real-time SXRD during thin film growth. The system is configured for using both laboratory excitation sources off-line and on-line operation at beamline 33-ID-D at the Advanced Photon Source. We illustrate the performance of the system by preliminary scattering and spectroscopic data on oxygen vacancy ordering induced perovskite-to-brownmillerite reversible phase transformation in La 2/3 Sr 1/3 MnO 3 films capped with oxygen deficient SrTiO 3–δ (100) layers of varying thickness.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

OpenRBC: A Fast Simulator of Red Blood Cells at Protein Resolution

We present OpenRBC, a coarse-grained molecular dynamics code, which is capable of performing an unprecedented in silico experiment—simulating an entire mammal red blood cell lipid bilayer and cytoskeleton as modeled by multiple millions of mesoscopic particles—using a single shared memory commodity workstation. To achieve this, we invented an adaptive spatial-searching algorithm to accelerate the computation of short-range pairwise interactions in an extremely sparse three-dimensional space. The algorithm is based on a Voronoi partitioning of the point cloud of coarse-grained particles, and is continuously updated over the course of the simulation. The algorithm enables the construction of the key spatial searching data structure in our code, i.e., a lattice-free cell list, with a time and space cost linearly proportional to the number of particles in the system. The position and the shape of the cells also adapt automatically to the local density and curvature. The code implements OpenMP parallelization and scales to hundreds of hardware threads. It outperforms a legacy simulator by almost an order of magnitude in time-to-solution and >40 times in problem size, thus providing, to our knowledge, a new platform for probing the biomechanics of red blood cells.

Biophysics↗
DOE OSTI2021

Porting hypre to heterogeneous computer architectures: Strategies and experiences

We report that linear systems are occurring in many applications, and solving them can take a large amount of the total simulation time. The high performance library hypre provides a variety of interfaces and linear solvers, including various multigrid methods, that have achieved good scalability on a variety of homogeneous parallel computer architectures. Heterogeneous architectures with nodes that have both CPUs and accelerators provide new challenges, since they require more fine-grained parallelism and reduced data movement between different memories on a single node as well as across nodes. We will discuss our experiences and strategies to port hypre to heterogeneous computers with accelerators, including the design of a new memory model, the use of abstractions, the BoxLoop macros in the structured and semi-structured interfaces, and the restructuring of algebraic multigrid (AMG) into modular components. We present numerical experiments comparing CPU and GPU performance for several test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Bit-GraphBLAS: Bit-Level Optimizations of Matrix-Centric Graph Processing on GPU

In the graph data structure like adjacency matrix, the connectivity of two nodes can be sufficiently represented using only 1 bit, but they are generally treated as 32-bit full-precision in state-of-the-art graph frameworks to adopt common sparse format such as CSR. Meanwhile, bit-level parallelism has recently be explored to have high-performance potential and low storage requirement on GPUs with dense bit-tiles. To fill the gap, our solution is a hierarchical storage format that contains the bit-indexing base and dense bit-tile units. Inherently, the granularity of the bit-tile is an essential factor in achieving both storage compression and GPU parallelism. How to find a sweet spot that trades off between avoiding sparsity and exploiting is comprehensively researched in this work. In the experiment, we evaluate the proposed storage format and algorithms on modern generation GPUs, including Pascal and Volta, to figure out critical software co-designs in conjunction with existing hardware-specific optimization.

Chen, Jou-An↗
DOE OSTI2021

Enabling particle applications for exascale computing platforms

The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four “sub-motifs”: short-range particle–particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle–particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or “apps”) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable “lessons learned” that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA’s ECP application partner codes are also addressed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Digital tracking loops for a programmable digital modem

In this paper, an analysis and hardware emulation of the tracking loops for a very flexible programmable digital modem (PDM) will be presented. The modem is capable of being programmed for 2, 4, 8, 16-PSK, 16-QAM, MSK, and Offset-QPSK modulation schemes over a range of data rates from 2.34 to 300 Mbps with programmable spectral occupancy from 1.2 to 1.8 times the symbol rate; these operational parameters are executable in burst or continuous mode. All of the critical processing in both the modulator and demodulator is done at baseband with very high-speed digital hardware and memory. Quadrature analog front-ends are used for translation between baseband and the IF center frequency. The modulator is based on a table lookup approach, where precomputed samples are stored in memory and clocked out according to the incoming data pattern. The sample values are predistorted to counteract the effects of the other filtering functions in the link as well as any transmission impairments. The demodulator architecture was adapted from a joint estimator-detector (JED) mathematical analysis. Its structure is applicable to most signalling formats that can be represented in a two-dimensional space. The JED realization uses interdependent, mutually aiding tracking loops with post-detection data feedback. To expedite and provide for more reliable synchronization, initial estimates for these loops are computed in a parallel acquisition processor. The cornerstone of the demodulator realization is the pre-averager received data filter which allows operation over a broad range of data rates without any hardware changes and greatly simplifies the implementation complexity. The emulation results confirmed tracking loop operation over the entire range of operational parameters listed above, as well as the capability of achieving and maintaining synchronization at BER's in excess of 10(exp -1). The emulation results also showed very close agreement with the tracking loop analysis, and validated the resolution apportionment of the various hardware elements in the tracking loops.

Poklemba, John J.↗
NASA NTRS2004

Simple Models of the Spatial Distribution of Cloud Radiative Properties for Remote Sensing Studies

This project aimed to assess the degree to which estimates of three-dimensional cloud structure can be inferred from a time series of profiles obtained at a point. The work was motivated by the desire to understand the extent to which high-frequency profiles of the atmosphere (e.g. ARM data streams) can be used to assess the magnitude of non-plane parallel transfer of radiation in thc atmosphere. We accomplished this by performing an observing system simulation using a large-eddy simulation and a Monte Carlo radiative transfer model. We define the 3D effect as the part of the radiative transfer that isn't captured by one-dimensional radiative transfer calculations. We assess the magnitude of the 3D effect in small cumulus clouds by using a fine-scale cloud model to simulate many hours of cloudiness over a continental site. We then use a Monte Carlo radiative transfer model to compute the broadband shortwave fluxes at the surface twice, once using the complete three-dimensional radiative transfer F(sup 3D), and once using the ICA F (sup ICA); the difference between them is the 3D effect given.

Source record↗