Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial load balancing”

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 19 records

Multitasking a three-dimensional Navier-Stokes algorithm on the Cray-2

A three-dimensional computational aerodynamics algorithm has been multitasked for efficient parallel execution on the Cray-2. It provides a means for examining the multitasking performance of a complete CFD application code. An embedded zonal multigrid scheme is used to solve the Reynolds-averaged Navier-Stokes equations for an internal flow model problem. The explicit nature of each component of the method allows a spatial partitioning of the computational domain to achieve a well-balanced task load for MIMD computers with vector-processing capability. Experiments have been conducted with both two- and three-dimensional multitasked cases. The best speedup attained by an individual task group was 3.54 on four processors of the Cray-2, while the entire solver yielded a speedup of 2.67 on four processors for the three-dimensional case. The multiprocessing efficiency of various types of computational tasks is examined, performance on two Cray-2s with different memory access speeds is compared, and extrapolation to larger problems is discussed.

Swisshelm, Julie M.↗

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗

Particle simulation in a multiprocessor environment

A parallel implementation of a particle simulation method that is portable between a wide class of multiprocessor computers is presented. A fine grain spatial decomposition is utilized where several subdomains having a regular structure are computed at each processing node. This leads directly to an efficient and straightforward load balancing scheme if the number of subdomains at each processor is permitted to vary in an appropriate manner. Three dimensional simulations incorporating full thermochemical nonequilibrium are possible using the resulting code. Vectorizable algorithms are retained from earlier work allowing efficient use of deeply pipelined node processors where available. Performance results are presented from three different machine architectures demonstrating the portability of the code. On a 128-node Intel iPSC/860, performance is twice that of a single Cray-Y/MP CPU running a highly vectorized simulation code. Speedup is linear over the full range of number of processors on all target machines, indicating scalability of the method to higher degrees of parallelism.

Mcdonald, Jeffrey D.↗

Investigation of Unsteady Pressure-Sensitive Paint (uPSP) and a Dynamic Loads Balance to Predict Launch Vehicle Buffet Environments

This NESC assessment examined the accuracy of estimating buffet loads on in-line launch vehicles without booster attachments using sparse unsteady pressure measurements. The buffet loads computed using sparse sensor data were compared with estimates derived using measurements with much higher spatial resolution. The current method for estimating launch vehicle buffet loads is through wind tunnel testing of models with approximately 400 unsteady pressure transducers. Even with this relatively large number of sensors, the coverage can be insufficient to provide reliable integrated unsteady loads on vehicles. In general, sparse sensor spacing requires the use of coherence-length-based corrections in the azimuthal and axial directions to integrate the unsteady pressures and obtain reasonable estimates of the buffet loads. Coherence corrections have been used to estimate buffet loads for a variety of launch vehicles with the assumption methodology results in reasonably conservative loads. For the Space Launch System (SLS), the first estimates of buffet loads exceeded the limits of the vehicle structure, so additional tests with higher sensor density were conducted to better define the buffet loads and possibly avoid expensive modifications to the vehicle design. Without the additional tests and improvements to the coherence-length analysis methods, there would have been significant impacts to the vehicle weight, cost, and schedule. If the load estimates turn out to be too low, there is significant risk of structural failure of the vehicle. This assessment used a combination of unsteady pressure-sensitive paint (uPSP), unsteady pressure transducers, and a dynamic force and moment balance to investigate the integration schemes used with limited unsteady pressure data by comparing them with direct integration of extremely dense fluctuating pressure measurements. An outfall of the assessment was to evaluate the potential of using the emerging uPSP technique in a production test environment for future launch vehicles. The results show that modifications to the current technique can improve the accuracy of buffet estimates. More importantly, the uPSP worked remarkably well and, with improvements to the frequency response, sensitivity, and productivity, will provide an enhanced method for measuring wind tunnel buffet forcing functions (BFFs).

Schuster, David M.↗

Sensitivity of acoustic predictions to variation of input parameters

The noise prediction code WOPWOP predicts the thickness and loading noise produced by a helicopter rotor, given the blade motion, rotor operating conditions, and fluctuating force distribution over the blade surface. However, the criticality of these various inputs, and their respective effects on the predicted acoustic field, have never been fully addressed. This paper examines the importance of these inputs, and the sensitivity of the acoustic predictions to a variation of each parameter. The effects of collective and cyclic pitch, as well as coning and flapping, are presented. Blade loading inputs are examined to determine the necessary spatial and temporal resolution, as well as the importance of the cordwise distribution. The acoustic predictions show regions in the acoustic field where significant errors occur when simplified blade motions or blade loadings are used. An assessment of the variation in the predicted acoustic field is balanced by a consideration of CPU time necessary for the various approximations.

