Abstract

High-fidelity simulations are used to conduct controlled numerical experiments to investigate the effect of periodically incoming wakes on profile and three-dimensional loss mechanisms. The present work considers the MTU-T161 cascade with spanwise diverging end-walls, representative of a high-lift, low-pressure turbine blade. All simulations are carried out at engine-relevant conditions, with exit Reynolds number of 90,000 and exit Mach number of 0.6. Upstream moving bars are used to generate incoming wakes which impinge on the blade and potentially alter its aerodynamic performance. Unlike previous studies, the incoming wakes are subjected to an additional axial pressure gradient when convecting through the passage, due to the divergence of the spanwise end-walls. The evolving secondary vortex systems around the bars periodically disturb the freestream end-wall boundary layer facing the blade leading edge. This ultimately influences the end-wall related losses downstream of the blade and governs the overall aerodynamic performance of the blade. Following validation against available experimental data, a systematic variation of flow coefficient and reduced frequency extends the parametric space studied to encompass engine-realistic operating conditions. The high-fidelity simulations reveal the impact of incoming wakes on blade boundary layer losses and wake-induced losses both at the mid-span and within the end-wall regions. Furthermore, by decomposing the total loss generation, the data-rich results shed light on the underlying physical mechanisms driving unsteady losses when applied to phase- and time-averaged flow fields. Secondary losses incurred in the end-wall region show little sensitivity toward unsteadiness associated with incoming wakes and are rather prone to the turbulence levels in the passage. On the other hand, profile losses show high dependency on bar wakes in the absence of wake fogging. While profile losses can be minimized by certain combinations of flow coefficients and reduced frequencies, they remain the dominant source of unsteady loss generation.

1 Introduction

In order to maximize the potential of modern geared turbofans, the low-pressure turbine (LPT) commonly features higher rotational speeds leading to a reduction in core engine size. In such scenarios, secondary flow effects can have increased impact on blade aerodynamics and overall performance. Furthermore, unsteadiness due to periodically incoming wakes, representative of LPT stage cases, increase the complexity for accurate time-averaged and time-resolved loss predictions. Loss generation mechanisms must be understood and well-predicted to better guide the development of more efficient designs and configurations in LPT design cycles. Generally speaking, the understanding of secondary flow phenomena in LPTs is well-documented for steady-state flow conditions, e.g., by Sieverding [1] and Goldstein and Spores [2], and their relation to loss generation is reported, e.g., by Langston [3]. The unsteady nature of blade-to-blade interactions and mechanisms driving unsteady losses, on the other hand, is challenging to predict. Most commonly, low-fidelity computational fluid dynamics (CFD) analyses and laboratory experiments have been standard industrial tools to yield detailed physical understanding. On the one hand, numerical studies are predominately carried out using unsteady Reynolds-averaged Navier–Stokes (URANS) analyses and compared to experimental campaigns, e.g., in Ref. [4]. The accuracy of URANS, driven by underlying turbulence modeling assumptions, however, is notoriously challenged by the complex flow conditions within LPTs, including strong pressure gradients, separated flow, and intermittent laminar to turbulent transition [5]. Experimental cascade testing, on the other hand, is time-consuming, costly, and notoriously difficult to instrument. Even though unsteady rotor-like wakes can be generated using upstream moving bars, engine-relevant flow coefficients and reduced frequencies often cannot be reproduced due to mechanical constraints [6,7].

With recent advancements in massively parallel and accelerated high-performance computing facilities, high-fidelity simulations are a promising additional tool for the detailed analysis of turbomachinery flows to URANS and experiments. They are not bound by mechanical constraints and, depending on resolution, only model the dissipation region while the dynamically relevant fluid flow scales are fully resolved. Engine-relevant operating conditions can now be simulated. Highly resolved, large eddy simulations (LES) and direct numerical simulations (DNS) have been validated against experimental data in the past and have been able to shed light on the physical understanding of unsteady LPT flows [8,–11]. Michelassi et al. [9,12] conducted DNS and LES of LPTs with upstream moving bars, discussing their effect on aerodynamic performance of the highly loaded T106 LPT blade. The studies comprised a set of cases with varying flow coefficient and reduced frequencies showing that the wake-to-wake width could be correlated to loss generation in the mid-span wake. To minimize computational requirements, spanwise periodic boundary conditions were employed. Taking spanwise end-wall effects into account, Cui et al. [13] and Koschichow et al. [14] studied the same cascade using low- and high-fidelity methods. Both cases reported time-dependent variations of mid-span profile losses due to intermittent separation of the suction side boundary layer, while deviations in secondary losses remained small compared to profile losses. Robison and Gross [15] concluded a stronger dependence if incoming wakes are generated by upstream rotor blades instead of cylindrical bars at a Mach number of 0.1. By contrast, Lopes et al. [16] reported a 10% increase of secondary losses due to moving bars at engine-relevant Mach numbers.

However, in all of these cases, the spanwise end-walls were perfectly parallel. Modern LPT configurations, on the other hand, feature divergent gas paths. The effect of tapering can lead to a significant increase in secondary losses [17] over a range of Mach number and Reynolds number and must be taken into account for modern LPTs. On top of that, axial-velocity-density-ratio effects further lead to an additional pressure gradient in the mid-span region which cannot be simulated using spanwise periodic boundary conditions or spanwise parallel end-walls [18].

In summary, the available literature is scarce when engine-like geometries are of interest at engine-relevant conditions subject to varying incoming wakes. The present investigation aims to close this gap, complementing existing knowledge and providing insights into unsteady loss generation mechanisms studied via high-fidelity simulations. The MTU-T161 LPT cascade featuring spanwise divergent end-walls at a tapering angle of 12 deg is chosen as a numerical test rig representative of modern LPT cascades. Its steady flow conditions have been widely discussed using high-fidelity simulations in the past, i.e., in Refs. [1923]. A brief summary of key aerodynamic parameters, numerically determined in the absence of incoming wakes [23], and geometrical parameters is given in Table 1. At engine-relevant conditions of Reynolds number of 90,000, based on isentropic exit velocity and chord length, and Mach number of 0.6, the turbine blade features a closed separation bubble at the mid-span and significant secondary flow structures. This configuration was also considered by Morsbach et al. [22], who compared steady and unsteady inflow conditions by means of upstream moving bars, although only at a single operating point. The time-averaged results revealed that the incoming wakes are able to completely suppress the suction side separation bubble. The total pressure losses downstream of the blade increased over the whole span taking into account the losses due to the mixing behind the moving bars.

