Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fixed grid”

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 577 records · Page 32

The Yale–Potsdam Stellar Isochrones

We introduce the Yale–Potsdam Stellar Isochrones (YaPSI), a new grid of stellar evolution tracks and isochrones of solar-scaled composition. In an effort to improve the Yonsei–Yale database, special emphasis is placed on the construction of accurate low-mass models (M{sub ∗}<0.6 M{sub ⊙}), and in particular on their mass–luminosity and mass–radius relations, both crucial for characterizing exoplanet-host stars, and, in turn, their planetary systems. The YaPSI models cover the mass range 0.15–5.0 M{sub ⊙} densely enough to permit detailed interpolation in mass, and the metallicity and helium abundance ranges [Fe/H] = −1.5 to +0.3 and Y {sub 0} = 0.25–0.37 are specified independently of each other (i.e., no fixed ΔY/ΔZ relation is assumed). The evolutionary tracks are calculated from the pre-main sequence up to the tip of the red giant branch. The isochrones, with ages between 1 Myr and 20 Gyr, provide UBVRI colors in the Johnson–Cousins system, and JHK colors in the homogenized Bessell and Brett system, derived from two different semi-empirical T{sub eff}–color calibrations from the literature. We also provide utility codes, such as an isochrone interpolator, in age, metallicity, and helium content, and an interface of the tracks with an open-source Monte Carlo Markov-Chain tool for the analysis of individual stars. Finally, we present comparisons of the YaPSI models with the best empirical mass–luminosity and mass–radius relations available to date, as well as isochrone fitting of well-studied stellar clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Validation Experiments for Turbulent Separation over an Axisymmetric Body of Revolution

Historically, the flow physics involved with most turbulent separated flows have presented fundamental challenges to validating numerical approaches. As recognized by the CFD Vision 2030 study commissioned by NASA, validation of Reynolds-averaged Navier-Stokes (RANS) models and other scale-resolving methods for turbulent separated flow requires data from advanced, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, a new test platform, referred to as the NASA Axisymmetric Afterbody, was designed to obtain detailed measurements of the flow field undergoing a smooth, adverse pressure gradient induced separation for a fixed Reynolds number, Re = 180,000. The parametric body offers a range of flow states progressing from fully attached, to incipient separation, and finally to small-scale separated flow based on variable afterbody geometries. In an initial effort to evaluate RANS turbulence model capabilities, the present configurations of the axisymmetric model host a mild adverse pressure gradient over the contoured boattail section, inducing incipient turbulent separation, as well as a slightly larger adverse pressure gradient, inducing a small-scale region of turbulent separation. Experiments include steady pressure measurements and 2-D PIV to provide the preliminary dataset for simulation studies, which examine the effect of variable grid domains and RANS turbulence models. This is done in an effort to understand and evaluate the critical variability between solutions for the present model configurations. Results indicate potential discrepancies may be due to the effect of the square tunnel test section walls, relatively large blockage ratio, and slight variability in reference parameters. Ongoing work will focus on higher fidelity experimental campaigns to obtain surface flow visualizations and Stereo-Particle Image Velocimetry (SPIV) to deliver higher spatial resolution of the three-dimensional flow field to aid turbulence modelers.

Validation↗

Synthetic, Realistic Transmission and Distribution Co-Simulation for Voltage Control Benchmarking: Preprint

As distributed energy resources play an increasing role in generation and grid services, control and market paradigms will need to be tested with models that incorporate transmission and distribution systems. The test systems will need to be realistic and accessible for researchers to form conclusions about realistic behaviors and maintain comparability and repeatability of results. These test systems also must have high resolution models of sufficient size for realistic distribution and transmission interaction. This paper demonstrates the tightly coupled co-simulation of realistic synthetic distribution and transmission models for the Austin metropolitan area and how this framework can be used to analyze resultant behaviors from distributed resources and their controls. The comparison of a distribution system run with fixed feed-in voltage is compared to the closely tied distribution-transmission co-simulation. It is determined that even without advanced voltage controls, the feed-in voltage response of the transmission system shows a lower level of voltage excursions. Therefore, control methods for voltage stability should be benchmarked on realistic transmission and distribution co-simulations for more accurate assessment of voltage profiles.

27 ARPA - Advanced Research Projects Agency-Energy↗

