Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiple scale”

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 325 records · Page 18

Numerical simulation of boundary layers. Part 3: Turbulence and relaminarization in sink flows

Direct numerical simulations of sink-flow boundary layers, with acceleration parameters K between 1.5 and 3.0 x 10 to the -6 power, are presented. The three-dimensional, time-dependent Navier-Stokes equations are solved numerically using a spectral method, with about 10 to the -6 degrees of freedom. The flow is assumed to be statistically steady, and self-similar. A multiple-scale approximation and periodic conditions are applied to the fluctuations. The turbulence is studied using instantaneous and statistical results. Good agreement with the experiments of Jones and Launder is observed. Two effects of the favorable pressure gradient are to extend the logarithmic layer, and to alter the energy balance of the turbulence near the edge of the boundary layer. At very low Reynolds number the logarithmic layer is shortened and slightly displaced, but wall-layer streaks are present even at the lowest values of the Reynolds number for which turbulence can be sustained. Large quiescent patches appear in the flow. Relaminarization occurs at K = 3.0 x 10 to the -6 power, corresponding to a Reynolds number, based on momentum thickness, of about 330.

Spalart, Philippe R.↗

Direct simulation of a turbulent boundary layer up to R sub theta = 1410

The turbulent boundary layer on a flat plate, with zero pressure gradient, is simulated numerically at four stations between R sub theta = 225 and R sub theta = 1410. The three-dimensional time-dependent Navier-Stokes equations are solved using a spectral method with up to about 10 to the 7th grid points. Periodic spanwise and streamwise conditions are applied, and a multiple-scale procedure is applied to approximate the slow streamwise growth of the boundary layer. The flow is studied, primarily, from a statistical point of view. The solutions are compared with experimental results. The scaling of the mean and turbulent quantities with Reynolds number is compared with accepted laws, and the significant deviations are documented. The turbulence at the highest Reynolds number is studied in detail. The spectra are compared with various theoretical models. Reynolds-stress budget data are provided for turbulence-model testing.

Spalart, Philippe R.↗

Computational simulation of turbulent flow in Space Shuttle Main Engine turnaround ducts

Two axisymmetric 180-deg turnaround ducts are used in the Space Shuttle Main Engine (SSME) to connnect the preburners with the main thrust chamber. The prediction of the flow field and heat transfer within the turnaround ducts is important to minimize the pressure drops, size of recirculation and stagnation zones, local overheating, etc. This paper presents the computational approach and selected results for the turbulent flow in the turnaround ducts of the fuel and oxidizer sides of the SSME. The time-averaged, Navier-Stokes equations for the viscous, compressible, turbulent flow are solved in body-fitted-coordinates by using a finite-volume approach. Two turbulence models, viz: the k-epsilon model and a multiple scale turbulence model, are used to examine the sensitivity of calculated flows. Both models produce almost identical solutions for the fuelside turnaround duct (which has no recirculation region). However, for the oxidizer-side duct which has a large recirculation region, the two models show quite different results.

Przekwas, Andrzej J.↗

A spectral multi-domain technique applied to high-speed chemically reacting flows

The first applications of a spectral multidomain method for viscous compressible flow is presented. The method imposes a global flux balance condition at the interface so that high-order continuity of the solution is preserved. The global flux balance is imposed in terms of a spectral integral of the discrete equations across adjoining domains. Since the discretized equations interior to each domain solved are uncoupled from each other, and since the interface relation has a block structure, the solution scheme can be adapted to the particular requirements of each subdomain. The spectral multidomain technique presented is well-suited for the multiple scales associated with the chemically reacting and transition flows in hypersonic research. A nonstaggered multidomain discretization is used for the chemically reacting flow calculation, and the first implementation of a staggered multidomain mesh is presented for accurately solving the stability equation for a viscous compressible fluid.

Macaraeg, Michele G.↗

On WKB expansions for Alfven waves in the solar wind

The WKB expansion for 'toroidal' Alfven waves in solar wind, which is described by equations of Heinemann and Olbert (1980), is examined. In this case, the multiple scales method (Nayfeh, 1981) is used to obtain a uniform expansion. It is shown that the WKB expansion used by Belcher (1971) and Hollweg (1973) for Alfven waves in the solar wind is nonuniformly convergent.

Hollweg, Joseph V.↗

