Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear Stability Theory”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The three-dimensional evolution of a plane mixing layer - Pairing and transition to turbulence

The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings has been simulated numerically. All simulations were begun from a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. Three-dimensional perturbations were used with amplitudes ranging from infinitesimal to large enough to trigger a rapid transition to turbulence. Pairing is found to inhibit the growth of infinitesimal three-dimensional disturbances, and to trigger the transition to turbulence in highly three-dimensional flows. The mechanisms responsible for the growth of three-dimensionality and onset of transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of spanwise vorticity, which undergo secondary rollups. The post-transitional simulated flow fields exhibit many properties characteristic of turbulent flows.

Moser, Robert D.↗
NASA NTRS1993

Spanwise scale selection in plane mixing layers

Direct numerical simulations of temporally evolving plane mixing layers undergoing as many as three pairings have been examined for evidence of spanwise scale change. All simulations were begun from a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. The amplitude of the initial three-dimensional disturbances varied from infinitesimal to large enough to trigger a rapid transition to turbulence. The time required for a change of characteristic spanwise scale with infinitesimal three-dimensional disturbances was found to be very long, requiring three or more pairings to complete a doubling of the spanwise scale. Stronger three-dimensionality can produce more rapid scale changes, but it is also likely to trigger transition to turbulence. No evidence was found for a change from an organized array of rib vortices at one spanwise scale to a similar array at a larger spanwise scale.

Rogers, Michael M.↗
NASA NTRS1993

Multigrid methods and high order finite difference for flow in transition - Effects of isolated and distributed roughness elements

The high order finite difference and multigrid methods have been successfully applied to direct numerical simulation (DNS) for flow transition in 3D channels and 3D boundary layers with 2D and 3D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semicoarsening multigrid method associated with line distributive relaxation scheme, and a new treatment of the outflow boundary condition, which needs only a very short buffer domain to damp all wave reflection, are developed. These approaches make the multigrid DNS code very accurate and efficient. This makes us not only able to do spatial DNS for the 3D channel and flat plate at low computational costs, but also able to do spatial DNS for transition in the 3D boundary layer with 3D single and multiple roughness elements. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments.

Liu, C.↗
NASA NTRS1992

Thermocapillary convection instability in microgravity crystal growth

Energy and linear stability theories are applied to study the stability properties of thermocapillary convection in a half zone model of the float zone crystal growth process. The analyses yield results which provide sufficient conditions for both stability and instability, respectively, and therefore, lower and upper bounds for the Marangoni numbers at which oscillatory thermocapillary convection will appear. The results of the analyses are compared with available experimental and numerical data.

Neitzel, G. P.↗
NASA NTRS1993

Transition prediction and control in subsonic flow over a hump

The influence of a surface roughness element in the form of a two-dimensional hump on the transition location in a two-dimensional subsonic flow with a free-stream Mach number up to 0.8 is evaluated. Linear stability theory, coupled with the N-factor transition criterion, is used in the evaluation. The mean flow over the hump is calculated by solving the interacting boundary-layer equations; the viscous-inviscid coupling is taken into consideration, and the flow is solved within the separation bubble. The effects of hump height, length, location, and shape; unit Reynolds number; free-stream Mach number, continuous suction level; location of a suction strip; continuous cooling level; and location of a heating strip on the transition location are evaluated. The N-factor criterion predictions agree well with the experimental correlation of Fage; in addition, the N-factor criterion is more general and powerful than experimental correlations. The theoretically predicted effects of the hump's parameters and flow conditions on transition location are consistent and in agreement with both wind-tunnel and flight observations.

Masad, Jamal A.↗
NASA NTRS1993

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗
NASA NTRS1993

Multigrid methods for flow transition in three-dimensional boundary layers with surface roughness

The efficient multilevel adaptive method has been successfully applied to perform direct numerical simulations (DNS) of flow transition in 3-D channels and 3-D boundary layers with 2-D and 3-D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semi-coarsening multigrid method associated with line distributive relaxation scheme, and an improved outflow boundary-condition treatment, which needs only a very short buffer domain to damp all order-one wave reflections, are developed. These approaches make the multigrid DNS code very accurate and efficient. This allows us not only to be able to do spatial DNS for the 3-D channel and flat plate at low computational costs, but also to do spatial DNS for transition in the 3-D boundary layer with 3-D single and multiple roughness elements, which would have extremely high computational costs with conventional methods. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments. The contribution of isolated and distributed roughness to transition is analyzed.

Liu, Chaoqun↗
NASA NTRS1993

Relationship between transition and modes of instability in supersonic boundary layers

The relationship between the predicted transition location and the first and second modes of instability in two-dimensional supersonic boundary-layer flow on a flat plate is examined. Linear stability theory and the N-factor criterion are used to predict transition location. The effect of heat transfer is also studied; the results demonstrate that the transition reversal phenomenon can be explained by the opposite effect of cooling on the first and second modes of instability. Compressibility of destabilizing at free-stream Mach numbers of 2 to 3.5. The predicted transition location is due to the oblique first modes of instability, up to free-stream Mach numbers between 6 and 6.5. At higher Mach numbers, the predicted transition location is due to a combination of two-dimensional first and second modes of instability.

