Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “algorithms and data structure”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

The Global Fire Atlas of Individual Fire Size, Duration, Speed and Direction

Natural and human-ignited fires affect all major biomes, altering ecosystem structure, biogeochemical cycles and atmospheric composition. Satellite observations provide global data on spatiotemporal patterns of biomass burning and evidence for the rapid changes in global fire activity in response to land management and climate. Satellite imagery also provides detailed information on the daily or sub-daily position of fires that can be used to understand the dynamics of individual fires. The Global Fire Atlas is a new global dataset that tracks the dynamics of individual fires to determine the timing and location of ignitions, fire size and duration, and daily expansion, fire line length, speed, and direction of spread. Here, we present the underlying methodology and Global Fire Atlas results for 2003–2016 derived from daily moderate-resolution (500 m) Collection 6 MCD64A1 burned-area data. The algorithm identified 13.3 million individual fires over the study period, and estimated fire perimeters were in good agreement with independent data for the continental United States. A small number of large fires dominated sparsely populated arid and boreal ecosystems, while burned area in agricultural and other human-dominated landscapes was driven by high ignition densities that resulted in numerous smaller fires. Long-duration fires in boreal regions and natural landscapes in the humid tropics suggest that fire season length exerts a strong control on fire size and total burned area in these areas. In arid ecosystems with low fuel densities, high fire spread rates resulted in large, short-duration fires that quickly consumed available fuels. Importantly, multiday fires contributed the majority of burned area in all biomass burning regions. A first analysis of the largest, longest and fastest fires that occurred around the world revealed coherent regional patterns of extreme fires driven by large-scale climate forcing. Global Fire Atlas data are publicly available through http://www.globalfiredata.org (last access: 9 August 2018) and https://doi.org/10.3334/ORNLDAAC/1642, and individual fire information and summary data products provide new information for benchmarking fire models within ecosystem and Earth system models, understanding vegetation–fire feedbacks, improving global emissions estimates, and characterizing the changing role of fire in the Earth system.

global fire atlas↗
NASA NTRSDate not supplied

Wall Modeled Large Eddy Simulations for NASA’s Jet Noise Consensus Database of Single-Flow, Round, Convergent Jets

A campaign of wall-modeled large-eddy simulations (WMLES) using structured curvilinear overlapping grids has been performed with the Launch Ascent and Vehicle Aerodynamics(LAVA) computational fluid dynamics (CFD) software to predict jet noise for single-stream axisymmetric round jets. The simulations address the new Prediction Uncertainty Reduction(PUR) technical challenge within the context of NASA’s Commercial Supersonic Technology(CST) project. The goal of PUR is to quantify and reduce uncertainties from scale-resolving simulations to assess noise characteristics of next generation quiet supersonic commercial jets during takeoff and landing conditions where the noise from the exhaust jet dominates. The focus of this effort is to generate a simulation database for single-stream axisymmetric round nozzles at several conditions both for static (no ambient co-flow), which is the focus of this article, and in-flight (flight stream co-flow) conditions, which are beyond the current scope. Nine different flow conditions ranging in jet exit Mach number from 0.38 to 1.0 with nozzle temperature ratios (NTR) from 0.84 to 2.7 have been conducted. Details of the structured overset grids, numerical discretization and wall-model are provided. Near-field comparisons to PIV data show great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the lip line shear-layer. A permeable Ffowcs Williams Hawkings (FWH) surface, enclosing the jet, is used to predict far-field noise from the simulated flow-field. Comparison of CFD predictions to microphone array measurements demonstrate excellent agreement within the resolved frequency range. A systematic under-prediction of far-aft observer angles larger than 150 degrees has been observed across all simulations. We achieved a cost reduction of an order of magnitude for these WMLES compared to an earlier study of this configuration due to algorithmic and software improvements. The accuracy of the results and short turnaround time demonstrate that WMLES within the LAVA framework is a cost-effective approach for jet noise predictions that could soon be incorporated into the design cycle of jet noise reduction technologies.

CST↗
DOE OSTI2021

Characterizing Signatures of Geothermal Exploration Data with Machine Learning Techniques: An Application to the Nevada Play Fairway Analysis

