Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid computing”

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 937 records · Page 52

ADI on staggered mesh - A method for the calculation of compressible convection

An alternating direction implicit (ADI) method has been applied to a staggered grid for the computation of convection in a highly stratified fluid. Since artificial viscosity is not needed, subtle effects like the onset of convection can be studied. These computations compare well with the 2-D results by Graham (1975) and also agree with standard Boussinesq results when taken to that limit. Good efficiency has been achieved with a time step hundreds of times larger than the stability limit imposed by the explicit treatment of diffusion and the Courant number is not restricted to be below 1. The Navier-Stokes equation contains cross spatial derivatives which are treated explicitly in most ADI schemes. The destabilizing effect of such a practice on a 2-D model system with second-order spatial derivative terms only was analyzed and found to be not excessive.

Chan, K. L.↗

Boundary conditions for the solution of compressible Navier-Stokes equations by an implicit factored method

A method is presented for formulating the boundary conditions in implicit finite-difference form needed for obtaining solutions to the compressible Navier-Stokes equations by the Beam and Warming implicit factored method. The usefulness of the method was demonstrated (a) by establishing the boundary conditions applicable to the analysis of the flow inside an axisymmetric piston-cylinder configuration and (b) by calculating velocities and mass fractions inside the cylinder for different geometries and different operating conditions. Stability, selection of time step and grid sizes, and computer time requirements are discussed in reference to the piston-cylinder problem analyzed.

Shih, T. I.-P.↗

Use of a hyperbolic grid generation scheme in simulating supersonic viscous flow about three-dimensional winged configuration

The present paper describes a numerical mesh generation technique to be used with an implicit finite difference method for simulating visous supersonic flow about low-aspect-ratio wing body configurations using a single grid strategy. The computational domain is segmented into multiple regions, with borders located in supersonic areas to avoid the otherwise costly interfacing procedure between adjacent segments. The numerical procedure is applied to calculate the turbulent flow around the shuttle orbiter and a canard projectile at supersonic free stream Mach number.

Rizk, Y. M.↗

Numerical study of three-dimensional separation and flow control at a wing/body junction

The problem of three-dimensional separation and flow control at a wing/body junction has been investigated numerically using a three-dimensional Navier-Stokes code. The numerical code employs an algebraic grid generation technique for generating the grid for unmodified junction and an elliptic grid generation technique for filleted fin junction. The results for laminar flow past a blunt fin/flat plate junction demonstrate that after grid refinement, the computations agree with experiment and reveal a strong dependency of the number of vortices at the junction on Mach number and Reynolds number. The numerical results for pressure distribution, particle paths and limiting streamlines for turbulent flow past a swept fin show a decrease in the peak pressure and in the extent of the separated flow region compared to the laminar case. The results for a filleted juncture indicate that the streamline patterns lose much of their vortical character with proper filleting. Fillets with a radius of three and one-half times the fin leading edge diameter or two times the incoming boundary layer thickness, significantly weaken the usual necklace interaction vortex for the Mach number and Reynolds number considered in the present study.

Ash, Robert L.↗

Aerodynamic applications of an efficient incompressible Navier-Stokes solver

A robust and accurate upwind scheme based on flux-difference splitting for three-dimensional, incompressible flows in general coordinates is described. Viscous effects are included through a thin-layer approximation to the Navier-Stokes equations. Turbulent flows are simulated with extensions to the Baldwin-Lomax (1978) turbulence model. The effects of accuracy and gridding on the computed flow field results are assessed. The results obtained from the code for vortical flows over a delta wing, a double delta wing, and a tangent-ogive forebody compare well with experimental data.

Hartwich, Peter M.↗

Parallel algorithms for interactive manipulation of digital terrain models

Interactive three-dimensional graphics applications, such as terrain data representation and manipulation, require extensive arithmetic processing. Massively parallel machines are attractive for this application since they offer high computational rates, and grid connected architectures provide a natural mapping for grid based terrain models. Presented here are algorithms for data movement on the massive parallel processor (MPP) in support of pan and zoom functions over large data grids. It is an extension of earlier work that demonstrated real-time performance of graphics functions on grids that were equal in size to the physical dimensions of the MPP. When the dimensions of a data grid exceed the processing array size, data is packed in the array memory. Windows of the total data grid are interactively selected for processing. Movement of packed data is needed to distribute items across the array for efficient parallel processing. Execution time for data movement was found to exceed that for arithmetic aspects of graphics functions. Performance figures are given for routines written in MPP Pascal.

Davis, E. W.↗

Hyperbolic conservation laws and numerical methods

The mathematical structure of hyperbolic systems and the scalar equation case of conservation laws are discussed. Linear, nonlinear systems and the Riemann problem for the Euler equations are also studied. The numerical methods for conservation laws are presented in a nonstandard manner which leads to large time steps generalizations and computations on irregular grids. The solution of conservation laws with stiff source terms is examined.

Leveque, Randall J.↗