Computational and Experimental Validation of Fractal-Fin, Dimpled Solar Heat Collector

A 3 foot x 3 foot x 3 foot aluminum solar collector was manufactured using computer numerical control. The interior of the device included six triangular dimpled fins for enhanced heat transfer. The interior vertical wall on the south side was also dimpled. The solar collector working fluid was based on water, and the collector consisted solely of passive heat transfer mechanisms (no moving parts), making it ideal for off-the-grid and rural applications. Two types of heat transfer experiments were conducted. One experiment had external flat heaters attached on the top and the front side, while the other four sides were insulated. Except for the bottom surface, the second experiment had all its exterior surfaces sprayed with black solar paint to collect as much solar heat as possible. Temperature data as a function of time was collected using 14 thermocouples spread strategically throughout the solar collector. In addition, computational fluid dynamics (CFD) simulations were conducted using the dynamic Smagorinsky large eddy simulation turbulence model. The first simulation considered that both the top and front surfaces were exposed to a fixed temperature of 313.7 K (105 °F), while the remaining four surfaces were insulated. For the second simulation, all conditions were the same, except that the temperature for both heated surfaces was raised to 350 K (170.3 °F). The two temperatures are expected to bound the solar collector operational temperature during the late- Spring, Summer, and early-Fall months. The solar collector design, experimental data, CFD output, and a discussion of five manufacturing approaches and costs are documented in this report.

14 SOLAR ENERGY↗

Weakly Supervised Event Classification Using Imperfect Real-world PMU Data with Scarce Labels

This paper studies event classification using imperfect real-world phasor measurement unit (PMU) data with scarce event types (labels). By investigating the real-world PMU data, it is observed that most real-world PMU data's event type is unknown, which makes it challenging to directly use such dataset to build event classifiers as existing classification techniques require high-quality training data with known event type (i.e., label). To address this challenge, a weakly supervised learning based event classification approach is developed, which can use noisy and low-quality PMU data for the training. First, data quality issues are fixed using data preprocessing techniques and then event features are constructed from the PMU data. Using these features, a series of labeling functions are learnt to generate initial estimates of the labels of large amounts of unlabeled PMU data. As the labeling functions are learnt using the same data with scarce labels, the label estimates from the labeling functions can be correlated, noisy, and bias. To enhance these initial estimates, a generative model is developed to characterize the dependencies among the estimated labels, based on which better labels are obtained for training event classifiers. Numerical experiments using the real-world dataset from the Western Interconnection of the U.S. power transmission grid show that the proposed weakly supervised event classifier trained using the dataset with only 5% labeled data can achieve 78.4% classification accuracy.

Liu, Yunchuan↗

Lithium in the most extreme halo stars - Trends with metallicity

A grid of extremely metal-poor stellar evolutionary tracks (Z = 0.00001) is constructed to compare the predicted lithium abundance depletion to observations of the lithium abundance in correspondingly metal-poor stars (Fe/H = -3.3), and, combined with an earlier work to investigate trends of lithium isochrones with metal abundance. Lithium abundances in the most extreme halo stars hint at the possibility of small downward slope toward lower T(eff) in the flat Li-T(eff) relation (the Spite lithium plateau) observed in slightly less extreme stars; the models also predict this. In the cool star region of the Li-T(eff) plane, higher metallicity lithium isochromes are cooler at fixed lithium abundance, consistent with the observations. Current data thus agree well with the lithium trends from standard model.

Deliyannis, Constantine P.↗

Identifying Potential Candidates for Renewable Energy Zones (REZs) in Bangladesh

Bangladesh faces several hurdles to achieving its renewable energy objectives, such as land availability and transmission congestion. Renewable Energy Zones (REZ), which are geographic areas with high-quality utility-scale renewable resources, suitable land topography, and commercial interest, can help address these challenges and support long-term generation and transmission planning. This study focuses on solar PV (fixed-tilt) and onshore wind, leveraging recently developed high temporal and spatial resolution resource data. In the moderate land exclusion scenario, large "study areas" are identified in Bangladesh with capacity factors in the top 25% for the entire country - 5 for wind and 8 for solar. Within these study areas, 19 candidates for REZ are identified based on the overlap between these study areas and upazilas (i.e., administrative subdivisions) containing economic development zones. Four of the identified candidate zones are opportunities for priority development, given the combination of strong wind and solar resources and the presence of economic zones. Furthermore, the geographic diversity of wind and solar resources in Bangladesh could help increase grid resilience by not concentrating all renewable energy development in the same region. Finally, pairing REZ with economic zones can bolster economic development, take advantage of large electricity demand, and leverage existing infrastructure investments.

