Search NASA⌕ Search

SEARCH · Search NASA

Results for “space 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 199 records · Page 11

Two-dimensional Navier-Stokes heat transfer analysis for rough turbine blades

A quasi-three-dimensional thin-layer Navier-Stokes analysis was used to predict heat transfer to rough surfaces. Comparisons are made between predicted and experimental heat transfer for turbine blades and flat plates of known roughness. The effect of surface roughness on heat transfer was modeled using a mixing length approach. The effect of near-wall grid spacing and convergence criteria on the accuracy of the heat transfer predictions are examined. An eddy viscosity mixing length model having an inner and outer layer was used. A discussion of the appropriate model for the crossover between the inner and outer layers is included. The analytic results are compared with experimental data for both flat plates and turbine blade geometries. Comparisons between predicted and experimental heat transfer showed that a modeling roughness effects using a modified mixing length approach results in good predictions of the trends in heat transfer due to roughness.

Boyle, R. J.↗

Magnetic effects in Venus/solar wind interaction

The overall objective of this research program is to better understand the interaction of a magnetized solar wind with the Venus atmosphere through the use of numerical solutions of the time-dependent, 2-D and 3-D magnetohydrodynamic (MHD) equations. Due to the more modest CPU requirements for the 2-D simulations, they were used for studies in which useful information not dependent on the third dimension could be obtained. The 2-D simulations served several purposes in addition to providing useful physical insight. They were used to determine the numerical parameters required in the 3-D studies, such as the grid spacing required to resolve particular features, and the damping that must be included to remove high-frequency oscillations. Among the specific studies performed with support from this grant that are discussed in this report are the following: comparison with other available models for purposes of testing the code and obtaining a baseline with which to evaluate the effects of additional physical processes, effects of a finite planet conductivity, bow shock standoff distance, and formation of the magnetic barrier and slippage of the magnetic field around the planet. A brief description of the methodology is presented before discussing the results.

Steinolfson, Richard S.↗

A typical pointed observation in IRAS designed to fit into a 15 minute window

A typical pointed observation in the Infrared Astronomical Satellite (IRAS) was designed to fit into a 15 minute window. Observing constraints, however, did not always permit targets from the scientific program to be observed on exactly 15 minute centers and a few months into the mission it was realized that the automatic scheduling program left times when no observations were being made. In order to use this potentially wasted observing time, a mode of pointed observations, the 'filler' mode, with shorter observation duration was initiated. Locations were picked on an arbitrary grid, spaced relatively uniformly in ecliptic coordinates, and the telescope was pointed to one of these sites whenever a gap potentially occurred. One of us (GXN) was the coordinator of this effort (the FL sub-group). There are a total of 22 FL fields, 13 of which have more than 50 component observations. These 13 fields cover about 30 square degrees. It was proposed to concentrate on one filler field, FL29, a field with 53 coverages which thus had among the deepest coverage of those observed in this mode. This field was also observed extensively at various non-infrared wavelengths by a variety of techniques. Because the IRAS observations of FL29 were made at the limits of the survey sensitivity and at different twist angles it was necessary to make significant technical advances in the methods used in the processing and analyzing the data. In particular, if normal IRAS coaddition procedures were used, the noise resulting from coaddition of the multiple frames did not decrease in the manner expected once significantly different twist angles were included in the coaddition. The technical effort needed to reduce the observations of FL29 was completed and now the noise behavior is understood and well behaved. The reduction of the field FL29 observations is done.

Neugebauer, G.↗

Diurnal tides in the Arctic Ocean

A 2D numerical model with a space grid of about 14 km is applied to calculate diurnal tidal constituents K(1) and O(1) in the Arctic Ocean. Calculated corange and cotidal charts show that along the continental slope, local regions of increased sea level amplitude, highly variable phase and enhanced currents occur. It is shown that in these local regions, shelf waves (topographic waves) of tidal origin are generated. In the Arctic Ocean and Northern Atlantic Ocean more than 30 regions of enhanced currents are identified. To prove the near-resonant interaction of the diurnal tides with the local bottom topography, the natural periods of oscillations for all regions have been calculated. The flux of energy averaged over the tidal period depicts the gyres of semitrapped energy, suggesting that the shelf waves are partially trapped over the irregularities of the bottom topography. It is shown that the occurrence of near-resonance phenomenon changes the energy flow in the tidal waves. First, the flux of energy from the astronomical sources is amplified in the shelf wave regions, and afterwards the tidal energy is strongly dissipated in the same regions.