Charge transfer and X-ray emission from supernova remnants

The X-ray line emission excited by charge transfer between neutral hydrogen and the hydrogenic or fully stripped ions of heavy elements has been computed for a grid of nonequilibrium ionization models for supernova remnants. The effects of charge transfer excitation on the X-ray spectra of supernova remnants are discussed. The charge transfer process depends very sensitively upon the abundance of neutral hydrogen behind the shock. Consequently, a careful treatment of the ionization of hydrogen was required. The model line fluxes are compared with those excited by conventional collisional processes involving electrons. Results are presented in the form of contour plots for both the flux excited by charge transfer and the ratio of the charge transfer flux to that associated solely with collisional processes involving electrons. Charge transfer never contributes more than 10 percent of the flux of the strongest electron collisionally excited lines, and its contribution is typically 0.001 to 10 to the -5th for a shock moving into a medium composed of neutral hydrogen.

Wise, Michael W.↗

Unsteady inviscid and viscous computations for vortex-dominated flows

The unsteady supersonic flow around a rigid sharp-edged delta wing is solved using the unsteady Euler and thin-layer Navier-Stokes equations. The problem is formulated relative to a moving frame of reference so that computation of the grid motion is not required. A three-dimensional, implicit, aproximately factored, central-differencing, finite-volume scheme is used to obtain the time-accurate, unsteady, locally conical flow. The inviscid and viscous results are compared during the periodic response of a wing undergoing forced rolling-oscillation motion. The results show substantial differences in the surface pressure results and other distributed flow characteristics. These differences are attributed to the inviscid and viscous flow predictions in the vortex-dominated regions. The results of the total loads do not show any significant differences.

Kandil, Osama A.↗

Three-dimensional simulation of slender delta wing rock and divergence

Computational simulation of three-dimensional flows around a delta wing undergoing rock and roll-divergence motions is presented. The problem is a multidisciplinary one where fluid-dynamics equations and rigid-body-dynamics equations are sequentially solved. For the fluid-dynamics part, the unsteady Euler equations, which are written relative to a moving frame of reference, are solved using an implicit, approximately-factored, central-difference, finite-volume scheme. For the rigid-body dynamics part, the Euler equation of rigid-body rolling motion is solved using a four-stage Runge-Kutta scheme. Since the applications do not include deforming wings or relative-rigid-body motions, the computational-fluid-dynamics grid, which is fixed in the moving frame of reference, does not need to be updated once it is generated.

Kandil, Osama A.↗

A global mean sea surface based upon GEOS 3 and Seasat altimeter data

A mean sea surface relative to the International Union of Geodesy 1980 Geodetic Reference System reference ellipsoid has been derived from Seasat and GEOS 3 altimeter measurements. This surface, called MSS-9012, has been computed on a grid of 1/8 deg resolution. The surface covers the global ocean between 70 deg N and 72 deg S. Comparisons are made between MSS-9012 and ocean bathymetry. The eastern extent of the Chain Fracture Zone in the Gulf of Guinea is more apparent in the altimetry than in the bathymetry data, as are many other features. The combination of data from the two satellites has successfully retrieved more information about the sea surface than was previously possible using only Seasat data.

Marsh, J. G.↗

Gigaflop (billion floating point operations per second) performance for computational electromagnetics

Accurate and rapid evaluation of radar signature for alternative aircraft/store configurations would be of substantial benefit in the evolution of integrated designs that meet radar cross-section (RCS) requirements across the threat spectrum. Finite-volume time domain methods offer the possibility of modeling the whole aircraft, including penetrable regions and stores, at longer wavelengths on today's gigaflop supercomputers and at typical airborne radar wavelengths on the teraflop computers of tomorrow. A structured-grid finite-volume time domain computational fluid dynamics (CFD)-based RCS code has been developed at the Rockwell Science Center, and this code incorporates modeling techniques for general radar absorbing materials and structures. Using this work as a base, the goal of the CFD-based CEM effort is to define, implement and evaluate various code development issues suitable for rapid prototype signature prediction.

Shankar, V.↗

Hybrid active vibration control of rotorbearing systems using piezoelectric actuators

The vibrations of a flexible rotor are controlled using piezoelectric actuators. The controller includes active analog components and a hybrid interface with a digital computer. The computer utilizes a grid search algorithm to select feedback gains that minimize a vibration norm at a specific operating speed. These gains are then downloaded as active stillnesses and dampings with a linear fit throughout the operating speed range to obtain a very effective vibration control.

Palazzolo, A. B.↗

Domain connectivity among systems of overset grids

This progress report describes research performed to achieve algorithmic improvements in the use of overset grids for the computation of unsteady flowfields about geometrically complex and moving component configurations. A hole-map computational method is detailed that reduces computational demand and corresponding demand on human resources. Results indicate that the method is both efficient and flexible for use in complex simulation cases. Tests were performed on shuttle orbiter/external tank, wing/store, and V-22 tilt-rotor configurations. The number of simulated points and the corresponding time on a SGI 4D-210 workstation are presented.