Bangladesh↗

Ordovician-Cambrian Units: Hierarchical Evaluation of Geologic Carbon Storage Resource Estimates

The Indiana Geological and Water Survey (IGWS) led subtask 1.1 to assess the regional distribution and estimate the storage capacity of Ordovician-Cambrian stratigraphic units located within the partnership region. A comprehensive data set of wireline logs and petrophysical information was used to generate these interpretations. These data include core analysis for porosity and permeability, mercury injection capillary pressure (MICP), and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for three potential storage reservoirs(limestone and dolostone from the Upper Ordovician Trenton Limestone/Black River Group and equivalent units, the Middle Ordovician St. Peter Sandstone, and primary target reservoir rocks of the Lower Ordovician and Upper Cambrian Knox Supergroup and equivalent units) calculated using six methodologies: (1) a fixed value of porosity of 10 percent in all units evaluated; (2) a unique average porosity (per well) from wireline-derived porosity (neutron, sonic, and/or density porosity for each unit); (3) porosity values from core analysis; (4) a depth-dependent porosity model (Knox Supergroup only); (5) porosity based on a model based on petrophysical facies; and (6) SREs using National Energy Technology Laboratory’s CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN beta V2). All methods used the same values for thickness for each unit. However, the areal extent of each assessment was limited by the data available for each method. Estimated volumes were calculated in 1-by-1 kilometer grid cells and summarized as county and total stratigraphic unit volumes. The resultant SREs mass are displayed using boxplots, which allow for comparing data statistics (mean values and variability) between methods. Differences observed in SRE results from the six methods are mainly attributable to differences in the data and conceptual models used to interpret or estimate porosity in each method. Based on this systematic variability between methods, it is inferred that methods 1, 4, and 6 are best used for regional-scale reconnaissance estimates of storage capacity while methods 2, 3, and 5 are more appropriate for local scales where more data is required. All estimates are data-density dependent and different methods require different amounts of data for reasonable assessments. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗

The compression–error trade-off for large gridded data sets

The netCDF-4 format is widely used for large gridded scientific data sets and includes several compression methods: lossy linear scaling and the non-lossy deflate and shuffle algorithms. Many multidimensional geoscientific data sets exhibit considerable variation over one or several spatial dimensions (e.g., vertically) with less variation in the remaining dimensions (e.g., horizontally). On such data sets, linear scaling with a single pair of scale and offset parameters often entails considerable loss of precision. We introduce an alternative compression method called "layer-packing" that simultaneously exploits lossy linear scaling and lossless compression. Layer-packing stores arrays (instead of a scalar pair) of scale and offset parameters. An implementation of this method is compared with lossless compression, storing data at fixed relative precision (bit-grooming) and scalar linear packing in terms of compression ratio, accuracy and speed. When viewed as a trade-off between compression and error, layer-packing yields similar results to bit-grooming (storing between 3 and 4 significant figures). Bit-grooming and layer-packing offer significantly better control of precision than scalar linear packing. Relative performance, in terms of compression and errors, of bit-groomed and layer-packed data were strongly predicted by the entropy of the exponent array, and lossless compression was well predicted by entropy of the original data array. Layer-packed data files must be "unpacked" to be readily usable. The compression and precision characteristics make layer-packing a competitive archive format for many scientific data sets.

58 GEOSCIENCES↗

Secure Information Exchange Gateway for Electric Grid Operations