We are introducing machine learning methods to the play fairway analysis to generate geothermal potential maps to support the evaluation of geothermal resource potential and the exploration for undiscovered blind geothermal systems in the Nevada Great Basin region. Our project aims to identify new ways to combine the play fairway data and empirically organize relationships between feature weights and labels in an improved workflow. As a means of doing this, we introduce machine learning methods to evaluate the influence of certain geological and geophysical features/feature sets in predicting geothermal favorability. This report highlights promising approaches based on supervised and unsupervised learning methods. First, we demonstrate a filter method applied to supervised classification modeling. The supervised filter method is based on permutation analysis to evaluate every possible feature combination/drop out scenario and rank feature influence based on the performance variance of supervised classification models. Additionally, we present an unsupervised factor analysis based on principal component analysis coupled with a semi-supervised kmeans clustering algorithm. This analysis allows us to identify the optimal number of groups/clusters for training sites and structural settings to identify feature patterns including correlation, variance, and latent and dominant feature relationships. The results from these methods offer a promising avenue for identifying favorable sources of predictive information to identify the locations of blind geothermal systems and furthering our understanding of complex geothermal feature and label relationships in the Great Basin region and beyond.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

A variant selection framework for genome graphs

Abstract Motivation Variation graph representations are projected to either replace or supplement conventional single genome references due to their ability to capture population genetic diversity and reduce reference bias. Vast catalogues of genetic variants for many species now exist, and it is natural to ask which among these are crucial to circumvent reference bias during read mapping. Results In this work, we propose a novel mathematical framework for variant selection, by casting it in terms of minimizing variation graph size subject to preserving paths of length α with at most δ differences. This framework leads to a rich set of problems based on the types of variants [e.g. single nucleotide polymorphisms (SNPs), indels or structural variants (SVs)], and whether the goal is to minimize the number of positions at which variants are listed or to minimize the total number of variants listed. We classify the computational complexity of these problems and provide efficient algorithms along with their software implementation when feasible. We empirically evaluate the magnitude of graph reduction achieved in human chromosome variation graphs using multiple α and δ parameter values corresponding to short and long-read resequencing characteristics. When our algorithm is run with parameter settings amenable to long-read mapping (α = 10 kbp, δ = 1000), 99.99% SNPs and 73% SVs can be safely excluded from human chromosome 1 variation graph. The graph size reduction can benefit downstream pan-genome analysis. Availability and implementation https://github.com/AT-CG/VF. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Large-scale structure catalogues for cosmological analysis

ABSTRACT We present large-scale structure catalogues from the completed extended Baryon Oscillation Spectroscopic Survey (eBOSS). Derived from Sloan Digital Sky Survey (SDSS) IV Data Release 16 (DR16), these catalogues provide the data samples, corrected for observational systematics, and random positions sampling the survey selection function. Combined, they allow large-scale clustering measurements suitable for testing cosmological models. We describe the methods used to create these catalogues for the eBOSS DR16 Luminous Red Galaxy (LRG) and Quasar samples. The quasar catalogue contains 343 708 redshifts with 0.8 < z < 2.2 over 4808 deg2. We combine 174 816 eBOSS LRG redshifts over 4242 deg2 in the redshift interval 0.6 < z < 1.0 with SDSS-III BOSS LRGs in the same redshift range to produce a combined sample of 377 458 galaxy redshifts distributed over 9493 deg2. Improved algorithms for estimating redshifts allow that 98 per cent of LRG observations result in a successful redshift, with less than one per cent catastrophic failures (Δz > 1000 km s−1). For quasars, these rates are 95 and 2 per cent (with Δz > 3000 km s−1). We apply corrections for trends between the number densities of our samples and the properties of the imaging and spectroscopic data. For example, the quasar catalogue obtains a χ2/DoF = 776/10 for a null test against imaging depth before corrections and a χ2/DoF= 6/8 after. The catalogues, combined with careful consideration of the details of their construction found here-in, allow companion papers to present cosmological results with negligible impact from observational systematic uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Distances to Prompt Effects for a Nuclear Device

