The design of low-pressure turbines (LPTs) must account for the losses generated by the unsteady interaction with the upstream blade row. The estimation of such unsteady wake-induced losses requires the accurate prediction of the incoming wake dynamics and decay. Existing linear turbulence closures (stress–strain relationships), however, do not offer an accurate prediction of the wake mixing. Therefore, machine-learnt, nonlinear turbulence closures (models) have been developed for LPT flows with unsteady inflow conditions using a zonal-based model development approach, with an aim to enhance the wake mixing prediction for unsteady Reynolds-averaged Navier–Stokes calculations. High-fidelity time-averaged and phase-lock averaged data at a realistic isentropic Reynolds number and two reduced frequencies, i.e., with discrete incoming wakes and with wake “fogging,” have been used as reference data for a machine learning algorithm based on gene expression programing to develop models. Models developed via phase-lock averaged data were able to capture the effect of certain prominent physical phenomena in LPTs such as wake–wake interactions, whereas models based on the time-averaged data could not. Correlations with the flow physics lead to a set of models that can effectively enhance the wake mixing prediction across the entire LPT domain for both cases. Based on a newly developed error metric, the developed models have reduced the a priori error over the Boussinesq approximation on average by 45%. This study thus aids blade designers in selecting the appropriate nonlinear closures capable of mimicking the physical mechanisms responsible for loss generation.

Turbofan rotor–stator aeromechanic and aeroacoustics modeling has been traditionally developed by considering separately the aerodynamic and acoustic response of the fan and the stator to inflow nonuniformities. The present paper develops a model for the coupled fan–stator response for a realistic 3D geometry. The coupling mechanism is assumed to be mainly carried by scattered waves bouncing back and forth between the fan and the stator. The relationship between sources and state elements at three regions (inlet (1), in-between fan and stator (2), and exit (3)) is derived in terms of a scattering matrix S. The model is applied to two fan configurations: (1) a fully subsonic fan and (2) a transonic tip speed fan. The scattering matrix terms, up to the 5th blade passing frequency (BPF), are calculated by using CAAT, an Euler-based code. Results show coupling adds about 6 dB to the sound power level (PWL) of the transonic fan configuration but has a small effect for the subsonic fan configuration. The analysis of the propagating modes shows that, for the transonic configuration, the inlet mode at 1BPF propagates in regions 2 and 3 but is cut off in region 1. This reinforces the coupling process by trapping the acoustic mode 1BPF in region 2. Although this trapped energy is mainly due to the 1BPF of the fan wake, the fan scatters this energy into higher order acoustic modes and thus produces redistribution toward higher frequency of the acoustic spectra. Finally, adding a liner in region 2 reduces the energy of 1BPF mode propagating upstream and impinging on the fan. This mitigates the effect of the fan–stator coupling.

In this paper, we investigated the effects of herringbone riblets, a type of bio-inspired micro-scale surface patterns, on the pressure losses and flow turning angles in a linear cascade over a range of low Reynolds numbers from 0.50×10 5 to 1.50×10 5 and at three different angles of attack. Our experiments showed that despite the micro-scale size of herringbone riblets they produced a significant reduction in pressure loss coefficient and a substantial increase in flow turning angle except at the low end of Reynolds numbers tested. In comparison to the baseline case without riblets, the highest reduction in mean pressure loss achieved was 36.4% which was accompanied by an increase of average turning angle of 4.1°. It was observed that the loss reduction is associated with a reduction in the maximum pressure loss at lower angles of attack and a narrower wake zone at higher angles of attack as a result of suppression of flow separation on the blade suction surface. It was also noted that such a significant improvement in performance was accompanied by the appearance of a serrated wake structure in the contours of pressure loss coefficient. The improvement in blade performance was attributed to the secondary flow motion produced by herringbone riblets which energises the boundary layer. Overall, this work has produced convincing experimental evidence that herringbone riblets can be used as passive flow control devices for reducing flow separation in compressors at low Reynolds numbers.