Non-linear dynamic analysis of geared systems

Under driving conditions, a typical geared system may be subjected to large dynamic loads. Also, the vibration level of the geared system is directly related to the noise radiated from the gear box. The steady state dynamic behavior of the system is examined in order to design reliable and quiet transmissions. The scope is limited to a system containing a spur gear pair with backlash and periodically time varying mesh stiffness, and rolling element bearings with clearance type nonlinearities. The internal static transmission error at the gear mesh, which is of importance from high frequency noise and vibration control view point, is considered in the formulation in sinusoidal or periodic form. A dynamic finite element model of the linear time invariant (LTI) system is developed. Effects of several system parameters, such as torsional and transverse flexibilities of the shafts and prime mover/load inertias, on free and forced vibration characteristics are investigated. Several reduced order LTI models are developed and validated by comparing their eigen solutions with the finite element model results. Using the reduced order formulations, a three degree of freedom dynamic model is developed which includes nonlinearities associated with radical clearances in the radial rolling element bearings, backlash between a spur gear pair and periodically varying gear mesh stiffness. As a limiting case, a single degree of freedom model of the spur gear pair with backlash is considered and mathematical conditions for tooth separation and back collision are defined. Both digital simulation technique and analytical models such as method of harmonic balance and the method of multiple scales were used to develop the steady state frequency response characteristics for various nonlinear and/or time varying cases.

Singh, Rajendra↗

Nonparallel instability of supersonic and hypersonic boundary layers

Multiple scaling technique is used to examine the nonparallel instability of supersonic and hypersonic boundary-layer flows to three dimensional (first mode) and two dimensional (second mode) disturbances. The method is applied to the flat plate boundary layer for a range of Mach numbers from 0 to 10. Growth rates of disturbances are calculated based on three different criteria: following the maximum of the mass-flow disturbance, using an integral of the disturbance kinetic energy, and using the integral of the square of the mass-flow amplitude. By following the maximum of the mass-flow disturbance, the calculated nonparallel growth rates are in good quantitative agreement with the experimental results at Mach number 4.5.

El-Hady, Nabil M.↗

Nonparallel instability of supersonic and hypersonic boundary layers

Multiple scaling technique is used to examine the nonparallel instability of supersonic and hypersonic boundary-layer flows to three-dimensional (first mode) and two-dimensional (second mode) disturbances. The method is applied to the flat plate boundary layer for a range of Mach numbers from 0 to 10. Growth rates of disturbances are calculated based on three different criteria: following the maximum of the mass-flow disturbance, using an integral of the disturbance kinetic energy, and using an integral of the square of the mass-flow amplitude. By following the maximum of the mass-flow dusturbance, the calculated nonparallel growth rates are in good quantitative agreement with the experimental results of Kendall (1967) at Mach number 4.5.

El-Hady, Nabil M.↗

Modulational stability of rotating-disk flow

The problem of modulational primary stability of rotating-disk boundary layer flow for stationary modes was studied by using the method of multiple scales. The system of equations for infinitesimal disturbances is solved by linear stability theory to determine the modulated solution which represents the time and azimuthal variations of the wave pattern. The solution for the amplitude of the most critical disturbance indicates that the disturbance spreads out as the flow time increases and grows exponentially at the center.

Riahi, Daniel N.↗

Waves produced from a harmonic point source in a supersonic boundary layer

The disturbance wave pattern produced by a harmonic point source in a compressible flat-plate boundary layer is computed using linear stability theory and direct numerical integration approach. Receptivity coefficients are computed for the spectrum of spanwise modes generated at the source. The effect of boundary layer growth on the development of linear waves is determined by using the method of multiple scales. Results are presented for Mach numbers of 0, 2, 4.5, and 7. It is found that disturbances spread in wedge-shaped regions behind the source and the wedge angle decreases with Mach number. The lateral spreading angle for the instability waves turns out to be quite close to the angle found experimentally for turbulence lateral contamination.

Balakumar, P.↗

Nonparallel instability of supersonic and hypersonic boundary layers

