SEARCH · Search NASA
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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).
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.
Hierarchical computational procedures are described to simulate the multiple scale thermal/mechanical behavior of high temperature metal matrix composites (HT-MMC) in the following three broad areas: (1) behavior of HT-MMC's from micromechanics to laminate via METCAN (Metal Matrix Composite Analyzer), (2) tailoring of HT-MMC behavior for optimum specific performance via MMLT (Metal Matrix Laminate Tailoring), and (3) HT-MMC structural response for hot structural components via HITCAN (High Temperature Composite Analyzer). Representative results from each area are presented to illustrate the effectiveness of computational simulation procedures and accompanying computer codes. The sample case results show that METCAN can be used to simulate material behavior such as the entire creep span; MMLT can be used to concurrently tailor the fabrication process and the interphase layer for optimum performance such as minimum residual stresses; and HITCAN can be used to predict the structural behavior such as the deformed shape due to component fabrication. These codes constitute virtual portable desk-top test laboratories for characterizing HT-MMC laminates, tailoring the fabrication process, and qualifying structural components made from them.
The main research activities at the Center for Modeling of Turbulence and Transition (CMOTT) are described. The research objective of CMOTT is to improve and/or develop turbulence and transition models for propulsion systems. The flows of interest in propulsion systems can be both compressible and incompressible, three dimensional, bounded by complex wall geometries, chemically reacting, and involve 'bypass' transition. The most relevant turbulence and transition models for the above flows are one- and two-equation eddy viscosity models, Reynolds stress algebraic- and transport-equation models, pdf models, and multiple-scale models. All these models are classified as one-point closure schemes since only one-point (in time and space) turbulent correlations, such as second moments (Reynolds stresses and turbulent heat fluxes) and third moments, are involved. In computational fluid dynamics, all turbulent quantities are one-point correlations. Therefore, the study of one-point turbulent closure schemes is the focus of our turbulence research. However, other research, such as the renormalization group theory, the direct interaction approximation method, and numerical simulations are also pursued to support the development of turbulence modeling.
This research brief contains the progress reports of the research staff of the Center for Modeling of Turbulence and Transition (CMOTT) from June 1992 to July 1993. It is also an annual report to the Institute for Computational Mechanics in Propulsion located at Ohio Aerospace Institute and NASA Lewis Research Center. The main objectives of the research activities at CMOTT are to develop, validate, and implement turbulence and transition models for flows of interest in propulsion systems. Currently, our research covers eddy viscosity one- and two-equation models, Reynolds-stress algebraic equation models, Reynolds-stress transport equation models, nonequilibrium multiple-scale models, bypass transition models, joint scalar probability density function models, and Renormalization Group Theory and Direct Interaction Approximation methods. Some numerical simulations (LES and DNS) have also been carried out to support the development of turbulence modeling. Last year was CMOTT's third year in operation. During this period, in addition to the above mentioned research, CMOTT has also hosted the following programs: an eighteen-hour short course on 'Turbulence--Fundamentals and Computational Modeling (Part I)' given by CMOTT at the NASA Lewis Research Center; a productive summer visitor research program that has generated many encouraging results; collaborative programs with industry customers to help improve their turbulent flow calculations for propulsion system designs; a biweekly CMOTT seminar series with speakers from within and without the NASA Lewis Research Center including foreign speakers. In addition, CMOTT members have been actively involved in the national and international turbulence research activities. The current CMOTT roster and organization are listed in Appendix A. Listed in Appendix B are the abstracts of the biweekly CMOTT seminar. Appendix C lists the papers contributed by CMOTT members.
In this paper the nonlinear stability of two-phase core-annular flow in a pipe is examined when the acting pressure gradient is modulated by time harmonic oscillations and viscosity stratification and interfacial tension is present. An exact solution of the Navier-Stokes equations is used as the background state to develop an asymptotic theory valid for thin annular layers, which leads to a novel nonlinear evolution describing the spatio-temporal evolution of the interface. The evolution equation is an extension of the equation found for constant pressure gradients and generalizes the Kuramoto-Sivashinsky equation with dispersive effects found by Papageorgiou, Maldarelli & Rumschitzki, Phys. Fluids A 2(3), 1990, pp. 340-352, to a similar system with time periodic coefficients. The distinct regimes of slow and moderate flow are considered and the corresponding evolution is derived. Certain solutions are described analytically in the neighborhood of the first bifurcation point by use of multiple scales asymptotics. Extensive numerical experiments, using dynamical systems ideas, are carried out in order to evaluate the effect of the oscillatory pressure gradient on the solutions in the presence of a constant pressure gradient.