Brentner, Kenneth S.↗

Sensitivity of acoustic predictions to variation of input parameters

Rotor noise prediction codes predict the thickness and loading noise produced by a helicopter rotor, given the blade motion, rotor operating conditions, and fluctuating force distribution over the blade surface. However, the criticality of these various inputs, and their respective effects on the predicted acoustic field, have never been fully addressed. This paper examines the importance of these inputs, and the sensitivity of the acoustic predicitions to a variation of each parameter. The effects of collective and cyclic pitch, as well as coning and cyclic flapping, are presented. Blade loading inputs are examined to determine the necessary spatial and temporal resolution, as well as the importance of the chordwise distribution. The acoustic predictions show regions in the acoustic field where significant errors occur when simplified blade motions or blade loadings are used. An assessment of the variation in the predicted acoustic field is balanced by a consideration of Central Processing Unit (CPU) time necessary for the various approximations.

Brentner, Kenneth S.↗

Parallel Adaptive High-Order CFD Simulations Characterizing Cavity Acoustics for the Complete SOFIA Aircraft

This paper presents one-of-a-kind MPI-parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting of the telescope. A temporally fourth-order Runge-Kutta, and spatially fifth-order WENO-5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32,000 cores and 4 billion cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregularities caused by the highly complex geometry. Limits to scaling beyond 32K cores are identified, and targeted code optimizations are discussed.

Acoustics↗

Parallel Adjective High-Order CFD Simulations Characterizing SOFIA Cavity Acoustics

This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A temporally fourth-order accurate Runge-Kutta, and spatially fth-order accurate WENO- 5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.

CFD↗

Parallel Adaptive High-Order CFD Simulations Characterizing SOFIA Cavitiy Acoustics

This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A tempo- rally fourth-order accurate Runge-Kutta, and a spatially fth-order accurate WENO-5Z scheme were used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.

Parallell↗

Research and test facilities

A description is given of each of the following Langley research and test facilities: 0.3-Meter Transonic Cryogenic Tunnel, 7-by 10-Foot High Speed Tunnel, 8-Foot Transonic Pressure Tunnel, 13-Inch Magnetic Suspension & Balance System, 14-by 22-Foot Subsonic Tunnel, 16-Foot Transonic Tunnel, 16-by 24-Inch Water Tunnel, 20-Foot Vertical Spin Tunnel, 30-by 60-Foot Wind Tunnel, Advanced Civil Transport Simulator (ACTS), Advanced Technology Research Laboratory, Aerospace Controls Research Laboratory (ACRL), Aerothermal Loads Complex, Aircraft Landing Dynamics Facility (ALDF), Avionics Integration Research Laboratory, Basic Aerodynamics Research Tunnel (BART), Compact Range Test Facility, Differential Maneuvering Simulator (DMS), Enhanced/Synthetic Vision & Spatial Displays Laboratory, Experimental Test Range (ETR) Flight Research Facility, General Aviation Simulator (GAS), High Intensity Radiated Fields Facility, Human Engineering Methods Laboratory, Hypersonic Facilities Complex, Impact Dynamics Research Facility, Jet Noise Laboratory & Anechoic Jet Facility, Light Alloy Laboratory, Low Frequency Antenna Test Facility, Low Turbulence Pressure Tunnel, Mechanics of Metals Laboratory, National Transonic Facility (NTF), NDE Research Laboratory, Polymers & Composites Laboratory, Pyrotechnic Test Facility, Quiet Flow Facility, Robotics Facilities, Scientific Visualization System, Scramjet Test Complex, Space Materials Research Laboratory, Space Simulation & Environmental Test Complex, Structural Dynamics Research Laboratory, Structural Dynamics Test Beds, Structures & Materials Research Laboratory, Supersonic Low Disturbance Pilot Tunnel, Thermal Acoustic Fatigue Apparatus (TAFA), Transonic Dynamics Tunnel (TDT), Transport Systems Research Vehicle, Unitary Plan Wind Tunnel, and the Visual Motion Simulator (VMS).

Source record↗

Vaporization characteristics of carbon heat shields under radiative heating.

Study of the vaporization characteristics of samples of ATJ graphite, a material that has been considered for use on a Jovian probe. These samples were subjected to radiative heating loads of approximately 2 kW/sq cm in argon atmospheres of pressures from 0.00046 to 1 atm. Surface temperatures, mass loss rates, and spatially resolved emission spectral data were recorded. These data are analyzed to determine carbon vapor pressure as a function of temperature and are compared with current models for the vapor pressure of carbon. The effects of finite vaporization (i.e., nonequilibrium) rates are considered and compared with experiment. Estimates of the heat of vaporization from an energy balance are also presented.

