SEARCH · Search NASA

Results for “Hydrodynamic stability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1985

A spectral collocation method for the Navier-Stokes equations

A Fourier-Chebyshev spectral method for the incompressible Navier-Stokes equations is described. It is applicable to a variety of problems including some with fluid properties which vary strongly both in the normal direction and in time. In this fully spectral algorithm, a preconditioned iterative technique is used for solving the implicit equations arising from semi-implicit treatment of pressure, mean advection and vertical diffusion terms. The algorithm is tested by applying it to hydrodynamic stability problems in channel flow and in external boundary layers with both constant and variable viscosity.

Malik, M. R.↗

NASA NTRS1986

Instability of elliptic jets

This paper presents a detailed account of the hydrodynamic stability characteristics of the initial mixing region of an elliptic jet. A realistic mean velocity profile is used. Calculations of growth rates, phase velocities and eigenfunctions are presented. The growth rates of all modes in the initial mixing region scale with the minimum momentum thickness. Pressure fluctuations are found to be greatest for all modes close to the major axis. An irregular normal mode is found at larger eccentricities. The flapping mode about the major axis becomes the least stable as the momentum thickness of the mixing layer increases.

Morris, P. J.↗

NASA NTRS1986

Numerical and physical aspects of aerodynamic flows III; Proceedings of the Third Symposium, California State University, Long Beach, CA, January 21-24, 1985

The present symposium considers hydrodynamic stability and turbulent transition, transition calculations in three-dimensional flows, a quasi-simultaneous finite difference approach for strongly interacting flows, the significance of the thin layer Navier-Stokes approximation, unsteady airfoil boundary layers, predictions and experiments on airfoils at low Reynolds numbers, and a comparison of interactive boundary layer and thin layer Navier-Stokes procedures. Also discussed are a viscous-inviscid interaction method for computing unsteady transonic separation, massive separation and dynamic stall on a cusped trailing edge airfoil, the computation of turbulent separated flows over wings, an iterative scheme for three-dimensional transonic flows, and the computation of three-dimensional flows with shock wave/boundary layer interaction. Many of the papers in this symposium were abstracted previously (cf., A85-42951).

Cebeci, T.↗

NASA NTRS1988

Inflow/outflow boundary conditions and global dynamics of spatial mixing layers

The numerical simulation of incompressible spatially-developing shear flows poses a special challenge to computational fluid dynamicists. The Navier-Stokes equations are elliptic and boundary equations need to be specified at the inflow and outflow boundaries in order to compute the fluid properties within the region of interest. It is, however, difficult to choose inflow and outflow conditions corresponding to a given experimental situation. Furthermore the effects that changes in the boundary conditions or in the size of the computational domain may induce on the global dynamics of the flow are presently unknown. These issues are examined in light of recent developments in hydrodynamic stability theory. The particular flow considered is the spatial mixing layer but it was expected that similar phenomena were bound to occur in other cases such as channel flow, the boundary layer, etc. A short summary of local/global and absolute/convective instability concepts is given. The results of numerical simulations are presented which strongly suggest that global resonances may be triggered in domains of finite streamwise extent although the evolution of the perturbation vorticity field is everywhere locally convective. A relationship between finite domains and pressure sources which might help in devising a scheme to eliminate these difficulties is discussed.

Buell, J. C.↗

NASA NTRS1989

Elimination of spurious eigenvalues in the Chebyshev tau spectral method

Spectral methods have been used to great advantage in hydrodynamic stability calculations; the concepts are described in Orszag's seminal application of the Chebyshev tau method to the Orr-Sommerfeld equation for plane Poiseuille flow in 1971. Orszag discusses both the Chebyshev Galerkin and the Chebyshev tau methods, but presents results for the tau method, which is easier to implement than the Galerkin method. The tau method has the disadvantage that two unstable eigenvalues are produced that are artifacts of the discretization. An extremely simple modification to the Chebyshev tau method is presented which eliminates the spurious eigenvalues. First a simplified model of the Orr-Sommerfeld equation discussed by Gottlieb and Orszag was studied. Then the Chebyshev tau method is considered, which has two spurious eigenvalues, and then a modification which eliminates them is described. Finally, results for the Orr-Sommerfeld equation are considered where the modified tau method also eliminates the spurious eigenvalues. The simplicity of the modification makes it a convenient alternative to other approaches to the problem.

Mcfadden, G. B.↗

NASA NTRS1990

Modern developments in shear flow control with swirl