In 2010, nuclear weapon effects experts at Sandia National Laboratories (SNL) were asked to provide a quick reference document containing estimated prompt nuclear effects. This report is an update to the 2010 document that includes updated model assumptions. This report addresses only the prompt effects associated with a nuclear detonation (e.g., blast, thermal fluence, and prompt ionizing radiation). The potential medium- and longer-term health effects associated with nuclear fallout are not considered in this report because, in part, of the impracticality of making generic estimates given the high dependency of fallout predictions on the local meteorological conditions at the time of the event. The results included in this report also do not consider the urban environment (e.g., shielding by or collapse of structures) which may affect the extent of prompt effects. It is important to note that any operational recommendations made using the estimates in this report are limited by the generic assumptions considered in the analysis and should not replace analyses made for a specific scenario/device. Furthermore, nuclear effects experts (John Hogan, SNL, and Byron Ristvet, Defense Threat Reduction Agency (DTRA)) have indicated that the accuracy of effects predictions below 0.5 kilotons (kT) or 500 tons nuclear yield have greater uncertainty because of the limited data available for the prompt effects in this regime. The Specialized Hazard Assessment Response Capability (SHARC) effects prediction tool was used for these analyses. Specifically, the NUKE model within SHARC 2021 Version 10.2 was used. NUKE models only the prompt effects following a nuclear detonation. The algorithms for predicting range-to-output data contained within the NUKE model are primarily based on nuclear test effects data. Probits have been derived from nuclear test data and the U.S. Environmental Protection Agency (EPA) protective action guides. Probits relate the probability of a hazard (e.g., fatality or injury) caused by a given insult (e.g., overpressure, thermal fluence, dose level). Several probits have been built into SHARC to determine the fatality and injury associated with a given level of insult. Some of these probits differ with varying yield. Such probits were used to develop the tables and plots in this report.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1999

Design Optimization and Analysis of a Composite Honeycomb Intertank

Intertanks, the structure between tanks of launch vehicles, are prime candidates for weight reduction of rockets. This paper discusses the optimization and detailed follow up analysis and testing of a 96 in. diameter, 77 in. tall intertank. The structure has composite face sheets with an aluminum honeycomb core. The ends taper to a thick built up laminate for a double lap bolted splice joint interface. It is made in 8 full length panels joined with bonded double lap joints. The nominal load is 4000 lb/in. Optimization is by Genetic Algorithm and minimizes weight by varying core thickness, number and orientation of acreage and buildup plies, and the size, number and spacing of bolts. A variety of design cases were run with populations up to 2000 and chromosomes as long as 150 bits. Constraints were buckling; face stresses (normal, shear, wrinkling and dimpling); bolt stress; and bolt hole stresses (bearing, net tension, wedge splitting, shear out and tension/shear out). Analysis is by a combination of elasticity solutions and empirical data. After optimization, a series of coupon tests were performed in conjunction with a rigorous analysis involving a variety of finite element models. This analysis and testing resulted in several small changes to the optimized design. The equation used for hole bearing strength was found to be inadequate, resulting in thicker ends. The core thickness increased 0.05", and potting compound was added in the taper to strengthen the facesheet bond. The intertank has undergone a 250,000 lb limit load test and been mated with a composite liquid hydrogen tank. The tank/intertank unit is being installed in a test stand where it will see 200 thermal/load cycles. Afterwards the intertank will be demated and loaded in compression to failure.

Finckenor, Jeff↗
NASA NTRS2000

A Search for Supernova in Star-Burst Galaxies

The rate at which massive supernovae occur in galaxies can be related directly to the rate of formation of massive stars this rate in turn can be computed from measurements and theoretical models of stellar mss distributions. Other measurements such as chemical composition and the dynamical structure of galaxies can also be used to infer the rate of massive supernovae. Various lines of reasoning using these measurements and theories of star formation have predicted a supernova rate which might be as much as ten times the rate of supernovae that had been observed using optical telescopes. The unseen supernovae were thought to be undetected due to a variety of circumstances, with the primary problem being extinction imposed by dust along the line of observation. This research used new infrared observations of the Infrared Space Observatory (ISO) space craft to explore active galaxies with a technique unhindered by dust extinction. New image processing algorithms for ISO point source analysis and cosmic ray deglitching were developed in order to search through a fairly complete set of the ISO data.

Rank, David M.↗
DOE OSTI2020

Multiphysics for nuclear energy applications using a cohesive computational framework