Masad, Jamal A.↗
NASA NTRS1993

Computational fluid mechanics

Two papers are included in this progress report. In the first, the compressible Navier-Stokes equations have been used to compute leading edge receptivity of boundary layers over parabolic cylinders. Natural receptivity at the leading edge was simulated and Tollmien-Schlichting waves were observed to develop in response to an acoustic disturbance, applied through the farfield boundary conditions. To facilitate comparison with previous work, all computations were carried out at a free stream Mach number of 0.3. The spatial and temporal behavior of the flowfields are calculated through the use of finite volume algorithms and Runge-Kutta integration. The results are dominated by strong decay of the Tollmien-Schlichting wave due to the presence of the mean flow favorable pressure gradient. The effects of numerical dissipation, forcing frequency, and nose radius are studied. The Strouhal number is shown to have the greatest effect on the unsteady results. In the second paper, a transition model for low-speed flows, previously developed by Young et al., which incorporates first-mode (Tollmien-Schlichting) disturbance information from linear stability theory has been extended to high-speed flow by incorporating the effects of second mode disturbances. The transition model is incorporated into a Reynolds-averaged Navier-Stokes solver with a one-equation turbulence model. Results using a variable turbulent Prandtl number approach demonstrate that the current model accurately reproduces available experimental data for first and second-mode dominated transitional flows. The performance of the present model shows significant improvement over previous transition modeling attempts.

Hassan, H. A.↗
NASA NTRS1994

Numerical simulation of two-dimensional spatially-developing mixing layers

Two-dimensional, incompressible, spatially developing mixing layer simulations are performed at Re = 10(exp 2) and 10(exp 4) with two classes of perturbations applied at the inlet boundary; combinations of discrete modes from linear stability theory, and a broad spectrum of modes derived from experimentally measured velocity spectra. The effect of the type and strength of inlet perturbations on vortex dynamics and time-averaged properties are explored. Two-point spatial velocity and autocorrelations are used to estimate the size and lifetime of the resulting coherent structures and to explore possible feedback effects. The computed time-averaged properties such as mean velocity profiles, turbulent statistics, and spread rates show good agreement with experimentally measured values. It is shown that by forcing with a broad spectrum of modes derived from an experimental energy spectrum many experimentally observed phenomena can be reproduced by a 2-D simulation. The strength of the forcing merely affected the length required for the dominant coherent structures to become fully-developed. Thus intensities comparable to those of the background turbulence in many wind tunnel experiments produced the same results, given sufficient simulation length.

Wilson, R. V.↗
NASA NTRS1994

A spectrally accurate boundary-layer code for infinite swept wings

This report documents the development, validation, and application of a spectrally accurate boundary-layer code, WINGBL2, which has been designed specifically for use in stability analyses of swept-wing configurations. Currently, we consider only the quasi-three-dimensional case of an infinitely long wing of constant cross section. The effects of streamwise curvature, streamwise pressure gradient, and wall suction and/or blowing are taken into account in the governing equations and boundary conditions. The boundary-layer equations are formulated both for the attachment-line flow and for the evolving boundary layer. The boundary-layer equations are solved by marching in the direction perpendicular to the leading edge, for which high-order (up to fifth) backward differencing techniques are used. In the wall-normal direction, a spectral collocation method, based upon Chebyshev polynomial approximations, is exploited. The accuracy, efficiency, and user-friendliness of WINGBL2 make it well suited for applications to linear stability theory, parabolized stability equation methodology, direct numerical simulation, and large-eddy simulation. The method is validated against existing schemes for three test cases, including incompressible swept Hiemenz flow and Mach 2.4 flow over an airfoil swept at 70 deg to the free stream.

Pruett, C. David↗
NASA NTRS1994

Effect of crossflow on Goertler instability in incompressible boundary layers

Linear stability theory is used to study the effect of crossflow on Goertler instability in incompressible boundary layers. The results cover a wide range of sweep angle, pressure gradient, and wall curvature parameters. It is shown that the crossflow stabilizes Goertler disturbances by reducing the maximum growth rate and shrinking the unstable band of spanwise wave numbers. On the other hand, the effect of concave wall curvature on crossflow instability is destabilizing. Calculations show that the changeover from Goertler to crossflow instabilities is a function of Goertler number, pressure gradient, and sweep angle. The results demonstrate that Goertler instability may still be relevant in the transition process on swept wings even at large angles of sweep if the pressure gradient is sufficiently small. The influence of pressure gradient and sweep can be combined by defining a crossflow Reynolds number. Thus, the changeover from Goertler to crossflow instability takes place at some critical crossflow Reynolds number whose value increases with Goertler number.

Zurigat, Y. H.↗
NASA NTRS1995

Three-dimensional Structure in a Supersonic Jet: Behavior near Centerline

