Experimental investigation of turbine stator well rim seal, reingestion and interstage seal flows using gas concentration techniques and displacement measurements

Gas turbine engine performance requires effective and reliable internal cooling over the duty cycle of the engine. Life predictions for rotating components subject to the main gas path temperatures are vital. This demands increased precision in the specification of the internal air system flows which provide turbine stator well cooling and sealing. This in turn requires detailed knowledge of the flow rates through rim seals and interstage labyrinth seals. Knowledge of seal movement and clearances at operating temperatures is of great importance when prescribing these flows. A test facility has been developed at the University of Sussex, incorporating a two stage turbine rated at 400 kW with an individual stage pressure ratio of 1.7:1. The mechanical design of the test facility allows internal cooling geometry to be rapidly re-configured, while cooling flow rates of between 0.71 CW, ENT and 1.46 C W, ENT, may be set to allow ingress or egress dominated cavity flows. The main annulus and cavity conditions correspond to in cavity rotational Reynolds numbers of 1.71×106< Reφ <1.93×106. Displacement sensors have been used to establish hot running seal clearances over a range of stator well flow conditions, allowing realistic flow rates to be calculated. Additionally, gas seeding techniques have been developed, where stator well and main annulus flow interactions are evaluated by measuring changes in gas concentration. Experiments have been performed which allow rim seal and re-ingestion flows to be quantified. It will be shown that this work develops the measurement of stator well cooling flows and provides data suitable for the validation of improved thermo-mechanical and CFD codes, beneficial to the engine design process. Copyright © 2011 by Rolls-Royce plc.

Non-reflective boundary conditions for non-uniform mean velocity inlet and outlet flows

The numerical solution of problems in unbounded physical space requires a truncation of the computational domain to a reasonable size. As a result, the conditions on the artificial boundaries are generally unknown. Assumptions like constant pressure or velocities are only valid in the far field and lead to spurious reflections if applied on the boundaries of the truncated domain. A number of attempts have been made over the past decades to design conditions that prevent such reflections. One approach is based on characteristics. The standard analysis assumes a spatially uniform mean flow field but this is often impractical. In the present paper we show how to extend the formulation to the more general case of a non-uniform mean velocity field. A number of test cases are provided and our results compare favourably with other boundary conditions. In principle the present approach can be extended to include non-uniformities in all variables.

Transient flows in enclosures: A generalised approach for modelling the effects of geometry, heat gains and wind

Transient flows in a confined ventilated space induced by a buoyancy source of time-varying strength and an external wind are examined. The space considered has varying cross-sectional area with height. A generalised theoretical model is proposed to investigate the flow dynamics following the activation of an external wind and an internal source of buoyancy. To investigate the effect of geometry, we vary the angle of the wall inclination of a particular geometry in which a point source of constant buoyancy is activated in the absence of wind. Counter-intuitively the ventilation is worse and lower airflow rates are established for geometries of increasing cross-sectional areas with height. We investigate the effect of the source buoyancy strength by comparing two cases: (1) when the buoyancy input is constant and (2) when the buoyancy input gradually increases over time so that after a finite time the total buoyancy inputs for (1) and (2) are identical. The rate at which the source heat gains are introduced has a significant role on the flow behaviour as we find that, in case (2), a warmer layer and a more pronounced overshoot are obtained than in case (1). The effect of assisting and opposing wind on the transient ventilation of an enclosure of constant cross-sectional area with height and constant heat gains is examined. A Froude number Fr is used to define the relative strengths of the buoyancy-induced and wind-induced velocities and five different transient states and their associated critical Fr are identified. © 2010 Elsevier Ltd.

Ventilation effectiveness measures based on heat removal: Part 2. Application to natural ventilation flows

The three effectiveness measures based on the ability of a flow to flush buoyancy from a ventilated space proposed by Coffey and Hunt [Ventilation effectiveness measures based on heat removal-part 1. Definitions. Building and Environment, in press, doi:10.1016/j.buildenv.2006.03.016.] are applied to assess and compare two fundamental natural ventilation flows. We focus on the limiting cases of passive displacement and passive mixing ventilation flows during transient conditions. These transient flows occur when, for example, heat is purged from a building at night. Whilst it is widely recognised that mixing flows are less efficient at purging heat than displacement flows, our results indicate that, when a particular zone of a room is considered, displacement ventilation can result in lower effectiveness than mixing ventilation. When a room is considered as a whole, displacement ventilation yields higher effectiveness than mixing ventilation and we quantify these differences in terms of the geometry of the space and opening area. The proposed theoretical predictions are compared with effectiveness deduced from measurements made during laboratory experiments and show good agreement. © 2006 Elsevier Ltd. All rights reserved.

Strategies for modeling turbulent flows in electronics

