Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear stability analysis”

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 55 records · Page 3

Gravitational Effects on Flow Instability and Transition in Low Density Jets

Recent experiments have shown that low-density gas jets injected into a high-density ambient gas undergo an instability mode, leading to highly-periodic oscillations in the flow-field for certain conditions. The transition from laminar to turbulent flow in these jets is abrupt, without the gradual change in scales. Even the fine scale turbulent structure repeats itself with extreme regularity from cycle to cycle. Similar observations were obtained in buoyancy-dominated and momentum-dominated jets characterized by the Richardson numbers, Ri = [gD(rho(sub a)-rho(sub j))/rho(sub j)U(sub j)(exp 2) ] where g is the gravitational acceleration, D is the jet diameter, rho(sub a) and rho(sub a) are, respectively, the free-stream and jet densities, and U(sub j) is the mean jet exit velocity. At high Richardson numbers, the instability is presumably caused by buoyancy since the flow-oscillation frequency (f) or the Strouhal number, St = [fD/U(sub j)] scales with Ri. In momentum-dominated jets, however, the Strouhal number of the oscillating flow is relatively independent of the Ri. In this case, a local absolute instability is predicted in the potential core of low-density jets with S [= rho(sub j)/rho(sub a)] < 0.7, which agrees qualitatively with experiments. Although the instability in gas jets of high Richardson numbers is attributed to buoyancy, direct physical evidence has not been acquired in experiments. If the instability is indeed caused by buoyancy, the near-field flow structure of the jet will change significantly when the buoyancy is removed, for example, in the microgravity environment. Thus, quantitative data on the spatial and temporal evolutions of the instability, length and time scale of the oscillating mode and its effects on the mean flow and breakdown of the potential core are needed in normal and microgravity to delineate gravitational effects in buoyant jets. In momentum dominated low-density jets, the instability is speculated to originate in the potential core. However, experiments have not succeeded in identifying the direct physical cause of the instability. For example, the theory predicts an oscillating mode for S<0.62 in the limit of zero momentum thickness, which contradicts with the experimental findings of Kyle and Sreenivasan. The analyses of momentum-dominated jets neglect buoyancy effects because of the small Richardson number. Although this assumption is appropriate in the potential core, the gravitational effects are important in the annular region surrounding the jet, where the density and velocity gradients are large. This reasoning provides basis for the hypothesis that the instability in low Richardosn number jets studied by Kyle and Sreenivasan and Monkewitz et al. is caused by buoyancy. The striking similarity in characteristics of the instability and virtually the identical conclusions reached by Subbarao and Cantwell in buoyant (Ri>0.5) helium jets on one hand and by Kyle and Sreenivasan in momentum-dominated (Ri<1x10(exp -3)) helium jets on the other support this hypothesis. However, quantitative experiments in normal and microgravity are necessary to obtain direct physical evidence of buoyancy effects on the flow instability and structure of momentum-dominated low-density jets. The primary objective of this new research project is to quantify how buoyancy affects the flow instability and structure in the near field of low-density jets. The flow will be described by the spatial and temporal evolutions of the instability, length and time scales of the oscillating mode, and the mean and fluctuating concentration fields. To meet this objective, concentration measurements will be obtained across the whole field using quantitative Rainbow Schlieren Deflectometry, providing spatial resolution of 0.1mm and temporal resolution of 0.017s to 1ms. The experimental effort will be supplemented with linear stability analysis of low-density jets by considering buoyancy. The first objective of this research is to investigate the effects of gravity on the flow instability and structure of low-density jets. The flow instability in these jets has been attributed to buoyancy. By removing buoyancy in our experiments, we seek to obtain the direct physical evidence of the instability mechanism. In the absence of the instability, the flow structure will undergo a significant change. We seek to quantify these changes by mapping the flow field (in terms of the concentration profiles) of these jets at non-buoyant conditions. Such information is presently lacking in the existing literature. The second objective of this research is to determine if the instability in momentum-driven, low-density jets is caused by buoyancy. At these conditions, the buoyancy effects are commonly ignored because of the small Richardson based on global parameters. By eliminating buoyancy in our experiments, globally as well as locally, we seek to examine the possibility that the instability mechanism in self-excited, buoyant or momentum-driven jets is the same. To meet this objective, we would quantify the jet flow in normal and microgravity, while systematically decreasing the Richardson number from buoyancy-driven to momentum driven flow regime. The third objective of this research is to perform a linear stability analysis of low-density gas jets by including the gravitational effects. The flow oscillations in these jets are attributed to an absolute instability, whereby the disturbance grows exponentially at the site to ultimately contaminate the entire flow field. We seek to study the characteristics of both convective and absolute instabilities and demarcate the boundary between them.

