Abstract

Compressible direct numerical simulations are conducted to investigate how surface roughness affects the aerothermal performance of a high-pressure turbine vane operating at an exit Reynolds number of 0.59 ×106 and exit Mach number of 0.92. The roughness under investigation here was synthesized with non-Gaussian statistical properties and an amplitude that varies over its chord length, representative of what truly occurs on an in-service vane. Particular attention is directed toward how systematically changing the axial extent of leading edge roughness affects convective heat transfer (Nusselt and Stanton numbers) and aerodynamic drag (skin friction coefficient) on the pressure and suction surfaces. The results of this investigation demonstrate that moving the larger amplitude roughness further along the suction surface can alter the blade boundary layer state. In fact, toward the trailing edge of one of the rough vanes investigated here, the local skin friction coefficient increases by a factor of 22 compared to smooth-vane levels, whereas the local Nusselt number increases by a factor 6. The disproportionate rise of drag compared to heat transfer is explored in further detail by quantifying the Reynolds’ analogy and by calculating the fractional contributions of pressure drag and viscous drag to the total drag force. The effect of varying the inlet turbulence intensity and integral length scale for a fixed roughness topography is also investigated, and the Reynolds number scaling of heat transfer and drag is examined in the context of the Chilton–Colburn analogy.

1 Introduction

Surface roughness arising through in-service wear and/or the manufacturing process plays a critical role in determining high-pressure turbine (HPT) aerothermal performance [1,2]. While it is well known that roughness tends to augment heat transfer at the expense of an even greater augmentation in drag [35], detailed accounts of the fundamental flow parameters that govern convective heat transfer (wall heat flux, Nusselt, and/or Stanton numbers) and boundary layer losses (wall shear stress, skin friction coefficient, momentum thickness Reynolds number) on rough vanes are still lacking, especially at engine-relevant conditions. Obtaining these data experimentally can be extremely challenging, and although there are some notable exceptions to this trend [6,7], the most detailed studies to date have focused on smooth vanes, e.g., see work by Thole et al. [8].

At the same time, high-fidelity simulations of flow past HPT vanes with realistic levels of roughness have only very recently become possible [9]. This is because of the tremendous computational resources required not only to resolve the full spectrum of hydrodynamic and thermodynamic scales in the flow field but also resolve the topographical scales of the surface roughness. The very first compressible direct numerical simulation (DNS) of flow past an HPT vane with realistic levels of roughness was conducted in Refs. [10,11], which built upon a precursory series of compressible large-eddy simulations performed by Ref. [12]. These past studies quantified the aerothermal performance of the (uncooled) VKI LS89 nozzle guide vane [13] covered with numerically generated multiscale surface roughness [14] at an exit Reynolds number of 590,000 and an exit Mach number of 0.92, relative to a smooth vane at matched conditions. While Refs. [1012] highlighted the profound impact that surface roughness can have upon HPT aerothermal performance, these studies mainly focused on near-Gaussian roughness (i.e., roughness with a symmetric height distribution) with a uniform roughness height on the suction and pressure surfaces.

However, as shown in measurements by Refs. [15,16], roughness is not homogeneously distributed over the surface of in-service vanes. In the context of first-stage HPT vanes, the highest roughness levels tend to occur in the leading edge region and along the pressure surface due to direct exposure to combustor exhaust gases, whereas the lowest roughness levels occur along the sheltered portion of the suction surface. Furthermore, it is well known that most turbine blade surface roughness is non-Gaussian [17], i.e., the height distribution is nonsymmetric and has nonzero skewness and kurtosis not equal to 3. As a result, systematic studies of how localized, non-Gaussian surface roughness—as it truly occurs on in-service vanes—affects HPT aerothermal performance are of significant practical interest.

To this end, a first-of-its-kind series of compressible DNSs with up to 14.6×109 grid points per simulation have been conducted, with the specific goal of performing a direct computational assessment of convective heat transfer, skin friction, and the Reynolds” analogy on realistically rough HPT vanes at engine-relevant conditions. Particular attention is directed toward understanding how systematically varying the axial extent of leading edge roughness affects the instantaneous (and averaged) quantities that govern HPT aerothermal performance, i.e., the skin friction coefficient and the Nusselt and Stanton numbers, and how these quantities are augmented (or reduced) compared to smooth-vane levels. The impact of varying the inflow turbulence intensity and integral length scale for a fixed roughness topography is also investigated, and, where appropriate, the Reynolds number scaling of the HPT data is compared against existing correlations for turbulent heat transfer and drag.

2 Computational Aspects

2.1 Blade Generation Procedure and Roughness Parameters.

The blade generation procedure adopted here is an extension of the method devised by Refs. [1012]. However, in those past studies, the main focus was on nonlocalized (i.e., uniformly distributed) roughness with a near-Gaussian distribution, as opposed to the localized non-Gaussian roughness under investigation here.

HPT vanes covered with localized non-Gaussian roughness were generated by taking the following steps. First, a surface generation algorithm [18] was used to create a doubly periodic planar roughness distribution, h(x,z), with zero mean height, specified root-mean-squared (RMS) roughness height, RqA1Ah2(x,z)dA, and non-Gaussian statistical properties, i.e., nonzero skewness SA1Rq3Ah3(x,z)dA0 and kurtosis not equal to 3, KA1Rq4Ah4(x,z)dA3, where A is the planform area of the heightmap. Here, x and z, respectively, denote the streamwise and spanwise directions of the planar roughness distribution. The statistical properties of h(x,z) were inferred from a set of roughness measurements provided by GE Aerospace. The key amplitude parameters of the planar roughness distribution are (Rq/c,S,K)=(4.2×103,0.42,3.75), where c is the axial blade chord. For a typical axial blade chord of, say, c=25mm, the physical RMS roughness height is therefore Rq100μm, which is typical of the most severely degraded portions of vanes, particularly in the leading edge region [17].

Second, the roughness distribution was wrapped around the reference smooth HPT coordinates by transforming the planar coordinate system of the roughness heightmap, (x,z), to the curvilinear coordinate system of the smooth blade, (s,z), where s is the curvilinear abscissa along the blade tangential direction. Mathematically, the coordinates of a generic point on the surface of the rough HPT vane, xr, can be expressed as follows:
(1)
where xs is the corresponding point on the smooth blade surface with an outward-pointing normal unit vector, ns. The term W(x) in Eq. (1) is a (spanwise homogeneous) weighting function used to control the roughness amplitude around the vane, whose shape was inferred from data provided by GE Aerospace. Here, the goal was to mimic the wear pattern of an in-service HPT vane, i.e., to prescribe larger amplitude roughness in the leading edge region and along the pressure surface, and smaller amplitude roughness along the sheltered portion of the suction surface.

The weighting function, W, is plotted against axial position normalized by the axial chord, x/c, in Fig. 1, where x/c<0 corresponds to the vane suction surface and x/c>0 corresponds to the vane pressure surface. Looking from left-to-right across Fig. 1, it can be observed that all the control points remain fixed with respect to x/c, with the exception of control point IV (white circles), which takes on one of four possible values along the suction surface, i.e., x/c={0.28,0.40,0.50,0.60}. In terms of amplitude, the weighting function increases from control point I at the leading edge (LE) (x/c=0.0), to a maximum at control point III (x/c0.15), before decreasing to a minimum at control point IV, beyond which the weighting function remains constant until control point V, which is located at the vane trailing edge (x/c=1.0). In contrast, the weighting function has a constant amplitude between control points I and VI, which covers most of the pressure surface. However, the weighting function alters the local roughness height, skewness and kurtosis remain unaffected, since standardized moments of a probability density function are scale invariant by definition (see Table 1).

Fig. 1
Roughness weighting function, W, versus x/c
Fig. 1
Roughness weighting function, W, versus x/c
Close modal
Table 1

Roughness amplitude parameters