The major objectives of the SIEGate project were to improve the security posture and minimize the cyber-attack surface of electric utility control centers and to reduce the cost of maintaining control-room-to-control-room information exchange. Major project goals included the design, development, testing, and commercialization of a single security-hardened appliance that could meet industry needs for resisting cyber-attacks while protecting the confidentiality and integrity of a growing volume of real-time information needed to ensure the reliability of the bulk electric system and interoperating with existing data formats and networking technologies. The SIEGate project has achieved its goals and objectives. The SIEGate Design Document, issued in March 2012, presented SIEGate use cases, provided SIEGate requirements, established SIEGate design principles, and prescribed design functionality of SIEGate as well as the components that make up SIEGate. SIEGate Release Version 1.0 was posted in January 2014. Release Version 1.0.83, which was posted on March 28, 2014, fixed many issues discovered by early adopters and added several new features. Release Candidate 1.1, which added additional improvements and bug fixes, was posted in June 2014. SIEGate executables have been downloaded more than 300 times. SIEGate has been tested at PJM, Entergy, TVA, and Southern. Security testing and analysis of SIEGate has been conducted at PNNL and PJM. Alstom has provided a summary of recommended steps for commercialization of the SIEGate Appliance and identified two deployment models with immediate commercial application.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Accurate derivative evaluation for any Grad–Shafranov solver

We present a numerical scheme that can be combined with any fixed boundary finite element based Poisson or Grad–Shafranov solver to compute the first and second partial derivatives of the solution to these equations with the same order of convergence as the solution itself. At the heart of our scheme is an efficient and accurate computation of the Dirichlet to Neumann map through the evaluation of a singular volume integral and the solution to a Fredholm integral equation of the second kind. Our numerical method is particularly useful for magnetic confinement fusion simulations, since it allows the evaluation of quantities such as the magnetic field, the parallel current density and the magnetic curvature with much higher accuracy than has been previously feasible on the affordable coarse grids that are usually implemented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A dedicated cone-beam CT system for musculoskeletal extremities imaging: Design, optimization, and initial performance characterization

Purpose: This paper reports on the design and initial imaging performance of a dedicated cone-beam CT (CBCT) system for musculoskeletal (MSK) extremities. The system complements conventional CT and MR and offers a variety of potential clinical and logistical advantages that are likely to be of benefit to diagnosis, treatment planning, and assessment of therapy response in MSK radiology, orthopaedic surgery, and rheumatology. Methods: The scanner design incorporated a host of clinical requirements (e.g., ability to scan the weight-bearing knee in a natural stance) and was guided by theoretical and experimental analysis of image quality and dose. Such criteria identified the following basic scanner components and system configuration: a flat-panel detector (FPD, Varian 3030+, 0.194 mm pixels); and a low-power, fixed anode x-ray source with 0.5 mm focal spot (SourceRay XRS-125-7K-P, 0.875 kW) mounted on a retractable C-arm allowing for two scanning orientations with the capability for side entry, viz. a standing configuration for imaging of weight-bearing lower extremities and a sitting configuration for imaging of tensioned upper extremity and unloaded lower extremity. Theoretical modeling employed cascaded systems analysis of modulation transfer function (MTF) and detective quantum efficiency (DQE) computed as a function of system geometry, kVp and filtration, dose, source power, etc. Physical experimentation utilized an imaging bench simulating the scanner geometry for verification of theoretical results and investigation of other factors, such as antiscatter grid selection and 3D image quality in phantom and cadaver, including qualitative comparison to conventional CT. Results: Theoretical modeling and benchtop experimentation confirmed the basic suitability of the FPD and x-ray source mentioned above. Clinical requirements combined with analysis of MTF and DQE yielded the following system geometry: a {approx}55 cm source-to-detector distance; 1.3 magnification; a 20 cm diameter bore (20 x 20 x 20 cm{sup 3} field of view); total acquisition arc of {approx}240 deg. The system MTF declines to 50% at {approx}1.3 mm{sup -1} and to 10% at {approx}2.7 mm{sup -1}, consistent with sub-millimeter spatial resolution. Analysis of DQE suggested a nominal technique of 90 kVp (+0.3 mm Cu added filtration) to provide high imaging performance from {approx}500 projections at less than {approx}0.5 kW power, implying {approx}6.4 mGy (0.064 mSv) for low-dose protocols and {approx}15 mGy (0.15 mSv) for high-quality protocols. The experimental studies show improved image uniformity and contrast-to-noise ratio (without increase in dose) through incorporation of a custom 10:1 GR antiscatter grid. Cadaver images demonstrate exquisite bone detail, visualization of articular morphology, and soft-tissue visibility comparable to diagnostic CT (10-20 HU contrast resolution). Conclusions: The results indicate that the proposed system will deliver volumetric images of the extremities with soft-tissue contrast resolution comparable to diagnostic CT and improved spatial resolution at potentially reduced dose. Cascaded systems analysis provided a useful basis for system design and optimization without costly repeated experimentation. A combined process of design specification, image quality analysis, clinical feedback, and revision yielded a prototype that is now awaiting clinical pilot studies. Potential advantages of the proposed system include reduced space and cost, imaging of load-bearing extremities, and combined volumetric imaging with real-time fluoroscopy and digital radiography.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