Table 1

Key aerodynamic parameters and geometrical properties of the MTU-T161 LPT cascade

AerodynamicValueGeometrical parameterValue
Re2th90,000lpitch/C0.956
Ma2th0.60lspan,LE/C2.467
α141 deglspan,TE/C2.937
α265 deg
Zw1.2
AerodynamicValueGeometrical parameterValue
Re2th90,000lpitch/C0.956
Ma2th0.60lspan,LE/C2.467
α141 deglspan,TE/C2.937
α265 deg
Zw1.2

In the present study, we complement the insights from Morsbach et al. by performing a parametric sweep study. A systematic variation of flow coefficient and reduced frequency is presented via high-fidelity methods comprising a set of eight highly resolved large eddy simulations. The data-rich results establish a database for three-dimensional losses given wake-induced losses to elucidate the driving mechanisms behind periodically incoming, rotor-like wakes impinging on the stator blade.

The paper is structured in the following way. First, the numerical setup and parameter space is outlined. Second, a grid convergence study is presented and validation against experimental data benchmarks the present high-fidelity simulations. Third, focus is put on the mid-span section highlighting the importance of additional pressure gradients due to the spanwise divergent end-walls on the unsteady (profile) losses. Afterwards, end-wall associated losses are considered and phase-averaged results are presented. Lastly, insights and understanding of the data-rich results are summarized.

2 Methodology

2.1 Numerical Setup.

The flow solver of choice is the in-house CFD solver high-performance solver for turbulence and aeroacoustic research (HiPSTAR), which has specifically been developed to undertake highly resolved LES and DNS efficiently on the latest high-performance computing architectures, with a detailed description given in Ref. [8]. In the context of turbomachinery flows, HiPSTAR has been extensively validated for compressor, LPT, and high-pressure turbine blades [5,2426]. Spatial discretization is achieved using a fourth-order accurate, finite-difference scheme. To advance the solution in time, an ultra-low storage, frequency optimized, fourth-order explicit Runge–Kutta method is used. For the present LES, the wall-adaptive, local eddy model is used as the subgrid-scale model [27]. Non-reflecting, zonal characteristic boundary conditions are employed at the outlet [28]. Periodicity is assumed in the pitchwise direction. The mesh configuration is based on a stable overset technique [29] with O-grids around the blades and bars to resolve the boundary layers adequately. These O-grids are embedded into two background H-grids which resolve the broader turbulent flow field as shown in Fig. 1 for multiple bar configurations. Each bar has a diameter of 2 mm to simulate rotor-like wakes. The sliding interface condition between the H-grids is based on the message passing interface halo exchange and subsequent fourth-order accurate Lagrange interpolation in the pitchwise direction. At the domain inlet, Riemann boundary conditions are prescribed with turbulent fluctuations based on a digital filter method [30] and have been calibrated in pre-curser simulations to match the decay of turbulent intensity in mid-span according to experimental data [31]. However, the experimentally determined end-wall boundary layer state cannot be reproduced numerically due to the nature of the test facility and the turbulent grid characteristics used in the measurement campaigns [32]. As the boundary layer state can influence the end-wall loss generation significantly, the simulations in this study use the same base inflow profile as Profile A in Ref. [23], which was found to result in good agreement with the experiments in the end-wall region of the blade wake.

Fig. 1
Mid-span overset grid setup showing sliding bar O-grid embedded into sliding H-grid (left-hand side of Section 1) and non-sliding blade O-grid embedded into non-sliding H-grid (right-hand side of Section 1). Every eighth grid line is shown.
Fig. 1
Mid-span overset grid setup showing sliding bar O-grid embedded into sliding H-grid (left-hand side of Section 1) and non-sliding blade O-grid embedded into non-sliding H-grid (right-hand side of Section 1). Every eighth grid line is shown.
Close modal

2.2 Validation.

First, a grid convergence and validation study is performed by comparing against experimental data of Entlesberger et al. [33] to justify the numerical setup and subsequent simulations. In order to minimize computational costs, the convergence study is based on a single blade with one upstream moving bar. The grid resolution of the steady configurations reported in Ref. [23] are taken as a baseline case. A further refined case, labeled fine, yields DNS resolution in the spanwise direction, necessary to converge the skin friction coefficient on the blade suction side of LPTs [18]. The final grid counts are given in Table 2. For reference purposes, the hours required to simulate one convective time unit tc=C/u¯2is on graphical processing units (GPUs), based on chord length C and isentropic exit velocity u¯2is, is appended alongside the total number of GPUs per simulation (NGPUs). The point spacings around the blade normalized on the local viscous length scale δν=ν¯/uτ¯ are denoted as Δn+ in the wall normal direction, Δs+ in the blade tangential direction, and Δz+ in the spanwise direction. Here, the wall friction velocity uτ=τw¯/ρ¯ is based on the wall shear stress τw=sgn(τw,x¯)τw,x¯2+τw,y¯2 and density ρ¯. Looking at the two top graphs of Fig. 2, Δn+<1.0 is satisfied and Δs+<20.0 ensures high quality streamwise resolution along the pressure side with s/Cax<0.0 and suction side s/Cax>0.0. The third graph outlines the differences between the two grids and reveals that the fine grid maintains DNS resolution of Δz+<10 over the entire blade. In terms of velocity gradients at the wall, the skin friction coefficient cf¯=τw¯/(pt1p2), based on total pressure at station 1, pt1, and static pressure at station 2, p2, shows grid sensitivity in the separation bubble region, where cf¯<0.0, which is expected and shown in the bottom plot. An adequate DNS-like resolution in the spanwise direction is therefore necessary to fully resolve the velocity gradients in the separated region along the blade [18]. The flow coefficient is commonly defined as

Fig. 2
Mid-span non-dimensional point spacing in viscous units and skin friction coefficient plotted over blade path length non-dimensionalized using chord length such that s/C<0 refers to pressure side and s/C>0 refers to suction side
Fig. 2
Mid-span non-dimensional point spacing in viscous units and skin friction coefficient plotted over blade path length non-dimensionalized using chord length such that s/C<0 refers to pressure side and s/C>0 refers to suction side
Close modal
Table 2

Grid convergence parameters for one bar setup