The unsteady structure of a supersonic jet is highly three dimensional, though the mean flow is axisymmetric. In simulating a circular jet, the centerline represents a computational boundary. As such, spurious modes can be generated near centerline, unless special attention is given to the behavior of the 3D structure near the centerline. Improper treatment of the dependent variables near the centerline results in the solution diverging or being suitable only for small amplitude excitation. With a careful treatment of the centerline formation, no spurious mode is generated. The results show that a near-linear disturbance growth is obtained, as the linear stability theory indicates. At high levels of excitation, nonlinear development of disturbances is evident and saturation is reached downstream.

Shih, S. H.↗
NASA NTRS1993

Numerical simulation of stability and stability control of high speed compressible rotating couette flow

The nonlinear temporal evolution of disturbances in compressible flow between infinitely long, concentric cylinders is investigated through direct numerical simulations of the full, three-dimensional Navier-Stokes and energy equations. Counter-rotating cylinders separated by wide gaps are considered with supersonic velocities of the inner cylinder. Initially, the primary disturbance grows exponentially in accordance with linear stability theory. As the disturbances evolve, higher harmonics and subharmonics are generated in a cascading order eventually reaching a saturation state. Subsequent highly nonlinear stages of the evolution are governed by the interaction of the disturbance modes, particularly the axial subharmonics. Nonlinear evolution of the disturbance field is characterized by the formation of high-shear layers extending from the inner cylinder towards the center of the gap in the form of jets similar to the ejection events in transitional and turbulent wall-bounded shear flows.

Biringen, Sedat↗
NASA NTRS1994

Transition correlations in three-dimensional boundary layers

The stability and transition characteristics of three-dimensional boundary-layer flows are examined. First, the flow over a rotating cone is considered computationally. An increase of stagnation temperature is found to be only slightly stabilizing. Parameter studies on the simple rotating-cone geometry provide a large database of three-dimensional boundary-layer profiles and associated stability characteristics. To determine the possibility of correlating transition location with parameters based purely on basic-state three-dimensional boundary-layer profile characteristics, an empirical transition location of N = 9 is assumed. Transition location does not correlate with the traditional crossflow Reynolds number. A more appropriate definition for crossflow Reynolds number is found and termed R(sub cf(new)). This new parameter appears to correlate for transition location when plotted against maximum crossflow velocity. Then, the flow over a yawed cone is considered experimentally. The correlation results obtained from the rotating-cone work are applied to the actual measured transition locations on two different yawed-cone models under various angle-of-attack conditions in two different experimental facilities and are verified. This correlation is only suggested as a tool for preliminary transition prediction and design in three-dimensional boundary layers; once a preliminary shape is selected, further linear stability theory or parabolized stability equation calculations are strongly urged.

Reed, Helen L.↗
NASA NTRS1996

Transition Studies on a Swept-Wing Model

The present investigation contributes to the understanding of boundary-layer stability and transition by providing detailed measurements of carefully-produced stationary crossflow vortices. It is clear that a successful prediction of transition in swept-wing flows must include an understanding of the detailed physics involved. Receptivity and nonlinear effects must not be ignored. Linear stability theory correctly predicts the expected wavelengths and mode shapes for stationary crossflow, but fails to predict the growth rates, even for low amplitudes. As new computational and analytical methods are developed to deal with three-dimensional boundary layers, the data provided by this experiment will serve as a useful benchmark for comparison.

Saric, William S.↗
NASA NTRS1997

Distributed-Roughness Effects on Stability and Transition In Swept-Wing Boundary Layers

Boundary-layer stability experiments are conducted in the Arizona State University Unsteady Wind Tunnel on a 45 deg swept airfoil. The pressure distribution and test conditions are designed to suppress Tollmien-Schlichting disturbances and provide crossflow-dominated transition. The surface of the airfoil is finely polished to a near mirror finish. Under these conditions, submicron surface irregularities cause the naturally occurring stationary crossflow waves to grow to nonuniform amplitudes. Spanwise-uniform stationary crossflow disturbances are generated through careful control of the initial conditions with full-span arrays of micron-high roughness elements near the attachment line. Detailed hot-wire measurements are taken to document the stationary crossflow structure and determine growth rates for the total and individual-mode disturbances. Naphthalene flow visualization provides transition location information. Roughness spacing and roughness height are varied to examine the effects on transition location and all amplified wavelengths. The measurements show that roughness spacings that do not contain harmonics equal to the most unstable wavelength as computed by linear stability theory effectively suppress the most unstable mode. Under certain conditions, subcritical roughness spacing delays transition past that of the corresponding smooth surface.

Carrillo, Ruben B., Jr.↗
NASA NTRS1997

A Review of Hypersonic Boundary Layer Stability Experiments in a Quiet Mach 6 Wind Tunnel

Three recent experimental studies of transition on cones with adverse pressure gradient produced by a flared afterbody and with the additive stability modifiers of wall cooling, angle of attack and bluntness are reviewed. All tests were conducted in a quiet Mach 6 wind tunnel. The dominant instability was found to be the second mode. For the cases examined with linear stability theory, the N factors at mode saturation were in the range of 8.5 to 11. Evidence of a combined second-mode/Gortler transition process was found. Mean, rms and spectral freestream data for the quiet facility is presented and the role of low frequency freestream noise is discussed.

Wilkinson, Stephen P.↗