The Dual Nature of Entrainment-Mixing Signatures Revealed through Large-Eddy Simulations of a Convection-Cloud Chamber

Abstract Entrainment of subsaturated air into a cloud can influence its optical and microphysical properties in various ways, depending on the droplet evaporation and turbulent mixing time scales. Previous experiments in the Pi convection-cloud chamber have revealed that, given a fixed entrained air property, the mixing of entrained subsaturated air results in complete evaporation of some cloud droplets, with the rest remaining unchanged. This is a signature of inhomogeneous mixing. While comparing the results of entrainment with varying air properties, the mixing signature appears as if the subsaturated air is well mixed with the cloud to evenly reduce the droplets’ size. In other words, taken together, the experiments appear to have the signature of homogeneous mixing. To explore these results in a greater depth, we conduct large-eddy simulations combined with a bin microphysics scheme. Our results reproduce the similar signatures of inhomogeneous and homogeneous mixing, implying that LES can resolve the inhomogeneous mixing when the grid spacing is smaller than the entrained air parcel. Additionally, we observe that increasing the aerosol injection rate enhances the signature of inhomogeneous mixing, while coarser grid spacing diminishes it. Finally, the change in wall fluxes in response to various entrained air properties confirms that the homogeneous signature seen in the analysis of an ensemble of simulations is the result of various equilibrium states. This further strengthens the suggestion that the homogeneous mixing signature found in aircraft observations near the cloud top may result from combining entrainment events of different intensities, possibly caused by various-sized eddies. Significance Statement Large-eddy simulation and size-resolved microphysics can resolve time scales for turbulent mixing and evaporation and, therefore, are well suited for reproducing, extending, and interpreting the entrainment experiment in the Pi convection-cloud chamber. Our simulation results confirm (i) the inhomogeneous mixing signature for an individual entrainment event and (ii) the appearance of homogeneous mixing in an ensemble of entrainment episodes. Furthermore, we demonstrate that the inhomogeneous mixing signature is more pronounced in a polluted cloud, but coarser grid spacing in simulations may compromise the accuracy of this signature. Last, the homogeneous mixing signature results from various equilibrium states established for different entrainment intensities and adjusted wall fluxes, which are challenging to measure experimentally but can be easily analyzed in the simulations.

54 ENVIRONMENTAL SCIENCES↗

Synthetic, Realistic Transmission and Distribution Co-Simulation for Voltage Control Benchmarking

As distributed energy resources (DERs) play an increasing role in generation and grid services, control and market paradigms will need to be tested with models that incorporate transmission and distribution systems. The test systems will need to be realistic and accessible for researchers to form conclusions about realistic behaviors and maintain comparability and repeatability of results. These test systems also must have high resolution models of sufficient size for realistic distribution and transmission interaction. This paper demonstrates the tightly coupled HELICS-based co-simulation of realistic synthetic distribution and transmission models for the Austin metropolitan area and how this framework can be used to analyze resultant behaviors from distributed resources and their controls. These results are compared to the same distribution system run without co-simulation using a fixed feed-in voltage. Results show that even without advanced voltage controls, the feed-in voltage response when coupled to the transmission system results a lower level of voltage excursions, suggesting the importance of using this type of realistic, tightly coupled, transmission and distribution cosimulation when assessing DERs.

27 ARPA - Advanced Research Projects Agency-Energy↗

Improving the Time Resolution of Hyper-Spectral Infrared Sounding Observations With Trajectory Enhancement

