Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 361 records · Page 20

Spatiotemporal chaos near the onset of cellular growth during thin-film solidification of a binary alloy

Thin-film solidification experiments with a succinonitrile-acetone alloy are used to observe the long time-scale dynamics of cellular crystal growth at growth rates only slightly above the critical value VC = Vc(lambda sub c) for the onset of morphological instability. Under these conditions only very small amplitude cells are observed with wavelengths near the value predicted by linear stability theory lambda = lambda sub c. At long times, microstructures with wavelengths significantly finer than lambda suc c form by nucleation at defects across the interface. These interfaces do not have a unique microstructure, but seem to exhibit spatiotemporal chaos on a long time scale caused by the continual birth and death of cells by tip splitting and cell annihilation in grooves.

Lee, J. T. C.↗

Crossflow Stability and Transition Experiments in a Swept-Wing Flow

An experimental examination of crossflow instability and transition on a 45 degree swept wing is conducted in the Arizona State University Unsteady Wind Tunnel. The stationary-vortex pattern and transition location are visualized using both sublimating-chemical and liquid-crystal coatings. Extensive hot-wire measurements are conducted at several measurement stations across a single vortex track. The mean and travelling-wave disturbances are measured simultaneously. Stationary-crossflow disturbance profiles are determined by subtracting either a reference or a span-averaged velocity profile from the mean-velocity data. Mean, stationary-crossflow, and travelling-wave velocity data are presented as local boundary-layer profiles and as contour plots across a single stationary-crossflow vortex track. Disturbance-mode profiles and growth rates are determined. The experimental data are compared to predictions from linear stability theory.

Dagenhart, John Ray↗

The three-dimensional evolution of a plane mixing layer. Part 2: Pairing and transition to turbulence

The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings was simulated numerically. Initial conditions for all simulations consisted of 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 both 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 as well as the initial phases of the transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of span wise vorticity, which undergo a secondary roll up. Transition also produces an increase in the degree of scalar mixing, in agreement with experimental observations of mixing transition. Simulations were also conducted to investigate changes in span wise length scale that may occur in response to the change in stream wise length scale during a pairing. The linear mechanism for this process was found to be very slow, requiring roughly three pairings to complete a doubling of the span wise scale. Stronger three-dimensionality can produce more rapid scale changes but is also likely to trigger transition to turbulence. No evidence was found for a change from an organized array of rib vortices at one span wise scale to a similar array at a larger span wise scale.

Moser, Robert D.↗

Direct numerical simulation of laminar breakdown in high-speed, axisymmetric boundary layers

The compressible Navier-Stokes equations are solved using spectral collocation and high-order compact-difference techniques to simulate the laminar breakdown in high-speed, axisymmetric boundary-layer flow. Mach 4.5 flow along a hollow cylinder and Mach 6.8 flow along a sharp cone are considered. Data obtained replicate two previously unexplained phenomena, namely, the appearance of so-called 'rope-like waves' and 'the precursor transition effect', in which transitional flow originates near the critical layer well upstream of the transition location at the wall. The numerical data also reveal that neither of these effects can be explained, even qualitatively, by linear stability theory alone. It is shown that rope-like appearance arises from secondary instability. Certain features of the precursor transition effect also emerge from secondary instability but its nature is revealed to be fundamentally nonlinear.

Pruett, C. D.↗

A high-frequency, secondary instability of crossflow vortices that leads to transition

Three-dimensional boundary-layer transition experiments are currently being conducted on a 45 deg swept wing in the Arizona State University Unsteady Wind Tunnel. Crossflow-dominated transition is produced via a model with contoured end liners to simulate infinite swept-wing flow. Fixed-wavelength, stationary and travelling crossflow vortices are observed. The frequencies of the most amplified travelling waves are in agreement with linear-stability theory; however travelling waves at frequencies an order of magnitude higher than predicted are also observed near transition. Near the transition location, the distorted boundary layer is due to the stationary crossflow vortex and is subject to a Rayleigh-type instability in the stream direction. As a result, a high-frequency secondary instability is detected in the transition region and spatial relations of the process are well documented by coupled use of flow visualization and hot-wire measurements.

Kohama, Yasuaki↗

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.↗

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.↗

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.↗

The effect of grid aspect ratio on convergence

