Mechanism of instabilities of exothermic hypersonic blunt-body flows
Exothermic hypersonic blunt body flow periodic instability mechanism, using ballistic range with schlieren photographic equipment
SEARCH · Search NASA
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.
Exothermic hypersonic blunt body flow periodic instability mechanism, using ballistic range with schlieren photographic equipment
The Joint Army, Navy, NASA, Air Force (JANNAF) Liquid Rocket Combustion Instability Panel was formed in 1988, drawing its members from industry, academia, and government experts. The panel was charted to address the needs of near-term engine development programs and to make recommendations whose implementation would provide not only sufficient data but also the analysis capabilities to design stable and efficient engines. The panel was also chartered to make long-term recommendations toward developing mechanistic analysis models that would not be limited by design geometry or operating regime. These models would accurately predict stability and thereby minimize the amount of subscale testing for anchoring. The panel has held workshops on acoustic absorbing devices, combustion instability mechanisms, instability test hardware, and combustion instability computational methods. At these workshops, research projects that would meet the panel's charter were suggested. The JANNAF Liquid Rocket Combustion Instability Panel's conclusions about the work that needs to be done and recommendations on how to approach it, based on evaluation of the suggested research projects, are presented.
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.
High-speed schlieren cinematography, combined with synchronized pressure transducer records, was used to investigate the mechanism of combustion instabilities leading to flashback. The combustion chamber had an oblong rectangular cross-section to model the essential features of planar flow, and was provided with a rearward facing step acting as a flameholder. As the rich limit was approached, three instability modes were observed: (1) humming - a significant increase in the amplitude of the vortex pattern; (2) buzzing - a large-scale oscillation of the flame; and (3) chucking - a cyclic reformation of the flame, which results in flashback. The mechanism of these phenomena is ascribed to the action of vortices in the recirculation zone and their interactions with the trailing vortex pattern of the turbulent mixing layer behind the step.
A refined Taylor instability model is developed to describe the surface morphology of rhyolite lava flows. The effect of the downslope flow of the lava on the structures resulting from the Taylor instability mechanism is considered. Squire's (1933) transformation is developed for this flow in order to extend the results to three-dimensional modes. This permits assessing why ridges thought to arise from the Taylor instability mechanism are preferentially oriented transverse to the direction of lava flow. Measured diapir and ridge spacings for the Little and Big Glass Mountain rhyolite flows in northern California are used in conjunction with the model in order to explore the implications of the Taylor instability for flow emplacement. The model suggests additional lava flow features that can be measured in order to test whether the Taylor instability mechanism has influenced the flows surface morphology.
Rotor dynamic instability problems in high performance turbomachinery are reviewed. Mechanical instability mechanisms are discussed. Seal forces and working fluid forces in turbomachinery are discussed. Control of rotor instability is also investigated.
There are many fluid flows where the onset of transition can be caused by different instability mechanisms which compete in the nonlinear regime. The interaction of a centrifugal instability mechanism with the viscous mechanism which causes Tollmein-Schlichting waves is discussed. The interaction between these modes can be strong enough to drive the mean state; the interaction is investigated in the context of curved channel flows so as to avoid difficulties associated with boundary layer growth. Essentially it was found that the mean state adjusts itself so that any modes present are neutrally stable even at finite amplitude. In the first instance, the mean state driven by a vortex of short wavelength in the absence of Tollmein-Schlichting waves is considered. It is shown that for a given channel curvature and vortex wavelength, there is an upper limit to the mass flow rate which the channel can support as the pressure gradient is increased. When the Tollmein-Schlichting waves are present, then the nonlinear differential equation to determine the mean state is modified.
Several experimental campaigns have been conducted across a number of Mach-6 facilities on the hypersonic flow around a CCF (CCF) geometry with a 5° half-angle cone and a 12° half-angle flare. These experiments were conducted as part of the NATO STO Research Task Group AVT-346, which is focused on predicting hypersonic boundary-layer transition on complex geometries. Two conventional wind tunnels (AFRL M6LT and ONERA R2Ch) and one quiet tunnel (Purdue BAM6QT) were used to test the same CCF geometry and to study the instabilities in both the boundary layer in the attached parts of the flow and the shear layer above the axisymmetric separation bubble near the cylinder-cone junction. Two nosetip radii (one nominally sharp and one blunt with a 5 mm radius) were tested. For the sharp nosetip case, there was a great deal of agreement between the measurements of both second-mode and shear-layer instabilities across the two conventional facilities. However, the measured spectra for the blunt nosetip case showed more significant differences between the two tunnels, potentially due to an alternate dominant instability mechanism coupled with the variations in the freestream noise spectra. The quiet facility resulted in a flow that remained laminar to much higher freestream unit Reynolds numbers, as well as in instability measurements that had more distinct spectral peaks for the sharp tip case and broadband rises for the blunt one. The instability mechanisms at play in the sharp quiet case were found to be the same as those in the conventional facilities.
Excitation of two miscible-viscous liquids inside a bounded enclosure in a microgravity environment has shown the evolution of quasi-stationary waves of various modes for a range of parameters. We examine computationally the nonlinear dynamics of the system as the interface breakup and bifurcates to resonance structures typified by the Rayleigh-Taylor instability mechanism. Results show that when the mean steady field is much smaller than the amplitude of the sinusoidal excitation, the system behaves linearly, and growth of quasi-stationary waves occurs through the Kelvin-Helmholtz instability mechanism. However, as the amplitude of excitation increases, nonlinearity occurs through subharmonic bifurcation prior to broadband chaos.
Numerical solutions of the thin-layer approximation of the compressible Navier-Stokes equations have been obtained for flows around an ogive-cylinder body with and without a small fixed disturbance placed upstream of the body tip. Locating the disturbance at positions in the flow field upstream of the tip provokes the same range of behavior of the asymmetric flow that was numerically produced earlier by use of a geometrical disturbance on the body tip. Results remain consistent with the presence of a convective instability mechanism, and demonstrate the potential for a precise mapping of the body's receptivity to fixed disturbances in the flow field. Numerical solutions were also obtained for the flow-field responses to impulsive upstream disturbances to determine whether there is a growing response to an asymmetric impulsive disturbance that is consistent with presence of a convective instability mechanism. Results for surface pressure are interpreted with the aid of a mathematical model. The model suggests that the observed growth of surface pressure gradient with time and distance along a ray in response to an asymmetnc impulsive disturbance is in accord with the solution of a Ginzberg-Landau equation, with distinguishing features of the solution being consistent with the convective instability mode of behavior.
Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.
Extensive measurements are made in a transitioning swept-wing boundary layer using hot-film, hot-wire and cross-wire anemometry. The crossflow-dominated flow contains stationary vortices that breakdown near mid-chord. The most amplified vortex wavelength is forced by the use of artificial roughness elements near the leading edge. Two-component velocity and spanwise surface shear-stress correlation measurements are made at two constant chord locations, before and after transition. Streamwise surface shear stresses are also measured through the entire transition region. Correlation techniques are used to identify stationary structures in the laminar regime and coherent structures in the turbulent regime. Basic techniques include observation of the spatial correlations and the spatially distributed auto-spectra. The primary and secondary instability mechanisms are identified in the spectra in all measured fields. The primary mechanism is seen to grow, cause transition and produce large-scale turbulence. The secondary mechanism grows through the entire transition region and produces the small-scale turbulence. Advanced techniques use Linear Stochastic Estimation (LSE) and Proper Orthogonal Decomposition (POD) to identify the spatio-temporal evolutions of structures in the boundary layer. LSE is used to estimate the instantaneous velocity fields using temporal data from just two spatial locations and the spatial correlations. Reference locations are selected using maximum RMS values to provide the best available estimates. POD is used to objectively determine modes characteristic of the measured flow based on energy. The stationary vortices are identified in the first laminar modes of each velocity component and shear component. Experimental evidence suggests that neighboring vortices interact and produce large coherent structures with spanwise periodicity at double the stationary vortex wavelength. An objective transition region detection method is developed using streamwise spatial POD solutions which isolate the growth of the primary and secondary instability mechanisms in the first and second modes, respectively. Temporal evolutions of dominant POD modes in all measured fields are calculated. These scalar POD coefficients contain the integrated characteristics of the entire field, greatly reducing the amount of data to characterize the instantaneous field. These modes may then be used to train future flow control algorithms based on neural networks.
The use of couplings for high speed turbocompressors or pumps is essential to transmit power from the driver. Typical couplings are either of the lubricated gear or dry diaphragm type design. Gear couplings have been the standard design for many years and recent advances in power and speed requirements have pushed the standard design criteria to the limit. Recent test stand and field data on continuous lube gear type couplings have forced a closer examination of design tolerances and concepts to avoid operational instabilities. Two types of mechanical instabilities are reviewed in this paper: (1) entrapped fluid, and (2) gear mesh instability resulting in spacer throw-out onset. Test stand results of these types of instabilities and other directly related problems are presented together with criteria for proper coupling design to avoid these conditions. An additional test case discussed shows the importance of proper material selection and processing and what can happen to an otherwise good design.
Theoretical developments enable suppression of chaotic structual motions. Theory enables prediction, avoidance, and suppression of chaotic vibrations in structures, especially using dynamic feedback stabilization. In new formulation, motion both repeatable and predictable.
An attempt is made to gain insight into the injection boundary instability mechanism as captured by a simple computational tool developed for liquid rocket engines. Considerable evidence is generated supporting a hypothesis that destructive tangential-mode instabilities in liquid-fuel combustors are driven by combustion processes occurring at the propellant-wetted injector face. The driving action is found to be directly analogous to a burning solid propellant surface mechanism, responsible for solid rocket motor instabilities. Since the proposed instability mechanism has not been fundamentally confirmed, experimental research is strongly recommended.
The post-instability can appear in the form of failure of hyperbolicity or in the form of cascade instability (stretching of vorticity in an inviscid fluid or instability of the Navier-Stokes equations). All the cases are characterized by an unlimited decrease of the scale of the motions, in the course of which the derivatives of the corresponding functions tend to infinity, though the functions themselves remain finite. Such an instability shows that the corresponding state of a continuum cannot be properly simulated by smooth functions. Hence, the original model must be corrected by giving up the requirement about differentiability and enlarging the class of functions. Applications to a model of an inviscid fluid and to a model of turbulence are discussed.
The aerodynamic forces on an oscillating airfoil or airfoil-aileron combination of three independent degrees of freedom were determined. The problem resolves itself into the solution of certain definite integrals, which were identified as Bessel functions of the first and second kind, and of zero and first order. The theory, based on potential flow and the Kutta condition, is fundamentally equivalent to the conventional wing section theory relating to the steady case. The air forces being known, the mechanism of aerodynamic instability was analyzed. An exact solution, involving potential flow and the adoption of the Kutta condition, was derived. The solution is of a simple form and is expressed by means of an auxiliary parameter k. The flutter velocity, treated as the unknown quantity, was determined as a function of a certain ratio of the frequencies in the separate degrees of freedom for any magnitudes and combinations of the airfoil-aileron parameters.
The aerodynamic forces on an oscillating airfoil or airfoil-aileron combination of three independent degrees of freedom have been determined. The problem resolves itself into the solution of certain definite integrals, which have been identified as Bessel functions of the first and second kind and of zero and first order. The theory, being based on potential flow and the Kutta condition, is fundamentally equivalent to the conventional wing-section theory relating to the steady case. The air forces being known, the mechanism of aerodynamic instability has been analyzed in detail. An exact solution, involving potential flow and the adoption of the Kutta condition, has been analyzed in detail. An exact solution, involving potential flow and the adoption of the Kutta condition, has been arrived at. The solution is of a simple form and is expressed by means of an auxiliary parameter K.