With the recent development of advanced numerical algorithms, software design, and low-cost high-performance computer hardware, reliance on coupled multiphysics to predict the behavior of complex physical systems is beginning to become standard practice. This is especially true in nuclear energy applications where strong nonlinear interdependencies exist between reactor physics, radiation transport, multi-scale nuclear fuels performance, thermal fluids, etc. Resolving these nonlinear dependencies requires choices in multiphysics software approaches. Two main multiphysics modeling and simulation approaches have emerged. The first is based upon "code coupling" where disparate physics codes of different software design, code languages, and spatial and temporal integration schemes are coupled together with relatively complex data passing interfaces. The second multiphysics software approach is to employ a "cohesive" framework where all physics applications are developed with a common software design, i.e., data structures, syntax, input format, integrated spatial and temporal discretization schemes, etc. In this paper we present the Multiphysics Object-Oriented Simulation Environment (MOOSE) development and runtime framework and describe the framework's cohesive modeling and simulation multiphysics approach. Then, a "cohesive-like" extension of the MOOSE framework is presented where MOOSE-based physics software applications are efficiently coupled to non-MOOSE (external) physics codes to form multiphysics applications using MOOSE's unique interface capabilities. Finally, several examples of MOOSE's cohesive and cohesive-like multiphysics applications will be demonstrated. These multiphysics demonstrations will incorporate both MOOSE-based applications and external codes, including Nek5000, RELAP-7, TRACE, BISON, and Pronghorn.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Deep learning symmetries and their Lie groups, algebras, and subalgebras from first principles

Abstract We design a deep-learning algorithm for the discovery and identification of the continuous group of symmetries present in a labeled dataset. We use fully connected neural networks to model the symmetry transformations and the corresponding generators. The constructed loss functions ensure that the applied transformations are symmetries and the corresponding set of generators forms a closed (sub)algebra. Our procedure is validated with several examples illustrating different types of conserved quantities preserved by symmetry. In the process of deriving the full set of symmetries, we analyze the complete subgroup structure of the rotation groups SO (2), SO (3), and SO (4), and of the Lorentz group S O ( 1 , 3 ) . Other examples include squeeze mapping, piecewise discontinuous labels, and SO (10), demonstrating that our method is completely general, with many possible applications in physics and data science. Our study also opens the door for using a machine learning approach in the mathematical study of Lie groups and their properties.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Dielectric function of InGaAs in the visible

Measurements are reported of the dielectric function of thermodynamically stable In(x)Ga(1-x)As in the composition range 0.3 equal to or less than X = to or less than 0.7. The optically thick samples of InGaAs were made by molecular beam epitaxy (MBE) in the range 0.4 = to or less than X = to or less than 0.7 and by metal-organic chemical vapor deposition (MOCVD) for X = 0.3. The MBE made samples, usually 1 micron thick, were grown on semi-insulating InP and included a strain release structure. The MOCVD sample was grown on GaAs and was 2 microns thick. The dielectric functions were measured by variable angle spectroscopic ellipsometry in the range 1.55 to 4.4 eV. The data was analyzed assuming an optically thick InGaAs material with an oxide layer on top. The thickness of this layer was estimated by comparing the results for the InP lattice matched material, i.e., X = 0.53, with results published in the literature. The top oxide layer mathematically for X = 0.3 and X = 0.53 was removed to get the dielectric function of the bare InGaAs. In addition, the dielectric function of GaAs in vacuum, after a protective arsenic layer was removed. The dielectric functions for X = 0, 0.3, and 0.53 together with the X = 1 result from the literature to evaluate an algorithm for calculating the dielectric function of InGaAs for an arbitrary value of X(0 = to or less than X = to or less than 1) were used. Results of the dielectric function calculated using the algorithm were compared with experimental data.

Alterovitz, S. A.↗
NASA NTRS1990

Dielectric function of InGaAs in the visible