Vector stability theory and computational experiments are used to examine the effects of high aspect ratio grids on algorithm convergence. Emphasis is placed on identifying the degree of convergence deterioration that ADI schemes experience as grid aspect ratio increases and on deriving methods for circumventing this problem. Approximately factored implicit schemes are found to experience convergence deterioration because of nonoptimum local time stepping procedures and the increased need for viscous preconditioning. An enhanced algorithm is devised on the basis of improved selection of the local time step, appropriate definition of viscous preconditioning matrix, and proper implementation of boundary conditions. The novel algorithm provides uniformly efficient convergence at all aspect ratios for both Euler and Navier-Stokes computations for a variety of test problems.

Buelow, P. E. O.↗

Stability analysis and control of tethered satellites

Issues of stability and stabilization of the Shuttle/Tethered Satellite System (TSS) are considered. The basic, nonlinear equations of motion of the TSS are derived using the system Lagrangian. A tension control strategy in feedback form is introduced to stabilize the TSS during station-keeping. In addition, the stability properties of a constant angle control method for deployment and retrieval are investigated using Liapunov stability theory.

Liaw, D.-C.↗

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.↗

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↗

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↗

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.↗

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.↗

Code validation for the simulation of supersonic viscous flow about the F-16XL

The viewgraphs and discussion on code validation for the simulation of supersonic viscous flow about the F-16XL are provided. Because of the large potential gains related to laminar flow on the swept wings of supersonic aircraft, interest in the applications of laminar flow control (LFC) techniques in the supersonic regime has increased. A supersonic laminar flow control (SLFC) technology program is currently underway within NASA. The objective of this program is to develop the data base and design methods that are critical to the development of laminar flow control technology for application to supersonic transport aircraft design. Towards this end, the program integrates computational investigations underway at NASA Ames-Moffett and NASA Langley with flight-test investigations being conducted on the F-16XL at the NASA Ames-Dryden Research Facility in cooperation with Rockwell International. The computational goal at NASA Ames-Moffett is to integrate a thin-layer Reynolds averaged Navier-Stokes flow solver with a stability analysis code. The flow solver would provide boundary layer profiles to the stability analysis code which in turn would predict transition on the F-16XL wing. To utilize the stability analysis codes, reliable boundary layer data is necessary at off-design cases. Previously, much of the prediction of boundary layer transition has been accomplished through the coupling of boundary layer codes with stability theory. However, boundary layer codes may have difficulties at high Reynolds numbers, of the order of 100 million, and with the current complex geometry in question. Therefore, a reliable code which solves the thin-layer Reynolds averaged Navier-Stokes equations is needed. Two objectives are discussed, the first in greater depth. The first objective is method verification, via comparisons of computations with experiment, of the reliability and robustness of the code. To successfully implement LFC techniques to the F-16XL wing, the flow about the leading edge must be maintained as laminar flow. Therefore, the second objective is to focus on a series of numerical simulations with different values of angle of attack, alpha, and Reynolds numbers. The purpose of the simulations is to study their effects on the two main factors which precipitate transition to turbulence at leading edges of highly swept wings (e.g., 'spanwise contamination' and 'crossflow instability').

Flores, Jolen↗

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.↗

Towards enhancing and delaying disturbances in free shear flows

The family of shear flows comprising the jet, wake, and the mixing layer are subjected to perturbations in an inviscid incompressible fluid. By modeling the basic mean flows as parallel with piecewise linear variations for the velocities, complete and general solutions to the linearized equations of motion can be obtained in closed form as functions of all space variables and time when posed as an initial value problem. The results show that there is a continuous as well as the discrete spectrum that is more familiar in stability theory and therefore there can be both algebraic and exponential growth of disturbances in time. These bases make it feasible to consider control of such flows. To this end, the possibility of enhancing the disturbances in the mixing layer and delaying the onset in the jet and wake is investigated. It is found that growth of perturbations can be delayed to a considerable degree for the jet and the wake but, by comparison, cannot be enhanced in the mixing layer. By using moving coordinates, a method for demonstrating the predominant early and long time behavior of disturbances in these flows is given for continuous velocity profiles. It is shown that the early time transients are always algebraic whereas the asymptotic limit is that of an exponential normal mode. Numerical treatment of the new governing equations confirm the conclusions reached by use of the piecewise linear basic models. Although not pursued here, feedback mechanisms designed for control of the flow could be devised using the results of this work.

Criminale, W. O.↗