Agrawal A. K.↗

Linear viscous stability analysis of the stratified Bickley jet

The stability of the stably stratified Bickley jet to small disturbances in a viscous, heat conducting fluid is explored by numerical integration of the governing differential equation. Curves of neutral stability are presented in the alpha-R plane and the parametric dependence of these curves upon Richardson number and Prandtl number is examined. The results show a curious bimodal behavior of the curve of neutral stability for small Richardson numbers. In addition, closed neutral curves are found for Richardson numbers in excess of the critical value for the inviscid flow. The unstable regions in the interior of these curves increase with decreasing Prandtl number. Energy balances are computed which show the relative magnitude of the Reynolds stress production, buoyancy dissipation, and viscous dissipation terms across the flow. As the thermal diffusivity becomes large compared with the kinematic viscosity, the buoyancy term is suppressed and the thermal stratification becomes irrelevant.

Gage, K. S.↗

Computational fluid dynamic analysis of liquid rocket combustion instability

The paper presents a computational analysis of liquid rocket combustion instability. Consideration is given to both a fully nonlinear unsteady calculation as well as a new CFD-based linearized stability analysis. An analytical solution for the linear stability problem in a constant area combustion chamber with uniform mean flow is developed to verify the numerical analyses.

Venkateswaran, Sankaran↗

Orion Entry, Descent, and Landing Simulation

The Orion Entry, Descent, and Landing simulation was created over the past two years to serve as the primary Crew Exploration Vehicle guidance, navigation, and control (GN&C) design and analysis tool at the National Aeronautics and Space Administration (NASA). The Advanced NASA Technology Architecture for Exploration Studies (ANTARES) simulation is a six degree-of-freedom tool with a unique design architecture which has a high level of flexibility. This paper describes the decision history and motivations that guided the creation of this simulation tool. The capabilities of the models within ANTARES are presented in detail. Special attention is given to features of the highly flexible GN&C architecture and the details of the implemented GN&C algorithms. ANTARES provides a foundation simulation for the Orion Project that has already been successfully used for requirements analysis, system definition analysis, and preliminary GN&C design analysis. ANTARES will find useful application in engineering analysis, mission operations, crew training, avionics-in-the-loop testing, etc. This paper focuses on the entry simulation aspect of ANTARES, which is part of a bigger simulation package supporting the entire mission profile of the Orion vehicle. The unique aspects of entry GN&C design are covered, including how the simulation is being used for Monte Carlo dispersion analysis and for support of linear stability analysis. Sample simulation output from ANTARES is presented in an appendix.

Hoelscher, Brian R.↗

A Stability Analysis of Cylindrical Panels Using a Finite Element Formulation

A cylindrical finite element suitable for the linear stability analysis of cylindrical shells is developed. Energy principles and variational methods lead to a problem formulation which lends itself to physical interpretations of the governing matrices of the finite element. By properly grouping the terms which result from taking the second variation of the potential energy of the element, it is possible to identify three distinct types of matrices. The first matrix is the conventional stiffness matrix; the second is an initial stress stiffness matrix; and the third is an initial displacement stiffness matrix. With the assumption of linearity, the buckling problem is stated in terms of the classical linear real eigenvalue equation. This problem formulation was programmed on the CDC 6600 series computer. The computer program is used to analyze the buckling of a variety of structures. Columns, arches, flat plates and curved panels with and without cutouts are considered.

Snyder, R. E.↗

Stability analysis of two coupled Lorenz lasers and the coupling-induced periodic to chaotic transition