Meakin, Robert L.↗

Semi-span model testing in the National Transonic Facility

A semi-span testing technique has been proposed for the NASA Langley Research Center's National Transonic Facility (NTF). Semi-span testing has several advantages including (1) larger model size, giving increased Reynolds number capability; (2) improved model fidelity, allowing ease of flap and slat positioning which ultimately improves data quality; and (3) reduced construction costs compared with a full-span model. In addition, the increased model size inherently allows for increased model strength, reducing aeroelastic effects at the high dynamic pressure levels necessary to simulate flight Reynolds numbers. The Energy Efficient Transport (EET) full-span model has been modified to become the EET semi-span model. The full-span EET model was tested extensively at both NASA LRC and NASA Ames Research Center. The available full-span data will be useful in validating the semi-span test strategy in the NTF. In spite of the advantages discussed above, the use of a semi-span model does introduce additional challenges which must be addressed in the testing procedure. To minimize the influence of the sidewall boundary layer on the flow over the semi-span model, the model must be off-set from the sidewall. The objective is to remove the semi-span model from the sidewall boundary layer by use of a stand-off geometry. When this is done however, the symmetry along the centerline of the full-span model is lost when the semi-span model is mounted on the wind tunnel sidewall. In addition, the large semi-span model will impose a significant pressure loading on the sidewall boundary layer, which may cause separation. Even under flow conditions where the sidewall boundary layer remains attached, the sidewall boundary layer may adversely effect the flow over the semi-span model. Also, the increased model size and sidewall mounting requires a modified wall correction strategy. With these issues in mind, the semi-span model has been well instrumented with surface pressure taps to obtain data on the expected complex flow field in the near wall region. This status report summarizes the progress to date on developing the semi-span geometry definition suitable for generating structured grids for the computational research. In addition, the progress on evaluating three state-of-the-art Navier-Stokes codes is presented.

Chokani, Ndaona↗

Computational geometry issues

Presented in viewgraph format are some of the benefits of computational geometry for grid generation, including benefits to grid quality and algorithm efficiency. Many efficient algorithms are available, but it is concluded that high aspect ratio triangulations are not among them.

Mary-Anne Posenau↗

On the Accuracy of CFD-Based Pressure Drop Predictions for Right-Angle Ducts

The predictive capability of computational fluid dynamics (CFD) codes for turbulent flow through curved ducts is of significant importance to the design and performance analysis of modern rocket engine flowpaths. Code calibration and validation studies for this class of flow are desireable to estimate the performance margin and operating range of components designed using Navier-Stokes methods. Parametric experimental studies such as that of Weske (NACA ARR W-39) provided a wealth of performance data for the design of single- and compound elbow configurations with various cross-sections, curvature and aspect ratios at varying Reynolds numbers. In that work, the majority of data is presented in the form of loss coefficients, characterizing pressure losses due to duct curvature, and including losses due to wall friction. Using measured friction coefficients, losses of equivalent straight lengths of duct are subtracted, resulting in performance curves useful for design computations. These data are currently used in a CFD-based parametric study covering a broad range of operating conditions. Of particular interest for the accuracy of CFD predictions are the effects on pressure loss due to inlet boundary layer thickness (dependent on upstream development length), and the wall treatment for the turbulence equations (conventional wall functions vs. wall integration using a two-layer model). The experimental data are reassessed in the form of an error analysis, and are compared with CFD predictions for 18 computational cases. Grid-independence, grid spacing, and convergence requirements of the cases are discussed. Conclusions regarding the relative importance of the parametric variables will be presented.

Brankovic, Andreja↗

Grains in ionized nebulae: Spectral line diagnostics

The depletion of condensable elements onto grains in gaseous nebulae can provide evidence that dust is well mixed with the ionized gas. Al and CA are two of the most depleted elements in the general interstellar medium, and it is therefore important to measure their abundances within the ionized region of nebulae. We compute a large grid of photoionization models and identify sets of line ratios which are relatively insensitive to stellar and nebular parameters, and are thus excellent diagnostics for determining relative abundances. Based on the absence of the ((Ca II) lambda lambda 291, 7324 doublet and the detection of Al II) lambda lambda 2660, 2669 in the ultraviolet, we determine the extent of aluminum and calcium depletion onto grains in NGC 7027 and the Orion Nebula. Our results show a approximately 0.3 dex depletion for Al, but a depletion of more than two and a half orders of magnitude for Ca. A similar calculation based on Mg II lambda 2798 yields roughly a 0.8 dex depletion for Mg. This reaffirms the discrepancy between depletion determined from high and low ionization Mg lines. We also find evidence for a 'depletion gradient' in Ca in NGC 7027, since the calcium depletion we infer for the outer, more neutral regions using (Ca II) is somewhat higher than that inferred for the inner high-ionization region, using (Ca v). This gradient can test current models of the survival of grains within hot ionized gas.

Kingdon, J.↗