Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow instability”

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 721 records · Page 40

On the secondary instability of the most dangerous Goertler vortex

Recent studies have demonstrated the most unstable Goertler vortex mode is found in flows, both two and three-dimensional, with regions of (moderately) large body curvature and these modes reside within a thin layer situated at the base of the conventional boundary layer. Further work concerning the nonlinear development of the most dangerous mode demonstrates that the flow results in a self induced flow reversal. However, prior to the point at which flow reversal is encountered, the total streamwise velocity profile is found to be highly inflectional in nature. Previous work then suggests that the nonlinear vortex state will become unstable to secondary, inviscid, Rayleigh wave instabilities prior to the point of flow reversal. Our concern is with the secondary instability of the nonlinear vortex states, which result from the streamwise evolution of the most unstable Goertler vortex mode, with the aim of determining whether such modes can induce a transition to a fully turbulent state before separation is encountered.

Otto, S. R.↗

Computation of multi-dimensional viscous supersonic flow

A method has been developed for two- and three-dimensional computations of viscous supersonic jet flows interacting with an external flow. The approach employs a reduced form of the Navier-Stokes equations which allows solution as an initial-boundary value problem in space, using an efficient noniterative forward marching algorithm. Numerical instability associated with forward marching algorithms for flows with embedded subsonic regions is avoided by approximation of the reduced form of the Navier-Stokes equations in the subsonic regions of the boundary layers. Supersonic and subsonic portions of the flow field are simultaneously calculated by a consistently split linearized block implicit computational algorithm. The results of computations for a series of test cases associated with supersonic jet flow is presented and compared with other calculations for axisymmetric cases. Demonstration calculations indicate that the computational technique has great promise as a tool for calculating a wide range of supersonic flow problems including jet flow. Finally, a User's Manual is presented for the computer code used to perform the calculations.

Buggeln, R. C.↗

Shear-induced instability and arch filament eruption - A magnetohydrodynamic (MHD) numerical simulation

A situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere undergoes symmetrical shear motion at the footpoints is investigated via a 2D (nonplanar) MHD simulation. It was found that the vertical plasma flow velocities grow exponentially, leading to a new type of global MHD instability. The growth rate increases almost linearly until it reaches the same order of magnitude as the Alfven speed. Then a nonlinear MHD instability occurs beyond this point. It was found that the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. The nonlinear dynamical shearing instability is illustrated by a numerical example that is given for three different values of the plasma beta that span several orders of magnitude.

Wu, S. T.↗

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the nonlocal, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, such as those in the vicinity of roughness elements, steps, gaps, or corners. To compute the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 10^(7) to 10^(8) even in relatively simple flows. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear systems. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE. Additionally, inflow-resolvent analysis based on the HLNSE for discrete-roughness-induced disturbances in the nose tip of a blunt cone at Mach 6 demonstrates the importance of including the disturbance amplification along the near vicinity of the roughness element and separation region.

Boundary Layer Stability↗

Tone excited jets. V - A theoretical model and comparison with experiment

A mathematical model of tone-excited jets is developed which consists of two major components: a mathematical description of the process by which the intrinsic instability waves of the jet are excited by the upstream tones; and the modeling of the nonlinear interaction between the mean flow of the jet, the excited large-scale instability waves or turbulence structure, and the fine scale turbulence. It is assumed that each of these jet flow components can be characterized by a few parameters, which are then related by conservation equations that are supplemented by closure models. This quasi-linear model's results are compared with experimental measurements, and good agreement is obtained over a wide range of excitation frequencies and excitation levels.

Tam, C. K. W.↗

On the interaction of stationary crossflow vortices and Tollmien-Schlichting waves in the boundary layer on a rotating disc