The aerodynamic interaction of upstream and downstream blade rows can have a significant impact on the forced response of the compressor. Previously, the authors carried out the forced response analysis of a three-row stator-rotor-stator (S1-R2-S2) configuration from a 3.5-stage compressor. However, since the stator vane counts in both the stators (S1 and S2) were the same, it was not possible to separate the excitations from both the rows as they excited the rotor at the same frequency. Hence, a new configuration was developed and tested in which the stator 1 blade count was changed to 38 and stator 2 blade count was maintained at 44 in order to study the individual influences of the stator on the embedded rotor. By using this method, the excitations from both rows can be determined, and the excitations can be quantified to determine the row having the maximum influence on the overall forcing. To achieve this, two sets of simulations were carried out. The three-row stator-rotor (S1-R2-S2) simulation was carried out at both the 38EO (engine order) and 44EO crossings at the peak efficiency (PE) operating condition. The two-row stator-rotor analysis (S1-R2) was carried out at the 38EO crossing, and the other two-Row (R2-S2) analyses were carried out at the 44EO crossing. The steady aerodynamics was preserved in both the cases. A study was done to determine the contribution of wave reflections from the stator inlet and exit planes to the forcing function. Two conclusions drawn from this study are as follows: (1) the modal force value decreased after the upstream stator was removed, which proved that wave reflections from this stator were significant and (2) the increase in modal force was in-line with experimental observations.

In this paper, we investigate the coupled interaction between a new short intake design with a modern fan in a high-bypass ratio civil engine, specifically under the off-design condition of high incidence. The interaction is expected to be much more significant than that on a conventional intake. The performance of both the intake-alone and rotor-alone configurations are examined under isolation. Subsequently, a comprehensive understanding on the two-way interaction between intake and fan is presented. This includes the effect of fan on intake angles of attack (AoA) tolerance (FoI) and the effect of circumferential and radial flow distortion induced by the intake on the fan performance (IoF). In the FoI scenario, the rotor effectively redistributes the mass flow at the fan-face. The AoA tolerance of the short-intake design has increased by ≈4 deg when compared with the intake-alone configuration. Dynamic nature of distortion due to shock unsteadiness has been quantified. ST plots and power spectral density (PSD) of pressure fluctuations show the existence of a spectral gap between the shock unsteadiness and blade passing, with almost an order of magnitude difference in the corresponding frequencies. In the IoF scenario, both the “large” (O(360 deg)) and “small” scale distortion (O(10–60 deg)) induced by the intake results in a non-uniform inflow to the rotor. Sector analysis reveals a substantial variation in the local operating condition of the fan as opposed to its steady characteristic. Streamline curvature, upwash, and wake thickening are identified to be the three key factors affecting the fan performance. These underlying mechanisms are discussed in detail to provide further insights into the physical understanding of the fan-intake interaction. In addition to the shock-induced separation on the intake lip, the current study shows that shorter intakes are much more prone to the upwash effect at higher AoA. Insufficient flow straightening along the engine axis is reconfirmed to be one of the limiting factors for the short-intake design.

The Boxprop is a novel, double-bladed, tip-joined propeller for high-speed flight. The concept draws inspiration from the box wing concept and could potentially decrease tip vortex strength compared with conventional propeller blades. Early Boxprop designs experienced significant amounts of blade interference. By performing a wake analysis and quantifying the various losses of the flow, it could be seen that these Boxprop designs produced 45% more swirl than a conventional reference blade. The reason for this was the proximity of the Boxprop blade halves to each other, which prevented the Boxprop from achieving the required aerodynamic loading on the outer parts of the blade. This paper presents an aerodynamic optimization of a 6-bladed Boxprop aiming at maximizing efficiency and thrust at cruise. A geometric parametrization has been adopted which decreases interference by allowing the blade halves to be swept in opposite directions. Compared with an earlier equal-thrust Boxprop design, the optimized design features a 7% percentage point increase in propeller efficiency and a lower amount of swirl and entropy generation. A vortex-like structure has also appeared downstream of the optimized Boxprop, but with two key differences relative to conventional propellers. (1) Its formation differs from a traditional tip vortex and (2) it is 46% weaker than the tip vortex of an optimized 12-bladed conventional propeller.

