Investigation on mechanism of critical cavitating flow in liquid jet pumps under operating limits

The critical cavitating flow in liquid jet pumps under operating limits is investigated in this paper. Measurements on the axial pressure distribution along the wall of jet pumps indicate that two-phase critical flow occurs in the throat pipe under operating limits. The entrained flow rate and the distribution of the wall pressure upstream lowest pressure section does not change when the outlet pressure is lower than a critical value. A liquid-vapor mixing shockwave is also observed under operating limits. The wave front moves back and forth in low frequency around the position of the lowest pressure. With the measured axial wall pressures, the Mach number of the two-phase cavitating flow is calculated. It's found that the maximum Mach number is very close to I under operating limits. Further analysis infers a cross-section where Mach number approaches to I near the wave front. Thus, the liquid-vapor mixture velocity should reach the local sound velocity and resulting in the occurrence of operating limits.

Compressible flows at small Reynolds numbers

The problem of the slow viscous flow of a gas past a sphere is considered. The fluid cannot be treated incompressible in the limit when the Reynolds number Re, and the Mach number M, tend to zero in such a way that Re ~ o(M^2 ). In this case, the lowest order approximation to the steady Navier-Stokes equations of motion leads to a paradox discovered by Lagerstrom and Chester. This paradox is resolved within the framework of continuum mechanics using the classical slip condition and an iteration scheme that takes into account certain terms in the full Navier-Stokes equations that drop out in the approximation used by the above authors. It is found however that the drag predicted by the theory does not agree with R. A. Millikan's classic experiments on sphere drag.

The whole question of the applicability of the Navier-Stokes theory when the Knudsen number M/Re is not small is examined. A new slip condition is proposed. The idea that the Navier-Stokes equations coupled with this condition may adequately describe small Reynolds number flows when the Knudsen number is not too large is looked at in some detail. First, a general discussion of asymptotic solutions of the equations for all such flows is given. The theory is then applied to several concrete problems of fluid motion. The deductions from this theory appear to interpret and summarize the results of Millikan over a much wider range of Knudsen numbers (almost up to the free molecular or kinetic limit) than hitherto Believed possible by a purely continuum theory. Further experimental tests are suggested and certain interesting applications to the theory of dilute suspensions in gases are noted. Some of the questions raised in the main body of the work are explored further in the appendices.

An investigation of detached shock waves

This investigation demonstrates an application of a flexible wall nozzle for testing in a supersonic wind tunnel. It is conservative to say that the versatility of this nozzle is such that it warrants the expenditure of time to carefully engineer a nozzle and incorporate it in the wind tunnel as a permanent part of the system. The gradients in the test section were kept within one percent of the calibrated Mach number, however, the gradients occurring over the bodies tested were only ± 0.2 percent in Mach number.

The conditions existing on a finite cone with a vertex angle of 75° were investigated by considering the pressure distribution on the cone and the shape of the shock wave. The pressure distribution on the surface of the 75° cone when based on upstream conditions does not show any discontinuities at the theoretical attachment Mach number.

Both the angle of the shock wave and the pressure distribution of the 75° cone are in very close agreement with the theoretical values given in the Kopal report, (Ref. 3).

The location of the intersection of the sonic line with the surface of the cone and with the shock wave are given for the cone. The blocking characteristics of the GALCIT supersonic wind tunnel were investigated with a series of 60° cones.

General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

Stability of Hypervelocity Boundary Layers

The early stage of laminar-turbulent transition in a hypervelocity boundary layer is studied using a combination of modal linear stability analysis, transient growth analysis, and direct numerical simulation. Modal stability analysis is used to clarify the behavior of first and second mode instabilities on flat plates and sharp cones for a wide range of high enthalpy flow conditions relevant to experiments in impulse facilities. Vibrational nonequilibrium is included in this analysis, its influence on the stability properties is investigated, and simple models for predicting when it is important are described.

