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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Single-mode projection filters for modal parameter identification for flexible structures

Single-mode projection filters are developed for eigensystem parameter identification from both analytical results and test data. Explicit formulations of these projection filters are derived using the orthogonal matrices of the controllability and observability matrices in the general sense. A global minimum optimization algorithm is applied to update the filter parameters by using the interval analysis method. The updated modal parameters represent the characteristics of the test data. For illustration of this new approach, a numerical simulation for the MAST beam structure is shown by using a one-dimensional global optimization algorithm to identify modal frequencies and damping. The projection filters are practical for parallel processing implementation.

Huang, Jen-Kuang↗
NASA NTRS2001

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗
DOE OSTI2017

Diagnosing collisionless energy transfer using field–particle correlations: gyrokinetic turbulence

Determining the physical mechanisms that extract energy from turbulent fluctuations in weakly collisional magnetized plasmas is necessary for a more complete characterization of the behaviour of a variety of space and astrophysical plasmas. Such a determination is complicated by the complex nature of the turbulence as well as observational constraints, chiefly that in situ measurements of such plasmas are typically only available at a single point in space. Recent work has shown that correlations between electric fields and particle velocity distributions constructed from single-point measurements produce a velocity-dependent signature of the collisionless damping mechanism. We extend this work by constructing field–particle correlations using data sets drawn from single points in strongly driven, turbulent, electromagnetic gyrokinetic simulations to demonstrate that this technique can identify the collisionless mechanisms operating in such systems. The velocity-space structure of the correlation between proton distributions and parallel electric fields agrees with expectations of resonant mechanisms transferring energy collisionlessly in turbulent systems. This work motivates the eventual application of field–particle correlations to spacecraft measurements in the solar wind, with the ultimate goal to determine the physical mechanisms that dissipate magnetized plasma turbulence.

Physics↗
DOE OSTI2020

A massively parallel and scalable multi-CPU material point method

Harnessing the power of modern multi-GPU architectures, we present a massively parallel simulation system based on the Material Point Method (MPM) for simulating physical behaviors of materials undergoing complex topological changes, self-collision, and large deformations. Our system makes three critical contributions. First, we introduce a new particle data structure that promotes coalesced memory access patterns on the GPU and eliminates the need for complex atomic operations on the memory hierarchy when writing particle data to the grid. Second, we propose a kernel fusion approach using a new Grid-to-Particles-to-Grid (G2P2G) scheme, which efficiently reduces GPU kernel launches, improves latency, and significantly reduces the amount of global memory needed to store particle data. Finally, we introduce optimized algorithmic designs that allow for efficient sparse grids in a shared memory context, enabling us to best utilize modern multi-GPU computational platforms for hybrid Lagrangian-Eulerian computational patterns. We demonstrate the effectiveness of our method with extensive benchmarks, evaluations, and dynamic simulations with elastoplasticity, granular media, and fluid dynamics. In comparisons against an open-source and heavily optimized CPU-based MPM codebase [Fang et al. 2019] on an elastic sphere colliding scene with particle counts ranging from 5 to 40 million, our GPU MPM achieves over 100x per-time-step speedup on a workstation with an Intel 8086K CPU and a single Quadro P6000 GPU, exposing exciting possibilities for future MPM simulations in computer graphics and computational science. Moreover, compared to the state-of-the-art GPU MPM method [Hu et al. 2019a], we not only achieve 2x acceleration on a single GPU but our kernel fusion strategy and Array-of-Structs-of-Array (AoSoA) data structure design also generalizes to multi-GPU systems. Our multi-GPU MPM exhibits near-perfect weak and strong scaling with 4 GPUs, enabling performant and large-scale simulations on a 10243 grid with close to 100 million particles with less than 4 minutes per frame on a single 4-GPU workstation and 134 million particles with less than 1 minute per frame on an 8-GPU workstation.

Wang, Xinlei↗
DOE OSTI2017

Extreme-Scale De Novo Genome Assembly

De novo whole genome assembly reconstructs genomic sequence from short, overlapping, and potentially erroneous DNA segments and is one of the most important computations in modern genomics. This work presents HipMER, a high-quality end-to-end de novo assembler designed for extreme scale analysis, via efficient parallelization of the Meraculous code. Genome assembly software has many components, each of which stresses different components of a computer system. This chapter explains the computational challenges involved in each step of the HipMer pipeline, the key distributed data structures, and communication costs in detail. We present performance results of assembling the human genome and the large hexaploid wheat genome on large supercomputers up to tens of thousands of cores.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1988