GridNxyNzNxyz/106NGPUsGPUh/tc
Baseline299,856661198.23263.66
Fine299,8561345403.364156.18
GridNxyNzNxyz/106NGPUsGPUh/tc
Baseline299,856661198.23263.66
Fine299,8561345403.364156.18
(1)
Here, u¯1ax denotes the pitchwise area-averaged axial velocity behind the bars at Sec. 1 (see Fig. 1) and ubar denotes the bar speed. However, the experimentally obtained data do not contain information about the axial velocity behind the bars and must therefore be iteratively approximated based on Φexp=uax/ubar using the freestream axial velocity upstream of the bars uax. Furthermore, reproducing the bar-to-blade ratio of the measurements would require simulating six blades which is too costly for a grid convergence study. As a result, the bar passing frequency f=ubar/lpitch,bar and subsequent reduced frequency
(2)
are about 20% higher in the numerical setup.

Lastly, experimentally traversing the bars along the entire span through slots results in a leakage flow at the spanwise end-walls altering the experimental operating conditions. The numerical operating conditions are therefore iteratively approximated via the baseline grid. The calibrated flow coefficient and reduced frequency are determined as Φcalibrateduax/ubar=2.49 and Fred,calibrated=0.20, respectively. The resulting blade loading and total pressure deficit in the wake are shown in Fig. 3. In the top graph, the blade pressure p¯ at the wall is non-dimensionalized using a reference pressure taken upstream of the moving bars in accordance with experimental data. Both grids are able to accurately capture the experimentally obtained data, giving confidence in the chosen numerical setup and subsequent parametric sweeps.

Fig. 3
Mid-span blade pressure distribution and total pressure loss coefficient for the baseline and fine grid simulations. Experimental values are additionally marked as circles.
Fig. 3
Mid-span blade pressure distribution and total pressure loss coefficient for the baseline and fine grid simulations. Experimental values are additionally marked as circles.
Close modal
The total pressure loss coefficient defined in the experimental test campaign follows
(3)
using static outlet pressure pout and the freestream total pressure pt. It is evaluated at 40% axial chord length downstream of the blade in accordance with the experimental setup and identical to Fig. 1 in Sec. 2. The wake peak is predicted more accurately by the baseline grid compared to the fine grid which might be fortuitous. On the other hand, the fine grid resolution agrees better with the wake onset and wake width which is a result of the DNS-like resolution of the separation bubble on the blade suction side. The simulations reported throughout this study are, therefore, using the fine grid resolution. After initial transients have been washed out, statistical quantities have been collected over 20 bar passes ensuring statistical convergence.

2.3 Parameter Space Investigation.

The present study investigates the effect of varying reduced frequency and flow coefficient on the aerodynamic performance of a high-lift LPT blade. Changing the bar speed results in a change of the flow coefficient ϕ. In the simulations reported here, the flow coefficient approximating that of the experiment is labeled as 1u, and is subsequently doubled, 2u, and tripled 3u. For fixed flow coefficients, the reduced frequency can then be modified by varying the bar-to-blade count. Numerical setups with one, two, and four bars are named 1b, 2b, and 4b, respectively. The resulting combinations of flow coefficient and reduced frequency are summarized in Table 3.

Table 3

Parameter space

SimulationNbar/NbladeFredΦGw/CRe2isMa2isα1 (deg)α2 (deg)Reθ1,z=0
1b1u10.212.250.59859530.5941.3266.25423
1b2u10.421.110.49858240.5941.8165.96394
1b3u10.630.730.42864610.5941.3766.50209
2b1u20.422.220.30851720.5945.6566.57318
2b2u20.851.100.25862090.5940.7766.70185
2b3u21.270.710.21842050.5838.6565.74130
4b1u40.852.170.15846740.5839.9666.23192
4b2u41.701.050.12835690.5835.7565.61126
SimulationNbar/NbladeFredΦGw/CRe2isMa2isα1 (deg)α2 (deg)Reθ1,z=0
1b1u10.212.250.59859530.5941.3266.25423
1b2u10.421.110.49858240.5941.8165.96394
1b3u10.630.730.42864610.5941.3766.50209
2b1u20.422.220.30851720.5945.6566.57318
2b2u20.851.100.25862090.5940.7766.70185
2b3u21.270.710.21842050.5838.6565.74130
4b1u40.852.170.15846740.5839.9666.23192
4b2u41.701.050.12835690.5835.7565.61126

The resulting wake-to-wake spacing non-dimensionalized by chord length, Gw/C, evenly covers the value range reported in Ref. [12] for the mid-span section of a T106A profile. Small values of Gw/C are associated with wake fogging, a phenomenon where the incoming wakes merge before entering the blade passage and the incoming disturbances, as a result, loose their discreteness and resemble more closely inflow turbulence with an elevated level [12]. On the other hand, large values of Gw/c are associated with discrete wakes entering the passage leading to profound unsteadiness around the blade. In summary, the chosen parameter space covers engine-relevant conditions for LPT designs which are typically around Fred0.6,ϕ0.75 with respect to the definitions given above spanning the simulations 1b3u, 2b1u, 2b2u, and 2b3u [6]. Less engine-realistic conditions achieved by extreme values of ϕ and Fred, such as 1b1u, 4b1u, and 4b2u, are necessary to determine trends as a function of Gw/C and validate the simulations against the experimental data approximated with 1b1u. In contrast to Ref. [22], the total pressure drop and flow turning associated with the moving bars are not corrected for, which is consistent with previous parametric sweep studies [12]. The resulting change of inlet flow angle behind the bars, α1, reported in Table 3, shows little deviation with increasing slide speed and bar count for almost all simulations, except the extreme case of four bars and the doubled slide speed (4b2u). Here, the absolute flow angle reduces about 6 deg from the nominal case (1b1u) and the isentropic exit Reynolds number Re2is decreases by less than 3% which is expected. However, for this case, due to the very small value of Gw/C, we expect significant wake fogging and therefore less contribution of the unsteady losses. On the other hand, the outlet flow angle α2 and isentropic exit Mach number Ma2is remain nearly constant for all cases. The resulting drift of operating conditions and inflow angle can therefore be assumed to have negligible effects. At the spanwise end-walls, the moving bars distort the prescribed boundary layer profile resulting in a drop of momentum thickness Reynolds number Reθ1,z=0=θρBLEuBLE/μBLE at Sec. 1. With more incoming wakes, i.e., lower values of Gw/C, the end-wall boundary layers get disturbed more drastically resulting in a further drop of Reθ1,z=0 summarized in the last row of Table 3.

3 Mid-Span Analysis

3.1 Blade Profiles.