A linear stability analysis of two Lorenz lasers coupled by their electric fields has been performed, and it is shown that the bad cavity condition becomes a function of coupling and that a good cavity instability may occur if the injected fields are inverted before injection. In addition, it is shown that the symmetrically coupled Lorenz system is isomorphic to the original Lorenz system with new parameters. The stability analysis also predicts a lowering of the second laser threshold with coupling for both the chaotic and self-pulsing regimes. Numerical integration of the equations is in agreement with these predictions and has revealed a coupling induced transition from self-pulsing to chaotic behavior. The classification of the behavior of the coupled system in the parameter space of the coupling constants has been investigated and shows that the results of symmetric coupling allow enough of a margin for an experimental test of the theory. This would allow experimentalists to observe the actual Lorenz instability at excitations as low as 4-5 times above threshold.

Lawandy, N. M.↗

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↗

In-flight detection of Tollmien-Schlichting instabilities in laminar flow

This paper presents the results of a recent flight experiment conducted to evaluate surface-mounted hot-film sensors designed to detect Tollmien-Schlichting instability in the boundary layer in flight. Flight data are presented to illustrate the design and operation of the sensors and to provide information on disturbance growth and transition mode. The results of calculations using a boundary-layer analysis code and boundary-layer linear stability analysis methods are included.

Croom, Cynthia C.↗

Effect of surface tension on the onset of convection in a double-diffusive layer

The effect of surface tension on the stability of a double-diffusive layer is considered using linear stability analysis. The surface tension is assumed to vary linearly with temperature and solute concentration. The eigenvalue problem is solved by the Galerkin method. Results show that the predicted stability boundary based on Marangoni effects alone is completely altered in the presence of buoyancy effects induced by low gravity levels (about 10 exp -5 g). At reduced gravity levels, salt-finger instability may onset in the overstable mode due to the stabilizing effect of surface tension. Fluid properties in terms of the Prandtl and the Lewis numbers have a profound effect on the stability conditions; opposite stability characteristics are found in salt solutions and in molten metals.

Chen, C. F.↗

Design optimization of natural laminar flow fuselages in compressible flow

An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions, and linear stability analysis of the compressible boundary layer. Boundary-layer transition is predicted by a 'hybrid' transition criterion based on Granville's transition criterion and a criterion using linear stability theory coupled with the e(n)-method. A tiptank of a business-jet is used as an example to illustrate that the method can be utilized to design an axisymmetric body shape with extensive natural laminar flow. Boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10 to the 6th on the original tiptank. On the designed body shape, a transition Reynolds number of 7.22 x 10 to the 6th is predicted using compressible linear stability theory coupled with e(n)-method.

Dodbele, Simha S.↗

Design of laminar flow bodies in compressible flow

An optimization method was developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions, and linear stability analysis of the compressible boundary-layer. Boundary-layer transition is predicted by a 'hybrid' transition criterion based on Granville's transition criterion and a criterion using linear stability theory coupled with the e(sup n)-method. A tiptank of a business jet is used as an example to illustrate that the method can be utilized to design an axisymmetric body shape with extensive natural laminar flow. On the original tiptank boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10 6 on the original tiptank. On the designed body shape a transition Reynolds number of 7.22 x 10 6 is predicted using compressible linear stability theory coupled with e n -method.

Simha S. Dodbele↗

On Time Delay Margin Estimation for Adaptive Control and Optimal Control Modification

This paper presents methods for estimating time delay margin for adaptive control of input delay systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent an adaptive law by a locally bounded linear approximation within a small time window. The time delay margin of this input delay system represents a local stability measure and is computed analytically by three methods: Pade approximation, Lyapunov-Krasovskii method, and the matrix measure method. These methods are applied to the standard model-reference adaptive control, s-modification adaptive law, and optimal control modification adaptive law. The windowing analysis results in non-unique estimates of the time delay margin since it is dependent on the length of a time window and parameters which vary from one time window to the next. The optimal control modification adaptive law overcomes this limitation in that, as the adaptive gain tends to infinity and if the matched uncertainty is linear, then the closed-loop input delay system tends to a LTI system. A lower bound of the time delay margin of this system can then be estimated uniquely without the need for the windowing analysis. Simulation results demonstrates the feasibility of the bounded linear stability method for time delay margin estimation.

Nguyen, Nhan T.↗

Design optimization of natural laminar flow bodies in compressible flow