Kowalik, Z.↗

A numerical analysis of tonal acoustics in rotor-stator interactions

In this study, the unsteady, thin-layer, Navier-Stokes equations are solved using a system of patched grids for a rotor-stator configuration of an axial turbine. The study focuses on the plurality of spinning modes that are present in such an interaction. The propagation of these modes in the upstream and downstream regions is analyzed and compared with numerical results. It was found that the numerically calculated tonal acoustics could be affected by the type of numerical boundary conditions employed at the inlet and exit of the computational boundaries and the grid spacing in the upstream and downstream regions. Results in the form of surface pressure amplitudes and the spectra of turbine tones and their far field behavior are presented. Numerical results and experimental data are compared wherever possible. The 'mode-content' for different harmonics of blade-passage frequency is shown to conform with that predicted by a kinematical analysis.

Rangwalla, A. A.↗

Instabilities orginating from suction holes used for Laminar Flow Control (LFC)

A small-scale wind tunnel previously used for turbulent boundary layer studies has been modified for experiments in laminar flow control. The facility incorporates suction through interchangeable porous test surfaces which are used to stabilize the boundary layer and delay transition to turbulent flow. The thin porous test surfaces are supported by a baffled plenum chamber box which also acts to gather the flow through the surface into tubes which are routed to a high pressure fan. An elliptic leading edge is attached to the assembly to establish a new layer on the test plate. A slot is used to remove the test section flow below the leading edge. The test section was lengthened and fitted with a new ceiling. Substantial modifications were also made to the 3D probe traverse. Detailed studies have been made using isolated holes to explore the underlying instability mechanisms. The suction is perturbed harmonically and data are averaged on the basis of the phase of the disturbance. Conditions corresponding to strong suction and without suction have been studied. In both cases, 3D contour surfaces in the vicinity of the hole show highly three-dimensional T-S waves that fan out away from the hole with streamwise distance. With suction, the perturbations on the centerline are much stronger and decay less rapidly, while the far field is similar to the case without suction. Downstream the contour surfaces of the bow-shaped TS waves develop spanwise irregularities which eventually form into clumps. The contours remain smooth when suction is not applied. Even without suction, the harmonic point source is challenging for CFD; e.g. DNS has been used for streamwise growth. With suction, grid resources are consumed by the hole and this makes DNS even more expensive. Limited DNS results so far indicate that the vortices which emanate from suction holes appear to be stable. The spanwise clumping observed in the experiment is evidence of a secondary instability that could be associated with suction vortices. A typical porous surface for LFC consists of 0.002 inch diameter holes with 0.020 inch grid spacing L, which is too small to resolve disturbances. A 20:1 scale porous test surface has been machined for improved spatial resolution while the L/d is still representative of flight conditions. Designers of porous surfaces use Goldsmith's criterion to minimize crossstream interaction. However nothing is known about the streamwise interactions. Results using two holes, aligned but displaced in the streamwise direction, indicate that partial TS wave cancellation is possible, depending on the hole spacing and disturbance frequency. Using DNS for streamwise interaction studies will be prohibitively expensive if linear superposition cannot be used for the multiple holes.

Watmuff, Jonathan H.↗

Aeroacoustic computation of gust-blade interaction