Pointx/cRq/cSK
I0.003.7 ×1040.423.75
II0.053.9 ×1040.423.75
III0.144.2 ×1040.423.75
IV{0.28,0.40,0.50,0.60}1.1 ×1040.423.75
V1.01.1 ×1040.423.75
VI0.883.9 ×1040.423.75
Pointx/cRq/cSK
I0.003.7 ×1040.423.75
II0.053.9 ×1040.423.75
III0.144.2 ×1040.423.75
IV{0.28,0.40,0.50,0.60}1.1 ×1040.423.75
V1.01.1 ×1040.423.75
VI0.883.9 ×1040.423.75

A zoomed-in view of the non-Gaussian surface roughness in the leading edge of each HPT vane under investigation here is provided in Fig. 2. While the relative changes in the roughness topography may at first appear too small to significantly affect HPT aerothermal performance, the results of this investigation suggest otherwise.

Fig. 2
Zoomed-in view of HPT leading edge colored by local roughness height normalized by axial blade chord, h/c: (a) case X28_0808, (b) case X40_0808, (c) case X50_0808, and (d) case X60_0808
Fig. 2
Zoomed-in view of HPT leading edge colored by local roughness height normalized by axial blade chord, h/c: (a) case X28_0808, (b) case X40_0808, (c) case X50_0808, and (d) case X60_0808
Close modal

2.2 Compressible Direct Numerical Simulations.

Transonic flow past rough HPT vanes was simulated using the HiPSTAR flow solver [19], which solves the compressible Navier–Stokes equations nondimensionalized by a reference length (c= axial blade chord), velocity (Uin= inflow velocity), density (ρin= inlet density), and temperature (Tin= inflow temperature). The governing equations for mass, momentum, and total energy solved by HiPSTAR can be written in the tensor form as follows:

(2)
(3)
(4)
where repeated i and j indices imply summation and δij is the Kronecker delta. In Eqs. (2)(4), ρ(x,t), ui(x,t), p(x,t), and E(x,t) are the density, velocity vector, pressure, and total energy, respectively, where x=(x,y,z) is the Cartesian coordinate system that defines the axial (x), pitchwise (y), and spanwise (z) directions, and t denotes time. The nondimensional groups in the governing Eqs. (3) and (4) include the inlet Reynolds number, Re(ρinUinc)/νin and the inlet Mach number, MUin/γRTin. The (molecular) Prandtl number, Pr(cpμin/k), is assumed to be constant and is set equal to 0.72. Here, νin is the inlet kinematic viscosity, R is the specific gas constant, cp is the specific heat, and k is the fluid thermal conductivity. The fluid is assumed to behave as a calorifically perfect gas and obeys the ideal gas law, which, after nondimensionalization, can be written as p=(ρT)/(γM2), with γ=1.4. The viscous stress tensor and the heat flux vector are defined as τijμ(uixj+ujxi23ukxkδij) and qjμ(γ1)M2RePrTxj, where μ(x,t) and T(x,t) are the instantaneous molecular viscosity and temperature, respectively, where the former variable is computed using Sutherland’s law based on a (dimensional) reference temperature of 300 K.

HiPSTAR solves Eqs. (2)(4) on structured overlapping grids using the fourth-order accurate compact finite differences for spatial discretization and an ultra-low storage frequency-optimized explicit fourth-order accurate Runge–Kutta method for time integration [19]. The computational domain is shown in Fig. 3(a) and is composed of three overlapping blocks: (i) a background Cartesian H-type grid (block 1); (ii) an outer curvilinear O-type grid (block 2), and (iii) an inner curvilinear O-type grid (block 3). At the overlapping block boundaries, continuity conditions are imposed as variables and are interpolated with a fourth-order Lagrangian method [20]. A zoomed-in view of the inner and outer O-meshes is shown in Fig. 3(b), which shows the typical mesh resolution around the roughness profile. Note that only every fourth in-plane point is shown.

Fig. 3
(a) Computational domain. Block 1: Cartesian H-grid; block 2: outer O-mesh and block 3: inner O-mesh. Note that only every 16th in-plane mesh point is shown. (b) Typical mesh resolution around the roughness profile, h(x,y) (black line), at z/c=0 in the leading edge region. Note that only every 4th in-plane mesh point is shown.
Fig. 3
(a) Computational domain. Block 1: Cartesian H-grid; block 2: outer O-mesh and block 3: inner O-mesh. Note that only every 16th in-plane mesh point is shown. (b) Typical mesh resolution around the roughness profile, h(x,y) (black line), at z/c=0 in the leading edge region. Note that only every 4th in-plane mesh point is shown.
Close modal

To represent the geometry of each rough HPT vane, a point cloud of the three-dimensional blade coordinates (Eq. (1)) was immersed within the inner O-grid of the computational domain (Fig. 3), and a compressible variant of the boundary data immersion method (BDIM) [21] was used to enforce no-slip boundary conditions on velocity (ui=0) and an isothermal boundary condition on temperature (T=Tw, where Tw is blade surface temperature).

To drive the flow through the passage, an outlet static pressure was prescribed and an inlet static pressure and temperature were imposed as Riemann conditions using standard isentropic flow relations. To mimic disturbances from an upstream combustor, synthetic eddies with a specified inlet turbulent intensity, Tu/Uin, and an (isotropic) integral length scale, L/c, were fed into the domain inlet using a compressible version of the digital filter technique [22]. A nonreflective zonal characteristic boundary condition was enforced at the outlet [23]. Periodic boundaries were imposed in the pitchwise (y) and spanwise (z) directions to simulate a linear HPT cascade without end-walls.

2.3 Description of Simulations.

As shown in Table 2, a total of seven different rough HPT simulations were undertaken and can be divided into two separate groups. Group 1 contains four cases where the axial extent of leading edge roughness was varied, while the inlet turbulence intensity and integral length scale constant were held fixed, along with a reference smooth-vane case (with matched inflow conditions) from the past studies [10,11]. Group 2 contains four cases (including one case from Group 1) and adopts the opposite approach, i.e., the inlet condition was varied for a fixed roughness topography. The seven cases investigated here are named using the identification code:
(5)
where the first three characters correspond to the axial position of control point IV of the roughness weighting function, e.g., X28 corresponds to x/c=0.28 (see Fig. 1 or Table 1), and the last four digits represent the values of the inlet turbulence intensity and integral length scale, e.g., 0820 corresponds to Tu/c=8% and L/c=20%, and so on. For cases X28_0820 and X28_2020, the width of the domain, Lz, had to be doubled from 0.4c to 0.8c to accommodate the wider length scale inflow disturbances, which doubled the total number of grid points from 7.3 billion to 14.6 billion. For all simulations, the Reynolds number and Mach number based on exit conditions were 590,000 and 0.92, respectively.
Table 2

Table of simulations

GroupCaseTu/Uin(%)L/c(%)Lz/c(%)(Nx×Ny×Nz)×109
Smooth8.08.00.43.8
1X28_08088.08.00.47.3
X40_08088.08.00.47.3
X50_08088.08.00.47.3
X60_08088.08.00.47.3
2X28_08088.08.00.47.3
X28_200820.08.00.47.3
X28_08208.020.00.814.6
X28_202020.020.00.814.6
GroupCaseTu/Uin(%)L/c(%)Lz/c(%)(Nx×Ny×Nz)×109
Smooth8.08.00.43.8
1X28_08088.08.00.47.3
X40_08088.08.00.47.3
X50_08088.08.00.47.3
X60_08088.08.00.47.3
2X28_08088.08.00.47.3
X28_200820.08.00.47.3
X28_08208.020.00.814.6
X28_202020.020.00.814.6

For groups 1 and 2, the background mesh (block 1) has dimensions of (Nx,Ny,Nz)=(1470,716,576) and the outer O-mesh (block 2) has dimensions of (Nx,Ny,Nz)=(8165,237,576), where Nx, Ny, and Nz are the number of grid points along the axial (x), pitchwise (y), and spanwise (z) directions, respectively. The inner O-mesh (block 3) has dimensions of (Nx,Ny,Nz)=(35595,230,1200), except for cases X28_0820 and X28_2020 in group 2, where (Nx,Ny,Nz)=(35595,230,2400) to account for the wider computational domain. For all cases presented here, the viscous-scaled grid resolutions in the blade tangential direction (s), blade normal direction (n), and spanwise direction (z) are such that Δs+<5.0, Δn+<1.5, and Δz+<8.0, where superscript + denotes nondimensional by the local (mean) viscous length scale. The viscous-scaled grid resolution adopted here is therefore commensurate with that in past work related to the present study [10,11]. The accuracy and reliability of the HiPSTAR flow solver and BDIM algorithm in simulating flow over rough (and smooth) HPTs have been validated extensively in the previous works related to the present study [1012].

