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 217 records · Page 12

On Transients in Detached Bridgman Growth

In detached Bridgman growth, a gap exists between the growing crystal and the crucible wall. According to crystal shape stability theory, only specific gap widths will be dynamically stable. Beginning with a crystal diameter that differs from stable conditions, the transient crystal growth process is analyzed. The transient shapes are calculated assuming that the growth angle is constant. Anisotropy and dynamic contact angle effects are considered. In microgravity, dynamic stability depends only on capillary effects and is decoupled from heat transfer. However, heat transfer will influence the crystal-melt interface shape. The local angles and the crystal-melt-vapor triple junction are analyzed and the applicability of the Herring formula is discussed. A potential microgravity experiment is proposed which would enhance our understanding of the detached growth dynamic stability problem.

Mazuruk, K.↗

Instability of a supersonic shock free elliptic jet

This paper presents a comparison of the measured and the computed spatial stability properties of an aspect ratio 2 supersonic shock free elliptic jet. The shock free nature of the elliptic jet provides an ideal test of validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with linear hydrodynamic stability theory. Both aerodynamic and acoustic data were measured. The data are used to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. A hybrid numerical scheme is applied to solve the Rayleigh problem governing the inviscid linear spatial stability of the jet. The measured mean velocity profiles are used to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles used in the stability analysis. Computational results are presented for several modes of instability at two jet cross sections. The acoustic measurements show that a varicose instability is the jet's perferred mode of motion. The stability analysis predicts that the Strouhal number varies linearly as a function of axial distance in the jet's initial mixing region, which is in good qualitative agreement with previous measurements.

Baty, Roy S.↗

Instability of a supersonic shock-free elliptic jet

A comparison of the measured and the computed spatial stability properties of an aspect ratio 2, supersonic, shock free, elliptic jet is presented. The shock free nature of the elliptic jet furnishes a perfect test of the validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with the linear hydrodynamic stability theory. Aerodynamic and acoustic data are measured to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. The measured mean velocity profiles are employed to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles utilized in the stability analysis. It is shown that the measured frequency associated with peak amplitude noise radiation agrees with that frequency predicted to dominate the large scale structure near the end of the potential core.

Baty, Roy S.↗

Boundary-Layer Transition Prediction Through Loose Coupling of OVERFLOW and LASTRAC

Transition prediction based on linear stability theory is expected to more accurately reflect the causality of transition onset than phenomenological transition models based on RANS-like transport equations. To help achieve the CFD vision 2030 aim of building a CFD tool chain with automated prediction of boundary layer transition, a technique to loosely tie the NASA OVERFLOW CFD solver with the LASTRAC stability analysis tool is described. The coupled solver is then used to compute transition over selected over a flat plate in a freestream with sufficiently low levels of turbulence, NLF(1)-0416 airfoil, the 6:1 prolate spheroid at an angle of attack, and a NASA juncture flow model with symmetric wing configuration. The findings show that the loosely coupled approach can reliably predict the transition location accurately in scenarios that are dominated by a single transition mechanism involving Tollmien-Schlichting instabilities, crossflow instabilities, or separation bubble-induced transition, or include a mixture of selected mechanisms. The toolset presents here appears to be robust to the prescription of the initial transition location, and it can lead to a converged solution in four or five rounds of the mean flow calculation and stability analysis, with minimal input from the user.

boundary layer transition↗

Instabilities and subharmonic resonances of subsonic heated round jets, volume 2