For the development of the latest generation of axial compressors, it is necessary to enlarge the design space by using advanced aerodynamic processes. This enables a further increase in efficiency and performance. The use of a tandem blade configuration in a transonic compressor row provides the possibility to enlarge the design space. It is necessary to address the design aspects a bit more in detail in order to efficiently apply this blading concept to turbomachinery. Therefore, in the current study, the known design aspects of tandem blading in compressors will be summed up under the consideration of the aerodynamic effects and construction characteristics of a transonic compressor tandem. Based on this knowledge, a new transonic compressor tandem cascade (DLR TTC) with an inflow Mach number of 0.9 is designed using modern numerical methods and a multi-objective optimization process. Three objective functions as well as three operating points are used in the optimization. Furthermore, both tandem blades and their arrangement are parameterized. From the resulting database of 1246 members, a final best member is chosen as the state-of-the-art design for further detailed investigation. The aim of the ensuing experimental and numerical investigation is to answer the question, whether the tandem cascade resulting from the modern design process fulfills the described design aspects and delivers the requested performance and efficiency criteria. The numerical simulations within the study are carried out with the DLR flow solver TRACE. The experiments are performed at the transonic cascade wind tunnel of DLR in Cologne. The inflow Mach number during the tests is 0.9, and the AVDR is adjusted to 1.3 (design value). Wake measurements with a three-hole probe are carried out in order to determine the cascade performance. The experimental results show an increase in losses and a reduction of the cascade deflection by about 2 deg compared to the design concept. Nevertheless, the experimental and numerical results allow a good understanding of the aerodynamic effects. In addition, planar PIV was applied in a single S1 plane located at midspan to capture the velocity field in the wake of blade 1 in order to analyze the wake flow in detail and describe its influence on the cascade deflection and loss behavior. Finally, an outlook will be given on what future tandem compressor research should be focused.

The experimental results reported in this contribution address the time-dependent impact of periodically unsteady wakes on the development of profile and end wall boundary layers and consequently on the secondary flow system. Experimental investigations are conducted on an annular 1.5 stage axial turbine rig at Ruhr-Universität Bochum’s Chair of Thermal Turbomachines and Aeroengines. The object under investigation is a modified T106 profile low-pressure turbine (LPT) stator row at a representative exit flow Reynolds number of 200,000. By making use of an annular geometry instead of a linear cascade, the influence of curvilinear end walls, nonuniform, increasing pitch across the span and radial flow migration can be represented. Incoming wakes are generated by a variable-speed driven rotor equipped with cylindrical bars. Special emphasis is put on the wake-induced recurrent formation, suppression, weakening, and displacement of individual vortices and separated flow regimes. For this, based on a comprehensive set of time-resolved measurement data, the interaction of impinging bar wakes and boundary layer flow and thus separation and its periodic manipulation along the passage end walls and on the blade suction surface are studied within the frequency domain.

Nonlinear turbulence closures were developed that improve the prediction accuracy of wake mixing in low-pressure turbine (LPT) flows. First, Reynolds-averaged Navier–Stokes (RANS) calculations using five linear turbulence closures were performed for the T106A LPT profile at isentropic exit Reynolds numbers 60,000 and 100,000. None of these RANS models were able to accurately reproduce wake loss profiles, a crucial parameter in LPT design, from direct numerical simulation (DNS) reference data. However, the recently proposed k v 2 ¯ ω transition model was found to produce the best agreement with DNS data in terms of blade loading and boundary layer behavior and thus was selected as baseline model for turbulence closure development. Analysis of the DNS data revealed that the linear stress–strain coupling constitutes one of the main model form errors. Hence, a gene-expression programming (GEP) based machine-learning technique was applied to the high-fidelity DNS data to train nonlinear explicit algebraic Reynolds stress models (EARSM), using different training regions. The trained models were first assessed in an a priori sense (without running any RANS calculations) and showed much improved alignment of the trained models in the region of training. Additional RANS calculations were then performed using the trained models. Importantly, to assess their robustness, the trained models were tested both on the cases they were trained for and on testing, i.e., previously not seen, cases with different flow features. The developed models improved prediction of the Reynolds stress, turbulent kinetic energy (TKE) production, wake-loss profiles, and wake maturity, across all cases.