To better understand and address the challenges faced in computing the acoustics of flow fields, test problems must be considered. In the present study, the sound radiated by the interaction of a flat plate with an oncoming gust containing a two component, mean velocity is computed. The gust has a uniform mean flow in x with Mach number M(infinity) equal to 0.5. The gust's mean velocity in y is of smaller amplitude and is given by: v = 0.1 sin(pi/8(x/M(sub infinity) - t)). This problem has been posed for an upcoming ICASE/LaRC workshop on benchmark problems in computational aeroacoustics. A plate with a length of 30 units in x is used. The plate is assumed to be infinitesimally thin and is centered at the origin. All variables are made dimensionless using the scales specified. Acoustic quantities are obtained by numerically integrating the linearized Euler equations. Integration is performed on the computational domain -100.0 less than or equal to x less than or equal to 100.0, -100.0 less than or equal to y less than or equal to 100.0, using unit length grid spacing in x and in y. An integration scheme is sought which will provide accurate solution to the small quantities of interest at a minimal computational expense. Results indicate that with the given discretization a scheme of minimal fourth order accuracy might be adequate to approximate the waves within the given flow. Thus, a variation of the MacCormack scheme with fourth order accuracy in space and second order accuracy in time was chosen. A scheme with sixth order accuracy in space has also been implemented and results compared with those of the fourth order accurate scheme. To ensure no mass flux, zero normal velocity is assigned at the plate. This condition will induce a discontinuity in the pressure across the plate location. Values for the perturbation pressure p' along the surface of the plate are obtained using a one-sided, third order Taylor expansion, such that p'(sub y) = O. In accordance with the Kutta condition, perturbation pressure at the trailing edge is assigned to zero. In the far field, radiation boundary conditions have been implemented. The effectiveness of the far field conditions are validated by computing in a larger computational domain and comparing the results. Early in time, sound waves begin to radiate from the plate. A Doppler effect is observed. After the initial transients disappear, the strongest waves leave the trailing edge at an approximate 45 degree angle. The intensity pattern of pressure fluctuations shows five lobes (of increasing magnitude with increasing downstream direction) emerging from the plate. Undesirable short wave contaminants are observed in the computed pressure distribution along the plate surface. For a more accurate solution at small scales, a more refined discretization will be required.

Martin, James E.↗

Numerical simulation of steady and unsteady viscous flow in turbomachinery using pressure based algorithm

The objective of this research is to simulate steady and unsteady viscous flows, including rotor/stator interaction and tip clearance effects in turbomachinery. The numerical formulation for steady flow developed here includes an efficient grid generation scheme, particularly suited to computational grids for the analysis of turbulent turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and convergency in the solution. The numerical model is used to investigate the structure of tip clearance flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to-blade variation of all three velocity components, pitch and yaw angles, losses and blade static pressures in the tip clearance region. The simulation also includes evaluation of such quantities of leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial dissipation, that the tip clearance flow field can be captured accurately. The above numerical formulation was modified to incorporate time accurate solutions. An inner loop iteration scheme is used at each time step to account for the non-linear effects. The computation of unsteady flow through a flat plate cascade subjected to a transverse gust reveals that the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the computation upstream of the stator, and the upstream rotor wake is specified from the experimental data. The results show that the stator potential effects have appreciable influence on the upstream rotor wake. The predicted unsteady wake profiles are compared with the available experimental data and the agreement is good. The numerical results are interpreted to draw conclusions on the unsteady wake transport mechanism in the blade passage.

Lakshminarayana, B.↗

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↗

Lattice gas methods for computational aeroacoustics

This paper presents the lattice gas solution to the category 1 problems of the ICASE/LaRC Workshop on Benchmark Problems in Computational Aeroacoustics. The first and second problems were solved for Delta t = Delta x = 1, and additionally the second problem was solved for Delta t = 1/4 and Delta x = 1/2. The results are striking: even for these large time and space grids the lattice gas numerical solutions are almost indistinguishable from the analytical solutions. A simple bug in the Mathematica code was found in the solutions submitted for comparison, and the comparison plots shown at the end of this volume show the bug. An Appendix to the present paper shows an example lattice gas solution with and without the bug.

Sparrow, Victor W.↗

Jet to freestream velocity ratio computations for a jet in a crossflow

The flowfield induced by a single, subsonic jet exhausting perpendicularly from a flat plate into a subsonic crossflow has been numerically investigated. Time-averaged solutions were obtained using the thin-layer Navier-Stokes equations and two overlapping grids. Test cases were chosen to match available experimental data where the jet Mach number was 0.78 and the freestream Mach number was varied to represent effective velocity ratios R from 4 to 12. Comparisons of the pressures induced on the flat plate are presented. The results show that the best agreement is obtained for R of 4 and 6, R of 8 is more difficult to resolve, and the R of 12 solution was not satisfactory. It is anticipated that the resolution could be improved by using either smaller time steps or finer grid spacing as R increases.