System architectures for PROLOG execution

The principal features of PROLOG are reviewed with particular reference to the use of the language in the development of expert systems. Attention is given to the language syntax, data structures and recursion, constructs that alter control, predicates that modify the database, predicates for input and output, and metalogical predicates. The parts of the Warren Abstract Machine are then described, and it is shown how it executes programs. Parallel architectures proposed for expert systems are examined, and the principal varieties of parallelism available in PROLOG and examples of architectures exploiting them are discussed.

Gooley, Markian M.↗
NASA NTRS1991

Lanczos eigensolution method for high-performance computers

The theory, computational analysis, and applications are presented of a Lanczos algorithm on high performance computers. The computationally intensive steps of the algorithm are identified as: the matrix factorization, the forward/backward equation solution, and the matrix vector multiples. These computational steps are optimized to exploit the vector and parallel capabilities of high performance computers. The savings in computational time from applying optimization techniques such as: variable band and sparse data storage and access, loop unrolling, use of local memory, and compiler directives are presented. Two large scale structural analysis applications are described: the buckling of a composite blade stiffened panel with a cutout, and the vibration analysis of a high speed civil transport. The sequential computational time for the panel problem executed on a CONVEX computer of 181.6 seconds was decreased to 14.1 seconds with the optimized vector algorithm. The best computational time of 23 seconds for the transport problem with 17,000 degs of freedom was on the the Cray-YMP using an average of 3.63 processors.

Bostic, Susan W.↗
NASA NTRS1976

Architecture and data processing alternatives for the TSE computer. Volume 3: Execution of a parallel counting algorithm using array logic (Tse) devices

A parallel algorithm for counting the number of logic-l elements in a binary array or image developed during preliminary investigation of the Tse concept is described. The counting algorithm is implemented using a basic combinational structure. Modifications which improve the efficiency of the basic structure are also presented. A programmable Tse computer structure is proposed, along with a hardware control unit, Tse instruction set, and software program for execution of the counting algorithm. Finally, a comparison is made between the different structures in terms of their more important characteristics.

Metcalfe, A. G.↗
DOE OSTI2018

Structural disorder and elementary magnetic properties of triangular lattice ErMgGaO 4 single crystals

The single crystal growth, structure, and basic magnetic properties of ErMgGaO 4 are reported. The structure consists of triangular layers of magnetic ErO 6 octahedra separated by a double layer of randomly occupied non-magnetic (Ga,Mg)O 5 bipyramids. The Er atoms are positionally disordered. Magnetic measurements parallel and perpendicular to the c axis of a single crystal reveal dominantly antiferromagnetic interactions, with a small degree of magnetic anisotropy. A weighted average of the directional data suggests an antiferromagnetic Curie Weiss temperature of approximately -30 K. Below 10 K the temperature dependences of the inverse susceptibilities in the in-plane and perpendicular-to plane directions are parallel, indicative of an isotropic magnetic moment at low temperatures. In conclusion, no sign of magnetic ordering is observed above 1.8 K, suggesting that ErMgGaO 4 is a geometrically frustrated magnet.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Parallel Grid Manipulations in Earth Science Calculations

The National Aeronautics and Space Administration (NASA) Data Assimilation Office (DAO) at the Goddard Space Flight Center is moving its data assimilation system to massively parallel computing platforms. This parallel implementation of GEOS DAS will be used in the DAO's normal activities, which include reanalysis of data, and operational support for flight missions. Key components of GEOS DAS, including the gridpoint-based general circulation model and a data analysis system, are currently being parallelized. The parallelization of GEOS DAS is also one of the HPCC Grand Challenge Projects. The GEOS-DAS software employs several distinct grids. Some examples are: an observation grid- an unstructured grid of points at which observed or measured physical quantities from instruments or satellites are associated- a highly-structured latitude-longitude grid of points spanning the earth at given latitude-longitude coordinates at which prognostic quantities are determined, and a computational lat-lon grid in which the pole has been moved to a different location to avoid computational instabilities. Each of these grids has a different structure and number of constituent points. In spite of that, there are numerous interactions between the grids, e.g., values on one grid must be interpolated to another, or, in other cases, grids need to be redistributed on the underlying parallel platform. The DAO has designed a parallel integrated library for grid manipulations (PILGRIM) to support the needed grid interactions with maximum efficiency. It offers a flexible interface to generate new grids, define transformations between grids and apply them. Basic communication is currently MPI, however the interfaces defined here could conceivably be implemented with other message-passing libraries, e.g., Cray SHMEM, or with shared-memory constructs. The library is written in Fortran 90. First performance results indicate that even difficult problems, such as above-mentioned pole rotation- a sparse interpolation with little data locality between the physical lat-lon grid and a pole rotated computational grid- can be solved efficiently and at the GFlop/s rates needed to solve tomorrow's high resolution earth science models. In the subsequent presentation we will discuss the design and implementation of PILGRIM as well as a number of the problems it is required to solve. Some conclusions will be drawn about the potential performance of the overall earth science models on the supercomputer platforms foreseen for these problems.