For each rough HPT case listed in Table 2, a series of extreme-scale compressible DNSs were performed on the Summit supercomputer based in the Oak Ridge National Laboratory, USA. Where possible, all seven cases were executed concurrently as a single job using 1152 Summit nodes, meaning that 6912 Nvidia Tesla V100 GPUs were required for each production-scale calculation. In terms of computational resources, each production scale simulation on Summit required approximately 168 h of wall-time to reach statistical convergence, meaning that 194,000 node-hours (or 1116,000 GPU-hours) were consumed in total. Some calculations were also performed on Summit’s successor, Frontier, using a comparable number of AMD Instinct M1250x Accelerators. For a typical job with 5.0×106 points per GPU, Frontier delivered a speed-up of three times compared to Summit. Herein, all statistical quantities were collected for a minimum of four nondimensional time units (based on the axial inlet velocity and axial blade chord) following an initial transient.

3 Results

3.1 Drag Partition Between Viscous and Pressure Forces.

Traditionally, roughness effects have been quantified in fully developed turbulent flows in canonical configurations, e.g., circular pipes [24], plane channels [25], and flat-plate boundary layers [26], often with the goal of mapping skin friction measurements to an equivalent value of Nikuradse’s sandgrain roughness, ks [27]. However, ks is only strictly valid for turbulent flows in the fully rough regime, i.e., a rough-wall flow where the skin friction coefficient is Reynolds number independent. Whether the fully rough regime is reached in the context of the present investigation is debatable, since detailed accounts of how skin friction drag varies as a function of the Reynolds number is typically limited to smooth HPT vanes [8,13,28]. Furthermore, despite having units of length, it is important to note that ks cannot be reliably determined from topographical data alone (e.g., from a digital scan of an ablated turbine blade [2]). This is because ks is a hydraulic length scale, i.e., a hydrodynamic property of the flow, not a topographical property of the surface, that must be determined on a surface-specific basis. In fact, as noted in a recent review article [29], determining ks remains the major bottleneck for accurate predictions of roughness effects in fully rough flows. Whether this bottleneck applies to the present study raises the following question: do any of the rough HPT vanes under investigation here attain fully rough conditions?

To answer this question, the fractional contribution of pressure drag, Fp, and viscous drag, Fν, to the total drag, Ftot=Fp+Fν, can be computed. To reach the fully rough regime, it is generally accepted that the tangential component of the (inviscid) pressure force acting on the roughness elements overwhelms its viscous counterpart, i.e., FpFν [30]. A previous experimental investigation of surface roughness effects on the fully developed turbulent pipe flow [31] suggests that the onset of the fully rough regime occurs when pressure drag accounts for at least 80% of the total drag force, i.e., Fp/Ftot0.80. Whether this threshold is exceeded in the context of the present study remains unclear, since very little is known regarding the relative magnitude of viscous and pressure forces acting against rough HPT vanes—with one recent exception being the work presented in Ref. [11]. Yet, knowledge of the drag partition over rough HPT vanes is of critical importance, especially when one considers that any ks-based modeling approach, e.g., surface-specific wall functions [32] or roughness-induced transition models [33], implicitly assume that fully rough conditions are attained.

To draw a fair comparison between the viscous and pressure forces exerted against the surface of the rough and smooth vanes, it is necessary to use a common tangential direction along which to quantify the drag force. Here, we follow [10,11] and take the reference as the tangential unit vector in the x-y plane for the smooth blade, txys, which defines the mean viscous force as follows:
(6)
and the mean pressure force as follows:
(7)
where Tint is the sampling period and t0 is the initial time. Considering Eqs. (6) and (7), it follows that Fp=0 for the smooth blade since its surface normal vector and surface tangential vector are orthogonal by definition, i.e., ntxys=0. On the other hand, the dot product between the normal vector on the rough surface, n, and smooth-vane tangential unit vector, txys, is typically nonzero, and, as a result, Fp0 for the rough vane.

The fractional contribution of pressure drag to the total drag, Fp/Ftot, is plotted against the axial position for groups 1 and 2 in Fig. 4. Here, 1<x/c<0 corresponds to the pressure surface, 0<x/c<1 corresponds to the suction surface, and the horizontal dashed line depicts the threshold of Fp/Ftot=0.8, which is required to attain fully rough conditions according to Ref. [31]. At first glance, the data in Fig. 4 bear a striking resemblance to the shape of the weighting function used to control the roughness amplitude around the vane (see Fig. 1). To be specific, group 1 data (Fig. 4(a)) show that the fractional contribution of pressure drag correlates with regions of high-amplitude roughness, i.e., in the forward half of the suction surface (0x/c0.5) and along the majority of the pressure surface (0.9x/c0.0). A simple explanation for this behavior is that the larger amplitude asperities have more surface normals with an axial component on which pressure drag can act (see integrand of Eq. (7)). This also explains why increasing the axial extent of larger amplitude roughness results in larger values of Fp/Ftot along a greater proportion of the vane suction surface. Quantitatively, the fractional contribution of pressure drag to the total drag reaches a peak value of Fp/Ftot0.6 at x/c0.25 and x/c0.85 on the suction and pressure surfaces, respectively, meaning that viscous drag still accounts for up to two fifths of the total drag, since Fν/Ftot1.0Fp/Ftot. On the other hand, in the aft portion of the suction surface, where the RMS roughness amplitude is four times smaller than elsewhere on the vane (see Fig. 1 and Table 1),– pressure drag accounts for just one tenth of the total drag force, i.e., Fp/Ftot0.1. Looking at Fig. 4(b), the variation of Fp/Ftot with respect to x/c is almost identical for all cases in Group 2, suggesting that pressure drag is relatively insensitive to changes in inflow conditions, although this is not necessarily the case since the magnitude of a force cannot be inferred from a ratio of forces. Nevertheless, beyond an axial distance on the suction surface of, say, x/c0.6, the viscous drag force for all cases in groups 1 and 2 is approximately nine times greater than the pressure force, i.e., FνFp, which is the opposite of what is required to reach fully rough conditions.

Fig. 4
Fp/Ftot versus x/c. (a) Effect of varying roughness and (b) effect of varying inlet conditions. The onset of the fully rough regime Fp/Ftot=0.8 [31] is also shown (dashed line).
Fig. 4
Fp/Ftot versus x/c. (a) Effect of varying roughness and (b) effect of varying inlet conditions. The onset of the fully rough regime Fp/Ftot=0.8 [31] is also shown (dashed line).
Close modal

Overall, the data in Fig. 4 imply that fully rough conditions are not achieved for any of the HPT vanes under investigation here. The fractional contribution of pressure drag to the total drag fails to exceed the threshold Fp/Ftot=0.8 of Ref. [31] at any point on the pressure or suction surface. As a result, the bottleneck of determining an equivalent value of Nikuradse’s sandgrain roughness [29] does not apply here, since it is not possible to ascribe a physically meaningful value (or values) of ks. The lack of fully rough conditions is most apparent in the aft portion of the suction surface, where viscous drag accounts for up to 90% of the total drag. As will be shown later, the flow in this region perceives the blade as an aerodynamically (and thermodynamically) smooth surface. Consequently, the Reynolds number scaling of the local skin friction and heat transfer coefficients agree well with the existing empirical correlations for (turbulent) flow over smooth surfaces.

3.2 Skin Friction Drag and Wall Heat Flux.

Attention is now turned toward evaluating HPT aerothermal performance. This is achieved by quantifying the sensitivity of the skin friction coefficient and wall heat flux with respect to changes in the axial extent of leading edge roughness (group 1) and changes in inlet conditions (group 2).