An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions and linear stability analysis of the compressible boundary-layer. In order to reduce the computer time, Granville's transition criterion is used to predict boundary-layer transition and to calculate the gradients of the objective function, and linear stability theory coupled with the e(exp n)-method is used to calculate the objective function at the end of each design iteration. Use of a method to design an axisymmetric body with extensive natural laminar flow is illustrated through the design of a tiptank of a business jet. For the original tiptank, boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10(exp 6). For the designed body shape, a transition Reynolds number of 7.22 x 10(exp 6) is predicted using compressible linear stability theory coupled with the e(exp n)-method.

Dodbele, Simha S.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper reports the results of a comprehensive linear stability analysis of the boundary layer flow over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien-Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. However, a significant scatter ( N ∈ (4,9)) was observed in the values of the logarithmic amplification factor along the measured transition front. This previous analysis was based on an approximate basic state (based on a boundary layer code with conical flow approximation) and parallel stability computations without surface curvature effects. Here, we examine the effects of these various approximations with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, both linear stability theory (LST) and the parabolized stability equations (PSE)are used in conjunction with an accurate definition of the laminar boundary layer flow as computed with a Navier-Stokes solver with a Reynolds-Averaged-Navier-Stokes (RANS) based turbulence model within the turbulent parts of the flow. Furthermore, the effects of instability wave propagation within a fully three-dimensional boundary layer are also evaluated by integrating the disturbance growth rates along suitably chosen, curvilinear (i.e., nonplanar) propagation trajectories. The results of this analysis are also used in an accompanying paper by Venkatachari et al. to develop improved, physics based transition predictions for the same CRM-NLF configuration.

Boundary layer transition↗

Single and double point modeling of homogeneous turbulence

Investigations carried out for revisiting homogeneous turbulent flows in the presence of mean shear, rotation, or external compression are summarized. The simplest and most concise RDT (Rapid Distortion Theory) formulation, which includes a comprehensive linear stability analysis, is used for this purpose. Such a linear approach could be extended by a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) to two point closure in order to model non-linear interactions, especially when pure Coriolis effects are present. The results are discussed in connection with databases obtained by DNS (Direct Numerical Simulations), including previous Center for Turbulence Research (CTR) results and new calculations in progress. The main goal is to contribute to and significantly improve on the rational one-point closure models in progress at the CTR and at Ecole Centrale de Lyon (ECL).

Cambon, C.↗

The three-dimensional evolution of a plane wake

In the past three decades, linear stability analysis has led to a comprehensive understanding of the linear stages of transition in plane wakes. Our understanding of the nonlinear and turbulent stages is less developed. Nonlinear theory developed by Papageorgiou and Smith was used to study the long-wavelength regime in wakes. The nonlinear and turbulent stages were investigated experimentally, and few numerical studies examined the early nonlinear stages of forced wakes. The evolution of three dimensional disturbances in an incompressible wake is investigated using direct numerical simulations. The instantaneous three-dimaensional structures and corresponding statistics are presented.

Maekawa, H.↗

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↗

Geometry of thin liquid sheet flows

Incompresible, thin sheet flows have been of research interest for many years. Those studies were mainly concerned with the stability of the flow in a surrounding gas. Squire was the first to carry out a linear, invicid stability analysis of sheet flow in air and compare the results with experiment. Dombrowski and Fraser did an experimental study of the disintegration of sheet flows using several viscous liquids. They also detected the formulation of holes in their sheet flows. Hagerty and Shea carried out an inviscid stability analysis and calculated growth rates with experimental values. They compared their calculated growth rates with experimental values. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. Brown experimentally investigated thin liquid sheet flows as a method of application of thin films. Clark and Dumbrowski carried out second-order stability analysis for invicid sheet flows. Lin introduced viscosity into the linear stability analysis of thin sheet flows in a vacuum. Mansour and Chigier conducted an experimental study of the breakup of a sheet flow surrounded by high-speed air. Lin et al. did a linear stability analysis that included viscosity and a surrounding gas. Rangel and Sirignano carried out both a linear and nonlinear invisid stability analysis that applies for any density ratio between the sheet liquid and the surrounding gas. Now there is renewed interest in sheet flows because of their possible application as low mass radiating surfaces. The objective of this study is to investigate the fluid dynamics of sheet flows that are of interest for a space radiator system. Analytical expressions that govern the sheet geometry are compared with experimental results. Since a space radiator will operate in a vacuum, the analysis does not include any drag force on the sheet flow.

Chubb, Donald L.↗