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 19 records

Transition prediction in external flows via linear stability theory

Linear stability theory results are presented for flight and low-disturbance wind tunnels, at low and high speeds, and for T-S, Gortler, and cross-flow modes. Results show that, in the absence of Morkovin 'bypasses', transition at low background disturbance levels corresponds to N-factors of the order of 9 to 11, provided that phenomena such as curvature effects are included. The results point to a wider range of applicability for the e exp N method than previously conjectured, and suggest that the e exp N method can be used to parametrize transition predictions as a function of parameters which affect the mean flow profiles for such purposes as LFC and minimum drag design.

Bushnell, D. M.↗

Application of linear stability theory in laminar flow design

The linear stability of fully three-dimensional supersonic boundary layers formed over swept-wing configurations is investigated using a modified version of the linear stability code COSAL. Configurations studied include a highly swept leading-edge model to be utilized for transition studies in the LARC Low-disturbance Mach 3.5 Pilot Tunnel. The model is a representation of the leading edge of a laminar flow control wing for the F-16XL aircraft. In addition, the region over a laminar flow control glove fitted on the midportion of an F-16XL wing was studied. For each configuration, estimates of the location of the onset of transition were computed using linear stability theory and the e exp N method. The effectiveness of suction in stabilizing the boundary layer over the F-16XL wing glove was also investigated.

Iyer, Venkit↗

On the application of linear stability theory to the problem of supersonic boundary-layer transition

Linear stability theory is used to calculate the amplitude ratio of constant-frequency disturbances as a function of Reynolds number for insulated and cooled-wall flat-plate boundary layers between Mach numbers 1.3 and 5.8. The growth curves are used to examine the consequences of using a fixed amplitude of the most unstable frequency as a transition criterion. The effect of free-stream Mach number on insulated-wall boundary layers is calculated, assuming that the initial disturbance level is constant, is proportional to the square of the free-stream Mach number and to the square root of the energy density of the one-dimensional power spectra of free-stream disturbances measured in supersonic wind tunnels.

Mack, L. M.↗

Boundary-layer linear stability theory

Most fluid flows are turbulent rather than laminar and the reason for this was studied. One of the earliest explanations was that laminar flow is unstable, and the linear instability theory was first developed to explore this possibility. A series of early papers by Rayleigh produced many notable results concerning the instability of inviscid flows, such as the discovery of inflectional instability. Viscosity was commonly thought to act only to stabilize the flow, and flows with convex velocity profiles appeared to be stable. The investigations that led to a viscous theory of boundary layer instability was reported. The earliest application of linear stability theory to transition prediction calculated the amplitude ratio of the most amplified frequency as a function of Reynolds number for a Blasius boundary layer, and found that this quantity had values between five and nine at the observed Ret. The experiment of Schubauer and Skramstad (1947) completely reversed the prevailing option and fully vindicated the Gottingen proponents of the theory. This experiment demonstrated the existence of instability waves in a boundary layer, their connection with transition, and the quantitative description of their behavior by the theory of Tollmien and Schlichting. It is generally accepted that flow parameters such as pressure gradient, suction and heat transfer qualitatively affect transition in the manner predicted by the linear theory, and in particular that a flow predicted to be stable by the theory should remain laminar. The linear theory, in the form of the e9, or N-factor is today in routine use in engineering studies of laminar flow. The stability theory to boundary layers with pressure gradients and suction was applied. The only large body of numerical results for exact boundary layer solutions before the advent of the computer age by calculating the stability characteristics of the Falkner-Skan family of velocity profiles are given. When the digital computer reached a stage of development which permit the direct solution of the primary differential equations, numerical results were obtained from the linear theory during the next 10 years for many different boundary layer flows: three dimensional boundary layers; free convention boundary layers; compressible boundary layers; boundary layers on compliant walls; a recomputation of Falkner-Skan flows; unsteady boundary layers; and heated wall boundary layers.

Leslie M Mack↗

Linear stability theory and three-dimensional boundary layer transition