Multiple scaling technique is used to examine the nonparallel instability of supersonic and hypersonic boundary-layer flows to three-dimensional (first mode) and two-dimensional (second mode) disturbances. The method is applied to the flat plate boundary layer for a range of Mach numbers from 0 to 10. Growth rates of disturbances are calculated based on three different criteria: following the maximum of the mass-flow disturbance, using an integral of the disturbance kinetic energy, and using the integral of the square of the mass-flow amplitude. By following the maximum of the mass-flow disturbance, the calculated nonparallel growth rates are in good quantitative agreement with the experimental results at Mach number 4.5.

El-Hady, Nabil M.↗

Turbine blade tip and seal clearance excitation forces

Experimental and theoretical work done as Phase 3 of a program sponsored by MSFC to investigate the magnitude, origin, and parametric variations of destabilizing forces which arise in high power turbines due to blade-tip leakage effects are described. The two facilities which were built for this purpose are first described. The larger one is a closed, 2 atm pressurized Freon-12 flow loop in which is installed a 1:1 replica of the SSME first stage hydrogen turbine, which can be driven by the flow, and which generates about 14 KW of power into a load-absorbing DC generator. The smaller facility is used to measure the forces on labyrinth seals of the same type as those used in our turbine tests with a shrouded turbine. The seals can be kinematically whirled and spun (independently), and the inlet swirl can be set to a variety of values. Air is the working fluid (with atmospheric discharge) and the data are real-time pressure distributions in the seal glands. The five different unshrouded turbine configurations were tested with static offsets, plus one with a shroud band and a two-ridge seal. Theoretical models of various degrees of complexity were developed to help interpreting and extrapolating the data. The notion of partial work done by the fluid leaking through the tip gaps was put on a quantitative basis by examining the leakage vortex roll-up dynamics. This was used to obtain a theory of the work loss due to a uniform gap. Perturbation and multiple scale arguments were then used to extend this to the case of an eccentric turbine. This yields an unsteady, 3-D theory which can predict the distribution of the approach flow, and its effect on work defect, cross-forces, pressure patterns, and dynamic damping. The predictions agree qualitatively with the data and exhibit the correct trends, but the cross-forces are generally under-predicted.

Martinez-Sanchez, M.↗

Nonlinear response of infinitely long circular cylindrical shells to subharmonic radial loads

The nonlinear response of infinitely long circular cylindrical shells (thin circular rings) in the presence of a two-to-one internal (autoparametric) resonance to a subharmonic excitation of order one-half of the higher mode is analyzed with the multiple-scale method. Four autonomous first-order ordinary differential equations are derived for the modulation of the amplitudes and phases of the interacting models. These modulation equations are used to determine the fixed points and their stability. The fixed points correspond to periodic oscillations of the shell, whereas the limit-cycle solutions of the modulation equations correspond to amplitude and phase-modulated oscillations of the shell. The force response curves exhibit saturation, jumps, and Hopf bifurcation. As excitation frequency changes, all limit cycles deform and lose stability through either pitchfork or cyclic-fold (saddle-node) bifurcations. Some of these saddle-node bifurcations cause a transition to chaos. The pitchfork bifurcations break the symmetry of the limit cycles.

Nayfeh, Ali H.↗

Wave diffraction in weak cosmic-ray-modified shocks

Weakly multidirectional, long-wavelength cosmic-ray-modified shocks are studied via multiple scale perturbation techniques. The effects of diffraction are discussed in terms of Green's function solutions of the linearized 1 + 3D Burgers and 1 + 3D KdVB equations, and also in terms of solutions with singular Dirac delta initial distributions. The solutions show a monotonic decrease of the wave-front curvature with increasing time owing to the effects of wave diffraction. The shape of the wave surface is discussed in terms of solutions S to the wave eikonal equation corresponding to singular initial conditions. For the fast magnetosonic wave propagating in the positive x-direction, the wave phase surface S = 0 has elliptic cross sections with the planes x = constant and has a convex paraboloidal shape. Plane-wave solutions of the 1 + 3D KdVB equation are discussed.

Webb, G. M.↗

Examination of various turbulence models for application in liquid rocket thrust chambers

There is a large variety of turbulence models available. These models include direct numerical simulation, large eddy simulation, Reynolds stress/flux model, zero equation model, one equation model, two equation k-epsilon model, multiple-scale model, etc. Each turbulence model contains different physical assumptions and requirements. The natures of turbulence are randomness, irregularity, diffusivity and dissipation. The capabilities of the turbulence models, including physical strength, weakness, limitations, as well as numerical and computational considerations, are reviewed. Recommendations are made for the potential application of a turbulence model in thrust chamber and performance prediction programs. The full Reynolds stress model is recommended. In a workshop, specifically called for the assessment of turbulence models for applications in liquid rocket thrust chambers, most of the experts present were also in favor of the recommendation of the Reynolds stress model.