The instantaneous skin friction coefficient is defined here as follows:
(8)
where subscript w denotes an on-surface quantity acquired using the data collection framework devised in Ref. [10]. Physically, this definition of Cf represents the ratio of instantaneous wall shear stress (i.e., the sum of viscous and pressure forces, τwFp+Fν) at a given point on the blade surface relative to the (mean) dynamic pressure at the inlet.

Snapshots of the instantaneous skin friction coefficient are shown on the suction and pressure surfaces in Fig. 5. Isovolumes of the Q-criterion [34] are also shown over the right-hand half of each vane. Looking at the data, it is clear that the suction-side (SS) Cf distributions exhibit a greater sensitivity with respect to changes in either leading edge roughness (Figs. 5(a) and 5(b)) or inlet conditions (Figs. 5(c) and 5(d)) compared to their pressure-side (PS) counterparts. Comparing the suction-side data for both groups, it is also clear that Cf is more sensitive to changes in the axial extent of leading edge roughness (Fig. 5(a)) than variations of inlet turbulence (Fig. 5(c)), at least for the seven cases considered here. For instance, although the Cf distributions and vortical structures for cases X28_0808 and X40_0808 show that the boundary layer breaks down into turbulence in the final 15% of the suction surface (Fig. 5(a)), the remaining group 1 cases show a different behavior. In particular, case X60_0808 exhibits a streaky Cf distribution and fine-scale vortical motions across its span and over most of the suction surface—flow characteristics that are synonymous with wall turbulence [35]. Overall, the data in Fig. 5 show that even a small change in the axial extend of leading edge roughness can trigger a large change in the suction-side boundary layer state, underlining the critical role that localized surface roughness plays in determining HPT aerothermal performance.

Fig. 5
Instantaneous skin friction coefficient, Cf, and Q-criterion: (a) and (b) effect of varying roughness and (c) and (d) effect of varying inlet conditions
Fig. 5
Instantaneous skin friction coefficient, Cf, and Q-criterion: (a) and (b) effect of varying roughness and (c) and (d) effect of varying inlet conditions
Close modal

Next, the mean skin friction coefficient, C¯f, was obtained by integrating the instantaneous Cf values (Eq. (8)) with respect to time and the spanwise direction. The variation of C¯f along the vane suction surface for groups 1 and 2 is plotted against the axial position in Fig. 6. Here, the gray and black lines represent the C¯f values corresponding to the smooth and rough vanes vane, respectively. The light and dark contours represent the Cf histogram for the smooth and rough vanes, respectively, and illustrate the frequency distribution of instantaneous Cf values at a given x/c location. Note that a vertical shift of 0.1 units has been applied to consecutive cases in groups 1 and 2 for brevity. Looking at the data in Fig. 6, it is clear that these statistical data support their instantaneous counterparts shown in Fig. 5 in the sense that: (i) Increasing the axial extent of leading edge roughness causes an abrupt upstream shift of the (mean) transition point, e.g., compare the C¯f curve for case X60_0808 against that of X40_0808, X28_0808, or the reference smooth vane. (ii) Compared to the effect of increasing the axial extent of leading edge roughness, increasing the inlet turbulence intensity or integral length scale has a relatively weaker effect upon the (mean) transition point, i.e., cases X28_0808 and X28_0820 closely resemble the smooth-vane data, the upward kink in the C¯f curves for cases X28_2008 and X28_2020 shows a more gradual upstream shift. Further understanding of the near-wall boundary layer state can be gained from comparing the shape of the Cf histograms included in Fig. 6. Focusing on the data in the trailing edge region, say, 0.7x/c1.0, it is clear that the presence of intermittent flow phenomena across the span (and in time) produce a wider distribution of instantaneous Cf values, compared to those that break down to turbulence closer to the vane leading edge. For instance, comparing the Cf histogram for case X60_0808 with that for X50_0808 or X28_2020, it is obvious that the broadband distribution of Cf values for the latter cases is driven by a rapidly varying boundary layer state due to the inception, growth, and coalescence of turbulent spots along the suction surface—the foot print of these structures are clearly visible on the suction surface of each vane (see Figs. 5(a) and 5(c)). Similar observations have been made in the context of smooth HPT vanes, e.g.. see Ref. [36].

Fig. 6
C¯f versus x/c on the vane suction surface: (a) effect of varying roughness and (b) effect of varying inlet conditions
Fig. 6
C¯f versus x/c on the vane suction surface: (a) effect of varying roughness and (b) effect of varying inlet conditions
Close modal
The instantaneous (nondimensional) wall heat flux is defined here as follows:
(9)

It also represents the product of the thermal conductivity and the temperature gradient at the wall in the direction perpendicular to the blade surface. Snapshots of the instantaneous wall heat flux on the vane suction and pressure surfaces are shown in Fig. 7. Isovolumes of the Q-criterion are once again included. In general, the spatial distribution of instantaneous heat flux data resembles their instantaneous Cf counterparts in Fig. 5, on both the pressure and suction surfaces. As expected, regions of high instantaneous heat flux correlate with regions of high instantaneous skin friction drag, and the presence of fine-scale vortical motions in the near-wall region enhance the rate at which heat is transferred from the hot fluid into the (relatively) cool blade. Next, the mean wall heat flux, Q¯, was obtained by integrating the instantaneous Q values (Eq. (9)) with respect to time and the spanwise direction. The Q¯-distribution along the vane suction surface is plotted against the axial position for groups 1 and 2 in Fig. 8. Note that a vertical shift of 0.5 units has been applied to consecutive cases in groups 1 and 2 for brevity. Again, the gray and black lines represent the Q¯ values corresponding to the smooth and rough vanes vane, respectively. The light and dark contours represent the Q-histogram for the smooth and rough vanes, respectively. Data points from a past study of a smooth HPT vane resolved using a body-fitted mesh [28] are also included (white squares) for comparison. Looking at the data shown in Fig. 8, these Q¯-curves support their instantaneous counterparts (Fig. 7) in the sense that: (i) increasing the axial extent of leading edge roughness can cause an abrupt change in the mean transition point, e.g., compare the Q¯-curves for cases X40_0808 and X50_0808 in the range of 0.5x/c1.0 in Fig. 8; (ii) changing the inlet condition also has a strong effect upon the levels of wall heat flux, particularly in the trailing edge region; and (iii) the presence of transitional flow phenomena (e.g., see localized “hot spots” on the suction surface of cases X50_0808 and X28_2020 in Fig. 7) results in the widest band of realizations of wall heat flux in their respective Q-histograms.

Fig. 7
Instantaneous wall heat flux, Q, and Q-criterion: (a) and (b) effect of varying roughness and (c) and (d) effect of varying inlet conditions
Fig. 7
Instantaneous wall heat flux, Q, and Q-criterion: (a) and (b) effect of varying roughness and (c) and (d) effect of varying inlet conditions
Close modal
Fig. 8
Q¯ versus x/c on the vane suction surface: (a) effect of varying roughness and (b) effect of varying inlet conditions
Fig. 8
Q¯ versus x/c on the vane suction surface: (a) effect of varying roughness and (b) effect of varying inlet conditions
Close modal

3.3 Augmentation of Heat Transfer and Drag.

The preceding subsection demonstrated that changes in the axial extent of leading edge roughness and inlet conditions can profoundly impact skin friction and wall heat flux (in both an instantaneous and averaged sense), particularly along the vane suction surface (see Figs. 58). However, what is not yet clear is how heat transfer and drag are augmented compared to smooth vane levels or to some other suitable reference case.