At mid-span, the time-averaged pressure coefficient defined as
(4)
using time-averaged pressure pw around the blade is shown in Fig. 4. The configurations with one bar and varying slide speed are displayed on the left plot. As expected, the effective incidence angle of the blade notably decreases with increasing slide speed at the stagnation point Cax=0 and further results in reduced loading on the aft portion of the suction side, i.e., for x/Cax<0.675. The pressure side loading remains unaltered. For increasing bar count at constant slide speed, one can see a further reduction of the blade loading in the aft portion of the suction side due to the change of incidence angle. This also explains the profound pressure drop on the pressure side in close vicinity of the stagnation point. However, toward the suction side trailing edge, all profiles collapse, indicating a closed separation bubble. This is confirmed by the positive values of the time-averaged skin friction coefficient at the trailing edge, displayed in Fig. 5 showing the rear section of the suction side over x/Cax=[0.75,1.0]. Here, each simulation is plotted in a separate subplot with cf¯ being the solid line. The time-averaged (or mean) on-surface flow separation is marked by the vertical black-dashed line and its reattachment is marked by the gray-dashed line, delineating the separation bubble for which cf¯<0.0 holds.
Fig. 4
Blade pressure coefficient in mid-span plotted over axial chord length
Fig. 4
Blade pressure coefficient in mid-span plotted over axial chord length
Close modal
Fig. 5
Blade skin friction coefficient at mid-span in the vicinity of the separation bubble on the blade suction side. Time-averaged distribution plotted as solid lines with vertical dashed lines indicating time-averaged flow separation (first vertical line) and flow reattachment (second vertical line). Added contours represent probability density function of the skin friction coefficient around its mean value. The mean separation bubble length non-dimensionalized with chord length over all cases as a function of the wake-to-wake width is displayed in bottom right graph.
Fig. 5
Blade skin friction coefficient at mid-span in the vicinity of the separation bubble on the blade suction side. Time-averaged distribution plotted as solid lines with vertical dashed lines indicating time-averaged flow separation (first vertical line) and flow reattachment (second vertical line). Added contours represent probability density function of the skin friction coefficient around its mean value. The mean separation bubble length non-dimensionalized with chord length over all cases as a function of the wake-to-wake width is displayed in bottom right graph.
Close modal

The probability density function (PDF) of the cf time-series is additionally indicated by the gray-scale map ranging from light-gray (low values of the PDF) to black (high values of the PDF) around the time-averaged mean. Across all simulations, the PDF shows high concentration upstream of the separation onset (i.e., small variation around the mean) while downstream of the separation a much broader spectrum of gray scales is observed. Here, positive and negative areas of cf are found, indicating the intermittent opening and closing of the separation bubble due to the incoming wakes impinging on the blade.

The time-averaged size of the separation bubble is computed as lbubble=(sreattachmentsseparation) where s denotes the blade path length. The bottom right graph in Fig. 5 plots lbubble/C over the wake-to-wake width Gw/C across the eight LES cases. The circle color represents the simulations titled in their respective cf graph, e.g., orange denotes the 1b2u case. One can observe that the bubble length as a function of Gw/C collapses onto a trend line indicated by the dashed black line with a local minimum around Gw/C0.4. This trend, first reported for spanwise periodic setups [9,12], therefore also appears to hold true in the presence of an additional axial pressure gradient introduced by the divergent case path and for a section (T161) with a different loading than the T106A.

3.2 Wake Profiles.

To evaluate the mid-span blade performance, the total pressure loss coefficient throughout is defined as
(5)
and taken at an axial location of 40% axial chord length downstream of the blade trailing edge coinciding with the experimental setup and Sec. 2. The time-averaged pitchwise distribution is given in Fig. 6 for the eight LES cases. Increasing values of y/lpitch, beyond the respective maximum of ζ¯, refer to the suction side related losses, and decreasing values to the pressure side related losses. The left hand plot in the top column shows the one bar configurations with increasing slide speed, resulting in different Φ and Fred combinations. Generally speaking, the wake width and wake maxima are decreasing with increasing slide speed. This can be explained by the smaller separation bubble on the blade suction side (see Fig. 5) for increasing bar speeds. Interestingly, a redistribution of the losses in the off-peak wake region can further be observed in the present simulations when increasing the bar count, highlighted in subplot (b). The suction side losses predominately shrink with increasing slide speed around y/lpitch=0, while pressure side losses grow marginally for y/lpitch0.3. Eventually, in the case of four bars, wake fogging occurs. The blade boundary layers do not experience the discrete, bar-induced, unsteadiness necessary to suppress the separation bubble. The bar-wake mixing results in elevated levels of background turbulence in the passage and, in a time-averaged sense, widening of the separated region compared to the two bar cases. As a result, the wake width and wake magnitude increase (see subplot (c)). For increasing bar count, at constant slide speed, the wake width narrows, as the suction side separation bubble decreases in size (see Fig. 5, bottom column), until wake fogging occurs. Despite the minimal separation bubble size for 1b3u, the wake loss can be minimized by finding the optimal combination of bar count and bar speed which is 2b3u in the chosen parameter space (see Fig. 6(f)) as the wake loss toward its suction side onset at y/lpitch=0 is minimized. This is caused by the presence of an additional axial pressure gradient due to the spanwise diverging end-walls which is contrasting the observations made for the T106A case in a spanwise periodic configuration [9,12]. Hence, increasing the unsteadiness by decreasing Gw/C has a positive effect on loss generation until wake fogging occurs.
Fig. 6
Mid-span total pressure loss coefficient at 40% axial chord length downstream of the blade trailing edge plotted over pitchwise coordinate. Fixed bar count and variable bar speed are shown in the top column, and variable bar count and fixed bar speed are shown in the bottom column.
Fig. 6
Mid-span total pressure loss coefficient at 40% axial chord length downstream of the blade trailing edge plotted over pitchwise coordinate. Fixed bar count and variable bar speed are shown in the top column, and variable bar count and fixed bar speed are shown in the bottom column.
Close modal
In order to isolate the unsteady loss generation associated with the incoming wakes, the loss contribution of the wake distortion can be defined as
(6)
where ζM denotes total pressure coefficient based on mixed-out quantities and ζDenton~ are the profile losses computed according to Denton [34] with the modifications introduced in Refs. [9,12]. Boxplots of the phase-averaged quantities, indicated by the overhead tilde, are given in Fig. 7. The mean values of all phases are denoted by the circles, on which the black-dashed trend line is based. The outlined black box spans the middle 50% around the mean, and phase-maxima and minima are the vertical bounds. The boxplots indicate large vertical ranges, associated with high levels of unsteadiness, for high values of Gw/C. The integrated total pressure loss can more than double depending on the phase. Low values of Gw/C, where wakes merge and wake fogging occurs, effectively result in overall increased levels of turbulence, but less discrete wake-related intermittency, and leads to smaller variations around the mean value. In summary, the unsteady loss trend with respect to the wake-to-wake width Gw/C is identical to that reported in previous work [9,12] even though the integrated total pressure loss coefficient trend is of opposite direction. This can be explained by investigating the extreme values of Δ~. First, the local maximum is for 1b3u as the suction side separation bubble is perturbed heavily leading to an intermittently fully attached boundary caused by a high flow coefficient and minimize the bubble size in a time-averaged sense (see Fig. 5). As a result, Denton’s loss correlation ζDenton~ is minimized for this Gw/C value. But, the Gw/C spacing between incoming wakes is high enough to allow full recovery of an open separation bubble, leading to profound unsteadiness and deviations from Denton’s correlations as highlighted in the boxplots. The local minimum is observed for 2b3u. In this case, ζDenton~ is close to the 1b3u case. However, the additional axial pressure gradient has a positive impact on the losses itself, as discussed above, and, on top of that, decrease the unsteadiness with respect to incoming wakes Gw/C such that Δ~ is minimized.
Fig. 7
Mixed-out losses, Denton losses, and wake distortion losses at mid-span
Fig. 7
Mixed-out losses, Denton losses, and wake distortion losses at mid-span
Close modal