Experimental studies have been conducted on a modified T106 low pressure turbine (LPT) profile in an annular 1.5 stage axial turbine rig at the Chair of Thermal Turbomachines and Aeroengines, Ruhr-Universität Bochum. The rig setup allows the highly resolved measurement of unsteady wake–stator flow interaction in both space and time. Incoming wakes are generated by a variable-speed driven rotor equipped with cylindrical bars. In the present paper, an experimental approach to the investigation of unsteady phenomena is proposed. Time-averaged and instantaneous measurement data from 2D flow field traverses at the stator exit are provided for the analysis of the periodically unsteady vortex formation, displacement, and suppression. Additional time-accurate blade pressure data are used to study the relationship between the flow structures downstream of the stator row and the immediate intermittent wake impact on the blades. Bar wake kinematics is also discussed in relation to the observations.

In this work, we present the results of the numerical investigations of periodic wake–secondary flow interaction carried out on a low pressure turbine (LPT) equipped with modified T106-profile blades. The numerical predictions obtained by means of unsteady Reynolds-averaged Navier–Stokes (URANS) simulations using a k-ω-model have been compared with measurements conducted in the same configuration and showed a good agreement. Based on the verified numerical data, the Q-criterion has been employed to characterize the secondary flow structures and accurately identify their origin. An analysis of the fundamental wake kinematics and the unsteady vortex migration revealed dominant interaction mechanisms such as the circumferential fluctuation of the pressure side horseshoe vortex (HSV) and its direct interaction with the passage vortex (PV) and the concentrated shed vortex (CSV). Finally, a correlation with the total pressure loss coefficient is provided and a link to the incoming wake structures is given.

Ultra-high bypass ratio (UHBR) engines are designed as compact as possible and are characterized by a short asymmetric air inlet and heterogeneous outlet guide vanes (OGVs). The flow close to the fan is therefore circumferentially nonuniform (or distorted) and the resulting noise might be impacted. This is studied here at take-off conditions by means of a simulation of the unsteady Reynolds-averaged Navier–Stokes (URANS) equations of a full-annulus fan stage. The model includes an asymmetric air inlet, a fan, heterogeneous OGVs, and homogeneous inlet guide vanes (IGVs). Direct acoustic predictions are given for both inlet and aft noises. A novel hydrodynamic/acoustic splitting method based on a modal decomposition is developed and is applied for the aft noise analysis. The noise mechanisms that are generally considered (i.e., interaction of fan-blade wakes with OGVs and fan self-noise) are shown to be impacted by the distortion. In addition, new sources caused by the interaction between the stationary distortion and the fan blades appear and contribute to the inlet noise.