Sawyer, W.↗
NASA NTRS1990

Disk mass densities in edge-on spiral galaxies

Very Large Array (VLA) observations of the neutral hydrogen (HI) gas in two nearby edge-on spirals (NGC 4565 and NGC 891) successfully resolve the thickness of the gas layers in both disks over a wide range in radii. The combination of B, C, and D array data produces a 4 arcsec (approx. 200 pc) beam and 21 km s(exp -1) velocity resolution, combined with sensitivity to structures as large as 18 arcmin (approx. 54 kpc). These observations directly constrain the mid-plane disk mass densities, under the assumption of an equilibrium between the thermal pressure of the gas and the gravitational attraction of the disk. The results of a preliminary analysis are given regarding the z-velocity dispersion of the gas, the mass-to-light ratio of the disk in NGC 4565, and the roles of atomic and molecular gases. The data also allow a detailed study of the HI in these galaxies; in general their brightness temperature distributions seem similar to that in the Milky Way. Both galaxies show asymmetric HI extensions beyond the optical disk. In NGC 4565 the extension is a surprisingly abrupt warp, which may bend back to parallel the galactic plane; the velocity structure implies the warp is continuous around the disk.

Rupen, Michael P.↗
DOE OSTI2022

Planar 2D wireframe DNA origami

Two-dimensional (2D) DNA origami is widely used for applications ranging from excitonics to single-molecule biophysics. Conventional, single-layer 2D DNA origami exhibits flexibility and curvature in solution; however, that may limit its suitability as a 2D structural template. In contrast, 2D wireframe DNA origami rendered with six-helix bundle edges offers local control over duplex orientations with enhanced in-plane rigidity. Here, we investigate the 3D structure of these assemblies using cryo–electron microscopy (cryo-EM). 3D reconstructions reveal a high degree of planarity and homogeneity in solution for polygonal objects with and without internal mesh, enabling 10-Å resolution for a triangle. Coarse-grained simulations were in agreement with cryo-EM data, offering molecular structural insight into this class of 2D DNA origami. Our results suggest that these assemblies may be valuable for 2D material applications and geometries that require high structural fidelity together with local control over duplex orientations, rather than parallel duplex assembly.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1979

Optical memories in digital computing

High capacity optical memories with relatively-high data-transfer rate and multiport simultaneous access capability may serve as basis for new computer architectures. Several computer structures that might profitably use memories are: a) simultaneous record-access system, b) simultaneously-shared memory computer system, and c) parallel digital processing structure.

Alford, C. O.↗
NASA NTRS1973

The Mariner 9 celestial mechanics experiment

There are two main objectives in the Mars gravity field analysis. The first is to generate a picture of the gravity field which can be used in conjunction with other types of data to elucidate the interior, surface structure, and history of Mars. The second is to aid the orbit determination aspect of the relativity experiment. In an effort to extract from the data a Mars gravity model of reasonable fidelity, various spherical harmonic models have been made. A parallel effort using a surface mass representation involves the differentiation of Doppler residuals for acceleration contours.

Lorell, J.↗
NASA NTRSDate not supplied

Application of the NASA Multiscale Analysis Tool: Multiscale Integration and Interoperability

The NASA Multiscale Analysis Tool (NASMAT) was developed recently to allow a wide variety of multiscale analysis problems to be effectively and efficiently solved. The architecture of NASMAT was established specifically to enable parallelized, “plug-and-play” functionality to reduce the complexity associated with adding new features to the code in the future and to allow end users to rapidly implement and evaluate user-defined capabilities. Additionally, the tool utilizes recursive data structures and subroutines to allow for an arbitrary number of length scales when performing multiscale analyses of heterogeneous materials. These features permit the rapid integration of user-defined capabilities (e.g., a material model, micromechanics approach, or failure theory) at all stages within a NASMAT calculation while leveraging built-in techniques where needed. Additionally, these features allow NASMAT to both be called from an external program as well as call an external program. This paper specifically focuses on the multiscale integration and interoperability of NASMAT with other analysis techniques through an illustrative, multiscale analysis of a 3D woven polymer matrix composite (PMC).