4 End-Wall Analysis

4.1 Phase-Averaged Loss Mechanisms.

A detailed discussion of the evolving secondary flow phenomena have been provided in previous work, for steady conditions in Ref. [23] and for a single set of Φ and Fred in Ref. [22], and will therefore not be repeated here. This paper focuses on the evolution of those structures with varying bar wake-to-wake width Gw/C. Looking at phase-averaged results for a fixed sliding speed and various bar configurations, Fig. 8 reveals multiple phase-dependent mechanisms. For the three bar configurations shown, four phases have been chosen from the interval ϕ[0;2π] spanning one single bar wake-to-wake phase. At mid-span, the phase-averaged turbulent kinetic energy k~=1/2(uiui~)=1/2((uiui~)(uiui~)~) highlights the bar wakes convecting through the passage. All surfaces are colored by the skin friction coefficient and streaklines are represented by a line integral convolution of the wall shear stress vector. White iso-lines of cf~=0 further indicate separated flow regions. Toward the spanwise end-walls, three axial constant slices show the phase-averaged streamwise vorticity ωsw~=uω/u~ in Cartesian coordinates ω.

Fig. 8
Phase-averaged, three-dimensional flow structures for varying bar configurations. Mid-span slice of turbulent kinetic energy k~ and end-wall slices of streamwise vorticity ωsw~. Surfaces are colored by friction coefficient cf~ with iso-contours of cf~=0 and streaklines are indicated by line integral convolution.
Fig. 8
Phase-averaged, three-dimensional flow structures for varying bar configurations. Mid-span slice of turbulent kinetic energy k~ and end-wall slices of streamwise vorticity ωsw~. Surfaces are colored by friction coefficient cf~ with iso-contours of cf~=0 and streaklines are indicated by line integral convolution.
Close modal

The passage vortex (PV), originating from the pressure-sided leg of the horse-shoe vortex, is the concentrated region of blue color in the -marked graph. The trailing shed vortex (TSV) is colored in red emerging from the trailing edge and rotating in the opposite direction to the passage vortex. Looking at phase ϕ=π/2 of the one bar configuration, the bar wake is close to the trailing edge. As a result, the separation bubble decreases in size visible by the smaller enclosed region of the white iso-lines. Toward the spanwise end-wall, the wake mixes with the emerging passage vortex.

The next phase, ϕ=π, shows locations of fully attached flow which is likely a result of spanwise-varying turbulent spots [22]. Compared to , the separation line of the passage vortex moves downstream, indicated by . Its overall structure remains intact despite the high level of fluctuations from the unsteady wakes.

At ϕ=3/2π, the wake has convected through the passage. In its absence, the separation bubble widens and the passage vortex recovers close to the trailing edge. Comparing at the same phase to the case with two bars shows a smaller separation region for the latter as the overall level of unsteadiness increases in the passage region.

The separation bubble grows and shrinks depending on the phase with locally fully attached flow at ϕ=3/2π while the passage vortex separation line remains nearly unaffected. However, due to increased levels of turbulent kinetic energy, the separation line is mostly attached for the two bar case. This results in a weakening of the passage vortex and subsequently the trailing shed vortex. Comparing the axial constant slices at the trailing edge of the 1b2u and 2b2u cases highlights lower values of streamwise rotation across all phases for the two bar configuration confirming the weakening of secondary flow structures for cases with more bars. Unsurprisingly, the four bar configuration, associated with wake fogging, shows little to no variation over the different phases. The high background turbulence level resulting from the interaction and mixing of the incoming bar wakes, clearly visible in mid-span, suppresses the detachment of the passage vortex from the blade surface indicated by . Therefore, and due to strong mixing in the passage, the passage vortex and trailing shed vortex are less profound downstream of the wake.

4.2 Time-Averaged Wake Profiles.

The evolving secondary vortex patterns and their relation toward loss generation in a time-averaged manner can be visualized by contour plots of the total pressure loss coefficient (see Eq. (5)) across all eight LES cases in Fig. 9. In all subplots, the spanwise coordinate is normalized over the local spanwise extend on the x-axis and the normalized pitchwise coordinate is plotted on the y-axis. In the top left plot, the 1b1u further highlights the PV, TSV, and the end-wall separated region (ES) for clarity. Increasing the bar speed for cases with one bar is shown in the left column. Looking from 1b1u to 1b3u reveals a gradual decrease in the size of the trailing shed vortex. Simultaneously, the passage vortex decreases in strength and appears to merge with the trailing shed vortex. This is caused by increased perturbations of the secondary vortex system, which forms as a consequence of the spanwise end-walls, due to the presence of more wakes, i.e., with decreasing Gw/C. A further decrease in Gw/C, by increasing the bar count, leads to an amplified disturbance of the secondary vortex system and a further reduction in size of the trailing shed vortex. Interestingly, in case of wake fogging, which occurs for the 4b2u case, the passage vortex starts to grow in strength again which can be seen when comparing to the 1b3u or 2b3u cases. This can be explained by the fact that the passage vortex and trailing shed vortex start to spatially separate further due to high background turbulence levels rather than discrete wakes. On the other hand, the end-wall region showing non-zero ζ, labeled ES, is continuously decreasing in size when decreasing Gw/C, i.e., comparing 1b1u and 4b2u. This is a result of a decreased inflow boundary layer momentum thickness at the spanwise end-walls upstream of the blade (see Table 3) due to the mixing of the boundary layer with the moving bars. Simulations in the absence of incoming wakes have shown that decreased levels of incoming end-wall boundary layer momentum thickness directly correlate to the end-wall associated losses in the wake [23]. This holds true for the unsteady cases as more incoming wakes continuously disturb the freestream boundary layer resulting in a drop of Reθ1,z=0 and, hence, less total pressure drop in the wake.