Margason, Richard J.↗

Thermal and chemical convection in planetary mantles

Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Dupeyrat, L.↗

An Extended Lagrangian Method

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Liou, Meng-Sing↗

Aeropropulsion Technology (APT). Task 23 - Stator Seal Cavity Flow Investigation

The focus of NASA Contract NAS3-25950 Task 23 was to numerically investigate the flow through an axial compressor inner-banded stator seal cavity. The Allison/NASA developed ADPAC code was used to obtain all flow predictions. Flow through a labyrinth stator seal cavity of a high-speed compressor was modeled by coupling the cavity flow path and the main flow path of the compressor. A grid resolution study was performed to guarantee adequate grid spacing was used. Both unsteady rotor-stator-rotor interactions and steady-state isolated blade calculations were performed with and without the seal cavity present. A parameterized seal cavity study of the high-speed stator seal cavity collected a series of solutions for geometric variations. The parameter list included seal tooth gap, cavity depth, wheel speed, radial mismatch of hub flowpath, axial trench gap, hub corner treatments, and land edge treatments. Solution data presented includes radial and pitchwise distributions of flow variables and particle traces describing the flow character.

Heidegger, N. J.↗

The Diagnostics of the External Plasma for the Plasma Rocket

The plasma rocket is located at NASA Johnson Space Center. To produce a thrust in space. an inert gas is ionized into a plasma and heated in the linear section of a tokamak fusion device to 1 x 10(exp 4) - 1.16 x 10(exp 6)K(p= 10(exp 10) - 10(exp 14)/cu cm ). The magnetic field used to contain the plasma has a magnitude of 2 - 10k Gauss. The plasma plume has a variable thrust and specific impulse. A high temperature retarding potential analyzer (RPA) is being developed to characterize the plasma in the plume and at the edge of the magnetically contained plasma. The RPA measures the energy and density of ions or electrons entering into its solid angle of collection. An oscilloscope displays the ion flux versus the collected current. All measurements are made relative to the facility ground. A RPA is being developed in a process which involves the investigation of several prototypes. The first prototype has been tested on a thermal plasma. The knowledge gained from its development and testing were applied to the development of a RPA for collimated plasma. The prototypes consist of four equally spaced grids and an ion collector. The outermost grid is a ground. The second grid acts as a bias to repel electrons. The third is a variable v voltage ion suppressor. Grid four (inner grid) acts to repel secondary electrons, being biased equal to the first. Knowledge gained during these two stages are being applied to the development of a high temperature RPA Testing of this device involves the determination of its output parameters. sensitivity, and responses to a wide range of energies and densities. Each grid will be tested individually by changing only its voltage and observing the output from the RPA. To verify that the RPA is providing proper output. it is compared to the output from a Langmuir or Faraday probe.

Karr, Gerald R.↗

Atmospheric Chemical Transport Based on High Resolution Model- Derived Winds: A Case Study

Flight 10 of NASA's Subsonic Assessment (SASS) Ozone and Nitrogen Oxide Experiment (SONEX) extended southwest of Lajes, Azores. A variety of chemical signatures were encountered. These signatures are examined in detail, relating them to meteorological data from a high resolution numerical model having horizontal grid spacing of 30 and 90 km and 26 vertical levels. The meteorological output at hourly intervals is used to create backward trajectories from the locations of the chemical signatures. Four major categories of chemical signatures are discussed-stratospheric, lightning, continental pollution, and a transition layer. The strong stratospheric signal is encountered just south of the Azores in a region of depressed tropopause height. Three chemical signatures at different altitudes in the upper troposphere are attributed to lightning. Backward trajectories arriving at locations of these signatures are related to locations of cloud-to-ground lightning. Results show that the trajectories pass through regions of lightning 1-2 days earlier over the eastern Gulf of Mexico and off the southeast coast of the United States. The lowest leg of the flight exhibits a chemical signature consistent with continental pollution. Trajectories arriving at this signature are found to pass over the highly populated Northeast Corridor of the United States. Surface based pollution apparently is lofted to the altitudes of the trajectories by convective clouds along the East Coast that did not contain lightning. Finally, a chemical transition layer is described. Its chemical signature is intermediate to those of lightning and continental pollution. Trajectories arriving in this layer pass between the trajectories of the lightning and pollution signatures. Thus, they probably are impacted by both sources.

