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 1,675 records · Page 93

Rocket Combustion Instability Suppression With Cavity Resonators

Properly tuned Helmholtz resonator damping devices are an effective method for eliminating high frequency combustion instability in rocket engine thrust chambers. A program was implemented with the objective of extending Helmholtz resonator theory to thrust chamber conditions to establish a standard design procedure. By using systematic cold-flow testing and hot firings, empirical design equations were formulated. An effort is currently underway to extend the design theory to other types of acoustical damping devices such as quarter-wave tubes and slots.

R H Counts↗

On the interaction of a sound pulse with the shear layer of an axisymmetric jet. II - Heated jets

The fluctuating field of a heated jet excited by a sound pulse is simulated numerically. The fluctuations in both the flow field and the far field are studied. The flow field results depend crucially on the stability properties of the heated jet. The altered stability properties due to the heating cause a reduction in the overall fluctuations and a shift of the instability waves into lower frequencies. These changes cause a similar downward shift in the far field spectrum and a reduction in the total far field noise. The results provide a partial explanation in terms of stability theory of experimentally observed properties of the noise of heated jets.

Bayliss, A.↗

Plasma transport in the auroral return current region

The classical and anomalous transport properties of a multifluid plasma consisting of H(+), O(+), and electron populations in the presence of auroral field-aligned return currents are investigated, using a multimoment fluid model with anomalous transport coefficients. The macroscopic effects of the electrostatic ion cyclotron (EIC) instability and of an EIC-related anomalous resistivity mechanism which heats the electrons are included in the present version of the model. The responses of the outflowing polar wind plasma to the application of current, with and without instabilities, are exhibited. The simulations show that the electron drift velocity corresponding to a return current of 0.65 micro-A/sq m is above the threshold for EIC waves. Downward electron heat flow competes with upward convection and adiabatic effects to determine the direction of the electron temperature anisotropy. Resistive electron heating lowers the critical drift velocity for marginal EIC stability and leads to enhanced ion heating.

Ganguli, Supriya B.↗

Ejection mechanisms in the sublayer of a turbulent channel

A possible model for the inception of vorticity ejections in the viscous sublayer of a turbulent rectangular channel is presented. It was shown that this part of the flow is dominated by protruding strong shear layers of z-vorticity, and it was proposed as a mechanism for their maintenance and reproduction which is essentially equivalent to that responsible for the instability of 2-D Tollmien-Schlichting waves. The efforts to isolate computationally a single structure for its study have failed up to now, since it appears that single structures decay in the absence of external forcing, but a convenient computation model was identified in the form of a long and narrow periodic computational box containing at each moment only a few structures. Further work in the identification of better reduced systems is in progress.

Jimenez, J.↗

Transonic adaptive flutter suppression using approximate unsteady time domain aerodynamics

A digital adaptive controller is applied to the active flutter suppression problem of a wing under time varying flight conditions in subsonic and transonic flow. Linear quadratic controller gain at each time step is obtained using an iterative Riccati solver. The digital adaptive optimal controller is robust with respect to the unknown external loads. Flutter and divergence instabilities are simultaneously suppressed using a trailing-edge control surface and displacement sensing. A new transonic unsteady aerodynamic approximation methodology is developed which enables one to carry out the rapid calculation required for transonic aeroservoelastic applications. This approximation is based on a combination of unsteady subsonic aerodynamics combined with a transonic correction procedure. Aeroservoelastic transient time response is obtained using Roger's approximation, state transition matrices and an iterative time marching algorithm. The aeroservoelastic system in the time domain is modeled using a deterministic ARMA model together with a parameter estimator. Transonic flutter boundaries of a wing structure are computed, in the time domain, using an estimated aeroelastic system matrix and are in good agreement with experimental data for the low transonic Mach number range.

Pak, Chan-Gi↗

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗

Performance of a low-power subsonic-arc-attachment arcjet thruster

A subsonic-arc-attachment thruster design was scaled from a 30 kW 1960's vintage thruster to operate at nominally 3 kW. Performance measurements were obtained over a 1-4 kW power range using hydrogen as the propellant. Several modes of operation were identified and were characterized by varying degrees of voltage instability. A stability map was developed showing that the voltage oscillations were brought upon by elevated current or propellant levels. At a given specific energy level the specific impulse increased asymptotically with increased flow rates. Comparisons of performance were made between radial and tangential propellant injection. When the vortex flow was eliminated using radial injection, the operating voltages were lower at a given current, and the specific impulse and efficiency decreased. Tests were also conducted to determine the effects of background pressure on operation. Performance data were obtained at pressures of 0.047 Pa and 18 Pa. For a given specific energy level, the performance increased with a decrease in facility background pressure. Lowering the background pressure also caused a dramatic change in the voltage-current characteristic and the voltage stability, a phenomenon not previously reported with conventional supersonic-arc-attachment thrusters.