Fig. 9
Time-averaged wake loss contours at 40% axial chord lengths downstream of the blade
Fig. 9
Time-averaged wake loss contours at 40% axial chord lengths downstream of the blade
Close modal

4.3 Integral Loss Parameters.

In order to quantify the three-dimensional integral losses across all simulations, the pitchwise and spanwise mass-averaged, total pressure loss coefficient can be computed as
(7)
and can be decomposed into profile losses ζp, end-wall boundary layer distortion losses ζbl, and secondary losses ζsec [11]. Once again, the phase-averaged quantities are chosen to highlight the unsteadiness around their mean values across all cases. The secondary losses ζsec are representative of the passage vortex and trailing shed vortex. In contrast to the profile losses ζp (top graph in Fig. 10), the black-dashed ζsec trend line does not show a local minimum when plotted against Gw/C, but instead steadily decreases when Gw/C tends to zero (second graph in Fig. 10). Extreme cases of wake fogging, such as the case 4b2u, are therefore able to reduce the secondary loss generation even given a re-growth of the passage vortex for this case. Hence, the trailing shed vortex is the dominant source of losses of the forming vortex system. For the geometry considered in this study, the losses associated with the variation of incoming boundary layers upstream of the blade are isolated to the first 10% relative span height [23], and are gradually decreasing as a function of Gw/C, as depicted in the third graph of Fig. 10 and as outlined above. The bottom graph in Fig. 10 plots ζint~, the sum of the individual loss terms. Its black-dashed trend line follows the trend line of the profile losses over Gw/C, highlighting the importance of blade design for LPTs. However, for the cases where total pressure loss is minimized, secondary losses can contribute up to 50% of the total loss generation, as seen for the 2b1u case.
Fig. 10
Phase-averaged, three-dimensional loss breakdown at 40% axial chord length downstream of the blade
Fig. 10
Phase-averaged, three-dimensional loss breakdown at 40% axial chord length downstream of the blade
Close modal

5 Conclusion

The present study reports on highly resolved large eddy simulations of a linear low-pressure turbine cascade in a diverging gas path. Rotor-like wakes are generated by upstream moving bars and periodically impinge on the blade leading edge. A grid convergence study and bench-marking against available experimental data validate the present numerical setup. In total, eight LES comprise a set of systematically varying flow coefficient and reduced frequency such that the wake-to-wake width Gw/C is evenly distributed over a wide range. While the total pressure loss and flow turning caused by the bar wakes are not corrected, their impact on operating conditions such as isentropic exit Reynolds number and isentropic exit Mach number remains small.

First, the mid-span section is analyzed. The separation bubble on the blade suction side is heavily influenced by the intermittent wakes impinging on the blade. Comparing phase-averaged results, high values of Gw/C result in strong, intermittent lengthening and shortening of the separated region due to the high-amplitude unsteadiness of the wakes periodically disturbing the blade boundary layers. On the other hand, wake fogging and merging is associated with low values of Gw/C. This leads to elevated levels of background turbulence in the passage, rather than discrete wakes, experienced by the blade boundary layers. The time-averaged bubble size therefore shows a minimum between the extreme cases of low and high wake-to-wake width. The total pressure loss coefficient in the wake downstream of the blade follows a similar trend as the separation bubble size with respect to Gw/C. However, its local minimum is at a decreased Gw/C value compared to the minimum of the bubble size. The trend line of the wake loss over Gw/C is opposite to those in previous spanwise periodic studies on the T106A profile which is caused by two distinct phenomena. First, the time-averaged suction side separation bubble decreases in size when Gw/C decreases, with a local minimum between the extreme values of Gw/C. And, second, the additional axial pressure gradient in mid-span impacts the incoming wakes positively by reducing the total pressure loss at the wake onset of its suction side related branch.

The profile losses calculated following Denton’s loss analysis are subtracted from the LES losses in order to isolate the wake distortion losses. For the current case, an identical trend of the wake distortion losses over wake-to-wake width is observed compared to the previous spanwise periodic studies on the T106A profile. The maximum wake distortion losses are associated with high unsteadiness of the blade boundary layer state along the suction side. In this scenario, the incoming discrete wakes are able to intermittently fully suppress the separated region. Yet, the wake-to-wake width is sufficiently high for the separation bubble to recover before the next wake affects it. On the other hand, the minimum of the wake distortion losses occurs at low unsteadiness around the separation bubble and the minimum of the total pressure loss coefficient. At low values of Gw/C, the wake distortion losses are also decreased as the incoming wakes merge. In summary, the wake distortion losses are characterized by their interaction with the blade profile and the axial pressure gradient rather than the unsteady nature of the incoming wakes themselves.

The secondary flow system toward the spanwise end-walls is found to remain mostly intact for all cases considered. The incoming secondary structures of the bars are not able to fully suppress the formation of the horse-shoe vortex which feeds the passage vortex and trailing shed vortex. However, the passage vortex and trailing shed vortex significantly reduce in strength with decreasing wake-to-wake width and are therefore more sensitive toward elevated background turbulence levels in the passage rather than discrete incoming wakes. The end-wall losses caused by the end-wall boundary layer deficit are decreasing in a similar manner due to the increased mixing of the freestream boundary layer with the moving bars upstream of the blade and a subsequent reduction in end-wall boundary layer momentum thickness. The trend lines of the secondary vortex system toward loss generation are therefore steadily decreasing with decreasing values of Gw/C. Comparing the total integrated loss coefficient across all simulations reveals a similar trend to those of the profile losses, with the profile losses being the predominant source of loss generation and highlighting the importance of mid-span blade design. However, for the 2b3u configuration showing the smallest total pressure loss coefficient, secondary losses amount to 50% of the total losses.

Acknowledgment