When a jet is perturbed by a periodic excitation of suitable frequency, a large-scale coherent structure develops and grows in amplitude as it propagates downstream. The structure eventually rolls up into vortices at some downstream location. The wavy flow associated with the roll-up of a coherent structure is approximated by a parallel mean flow and a small, spatially periodic, axisymmetric wave whose phase velocity and mode shape are given by classical (primary) stability theory. The periodic wave acts as a parametric excitation in the differential equations governing the secondary instability of a subharmonic disturbance. The (resonant) conditions for which the periodic flow can strongly destabilize a subharmonic disturbance are derived. When the resonant conditions are met, the periodic wave plays a catalytic role to enhance the growth rate of the subharmonic. The stability characteristics of the subharmonic disturbance, as a function of jet Mach number, jet heating, mode number and the amplitude of the periodic wave, are studied via a secondary instability analysis using two independent but complementary methods: (1) method of multiple scales, and (2) normal mode analysis. It is found that the growth rates of the subharmonic waves with azimuthal numbers beta = 0 and beta = 1 are enhanced strongly, but comparably, when the amplitude of the periodic wave is increased. Furthermore, compressibility at subsonic Mach numbers has a moderate stabilizing influence on the subharmonic instability modes. Heating suppresses moderately the subharmonic growth rate of an axisymmetric mode, and it reduces more significantly the corresponding growth rate for the first spinning mode. Calculations also indicate that while the presence of a finite-amplitude periodic wave enhances the growth rates of subharmonic instability modes, it minimally distorts the mode shapes of the subharmonic waves.

Ng, Lian Lai↗

Analogy Between the Collatz Conjecture and Sliding Mode Control

The Collatz Conjecture is a famous mathematics problem that is simple to state and understand but remains unsolved. In this paper, the problem is recast as a discrete-time nonlinear system. and viewed in a new light from the perspective of nonlinear systems and feedback control theory. In particular, connections are made between the Collatz sequence of numbers and the behavior of closed-loop dynamical systems designed using a feedback control method called sliding mode control. Trajectories of such systems are characterized by a reaching phase and a sliding phase, the latter of which exhibits exponential convergence. As sliding mode control design is rooted in Lyapunov stability theory, the analogy suggests a new direction for proving or disproving the Collatz Conjecture. Although several possible approaches are discussed, no formal proof is given in this paper.

Collatz conjecture↗

Theory of Type 3 and Type 2 Solar Radio Emissions

The main features of some current theories of type III and type II bursts are outlined. Among the most common solar radio bursts, type III bursts are produced at frequencies of 10 kHz to a few GHz when electron beams are ejected from solar active regions, entering the corona and solar wind at typical speeds of 0.1c. These beams provide energy to generate Langmuir waves via a streaming instability. In the current stochastic-growth theory, Langmuir waves grow in clumps associated with random low-frequency density fluctuations, leading to the observed spiky waves. Nonlinear wave-wave interactions then lead to secondary emission of observable radio waves near the fundamental and harmonic of the plasma frequency. Subsequent scattering processes modify the dynamic radio spectra, while back-reaction of Langmuir waves on the beam causes it to fluctuate about a state of marginal stability. Theories based on these ideas can account for the observed properties of type III bursts, including the in situ waves and the dynamic spectra of the radiation. Type 11 bursts are associated with shock waves propagating through the corona and interplanetary space and radiating from roughly 30 kHz to 1 GHz. Their basic emission mechanisms are believed to be similar to those of type III events and radiation from Earth's foreshock. However, several sub-classes of type II bursts may exist with different source regions and detailed characteristics. Theoretical models for type II bursts are briefly reviewed, focusing on a model with emission from a foreshock region upstream of the shock for which observational evidence has just been reported.

Robinson, P. A.↗

Aeroacoustic and aerodynamic applications of the theory of nonequilibrium thermodynamics

Recent developments in the field of nonequilibrium thermodynamics associated with viscous flows are examined and related to developments to the understanding of specific phenomena in aerodynamics and aeroacoustics. A key element of the nonequilibrium theory is the principle of minimum entropy production rate for steady dissipative processes near equilibrium, and variational calculus is used to apply this principle to several examples of viscous flow. A review of nonequilibrium thermodynamics and its role in fluid motion are presented. Several formulations are presented of the local entropy production rate and the local energy dissipation rate, two quantities that are of central importance to the theory. These expressions and the principle of minimum entropy production rate for steady viscous flows are used to identify parallel-wall channel flow and irrotational flow as having minimally dissipative velocity distributions. Features of irrotational, steady, viscous flow near an airfoil, such as the effect of trailing-edge radius on circulation, are also found to be compatible with the minimum principle. Finally, the minimum principle is used to interpret the stability of infinitesimal and finite amplitude disturbances in an initially laminar, parallel shear flow, with results that are consistent with experiment and linearized hydrodynamic stability theory. These results suggest that a thermodynamic approach may be useful in unifying the understanding of many diverse phenomena in aerodynamics and aeroacoustics.