Hannan, John R.↗

A Global Magnetohydrodynamic Model of Jovian Magnetosphere

The goal of this project was to develop a new global magnetohydrodynamic model of the interaction of the Jovian magnetosphere with the solar wind. Observations from 28 orbits of Jupiter by Galileo along with those from previous spacecraft at Jupiter, Pioneer 10 and 11, Voyager I and 2 and Ulysses, have revealed that the Jovian magnetosphere is a vast, complicated system. The Jovian aurora also has been monitored for several years. Like auroral observations at Earth, these measurements provide us with a global picture of magnetospheric dynamics. Despite this wide range of observations, we have limited quantitative understanding of the Jovian magnetosphere and how it interacts with the solar wind. For the past several years we have been working toward a quantitative understanding of the Jovian magnetosphere and its interaction with the solar wind by employing global magnetohydrodynamic simulations to model the magnetosphere. Our model has been an explicit MHD code (previously used to model the Earth's magnetosphere) to study Jupiter's magnetosphere. We continue to obtain important insights with this code, but it suffers from some severe limitations. In particular with this code we are limited to considering the region outside of 15RJ, with cell sizes of about 1.5R(sub J). The problem arises because of the presence of widely separated time scales throughout the magnetosphere. The numerical stability criterion for explicit MHD codes is the CFL limit and is given by C(sub max)(Delta)t/(Delta)x less than 1 where C(sub max) is the maximum group velocity in a given cell, (Delta)x is the grid spacing and (Delta)t is the time step. If the maximum wave velocity is C(sub w) and the flow speed is C(sub f), C(sub max) = C(sub w) + C(sub f). Near Jupiter the Alfven wave speed becomes very large (it approaches the speed of light at one Jovian radius). Operating with this time step makes the calculation essentially intractable. Therefore under this funding we have been designing a new MHD model that will be able to compute solutions in the wide parameter regime of the Jovian magnetosphere.

Walker, Raymond J.↗

A Cloud-Resolving Simulation of Hurricane Bob (1991): Storm Structure and Eyewall Buoyancy

A numerical simulation of Hurricane Bob (1991) is conducted using the Penn State University-National Center for Atmospheric Research mesoscale model MM5 with a horizontal grid spacing of 1.3 Km on the finest nested mesh The model produces a realistic hurricane that intensifies slowly during the period of fine-scale simulation. Time-averaged results reveal the effects of storm motion. vertical shear, beta gyres and deformation forcing on the structure of radial inflow, vertical motion, and precipitation. Instantaneous model fields show that radial inflow in the eyewall is very intense near the surface but transitions to strong low-level outflow near the top of the boundary layer. The low-level structure is modulated by a wavenumber 2 disturbance that rotates around the eyewall at half the speed of the maximum tangential winds and is consistent with a vortex Rossby edge wave. The statistical distribution of vertical velocity in the eyewall indicates that the eyewall is composed of a small number of intense updrafts that account for the majority of the upward mass flux rather than a more gradual and symmetric eyewall circulation, consistent with the concept of hot towers. Tongues of high equivalent potential temperature, Theta(sub e), are seen along the inner edge of the eyewall updraft and within the low-level outflow. This air originates from outside of the eyewall with the highest theta(sub e) air coming from the layer closest to the surface after penetrating closest to the center. Occasionally, high Theta(sub e), air within the eye is drawn into the eyewall updrafts. The high Theta(sub e), air rising within the eyewall is shown to be associated with positive eyewall buoyancy with sufficient convective available potential energy along its path to produce relatively strong convective updrafts. Although the requirements for conditional symmetric instability are met within the eyewall and the air parcel trajectories follow slanted paths, the radial displacement of air parcels in the low-level outflow moves the air parcel sufficiently far away from the upper- warm core that the air becomes unstable to vertical displacements. Hence, convective instability rather than symmetric instability accounts for the stronger updrafts in the eyewall.

Braun, Scott A.↗