The viewgraphs and discussion of linear stability theory and three dimensional boundary layer transition are provided. The ability to predict, using analytical tools, the location of boundary layer transition over aircraft-type configurations is of great importance to designers interested in laminar flow control (LFC). The e(sup N) method has proven to be fairly effective in predicting, in a consistent manner, the location of the onset of transition for simple geometries in low disturbance environments. This method provides a correlation between the most amplified single normal mode and the experimental location of the onset of transition. Studies indicate that values of N between 8 and 10 correlate well with the onset of transition. For most previous calculations, the mean flows were restricted to two-dimensional or axisymmetric cases, or have employed simple three-dimensional mean flows (e.g., rotating disk, infinite swept wing, or tapered swept wing with straight isobars). Unfortunately, for flows over general wing configurations, and for nearly all flows over fuselage-type bodies at incidence, the analysis of fully three-dimensional flow fields is required. Results obtained for the linear stability of fully three-dimensional boundary layers formed over both wing and fuselage-type geometries, and for both high and low speed flows are discussed. When possible, transition estimates form the e(sup N) method are compared to experimentally determined locations. The stability calculations are made using a modified version of the linear stability code COSAL. Mean flows were computed using both Navier Stokes and boundary-layer codes.

Spall, Robert E.↗

Linear-stability theory of thermocapillary convection in a model of float-zone crystal growth

Linear-stability theory has been applied to a basic state of thermocapillary convection in a model half-zone to determine values of the Marangoni number above which instability is guaranteed. The basic state must be determined numerically since the half-zone is of finite, O(1) aspect ratio with two-dimensional flow and temperature fields. This, in turn, means that the governing equations for disturbance quantities will remain partial differential equations. The disturbance equations are treated by a staggered-grid discretization scheme. Results are presented for a variety of parameters of interest in the problem, including both terrestrial and microgravity cases.

Neitzel, G. P.↗

Eigenvalue Sensitivity Computations for Linear Stability Theory

To realize the drag reduction benefit of boundary-layer transition control strategies, it is crucial to integrate transition prediction into the vehicle design through an optimization process. The integration of transition prediction based on linear stability analysis into adjoint d design optimization requires coupling an adjoint enabled computational fluid dynamics (CFD) solver with an adjoint enabled linear stability code. In particular, the boundary-layer transition location is often predicted using the N-factor method based on linear stability theory (LST). Thus, sensitivity of the linear-stability eigenvalues constitute an essential building block for optimizing the laminar flow performance. The present paper describes an implementation of LST eigenvalue sensitivity analysis that can be easily coupled with a CFD solver. Specifically, we describe a discrete adjoint formulation for the transition location prediction based on the N-factor method. The verification of this formulation is carried out by comparing the adjoint-based sensitivity of the local growth rate of a given instability mode with respect to the disturbance frequency, and the adjoint-based sensitivity of the transition location with respect to spanwise wavenumber with those sensitivities computed using a finite-difference approximation. Finally, the adjoint LST formulation is applied to flat-plate boundary-layer flows at transonic, supersonic, and hypersonic conditions, to determine the behavior and sensitivities of the transition location with respect to a range of disturbance spanwise wavenumbers.

Boundary Layer Transition↗

Morphological instability in epitaxially strained dislocation-free solid films - Linear stability theory

The morphological instability of a growing epitaxially strained dislocation-free solid film is analyzed. An evolution equation for the film surface is derived in the dilute limit of vacancies based on surface diffusion driven by a stress-dependent chemical potential. From the time-dependent linear stability problem the conditions for which a growing film is unstable are determined. It is found that the instability is driven by the lattice mismatch between the film and the substrate; however, low temperatures as well as elastically stiff substrates are stabilizing influences. The results also reveal that the critical film thickness for instability depends on the growth rate of the film itself. Detailed comparison with experimental observations indicates that the instability described exhibits many of the observed features of the onset of the 'island instability'.

Spencer, B. J.↗

Some comparisons of linear stability theory with experiment at supersonic and hypersonic speed

