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 199 records · Page 11

Parallel Adaptive Mesh Refinement Library

Parallel Adaptive Mesh Refinement Library (PARAMESH) is a package of Fortran 90 subroutines designed to provide a computer programmer with an easy route to extension of (1) a previously written serial code that uses a logically Cartesian structured mesh into (2) a parallel code with adaptive mesh refinement (AMR). Alternatively, in its simplest use, and with minimal effort, PARAMESH can operate as a domain-decomposition tool for users who want to parallelize their serial codes but who do not wish to utilize adaptivity. The package builds a hierarchy of sub-grids to cover the computational domain of a given application program, with spatial resolution varying to satisfy the demands of the application. The sub-grid blocks form the nodes of a tree data structure (a quad-tree in two or an oct-tree in three dimensions). Each grid block has a logically Cartesian mesh. The package supports one-, two- and three-dimensional models.

Mac-Neice, Peter↗

Automatic Incorporation of Flaws into Cad and FEA Models Based On Measured NDE Data

Automatic incorporation of non-destructive evaluation (NDE) data into CAD and finite element (FE) models for structural analysis is essential for accurate residual-life and strength analyses of composite structures. Currently, two significant but independent challenges exist for incorporating NDE data in such analyses: interpretation of the NDE signals to locate and characterize the structural defects and the incorporation of that information into structural analysis tools. The first challenge will be addressed through linear inversion of calculated synthetic responses generated by a simulation engine to evaluate flaw size, location, and other characteristics. The second challenge will be addressed by building on prior and parallel defect modeling efforts to integrate suitable flaw models into the CAD/finite-element representation of a specimen. Ultrasonic NDE information is usually represented using C-Scans of time-of-flight or intensity information. Other modalities such as thermography and computed tomography (CT) could also be used to characterize the structure. This information will first be projected into material coordinates and then aligned with the specimen’s finite element mesh. The NDE data can then be incorporated by inserting flaw geometry into the CAD and FE analysis model and automatic modification of local material properties and boundary conditions. This process will enable the use of structural analysis tools even for practical modalities such as ultrasound and thermography which may not directly give accurate-enough geometric information to fully define the flaw. The full flaw characteristics will be inferred using prior expert and modeling knowledge to categorize flaw signatures. Outputs for this project will include demonstrated methods for creating structural analysis geometric models that incorporate flaw characterization information; these models can then be used directly in popular commercial FEA software. In addition, the primary deliverable for the project will include a software that will be able to automatically incorporate different kinds of flaws into an existing FE models. The models will be generated in standard file formats to enable direct import into the analysis software.

Adarsh Krishnamurthy↗

Internal structure of shock waves in disparate mass mixtures

The detailed flow structure of a normal shock wave for a gas mixture is investigated using the direct-simulation Monte Carlo method. A variable diameter hard-sphere (VDHS) model is employed to investigate the effect of different viscosity temperature exponents (VTE) for each species in a gas mixture. Special attention is paid to the irregular behavior in the density profiles which was previously observed in a helium-xenon experiment. It is shown that the VTE can have substantial effects in the prediction of the structure of shock waves. The variable hard-sphere model of Bird shows good agreement, but with some limitations, with the experimental data if a common VTE is chosen properly for each case. The VDHS model shows better agreement with the experimental data without adjusting the VTE. The irregular behavior of the light-gas component in shock waves of disparate mass mixtures is observed not only in the density profile, but also in the parallel temperature profile. The strength of the shock wave, the type of molecular interactions, and the mole fraction of heavy species have substantial effects on the existence and structure of the irregularities.

Chung, Chan-Hong↗

A Rover Deployed Ground Penetrating Radar on Mars

Radar is a fundamental tool capable of addressing a variety of geological problems on Mars via collection of data suitable for interpreting variations in surface morphology and reflectivity. Surface-deployed ground penetrating radar (GPR) can help further constrain the geology and structure of the near surface of Mars by directly measuring the range and character of in situ radar properties. In recognition of this potential, a miniaturized, easily modified GPR is being developed for possible deployment on a future Mars rover and will enable definition of radar stratigraphy at high spatial resolution to depths of 10-20 meters. Ongoing development of a Mars impulse GPR with industry partners at Geophysical Survey Systems, Inc., focuses on design and testing of a prototype transducer array (with both high frequency bistatic and low frequency monostatic components) in parallel with fabrication of a low power, mass, and volume control unit. The operational depth of 10-20 meters is geared towards definition of stratigraphy, subsurface blocks, and structure at the decimeter to meter scale that is critical for establishing the geologic setting of the rover. GPR data can also be used to infer the degree of any post-depositional pedogenic alteration or weathering that has subsequently taken place, thereby enabling assessment of pristine versus secondary morphology at the landing site. As is the case for most remote sensing instruments, a GPR may not detect water unambiguously. Nevertheless, any local, near-surface occurrence of liquid water will lead to large, easily detected dielectric contrasts. Moreover, definition of stratigraphy and setting will help in evaluating the history of aqueous activity and where any water might occur and be accessible. Most importantly perhaps, GPR can provide critical context for other rover and orbital instruments/data sets. Hence, GPR deployment along well positioned transects in the vicinity of a lander should enable 3-D mapping of stratigraphy and could serve to guide direct subsurface sampling. Additional information is contained in the original extended abstract.