Hybrid methods based on the Reynolds Averaged Navier Stokes (RANS) equations and the Large Eddy Simulation (LES) formulation are investigated to try and improve the accuracy of heat transfer and surface temperature predictions for electronics systems and components. Two relatively low Reynolds number flows are studied using hybrid RANS-LES, RANS-Implicit-LES (RANS-ILES) and non-linear LES models. Predictions using these methods are in good agreement with each other, even using different grid resolutions. © 2008 IEEE.

Mixed nonlinear LES for DES suitable flows

It is shown that a new mixed nonlinear/eddy viscosity LES model reproduces profiles better than a number of competing nonlinear and mixed models for plane channel flow. The objective is an LES method that produces a fully resolved turbulent boundary layer and could be applied to a variety of aerospace problems that are currently studied with RANS, RANS-LES, or DES methods that lack a true turbulent boundary layer. There are two components to the new model. One an eddy viscosity based upon the advected subgrid scale energy and a relatively small coefficient. Second, filtered nonlinear terms based upon the Leray regularization. Coefficients for the eddy viscosity and nonlinear terms come from LES tests in decaying, isotropic turbulence. Using these coefficients, the velocity profile matches measurements data at Reτ ≈ 1000 exactly. Profiles of the components of kinetic energy have the same shape as in the experiment, but the magnitudes differ by about 25%. None of the competing LES gets the shape correct. This method does not require extra operations at the transition between the boundary layer and the interior flow.

Predictive techniques for turbulent oscillatory flows

The prediction of turbulent oscillatory flow at around transitional Reynolds numbers is considered for an idealized electronics system. To assess the accuracy of turbulence models, comparison is made with measurements. A stochastic procedure is used to recover instantaneous velocity time traces from predictions. This procedure enables more direct comparison with turbulence intensity measurements which have not been filtered to remove the oscillatory flow component. Normal wall distances, required in some turbulence models, are evaluated using a modified Poisson equation based technique. A range of zero, one and two equation turbulence models are tested, including zonal and a non-linear eddy viscosity models. The non-linear and zonal models showed potential for accuracy improvements.

Computation of transient and steady turbulent isothermal flows with an impinging rectangular jet

The computation of both transient and steady turbulent incompressible isothermal flows is studied. The flow is very complex, having streamline curvature, large vortex structures and stagnation resulting from an impinging rectangular jet. For transient computations, the standard k-ε model is adopted. For steady flows, the k-ε, high and low Reynolds number k-l and mixing length models are tried. Zonal approaches combining the above turbulence models are also investigated. None of the models are found to give satisfactory agreement with velocity measurements.

Multi-objective optimization for multiphase orbital rendezvous missions

A multi-objective optimization approach was proposed for multiphase orbital rendezvous missions and validated by application to a representative numerical problem. By comparing the Pareto fronts obtained using the proposed method, the relationships between the three objectives considered are revealed, and the influences of other mission parameters, such as the sensors' field of view, can also be analyzed effectively. For multiphase orbital rendezvous missions, the tradeoff relationships between the total velocity increment and the trajectory robustness index as well as between the total velocity increment and the time of flight are obvious and clear. However, the tradeoff relationship between the time of flight and the trajectory robustness index is weak, especially for the four- and five-phase missions examined. The proposed approach could be used to reorganize a stable rendezvous profile for an engineering rendezvous mission, when there is a failure that prevents the completion of the nominal mission.

Turbine blade tip heat transfer in low speed and high speed flows

In this paper, high and low speed tip flows are investigated for a high-pressure turbine blade. Previous experimental data are used to validate a CFD code, which is then used to study the tip heat transfer in high and low speed cascades. The results show that at engine representative Mach numbers the tip flow is predominantly transonic. Thus, compared to the low speed tip flow, the heat transfer is affected by reductions in both the heat transfer coefficient and the recovery temperature. The high Mach numbers in the tip region (M>1.5) lead to large local variations in recovery temperature. Significant changes in the heat transfer coefficient are also observed. These are due to changes in the structure of the tip flow at high speed. At high speeds, the pressure side corner separation bubble reattachment occurs through supersonic acceleration which halves the length of the bubble when the tip gap exit Mach number is increased from 0.1 to 1.0. In addition, shock/boundary-layer interactions within the tip gap lead to large changes in the tip boundary-layer thickness. These effects give rise to significant differences in the heat-transfer coefficient within the tip region compared to the low-speed tip flow. Compared to the low speed tip flow, the high speed tip flow is much less dominated by turbulent dissipation and is thus less sensitive to the choice of turbulence model. These results clearly demonstrate that blade tip heat transfer is a strong function of Mach number, an important implication when considering the use of low speed experimental testing and associated CFD validation in engine blade tip design. Copyright © 2009 by ASME.