Measurements are reported of the dielectric function of thermodynamically stable In(x)Ga(1-x)As in the composition range 0.3 equal to or less than X = to or less than 0.7. The optically thick samples of InGaAs were made by molecular beam epitaxy (MBE) in the range 0.4 = to or less than X = to or less than 0.7 and by metal-organic chemical vapor deposition (MOCVD) for X = 0.3. The MBE made samples, usually 1 micron thick, were grown on semi-insulating InP and included a strain release structure. The MOCVD sample was grown on GaAs and was 2 microns thick. The dielectric functions were measured by variable angle spectroscopic ellipsometry in the range 1.55 to 4.4 eV. The data was analyzed assuming an optically thick InGaAs material with an oxide layer on top. The thickness of this layer was estimated by comparing the results for the InP lattice matched material, i.e., X = 0.53, with results published in the literature. The top oxide layer mathematically for X = 0.3 and X = 0.53 was removed to get the dielectric function of the bare InGaAs. In addition, the dielectric function of GaAs in vacuum, after a protective arsenic layer was removed. The dielectric functions for X = 0, 0.3, and 0.53 together with the X = 1 result from the literature to evaluate an algorithm for calculating the dielectric function of InGaAs for an arbitrary value of X (0 = to or less than X = to or less than 1) were used. Results of the dielectric function calculated using the algorithm were compared with experimental data.

Alterovitz, S. A.↗
NASA NTRS1991

Distributed parameter estimation for NASA Mini-Mast truss through displacement measurements

Most methods of system identification of large flexible structures by far are based on the lumped parameter approach. Because of the considerable computational burden due to the large number of unknown parameters, distributed parameter approach, which greatly decreases the number of unknowns, has being investigated. In this paper a distributed parameter model for the estimation of modal characteristics of NASA Mini-Mast truss has been formulated. Both Bernoulli-Euler beam and Timoshenko beam equations are used to characterize the lateral bending vibrations of the truss. The measurement of the lateral displacement at the tip of the truss is provided to the maximum likelihood estimator. Closed-form solutions of the partial differential equations and closed-form expressions of the sensitivity functions are derived so that the estimation algorithm is highly efficient. The resulting estimates from test data by using Timoshenko beam model are found to be comparable to those derived from finite element analysis.

Huang, Jen-Kuang↗
NASA NTRS2000

A LEO Satellite Navigation Algorithm Based on GPS and Magnetometer Data

The Global Position System (GPS) has become a standard method for low cost onboard satellite orbit determination. The use of GPS as an attitude and rate sensor has also been developed in the recent past. Additionally, focus has been given to attitude and orbit estimation using the magnetometer, a low cost, reliable sensor. Combining measurements from both GPS and a magnetometer can provide a robust navigation system which takes advantage of the estimation qualities of both measurements. Ultimately a low cost, accurate navigation system can result, potentially eliminating the need for more costly sensors, including gyroscopes. This work presents the development and preliminary testing of a unified navigation algorithm which produces estimates of attitude, angular rate, position, and velocity for a low earth orbit (LEO) spacecraft. The system relies on GPS phase, range, and range rate data as well as magnetometer data. The algorithm used is an extended Kalman filter (EKF) developed to provide LEO attitude, orbit, and rate estimates using magnetometer and sun sensor data. Incorporating sun sensor data into the EKF improved the attitude and rate estimates. For many LEO spacecraft the sun data is available during only a portion of the orbit. However, GPS data is available continuously throughout the orbit. GPS can produce accurate orbit estimates and combining GPS and magnetometer data improves the attitude and rate estimates. The magnetometer based EKF can converge from large initial errors in position, velocity, and 3 attitude . Combining the magnetometer and GPS data into a single EKF will provide a more robust and accurate system. The EKF is based on an existing EKF. The GPS measurement models for phase, range, and range rate are incorporated into the existing structure of the filter. The original EKF produced the orbit estimates in terms of Keplerian elements. Due to the nature of the GPS measurements and ease of computation, the orbit estimates are converted to the Cartesian position and velocity. The measurement model for the magnetometer is adjusted for this change in the state and the measurements. Preliminary test results based on simulated GPS and magnetometer data are included.

Deutschmann, Julie↗
NASA NTRS2006

DCS-Neural-Network Program for Aircraft Control and Testing

