Application of a laminar flamelet model to confined explosion hazards

The majority of computational studies of confined explosion hazards apply simple and inaccurate combustion models, requiring ad hoc corrections to obtain realistic flame shapes and often predicting an order of magnitude error in the overpressures. This work describes the application of a laminar flamelet model to a series of two-dimensional test cases. The model is computationally efficient applying an algebraic expression to calculate the flame surface area, an empirical correlation for the laminar flame speed and a novel unstructured, solution adaptive numerical grid system which allows important features of the solution to be resolved close to the flame. Accurate flame shapes are predicted, the correct burning rate is predicted near the walls, and an improvement in the predicted overpressures is obtained. However, in these fully turbulent calculations the overpressures are still too high and the flame arrival times too low, indicating the need for a model for the early laminar burning phase. Due to the computational expense, it is unrealistic to model a laminar flame in the complex geometries involved and therefore a pragmatic approach is employed which constrains the flame to propagate at the laminar flame speed. Transition to turbulent burning occurs at a specified turbulent Reynolds number. With the laminar phase model included, the predicted flame arrival times increase significantly, but are still too low. However, this has no significant effect on the overpressures, which are predicted accurately for a baffled channel test case where rapid transition occurs once the flame reaches the first pair of baffles. In a channel with obstacles on the centreline, transition is more gradual and the accuracy of the predicted overpressures is reduced. However, although the accuracy is still less than desirable in some cases, it is much better than the order of magnitude error previously expected.

Wake, shock and potential field interactions in a 1.5 stage turbine: Part I: Vane-rotor and rotor-vane interaction

The composition of the time-resolved surface pressure field around a high-pressure rotor blade caused by the presence of neighboring blade rows was studied, with the individual effects of wake, shock and potential field interaction being determined. Two test geometries were considered: first, a high-pressure turbine stage coupled with a swan-necked diffuser exit duct; secondly, the same high-pressure stage but with a vane located in the downstream duct. Both tests were carried out at engine-representative Mach and Reynolds numbers. By comparing the results to time-resolved computational predictions of the flowfield, the accuracy with which the computation predicts blade interaction was determined. It was found that in addition to upstream vane-rotor and rotor-downstream vane interactions, a new interaction mechanism was found resulting from the interaction between the downstream vane's potential field and the upstream vane's trailing edge potential field and shock.

Wake, shock and potential field interactions in a 1.5 stage turbine: Part II: Vane-vane interaction and discussion of results

The composition of the time-resolved surface pressure field around a high-pressure rotor blade caused by the presence of neighboring blade rows was studied, with the individual effects of wake, shock and potential field interaction being determined. Two test geometries were considered: first, a high-pressure turbine stage coupled with a swan-necked diffuser exit duct; secondly, the same high-pressure stage but with a vane located in the downstream duct. Both tests were carried out at engine-representative Mach and Reynolds numbers. By comparing the results to time-resolved computational predictions of the flowfield, the accuracy with which the computation predicts blade interaction was determined. Evidence was obtained that for a large downstream vane, the flow conditions in the rotor passage, at any instant in time, are close to being steady state.

The effect of wake passing on a flow separation in a low pressure turbine cascade

This paper describes an investigation into the effect that passing wakes have on a separation bubble that exists on the pressure surface and near the leading edge of a low pressure turbine blade. Previous experimental studies have shown that the behaviour of this separation is strongly incidence dependent and that it responds to its disturbance environment. The results presented in this paper examine the effect of wake passing in greater detail. Two dimensional, Reynolds averaged, numerical predictions are first used to examine qualitatively the unsteady interaction between the wakes and the separation bubble. The separation is predicted to consist of spanwise vortices whose development is in phase with the wake passing. However, comparison with experiments shows that the numerical predictions exaggerate the coherence of these vortices and also overpredict the time-averaged length of the separation. Nonetheless, experiments strongly suggest that the predicted phase locking of the vortices in the separation onto the wake passing is physical.

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 an upstream turbine on a low-aspect ratio vane