There are many fluid flows where the onset of transition can be caused by different instability mechanisms which compete among themselves. The interaction is considered of two types of instability mode (at an asymptotically large Reynolds number) which can occur in the flow above a rotating disc. In particular, the interaction is examined between lower branch Tollmien-Schlichting (TS) waves and the upper branch, stationary, inviscid crossflow vortex whose asymptotic structure has been described by Hall (1986). This problem is studied in the context of investigating the effect of the vortex on the stability characteristics of a small TS wave. Essentially, it is found that the primary effect is felt through the modification to the mean flow induced by the presence of the vortex. Initially, the TS wave is taken to be linear in character and it is shown (for the cases of both a linear and a nonlinear stationary vortex) that the vortex can exhibit both stabilizing and destabilizing effects on the TS wave and the nature of this influence is wholly dependent upon the orientation of this latter instability. Further, the problem is examined with a larger TS wave, whose size is chosen so as to ensure that this mode is nonlinear in its own right. An amplitude equation for the evolution of the TS wave is derived which admits solutions corresponding to finite amplitude, stable, traveling waves.

Bassom, Andrew P.↗

On the interaction of stationary crossflow vortices and Tollmien-Schlichting waves in the boundary layer on a rotating disc

There are many fluid flows where the onset of transition can be caused by different instability mechanisms which compete among themselves. The interaction is considered of two types of instability mode (at an asymptotically large Reynolds number) which can occur in the flow above a rotating disc. In particular, the interaction is examined between lower branch Tollmien-Schlichting (TS) waves and the upper branch, stationary, inviscid crossflow vortex whose asymptotic structure has been described by Hall (1986). This problem is studied in the context of investigating the effect of the vortex on the stability characteristics of a small TS wave. Essentially, it is found that the primary effect is felt through the modification to the mean flow induced by the presence of the vortex. Initially, the TS wave is taken to be linear in character and it is shown (for the cases of both a linear and a nonlinear stationary vortex) that the vortex can exhibit both stabilizing and destabilizing effects on the TS wave and the nature of this influence is wholly dependent upon the orientation of this latter instability. Further, the problem is examined with a larger TS wave, whose size is chosen so as to ensure that this mode is nonlinear in its own right. An amplitude equation for the evolution of the TS wave is derived which admits solutions corresponding to finite amplitude, stable, traveling waves.

Bassom, Andrew P.↗

Transonic small disturbances equation applied to the solution of two-dimensional nonsteady flows

Transonic nonsteady flows are of large practical interest. Aeroelastic instability prediction, control figured vehicle techniques or rotary wings in forward flight are some examples justifying the effort undertaken to improve knowledge of these problems is described. The numerical solution of these problems under the potential flow hypothesis is described. The use of an alternating direction implicit scheme allows the efficient resolution of the two dimensional transonic small perturbations equation.

Couston, M.↗

Status of NASA full-scale engine aeroelasticity research

The paper presents data relevant to several types of aeroelastic instabilities which have been obtained using several types of turbojet and turbofan engines. Special attention is given to data relative to separated flow (stall) flutter, choke flutter, and system mode instabilities. The discussion covers the characteristics of these instabilities, and a number of correlations are presented that help identify the nature of the phenomena.

Lubomski, J. F.↗

Supercritical transition in plane channel flow with spatially periodic perturbations

Laboratory experiments and numerical simulations have been conducted for plane channel flow with a streamwise-periodic array of cylinders. The primary transition in this open flow occurs as a convective rather than absolute instability and leads to traveling-wave packets, which advect out of the system. The ordered secondary state is characteristic of closed flows, in contrast with other open flows where the primary transition often leads directly to turbulence.

Schatz, Michael F.↗

Stability of oscillatory two phase Couette flow

We investigate the stability of two phase Couette flow of different liquids bounded between plane parallel plates. One of the plates has a time dependent velocity in its own plane, which is composed of a constant steady part and a time harmonic component. In the absence of time harmonic modulations, the flow can be unstable to an interfacial instability if the viscosities are different and the more viscous fluid occupies the thinner of the two layers. Using Floquet theory, we show analytically in the limit of long waves, that time periodic modulations in the basic flow can have a significant influence on flow stability. In particular, flows which are otherwise unstable for extensive ranges of viscosity ratios, can be stabilized completely by the inclusion of background modulations, a finding that can have useful consequences in many practical applications.

Coward, Adrian V.↗

Active Control of High Frequency Combustion Instability in Aircraft Gas-Turbine Engines