Davy, W. C.↗

Design and implementation of scalable tape archiver

In order to reduce costs, computer manufacturers try to use commodity parts as much as possible. Mainframes using proprietary processors are being replaced by high performance RISC microprocessor-based workstations, which are further being replaced by the commodity microprocessor used in personal computers. Highly reliable disks for mainframes are also being replaced by disk arrays, which are complexes of disk drives. In this paper we try to clarify the feasibility of a large scale tertiary storage system composed of 8-mm tape archivers utilizing robotics. In the near future, the 8-mm tape archiver will be widely used and become a commodity part, since recent rapid growth of multimedia applications requires much larger storage than disk drives can provide. We designed a scalable tape archiver which connects as many 8-mm tape archivers (element archivers) as possible. In the scalable archiver, robotics can exchange a cassette tape between two adjacent element archivers mechanically. Thus, we can build a large scalable archiver inexpensively. In addition, a sophisticated migration mechanism distributes frequently accessed tapes (hot tapes) evenly among all of the element archivers, which improves the throughput considerably. Even with the failures of some tape drives, the system dynamically redistributes hot tapes to the other element archivers which have live tape drives. Several kinds of specially tailored huge archivers are on the market, however, the 8-mm tape scalable archiver could replace them. To maintain high performance in spite of high access locality when a large number of archivers are attached to the scalable archiver, it is necessary to scatter frequently accessed cassettes among the element archivers and to use the tape drives efficiently. For this purpose, we introduce two cassette migration algorithms, foreground migration and background migration. Background migration transfers cassettes between element archivers to redistribute frequently accessed cassettes, thus balancing the load of each archiver. Background migration occurs the robotics are idle. Both migration algorithms are based on access frequency and space utility of each element archiver. To normalize these parameters according to the number of drives in each element archiver, it is possible to maintain high performance even if some tape drives fail. We found that the foreground migration is efficient at reducing access response time. Beside the foreground migration, the background migration makes it possible to track the transition of spatial access locality quickly.

Nemoto, Toshihiro↗

An implementation of a tree code on a SIMD, parallel computer

We describe a fast tree algorithm for gravitational N-body simulation on SIMD parallel computers. The tree construction uses fast, parallel sorts. The sorted lists are recursively divided along their x, y and z coordinates. This data structure is a completely balanced tree (i.e., each particle is paired with exactly one other particle) and maintains good spatial locality. An implementation of this tree-building algorithm on a 16k processor Maspar MP-1 performs well and constitutes only a small fraction (approximately 15%) of the entire cycle of finding the accelerations. Each node in the tree is treated as a monopole. The tree search and the summation of accelerations also perform well. During the tree search, node data that is needed from another processor is simply fetched. Roughly 55% of the tree search time is spent in communications between processors. We apply the code to two problems of astrophysical interest. The first is a simulation of the close passage of two gravitationally, interacting, disk galaxies using 65,636 particles. We also simulate the formation of structure in an expanding, model universe using 1,048,576 particles. Our code attains speeds comparable to one head of a Cray Y-MP, so single instruction, multiple data (SIMD) type computers can be used for these simulations. The cost/performance ratio for SIMD machines like the Maspar MP-1 make them an extremely attractive alternative to either vector processors or large multiple instruction, multiple data (MIMD) type parallel computers. With further optimizations (e.g., more careful load balancing), speeds in excess of today's vector processing computers should be possible.

Olson, Kevin M.↗

Climatology of Aerosol Optical Properties in Southern Africa

A thorough regionally dependent understanding of optical properties of aerosols and their spatial and temporal distribution is required before we can accurately evaluate aerosol effects in the climate system. Long term measurements of aerosol optical depth, Angstrom exponent and retrieved single scattering albedo and size distribution, were analyzed and compiled into an aerosol optical properties climatology for southern Africa. Monitoring of aerosol parameters have been made by the AERONET program since the middle of the last decade in southern Africa. This valuable information provided an opportunity for understanding how aerosols of different types influence the regional radiation budget. Two long term sites, Mongu in Zambia and Skukuza in South Africa formed the core sources of data in this study. Results show that seasonal variation of aerosol optical thicknesses at 500 nm in southern Africa are characterized by low seasonal multi-month mean values (0.11 to 0.17) from December to May, medium values (0.20 to 0.27) between June and August, and high to very high values (0.30 to 0.46) during September to November. The spatial distribution of aerosol loadings shows that the north has high magnitudes than the south in the biomass burning season and the opposite in none biomass burning season. From the present aerosol data, no long term discernable trends are observable in aerosol concentrations in this region. This study also reveals that biomass burning aerosols contribute the bulk of the aerosol loading in August-October. Therefore if biomass burning could be controlled, southern Africa will experience a significant reduction in total atmospheric aerosol loading. In addition to that, aerosol volume size distribution is characterized by low concentrations in the non biomass burning period and well balanced particle size contributions of both coarse and fine modes. In contrast high concentrations are characteristic of biomass burning period, combined with significant dominance of fine mode particles.