To quantify the overall augmentation of heat transfer and drag, the ratio of the mean Nusselt number, Nu¯/Nu¯0, and the mean skin friction coefficient, C¯f/C¯f,0, were integrated over the suction surface and the pressure surface to obtain “global mean” quantities, herein denoted using angular brackets, i.e., Nu/Nu0 and Cf/Cf0, where subscript 0 denotes a specified reference state. Group 1 reference state is taken as the smooth-vane case, whereas Group 2 reference state is taken as case X28_0808. Hence, any augmentation in heat transfer or drag in group 1 is due to the effect of varying the axial extent of leading edge surface roughness (for a fixed inlet condition), whereas changes in group 2 are due to the effect of varying the inlet conditions (for a fixed roughness topography). If Cf/Cf0>Nu/Nu0, then the fractional increase in heat transfer comes at the expense of a greater fractional increase in drag compared to the reference state, whereas Cf/Cf0<Nu/Nu0 indicates the opposite. Otherwise, if Cf/Cf0=Nu/Nu0, then the fractional increase in heat transfer and drag matches those of the reference state. Note the (instantaneous) Nusselt number is defined as Nu(xw,zw,t)(c/ΔT)/T(x,z,t)n|w, which represents the ratio of convective to conductive heat transfer at a given point on the blade surface. Here, thermal conduction is driven by the difference between the (constant) total inlet temperature and the (isothermal) blade surface temperature, i.e., ΔT(TintotTw)>0, whereas thermal convection is solely determined by the instantaneous temperature gradient at the wall in the direction perpendicular to the blade surface.

The relationship between Nu/Nu0 and Cf/Cf0 for cases in groups 1 and 2 is shown in Fig. 9, a plot that is inspired from the past work of Bunker, e.g., see Fig. 3 in [37]. Looking at the data, it is clear that all of the suction-side data points (squares) fall below the diagonal dashed line, i.e., Cf/Cf0=Nu/Nu0, meaning that the fractional increase in (global mean) heat transfer comes at the expense of an even greater fractional increase in (global mean) drag, i.e., Cf/Cf0>Nu/Nu0. Data points from Ref. [11] that quantified the aerothermal performance of HPT vanes covered with uniformly distributed near-Gaussian roughness and the same inlet condition as group 1 (see Table 2) is included in Fig. 9 (white squares) and show a similar trend. The same behavior has been reported in the past studies of fully developed channel-flow turbulence over irregular roughness [32] and streamwise-aligned riblets [38], and as result, the breakdown in Reynolds analogy due to roughness observed here in some sense is not surprising. On the pressure surface, the data points (circles) fall closer to the diagonal dashed line, indicating a relatively small fractional increase in heat transfer and drag compared to the specified reference state, i.e., the smooth vane for group 1 and case X28_2020 for group 2.

Fig. 9
⟨Nu⟩/⟨Nu⟩0 versus ⟨Cf⟩/⟨Cf⟩0. (a) Effect of varying roughness. (b) Effect of varying inlet conditions. Circular and square symbols correspond to SS and PS data, respectively. Data obtained from Nardini et al. [11] are shown as white symbols.
Fig. 9
⟨Nu⟩/⟨Nu⟩0 versus ⟨Cf⟩/⟨Cf⟩0. (a) Effect of varying roughness. (b) Effect of varying inlet conditions. Circular and square symbols correspond to SS and PS data, respectively. Data obtained from Nardini et al. [11] are shown as white symbols.
Close modal

To understand how the augmentation of heat transfer and drag vary along the pressure and suction surfaces, the ratios of the mean skin friction coefficient, C¯f/C¯f,0, and the mean Nusselt number, Nu¯/Nu¯0, for groups 1 and 2 are plotted against the axial position in Fig. 10. Note that a vertical offset of 10 units has been applied to the former ratio for clarity. Looking from left-to-right across Fig. 10(a), it is clear that the fractional increases in heat transfer and drag are (approximately) in balance across the entire pressure side and for axial positions less than x/c0.4 on the suction side. It is also clear that changes in the axial extent of leading edge roughness can lead to very significant augmentations of heat transfer and drag in the aft portion of the suction surface. To be specific, at an axial position of x/c0.90, the mean skin friction drag for cases X50_0808 and X60_0808 is augmented by up to a factor of 22 compared to smooth-vane levels, whereas Nu¯ is augmented by just a factor of six at the same location. Quantitatively, the mean Nusselt number increases from Nu¯088 to Nu¯527 at x/c0.90 on the smooth vane and case X50_0808, respectively. The remaining group 1 cases exhibit a much weaker augmentation of heat transfer and drag—consistent with the previous analysis of the skin friction coefficient and wall heat flux, e.g., see Figs. 5(a)5(c) and Figs. 7(a)7(c). Looking from left-to-right across Fig. 10(b), it is clear that changes in inlet turbulence can also lead to appreciable augmentations of trailing edge heat transfer and drag, although the fractional increases are less significant compared to the group 1 data. Nevertheless, at an axial position of x/c=0.90, skin friction drag and convective heat transfer are, respectively, augmented by a factor of 12 and 4 for case X28_2020 compared to the reference rough-vane state (case X28_0808). The integrated effect of these large, local augmentations of heat transfer and drag in the trailing-edge region explain why the global mean ratios, Nu/Nu0 and Cf/Cf0 fall below the diagonal dashed line in Fig. 9. Compared to the suction surface, the peak augmentation of heat transfer and drag on the pressure surface is at least one order of magnitude smaller for all cases in groups 1 and 2, which explains why the global mean ratios closely follow the diagonal dashed line in Fig. 9.

Fig. 10
Cf¯/Cf¯0 and Nu¯/Nu¯0, versus x/c: (a) effect of varying roughness and (b) effect of varying inlet conditions
Fig. 10
Cf¯/Cf¯0 and Nu¯/Nu¯0, versus x/c: (a) effect of varying roughness and (b) effect of varying inlet conditions
Close modal

To illustrate the differences in boundary layer state responsible for the augmentation of drag and heat transfer, profiles of the mean blade tangential velocity and mean temperature are plotted against blade normal distance at an axial position of x/c=0.90 in Fig. 11. Here, velocity is normalized by its corresponding mean value at the boundary layer edge, i.e., U¯t/U¯e, temperature is shown in “reduced” form, i.e., (T¯T¯w)/(T¯eT¯w), where T¯e is the mean temperature at the edge of the thermal boundary layer, and blade normal distance is normalized by the local boundary layer thickness, i.e., n/δe. Looking at group 1 data first (Fig. 11(a)), it is clear that the thermal and hydrodynamic boundary layer slopes at the wall are far steeper for cases X50_0808 and X60_0808 compared to reference smooth-vane profile, which explains why skin friction and heat transfer are augmented so dramatically for these two particular cases (see Fig. 10(a)). Similar behavior for cases X28_0820 and X28_0820 (relative to case X28_0808) can be observed in Fig. 11(b). Also included in Fig. 11 are DNS data points for the mean streamwise velocity profile of an incompressible zero-pressure gradient (ZPG) turbulent boundary layer (TBL) over a smooth flat plate [39] at a momentum thickness Reynolds number of Reθ1400. Remarkably, cases X50_0808 and X60_0808 show good levels of agreement with the DNS data points at all blade normal positions, whereas X28_2020 shows the best overall agreement among group 2 cases. For the remaining cases in groups 1 and 2, the blade normal velocity profiles show a point of inflection at around n/δe0.3 (a feature that is absent in the corresponding temperature profiles), which is consistent with the sharp drop in C¯f over the range 0.85x/c0.95 (see Fig. 6) and indicative of a flow approaching a (mean) separation as a result of adverse pressure gradient forcing.

Fig. 11
Profiles of U¯t/U¯e and (T¯−T¯w)/(T¯e−T¯w) versus n/δe at x/c=0.90 on the vane suction surface: (a) effect of varying roughness and (b) effect of varying inlet conditions. Data from the study by Schlatter and Örlü [39] are also included (white squares).
Fig. 11
Profiles of U¯t/U¯e and (T¯−T¯w)/(T¯e−T¯w) versus n/δe at x/c=0.90 on the vane suction surface: (a) effect of varying roughness and (b) effect of varying inlet conditions. Data from the study by Schlatter and Örlü [39] are also included (white squares).
Close modal