Measurements of waves initiated by wind-tunnel sounds are obtained for both supersonic and hypersonic speed. Flat-plate experiments carried out at Mach numbers of 4.5, 3.7, and 2.4 under the condition of laminar tunnel-wall boundary layers whereby a low background fluctuation level prevailed are outlined, as well as flat-plate experiments with sound-induced waves at supersonic speeds. Experimental studies of 'natural' wave growth on sharp cones at zero angle of attack and Mach numbers near 8 are discussed. Attention is focused on experimental research needed on boundary-layer receptivity to eddy Mach-wave radiation, hypersonic speed, and the stability-transition relationship in zero-pressure-gradient layers.

Kendall, James M., Jr.↗

Convolutional Neural Network for Transition Modeling Based on Linear Stability Theory

Transition prediction is an important aspect of aerodynamic design because of its impact on skin friction and potential coupling with flow separation characteristics. Traditionally, the modeling of transition has relied on correlation-based empirical formulas based on integral quantities such as the shape factor of the boundary layer. However, in many applications of computational fluid dynamics, the shape factor is not straightforwardly available or not well-defined. We propose using the complete velocity profile along with other quantities (e.g., frequency, Reynolds number) to predict the perturbation amplification factor. While this can be achieved with regression models based on a classical fully connected neural network, such a model can be computationally more demanding. We propose a novel convolutional neural network inspired by the underlying physics as described by the stability equations. Specifically, convolutional layers are first used to extract integral quantities from the velocity profiles, and then fully connected layers are used to map the extracted integral quantities, along with frequency and Reynolds number, to the output (amplification ratio). Numerical tests on classical boundary layers clearly demonstrate the merits of the proposed method. More importantly, we demonstrate that, for Tollmien-Schlichting instabilities in two-dimensional, low-speed boundary layers, the proposed network encodes information in the boundary layer profiles into an integral quantity that is strongly correlated to a well-known, physically defined parameter – the shape factor.

Laminar-turbulent transition↗

Review of linear compressible stability theory

Aspects of the linear compressible stability theory are considered, with special attention given to the inviscid theory and the additional solutions that arise when there is a region of supersonic flow relative to the phase velocity. For the case of highly cooled flat-plate boundary layers at Mach number 5.8, the unstable region is found to include supersonic outgoing waves. A previously unknown neutral incoming wave has also been identified. An example of viscous multiple solutions is provided, and calculations of higher viscous discrete modes and the compressible counterpart of the Squire mode are presented. It is noted that compressible stability theory is missing a firm connection with boundary-layer transition.

Mack, Leslie M.↗

On the role of secondary instabilities in laminar-turbulent transition of 2D and 3D boundary layers

The accurate prediction of laminar/turbulent transition is one of the fundamental problems in engineering fluid mechanics. There is almost unanimous consent that such a transition criterion should come from stability theory. Linear primary stability theory describes the initial stage of transition, but falls short of predicting transition. Only in conjunction with empirical correlations, the widely used e(sup n) method is obtained, which, however, lacks a solid physical base. Three-dimensional secondary instabilities are known to play an important role in the transition process. However, no use has been made so far of secondary instabilities, instability interactions or wave resonances to define a 'transition location'. The paper summarizes new attempts to identify certain interaction and resonance phenomena within the laminar-turbulent transition regime in two and three-dimensional boundary layers which are associated with rapid structural and temporal changes of fluctuations beyond their exponential growths.

Dallman, Uwe↗

Langley Stability and Transition Analysis Code (LASTRAC) Version 1.2 User Manual

LASTRAC is a general-purposed, physics-based transition prediction code released by NASA for Laminar Flow Control studies and transition research. The design and development of the LASTRAC code is aimed at providing an engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. It was written from scratch based on the state-of-the-art numerical methods for stability analysis and modern software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory or linear parabolized stability equations method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. This document describes the governing equations, numerical methods, code development, detailed description of input/output parameters, and case studies for the current release of LASTRAC.

Chang, Chau-Lyan↗

Hypersonic Boundary Layer Stability Experiments in a Quiet Wind Tunnel with Bluntness Effects