Queface, Antonio J.↗

Coupling Mars' Dust and Water Cycles: Effects on Dust Lifting Vigor, Spatial Extent and Seasonality

The dust cycle is an important component of Mars' current climate system. Airborne dust affects the radiative balance of the atmosphere, thus greatly influencing the thermal and dynamical state of the atmosphere. Dust raising events on Mars occur at spatial scales ranging from meters to planet-wide. Although the occurrence and season of large regional and global dust storms are highly variable from one year to the next, there are many features of the dust cycle that occur year after year. Generally, a low-level dust haze is maintained during northern spring and summer, while elevated levels of atmospheric dust occur during northern autumn and winter. During years without global-scale dust storms, two peaks in total dust loading were observed by MGS/TES: one peak occurred before northern winter solstice at Ls 200-240, and one peak occurred after northern winter solstice at L(sub s) 305-340. These maxima in dust loading are thought to be associated with transient eddy activity in the northern hemisphere, which has been observed to maximize pre- and post-solstice. Interactive dust cycle studies with Mars General Circulation Models (MGCMs) have included the lifting, transport, and sedimentation of radiatively active dust. Although the predicted global dust loadings from these simulations capture some aspects of the observed dust cycle, there are marked differences between the simulated and observed dust cycles. Most notably, the maximum dust loading is robustly predicted by models to occur near northern winter solstice and is due to dust lifting associated with down slope flows on the flanks of the Hellas basin. Thus far, models have had difficulty simulating the observed pre- and post- solstice peaks in dust loading. Interactive dust cycle studies typically have not included the formation of water ice clouds or their radiative effects. Water ice clouds can influence the dust cycle by scavenging dust from atmosphere and by interacting with solar and infrared radiation, thereby modifying the thermal structure of the atmosphere and its circulation. Results presented in other papers at this workshop show that including the radiative effects of water ice clouds greatly influence the water cycle and the vigor of weather systems in both the northern and southern hemispheres. Our goal is to investigate the effects of fully coupling the dust and water cycles on the dust cycle. We show that including water ice clouds and their radiative effects greatly affect the magnitude, spatial extent and seasonality of dust lifting and the season of maximum atmospheric dust loading.

Kahre, M. A.↗

Sensitivity analysis for the dynamic response of thermoviscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric, fully-coupled, thermoviscoplastic response of shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of large-strain, transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be viscoplasticity with strain hardening and thermal hardening, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used for the shell equations with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. The energy-balance equation is solved using a Galerkin procedure, with the temperature as the fundamental unknown. Spatial discretization is performed in one dimension (meridional direction) for the momentum and constitutive equations of the shell, and in two dimensions (meridional and thickness directions) for the energy-balance equation. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) for the momentum equation; a predictor-corrector version of the trapezoidal rule is used for the energy-balance equation; and an explicit scheme consistent with the central difference method is used to integrate the constitutive equations. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading.

Kulkarni, Makarand↗

Measurement and evaluation of swirl-type flow in labyrinth seals of conventional turbine stages

The effects of load factor and rotor eccentricity were determined on flow conditions in test series for two stages. The results indicate that swirl-type entry flow follows the rules of potential swirl. Within the labyrinth cavities two spatial separated flow areas are considered. A dominating flow in periphera direction nearly fills the space between the sealing strips and the ceiling of the cavity. Below this flow, an area of axial mass transport is situated, with a slight peripheral component, limited on the nearest surroundings of the seals gap and the rotor surface. Between both flows, an exchange of energy takes place. Within the gaps, flow direction depends on axial velocity and therefore on variable flow contraction. A balance of energy within the seal and the cavities interprets the results of lateral force measurements as an influence of friction at the sealing strips surface and the rotating shaft surface. Stages with their blades put together in buckets by means of shrouding segments are particularly influenced by the rotating speed of the shrouding.

Hauck, L.↗

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