To further characterize the velocity and temperature profiles plotted in Fig. 11, some nondimensional boundary layer parameters are given in Table 3. The parameters include the following: the friction Reynolds number, ReτUτδe/νw; the momentum thickness Reynolds number, ReθUeθ/νw; the boundary layer shape factor, Hδ*/θ; the mean skin friction coefficient, C¯f; the mean Stanton number, C¯hNu¯/(RePr); and twice the mean Stanton number divided by the mean skin friction coefficient, i.e., the Reynolds analogy factor, 2C¯h/C¯f. Looking at Table 3, stark differences between the various cases can be observed. For instance, although the smooth-vane friction Reynolds number reaches a value of Reτ60, this quantity is a factor of nine greater for case X60_0808. Likewise, although the smooth-vane shape factor attains a value of H3.0, a value of H1.4 is attained for case X60_0808—the latter typical of turbulent flow and the former for a near-separated flow. In terms the Reynolds’ analogy, which, according to Bons [40] and others, is preserved when 2C¯h/C¯f1, the data in Table 3 indicate that this approximation holds, but only for three out of the seven cases investigated here (X50_0808, X60_0808, and X20_2020). Two common features among these three cases are their relatively high friction Reynolds numbers (Reτ400) and relatively low shape factors (H1.4), indicating the flow has attained fully turbulent conditions at x/c=0.90. A further common feature among these cases (and all other cases) is that the fractional contribution of pressure drag to the total drag is far smaller than the viscous drag force (see Fig. 4). As a result, one may conclude that the preservation of Reynolds’ analogy, i.e., 2C¯h/C¯f1, for cases X50_0808 and X60_0808 from group 1 and case X20_2020 from group 2 is due to the presence of fully turbulent flow and the absence of fully rough conditions. Note that the Reynolds’ analogy would likely not be preserved if fully rough conditions were reached. This is because although C¯f would asymptote a surface-specific constant, C¯h would continue to decrease with the increasing Reynolds number and, as a result, the Reynolds’ analogy factor 2C¯h/C¯f<1, e.g., see rough-wall channel-flow studies of Refs. [32,41]. Put in other words, the breakdown of the Reynolds’ analogy in a fully rough flow is due to an absence of viscous-dominated transport in the near-wall region—the opposite of what is observed here (see Fig. 4)

Table 3

Boundary layer parameters evaluated at x/c=0.90 on the vane suction surface

CaseReτReθHC¯f×103C¯h×1032C¯hC¯f
Smooth60.1649.12.980.1850.3473.76
X28_080881.5687.03.050.3200.4002.49
X40_0808105.2727.92.750.4610.4892.12
X50_0808483.21236.41.453.9912.1041.05
X60_0808529.61485.21.473.7211.9981.07
X28_080881.5687.03.050.3200.4002.49
X28_2008250.9854.41.951.9361.2341.28
X28_082083.7704.62.990.3280.4052.47
X28_2020422.6969.21.483.7341.9921.07
TBL [39]492.21420.91.433.880
CaseReτReθHC¯f×103C¯h×1032C¯hC¯f
Smooth60.1649.12.980.1850.3473.76
X28_080881.5687.03.050.3200.4002.49
X40_0808105.2727.92.750.4610.4892.12
X50_0808483.21236.41.453.9912.1041.05
X60_0808529.61485.21.473.7211.9981.07
X28_080881.5687.03.050.3200.4002.49
X28_2008250.9854.41.951.9361.2341.28
X28_082083.7704.62.990.3280.4052.47
X28_2020422.6969.21.483.7341.9921.07
TBL [39]492.21420.91.433.880

3.4 Reynolds Number Scaling of Drag and Heat Transfer.

To understand the Reynolds number dependence of heat transfer and drag, HPT Moody diagrams can be constructed. Here, the Re-scaling is examined in the context of the Chilton-Colburn analogy, which states that the (mean) friction factor coincides with the product of the (mean) Stanton number and two-third power of the (molecular) Prandtl number, i.e., C¯f/2=C¯hPr2/3. While HPT Moody diagrams have been constructed using measurements acquired at mid-span over smooth HPT vanes, e.g., see Figs. 7 and 8 in the study by Thole et al. [8], detailed information in the context of rough vanes remains rare.

The Chilton–Colburn form of the mean skin friction coefficient based on the local dynamic pressure at the boundary layer edge, C¯fτw/(ρ¯eU¯e2/2), and mean Stanton number, C¯h, are plotted against the momentum thickness Reynolds number, Reθ, on the suction surface in Fig. 12. Also included on this plot are the Coles–Fernholz correlation [42] for the mean skin friction coefficient of an incompressible smooth-wall ZPG TBL, C¯f=2(κ1log(Reθ)+C)2 (black dashed line), with von Kármán coefficient of κ=0.387 and offset of C=4.1, as well as TBL data points on the range 670Reθ2000 [39]. Looking again at Fig. 12, it is clear that all of the HPT data points (in both groups 1 and 2) collapse onto the Coles–Fernholz correlation and the DNS data beyond Reθ1000. This behavior not only indicates that the (viscous) skin friction losses scale in a manner similar to that of a smooth-wall ZPG TBL but also indicates that the mean Stanton number can be approximated as C¯h4(κ1log(Reθ)+C)2Pr2/3 provided that the momentum thickness Reynolds number exceeds Reθ1000. The collapse of the HPT data on to the Coles–Fernholz correlation occurs at increasingly small Reθ values for cases X50_0808 (Reθ800) and X60_0808 (Reθ500) since increasing the axial extent of leading edge roughness causes the transition point to shift upstream (see Figs. 5 and 7). Ultimately, earlier transition leads to a larger value of Reθ at the vane trailing edge since more (turbulent) momentum loss is accumulated along a greater proportion of the suction surface, relative to the other cases. Increasing the axial extent of leading edge roughness also causes the C¯f and C¯h curves to shift upward (relative to the smooth-vane level) on the range 400Reθ700, which is consistent with the remarks of Denton [43], who predicted greater losses in the transition region at higher values of roughness. The C¯f curves for cases X50_0808 and X60_0808 also highlight how an early transition produces a boundary layer that is more resistant to the suction surface adverse pressure gradient, whereas C¯f curves for the other (nonturbulent) cases exhibit a rapid reduction within range 400Reθ700, before finally transitioning to turbulence at around Reθ1,000. A similar behavior is evident in the HPT Moody diagrams constructed in a past experimental study by Thole et al. [8], although those authors preferred to examine the C¯f/2=C¯hPr2/3 scaling as a function of Reynolds number based on distance from the leading edge of a smooth vane, as opposed to Reθ scaling presented here.

Fig. 12
C¯f/2, C¯hPr2/3 versus Rθ on the vane suction surface: (a) effect of varying roughness and (b) effect of varying inlet conditions
Fig. 12
C¯f/2, C¯hPr2/3 versus Rθ on the vane suction surface: (a) effect of varying roughness and (b) effect of varying inlet conditions
Close modal

4 Conclusions

The impact of localized, non-Gaussian surface roughness on HPT aerothermal performance was quantified at an exit Reynolds number of 590,000 and exit Mach number of 0.92. This was achieved by undertaking a series of extreme-scale compressible DNSs for a range of HPT vanes covered with systematically varied surface roughness for a range of realistic turbulence inflow conditions using up to 14.6 billion grid points per simulation. The key results of this study are as follows:

  • Drag partition between pressure and viscous forces. The present study shows that the fractional contribution of pressure drag to the total drag falls short of the fully rough regime. In the aft portion of the suction surface, the (inviscid) pressure force is nine times smaller than its viscous counterpart, and, as a result, the (turbulent) near-wall flow perceives the blade surface as an aerodynamically (and thermodynamically) smooth surface. As a result, roughness effects cannot be modeled using an equivalent value of Nikuradse’s sandgrain roughness, since ks is only meaningful if pressure drag overwhelms viscous drag, i.e., the opposite of what is observed here.

  • Augmentation of heat transfer and drag. While the augmentation of convective heat transfer comes at the expense of an even greater augmentation in drag—in both the global and local sense—significant augmentations only occur on the suction surface. For one of the rough HPT vanes investigated here, the local skin friction coefficient is a factor of 22 greater than the reference smooth-vane level, whereas the local Nusselt number is just six times higher in comparison. The rise in skin friction drag and heat transfer occurs as a direct result of contrasting boundary layer states over the smooth and rough vanes, which are driven by differing transition mechanisms and vortical structures along the suction surface.

  • Reynolds’ analogy and Reynolds number scaling. The present investigation demonstrates that the Reynolds’ analogy, i.e., 2C¯h/C¯f1, can hold over the HPT vanes under investigation here, but only in regions (i) where the flow is fully turbulent and (ii) where viscous drag dominates pressure drag. Here, conditions (i) and (ii) are satisfied for three cases that undergo the earliest transition to turbulence on the vane suction surface. The Re scaling of heat transfer and drag was examined in the context of the Chilton–Colburn analogy. The HPT Moody diagrams demonstrate that the C¯f/2=C¯hPr2/3 scaling holds, but only for regions of turbulent flow where the momentum thickness Reynolds number was sufficiently large, i.e., Reθ1000.