Sankovic, John M.↗

On the Resistance of the Air at High Speeds and on the Automatic Rotation of Projectiles

Here, the laws governing the flow of a compressible fluid through an opening in a thin wall are applied to the resistance of the air at high speeds, especially as applied to the automatic rotation of projectiles. The instability which we observe in projectiles shot into the air without being given a moment of rotation about their axis of symmetry, or without stabilizing planes, is a phenomenon of automatic rotation. It is noted that we can prevent this phenomenon of automatic rotation by bringing the center of gravity sufficiently near one end, or by fitting the projectile with stabilizing planes or a tail. The automatic rotation of projectiles is due to the suction produced by the systematic formation of vortices behind the extremity of the projectile moving with the wind.

Riabouchinski, D↗

Sound radiation due to boundary layer transition

This report describes progress made to date towards calculations of noise produced by the laminar-turbulence transition process in a low Mach number boundary layer formed on a rigid wall. The primary objectives of the study are to elucidate the physical mechanisms by which acoustic waves are generated, to clarify the roles of the fluctuating Reynolds stress and the viscous stress in the presence of a solid surface, and to determine the relative efficiency as a noise source of the various transition stages. In particular, we will examine the acoustic characteristics and directivity associated with three-dimensional instability waves, the detached high-shear layer, and turbulent spots following a laminar breakdown. Additionally, attention will be paid to the unsteady surface pressures during the transition, which provide a source of flow noise as well as a forcing function for wall vibration in both aeronautical and marine applications.

Wang, Meng↗

Performance of a low-power subsonic-arc-attachment arcjet thruster

A subsonic-arc-attachment thruster design was scaled from a 30 kW 1960's vintage thruster to operate at nominally 3 kW. Performance measurements were obtained over a 1-4 kW power range using hydrogen as the propellant. Several modes of operation were identified and were characterized by varying degrees of voltage instability. A stability map was developed showing that the voltage oscillations were brought upon by elevated current or propellant levels. At a given specific energy level the specific impulse increased asymptotically with increased flow rates. Comparisons of performance were made between radial and tangential propellant injection. When the vortex flow was eliminated using radial injection, the operating voltages were lower at a given current, and the specific impulse and efficiency decreased. Tests were also conducted to determine the effects of background pressure on operation, and performance data were obtained at pressures of 0.047 Pa and 18 Pa. For a given specific energy level, the performance increased with a decrease in facility background pressure. Lowering the background pressure also caused a dramatic change in the voltage-current characteristic and the voltage stability, a phenomenon not previously reported with conventional supersonic-arc-attachment thrusters.

Sankovic, John M.↗

Convection in the Physical Vapor Transport Process: Thermal - 1

The effects of convection on diffusive-convective physical vapor transport process are examined computationally. We analyze conditions ranging from typical laboratory conditions to conditions achievable only in a low gravity environment. This corresponds to thermal Rayleigh numbers Ra(sub tau) ranging from 1.80 x 10 to 1.92 x 10(exp 6). Our results indicate that the effect of the sublimation and condensation fluxes at the boundaries is to increase the threshold of instability. For typical ground based conditions, time dependent oscillatory convection can occur. This results in unsteady transport, and non-uniform temperature and concentration gradients at the crystal interface. Spectral analysis of the flow field shows parametric regions exhibiting both an oscillatory approach to steady state and a chaotic transient to a periodic state. Low gravity conditions stabilize the flow field. Convective effects are effectively reduced, thus resulting in uniform temperature and concentration gradients at the interface, a desirable condition for crystal growth.

Duval, Walter M. B.↗

Rotor Wake Development During the First Revolution

The wake behind a two-bladed model rotor in light climb was measured using particle image velocimetry, with particular emphasis on the development of the trailing vortex during the first revolution of the rotor. The distribution of vorticity was distinguished from the slightly elliptical swirl pattern. Peculiar dynamics within the void region may explain why the peak vorticity appeared to shift away from the center as the vortex aged, suggesting the onset of instability. The swirl and axial velocities (which reached 44 and 12 percent of the rotor-tip speed, respectively) were found to be asymmetric relative to the vortex center. In particular, the axial flow was composed of two concentrated zones moving in opposite directions. The radial distribution of the circulation rapidly increased in magnitude until reaching a point just beyond the core radius, after which the rate of growth decreased significantly. The core-radius circulation increased slightly with wake age, but the large-radius circulation appeared to remain relatively constant. The radial distributions of swirl velocity and vorticity exhibit self-similar behaviors, especially within the core. The diameter of the vortex core was initially about 10 percent of the rotor-blade chord, but more than doubled its size after one revolution of the rotor. According to vortex models that approximate the measured data, the core-radius circulation was about 79 percent of the large-radius circulation, and the large-radius circulation was about 67 percent of the maximum bound circulation on the rotor blade. On average, about 53 percent of the maximum bound circulation resides within the vortex core during the first revolution of the rotor.