Grant, J. A.↗

Unsteady interaction of a rotor with a vortex

The unsteady, three-dimensional flow field of a helicopter rotor blade encountering a passing vortex is calculated by solving the Euler/thin layer Navier-Stokes equations by a finite-difference numerical procedure. A prescribed vortex method is adopted to preserve the structure of the interacting vortex. The cases considered for computation correspond to the experimental model rotor test conditions of Caradonna, et al. and consist of parallel and oblique interactions. Comparison of the numerical results with test data show good agreement for both parallel and oblique interactions at subsonic and transonic tip speeds.

Srinivasan, G. R.↗

Mid-Infrared Observational and Theoretical Studies of Star Formation and Early Solar Systems

The first 2 years of this program were used to make mid-IR observations of regions of star formation in the Orion nebula with the UCSD mid-IR camera at the UCSD/University of Minnesota telescope at Mt. Lemmon. These observations attempted to make the first systematic study of an extended region, known to have newly forming stars, and expected to have complex mid-IR emission. We discovered, to our surprise, that most of the thermal emission originated from extended sources rather than from point sources. This interesting observation made the analysis of the data much more complex, since the chop/nod procedures used at these wavelengths produce a differential measurement of the emission in one region compared to that in the adjacent region. Disentangling complex extended emission in such a situation is very difficult. In parallel with this work we were also observing comets in the thermal infrared, the other component of the original proposal. Some spectacular data on the comet Swift-Tuttle was acquired and published. A changing jet structure observed over a 2 week period is described. The rotation period of the comet can be measured at 66 hours. The size of the nucleus can also be estimated (at 30 km) from the observed excess flux from the nucleus. These data have lead to the development of models describing the action of dust particles of differing sizes and composition leaving the nucleus. The spatial distribution of the predicted IR emission has been compared to the observed jet structures, leading to estimates of both particles sizes, relative amounts of silicate vs organic grains, and the amounts of dust emitted in the jets vs isotopic emission.

Jones, Barbara↗

Observations of Low-Latitude Plasma Density Enhancements and their Associated Plasma Drifts

Plasma density structures are frequently encountered in the nighttime low-latitude ionosphere by probes on the Communication/Navigation Outage Forecasting System (C/NOFS) satellite. Of particular interest to us here are plasma density enhancements, which are typically observed +/- 15 deg away from the magnetic equator. The low inclination of the C/NOFS satellite offers an unprecedented opportunity to examine these structures and their associated electric fields and plasma velocities, including their field-aligned components, along an east-west trajectory. Among other observations, the data reveal a clear asymmetry in the velocity structure within and around these density enhancements. Previous observations have shown that the peak change in drift velocity associated with a density enhancement occurs simultaneously both perpendicular and parallel to the magnetic field, while the 1results in this paper show that the peak change in parallel fl ow typically occurs 25-100 km to the east of the peak perpendicular ow. We discuss this and other aspects of the observations in relation to the characteristics of the plasma depletions formed near the magnetic equator detected by the same probes on the C/NOFS satellite and to previous observations and theories.

Klenzing, J. H.↗

RAMP: A fault tolerant distributed microcomputer structure for aircraft navigation and control

RAMP consists of distributed sets of parallel computers partioned on the basis of software and packaging constraints. To minimize hardware and software complexity, the processors operate asynchronously. It was shown that through the design of asymptotically stable control laws, data errors due to the asynchronism were minimized. It was further shown that by designing control laws with this property and making minor hardware modifications to the RAMP modules, the system became inherently tolerant to intermittent faults. A laboratory version of RAMP was constructed and is described in the paper along with the experimental results.

Dunn, W. R.↗

Study on multi-satellite, multi-measurement of the structure of the earth's bow shock

The pulsation model of the earth's bow shock proposed a nonuniform shock having both perpendicular (abrupt, monotonic) and oblique (oscillatory, multigradient) properties simultaneously, depending on local orientation of the shock surface to the interplanetary field B sub sw in parallel planes defined by B sub sw and solar wind velocity. Multiple, concurrent, satellite observations of the shock and solar wind conditions were used. Twenty-six potentially useful intervals of concurrent Explorer 33 and 35 data acquisition were examined, of which six were selected for closer study. In addition, two years of OGO-5 and HEOS 1 magnetometer data were examined for possible conjunctions to these spacecraft having applicable data. One case of clear nonuniformity and several of field-dependent structure were documented. A computational aid, called pulsation index, was developed.

Source record↗

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↗

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↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