Acknowledgment

We are grateful for the permission of GE Aerospace to publish results from this study. Support from the ARC is acknowledged. This research used data generated on the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC05-00OR22725.

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

c =

axial blade chord (m)

h =

roughness height (m)

k =

thermal conductivity (W/m/K)

n =

blade normal distance (m)

p =

pressure (kg/m/s2)

s =

blade tangential distance (m)

t =

time (s)

x =

axial distance (m)

y =

pitchwise distance (m)

z =

spanwise distance (m)

n =

surface normal unit vector (–)

A =

planform area (m2)

C =

log-law offset (–)

E =

total energy (J)

Q =

wall heat flux (W/m2)

R =

specific gas constant (J/K/mol)

T =

temperature (K)

L =

turbulence length scale (m)

cp =

specific heat capacity (J/K/kg)

ks =

equivalent sandgrain height (m)

qj =

heat flux vector (W/m2)

t0 =

initial time (s)

ui =

velocity (m/s)

Fp =

pressure drag force (kgm/s2)

Fν =

viscous drag force (kgm/s2)

Ftot =

total drag force (kgm/s2)

Lz =

spanwise domain width (m)

Nx =

number of axial grid points (–)

Ny =

number of pitchwise grid points (–)

Nz =

number of spanwise grid points (–)

Rq =

root-mean-square roughness height (m)

Tin =

inlet temperature (K)

Tint =

sampling period (s)

Uin =

inlet velocity (m/s)

U¯τ =

mean friction velocity (m/s)

txys =

smooth-vane tangential unit vector (–)

Tu =

turbulence intensity (m/s)

Greek Symbols

γ =

heat capacity ratio

δe =

boundary layer thickness (m)

δij =

Kronecker delta (–)

δ* =

displacement thickness (m)

ΔT =

temperature difference (K)

θ =

momentum thickness (m)

κ =

von Kármán coefficient

μ =

dynamic viscosity (kg/m/s)

ν =

kinematic viscosity (m2/s)

ρ =

density (kg/m3)

τij =

viscous stress tensor (kg/m/s2)

τw =

wall shear stress (kg/m/s2)

Dimensionless Groups

M =

inlet Mach number, MUin/γRTin

H =

shape factor, δ*/θ

K =

Kurtosis

S =

skewness

Cf =

instantaneous skin friction coefficient, 2τw/(ρinUin2)

C¯f =

mean skin friction coefficient, 2τw/(ρ¯eU¯e2)

C¯h =

mean Stanton number, Nu¯/(RePr)

Reτ =

friction Reynolds number, U¯τδe/νw

Reθ =

momentum thickness Reynolds number, U¯eθ/νw

Nu =

instantaneous Nusselt number, (c/ΔT)/T(x,z,t)n|w

Pr =

Prandtl number, cpμin/k

Re =

inlet Reynolds number, ρinUinc/μin

Nu¯ =

mean Nusselt number, (c/ΔT)/T¯n|w

Superscripts and Subscripts

+ =

viscous-scaled quantity

0 =

reference quantity

e =

boundary layer edge quantity

in =

inlet quantity

r =

rough-vane quantity

s =

smooth-vane quantity

t =

tangential component

tot =

total quantity

w =

wall quantity

Averaging Notations

()¯ =

time-then-spanwise average

() =

time-then-surface-area average

References