The interaction between a high-pressure rotor and a downstream vane is dominated by vortex-blade interaction. Each rotor blade passing period two co-rotating vortex pairs, the tip-leakage and upper passage vortex and the lower passage and trailing shed vortex, impinge on, and are cut by, the vane leading edge. In addition to the streamwise vortex the tip-leakage flow also contains a large velocity deficit. This causes the interaction of the tip-leakage flow with a downstream vane to differ from typical vortex blade interaction. This paper investigates the effect these interaction mechanisms have on a downstream vane. The test geometry considered was a low aspect ratio second stage vane located within a S-shaped diffuser with large radius change mounted downstream of a shroudless high-pressure turbine stage. Experimental measurements were conducted at engine-representative Mach and Reynolds numbers, and data was acquired using a fast-response aerodynamic probe upstream and downstream of the vane. Time-resolved numerical simulations were undertaken with and without a rotor tip gap in order to investigate the relative magnitude of the interaction mechanisms. The presence of the upstream stage is shown to significantly change the structure of the secondary flow in the vane and to cause a small drop in its performance.

Multi-objective design and optimisation of bypass outlet-guide vanes

This paper describes the development of an automated design optimization system that makes use of a high fidelity Reynolds-Averaged CFD analysis procedure to minimize the fan forcing and fan BOGV (bypass outlet guide vane) losses simultaneously taking into the account the down-stream pylon and RDF (radial drive fairing) distortions. The design space consists of the OGV's stagger angle, trailing-edge recambering, axial and circumferential positions leading to a variable pitch optimum design. An advanced optimization system called SOFT (Smart Optimisation for Turbomachinery) was used to integrate a number of pre-processor, simulation and in-house grid generation codes and postprocessor programs. A number of multi-objective, multi-point optimiztion were carried out by SOFT on a cluster of workstations and are reported herein.

The effects of a trip wire and unsteadiness on a high speed highly loaded low-pressure turbine blade

This paper presents the effect of a single spanwise 2D wire upon the downstream position of boundary layer transition under steady and unsteady inflow conditions. The study is carried out on a high turning, high-speed, low pressure turbine (LPT) profile designed to take account of the unsteady flow conditions. The experiments were carried out in a transonic cascade wind tunnel to which a rotating bar system had been added. The range of Reynolds and Mach numbers studied includes realistic LPT engine conditions and extends up to the transonic regime. Losses are measured to quantify the influence of the roughness with and without wake passing. Time resolved measurements such as hot wire boundary layer surveys and surface unsteady pressure are used to explain the state of the boundary layer. The results suggest that the effect of roughness on boundary layer transition is a stability governed phenomena, even at high Mach numbers. The combination of the effect of the roughness elements with the inviscid Kelvin-Helmholtz instability responsible for the rolling up of the separated shear layer (Stieger [1]) is also examined. Wake traverses using pneumatic probes downstream of the cascade reveal that the use of roughness elements reduces the profile losses up to exit Mach numbers of 0.8. This occurs with both steady and unsteady inflow conditions.

Heat transfer and aerodynamics of turbine blade tips in a linear cascade

Local measurements of the heat transfer coefficient and pressure coefficient were conducted on the tip and near tip region of a generic turbine blade in a five-blade linear cascade. Two tip clearance gaps were used: 1.6% and 2.8% chord. Data was obtained at a Reynolds number of 2.3 × 10 5 based on exit velocity and chord. Three different tip geometries were investigated: a flat (plain) tip, a suction-side squealer, and a cavity squealer. The experiments reveal that the flow through the plain gap is dominated by flow separation at the pressure-side edge and that the highest levels of heat transfer are located where the flow reattaches on the tip surface. High heat transfer is also measured at locations where the tip-leakage vortex has impinged onto the suction surface of the aerofoil. The experiments are supported by flow visualisation computed using the CFX CFD code which has provided insight into the fluid dynamics within the gap. The suction-side and cavity squealers are shown to reduce the heat transfer in the gap but high levels of heat transfer are associated with locations of impingement, identified using the flow visualisation and aerodynamic data. Film cooling is introduced on the plain tip at locations near the pressure-side edge within the separated region and a net heat flux reduction analysis is used to quantify the performance of the successful cooling design. copyright © 2005 by ASME.

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.

LES and hybrid LES/RANS simulations for conjugate heat transfer over a matrix of cubes