Passive and active control of swirling turbulent jets is experimentally investigated. Initial swirl distribution is shown to dominate the free jet evolution in the passive mode. Vortex breakdown, a manifestation of high intensity swirl, was achieved at below critical swirl number (S = 0.48) by reducing the vortex core diameter. The response of a swirling turbulent jet to single frequency, plane wave acoustic excitation was shown to depend strongly on the swirl number, excitation Strouhal number, amplitude of the excitation wave, and core turbulence in a low speed cold jet. A 10 percent reduction of the mean centerline velocity at x/D = 9.0 (and a corresponding increase in the shear layer momentum thickness) was achieved by large amplitude internal plane wave acoustic excitation. Helical instability waves of negative azimuthal wave numbers exhibit larger amplification rates than the plane waves in swirling free jets, according to hydrodynamic stability theory. Consequently, an active swirling shear layer control is proposed to include the generation of helical instability waves of arbitrary helicity and the promotion of modal interaction, through multifrequency forcing.

Farokhi, Saeed↗

NASA NTRS1988

Contributions to the understanding of large-scale coherent structures in developing free turbulent shear flows

The physical problem of large-scale coherent structures in real, developing free turbulent shear flows are discussed from the point of view of a broader interpretation of the nonlinear aspects of hydrodynamic stability. Variations on the Amsden and Harlow problem are considered, and the role of linear theory in nonlinear problems is addressed. Spatially developing two-dimensional coherent structures and three-dimensional nonlinear effects in large-scale coherent mode interactions are considered.

Liu, J. T. C.↗

NASA NTRS1990

Theory of electrohydrodynamic instabilities in electrolytic cells

The paper develops the theory of the hydrodynamic stability of an electrolytic cell as a function of the imposed electric current. A new electrohydrodynamic instability is encountered when the current is forced to exceed the Nernst limit. The convection is driven by the volume force exerted by the electric field on space charges in the electrolyte. This intrinsic instability is found to be easily masked by extrinsic convection sources such as gravity or stirring. A linear stability analysis is performed and a dimensionless number Le is derived whose value determines the convection pattern.

Bruinsma, R.↗

NASA NTRS1990

The vibrating ribbon problem revisited

A revised formal solution of the vibrating ribbon problem of hydrodynamic stability is presented. The initial formulation of Gaster (1965) is modified by application of the Briggs method and a careful treatment of the complex double Fourier transform inversions. Expressions are obtained in a natural way for the discrete spectrum as well as for the four branches of the continuous spectra. These correspond to discrete and branch-cut singularities in the complex wavenumber plane. The solutions from the continuous spectra decay both upsteam and downstream of the ribbon, with the decay in the upstream direction being much more rapid than that in the downstream direction. Comments and clarification of related prior work are made.

Ashpis, David E.↗

NASA NTRS1989

Computational fluid dynamics analysis of space shuttle main propulsion feed line 17-inch disconnect valves

A steady incompressible three-dimensional (3-D) viscous flow analysis was conducted for the Space Shuttle Main Propulsion External Tank (ET)/Orbiter (ORB) propellant feed line quick separable 17-inch disconnect flapper valves for liquid oxygen (LO2) and liquid hydrogen (LH2). The main objectives of the analysis were to predict and correlate the hydrodynamic stability of the flappers and pressure drop with available water test data. Computational Fluid Dynamics (CFD) computer codes were procured at no cost from the public domain, and were modified and extended to carry out the disconnect flow analysis. The grid generator codes SVTGD3D and INGRID were obtained. NASA Ames Research Center supplied the flow solution code INS3D, and the color graphics code PLOT3D. A driver routine was developed to automate the grid generation process. Components such as pipes, elbows, and flappers can be generated with simple commands, and flapper angles can be varied easily. The flow solver INS3D code was modified to treat interior flappers, and other interfacing routines were developed, which include a turbulence model, a force/moment routine, a time-step routine, and initial and boundary conditions. In particular, an under-relaxation scheme was implemented to enhance the solution stability. Major physical assumptions and simplifications made in the analysis include the neglect of linkages, slightly reduced flapper diameter, and smooth solid surfaces. A grid size of 54 x 21 x 25 was employed for both the LO2 and LH2 units. Mixing length theory applied to turbulent shear flow in pipes formed the basis for the simple turbulence model. Results of the analysis are presented for LO2 and LH2 disconnects.

Kandula, Max↗

NASA NTRS1990

Reynolds stress closure in jet flows using wave models

A collection of papers is presented. The outline of this report is as follows. Chapter three contains a description of a weakly nonlinear turbulence model that was developed. An essential part of the application of such a closure scheme to general geometry jets is the solution of the local hydrodynamic stability equation for a given jet cross-section. Chapter four describes the conformal mapping schemes used to map such geometries onto a simple computational domain. Chapter five describes a solution of a stability problem for circular, elliptic, and rectangular geometries. In chapter six linear models for the shock shell structure in non-circular jets is given. The appendices contain reprints of papers also published during this study including the following topics: (1) instability of elliptic jets; (2) a technique for predicting the shock cell structure in non-circular jets using a vortex sheet model; and (3) the resonant interaction between twin supersonic jets.