Hung, R. J.↗

Far-infrared emission, gas, and ultraluminous HII regions in M101

We have made deconvolved maps of M101 in each of the 4 IRAS bands using a new maximum-entropy based model for reconstructing distributions with correlated structure on multiple scales. This new deconvolution procedure, developed by Weir (1992, J.Opt.Soc.Am., submitted), is superior to other maximum entropy-based techniques for several reasons. For our purposes, an important advantage is that it has less artifacts and greatly reduced systematic biases compared with the 'HiRes' maps of M101 produced by Aumann, Fowler, and Melnyk (1990, AJ, 99, 1674), thus it is accurate enough to perform aperture photometry. The deconvolved maps have a maximum resolution of approximately 30 sec, sufficient to resolve the brightest HII complexes and much of the spiral structure. The new maps and our comparison of the far-infrared, gas, and optical distributions are being written up in Kenney, Weir & Scoville (1992, in preparation). The quality of the reconstructed far-infrared maps is good enough for us to carry out most of the analyses outlined in the original proposal. From a ratio of the 60 micron and 100 micron maps, we have found that the ultraluminous HII complexes in the outer galaxy have the hottest dust temperatures with T approximately 50-60 K, which is twice as hot as most of the disk. Their extraordinary luminosity in the far-infrared is due to a relatively small amount of dust being heated to high temperatures, rather than a large concentration of dust and gas. A map of the dust opacity at 60 microns shows good overall agreement with a map of cold gas (HI+H2), indicating that throughout most of the galaxy only approximately 20 percent of the dust is warm enough to be detected by IRAS. The ultraluminous HII complexes have a high luminosity-to-gas mass ratio, independent of whether the far-infrared or the H-alpha emission line is used to measure the luminosity, which implies that gas is being converted into high mass stars more rapidly in these complexes compared to other locations in the disk.

Kenney, J. D. P.↗

An unusual spectral unit in West Candor Chasma: Evidence for hydrothermal or aqueous alteration?

A spectrally distinctive unit on the floor of W. Candor Chasma (6 S, 76 W) in the central Valles Marineris may be a likely candidate for hydrothermal or aqueous alteration. This unusual material is noticeably redder than nearby plains and canyon floor-covering deposits of similar brightness in several Viking Orbiter color composite images calibrated using PICS Level I procedures. The surrounding plains and canyon floor units have colors that are typical for much of the weathered soil on Mars. Relative to adjacent materials, the West Candor unit has lower green-filter reflectance and higher red-filter reflectance. While subtle, these spectral characteristics were observed for this unit in a number of multi-spectral images acquired at different seasons and phase angles. When the color image data are transformed to hue, saturation, and value coordinates, the West Candor material stands out prominently among the rocks and soils of Coprates Quadrangle as a spatially coherent unit with a unique hue. Physically, this means that the unit is compositionally distinct (unlike most of the bright materials in the region), since its spectral reflectance cannot be obtained through simple multiplicative scaling of the reflectance of the surrounding bright materials or by altering the spectral reflectance of the surrounding materials with an additive constant that is independent of wavelength. Recent studies of the directional reflectance properties of the surface in this region indicate that the photometric phase function of the unit is similar to that of surrounding bright materials (isotropic to slightly back scattering).

Geissler, P. E.↗

Pulsatile instability in rapid directional solidification: Strongly-nonlinear analysis

In models of rapid directional solidification, non-equilibrium interfacial conditions are employed. As a result, there is an oscillatory mode of instability, as well as the steady cellular mode, found in the equilibrium model of Mullins and Sekerka. When the temperature field is decoupled from the problem, the preferred wave number for the oscillatory mode is zero, and the interface pulsates in time while remaining spatially uniform. Results from multiple scale analyses in the two limiting cases of the parameters are reported. In these limits, it is found that the instability is a bifurcation to relaxation oscillations; these nonlinear oscillations may be related to the observed microstructure that results from rapid solidification processes such as laser surface remelting.

Braun, Richard J.↗