Large Eddy Simulation (LES) and a novel k -l based hybrid LES/RANS approach have been applied to simulate a conjugate heat transfer problem involving flow over a matrix of surface mounted cubes. In order to assess the capability and reliability of the newly developed k -l based hybrid LES/RANS, numerical results are compared with new LES and existing RANS results. Comparisons include mean velocity profiles, Reynolds stresses and conjugate heat transfer. As well as for hybrid LES/RANS validation purposes, the LES results are used to gain insights into the complex flow physics and heat transfer mechanisms. Numerical simulations show that the hybrid LES/RANS approach is effective. Mean and instantaneous fluid temperatures adjacent to the cube surface are found to strongly correlate with flow structure. Although the LES captures more mean velocity field complexities, broadly time averaged wake temperature fields are found similar for the LES and hybrid LES/RANS. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Vortex generators near shock/ boundary layer interactions

Experiments have been performed in a blowdown supersonic wind tunnel to investigate the effect of arrays of sub-boundary layer vortex generators placed upstream of a normal shock/ boundary layer interaction. The investigation makes use of a recovery shock wave and the naturally grown turbulent boundary layer on the wind tunnel floor. Experiments were performed at Mach numbers of 1.5 and 1.3 and a freestream Reynolds number of 28 × 106. Two types of vortex generators were investigated - wedge-shaped and arrays of counter-rotating vanes. It was found that at Mach 1.5 the vane-type VGs eliminated and the wedge-type VGs greatly reduced the separation bubble under the shock. When placed in the supersonic part of the flow both VGs caused a wave pattern consisting of a shock, re-expansion and shock. The re-expansion and double shocks are undesirable features since they equate to increased total pressure losses and hence increased -wave drag. Furthermore there are indications that the vortex intensity is reduced by the normal shock/ boundary layer interaction. When the shock was located directly over the VGs there was no re-expansion present, but the 'damping' effect of the shock on the vortex persisted. It appears that the vortices produced by the wedge-shaped VGs lift off the surface more rapidly. Similar results were observed at Mach 1.3, where the flow was unseparated.

Secondary flow near an undulatory surface induced by wall blocking effect

Prandtl's secondary mean motions of the second kind near an undulating surface were explained in terms of turbulent blocking effect and kinematic boundary conditions at the surface, and its order of magnitude was estimated. Isotropic turbulence is distorted by the undulating surface of wavelength λ and amplitude h with a low slope, so that h « λ. The prime mechanism for generating the mean flow is that the far-field Isotropic turbulence is distorted by the non-local blocking effect of the surface to become anisotropic axisymmetric turbulence near the surface with principal axis that is not aligned with the local curvature of the undulation. Then the local analysis can be applied and the mechanism is similar to the mean flow generation mechanism for homogeneous axisymmetric turbulence over a planer surface, i.e. gradients of the Reynolds stress caused by the turbulent blocking effect generate the mean motions. The results from this simple analysis are consistent with previous exact analysis in which the effects of curvature are strictly taken into account. The results also qualitatively agree with flow visualization over an undulating surface in a mixing-box.

Adjoint linearised Euler solver in the frequency domain for jet noise modelling

The paper is devoted to extending the new efficient frequency-domain method of adjoint Green's function calculation to curvilinear multi-block RANS domains for middle and farfield sound computations. Numerical details of the method such as grids, boundary conditions and convergence acceleration are discussed. Two acoustic source models are considered in conjunction with the method and acoustic modelling results are presented for a benchmark low-Reynolds-number jet case.

The prediction of hot streak migration in a high-pressure turbine

Accurate predictions of combustor hot streak migration enable the turbine designer to identify high-temperature regions that can limit component life. It is therefore important that these predictions are achieved within the short time scales of a design process. This article compares temperature measurements of a circular hot streak through a turning duct and a research turbine with predictions using a three-dimensional Reynolds-averaged Navier-Stokes solver. It was found that the mixing length turbulence model did not predict the hot streak dissipation accurately. However, implementation of a very simple model of the free stream turbulence (FST) significantly improved the exit temperature predictions on both the duct and research turbine. One advantage of the simple FST model described over more complex alternatives is that no additional equations are solved. This makes the method attractive for design purposes, as it is not associated with any increase in computational time.