McAlister, Kenneth W.↗

Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods

A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.

Atkins, Harold L.↗

Stability of absorption phenomena in laser-thermal propulsion

The mean flow and stability characteristics of laser absorption phenomena in a choked converging-diverging nozzle are considered. Calculations are presented for a given nozzle geometry and a series of laser intensities. Gas absorptivities corresponding to a hydrogen-cesium mixture are used with different initial temperatures being selected to investigate the effects of changes in the shape of the k-T curve. Both stability and mean flow calculations are limited to the one-dimensional case. The mean flow results show a decrease in mass flow as laser power is increased, along with increasingly steep temperature profiles. Calculations span regions of partial and complete absorption. One region is found where multiple solutions exist. Local stability results indicate the u-c characteristic is the only unstable mode in the unheated case. Laser heat addition makes this mode more unstable and also destabilizes the u-characteristic. Numerical calculations of disturbance propagation show that the instability of the u-c disturbances is counteracted by their reflection to u + c disturbances at the upstream end. The growth of the u-disturbances is localized in regions where the temperature profile is steep and they are damped in other regions. The increasing destabilization that is observed with increased laser power is probably the reason for difficulty in obtaining converged mean flow solutions at high laser intensities.

Merkle, C. L.↗

Rotordynamic Instability Problems in High-Performance Turbomachinery, 1993

The first rotordynamics workshop proceedings (NASA CP-2133, 1980) emphasized a feeling of uncertainty in predicting the stability of characteristics of high-performance turbomachinery. In the second workshop proceedings (NASA CP-2250, 1982) these uncertainties were reduced through programs established to systematically resolve problems, with emphasis on experimental validation of the forces that influence rotordynamics. In the third proceedings (NASA CP-2338, 1984) many programs for predicting or measuring forces and force coefficients in high-performance turbomachinery produced results. Data became available for designing new machines with enhanced stability characteristics or for upgrading existing machines. In the fourth proceedings (NASA CP-2443, 1986) there emerged trends toward a more unified view of rotordynamic instability problems and several encouraging new analytical developments. The fifth workshop (NASA CP-3026, 1988) supported the continuing trend toward a unified view with several new developments in the design and manufacture of new turbomachineries with enhanced stability characteristics along with new data and associated numerical/theoretical results. The sixth workshop report (NASA CP-3 I 22, 1990) field experience and experimental results, and expanded the use of computational and control techniques with integration of damper, bearing, and eccentric seal operation results. The present workshop continues to report field experiences, numerical, theoretical, and experimental results and control methods for seals, bearings, and dampers with some attention given to variable thermo physical properties and turbulence measurements, and introduction of two-phase flow results. The intent of the workshop and this proceedings is to provide a continuing impetus for an understanding and resolution of these problems.

Instabilities↗

Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.

Culick, F. E. C.↗

Unstable vortices in the near region of an internal flow cavity

This paper presents experimental data taken in the forward region of a separated internal free-shear layer produced in an internal cavity flowfield. It has been found that in the region very near the forward restrictor, experimental velocity profiles agree closely with the exact Stuart instability velocity profile with various values of a steepness parameter. Reynolds shear-stress profiles suggest the presence of counter-rotating longitudinal vortices. Spectral analysis by the maximum entropy method of the time samples within the vortices indicates subharmonic and harmonic components of the fundamental frequency, with a weak indication of the fundamental frequency itself.

Isaacson, L. K.↗

Excitation of Crossflow Instabilities in a Swept Wing Boundary Layer

The problem of crossflow receptivity is considered in the context of a canonical 3D boundary layer (viz., the swept Hiemenz boundary layer) and a swept airfoil used recently in the SWIFT flight experiment performed at Texas A&M University. First, Hiemenz flow is used to analyze localized receptivity due to a spanwise periodic array of small amplitude roughness elements, with the goal of quantifying the effects of array size and location. Excitation of crossflow modes via nonlocalized but deterministic distribution of surface nonuniformity is also considered and contrasted with roughness induced acoustic excitation of Tollmien-Schlichting waves. Finally, roughness measurements on the SWIFT model are used to model the effects of random, spatially distributed roughness of sufficiently small amplitude with the eventual goal of enabling predictions of initial crossflow disturbance amplitudes as functions of surface roughness parameters.

Carpenter, Mark H.↗