Horne, W. Clifton↗

Absolute and convective instability of a cylinder wake

The stability behavior of a circular cross-section cylinder's wake at Reynolds number values of up to 45 is presently investigated by means of local linear stability theory. The steady-wake profiles computed are Navier-Stokes solutions of a uniform, incompressible viscous flow around a cylinder obtained by a spectral method. An absolutely unstable region is found to begin to form at a Reynolds number of about 20, and grows with incresing Reynolds number. The onset of global instability response must be characterized by a critical length of an absolutely unstable region; a critical Reynolds number criterion and preferred frequency are proposed based on linear stability analysis.

Yang, X.↗

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗

Stability of time dependent and spatially varying flows; Proceedings of the Symposium, Hampton, VA, Aug. 19-23, 1985

Papers are presented on the application of stability theory to laminar flow control, secondary instabilities in boundary layers, a Floquet analysis of secondary instability in shear flows, and the generation of Tollmien-Schlichting waves by long wavelength free stream disturbances. Also considered are numerical experiments on boundary-layer receptivity, short-scale inviscid instabilities in the flow past surface-mounted obstacles, wave phenomena in a high Reynolds number compressible boundary layer, and instability of time-periodic flows. Other topics include high frequency Rayleigh instability of Stokes layers, stability and resonance in grooved-channel flows, finite length Taylor Couette flow, and vortical structures in the breakdown stage of transition.

Dwoyer, Douglas L.↗

Task II: Three-dimensional Rotating Stall Inception and Effects of Rotating Tip Clearance Asymmetry in Axial Compressors

The effects of two types of flow non-uniformity on stall inception behavior were assessed with linearized stability analyses of two compressor flow models. Response to rotating tip clearance asymmetries induced by a whirling rotor shaft or rotor height variations were investigated with a two-dimensional flow model. A 3-D compressor model was also developed to study the stability of both full-span and part-span rotating stall modes in annular geometries with radial flow variations. The studies focussed on (1) understanding what compressor designs were sensitive to these types of circumferential and spanwise flow non-uniformities, and (2) situations where 2-D stability theories were inadequate because of 3-D flow effects. Rotating tip clearance non-uniformity caused the greatest performance loss for shafts whirling at the rotating stall frequency. A whirling shaft displacement of 1 percent chord caused the stalling mass flow to rise by as much as 10 percent and the peak pressure rise to decrease by 6 percent. These changes were an order of magnitude larger than for equivalent-sized stationary or rotor-locked clearance asymmetries. Spanwise flow non-uniformity always destabilized the compressor, so that 2-D models over-predicted that stall margin compared to 3-D theory. The difference increased for compressors with larger spanwise variations of characteristic slope and reduced characteristic curvature near the peak. Differences between 2-D and 3-D stall point predictions were generally unacceptable (2 - 4 percent of flow coefficient) for single-stage configurations, but were less than 1 percent for multistage compressors. 2-D analyses predicted the wrong stall mode for specific cases of radial inlet flow distortion, mismatching and annulus area contraction, where higher-order radial modes led to stall. The stability behavior of flows with circumferential or radial non-uniformity was unified through a single stability criterion. The stall point for both cases was set by the integral around the annulus of the pressure rise characteristic slope, weighted by the amplitude of the mode shape. For the case of steady circumferential variations, this criterion reduced to the integrated mean slope (IMS) condition associated with steady inlet distortions. The rotating tip clearance asymmetry model was also used to demonstrate the feasibility of actively controlling the shaft position to suppress rotating stall. In axisymmetric mean flow, this method only stabilized the first harmonic mode, increasing the operating range until surge or higher harmonic modes became unstable.