Active control of high-frequency (greater than 500 Hz) combustion instability has been demonstrated in the NASA single-nozzle combustor rig at United Technologies Research Center. The combustor rig emulates an actual engine instability and has many of the complexities of a real engine combustor (i.e. actual fuel nozzle and swirler, dilution cooling, etc.) In order to demonstrate control, a high-frequency fuel valve capable of modulating the fuel flow at up to 1kHz was developed. Characterization of the fuel delivery system was accomplished in a custom dynamic flow rig developed for that purpose. Two instability control methods, one model-based and one based on adaptive phase-shifting, were developed and evaluated against reduced order models and a Sectored-1-dimensional model of the combustor rig. Open-loop fuel modulation testing in the rig demonstrated sufficient fuel modulation authority to proceed with closed-loop testing. During closed-loop testing, both control methods were able to identify the instability from the background noise and were shown to reduce the pressure oscillations at the instability frequency by 30%. This is the first known successful demonstration of high-frequency combustion instability suppression in a realistic aero-engine environment. Future plans are to carry these technologies forward to demonstration on an advanced low-emission combustor.

Corrigan, Bob↗

Mechanisms of Rotating Instability in Axial Compressors Investigated

Rotating instability is a phenomenon that occurs in the tip flow region of axial compressor stages during stable operation. It can be observed in highly staggered rotors with significant tip clearance and is strongest at high-load operating points where the characteristic levels off. In this condition, the single-stage fan under investigation radiates an audible, whistling tone, and wall pressure spectra in the vicinity of the rotor disk exhibit nonrotational components. The graph shows the spectrum of static pressure at a point on the endwall near the leading edge. A hump appears at roughly half of the blade passing frequency (BPF) and is characteristic of rotating instability. A computational model was developed at the NASA Glenn Research Center to investigate the mechanism behind this phenomenon. A three-dimensional steady Navier-Stokes code that has been successfully tested for a wide range of turbomachinery flows was modified to execute a time-accurate simulation of the full annulus of the compressor. At the inlet of the computational domain, the total pressure, total temperature, and two velocity components are specified. Since no unsteady measurements of static pressure or other flow variables were available downstream of the rotor, circumferentially averaged static pressure was specified on the shroud at the outlet of the computational domain. A three-dimensional view of the vortex from the numerical model is shown. Particle traces released near the leading edge tip have rolled up to illustrate the tip clearance vortex. Flow near the trailing edge is pushed forward by the axially reversed flow. It then interacts with the tip clearance flow and the incoming flow and results in the rotating instability vortex, the core of which is illustrated by total pressure shading on planes located successively downstream. The rotating instability vortex is formed periodically midway between the blades and moves toward the pressure side of the passage. The unsteady behavior of this vortex structure is the main mechanism of the rotating instability is shown. The numerical model can be used to detect any possible occurrence of rotating instability when the tip clearance increases during engine service.

Hah, Chunill↗

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

The BOLT-II configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. The present paper is related to the NASA-led roughness experiment on the secondary side of the BOLT-II configuration. First, computations are performed to aid mean-flow based correlations for the transition measurements acquired during an experimental campaign in the Ludwieg Tube facility at the United States Air Force Academy. To that end, both laminar and turbulent flow computations are performed for selected conditions of the experiment with a 1:3 scale model of the slightly shorter BOLT configuration. The variation in well-known correlation parameters for transition due to isolated roughness elements with respect to element position, height, and flow Reynolds number is delineated. Furthermore, the wake instability theory is used to examine the amplification of unstable disturbances behind the centerline trip at the design conditions of the BOLT-II flight experiment and for its dynamically scaled version on the 1:3 scale model of the wind tunnel experiment. To our knowledge, this study represents the first analysis of the wake instabilities in the presence of a prominent entropy layer outside the boundary layer. At the nominal flight design condition of Re ∞ ) = 5.44 and Re ∞ = 2.5 x 10 6 /m, a diamond planform (“pizza-box”) trip with k/δ = 0.70 and planform-halfwidth-to-height ratio of b/k = 3.0 at X/L = 0.254 induces a strong distortion of the boundary-layer flow in the vicinity of the symmetry plane, such that the wake structure includes a centerline streak produced by the flow retardation behind the trip along with multiple additional, prominent vorticity structures representing the necklace vortex system. Stability results reveal that each of these streaks can support a strong amplification of wake instabilities within a frequency band of 95 kHz–255 kHz, with a peak N-factor of greater than 15 beyond X/L ≈ 0.42. The dynamically scaled wind tunnel configuration with Re ∞ = 20.25 x 10 6 /m is also found to support N > 15 at similar axial stations.