The bulk of the computing time has been provided by the Pawsey Supercomputing Center in Perth, Australia, awarded by the National Computational Merit Allocation Scheme. Furthermore, this research was undertaken using the LIEF HPC-GPGPU Facility hosted at the University of Melbourne. This Facility was established with the assistance of LIEF Grant LE170100200. Experimental data were provided under a non-disclosure agreement between MTU Aero Engines, A.G. and The University of Melbourne.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

Nomenclature

l =

length

p =

pressure

u =

velocity vector in Cartesian coordinates

C =

chord

cf =

skin friction coefficient

cp =

pressure coefficient

x,y,z =

Cartesian coordinates

Greek Symbols

α =

absolute flow angle

Δ =

wake distortion losses

Δs+ =

viscous point spacing in the blade tangential direction

Δn+ =

viscous point spacing in the blade normal direction

Δz+ =

viscous point spacing in the blade spanwise direction

ζ =

total pressure loss coefficient

ν =

kinematic viscosity

ρ =

density

τw =

wall shear stress

uτ =

wall friction velocity

Φ =

flow coefficient

ϕ =

bar passing phase

Superscripts and Subscripts

1 =

quantity at section 1

2 =

quantity at section 2

m =

mixed-out value

t =

total

ax =

axial

is =

isentropic

BLE =

boundary layer edge

LE =

leading edge

TE =

trailing edge

¯ =

time-averaged value

~ =

phase-averaged value

=

freestream value

Dimensionless Groups

Fred =

reduced frequency

Gw =

wake-to-wake width

Ma =

Mach number

Re =

Reynolds number

Zw =

Zweifel number

References