A computer program implements a dynamic-cell-structure (DCS) artificial neural network that can perform such tasks as learning selected aerodynamic characteristics of an airplane from wind-tunnel test data and computing real-time stability and control derivatives of the airplane for use in feedback linearized control. A DCS neural network is one of several types of neural networks that can incorporate additional nodes in order to rapidly learn increasingly complex relationships between inputs and outputs. In the DCS neural network implemented by the present program, the insertion of nodes is based on accumulated error. A competitive Hebbian learning rule (a supervised-learning rule in which connection weights are adjusted to minimize differences between actual and desired outputs for training examples) is used. A Kohonen-style learning rule (derived from a relatively simple training algorithm, implements a Delaunay triangulation layout of neurons) is used to adjust node positions during training. Neighborhood topology determines which nodes are used to estimate new values. The network learns, starting with two nodes, and adds new nodes sequentially in locations chosen to maximize reductions in global error. At any given time during learning, the error becomes homogeneously distributed over all nodes.

Jorgensen, Charles C.↗
DOE OSTI2023

Differentiable physics-enabled closure modeling for Burgers’ turbulence

Abstract Data-driven turbulence modeling is experiencing a surge in interest following algorithmic and hardware developments in the data sciences. We discuss an approach using the differentiable physics paradigm that combines known physics with machine learning to develop closure models for Burgers’ turbulence. We consider the one-dimensional Burgers system as a prototypical test problem for modeling the unresolved terms in advection-dominated turbulence problems. We train a series of models that incorporate varying degrees of physical assumptions on an a posteriori loss function to test the efficacy of models across a range of system parameters, including viscosity, time, and grid resolution. We find that constraining models with inductive biases in the form of partial differential equations that contain known physics or existing closure approaches produces highly data-efficient, accurate, and generalizable models, outperforming state-of-the-art baselines. Addition of structure in the form of physics information also brings a level of interpretability to the models, potentially offering a stepping stone to the future of closure modeling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Via Machinae 2.0 : Full-sky, model-agnostic search for stellar streams in Gaia DR2

ABSTRACT We present an update to Via Machinae, an automated stellar stream-finding algorithm based on the deep learning anomaly detector ANODE. Via Machinae identifies stellar streams within Gaia, using only angular positions, proper motions, and photometry, without reference to a model of the Milky Way potential for orbit integration or stellar distances. This new version, Via Machinae 2.0, includes many improvements and refinements to nearly every step of the algorithm, that altogether result in more robust and visually distinct stream candidates than our original formulation. In this work, we also provide a quantitative estimate of the false positive rate of Via Machinae 2.0 by applying it to a simulated Gaia-mock catalogue based on galaxia, a smooth model of the Milky Way that does not contain substructure or stellar streams. Finally, we perform the first full-sky search for stellar streams with Via Machinae 2.0, identifying 102 streams at high significance within the Gaia Data Release 2, of which only 10 have been previously identified. While follow-up observations for further confirmation are required, taking into account the false positive rate presented in this work, we expect approximately 90 of these stream candidates to correspond to real stellar structures.

Astronomy & Astrophysics↗
NASA NTRS1992

Mapping unstructured grid computations to massively parallel computers

Investigated here is this mapping problem: assign the tasks of a parallel program to the processors of a parallel computer such that the execution time is minimized. First, a taxonomy of objective functions and heuristics used to solve the mapping problem is presented. Next, we develop a highly parallel heuristic mapping algorithm, called Cyclic Pairwise Exchange (CPE), and discuss its place in the taxonomy. CPE uses local pairwise exchanges of processor assignments to iteratively improve an initial mapping. A variety of initial mapping schemes are tested and recursive spectral bipartitioning (RSB) followed by CPE is shown to result in the best mappings. For the test cases studied here, problems arising in computational fluid dynamics and structural mechanics on unstructured triangular and tetrahedral meshes, RSB and CPE outperform methods based on simulated annealing. Much less time is required to do the mapping and the results obtained are better. Compared with random and naive mappings, RSB and CPE reduce the communication time two fold for the test problems used. Finally, we use CPE in two applications on a CM-2. The first application is a data parallel mesh-vertex upwind finite volume scheme for solving the Euler equations on 2-D triangular unstructured meshes. CPE is used to map grid points to processors. The performance of this code is compared with a similar code on a Cray-YMP and an Intel iPSC/860. The second application is parallel sparse matrix-vector multiplication used in the iterative solution of large sparse linear systems of equations. We map rows of the matrix to processors and use an inner-product based matrix-vector multiplication. We demonstrate that this method is an order of magnitude faster than methods based on scan operations for our test cases.

Hammond, Steven Warren↗