Hypersonic boundary layer measurements over a flared cone were conducted in a Mach 6 quiet wind tunnel at a freestream unit Reynolds number of 2.82 million/ft. This Reynolds number provided laminar-to-transitional flow over the cone model in a low-disturbance environment. Four interchangeable nose-tips, including a sharp-tip, were tested. Point measurements with a single hot-wire using a novel constant voltage anemometer were used to measure the boundary layer disturbances. Surface temperature and schlieren measurements were also conducted to characterize the transitional state of the boundary layer and to identify instability modes. Results suggest that second mode disturbances were the most unstable and scaled with the boundary layer thickness. The second mode integrated growth rates compared well with linear stability theory in the linear stability regime. The second mode is responsible for transition onset despite the existence of a second mode subharmonic. The subharmonic disturbance wavelength also scales with the boundary layer thickness. Furthermore, the existence of higher harmonics of the fundamental suggests that nonlinear disturbances are not associated with 'high' free stream disturbance levels. Nose-tip radii greater than 2.7% of the base radius completely stabilized the second mode.

Lachowicz, Jason T.↗

Hypersonic Boundary Layer Stability over a Flared Cone in a Quiet Tunnel

Hypersonic boundary layer measurements were conducted over a flared cone in a quiet wind tunnel. The flared cone was tested at a freestream unit Reynolds number of 2.82x106/ft in a Mach 6 flow. This Reynolds number provided laminar-to-transitional flow over the model in a low-disturbance environment. Point measurements with a single hot wire using a novel constant voltage anemometry system were used to measure the boundary layer disturbances. Surface temperature and schlieren measurements were also conducted to characterize the laminar-to-transitional state of the boundary layer and to identify instability modes. Results suggest that the second mode disturbances were the most unstable and scaled with the boundary layer thickness. The integrated growth rates of the second mode compared well with linear stability theory in the linear stability regime. The second mode is responsible for transition onset despite the existence of a second mode sub-harmonic. The sub-harmonic wavelength also scales with the boundary layer thickness. Furthermore, the existence of higher harmonics of the fundamental suggests that non-linear disturbances are not associated with high free stream disturbance levels.

Lachowicz, Jason T.↗

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗

Linear viscous stability theory for stably stratified shear flow - A review.

In this paper the linear viscous stability theory for stably stratified parallel shear flow is reviewed and some new results are presented. Attention is focused on recent work on unbounded flows with emphasis placed on results which demonstrate apparent destabilizing effects of diffusivity which lie beyond the scope of inviscid theory. In particular it is shown that self-excited disturbances may exist with phase speeds lying outside the range of basic flow speed and that with strong thermal diffusivity, instability may occur even though the local Richardson number exceeds 0.25 throughout the flow.

Gage, K. S.↗

Boundary-layer receptivity to oblique freestream vorticity waves for a high-enthalpy hypersonic flow

The receptivity of a Mach 15 straight-cone boundary layer to oblique freestream vorticity waves is investigated using direct numerical simulation (DNS) alongside linear stability theory and the linear parabolized stability equations. A thermochemical nonequilibrium gas model is used. Oblique freestream vorticity waves at frequencies of 400, 800, and 1200 kHz are considered, with incident angles ranging from 0° to 29.4° at 400 kHz, 0° to 15.7° at 800 kHz, and 0° to 20.6° at 1200 kHz. The 400 kHz case is of primary interest due to the strong second-mode amplification at this frequency. Although the underlying base flow is axisymmetric, the oblique vorticity waves lead to a fully three-dimensional boundary-layer disturbance whose characteristics vary depending on the azimuthal ray relative to the freestream wave. Moving downstream within the second-mode instability region, some clear trends emerge in terms of the boundary-layer disturbance amplitudes; that is, disturbance amplitudes are highest at the leeward ray (relative to the freestream wave), but weakest about halfway between the windward and leeward rays. Increasing the incident angle causes the amplitudes to increase on the leeward ray and decrease on the windward ray. Moreover, the boundary-layer disturbance throughout contains a wide spectrum of azimuthal wavenumbers in which the disturbance energy falls off at higher wavenumbers. Increasing the incident angle causes the azimuthal spectrum of the boundary-layer disturbance to broaden overall. Qualitatively similar results are found for the two higher frequencies leading up to the peak-amplitude locations corresponding to the second-mode instability.

42 ENGINEERING↗