Development of a turning mid turbine frame with embedded design-part II: Unsteady measurements

© 2014 by ASME. This paper, the second of two parts, presents a new setup for the two-stage two-spool facility located at the Institute for Thermal Turbomachinery and Machine Dynamics (ITTM) of Graz University of Technology. The rig was designed to reproduce the flow behavior of a transonic turbine followed by a counter-rotating low pressure stage such as those in high bypass aero-engines. The meridional flow path of the machine is characterized by a diffusing S-shaped duct between the two rotors. The role of wide chord vanes placed into the mid turbine frame is to lead the flow towards the low pressure (LP) rotor with appropriate swirl. Experimental and numerical investigations performed on this setup showed that the wide chord struts induce large wakes and extended secondary flows at the LP inlet flow. Moreover, large deterministic fluctuations of pressure, which may cause noise and blade vibrations, were observed downstream of the LP rotor. In order to minimize secondary vortices and to damp the unsteady interactions, the mid turbine frame was redesigned to locate two zero-lift splitters into each vane passage. While in the first part of the paper the design process of the splitters and the time-averaged flow field were presented, in this second part the measurements performed by means of a fast response probe will support the explanation of the time-resolved field. The discussion will focus on the comparison between the baseline case (without splitters) and the embedded design.

Development of a turning mid turbine frame with embedded design-part I: Design and steady measurements

© 2014 by ASME. The paper presents a new setup for the two-stage two-spool facility located at the Institute for Thermal Turbomachinery and Machine Dynamics (ITTM) of Graz University of Technology. The rig was designed in order to simulate the flow behavior of a transonic turbine followed by a counter-rotating low pressure (LP) stage like the spools of a modern high bypass aeroengine. The meridional flow path of the machine is characterized by a diffusing S-shaped duct between the two rotors. The role of turning struts placed into the mid turbine frame is to lead the flow towards the LP rotor with appropriate swirl. Experimental and numerical investigations performed on the setup over the last years, which were used as baseline for this paper, showed that wide chord vanes induce large wakes and extended secondary flows at the LP rotor inlet flow. Moreover, unsteady interactions between the two turbines were observed downstream of the LP rotor. In order to increase the uniformity and to decrease the unsteady content of the flow at the inlet of the LP rotor, the mid turbine frame was redesigned with two zero-lifting splitters embedded into the strut passage. In this first part of the paper the design process of the splitters and its critical points are presented, while the time-averaged field is discussed by means of five-hole probe measurements and oil flow visualizations. The comparison between the baseline case and the embedded design configuration shows that the new design is able to reduce the flow gradients downstream of the turning struts, providing a more suitable inlet condition for the low pressure rotor. The improvement in the flow field uniformity is also observed downstream of the turbine and it is, consequently, reflected in an enhancement of the LP turbine performance. In the second part of this paper the influence of the embedded design on the time-resolved field is investigated.

Continuous spectrum growth and modal instability in swirling duct flow

We present results on the stability of compressible inviscid swirling flows in an annular duct. Such flows are present in aeroengines, for example in the by-pass duct, and there are also similar flows in many aeroacoustic or aeronautical applications. The linearised Euler equations have a ('critical layer') singularity associated with pure convection of the unsteady disturbance by the mean flow, and we focus our attention on this region of the spectrum. By considering the critical layer singularity, we identify the continuous spectrum of the problem and describe how it contributes to the unsteady field. We find a very generic family of instability modes near to the continuous spectrum, whose eigenvalue wavenumbers form an infinite set and accumulate to a point in the complex plane. We study this accumulation process asymptotically, and find conditions on the flow to support such instabilities. It is also found that the continuous spectrum can cause a new type of instability, leading to algebraic growth with an exponent determined by the mean flow, given in the analysis. The exponent of algebraic growth can be arbitrarily large. Numerical demonstrations of the continuous spectrum instability, and also the modal instabilities are presented.

Electroosmotic flow analysis of a branched U-turn nanofluidic device.

In this paper, we present the analysis of electroosmotic flow in a branched -turn nanofluidic device, which we developed for detection and sorting of single molecules. The device, where the channel depth is only 150 nm, is designed to optically detect fluorescence from a volume as small as 270 attolitres (al) with a common wide-field fluorescent setup. We use distilled water as the liquid, in which we dilute 110 nm fluorescent beads employed as tracer-particles. Quantitative imaging is used to characterize the pathlines and velocity distribution of the electroosmotic flow in the device. Due to the device's complex geometry, the electroosmotic flow cannot be solved analytically. Therefore we use numerical flow simulation to model our device. Our results show that the deviation between measured and simulated data can be explained by the measured Brownian motion of the tracer-particles, which was not incorporated in the simulation.

Evolution of localized blobs of swirling or buoyant fluid with and without an ambient magnetic field.