BOLT↗

Interface Morphology During Crystal Growth: Effects of Anisotropy and Fluid Flow

The effect of a parallel shear flow and anisotropic interface kinetics on the onset of instability during growth from a supersaturated solution is analyzed. The model used for anisotropy is based on the microscopic picture of step motion. A shear flow (linear Couette flow or asymptotic suction profile) parallel to the crystal solution interface in the same direction as the step motion decreases interface stability. A shear flow counter to the step motion enhances stability and for sufficiently large shear rates the interface is absolutely morphologically stable. For large wave numbers, the perturbed flow field can be neglected and a simple analytic approximation for the stability-instability demarcation is found.

Coriell, S. R.↗

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗

Investigation of Instabilities and Heat Transfer Phenomena in Supercritical Fuels at High Heat Flux and Temperatures

A series of heated tube experiments was performed to investigate fluid instabilities that occur during heating of supercritical fluids. In these tests, JP-7 flowed vertically through small diameter tubes at supercritical pressures. Test section heated length, diameter, mass flow rate, inlet temperature, and heat flux were varied in an effort to determine the range of conditions that trigger the instabilities. Heat flux was varied up to 4 BTU/sq in./s, and test section wall temperatures reached as high as 1950 F. A statistical model was generated to explain the trends and effects of the control variables. The model included no direct linear effect of heat flux on the occurrence of the instabilities. All terms involving inlet temperature were negative, and all terms involving mass flow rate were positive. Multiple tests at conditions that produced instabilities provided inconsistent results. These inconsistencies limit the use of the model as a predictive tool. Physical variables that had been previously postulated to control the onset of the instabilities, such as film temperature, velocity, buoyancy, and wall-to-bulk temperature ratio, were evaluated here. Film temperatures at or near critical occurred during both stable and unstable tests. All tests at the highest velocity were stable, but there was no functional relationship found between the instabilities and velocity, or a combination of velocity and temperature ratio. Finally, all of the unstable tests had significant buoyancy at the inlet of the test section, but many stable tests also had significant buoyancy forces.

Linne, Diane L.↗

Active Suppression of Instabilities in Engine Combustors

A method of feedback control has been proposed as a means of suppressing thermo-acoustic instabilities in a liquid- fueled combustor of a type used in an aircraft engine. The basic principle of the method is one of (1) sensing combustor pressure oscillations associated with instabilities and (2) modulating the rate of flow of fuel to the combustor with a control phase that is chosen adaptively so that the pressure oscillations caused by the modulation oppose the sensed pressure oscillations. The need for this method arises because of the planned introduction of advanced, lean-burning aircraft gas turbine engines, which promise to operate with higher efficiencies and to emit smaller quantities of nitrogen oxides, relative to those of present aircraft engines. Unfortunately, the advanced engines are more susceptible to thermoacoustic instabilities. These instabilities are hard to control because they include large dead-time phase shifts, wide-band noise characterized by amplitudes that are large relative to those of the instabilities, exponential growth of the instabilities, random net phase walks, and amplitude fluctuations. In this method (see figure), the output of a combustion-pressure sensor would be wide-band-pass filtered and then further processed to generate a control signal that would be applied to a fast-actuation valve to modulate the flow of fuel. Initially, the controller would rapidly take large phase steps in order to home in, within a fraction of a second, to a favorable phase region within which the instability would be reduced. Then the controller would restrict itself to operate within this phase region and would further restrict itself to operate within a region of stability, as long as the power in the instability signal was decreasing. In the phase-shifting scheme of this method, the phase of the control vector would be made to continuously bounce back and forth from one boundary of an effective stability region to the other. Computationally, this scheme would be implemented by the adaptive sliding phaser averaged control (ASPAC) algorithm, which requires very little detailed knowledge of the combustor dynamics. In the ASPAC algorithm, the power of the instability signal would be calculated from the wide-bandpass- filtered combustion-pressure signal and averaged over a period of time (typically of the order of a few hundredths of a second) corresponding to the controller updating cycle [not to be confused with the controller sampling cycle, which would be much shorter (typically of the order of 10(exp -4) second)].

Kopasakis, George↗