1.
Bons
,
J. P
,
2002
, “
St and Cf Augmentation for Real Turbine Roughness With Elevated Freestream Turbulence
,”
ASME J. Turbomach.
,
124
(
4
), pp.
632
644
.
2.
Bons
,
J. P.
,
2010
, “
A Review of Surface Roughness Effects in Gas Turbines
,”
ASME J. Turbomach.
,
132
(
2
), p.
021004
.
3.
Bammert
,
K.
, and
Sandstede
,
H.
,
1980
, “
Measurements of the Boundary Layer Development Along a Turbine Blade With Rough Surfaces
,”
ASME J. Eng. Power
,
102
(
4
), pp.
978
983
.
4.
Stripf
,
M.
,
Schulz
,
A.
, and
Wittig
,
S
,
2005
, “
Surface Roughness Effects on External Heat Transfer of a HP Turbine Vane
,”
ASME J. Turbomach.
,
127
(
1
), pp.
200
208
.
5.
Bai
,
T.
,
Liu
,
J.
,
Zhang
,
W.
, and
Zou
,
Z
,
2014
, “
Effect of Surface Roughness on the Aerodynamic Performance of Turbine Blade Cascade
,”
Propuls. Power Res.
,
3
(
2
), pp.
82
89
.
6.
Abuaf
,
N.
,
Bunker
,
R. S.
, and
Lee
,
C. P
,
1997
, “
Effects of Surface Roughness on Heat Transfer and Aerodynamic Performance of Turbine Airfoils
,”
ASME J. Turbomach.
,
102
(
3
), pp.
522
529
.
7.
Bacci
,
T.
,
Picchi
,
A.
,
Lenzi
,
T.
,
Facchini
,
B.
, and
Innocenti
,
L.
,
2021
, “
Effect of Surface Roughness and Inlet Turbulence Intensity on a Turbine Nozzle Guide Vane External Heat Transfer: Experimental Investigation on a Literature Test Case
,”
ASME J. Turbomach.
,
143
(
4
), p.
041006
.
8.
Thole
,
K. A.
,
Radomsky
,
R. W.
,
Kang
,
M. B.
, and
Kohli
,
A
,
2002
, “
Elevated Freestream Turbulence Effects on Heat Transfer for a Gas Turbine Vane
,”
Int. J. Heat Fl. Flow
,
23
(
2
), pp.
137
147
.
9.
Sandberg
,
R. D.
, and
Michelassi
,
V.
,
2022
, “
Fluid Dynamics of Axial Turbomachinery: Blade- and Stage-Level Simulations and Models
,”
Annu. Rev. Fluid Mech.
,
54
(
1
), pp.
255
285
.
10.
Nardini
,
M.
,
Jelly
,
T. O.
,
Kozul
,
M.
, and
Sandberg
,
R. D
,
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
.
11.
Nardini
,
M.
,
Jelly
,
T. O.
,
Kozul
,
M.
,
Sandberg
,
R. D.
,
Vitt
,
P.
, and
Sluyter
,
R. 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 Vane
,”
ASME J. Turbomach.
,
146
(
3
), p.
031009
.
12.
Jelly
,
T. O.
,
Nardini
,
M.
,
Rosenzweig
,
M.
,
Leggett
,
J.
,
Marusic
,
I.
, and
Sandberg
,
R. D.
,
2023
, “
High-Fidelity Computational Study of Roughness Effects on High Pressure Turbine Performance and Heat Transfer
,”
Int. J. Heat Fl. Flow
,
101
, p.
109134
.
13.
Arts
,
T.
, and
Lambert de Rouvroit
,
M.
,
1992
, “
Aero-Thermal Performance of a Two-Dimensional Highly Loaded Transonic Turbine Nozzle Guide Vane: A Test Case for Inviscid and Viscous Flow Computations
,”
VKI Training Center for Experimental Aerodynamics
,
Technical Report 174
, von Karman Institute for Fluid Dynamics, Belgium.
14.
Jelly
,
T. O.
, and
Busse
,
A.
,
2019
, “
Multi-Scale Anisotropic Rough Surface Algorithm: Technical Documentation and User Guide
,”
Technical Report
,
University of Glasgow
, Glasgow, UK.
15.
Taylor
,
R. P.
,
1990
, “
Surface Roughness Measurements on Gas Turbine Blades
,”
ASME J. Turbomach.
,
112
(
2
), pp.
175
180
.
16.
Bogard
,
D. G.
,
Schmidt
,
D. L.
, and
Tabbita
,
M
,
1996
, “
Characterization and Laboratory Simulation of Turbine Airfoil Surface Roughness and Associated Heat Transfer
,”
ASME J. Turbomach.
,
120
(
2
), pp.
337
342
.
17.
Bons
,
J. P.
,
Taylor
,
R. P.
,
McClain
,
S. T.
, and
Rivir
,
R. B
,
2001
, “
The Many Faces of Turbine Surface Roughness
,”
ASME J. Turbomach.
,
123
(
3
), pp.
739
748
.
18.
Busse
,
A.
, and
Jelly
,
T. O.
,
2023
, “
Effect of High Skewness and Kurtosis on Turbulent Channel Flow Over Irregular Rough Walls
,”
J. Turbul.
,
24
(
1–2
), p.
e2173761
.
19.
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
.
20.
Deuse
,
M.
, and
Sandberg
,
R. D.
,
2020
, “
Implementation of a Stable High-Order Overset Grid Method for High-Fidelity Simulations
,”
Comp. Fluids
,
211
, p.
104449
.
21.
Schlanderer
,
S. C.
,
Weymouth
,
G. D.
, and
Sandberg
,
R. D
,
2017
, “
The Boundary Data Immersion Method for Compressible Flows With Application to Aeroacoustics
,”
J. Comp. Phys.
,
333
, pp.
440
461
.
22.
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. Comp. Phys.
,
186
(
2
), pp.
652
665
.
23.
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
.
24.
Moody
,
L. F.
,
1944
, “
Friction Factors for Pipe Flow
,”
Trans. ASME
,
66
, pp.
671
684
.
25.
Busse
,
A.
,
Thakkar
,
M.
, and
Sandham
,
N. D.
,
2017
, “
Reynolds-Number Dependence of the Near-Wall Flow Over Irregular Rough Surfaces
,”
J. Fluid Mech.
,
810
, pp.
196
224
.
26.
Flack
,
K. A.
,
Schultz
,
M. P.
, and
Connelly
,
J. S.
,
2007
, “
Examination of a Critical Roughness Height for Outer Layer Similarity
,”
Phys. Fluids
,
19
(
9
), p.
095104
.
27.
Nikuradse
,
J.
,
1933
,
Laws of Flow in Rough Pipes
, VDI Forschungsheft, Vol.
361
, von Karman Institute for Fluid Dynamics, Belgium.
28.
Zhao
,
Y.
, and
Sandberg
,
R. D.
,
2021
, “
High-Fidelity Simulations of a High-Pressure Turbine Vane Subject to Large Disturbances: Effect of Exit Mach Number on Losses
,”
ASME J. Turbomach.
,
143
(
9
), p.
091002
.
29.
Chung
,
D. C.
,
Hutchins
,
N.
,
Schultz
,
M. P.
, and
Flack
,
K. A.
,
2021
, “
Predicting the Drag of Rough Surfaces
,”
Annu. Rev. Fluid Mech
,
53
, pp.
439
471
.
30.
Flack
,
K. A.
, and
Schultz
,
M. P.
,
2010
, “
Review of Hydraulic Roughness Scales in the Fully Rough Regime
,”
ASME J. Fluids Eng.
,
132
(
4
), p.
041203
.
31.
Scaggs
,
W. F.
,
Taylor
,
R. P.
, and
Coleman
,
H. W.
,
1988
, “
Measurement and Prediction of Rough Wall Effects on Friction Factor—Uniform Roughness Results
,”
ASME J. Fluids Eng.
,
110
(
4
), pp.
385
391
.
32.
Jelly
,
T. O.
,
Abu Rowin
,
W.
,
Hutchins
,
N.
,
Chung
,
D.
,
Tanimoto
,
K.
,
Oda
,
T.
, and
Sandberg
,
R. D.
,
2023
, “
High-Fidelity Computational Assessment of Aero-Thermal Performance and the Reynolds’ Analogy for Additively Manufactured Anisotropic Surface Roughness
,”
ASME J. Turbomach.
,
145
(
11
), p.
111005
.
33.
Kožulović
,
D.
,
Juchmann
,
N.
,
Führing
,
A.
, and
Bode
,
C.
,
2022
, “
Transition Model Extension for Roughness Effects
,” T
urbo Expo: Power for Land, Sea, and Air
,
Rotterdam, The Netherlands
,
June 13–17
,
American Society of Mechanical Engineers
, p.
V10CT32A042
.
34.
Hunt
,
J. C. R.
,
Wray
,
A. A.
, and
Moin
,
P.
,
1988
, “
Eddies, Streams, and Convergence Zones in Turbulent Flows
,”
Technical Report CTR-S88
,
Center for Turbulence Research Report
, von Karman Institute for Fluid Dynamics, Belgium.
35.
Robinson
,
S. K.
,
1991
, “
Coherent Motions in the Turbulent Boundary Layer
,”
Annu. Rev. Fluid Mech.
,
23
(
1
), pp.
601
639
.
36.
Pichler
,
R.
,
Sandberg
,
R. D.
,
Laskowski
,
G.
, and
Michelassi
,
V.
,
2017
, “
High-Fidelity Simulations of a Linear HPT Vane Cascade Subject to Varying Inlet Turbulence
,”
Turbo Expo: Power for Land, Sea, and Air, Vol. 50787
,
Chartlotte, NC
,
June 26–30
, p.
V02AT40A001
.
37.
Bunker
,
R. S.
,
2013
, “
Gas Turbine Cooling: Moving From Macro to Micro Cooling
,”
Turbo Expo: Power for Land, Sea, and Air, Volume 55164
,
San Antonio, TX
,
June 3–7
, p.
V03CT14A002
.
38.
Rouhi
,
A.
,
Endrikat
,
S.
,
Modesti
,
D.
,
Sandberg
,
R. D.
,
Oda
,
T.
,
Tanimoto
,
K.
,
Hutchins
,
N.
, and
Chung
,
D.
,
2022
, “
Riblet-Generated Flow Mechanisms that Lead to Local Breaking of Reynolds Analogy
,”
J. Fluid Mech.
,
951
, p.
A45
.
39.
Schlatter
,
P.
, and
Örlü
,
R.
,
2010
, “
Assessment of Direct Numerical Simulation Data of Turbulent Boundary Layers
,”
J. Fluid Mech.
,
659
, pp.
116
126
.
40.
Bons
,
J. P.
,
2005
, “
A Critical Assessment of Reynolds Analogy for Turbine Flows
,”
ASME J. Heat Transfer
,
127
(
5
), pp.
472
485
.
41.
Belnap
,
B. J.
,
Van Rij
,
J. A.
, and
Ligrani
,
P. M.
,
2002
, “
A Reynolds Analogy for Real Component Surface Roughness
,”
Int. J. Heat Mass Transf.
,
45
(
15
), pp.
3089
3099
.
42.
Fernholz
,
H. H.
, and
Finleyt
,
P. J.
,
1996
, “
The Incompressible Zero-Pressure-Gradient Turbulent Boundary Layer: An Assessment of the Data
,”
Prog. Aerosp. Sci.
,
32
(
4
), pp.
245
311
.
43.
Denton
,
J. D.
,
1993
, “
Loss Mechanisms in Turbomachines
,”
ASME J. Turbomach.
,
115
(
4
), pp.
621
656
.