Transient growth analysis is used to determine the optimal initial conditions that lead to the largest possible energy amplification within the flow. Such analysis is performed for both spatially and temporally evolving disturbances. The analysis again targets flows that have large stagnation enthalpy, such as those found in shock tunnels, expansion tubes, and atmospheric flight at high Mach numbers, and clarifies the effects of Mach number and wall temperature on the amplification achieved. Direct comparisons between modal and non-modal growth are made to determine the relative importance of these mechanisms under different flow regimes.

Conventional stability analysis employs the assumption that disturbances evolve with either a fixed frequency (spatial analysis) or a fixed wavenumber (temporal analysis). Direct numerical simulations are employed to relax these assumptions and investigate the downstream propagation of wave packets that are localized in space and time, and hence contain a distribution of frequencies and wavenumbers. Such wave packets are commonly observed in experiments and hence their amplification is highly relevant to boundary layer transition prediction. It is demonstrated that such localized wave packets experience much less growth than is predicted by spatial stability analysis, and therefore it is essential that the bandwidth of localized noise sources that excite the instability be taken into account in making transition estimates. A simple model based on linear stability theory is also developed which yields comparable results with an enormous reduction in computational expense. This enables the amplification of finite-width wave packets to be taken into account in transition prediction.

Transonic wind tunnel for investigation of turbomachinery blades

The design of a two-stream wind tunnel was undertaken to allow the simulation and study of certain features of the flow field around the blades of high-speed axial-flow turbomachineries. The mixing of the two parallel streams with designed Mach numbers respectively equal to 1.4 and 0.7 will simulate the transonic Mach number distribution generally obtained along the tips of the first stage blades in large bypass-fan engines.

The GALCIT hypersonic compressor plant will be used as an air supply for the wind tunnel, and consequently the calculations contained in the first chapter are derived from the characteristics and the performance of this plant.

The transonic part of the nozzle is computed by using a method developed by K. O. Friedrichs. This method consists essentially of expanding the coordinates and the characteristics of the flow in power series. The development begins with prescribing, more or less arbitrarily, a Mach number distribution along the centerline of the nozzle. This method has been programmed for an IBM 360 computer to define the wall contour of the nozzle.

A further computation is carried out to correct the contour for boundary layer buildup. This boundary layer analysis included geometry, pressure gradient, and Mach number effects. The subsonic nozzle is calculated {including boundary layer buildup) by using the same computer programs. Finally, the mixing zone downstream of the splitter plate was investigated to prescribe the wall contour correction necessary to ensure a constant-pressure test section.

Simulation of three-dimensional viscous flow in turbomachinery geometries using a solution-adaptive unstructured mesh methodology

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodology. The fully three-dimensional, compressible, Reynolds-averaged Navier-Stokes equations with k-ε turbulence modeling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.

The development of turbine exit flow in a swan-necked inter-stage diffuser

This paper describes both the migration and dissipation of flow phenomena downstream of a transonic high-pressure turbine stage. The geometry of the HP stage exit duct considered is a swan-necked diffuser similar to those likely to be used in future engine designs. The paper contains results both from an experimental programme in a turbine test facility and from numerical predictions. Experimental data was acquired using three fast-response aerodynamic probes capable of measuring Mach number, whirl angle, pitch angle, total pressure and static pressure. The probes were used to make time-resolved area traverses at two axial locations downstream of the rotor trailing edge. A 3D time-unsteady viscous Navier-Stokes solver was used for the numerical predictions. The unsteady exit flow from a turbine stage is formed from rotordependent phenomena (such as the rotor wake, the rotor trailing edge recompression shock, the tip-leakage flow and the hub secondary flow) and vane-rotor interaction dependant phenomena. This paper describes the time-resolved behaviour and three-dimensional migration paths of both of these phenomena as they convect downstream. It is shown that the inlet flow to a downstream vane is dominated by two corotating vortices, the first caused by the rotor tip-leakage flow and the second by the rotor hub secondary flow. At the inlet plane of the downstream vane the wake is extremely weak and the radial pressure gradient is shown to have caused the majority of the high loss wake fluid to be located between the mid-height of the passage and the casing wall. The structure of the flow indicates that between a high pressure stage and a downstream vane simple two-dimensional blade row interaction does not occur. The results presented in this paper indicate that the presence of an upstream stage is likely to significantly alter the structure of the secondary flow within a downstream vane. The paper also shows that vane-rotor interaction within the upstream stage causes a 10° circumferential variation in the inlet flow angle of the 2nd stage vane.