We investigate the evolution of localized blobs of swirling or buoyant fluid in an infinite, inviscid, electrically conducting fluid. We consider the three cases of a strong imposed magnetic field, a weak imposed magnetic field, and no magnetic field. For a swirling blob in the absence of a magnetic field, we find, in line with others, that the blob bursts radially outward under the action of the centrifugal force, forming a thin annular vortex sheet. A simple model of this process predicts that the vortex sheet thins exponentially fast and that it moves radially outward with constant velocity. These predictions are verified by high-resolution numerical simulations. When an intense magnetic field is applied, this phenomenon is suppressed, with the energy and angular momentum of the blob now diffusing axially along the magnetic field lines, converting the blob into a columnar structure. For modest or weak magnetic fields, there are elements of both types of behavior, with the radial bursting dominating over axial diffusion for weak fields. However, even when the magnetic field is very weak, the flow structure is quite distinct to that of the nonmagnetic case. In particular, a small but finite magnetic field places a lower bound on the thickness of the annular vortex sheet and produces an annulus of counter-rotating fluid that surrounds the vortex core. The behavior of the buoyant blob is similar. In the absence of a magnetic field, it rapidly develops the mushroomlike shape of a thermal, with a thin vortex sheet at the top and sides of the mushroom. Again, a simple model of this process predicts that the vortex sheet at the top of the thermal thins exponentially fast and rises with constant velocity. These predictions are consistent with earlier numerical simulations. Curiously, however, it is shown that the net vertical momentum associated with the blob increases linearly in time, despite the fact that the vertical velocity at the front of the thermal is constant. As with the swirling blob, an imposed magnetic field inhibits the formation of a vortex sheet. A strong magnetic field completely suppresses the phenomenon, replacing it with an axial diffusion of momentum, while a weak magnetic field allows the sheet to form, but places a lower bound on its thickness. The magnetic field does not, however, change the net vertical momentum of the blob, which always increases linearly with time.

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.

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.

Towards jet flow LES of conceptual nozzles for acoustic predictions

Progress in simulating chevron nozzle jet flows using ILES/RANS-ILES approaches and using the Ffowcs Williams and Hawkings (FW-H) surface integral method to predict the radiated far field sound is presented in this paper. With the focus on the realistic chevron geometries, SMC001 and SMC006, coarse and fine meshes are generated in the range of 3∼13 million mesh cells. Throughout this work, to minimize numerical dissipation introduced by mesh quality issues, the hexahedral cell type is used. Numerical simulations are then carried out with cell-vertex and cell-centered codes. Despite the modest grids, mean velocities and turbulent statistics are found to be in reasonable accord with measurements. Also, far field sound levels predicted by the FW-H post processor are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc.

Finite volume distance field solution applied to medial axis transform

Accurate and efficient computation of the nearest wall distance d (or level set) is important for many areas of computational science/engineering. Differential equation-based distance/ level set algorithms, such as the hyperbolic-natured Eikonal equation, have demonstrated valuable computational efficiency. Here, in the context, as an 'auxiliary' equation to the main flow equations, the Eikonal equation is solved efficiently with two different finite volume approaches (the cell vertex and cell-centered). Application of the distance solution is studied for various geometries. Moreover, a procedure using the differential field to obtain the medial axis transform (MAT) for different geometries is presented. The latter provides a skeleton representation of geometric models that has many useful analysis properties. As an alternative approach to the pure geometric methods (e.g. the Voronoi approach), the current d-MAT procedure bypasses many difficulties that are usually encountered by pure geometric methods, especially in three dimensional space. It is also shown that the d-MAT approach provides the potential to sculpt/control the MAT form for specialized solution purposes. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Rain igestion in axial flow compressors at part speed

New experimental work is reported on the effects of water ingestion on the performance of an axial flow compressor. The background to the work is the effect that heavy rain has on an aeroengine compressor when operating in a "descent idle" mode, i.e., when the compressor is operating at part speed and when the aeromechanical effects of water ingestion are more important than the thermodynamic effects. Most of our existing knowledge in this field comes from whole engine tests. The current work provides the first known results from direct measurements on a stand-alone compressor. The influence of droplet size on path trajectory is considered both computationally and experimentally to show that most rain droplets will collide with the first row of rotor blades. The water on the blades is then centrifuged toward the casing where the normal airflow patterns in the vicinity of the rotor tips are disrupted. The result of this disruption is a reduction in compressor delivery pressure and an increase in the torque required to keep the compressor speed constant. Both effects reduce the efficiency of the machine. The behavior of the water in the blade rows is examined in detail, and simple models are proposed to explain the loss of pressure rise and the increase in torque. The measurements were obtained in a low speed compressor, making it possible to study the mechanical (increase in torque) and aerodynamic (reduction in pressure rise) effects of water ingestion without the added complication of thermodynamic effects. Copyright © 2008 by ASME.

2D PIV measurements in the near field of grid turbulence using stitched fields from multiple cameras

We present measurements of grid turbulence using 2D particle image velocimetry taken immediately downstream from the grid at a Reynolds number of Re M = 16500 where M is the rod spacing. A long field of view of 14M x 4M in the down- and cross-stream directions was achieved by stitching multiple cameras together. Two uniform biplanar grids were selected to have the same M and pressure drop but different rod diameter D and crosssection. A large data set (10 4 vector fields) was obtained to ensure good convergence of second-order statistics. Estimations of the dissipation rate ε of turbulent kinetic energy (TKE) were found to be sensitive to the number of meansquared velocity gradient terms included and not whether the turbulence was assumed to adhere to isotropy or axisymmetry. The resolution dependency of different turbulence statistics was assessed with a procedure that does not rely on the dissipation scale η. The streamwise evolution of the TKE components and ε was found to collapse across grids when the rod diameter was included in the normalisation. We argue that this should be the case between all regular grids when the other relevant dimensionless quantities are matched and the flow has become homogeneous across the stream. Two-point space correlation functions at x/M = 1 show evidence of complex wake interactions which exhibit a strong Reynolds number dependence. However, these changes in initial conditions disappear indicating rapid cross-stream homogenisation. On the other hand, isotropy was, as expected, not found to be established by x/M = 12 for any case studied. © Springer-Verlag 2012.