Large eddy simulation (LES) is used to explore the boundary layer transition mechanisms in two rectilinear compressor cascades. To reduce numerical dissipation, a novel locally adaptive smoothing (LAS) scheme is added to an unstructured finite volume solver. The performance of a number of subgrid scale (SGS) models is explored. With the first cascade, numerical results at two different freestream turbulence intensities (Ti's), 3.25% and 10%, are compared. At both Ti's, time-averaged skin-friction and pressure coefficient distributions agree well with previous direct numerical simulations (DNS). At Ti = 3.25%, separation-induced transition occurs on the suction surface, while it is bypassed on the pressure surface. The pressure surface transition is dominated by modes originating from the convection of Tollmien–Schlichting waves by Klebanoff streaks. However, they do not resemble a classical bypass transition. Instead, they display characteristics of the “overlap” and “inner” transition modes observed in the previous DNS. At Ti = 10%, classical bypass transition occurs, with Klebanoff streaks incepting turbulent spots. With the second cascade, the influence of unsteady wakes on transition is examined. Wake-amplified Klebanoff streaks were found to instigate turbulent spots, which periodically shorten the suction surface separation bubble. The celerity line corresponding to 70% of the free-stream velocity, which is associated with the convection speed of the amplified Klebanoff streaks, was found to be important.

In low-pressure turbines (LPT) at design point, around 60–70% of losses are generated in the blade boundary layers far from end walls, while the remaining 30–40% is controlled by the interaction of the blade profile with the end-wall boundary layer. Increasing attention is devoted to these flow regions in industrial design processes. This paper discusses the end-wall flow characteristics of the T106 profile with parallel end walls at realistic LPT conditions, as described in the experimental setup of Duden, A., and Fottner, L., 1997, “Influence of Taper, Reynolds Number and Mach Number on the Secondary Flow Field of a Highly Loaded Turbine Cascade,” Proc. Inst. Mech. Eng., Part A, 211(4), pp.309–320. Calculations are carried out by both Reynolds-averaged Navier–Stokes (RANS), due to its continuing role as the design verification workhorse, and highly resolved large eddy simulation (LES). Part II of this paper focuses on the loss generation associated with the secondary end-wall vortices. Entropy generation and the consequent stagnation pressure losses are analyzed following the aerodynamic investigation carried out in the companion paper (GT2018-76233). The ability of classical turbulence models generally used in RANS to discern the loss contributions of the different vortical structures is discussed in detail and the attainable degree of accuracy is scrutinized with the help of LES and the available test data. The purpose is to identify the flow features that require further modeling efforts in order to improve RANS/unsteady RANS (URANS) approaches and make them able to support the design of the next generation of LPTs.

The present paper describes the application of proper orthogonal decomposition (POD) to large eddy simulation (LES) of the T106A low-pressure-turbine profile with unsteady incoming wakes at two different flow conditions. Conventional data analysis applied to time averaged or phase-locked averaged flow fields is not always able to identify and quantify the different sources of losses in the unsteady flow field as they are able to isolate only the deterministic contribution. A newly developed procedure allows such identification of the unsteady loss contribution due to the migration of the incoming wakes, as well as to construct reduced order models that are able to highlight unsteady losses due to larger and/or smaller flow structures carried by the wakes in the different parts of the blade boundary layers. This enables a designer to identify the dominant modes (i.e., phenomena) responsible for loss, the associated generation mechanism, their dynamics, and spatial location. The procedure applied to the two cases shows that losses in the fore part of the blade suction side are basically unaffected by the flow unsteadiness, irrespective of the reduced frequency and the flow coefficient. On the other hand, in the rear part of the suction side, the unsteadiness contributes to losses prevalently due to the finer scale (higher order POD modes) embedded into the bulk of the incoming wake. The main difference between the two cases has been identified by the losses produced in the core flow region, where both the largest scale structures and the finer ones produces turbulence during migration. The decomposition into POD modes allows the quantification of this latter extra losses generated in the core flow region, providing further inputs to the designers for future optimization strategies.

Accurate estimation of slip factor is of paramount importance to ensure centrifugal compressor work input is adequately predicted during the preliminary design process. However, variations in the flow field at impeller exit in both the pitchwise and spanwise directions complicate the evaluation procedure considerably. With the increasing implementation of engine downsizing technologies in the automotive sector, achieving a wide operating range has become a factor of prime importance for centrifugal compressors used in automotive turbocharging applications. As a result of the design features required to achieve this aim, modern impeller geometries have been shown to exhibit an approximately parabolic variation in slip factor across their respective operating maps. By comparison, traditional slip correlations typically exhibit a constant, or at best monotonic, relationship between slip factor and impeller exit flow coefficient. It is this lack of modeling fidelity which the current work seeks to address. In order to tackle these shortcomings, it is proposed that the impeller exit flow should be considered as being made up of three distinct regions: a region of recirculation next to the shroud providing aerodynamic blockage to the stage active flow, and a pitchwise subdivision of the active flow region into jet and wake components. It is illustrated that this hybrid approach in considering both spanwise and pitchwise stratification of the flow permits a better representation of slip factor to be achieved across the operating map. The factors influencing the relative extent of each of these three distinct regions of flow are numerous, requiring detailed investigations to successfully understand their sources and to characterize their extent. A combination of 3D computational fluid dynamics (CFD) data and gas stand test data for six automotive turbocharger compressor stages was employed to achieve this aim. Through application of the extensive interstage static pressure data gathered during gas stand testing at Queen's University Belfast, the results from the 3D CFD models were validated, thus permitting a more in-depth evaluation of the flow field in terms of locations and parameters that could not easily be measured under gas stand test conditions. Building on previous knowledge gained about the variation in shroud side recirculation with geometry and operating condition, the characteristic jet/wake flow structure emanating from the active flow region of the impeller was represented in terms of area and mass flow components. This knowledge allowed individual slip factor values for the jet and wake to be calculated and combined to give an accurate passage average value which exhibited the distinctive nonlinear variation in slip across the operating map which is frequently absent from existing modeling methods. Fundamental considerations of the flow phenomena in each region provided explanation of the results and permitted a modeling approach to be derived to replicate the trends observed in both the experimental data and the CFD simulations.

The cooling performance of sweeping jet film cooling was studied on a turbine vane suction surface in a low-speed linear cascade wind tunnel. The sweeping jet holes consist of fluidic oscillators with an aspect ratio (AR) of unity and a hole spacing of P d /D = 6. Infrared (IR) thermography was used to estimate the adiabatic film effectiveness at several blowing ratios and two different freestream turbulence levels (Tu = 0.3% and 6.1%). Convective heat transfer coefficient was measured by a transient IR technique, and the net heat flux benefit was calculated. The total pressure loss due to sweeping jet film cooling was characterized by traversing a total pressure probe at the exit plane of the cascade. Tests were performed with a baseline shaped hole (SH) (777-shaped hole) for comparison. The sweeping jet hole showed higher adiabatic film effectiveness than the 777-shaped hole in the near hole region. Although the unsteady sweeping action of the jet augments heat transfer, the net positive cooling benefit is higher for sweeping jet holes compared to 777 hole at particular flow conditions. The total pressure loss measurement showed a 12% increase in total pressure loss at a blowing ratio of M = 1.5 for sweeping jet hole, while 777-shaped hole showed a 8% total pressure loss increase at the corresponding blowing ratio.

The influences of blade loading profile on wake convection and wake/wake interaction were studied in two different blade designs for high-pressure (HP) turbines (front-loaded (FL) and aft-loaded (AL)), installed in linear cascades. A high-speed moving bar (HSMB) apparatus replicated wake shedding, and a closed loop wind tunnel produced engine-relevant Mach numbers (Ma = 0.7) and Reynolds numbers (Re = 3 × 10 5 ). The FL blades had approximately 10% greater total pressure loss when operated with unsteady wake passage. Phase conditioned particle image velocimetry (PIV) measurements were made in the aft portion of the blade channel and downstream of the blade trailing edge. The turbulence kinetic energy (TKE) in the wake was approximately 30% higher for the FL blades when the wake entered the measurement field-of-view. The pressure field in the upstream region of the FL blade design is believed to induce high magnitude strain rates—leading to increased TKE production—and more aggressively turn and dilate the unmixed wake—leading to increased mixing related losses. The higher TKE for the FL blades largely dissipated, being approximately equal to the AL wake by the time the wake reached the end of the blade passage. The interaction of the convected wake with the wake from the blade trailing edge caused periodic vortex shedding at the second harmonic of the convected wake frequency. This interaction also modulated the strength of the trailing edge wake. However, little difference was found in the modulation amplitudes between different cases due to similar strengths of the convected wakes in this region. The higher wake TKE in the upstream portion of the blade channel for the FL blades, therefore, is expected to be the cause of the higher total pressure loss.

Computational fluid dynamics (CFD) has been widely adopted in the compressor design process, but it remains a challenge to predict the flow details, performance, and stage matching for multistage, high-speed machines accurately. The Reynolds Averaged Navier-Stokes (RANS) simulation with mixing plane for bladerow coupling is still the workhorse in the industry and the unsteady bladerow interaction is discarded. This paper examines these discarded unsteady effects via deterministic fluxes using semi-analytical and unsteady RANS (URANS) calculations. The study starts from a planar duct under periodic perturbations. The study shows that under large perturbations, the mixing plane produces dubious values of flow quantities (e.g., whirl angle). The performance of the mixing plane can be considerably improved by including deterministic fluxes into the mixing plane formulation. This demonstrates the effect of deterministic fluxes at the bladerow interface. Furthermore, the front stages of a 19-blade row compressor are investigated and URANS solutions are compared with RANS mixing plane solutions. The magnitudes of divergence of Reynolds stresses (RS) and deterministic stresses (DS) are compared. The effect of deterministic fluxes is demonstrated on whirl angle and radial profiles of total pressure and so on. The enhanced spanwise mixing due to deterministic fluxes is also observed. The effect of deterministic fluxes is confirmed via the nonlinear harmonic (NLH) method which includes the deterministic fluxes in the mean flow, and the study of multistage compressor shows that unsteady effects, which are quantified by deterministic fluxes, are indispensable to have credible predictions of the flow details and performance of compressor even at its design stage.