NASMAT↗
DOE OSTI2022

Discovering causal structure with reproducing-kernel Hilbert space ε -machines

We merge computational mechanics’ definition of causal states (predictively equivalent histories) with reproducing-kernel Hilbert space (RKHS) representation inference. The result is a widely applicable method that infers causal structure directly from observations of a system’s behaviors whether they are over discrete or continuous events or time. A structural representation—a finite- or infinite-state kernel ϵ-machine—is extracted by a reduced-dimension transform that gives an efficient representation of causal states and their topology. In this way, the system dynamics are represented by a stochastic (ordinary or partial) differential equation that acts on causal states. We introduce an algorithm to estimate the associated evolution operator. Paralleling the Fokker–Planck equation, it efficiently evolves causal-state distributions and makes predictions in the original data space via an RKHS functional mapping. We demonstrate these techniques, together with their predictive abilities, on discrete-time, discrete-value infinite Markov-order processes generated by finite-state hidden Markov models with (i) finite or (ii) uncountably infinite causal states and (iii) continuous-time, continuous-value processes generated by thermally driven chaotic flows. The method robustly estimates causal structure in the presence of varying external and measurement noise levels and for very high-dimensional data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Parallel Tensor Decompositions for Massive Heterogeneous, Incomplete Data (Final Report: LDRD Project 199986).

We extended the fundamental capabilities of tensor decomposition to a broader range of problems - handling new data types and larger problems. This has implications for data analysis across a range of applications in sensor monitoring, cybersecurity, treaty verification, signal processing, etc. Identifies latent structure within data, enabling anomaly detection, process monitoring, scientific discovery.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Soil and Water Chemistry and Trace Metal Extractability and Speciation in Wetland Soils from Illinois and South Carolina and Stream Sediments from Tennessee

Dataset revised on October 15, 2021. This revision adds sulfur and iron X-ray absorption near-edge structure spectra for the wetland soils and stream sediments from the field areas. It also renames the sample locations in a way that is more intuitive to readers of the companion paper that is under review. Finally, the data filenames and organization have been updated in their labeling to parallel the data sources in the associated paper. The abstract text and methods were also revised to reflect the data that was added to the dataset.Trace metals are essential for microbially-mediated biogeochemical processes occurring in anoxic wetland soils and stream bed sediments, such as denitrification, methanogenesis, and mercury methylation. Low availability of these elements may potentially inhibit key components of anaerobic carbon and nitrogen cycling and contaminant transformation. The solid-phase speciation of trace metals likely plays an important role in controlling their bioavailability. Metal speciation is well studied in contaminated soils and sediments as well as those naturally elevated in trace metals. However, less is known regarding the chemical forms of trace metals in systems having concentrations similar to geological background levels, the very settings where metal limitations may be most prevalent. We have investigated trace metal concentrations, extractability, and solid-phase speciation in three freshwater subsurface aquatic systems: marsh wetland soils, riparian wetland soils, and the sediments of a streambed.Data are provided for marsh wetland soils at Argonne National Laboratory, riparian wetland soils in the Tims Branch watershed at Savannah River National Laboratory, and stream bed sediments from East Fork Poplar Creek near Oak Ridge National Laboratory. Soil and sediment elemental abundances, mineralogy, and extractable nutrients as well as dissolved major elements, anions, trace metals, and nutrients in the overlying surface waters are provided. In addition, the results of sequential chemical extraction for the trace metals cobalt, nickel, copper, and zinc from the soils and sediment are reported as well as X-ray absorption near-edge structure (XANES) spectra in these materials are reported. To aid interpretation of these data, XANES spectra of sulfur in the soils and sediments as well as both XANES and extended X-ray absorption fine structure (EXAFS) spectra of iron in these materials are reported. The data package also includes the XANES spectra of reference standards and a potential interferent in the measurements. All data are provided in text-based CSV format with header sections indicating the data contained in each file and the corresponding units. Note that "u" is used in place of Greek lower case mu to indicate the micro prefix on units.

54 ENVIRONMENTAL SCIENCES↗