The development of aerodynamic performance measurements in a transient test facility

Transient test facilities offer the potential for the simultaneous study of turbine aerodynamic performance, unsteady flow phenomena and the heat transfer characteristics of a turbine stage. This paper describes the development of aerodynamic performance measurement techniques in the Oxford Rotor Facility (ORF). The solutions to the technological issues involved with transient testing presented in this paper are expected to achieve levels of precision uncertainty comparable with traditional steady flow test rigs. The theoretical background to the measurement of aerodynamic performance is presented together with a comprehensive pre-test uncertainty analysis. The instrumentation scheme for the measurement of stage mass flow rate is discussed in detail, the measurements of shaft power, total inlet enthalpy, and stage pressure ratio are also outlined. The current working section features a 62% scale, 1-1/2 stage, high-pressure shroudless transonic turbine. The required inlet flow conditions are provided by an Isentropic Light Piston Tunnel (ILPT) with a quasi-steady state run time of approximately 70ms. The testing is conducted at engine representative specific speed, pressure ratio, gas-to-wall temperature ratio, Mach number and Reynolds number.

The effect of work processes on the casing heat transfer of a transonic turbine

This paper considers the effect of the rotor tip on the casing heat load of a transonic axial flow turbine. The aim of the research is to understand the dominant causes of casing heat-transfer. Experimental measurements were conducted at engine-representative Mach number, Reynolds number and stage inlet to casing wall temperature ratio. Time-resolved heat-transfer coefficient and gas recovery temperature on the casing were measured using an array of heat-transfer gauges. Time-resolved static pressure on the casing wall was measured using Kulite pressure transducers. Time-resolved numerical simulations were undertaken to aid understanding of the mechanism responsible for casing heat load. The results show that between 35% and 60% axial chord the rotor tip-leakage flow is responsible for more than 50% of casing heat transfer. The effects of both gas recovery temperature and heat transfer coefficient were investigated separately and it is shown that an increased stagnation temperature in the rotor tip gap dominates casing heat-transfer. In the tip gap the stagnation temperature is shown to rise above that found at stage inlet (combustor exit) by as much as 35% of stage total temperature drop. The rise in stagnation temperature is caused by an isentropic work input to the tip-leakage fluid by the rotor. The size of this mechanism is investigated by computationally tracking fluid path-lines through the rotor tip gap to understand the unsteady work processes that occur. Copyright © 2005 by ASME.

Unsteady shock motion on a NACA0012 aerofoil at low reduced frequencies

