SEARCH · Engineering Papers
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 higher spatial-stability eigenmodes for the Blasius boundary layer are examined by using asymptotic theory, and an infinite number of modes are found. The asymptotic results are shown to be in close agreement with results from a direct numerical solution of the Orr-Sommerfeld problem. The asymptotic theory would therefore provide an efficient tool in exploratory searches for the eigenvalues.
In this paper numerical calculation of the spatial stability of disturbances in the parallel and nonparallel Blasius boundary layers is considered. Chebyshev polynomials are used for discretization. The problem with the boundary condition at infinity is overcome, and the resulting nonlinear matrix eigenvalue problem is attacked directly. The secondary eigenvalue problem for three-dimensional disturbances is shown to be uniformly stable, and particular solutions of this problem generated by the Orr-Sommerfeld equation are shown. A numerical solution of the nonparallel problem is considered using Chebyshev polynomials. The matrix equations are analyzed directly and the problem of uniqueness of the nonparallel correction is settled by careful application of the Fredholm alternative. Nonparallel corrections to the streamwise eigenfunction are shown.
In this paper, we apply two theoretical turbulence models, DIA and the recent GISS model, to study properties of a turbulent channel flow. Both models provide a turbulent kinetic energy spectral function E(k) as the solution of a non-linear equation; the two models employ the same source function but different closures. The source function is characterized by a rate n sub s (k) which is derived from the complex eigenvalues of the Orr--Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. The O--S equation is solved for a variety of Reynolds numbers corresponding to available experimental data. A physical argument is presented whereby the central line velocity characterizing the basic flow, U0 sup L, is not to be identified with the U0 appearing in the experimental Reynolds number. The theoretical results are compared with two types of experimental data: (1) turbulence bulk properties, and (2) properties that depend stongly on the structure of the turbulence spectrun at low wave numbers. The only existing analytical expression for Pi (k) cannot be used in the present case because it applies to the case of a flat plate, not a finite channel.
The preferred mode of instability was investigated in an axisymmetric air jet of moderate Reynolds number. Natural instabilities are shown to scale with local shear-layer thickness and the preferred mode is shown to be a shear-layer instability. The spatial evolution of the preferred mode was examined by exciting the flow acoustically and then mapping the phase-locked velocity fluctuations. Throughout the potential core region the phase-locked profiles are shown to agree with the eigensolutions of the Orr-Sommerfeld stability equations provided the calculations are based on measured, mean velocity profiles. The excitation intensity was varied from low levels, where the flow was merely tagged, to high levels where the mean flow was substantially distorted, and over that range of excitation there was no apparent deterioration in the agreement with stability predictions.
The behavior of air pressure fluctuations induced by progressive water waves generated mechanically in a laboratory tank was simulated by solving a modified Orr-Sommerfeld equation in a transformed Eulerian wave-following frame of reference. Solution is obtained by modeling the mean and wave-coherent turbulent Reynolds stresses, the behavior of which in the turbulent boundary layer above the waves was simulated using a turbulent kinetic energy-dissipation model, properly modified to account for free-surface proximity and favorable pressure gradient effects. The distribution of both the wave-coherent turbulent Reynolds stress and pressure amplitudes and their corresponding phase lags was found to agree reasonably well with available laboratory data.
Numerical experiments were performed to clarify apparent differences between experimental observations and a theoretical prediction of the secondary instability in plane Poiseuille flow. It is shown that subharmonic breakdown is unlikely in natural transition as a result of the initial growth of what we call the 'minus' modes and consequent forcing of Orr-Sommerfeld modes present in the background noise. Subharmonic breakdown was achieved only when these minus modes were continuously suppressed.
Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes (NS) and transport equations. Here, practical aspects of these methods and their efficient implementation are examined, and several examples of flows in truly complex geometries are presented. The spectral element discretization for NS equations is introduced, and the convergence of the method is addressed. An efficient data management scheme is discussed in the context of parallel processing computations. The method is validated by comparing the spectral element solutions with the exact eigensolutions for the Orr-Sommerfeld equations in two and three dimensions. Computer-aided flow visualizations are presented for an impulsive flow past a sharp edge wedge. Three-dimensional states of channel flow disrupted by an array of cylindrical eddy promoters are studied, and the results of a direct simulation of the turbulent flow in a plane channel are presented.
In this paper, we apply two theoretical turbulence models, DIA and the recent GISS model, to study properties of a turbulent channel flow. Both models provide a turbulent kinetic energy spectral function E(k) as the solution of a non-linear equation; the two models employ the same source function but different closures. The source function is characterized by a rate n sub s (k) which is derived from the complex eigenvalues of the Orr-Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. The O-S equation is solved for a variety of Reynolds numbers corresponding to available experimental data. A physical argument is presented whereby the central line velocity characterizing the basic flow, U0 sup L, is not to be identified with the U0 appearing in the experimental Reynolds number. The theoretical results are compared with two types of experimental data: (1) turbulence bulk properties, and (2) properties that depend strongly on the structure of the turbulence spectrum at low wave numbers. The only existing analytical expression for Pi (k) cannot be used in the present case because it applies to the case of a flat plate, not a finite channel.
An analysis is made based upon the concept that the velocity fluctuations, and therefore, the Reynolds stresses, driven by the instability of the original flow grow until a new stable state is approached. The Reynolds stresses incorporated into the Orr-Sommerfeld equation are coupled with the main flow such that all the imaginary parts of the complex eigenvalues vanish, i.e., the original instability is eliminated. Using this stabilization principle, it is possible to find the Reynolds stresses as well as the mean velocity for plane Poiseuille flow with the Reynolds number slightly higher than the critical.
A very simple modification is presented for the Chebyshev tau method which can eliminate spurious eigenvalues, proceeding from a consideration of the vorticity-streamfunction reformulation of the Chebyshev tau method and the Chebyshev-Galerkin method, which have no spurious modes. Consideration of a model problem indicates that these two approaches are equivalent, and that they reduce to the present modification of the tau method. This modified tau method also eliminates spurious eigenvalues from the Orr-Sommerfeld equation.
The state of transition calculations in 1989 and desirable future trends are analyzed within the framework of aeronautical applications with emphasis on wall-bounded flows. A shift to studies of the early stages of transition, particularly to receptivity, from studies of the late stages in excessively confined periodic domains is proposed. It is pointed out that early-stage studies require a larger computational domain and more difficult boundary conditions but not a very fine grid. Intermediate strategies, more involved than the classical Orr-Sommerfeld equation but much more economical than full Navier-Stokes equations, are emphasized, along with the need for compressible calculations.
A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.
Pulsed field ionization (PFI) has been used in conjunction with a coherent VUV source to obtain high-resolution threshold photoelectron spectra for the (000), (010), (020), and (100) vibrational states of the N2O(+) cation. Simulations for the rotational profiles of each vibronic level were obtained by fitting the Buckingham-Orr-Sichel equations using accurate spectroscopic constants for the ground states of the neutral and the ion. The relative branch intensities are interpreted in terms of the partial waves of the outgoing photoelectron to which the ionic core is coupled and in terms of the angular momentum transferred to the core.
A model problem related to distributed receptivity to free-stream acoustic waves in laminar flow control (LFC) configurations is studied, within the Orr-Sommerfield framework, by a suitable extension of the Goldstein-Ruban theory for receptivity due to localized disturbances on the airfoil surface. The results, thus, complement the earlier work on the receptivity produced by local variations in the surface suction and/or surface admittance. In particular, we show that the cumulative effect of the distributed receptivity can be substantially larger than that of a single, isolated suction strip or slot. Furthermore, even if the receptivity is spread out over very large distances, the most effective contributions come from a relatively short region in vicinity of the lower branch of the neutral stability curve. The length scale of this region is intermediate to that of the mean of these two length scales. Finally, it is found that the receptivity is effectively dominated by a narrow band of Fourier components from the wall-suction and admittance distributions, roughly corresponding to a detuning of less than ten percent with respect to the neutral instability wavenumber at the frequency under consideration. The results suggest that the drop-off in receptivity magnitudes away from the resonant wavenumber is nearly independent of the frequency parameter.
The ultimate goal of this work is to study the flow inside a channel with riblets on one of the two walls. The method has been tested for two dimensional flows in the presence of bodies with a geometrical singularity and for three dimensional flows inside domains described by Cartesian coordinates. The results have been compared with previous numerical simulations and with experimental results. The cases considered are: (1) the growth of Orr-Sommerfield waves in plane Poiseulle flow; (2) the flow over a backward facing step; (3) the flow past a wedge; and (4) the flow inside a narrow channel. Finally, the case of a channel with two large riblets on a wall has been simulated. In this case a limited number of grid points is sufficient, and in spite of the slow convergence for the pressure solver, one is able to obtain solutions with a reasonable amount of computer time. At present solutions with very fine grids in all three directions can not be obtained, due to lack of a fast pressure solver for general curvilinear coordinates.
In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.
An analysis and calculations on the stability of the free laminar boundary layer between parallel streams were made for an incompressible fluid using the Tollmien-Schlichting theory of small disturbances. Because the boundary conditions are at infinity, two solutions of the Orr-Sommerfeld stability equations need not be considered, and the remaining two solutions are exponential in character at the infinite boundaries. The calculations show that the flow is unstable except for very low Reynolds numbers.
Linear receptivity studies in parallel flows are aimed at understanding how external forcing couples to the natural unstable motions which a flow can support. The vibrating ribbon problem models the original Schubauer and Skramstad boundary layer experiment and represents the classic boundary layer receptivity problem. The process by which disturbances are initiated in convectively-unstable jets and shear layers has also received attention. Gaster was the first to handle the boundary layer analysis with the recognition that spatial modes, rather than temporal modes, were relevant when studying convectively-unstable flows that are driven by a time-harmonic source. The amplitude of the least stable spatial mode, far downstream of the source, is related to the source strength by a coupling coefficient. The determination of this coefficient is at the heart of this type of linear receptivity study. The first objective of the present study was to determine whether the various wave number derivative factors, appearing in the coupling coefficients for linear receptivity problems, could be reexpressed in a simpler form involving adjoint eigensolutions. Secondly, it was hoped that the general nature of this simplification could be shown; indeed, a rather elegant characterization of the receptivity properties of spatial instabilities does emerge. The analysis is quite distinct from the usual Fourier-inversion procedures, although a detailed knowledge of the spectrum of the Orr-Sommerfeld equation is still required. Since the cylinder wake analysis proved very useful in addressing control considerations, the final objective was to provide a foundation upon which boundary layer control theory may be developed.