1.
Sieverding
,
C. H.
,
1985
, “
Recent Progress in the Understanding of Basic Aspects of Secondary Flows in Turbine Blade Passages
,”
ASME J. Eng. Gas Turbines Power
,
107
(
2
), pp.
248
257
.
2.
Goldstein
,
R. J.
, and
Spores
,
R. A.
,
1988
, “
Turbulent Transport on the Endwall in the Region Between Adjacent Turbine Blades
,”
ASME J. Heat Transfer-Trans. ASME
,
110
(
4a
), pp.
862
869
.
3.
Langston
,
L.
,
2001
, “
Secondary Flows in Axial Turbines—A Review
,”
Ann. N. Y. Acad. Sci.
,
934
(
1
), pp.
11
26
.
4.
Stadtmueller
,
Peter
,
Fottner
,
Leonhard
, and
Fiala
,
Andreas
,
2000
, “
Experimental and Numerical Investigation of Wake-Induced Transition on a Highly Loaded LP Turbine at Low Reynolds Numbers
,”
Turbo Expo: Power for Land, Sea, and Air
,
78569
.
5.
Sandberg
,
R. D.
, and
Michelassi
,
V.
,
2022
, “
Fluid Dynamics of Axial Turbomachinery: Blade- and Stage-Level Simulations and Models
,”
Annu. Rev. Fluid Mech.
,
54
, pp.
255
285
.
6.
Martinstetter
,
M.
,
2009
, “
Experimentelle Untersuchungen zur Aerodynamik hoch belasteter Niederdruckturbinen-Beschaufelungen
,” Dissertation, University of Armed Forces Munich.
7.
Homeier
,
L.
,
2006
, “
Experimentelle Untersuchungen an hochbelasteten Hochdruckturbinengittern mit Filmkuhlausblasung in Gebieten mit positiven Druckgradienten
,” Dissertation, University of Armed Forces Munich.
8.
Sandberg
,
R. D.
,
Michelassi
,
V.
,
Pichler
,
R.
,
Chen
,
L.
, and
Johnstone
,
R.
,
2015
, “
Compressible Direct Numerical Simulation of Low-Pressure Turbines—Part I: Methodology
,”
ASME J. Turbomach.
,
137
(
5
), p.
051011
.
9.
Michelassi
,
V.
,
Chen
,
L.-W.
,
Pichler
,
R.
, and
Sandberg
,
R.
,
2014
, “
Compressible Direct Numerical Simulation of Low-Pressure Turbines—Part II: Effect of Inflow Disturbances
,”
ASME J. Turbomach.
,
137
(
5
), p.
071005
.
10.
Pichler
,
R.
,
Zhao
,
Y.
,
Sandberg
,
R. D.
,
Michelassi
,
V.
,
Pacciani
,
R.
,
Marconcini
,
M.
, and
Arnone
,
A.
,
2019
, “
Large-Eddy Simulation and RANS Analysis of the End-Wall Flow in a Linear Low-Pressure Turbine Cascade, Part I: Flow and Secondary Vorticity Fields Under Varying Inlet Condition
,”
ASME J. Turbomach.
,
141
(
12
), p.
121005
.
11.
Marconcini
,
M.
,
Pacciani
,
R.
,
Arnone
,
A.
,
Michelassi
,
V.
,
Pichler
,
R.
,
Zhao
,
Y.
, and
Sandberg
,
R.
,
2019
, “
Large Eddy Simulation and RANS Analysis of the End-Wall Flow in a Linear Low-Pressure-Turbine Cascade—Part II: Loss Generation
,”
ASME J. Turbomach.
,
141
(
5
), p.
051004
.
12.
Michelassi
,
V.
,
Chen
,
L.
,
Pichler
,
R.
,
Sandberg
,
R.
, and
Bhaskaran
,
R.
,
2016
, “
High-Fidelity Simulations of Low-Pressure Turbines: Effect of Flow Coefficient and Reduced Frequency on Losses
,”
ASME J. Turbomach.
,
138
(
11
), p.
111006
.
13.
Cui
,
J.
,
Rao
,
V.
, and
Tucker
,
P.
,
2016
, “
Numerical Investigation of Secondary Flows in a High-Lift Low Pressure Turbine
,”
Int. J. Heat Fluid Flow
,
63
, pp.
149
157
.
14.
Koschichow
,
D.
,
Fröhlich
,
J.
,
Ciorciari
,
R.
, and
Niehuis
,
R.
,
2015
, “
Analysis of the Influence of Periodic Passing Wakes on the Secondary Flow Near the Endwall of a Linear LPT Cascade Using DNS and U-RANS
,”
11th European Conference on Turbomachinery Fluid Dynamics and Thermodynamics
,
Madrid, Spain
, Paper No. WETC2015-151.
15.
Robison
,
Z.
, and
Gross
,
A.
,
2022
, “
Comparative Numerical Investigation of Wake Effect on Low-Pressure Turbine Endwall Flow
,”
Aerosp. Sci. Technol.
,
131
, p.
107970
.
16.
Lopes
,
G.
,
Simonassi
,
L.
, and
Lavagnoli
,
S.
,
2023
, “
Impact of Unsteady Wakes on the Secondary Flows of a High-Speed Low-Pressure Turbine Cascade
,”
Int. J. Turbomach. Propul. Power
,
8
, p.
36
. DOI: 10.3390/ijtpp8040036
17.
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: J. Power Energy
,
211
(
4
), pp.
309
320
.
18.
Rosenzweig
,
M.
,
Giaccherini
,
S.
,
Pinelli
,
L.
,
Kozul
,
M.
,
Sandberg
,
R. D.
,
Marconcini
,
M.
, and
Pacciani
,
R.
,
2023
, “
Best-practice Guidelines for High-Fidelity Simulations Based on Detailed Analysis of a Highly-Loaded Low-Pressure Turbine Cascade
,”
Turbo Expo: Power for Land, Sea, and Air
, Vol. 138, p.
V13BT30A021
.
19.
Schwarzbach
,
F.
,
Müller-Schindewolffs
,
S.
,
Bode
,
C.
, and
Herbst
,
F.
,
2018
, “
The Effect of Turbulent Scales on Low-Pressure Turbine Aerodynamics: Part B - Scale Resolving Simulations
,”
Proceedings of the ASME Turbo Expo 2018: Turbomachinery Technical Conference and Exposition
, p. V02BT41A005.
20.
Afshar
,
N.
,
Kozulovic
,
D.
,
Henninger
,
S.
,
Deutsch
,
J.
, and
Bechlars
,
P.
,
2023
, “
Turbulence Anisotropy Analysis at the Middle Section of a Highly Loaded 3D Linear Turbine Cascade Using Large Eddy Simulation
,”
J. Global Power Propul. Soc.
,
7
, pp.
71
84
. DOI: 10.33737/jgpps/159784
21.
Iyer
,
A.
,
Abe
,
Y.
,
Vermeire
,
B. C.
,
Bechlars
,
P.
,
Baier
,
R. D.
,
Jameson
,
A.
,
Witherden
,
R.
, and
Vincent
,
P. E.
,
2021
, “
High-Order Accurate Direct Numerical Simulation of Flow Over a MTU-T161 Low Pressure Turbine Blade
,”
Comput. Fluids
,
226
(
10
), p.
104989
.
22.
Morsbach
,
C.
,
Bergmann
,
M.
,
Tosun
,
A.
,
Klose
,
B. F.
,
Bechlars
,
P.
, and
Kügeler
,
E.
,
2024
, “
A Numerical Test Rig for Turbomachinery Flows Based on Large Eddy Simulations With a High-Order Discontinuous Galerkin Scheme—Part III: Secondary Flow Effects
,”
ASME J. Turbomach.
,
146
(
2
), p.
021007
.
23.
Rosenzweig
,
M.
,
Kozul
,
M.
,
Sandberg
,
R.
,
Marconcini
,
M.
, and
Pacciani
,
R.
,
2023
, “
High-Fidelity Simulations and RANS Analysis of a Low-Pressure Turbine Cascade in a Spanwise-Diverging Gas Path - End-Wall Flow Analysis and Loss Generation Mechanisms
,”
Turbo Expo: Power for Land, Sea, and Air
, p.
V13BT30A003
.
24.
Nardini
,
M.
,
Kozul
,
M.
,
Jelly
,
T.
, and
Sandberg
,
R.
,
2024
, “
Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges
,”
ASME J. Turbomach.
,
146
(
3
), p.
031008
.
25.
Nardini
,
M.
,
Jelly
,
T.
,
Kozul
,
M.
,
Sandberg
,
R.
,
Vitt
,
P.
, and
Sluyter
,
G.
,
2024
, “
Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-scale Surface Roughness—Part II: The Effect of Roughness on the Performance of a High-Pressure Turbine Blade
,”
ASME J. Turbomach.
,
146
(
3
), p.
031009
.
26.
Jelly
,
T.
,
Nardini
,
M.
,
Rosenzweig
,
M.
,
Leggett
,
J.
,
Marusic
,
I.
, and
Sandberg
,
R.
,
2023
, “
High-Fidelity Computational Study of Roughness Effects on High Pressure Turbine Performance and Heat Transfer
,”
Int. J. Heat Fluid Flow
,
101
, p.
109134
.
27.
Nicoud
,
F.
, and
Ducros
,
F.
,
1999
, “
Subgrid-Scale Stress Modelling Based on the Square of the Velocity Gradient Tensor
,”
Flow Turbulence Combust.
,
62
(
3
), pp.
183
200
.
28.
Sandberg
,
R. D.
, and
Sandham
,
N. D.
,
2006
, “
Nonreflecting Zonal Characteristic Boundary Condition for Direct Numerical Simulation of Aerodynamic Sound
,”
AIAA J.
,
44
(
2
), pp.
402
405
.
29.
Deuse
,
M.
, and
Sandberg
,
R. D.
,
2020
, “
Implementation of a Stable High-Order Overset Grid Method for High-Fidelity Simulations
,”
Comput. Fluids
,
211
, p.
104449
.
30.
Klein
,
M.
,
Sadiki
,
A.
, and
Janicka
,
J.
,
2003
, “
A Digital Filter Based Generation of Inflow Data for Spatially Developing Direct Numerical or Large Eddy Simulations
,”
J. Comput. Phys.
,
186
(
2
), pp.
652
665
.
31.
Rosenzweig
,
M.
,
Giannini
,
S.
,
Jelly
,
T. O.
, and
Sandberg
,
R. D.
,
2022
, “
High-Fidelity Simulation and Flow Analysis of a Low-Pressure Turbine Cascade in a Spanwise-Diverging Gas Path—Calibration of Inlet Boundary Conditions
,”
Australasian Fluid Mechanics Conference
,
Sydney, Australia
, Vol. 23.
32.
MTU Aero Engines, A.G.
,
2023
, personal communication.
33.
Entlesberger
,
R. G.
,
Martinstetter
,
M.
, and
Staudacher
,
W.
,
2005
, “Untersuchungen am Turbinengitter T161 zur Bestimmung der Profildruckverteilung und der Gittercharakteristik,” Vol. LRT-WE12-05/02, Institutsbericht der Universität der Bundeswehr München, Neubiberg, Germany.
34.
Denton
,
J. D.
,
1993
, “
Loss Mechanisms in Turbomachines
,”
Turbo Expo: Power for Land, Sea, and Air
, p.
V002T14A001
.