Suder, Kenneth↗

Prediction and control of transition in hypersonic boundary layers

In this paper, the role of compressible linear stability theory in prediction of boundary layer transition at supersonic and hypersonic speeds is investigated. Computations for sharp cones, using the e exp N method with N = 10, show that the first oblique Tollmien-Schlichting mode is responsible for transition at adiabatic wall conditions for freestream Mach numbers up to 7. For cold walls, the two-dimensional second mode dominates the transition process at lower hypersonic Mach numbers due to the well-known destabilizing effect of cooling on the second mode. It is shown that pressure gradient and suction may be used to stabilize this mode. Some results on the real gas effects on hypersonic boundary-layer stability are presented.

Malik, Mujeeb R.↗

Unstable transition properties of the driven magnetohydrodynamic sheet pinch

The unstable transition behavior of a bounded current-carrying two-dimensional magnetofluid is explored, using the hydrodynamic theory developed for parallel shear flows as a guide. Nonlinear excitation of the higher wavenumbers results in the development of electric current sheets of finite extent, as well as the formation of 'attraction currents' centered at the magnetic O-points. A secondary instability mechanism, the dynamic tearing of the electric current sheet, is also observed. This dynamic tearing leads to sawtoothlike temporal oscillations in certain global quantities. The long-time state of the system resembles a nonlinearly saturated state with significant excitation of many wavenumbers. Some features of this state can be understood by means of a Landau nonlinear stability theory based on certain assumptions about the perturbation energy balance.

Dahlburg, R. B.↗

Near minimum-time maneuvers of the advanced space structures technology research experiment (ASTREX) test article: Theory and experiments

The Phillips Laboratory at the Edwards Air Force Base has developed the Advanced Space Structures Technology Research Experiment (ASTREX) facility to serve as a testbed for demonstrating the applicability of proven theories to the challenges of spacecraft maneuvers and structural control. This report describes the work performed on the ASTREX test article by Texas A&M University under contract NAS119373 as a part of the Control-Structure Interaction (CSI) Guest Investigator Program. The focus of this work is on maneuvering the ASTREX test article with compressed air thrusters that can be throttled, while attenuating structural excitation. The theoretical foundation for designing the near minimum-time thrust commands is based on the generation of smooth, parameterized optimal open-loop control profiles, and the determination of control laws for final position regulation and tracking using Lyapunov stability theory. Details of the theory, mathematical modeling, model updating, and compensation for the presence of 'real world' effects are described and the experimental results are presented. The results show an excellent match between theory and experiments.

Vadali, Srinivas R.↗

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗

Development of a two-dimensional zonally averaged statistical-dynamical model. I - The parameterization of moist convection and its role in the general circulation

The moist convection parameterization used in the GISS 3-D GCM is adapted for use in a two-dimensional (2-D) zonally averaged statistical-dynamical model. Experiments with different versions of the parameterization show that its impact on the general circulation in the 2-D model does not parallel its impact in the 3-D model unless the effect of zonal variations is parameterized in the moist convection calculations. A parameterization of the variations in moist static energy is introduced in which the temperature variations are calculated from baroclinic stability theory, and the relative humidity is assumed to be constant. Inclusion of the zonal variations of moist static energy in the 2-D moist convection parameterization allows just a fraction of a latitude circle to be unstable and enhances the amount of deep convection. This leads to a 2-D simulation of the general circulation very similar to that in the 3-D model. The experiments show that the general circulation is sensitive to the parameterized amount of deep convection in the subsident branch of the Hadley cell. The more there is, the weaker are the Hadley cell circulations and the westerly jets. The experiments also confirm the effects of momentum mixing associated with moist convection found by earlier investigators and, in addition, show that the momentum mixing weakens the Ferrel cell. An experiment in which the moist convection was removed while the hydrological cycle was retained and the eddy forcing was held fixed shows that moist convection by itself stabilizes the tropics, reduces the Hadley circulation, and reduces the maximum speeds in the westerly jets.

Yao, Mao-Sung↗

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