Low Earth orbit (LEO) hyper-spectral infrared (IR) sounders have significant potential for characterizing the complex evolution of thermodynamic environments favorable for convective initiation and ongoing convection. The snapshots at fixed local times are unable to provide the temporal resolution needed to resolve the rapidly evolving convective environment. A novel methodology using trajectory modeling coupled with satellite soundings was developed to create proximity soundings near NCEI Storm Events to investigate differences in severe weather environments (Kalmus et al., 2019, https://doi.org/10.1175/MWR-D-18-0055.1). This methodology was recently extended to entire satellite swaths forward in time up to six hours into the future and was evaluated during NOAA’s Hazardous Weather Testbed (HWT) in a quasi-operational setting (Kahn et al., 2023, https://doi.org/10.1175/WAF-D-22-0204.1). This methodology is based on parcel forward-trajectory calculations from the satellite observing time to recreate future soundings of temperature and moisture at regularly gridded intervals in space and time. Using coincident Multi-Radar Multi-Sensor (MRMS) rainfall estimates, we show that convective available potential energy (CAPE) is increased and convective inhibition (CIN) is decreased at times and locations where convection initiated. This methodology was evaluated with ERA5 data sampled to mimic the satellite observing swaths (Richardson et al., 2023, https://doi.org/10.5194/egusphere-2023-97). Approximately 60–90% of the temporal and spatial variability in temperature and humidity is explained by parcel advection using trajectory modeling. This method will be applied to the multi-decadal LEO satellite record and could be used to examine capabilities of future satellite missions such as the GeoXO mission.

Emily Berndt↗

Synthetic, Realistic Transmission and Distribution Co-Simulation for Voltage Control Benchmarking

As distributed energy resources (DERs) play an increasing role in generation and grid services, control and market paradigms will need to be tested with models that incorporate transmission and distribution systems. The test systems will need to be realistic and accessible for researchers to form conclusions about realistic behaviors and maintain comparability and repeatability of results. These test systems also must have high resolution models of sufficient size for realistic distribution and transmission interaction. This paper demonstrates the tightly coupled HELICS-based co-simulation of realistic synthetic distribution and transmission models for the Austin metropolitan area and how this framework can be used to analyze resultant behaviors from distributed resources and their controls. These results are compared to the same distribution system run without co-simulation using a fixed feed-in voltage. Results show that even without advanced voltage controls, the feed-in voltage response when coupled to the transmission system results a lower level of voltage excursions, suggesting the importance of using this type of realistic, tightly coupled, transmission and distribution co-simulation when assessing DERs.

27 ARPA - Advanced Research Projects Agency-Energy↗

An approach for drag correction based on the local heterogeneity for gas-solid flows

The drag models typically used for gas-solids interaction are mainly developed based on homogeneous systems of flow passing fixed particle assembly. It has been shown that the heterogeneous structures, i.e., clusters and bubbles in fluidized beds, need to be resolved to account for their effect in the numerical simulations. Since the heterogeneity is essentially captured through the local concentration gradient in the computational cells, this study proposes a simple approach to account for the non-uniformity of solids spatial distribution inside a computational cell and its effect on the interaction between gas and solid phases. Finally, to validate this approach, the predicted drag coefficient has been compared to the results from direct numerical simulations. In addition, the need to account for this type of heterogeneity is discussed for a periodic riser flow simulation with highly resolved numerical grids and the impact of the proposed correction for drag is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of multigrid methods on massively parallel computers: Architectural implications

We study the potential performance of multigrid algorithms running on massively parallel computers with the intent of discovering whether presently envisioned machines will provide an efficient platform for such algorithms. We consider the domain parallel version of the standard V cycle algorithm on model problems, discretized using finite difference techniques in two and three dimensions on block structured grids of size 10(exp 6) and 10(exp 9), respectively. Our models of parallel computation were developed to reflect the computing characteristics of the current generation of massively parallel multicomputers. These models are based on an interconnection network of 256 to 16,384 message passing, 'workstation size' processors executing in an SPMD mode. The first model accomplishes interprocessor communications through a multistage permutation network. The communication cost is a logarithmic function which is similar to the costs in a variety of different topologies. The second model allows single stage communication costs only. Both models were designed with information provided by machine developers and utilize implementation derived parameters. With the medium grain parallelism of the current generation and the high fixed cost of an interprocessor communication, our analysis suggests an efficient implementation requires the machine to support the efficient transmission of long messages, (up to 1000 words) or the high initiation cost of a communication must be significantly reduced through an alternative optimization technique. Furthermore, with variable length message capability, our analysis suggests the low diameter multistage networks provide little or no advantage over a simple single stage communications network.

Matheson, Lesley R.↗