Using transonic blowdown windtunnel experiments, the 2D unsteady shock motion on a NACA0012 aerofoil is examined at various frequencies typical for helicopter blades in forward flight. The aerofoil is subjected to freestream velocities oscillating periodically between M = 0.66 and M = 0.77. Unsteady pressure traces and schlieren images are analyzed over a range of low reduced frequencies to provide information on shock location and strength throughout the cycle. Unsteady effects were noticeable even at very low reduced frequencies (down to O(0.01). However, through the range of frequencies investigated, and within experimental error, the unsteady shock location showed no discernible lag compared to the quasi-steady behaviour. On the other hand, significant variations were observed in shock strengths with the upstream running part of the cycle (decreasing Mach number) displaying considerably stronger shocks than during the accelerating part of the cycle. It could be shown that this variation in shock strength is primarily caused by the shock motion modifying the relative shock Mach number. As a result is was possible to use the quasi-steady results to predict the unsteady shock behaviour at the frequencies investigated here (below 0(0.1)).

Unsteady shock motion on a NACA0012 aerofoil at low reduced frequencies

Using transonic blowdown windtunnel experiments, the 2D unsteady shock motion on a NACA0012 aerofoil is examined at various frequencies typical for helicopter blades in forward flight. The aerofoil is subjected to freestream velocities oscillating periodically between M = 0.66 and M = 0.77. Unsteady pressure traces and schlieren images are analyzed over a range of low reduced frequencies to provide information on shock location and strength throughout the cycle. Unsteady effects were noticeable even at very low reduced frequencies (down to O(0.01). However, through the range of frequencies investigated, and within experimental error, the unsteady shock location showed no discernible lag compared to the quasi-steady behaviour. On the other hand, significant variations were observed in shock strengths with the upstream running part of the cycle (decreasing Mach number) displaying considerably stronger shocks than during the accelerating part of the cycle. It could be shown that this variation in shock strength is primarily caused by the shock motion modifying the relative shock Mach number. As a result is was possible to use the quasi-steady results to predict the unsteady shock behaviour at the frequencies investigated here (below 0(0.1)).

An experimental study of unsteady separated shock wave boundary layer interactions

An experimental investigation of the unsteady interaction between a turbulent boundary layer and a normal shock wave of strength M∞ = 1.4 subject to periodic forcing in a parallel walled duct has been conducted. Emphasis has been placed on the mechanism by which changes in the global flow field influence the local interaction structure. Static pressure measurements and high speed Schlieren images of the unsteady interaction have been obtained. The pressure rise across the interaction and the appearance of the local SBLI structure have been observed to vary during the cycle of periodic shock wave motion. The magnitude of the pressure rise across the interaction is found to be related to the relative Mach number of the unsteady shock wave as it undergoes periodic motion. Variations in the upstream Influence of the interaction are sensitive to the magnitude and direction of shock wave velocity and acceleration and it is proposed that a viscous lag exists between the point of boundary layer separation and the shock wave position. Further work exploring the implications of these findings is proposed, including studies of the variation in position of the points of boundary layer separation and reattachment.

Micro-vortex generator flow control for supersonic engine inlets

To investigate the flow control potential of micro-vortex generators for supersonic mixed-compression inlets, a basic model experiment has been designed which combines a normal shock wave with a subsonic diffuser. The diffuser is formed by a simple expansion corner, with a divergence angle of 6 degrees. The diffuser entry Mach numbers were M=1.3 and M=1.5 and a number of shock locations relative to the corner position were tested. Flow control was applied in the form of counter-rotating micro-vanes with heights of approximately 20% of boundary layer thickness. Furthermore, corner fences where employed to reduce sidewall effects. It was found that micro-vortex generators were able to significantly reduce the extent of flow separation under all conditions, but could not eliminate it altogether. Corner fences also demonstrated potential for improving the flow in rectangular cross section channels and the combination of corner fences with micro-vortex generators was found to give the greatest benefits. At M=1.3 the combination of corner fences and micro-vanes placed close to the diffuser entry could prevent separation for a wide range of conditions. At the higher diffuser entry Mach number the benefits of flow control were less significant although a reduction of separation size and an improved pressure recovery was observed. It is thought that micro-vortex generators can have significant flow control potential if they are placed close to the expected separation onset and when the adverse pressure gradient is not too far above the incipient separation level. The significant beneficial effects of corner fences warrant a more comprehensive further investigation. It is thought that the control methods suggested here are capable of reducing the bleed requirement on an inlet, which could provide significant performance advantages.

The performance of ejectors driven by sinusoidally unsteady jets

The operation of ejectors driven by a low-speed, sinusoidally unsteady jet has been studied. The thrust augmentation is shown to be highly dependent on the non-dimensional frequency of the driver jet, but independent of its Mach Number. Convective rather than acoustically propagated phenomena dominate the ejector flowfield. Unsteady pressure measurements on the internal surfaces of the ejector have enabled convecting ring vortices to be identified. The impingement of a ring vortex on the leading edge of the ejector causes the peak unsteady body force. The non-dimensional diameter of the ejector is shown to be the only geometric variable that affects the optimum non-dimensional frequency for thrust augmentation. An experimentally optimised geometry is presented. An expression relating the mechanical efficiency and thrust augmentation of the ejector is developed, and shown to be crucially dependent on the degree of unsteadiness in the ejector exit plane.