Morris, Philip J.↗

NASA NTRS1991

A comparison of numerical methods for the Rayleigh equation in unbounded domains

A second-order finite difference and two spectral methods, including a Chebyshev tau and a Chebyshev collocation method were implemented to determine the linear hydrodynamic stability of an unbounded shear flow. The velocity profile of the basic flow in the stability analysis mimicks that of a one-stream free mixing layer. Local and global eigenvalue solution methods are used to determine individual eigenvalues and the eigenvalue spectrum, respectively. The calculated eigenvalue spectrum includes a discrete mode, a continuous spectrum associated with the equation singularity and a continuous spectrum associated with the domain unboundedness. The efficiency and the accuracy of these discretization methods in the prediction of the eigensolutions of the discrete mode were evaluated by comparison with a conventional shooting procedure. Their capabilities in mapping out the continuous eigenvalue spectra are also discussed.

Liou, W. W.↗

NASA NTRS1990

Modern developments in shear flow control with swirl

Farokhi, S.↗

NASA NTRS1991

On spatial instability of a swirling jet - Theory and experiment

The spatial instability of a swirling jet is investigated both experimentally and theoretically. A hydrodynamic stability analysis is applied to an inviscid incompressible top-hat jet, with a swirl distribution of solid-body rotation and free vortex in and outside the vortex core, respectively. Both plane and helical instability modes are examined. It is found that the top-hat jet with swirl distribution of Rankine vortex type is unstable in all the modes studied. The higher the positive helicity, the less spatially unstable the jet behavior; the higher the negative helicity, the more spatially unstable this behavior becomes. A comparison is made between theoretical results and experimental data on a low-intensity swirling jet. The trend of the initial growth of the instability waves in the near field is captured.

Wu, C.↗

NASA NTRS1993

Linear instability of curved free shear layers

The linear inviscid hydrodynamic stability of slightly curved free mixing layers is studied in this paper. The disturbance equation is solved numerically using a shooting technique. Two mean velocity profiles that represent stably and unstably curved free mixing layers are considered. Results are shown for cases of five curvature Richardson numbers. The stability characteristics of the shear layer are found to vary significantly with the introduction of the curvature effects. The results also indicate that, in a manner similar to the Goertler vortices observed in a boundary layer along a concave wall, instability modes of spatially developing streamwise vortex pairs may appear in unstable curved mixing layers.

Liou, William W.↗

NASA NTRS1993

Linear instability of curved free shear layers

Liou, William W.↗

NASA NTRS1995

The Benard problem: A comparison of finite difference and spectral collocation eigen value solutions

The application of spectral methods, using a Chebyshev collocation scheme, to solve hydrodynamic stability problems is demonstrated on the Benard problem. Implementation of the Chebyshev collocation formulation is described. The performance of the spectral scheme is compared with that of a 2nd order finite difference scheme. An exact solution to the Marangoni-Benard problem is used to evaluate the performance of both schemes. The error of the spectral scheme is at least seven orders of magnitude smaller than finite difference error for a grid resolution of N = 15 (number of points used). The performance of the spectral formulation far exceeded the performance of the finite difference formulation for this problem. The spectral scheme required only slightly more effort to set up than the 2nd order finite difference scheme. This suggests that the spectral scheme may actually be faster to implement than higher order finite difference schemes.

Skarda, J. Raymond Lee↗

NASA NTRS1995

Combustion of two-component miscible droplets in reduced gravity

This research focuses on the combustion of binary miscible droplets initially in the mm size range. Experiments are performed using the NASA Lewis 2.2 sec drop tower in Cleveland, Ohio, where mixtures of alkanes and/or alcohols are studied. The fuel components are selected to have significantly different volatilities. Initial oxygen mole fractions from about 0.15-0.5 and initial pressures from 0.2-2 atm are employed. Different inerts are used (He, CO2, Ar, N2) to change burning rates and sooting behaviors. Objectives are to observe the following: (1) Transient droplet diameters (including three-staged combustion behaviors and microexplosion; (2) Transient flow behaviors (sudden flame contraction, luminosity, extinction); and (3) Behaviors of observable soot particles. theoretical and computational research in support of this program has also been undertaken. This research includes analytical studies to determine the effects of small but nonzero gravitational levels on droplet gasification, analytical studies of hydrodynamic stability of spherically-symmetrical droplet gasification (to address the question as to whether spherically-symmetrical droplet gasification may be destabilized from capillary, i.e., Marangoni effects), and computational modeling of effects of capillary stresses on droplet gasification.

Shaw, Benjamin D.↗