4.1. Systems of the equations and the spectral method

In this study, linear stability analysis is performed for both incompressible and compressible flows to evaluate the compressibility effects. The procedure for incompressible flows was documented in the aforementioned study (Suzuki et al. Reference Suzuki, Shur, Strelets and Travin2022); accordingly, the procedure for compressible flows is mainly formulated in the rest of this section. For both problems, the approach taken by Takeuchi & Jankowski (Reference Takeuchi and Jankowski1981) is essentially followed: the Navier–Stokes equations are linearized by assuming axisymmetry of the mean flow as well as homogeneity in the axial direction, and eigenvalue problems are solved using a spectral method with the Chebyshev polynomials (Orszag Reference Orszag1971). The axial and azimuthal mean-velocity profiles, expressed as $U(r)$ and $V(r)$ respectively, are taken from the aforementioned IDDES (Suzuki et al. Reference Suzuki, Shur, Strelets and Travin2022) approximately at the middle of the interstage, $x/D \approx 0.18$ of the SDT geometry. Here, $x = 0$ is taken at the fan-blade pitch axis, as shown in figure 3. These velocity profiles at three engine speeds are compared in figure 4, highlighting that the azimuthal-velocity peak is most pronounced at the approach condition for the SDT. For compressible flows, the mean density and temperature profiles, expressed as $\rho (r)$ and $T(r)$ respectively, are additionally taken from the IDDES. Since all these profiles including the thermodynamic quantities are assumed to be invariant in the axial direction, the axial pressure gradient is accordingly neglected; namely, the viscous effect generating the axial pressure gradient is assumed to be much smaller than the radial variation and other terms.

Figure 4.

Mean-velocity profiles in the interstage of SDT at three engine speeds (normalized by the ambient speed of sound) at $x/D \approx 0.18$ computed by the IDDES. (a) Azimuthal velocity. (b) Axial velocity.

4.2. Temporal analysis and neutrally stable boundaries

The compressible Navier–Stokes equations are non-dimensionalized by using the ambient density and speed of sound for the air, expressed as $\rho _{\infty }$ and $a_{\infty } = \sqrt {\gamma R T_{\infty }}$ ($= 1.226\,{\rm kg}\,{\rm m}^{-3}$ and $340.2\,{\rm m}\,{\rm s}^{-1}$ in this study) respectively, as well as the fan diameter $D$ ($= 0.5588$ m). By expressing small flow perturbations as $\hat {f} ( r ) \exp [ - {\rm i} ( \omega t - k_x x - m \phi ) ]$, where $\omega$, $k_x$ and $m$ respectively denote the angular frequency, the axial wavenumber and the azimuthal-mode number, a system of the non-dimensional linearized equations in the Fourier domain is written in the cylindrical coordinates as (the corresponding system for an incompressible flow can be found in Suzuki et al. Reference Suzuki, Shur, Strelets and Travin2022)

(4.1a)\begin{align} &{\rm i} k_x \hat{u} + {\rm i} \frac{m}{r} \hat{v} + \left( \frac{\partial}{\partial r} + \frac{1}{r} + \frac{1}{\rho} \frac{\partial \rho}{\partial r} \right) \hat{w} + {\rm i} \frac{1}{\rho} \left( \frac{m V}{r} + k_x U \right) \hat{\rho} = {\rm i} \omega \frac{1}{\rho} {\hat \rho},\\ \end{align}

(4.1b)\begin{align} &\left[ {\rm i} \left( \frac{m V}{r} + k_x U \right) - \frac{\mu Ma}{\rho Re} \left( \Delta - \frac{k_x^2}{3} + \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \frac{\partial}{\partial r} \right) \right] \hat{u} \nonumber\\ &\qquad + \frac{k_x m \mu Ma}{3 r \rho Re} \hat{v} + \left[ \frac{\partial U}{\partial r} - {\rm i} \frac{k_x \mu Ma}{\rho Re} \left( \frac{1}{3} \frac{\partial }{\partial r} + \frac{1}{3r} + \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \right) \right] \hat{w} \nonumber\\ &\qquad + {\rm i} \frac{k_x T}{\gamma \rho} \hat{\rho} + \left[ {\rm i} \frac{k_x}{\gamma} - \frac{\mu Ma}{\rho Re} \frac{{\rm d} \mu}{{\rm d} T} \left( \frac{\partial U}{\partial r} \frac{\partial}{\partial r} + \frac{\partial^2 U}{\partial r^2} + \frac{1}{r} \frac{\partial U}{\partial r} \right) \right] \hat{T} = {\rm i} \omega \hat{u},\\ \end{align}

(4.1c)\begin{align} & \frac{k_x m \mu Ma}{3 r \rho Re} \hat{u} \nonumber\\ &\qquad + \left[ {\rm i} \left(\frac{m V}{r} + k_x U\right) - \frac{\mu Ma}{\rho Re} \left\{ \Delta - \frac{1}{r^2} - \frac{m^2}{3 r^2} + \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \left( \frac{\partial}{\partial r} - \frac{1}{r} \right) \right\} \right] \hat{v} \nonumber\\ &\qquad + \left[ \left( \frac{\partial V}{\partial r} + \frac{V}{r} \right) - {\rm i} \frac{m \mu Ma}{r \rho Re} \left( \frac{7}{3r} + \frac{1}{3} \frac{\partial}{\partial r} + \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \right) \right] \hat{w} + {\rm i} \frac{m T}{\gamma \rho r} {\hat \rho} \nonumber\\ &\qquad + \left[ {\rm i} \frac{m}{\gamma r} - \frac{\mu Ma}{\rho Re} \frac{{\rm d} \mu}{{\rm d} T} \left\{ \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right) \frac{\partial}{\partial r} + \left( \frac{\partial^2 V}{\partial r^2} + \frac{1}{r} \frac{\partial V}{\partial r} \right) \right\} \right] \hat{T} = {\rm i} \omega \hat{v},\\ \end{align}

(4.1d)\begin{align} &- {\rm i} \frac{k_x \mu Ma}{3 \rho Re} \left( \frac{\partial}{\partial r} - 2 \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \right) \hat{u} \nonumber\\ &\qquad - \left[ \frac{2~V}{r} - {\rm i} \frac{m \mu Ma}{3 r \rho Re} \left( \frac{7}{r} - \frac{\partial }{\partial r} + 2 \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \right) \right] \hat{v} \nonumber\\ &\qquad + \left[ {\rm i} \left( \frac{m V}{r} + k_x U \right) \right.\nonumber\\ &\left.\qquad - \frac{\mu Ma}{\rho Re} \left\{ \Delta + \frac{1}{3} \frac{\partial^2}{\partial r^2} + \frac{1}{3r} \frac{\partial}{\partial r} - \frac{4}{3 r^2} + \frac{2}{3} \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \left( 2 \frac{\partial}{\partial r} - \frac{1}{r} \right) \right\} \right] \hat{w} \nonumber\\ &\qquad + \frac{1}{\gamma \rho} \left( T \frac{\partial }{\partial r} + \frac{\partial T }{\partial r} \right) \hat{\rho} \nonumber\\ &\qquad + \left[ \frac{1}{\gamma} \left( \frac{\partial }{\partial r} + \frac{1}{\rho} \frac{\partial \rho}{\partial r} \right) - \frac{\mu Ma}{\rho Re} \frac{{\rm d} \mu}{{\rm d} T} \left\{ {\rm i} \frac{m}{r} \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right) + {\rm i} k_x \frac{\partial U}{\partial r} \right\} \right] \hat{T} = {\rm i} \omega \hat{w},\\ \end{align}

(4.1e)\begin{align} & -\frac{2 \gamma \mu Ma}{\rho Re} \frac{\partial U}{\partial r} \frac{\partial}{\partial r} \hat{u} - \frac{2 \gamma \mu Ma}{\rho Re} \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right) \left( \frac{\partial }{\partial r} - \frac{1}{r} \right) \hat{v} \nonumber\\ &\qquad + \left[ \frac{1}{\gamma - 1} \frac{\partial T}{\partial r} - \frac{T}{\rho} \frac{\partial \rho}{\partial r} - {\rm i} \frac{2 \gamma \mu Ma}{\rho Re} \left\{ k_x \frac{\partial U}{\partial r} + \frac{m}{r} \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right) \right\} \right] \hat{w} \nonumber\\ &\qquad - {\rm i} \frac{T}{\rho} \left( \frac{m V}{r} + k_x U \right) \hat{\rho} \nonumber\\ &\qquad + \frac{{\rm i}}{\gamma - 1} \left[ \left( \frac{m V}{r} + k_x U \right) + \frac{{\rm i} \gamma \mu Ma}{\rho Re Pr} \left\{ \Delta +\frac{{\rm d} \mu}{{\rm d} T} \left( \frac{\partial^2~T}{\partial r^2} + \frac{1}{r} \frac{\partial T}{\partial r} + 2 \frac{\partial T}{\partial r} \frac{\partial}{\partial r} \right) \right.\right.\nonumber \\ &\left.\left.\qquad + \left( \gamma - 1 \right) Pr \frac{{\rm d} \mu}{{\rm d} T} \left( \left( \frac{\partial U}{\partial r} \right)^2 + \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right)^2 \right) \right\} \right] \hat{T} = {\rm i} \omega \left[ \frac{\hat{T}}{\gamma - 1} - \frac{T}{\rho} \hat{\rho} \right], \end{align}

where $\hat {u}$, $\hat {v}$ and $\hat {w}$ denote the axial-, azimuthal- and radial-velocity components of the perturbation, respectively, and $\hat {\rho }$ and $\hat {T}$ the density and temperature perturbations, respectively. The Reynolds number, the Prandtl number and the Mach number are defined as $Re \equiv \rho _{\infty } \bar {U} D / \bar {\mu }$, $Pr \equiv \gamma R \bar {\mu } / (\gamma - 1 ) \bar {\kappa } = 0.72$ and $Ma \equiv \bar {U} / a_{\infty }$, where the gas constant $R$ is assumed to be constant; moreover, the quantities with the top bar are those averaged over the interstage station (refer to table 1 for $\bar {U}$ and $\bar {T}$ in each condition) only to characterize the Reynolds number. To include the temperature dependence of the dynamic viscosity $\mu$ and the heat conductivity $\kappa$, the Sutherland formula (Sutherland Reference Sutherland1893) is adopted (the dimensional formula is given below)

(4.2)\begin{equation} \mu = \frac{ 0.72 \left( \gamma - 1 \right)}{\gamma R} \kappa \equiv \mu_0 \left( \frac{T}{T_0} \right)^{1.5} \frac{T_0 + C}{T + C} , \end{equation}

where $\mu _0 = 1.716 \times 10^{-5}$ kg (s m)$^{-1}$, $T_0 = 273.15$ K, $C = 110.5$ K, respectively. The same constants are used for the heat conductivity, as it is related by the Prandtl number defined above. The Laplacian in (4.1) is also defined as

(4.3)\begin{equation} \Delta \equiv \frac{\partial^2}{\partial r^2} + \frac{1}{r} \frac{\partial }{\partial r} - \frac{m^2}{r^2} - k_x^2 , \end{equation}

which is different from the notation by Takeuchi & Jankowski (Reference Takeuchi and Jankowski1981). On both sides of the walls, the no-slip boundary conditions are imposed for the three velocity components.

The equations above are cast such that the angular frequency $\omega$ is linearly placed on the right-hand side and all the other terms on the left-hand side. Accordingly, temporal instability can be solved as a generalized eigenvalue problem using a spectral method developed by Orszag (Reference Orszag1971): using the Chebyshev polynomials, the system of the equation (4.1) is discretized in the radial direction into $N_p$ ($= 256$ in most cases) cells in a staggered manner (Macaraeg, Streett & Hussaini Reference Macaraeg, Streett and Hussaini1988). Namely, the three velocity components are located at $N_p+1$ grid points including the boundary points on the wall, while the density and the temperature are assigned at $N_p$ mid points; accordingly, $3 \times ( N_p+1 )$ equations for the momentum and $2 \times N_p$ equations for continuity and energy are solved in essence. Yet, six equations are sacrificed in the highest orders of these equations for the no-slip boundary conditions on both walls. Consequently, the total of $5N_p + 3$ equations are simultaneously solved for a compressible flow, as opposed to $4N_p + 3$ equations for an incompressible flow. It should be noted that the incompressible solver has been verified by comparing with the results reported by Takeuchi & Jankowski (Reference Takeuchi and Jankowski1981).

We begin with comparing the neutrally stable boundaries (the lines at which the most unstable growth rate vanishes) for incompressible and compressible solutions by solving this temporal problem (as the boundaries are common to a spatial problem). Figure 5 plots those for $m = 22$ (which is the same as the fan-blade count as an example) at the approach condition up to the Reynolds number at which the current spectral code can resolve the eigenvalue problem (the Reynolds number of the SDT/IDDES is $Re \approx 4 \times 10^6$). The solutions for $m = -22$ become the mirror image of these figures (i.e. the signs of $\omega$ and $k_x$ are flipped). First, the boundaries for incompressible and compressible flows almost overlap each other, indicating negligible compressibility effects. Second, we find three distinct unstable regimes, as labelled in figure 5(b): the first regime becomes unstable over a relatively wide frequency range, and its boundary spreads with the Reynolds number. This mode can destabilize disturbances moving with the rotor speed, indicated by the dashed vertical line in figure 5(a), i.e. the rotating speed of the wakes given by $\omega = m \varOmega$. The second regime resides only within a narrow range of small axial wavenumbers at low frequencies, corresponding to the local peak azimuthal velocity or lower near the outer wall, i.e. $\omega \approx mV/r$, indicated by the dotted line. These two modes rotate in the direction of the rotor. The third regime becomes unstable at relatively higher Reynolds numbers but spreads its boundary very rapidly. This mode physically appears as a counter-rotating mode (i.e. $m = -22$ in this case) with the opposite sign for $\omega$ and $k_x$; hence, it can still convect downstream. As discussed later, the rotational direction of a disturbance is indistinguishable from snapshots of flow fields in post-processing unless we track the evolution in time.

Figure 5.

Neutrally stable boundaries in the interstage of SDT for $m = 22$ at $61.7\,\%$ RPM. (a) Frequency dependence. (b) Wavenumber dependence.

Figure 6 similarly compares the neutrally stable boundaries at the cut-back condition. It is clear that the compressibility is still negligible. We find similar three unstable regimes; to be precise, the second and the third regimes merge at higher Reynolds numbers, and the boundary of the first regime depicts a gentle double-hump shape. More importantly, the first unstable regime spreads over wider frequencies but migrates to a higher-frequency range beyond the rotor speed; in addition, the critical Reynolds number of the first regime becomes noticeably higher than that at the approach. This implies that disturbances in the wakes are less amplified compared with the lower speed. On the other hand, the second unstable regime hardly shifts the frequency range relative to the local azimuthal velocity, but its wavenumber range migrates to the negative side in contrast to the approach condition.

Figure 6.

Neutrally stable boundaries in the interstage of SDT for $m = 22$ at $87.5\,\%$ RPM. (a) Frequency dependence. (b) Wavenumber dependence.

When the engine speed is increased to the take-off speed, the first and third regimes cannot be found up to the Reynolds-number range which this code can handle ($Re \lesssim 5 \times 10^5$), even though unstable modes are searched over a wider range than the two previous lower engine speeds (i.e. $-600 \lesssim k_x \lesssim 500$). Presumably, the critical Reynolds numbers of these two regimes are even higher than the resolvable range, possibly because the peak of the azimuthal velocity is smeared at the take-off speed, as shown in figure 4. Figure 7 shows that the second unstable regime stays almost in the same $\omega$ and $k_x$ ranges as those at the cut-back condition. We also notice that the boundaries for incompressible and compressible flows still agree very well even at the full take-off speed (the tip Mach number of $M_{{ tip}} \approx 1.2$). Therefore, it is safe to conclude that the compressibility effect itself on the hydrodynamic instability is negligible in the interstage, while the velocity profiles, particularly in the azimuthal component, change with the engine speed, and the profile change has a greater impact on the instability characteristics.

Figure 7.

Neutrally stable boundaries in the interstage of SDT for $m = 22$ at $100\,\%$ RPM. (a) Frequency dependence. (b) Wavenumber dependence.

Figure 8 also plots the neutrally stable boundaries for $m =39$ at the approach condition. This azimuthal-mode number is close to the upper threshold of the stripes associated with the fan broadband noise in this condition (refer to the discussion in Suzuki et al. Reference Suzuki, Shur, Strelets and Travin2022). While the critical Reynolds numbers for $m =39$ are higher in all the regimes, the three unstable regimes depict similar patterns to those for $m =22$ shown in figure 5 (interestingly, the first unstable regime depicts a clear double-hump shape even though processing the same velocity profiles as $m =22$). Thus, we still anticipate that the flow is unstable over the range of azimuthal modes mainly contained in the wake turbulence.

Figure 8.

Neutrally stable boundaries in the interstage of SDT for $m = 39$ at $61.7\,\%$ RPM. (a) Frequency dependence. (b) Wavenumber dependence.

We should remark that the second unstable regime stays close to the local peak velocity regardless of the engine speed or the azimuthal mode (refer to the $\omega - Re$ plots from figures 5–8). This suggests that the pointy velocity profile near the outer wall drives the second unstable mode. On the other hand, the first regime becomes highly unstable at the rotor speed for both $m = 22$ and $39$, but only at the approach condition. Referring to the wave equation in the cylindrical coordinates (i.e. the Bessel equation), it may be possible to scale $\omega$ and $k_x$ proportionally with $m$ as a limit of a narrow annulus. Namely, the ranges of the unstable frequency and wavenumber shift approximately proportionally; as a result, these unstable ranges stay about the same relative to the rotor speed and the local azimuthal velocity at the approach condition. As the engine speed increases, the dispersion relation of $\mbox {Re}[ \omega ] - \mbox {Re}[ k_x ]$ retains the same similarity, but the unstable regime of the first instability shifts to higher frequencies relative to the rotor speed. Thus, the first unstable mode amplifies the wake disturbances less at higher engine speeds. We will closely observe the azimuthal-mode dependence in § 4.4.

4.3. Spatial analysis of instability waves and acoustic waves

Next, a spatial problem is solved to observe the dispersion relations and the eigenfunctions representative of a turbomachinery problem. Here, (4.1) is recast so that the axial wavenumber $k_x$ is linearly placed on the right-hand side as

(4.4a)\begin{align} & {\rm i} \frac{m}{r} \hat{v} + \left( \frac{\partial }{\partial r} + \frac{1}{r} + \frac{1}{\rho} \frac{\partial \rho}{\partial r} \right) \hat{w} + {\rm i} \frac{1}{\rho} \left( \frac{m V}{r} - \omega \right) \hat{\rho} ={-} {\rm i} k_x \left[ \hat{u} + \frac{U}{\rho} \hat{\rho} \right], \\\end{align}

(4.4b)\begin{align} & \left[ {\rm i} \left( \frac{m V}{r} - \omega \right) - \frac{\mu Ma}{\rho Re} \left( \Delta - \frac{k_x^2}{3} + \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \frac{\partial}{\partial r} \right) \right] \hat{u} \nonumber\\ &\qquad + \frac{\partial U}{\partial r} \hat{w} - \frac{\mu Ma}{\rho Re} \frac{{\rm d} \mu}{{\rm d} T} \left( \frac{\partial U}{\partial r} \frac{\partial}{\partial r} + \frac{\partial^2~U}{\partial r^2} + \frac{1}{r} \frac{\partial U}{\partial r} \right) \hat{T} ={-} {\rm i} k_x \nonumber\\ & \left[U \hat{u} - {\rm i} \frac{m \mu Ma}{3 r \rho Re} \hat{v} - \frac{\mu Ma}{\rho Re} \left( \frac{1}{3} \frac{\partial }{\partial r} + \frac{1}{3r} + \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \right) \hat{w} + \frac{T}{\gamma \rho} \hat{\rho} + \frac{1}{\gamma} \hat{T} \right], \\\end{align}

(4.4c)\begin{align} & \left[ {\rm i} \left( \frac{m V}{r} - \omega \right) - \frac{\mu Ma}{\rho Re} \left\{ \Delta - \frac{1}{r^2} - \frac{m^2}{3 r^2} + \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \left( \frac{\partial}{\partial r} - \frac{1}{r} \right) \right\} \right] \hat{v} \nonumber\\ &\qquad+ \left[ \left( \frac{\partial V}{\partial r} + \frac{V}{r} \right) - \frac{{\rm i} m \mu Ma}{r \rho Re} \left( \frac{7}{3r} + \frac{1}{3} \frac{\partial}{\partial r} + \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \right) \right] \hat{w} + {\rm i} \frac{m T}{\gamma r \rho} {\hat \rho} \nonumber\\ &\qquad+ \left[ {\rm i} \frac{m}{\gamma r} - \frac{\mu Ma}{\rho Re} \frac{{\rm d} \mu}{{\rm d} T} \left\{ \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right) \frac{\partial}{\partial r} + \left( \frac{\partial^2~V}{\partial r^2} + \frac{1}{r} \frac{\partial V}{\partial r} \right) \right\} \right] \hat{T} ={-} {\rm i} k_x \nonumber\\ &\qquad\times \left[ -{\rm i} \frac{m \mu Ma}{3 r \rho Re} \hat{u} + U \hat{v} \right], \\\end{align}

(4.4d)\begin{align} & - \left[ \frac{2V}{r} - {\rm i} \frac{m \mu Ma}{3 r \rho Re} \left( \frac{7}{r} - \frac{\partial }{\partial r} + 2 \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \right) \right] \hat{v} + \left[ {\rm i} \left( \frac{m V}{r} - \omega \right) \right.\nonumber\\ &\left.\qquad- \frac{\mu Ma}{\rho Re} \left\{ \Delta + \frac{1}{3} \frac{\partial^2}{\partial r^2} + \frac{1}{3r} \frac{\partial}{\partial r} - \frac{4}{3 r^2} + \frac{2}{3} \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \left( 2 \frac{\partial}{\partial r} - \frac{1}{r} \right) \right\} \right] \hat{w} \nonumber\\ &\qquad+ \frac{1}{\gamma \rho} \left( T \frac{\partial }{\partial r} + \frac{\partial T }{\partial r} \right) \hat{\rho} \nonumber\\ &\qquad+ \left[ \frac{1}{\gamma} \left( \frac{\partial }{\partial r} + \frac{1}{\rho} \frac{\partial \rho}{\partial r} \right) - {\rm i} \frac{m \mu Ma}{r \rho Re} \frac{{\rm d} \mu}{{\rm d} T} \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right) \right] \hat{T} ={-} {\rm i} k_x\nonumber\\ &\qquad\times \left[ - \frac{\mu Ma}{3 \rho Re} \left( \frac{\partial}{\partial r} - 2 \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial T}{\partial r} \right) \hat{u} + U \hat{w} - \frac{\mu Ma}{\rho Re} \frac{{\rm d} \mu}{{\rm d} T} \frac{\partial U}{\partial r} \hat{T} \right], \\\end{align}

(4.4e)\begin{align} & -\frac{2 \gamma \mu Ma}{\rho Re} \frac{\partial U}{\partial r} \frac{\partial}{\partial r} \hat{u} - \frac{2 \gamma \mu Ma}{\rho Re} \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right) \left( \frac{\partial }{\partial r} - \frac{1}{r} \right) \hat{v} \nonumber\\ &\qquad+ \left[ \frac{1}{\gamma - 1} \frac{\partial T}{\partial r} - \frac{T}{\rho} \frac{\partial \rho}{\partial r} - {\rm i} \frac{2 \gamma m \mu Ma}{r \rho Re} \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right) \right] \hat{w} \nonumber\\ &\qquad- {\rm i} \frac{T}{\rho} \left( \frac{m V}{r} - \omega \right) \hat{\rho} + \frac{{\rm i}}{\gamma - 1} \left[ \left( \frac{m V}{r} - \omega \right) \right.\nonumber\\ &\qquad+ {\rm i} \frac{\gamma \mu Ma}{\rho Re Pr} \left\{ \Delta + \frac{{\rm d} \mu}{{\rm d}T} \left( \frac{\partial^2 T}{\partial r^2} + \frac{1}{r} \frac{\partial T}{\partial r} + 2 \frac{\partial T}{\partial r} \frac{\partial}{\partial r} \right) \right. \nonumber\\ &\left.\left.\qquad+ \left( \gamma - 1 \right) Pr \frac{{\rm d} \mu}{{\rm d}T} \left( \left( \frac{\partial U}{\partial r} \right)^2 + \left( \frac{\partial V}{\partial r} - \frac{V}{r} \right)^2 \right) \right\} \right] \hat{T} ={-} {\rm i} k_x \nonumber \\ &\qquad\times \left[ - \frac{2 \gamma \mu Ma}{\rho Re} \frac{\partial U}{\partial r} \hat{w} + \frac{U}{\gamma - 1} \hat{T} - \frac{T U}{\rho} \hat{\rho} \right]. \end{align}

Namely, it is expressed as ${\boldsymbol A} ( k_x^{n} ) {\boldsymbol u} = k_x^{n+1} {\boldsymbol B} {\boldsymbol u}$, where ${\boldsymbol u}$ represents a column vector consisting of the density, temperature and three velocity perturbations, and the superscript $n$ for $k_x$ denotes the iteration counter. A general eigenvalue problem can then be solved for given $\omega$ and $m$ by iteratively updating $k_x$ for a few times (refer to the numerical technique proposed by Bridges & Morris Reference Bridges and Morris1984). Thus, an eigenvalue $k_x$ for only one mode is accurately computed at a time for a spatial problem. The system of equations for an incompressible flow is formulated in (5) of Suzuki et al. (Reference Suzuki, Shur, Strelets and Travin2022).

We should remark that at greater $|\omega |$, solutions for the system of the compressible equations include acoustic modes. This section also reviews duct acoustic modes for the Euler equations and explores their potential interaction with unstable modes. By assuming the axial mean velocity to be uniform (the azimuthal mean velocity and the thermodynamic mean quantities are still functions of radius), Golubev & Atassi (Reference Golubev and Atassi1998) derived the wave equation in an infinitely extended co-axial cylinder as

(4.5)\begin{align} \frac{\partial^2 \hat{\varphi}}{\partial r^2} + \left( \frac{1}{r} + \frac{\partial \log {\rho}}{\partial r} \right) \frac{\partial \hat{\varphi}}{\partial r} + \left[ \left( \frac{\omega - k_x {U} - m {V}/r}{a} \right)^2 - k_x^2 - \frac{m^2}{r^2} \right] \hat{\varphi} ={-} \frac{1}{\rho} \boldsymbol{\nabla} \boldsymbol{\cdot} \left( \rho \hat{\boldsymbol u}_{{ rot}} \right) , \end{align}

where $\hat {\varphi }$ and $\hat {\boldsymbol u}_{{rot}}$ represent the velocity potential function and the rotational part of the velocity perturbation in the Fourier domain, respectively, and $a$ denotes the mean speed of sound. Although the density gradient is created in the radial direction due to the swirl velocity, the variation is mostly within $\pm$5 % in the interstage even at the cut-back condition; therefore, the mean density gradient in (4.5) is neglected below, and the mean speed of sound is assumed to be constant ($= a_{\infty }$) as well.

Since the rotational-velocity perturbation is related to the velocity potential as

(4.6)\begin{equation} -{\rm i} \left( \omega - k_x U - \frac{m V}{r} \right) \hat{\boldsymbol u}_{{ rot}} + \hat{\boldsymbol u}_{{ rot}} \boldsymbol{\cdot} \boldsymbol{\nabla} {\boldsymbol U} ={-} \left( \boldsymbol{\nabla} \times {\boldsymbol U} \right) \times \boldsymbol{\nabla} \varphi, \end{equation}

where ${\boldsymbol U}$ denotes the mean-velocity vector, it can be estimated as $\omega | \hat {\boldsymbol u}_{{rot}} | \sim | \partial V / \partial r | | \boldsymbol {\nabla } \varphi |$ by assuming $|\omega | \gtrsim | k_x U | + | m V/r | + | \partial V / \partial r |$. In the asymptote of $\omega \gg | \partial V / \partial r |$, (4.5) can be simplified as

(4.7)\begin{equation} \frac{\partial^2 \hat{\varphi}}{\partial r^2} + \frac{1}{r} \frac{\partial \hat{\varphi}}{\partial r} + \left[ \left( \frac{\omega - k_x {U} - m {V}/r}{a_{\infty}} \right)^2 - k_x^2 - \frac{m^2}{r^2} \right] \hat{\varphi} = 0 . \end{equation}

Golubev & Atassi (Reference Golubev and Atassi1996) derived this equation by assuming that the mean velocity has no vorticity, while Heaton & Peake (Reference Heaton and Peake2005) expanded the right-hand side by assuming rigid-body rotation. In reality, the azimuthal velocity in the interstage is in between them (i.e. the slope is rather flat). If the azimuthal velocity can be further characterized by rigid-body rotation, such as $V \approx \chi \varOmega r$, (4.7) can be regarded as the Bessel function. For a given frequency, a pair of real-value solutions can be found when

(4.8)\begin{equation} \left( \frac{\omega - m \chi \varOmega}{a_{\infty}} \right)^2 - \kappa^2 \left( 1- \frac{U^2}{a_{\infty}^2} \right) \ge 0 ,\quad \therefore \omega \ge \kappa a_{\infty} \sqrt{1- \frac{U^2}{a_{\infty}^2}} + m \chi \varOmega , \end{equation}

where $\kappa ^2 a_{\infty }^2 \equiv ( \omega - k_x U -m V/r )^2 - k_x^2 a_{\infty }^2$ scales the independent variable $r$ of the Bessel function. Thus, (4.8) crudely determines the cut-on boundaries of the duct acoustic modes in the interstage: the first term governs the cut-on duct acoustic modes in a uniform flow, and the second term implies the motion associated with the frame of reference spinning at $\chi \varOmega$. We will evaluate this relation in § 5.3 later.

Figure 9 plots the dispersion relations solved as a spatial problem (4.4) near unstable ranges for the same condition as figure 5 (i.e. $m = 22$ at the approach condition). Here, the figure compares three cases: the incompressible solutions at $Re = 10^5$ and $Re = 10^6$ and the compressible solution at $Re = 10^5$. The three solutions for the phase relation in figure 9(a) virtually overlap, indicating no impacts of the Reynolds number and the compressibility effect. The slope of the second unstable regime is slightly higher than the other two regimes. The growth rates in figure 9(b) somewhat depend on the Reynolds number; namely, the growth rates increase and the unstable ranges spread at the higher Reynolds number, particularly for counter-rotating modes (i.e. the third regime). The second mode always takes the highest growth rate, likely even at the Reynolds number of the real fan rig ($Re \approx 4 \times 10^6$).

Figure 9.

Dispersion relations as a spatial problem of SDT for $m = 22$ at $61.7\,\%$ RPM. (a) Axial wavenumber. (b) Growth rate.

For $m = 22$, acoustic modes start to branch out as a pair of highly decaying modes, one propagating downstream and the other upstream, somewhere at the non-dimensional angular frequencies of $\omega \gtrsim 60$ ($\approx 2.03 f_{{ BPF}}$) and $\omega \lesssim -40$ ($\approx -1.35 f_{{ BPF}}$) at the approach condition. These frequencies act as the threshold cut-on frequencies, which are skewed due to the swirl velocity, as explained in (4.8). As the frequency increases, the second pair of acoustic modes appears, and their decay rates become much less, say $\mbox {Im} [ k_x ] < 0.2$, at $\omega \gtrsim 72$ ($\approx 1.64 f_{{ BPF}}$) and $\omega \lesssim -49$ ($\approx -1.11 f_{{ BPF}}$). At these two frequencies, the eigenvalue distributions of the compressible analysis are plotted in figure 10 in the low $\mbox {Im} [ k_x ]$ range on the complex $k_x$ plane: as noted above, because the spatial problem is iteratively solved for a specific $k_x$ due to the nonlinear $k_x$ terms on the left-hand side of (4.4), values of $k_x$ are approximate except for the least-decaying acoustic mode in these plots.

Figure 10.

Eigenvalue distribution of $k_x$ at $61.7\,\%$ RPM ($Re = 10^5$); (a) $m = 22$ at $\omega = 72$, (b) $m = -22$ at $\omega = 49$.

In figure 10(a), there are no modes in the fourth quadrant (i.e. downstream-propagating growing mode) and in the second quadrant (i.e. upstream-propagating growing mode). In general, for co-rotating modes (i.e. $m > 0$), the hydrodynamic mode is stable or very weakly unstable at such high frequencies; therefore, we expect that no acoustic modes excite hydrodynamic modes significantly at the same frequency and azimuthal mode. For $\omega < 0$, figure 10(b) plots $k_x$ by setting $m = -22$ so that $\mbox {Re} [ k_x ]$ takes the opposite sign relative to figure 9. We now find an unstable mode in the fourth quadrant. Namely, for counter-rotating modes (i.e. $m < 0$), acoustic modes in both directions can excite the hydrodynamic mode in the third regime, but magnitudes of counter-rotating hydrodynamic disturbances may be insignificant in the interstage.

Figure 11(a) compares the dispersion relations at the cut-back condition. The phase relations for the three cases again overlap each other and retain the features found at the approach condition in figure 9(a). The comparison of the growth rates in figure 11(b) reveals that the Reynolds number has a significant impact, while the compressibility effect is still negligible at the cut-back condition. Compared with the approach condition, the first unstable regime spreads over a wider high-frequency range with a double-hump shape, but the peak growth rate at $Re = 10^5$ becomes lower; in contrast, the other two unstable regimes tend to remain in similar frequency ranges. These trends are consistent with those derived by the temporal analysis shown in figures 5 and 6.

Figure 11.

Dispersion relations as a spatial problem of SDT for $m = 22$ at $87.5\,\%$ RPM. (a) Axial wavenumber. (b) Growth rate.

Such a large high-frequency shift of the first unstable regime can keep the hydrodynamic mode unstable beyond the cut-on frequency of the acoustic mode mentioned above. At the cut-back condition, the acoustic modes start to brunch out at $\omega \gtrsim 65$ ($\approx 1.55 f_{{ BPF}}$) as a co-rotating mode of $m = 22$, and $\omega \lesssim -35$ ($\approx -0.84 f_{{ BPF}}$) as a counter-rotating mode of $m = -22$. Figure 12 similarly plots the complex $k_x$ distributions of both $m = \pm 22$ when the second acoustic modes appear on both directions. We confirm from figure 12(a) that the acoustic mode and the unstable hydrodynamic mode in the first regime now coexist with the same $m$ at the same $\omega$; however, this frequency ($\omega \gtrsim 77$) is much greater than the rotor speed ($m \varOmega \approx 41.9$). Therefore, it is unlikely that the amplification of hydrodynamic disturbances in the first unstable regime by acoustic modes has an appreciable impact on the rotor-related interaction noise.

Figure 12.

Eigenvalue distribution of $k_x$ at $87.5\,\%$ RPM ($Re = 10^5$); (a) $m = 22$ at $\omega = 77$, (b) $m = -22$ at $\omega = 43$.

At the take-off condition, we are unable to find unstable modes in the first and third regimes up to $Re = 10^6$ (the first regime may exist beyond $\omega > 300$, which is yet unresolvable with the current polynomial order). This is consistent with the observation in figure 7. This suggests that the interaction between the first and second unstable modes becomes weaker at higher engine speed. Accordingly, figure 13 plots the dispersion relations only in the second unstable regime. The growth rate is nearly doubled by increasing the Reynolds number by one order of magnitude. Across the three engine speeds, the peak growth rates of the second mode are found to be comparable at $Re = 10^6$, indicating the existence of the second unstable mode regardless of the engine speed (more precisely, the velocity profiles of these engine speeds). On the other hand, the compressibility effect remains to be negligible even at the take-off condition.

Figure 13.

Dispersion relations as a spatial problem of SDT for $m = 22$ at $100\,\%$ RPM. (a) Axial wavenumber. (b) Growth rate.

We now observe the eigenfunctions in the radial direction based on the compressible analysis. Figure 14 plots the eigenfunctions of the first and second unstable modes near the most unstable frequency in each regime at the approach condition. Both modes largely fluctuate near the outer wall; to be precise, the first mode forms a peak approximately above the outer-wall boundary layer, while the second mode resides inside of it. In contrast, figure 15, which similarly compares the eigenfunctions at the cut-back condition, reveals that the first mode is peaked away from the outer wall, while the second mode is still confined within the outer-wall boundary layer. These features are consistent with the radial profiles of turbulent kinetic energy in the interstage compared later in figure 18(b). This may also imply that the two modes can readily couple each other at $61.7\,\%$ RPM, while they can hardly interact at $87.5\,\%$ RPM. If the interaction between the two modes contributes to the fan broadband noise, as hypothesized by Suzuki et al. (Reference Suzuki, Shur, Strelets and Travin2022), this noise-generation mechanism would be less important at high engine speeds. This may be one of the reasons why the contribution of the stripe patterns to the broadband noise is diminished in the azimuthal-mode map at higher speeds, as discussed in § 5.4 later, as well as from the mode-ring measurement by Premo & Joppa (Reference Premo and Joppa2002).

Figure 14.

Eigenfunctions of the two unstable regimes for $m = 22$ at $61.7\,\%$ RPM ($Re = 10^5$). (a) First mode ($\omega = 33$ in figure 9). (b) Second mode ($\omega = 9$). (i) Axial-velocity fluctuation. (ii) Azimuthal-velocity fluctuation. (iii) Pressure fluctuation.

Figure 15.

Eigenfunctions of the two unstable regimes for $m = 22$ at $87.5\,\%$ RPM ($Re = 10^5$). (a) First mode ($\omega = 119$ in figure 11). (b) Second mode ($\omega = 7$). (i) Axial-velocity fluctuation. (ii) Azimuthal-velocity fluctuation. (iii) Pressure fluctuation.

4.4. Azimuthal-mode dependence

We observe the azimuthal-mode dependence, which clarifies the unstable range relevant to the rotor wakes. Figure 16 compares the dispersion relations of three different azimuthal modes near unstable frequency ranges at the approach condition by solving the incompressible problems at $Re = 10^5$. As discussed in figure 8, $\omega$ and $k_x$ should reach the linear dependence on $m$ for a large asymptote of $m$ in (4.4); accordingly, figure 16 normalizes $\omega$ and $\mbox {Re} [ k_x ]$ with $m$ by taking the reference as the fan-blade count $B = 22$, except for $-\mbox {Im} [ k_x ]$. The phase relations collapse very well across all three regimes in figure 16(a). The unstable frequency ranges of the first and second regimes tend to align with the rotor speed and the local azimuthal velocity, respectively, in figure 16(b), but the peak growth rates depend on the azimuthal-mode number.

Figure 16.

Dispersion relations of different azimuthal modes at $61.7\,\%$ RPM ($Re = 10^5$). (a) Axial wavenumber. (b) Growth rate.

Figure 17 compares the peak growth rates of the first and second unstable regimes as a function of the azimuthal mode at $Re = 10^5$. At the approach condition, the first and second regimes take the maximum values at similar azimuthal-mode numbers; accordingly, the first mode tends to amplify the wake disturbances in the vicinity of the fan-blade count. Moreover, the peak growth rate vanishes for $m \gtrsim 40$, and also other types of modes seem to take a higher peak growth rate for $m \lesssim 6$ (not shown) in both regimes. At the cut-back condition, in contrast, the unstable range of the first mode substantially shifts to the lower $m$ side, while the that of the second mode extends on the higher $m$ side. In fact, the unstable ranges of the second regime at both engine speeds coincide with the bundle of the stripe pattern in the azimuthal-decomposition map, shown later in figures 24, 29 and 30 (we discuss this point in §§ 5.3 and 5.4).

Figure 17.

Peak growth rates as a function of the azimuthal mode at $Re = 10^5$; (a) $61.7\,\%$ RPM, (b) $87.5\,\%$ RPM.

5.1. Post-processing of the IDDES database

To analyse the mode structures in the interstage and minimize artifacts in post-processing of the IDDES data for a limited time duration, the computational time-step was set to be $\Delta t_{{ CFD}} \equiv 2 {\rm \pi}/ ( 2112 \varOmega )$, and an instantaneous velocity field was recorded in the rotating frame at every $16$ and $6$ time steps, i.e. $\Delta t_{{ DMD}} \equiv 2 {\rm \pi}/ ( 132 \varOmega )$ and $2 {\rm \pi}/ ( 352 \varOmega )$, at the approach and cut-back conditions, respectively. During two and one revolutions, a total of $N_t = 264$ and $352$ frames were recorded, respectively, while the azimuthal grid count in the interstage block was set to be $N_{\phi } = 792$. At the approach condition, the flow quantities in the rotating frame can then be projected onto the stationary frame by just clocking the mesh by six grids in the uniform grid region (in the upstream non-uniform grid region, they are interpolated into the same uniform grid), and also the same fan-blade position is repeated at every six recorded frames. At the cut-back condition, the flow quantities can be similarly projected on the stationary mesh by clocking $2.25$ grids with interpolation; consequently, the same fan-blade position is repeated at every $16$ recorded frames. In this process, the steady-wake component can be subtracted by averaging the flow quantities in the rotating frame, and the fluctuating components can be post-processed in the stationary frame to track their time evolution.

To extract coherent structures and compare them with the eigenfunctions from the linear stability analysis, the velocity fields were recorded on two planes, one plane being on a mid-axial station at $x/D \approx 0.18$ and the other on a ‘quasi’-cylindrical surface (i.e. along iso-radial grid points) near the peak turbulent kinetic energy above the outer wall, as plotted in figure 18. On each plane, an approach equivalent to DMD formulated by Schmid (Reference Schmid2010) is applied to turbulent kinetic energy, which dictates hydrodynamic disturbances (details of the procedure are given in Suzuki et al. Reference Suzuki, Shur, Strelets and Travin2022): after velocity fluctuations are azimuthally decomposed up to $m = 44$ modes (i.e. twice the fan-blade count) on each plane, snapshot proper orthogonal decomposition (POD) modes (Sirovich Reference Sirovich1987) are obtained by optimizing the turbulent kinetic energy, and the time histories of the POD coefficients are then computed. Subsequently, the linear evolution among limited POD coefficients is calculated using a generalized-inverse technique (Perret, Collin & Delville Reference Perret, Collin and Delville2006). If the linear transfer matrix is eigen-decomposed, the sets of the eigenvalue and eigenvector can be regarded as DMD modes except that higher-order POD modes, which are less coherent, have been excluded.

Figure 18.

Surfaces on which DMD modes are extracted. (a) Mid-station plane and quasi-cylindrical surface. (b) Position of the quasi-cylindrical surface relative to turbulent kinetic energy at $x \approx 0.18$.

In addition to the velocity fields, pressure fields on the interstage outer wall are also post-processed mainly to analyse acoustic modes. The wall-pressure fluctuations in the stationary frame are obtained by similarly clocking and interpolating the azimuthal grids, and those at $x/D \approx 0.18$ are azimuthally decomposed using the full $N_{\phi } = 792$ grids over the frequencies. In addition, the snapshot POD analysis is applied to extract coherent modes dominating the wall pressure, and the time histories of the leading POD-mode coefficients are Fourier transformed to calculate their frequency spectra.

5.2. The DMD modes extracted from velocity fields on two planes

We first observe mode structures on the cylindrical surface near the outer wall by taking $m = 23$ at the approach condition, which is considered to be the peak of the azimuthal-mode power in figure 31(a) shown later. Figure 19 exhibits the velocity fluctuations of typical DMD modes that represent eigenfunctions in the first and second unstable regimes. In this study, we count the order of DMD modes so that the first mode contains the greatest turbulent kinetic energy. The first unstable mode, which is the fifth DMD mode in this case, changes the phase rapidly due to a high axial wavenumber. In contrast, the second unstable mode, which is the second DMD mode in this case, appears to grow monotonically with a very small negative wavenumber. For $m = 23$, the first four DMD modes all belong to the second unstable regime, and many higher modes depict wavepacket structures, growing and decaying weakly in the axial direction, relevant to the first unstable regime, as shown in figure 20(b i) below.

Figure 19.

Typical DMD modes near the wall for $m = 23$ at $61.7\,\%$ RPM. Colour schemes are common to all plots. (a) Second DMD mode representing the second unstable regime. (b) Fifth DMD mode representing the first unstable regime. (i) Amplitude. Line patterns are common to all components. (ii) Phase. Negative slope denotes an upstream-propagating mode and vice versa.

Figure 20.

Typical DMD modes near the wall for $m = 18$ at $87.5\,\%$ RPM. Colour schemes and line patterns are the same as figure 19. (a) First DMD mode representing the second unstable regime. (b) Third DMD mode representing the first unstable regime. (i) Amplitude. (ii) Phase.

Likewise, figure 20 shows typical DMD modes at the cut-back condition. Here, we take $m = 18$, which is a local peak in the azimuthally decomposed PWL contour shown later in figure 31(b). The first DMD mode seems to be close to an eigenfunction in the second unstable regime (or a slightly higher frequency/wavenumber). The third DMD mode represents the characteristics of the first unstable regime; in particular, this DMD mode indicates a clear wavepacket structure. Compared with the approach condition, the number of DMD modes in the first regime becomes more dominant than that in the second regime for the cut-back condition, which can be confirmed in figure 21 below.

Figure 21.

Dispersion relations of the first $10$ DMD modes overlaid from $m = 14$ to $41$. $\omega$ and $k_x$ are multiplied by $\times B/m$ (i.e. normalized to $m = 22$); (a) $61.7\,\%$ RPM, (b) $87.5\,\%$ RPM.

Figure 21 overlays the dispersion relations of the first ten DMD modes for each azimuthal mode from $m = 14$ to $41$ extracted from the IDDES at the approach and cut-back conditions. Here, we take the frequency directly from the DMD analysis and calculate the axial wavenumber from the slope of the phase averaged among the three velocity components at $x/D = 0.18$. If there is either a discontinuity in the phase over the axial extent or a discrepancy between the three phase slopes, such a DMD mode is declared as an outlier. To also compare them with the dispersion relation predicted by the compressible linear stability analysis in figures 9 and 11, $\omega$ and $k_x$ are scaled by multiplying them by $B/m$. For reference, the local peak of the azimuthal velocity and the rotor speed are also drawn.

Figure 21 reveals that the IDDES data collapse very well at both approach and cut-back conditions. The linear stability analysis predicts the extracted dispersion relations generally well except for outliers, although it slightly under-estimates the wavenumber for the first unstable regime for an unknown reason. For both engine speeds, we can find one cluster in the second unstable regime (i.e. along the straight line near the origin), and another cluster near the rotor speed, which represents modes convecting with the rotor wakes. In particular, the separation between the first and second regimes is discernible in figure 21(a). In figure 21(b), in contrast, the biggest cluster appears between the first and second regimes. We should recall that the eigenfunction for the first regime forms a peak away from the outer wall at $87.5\,\%$ RPM, as plotted in figure 15(a); hence, the velocity fluctuations on this cylindrical surface may not be exactly representative of the first unstable mode. Nonetheless, figures 20(b) and 21(b) suggest that the first unstable mode carries coherent structures in the rotor wakes without decaying the amplitude.

Next, we observe the radial dependence of the coherent structures. Figure 22 plots the velocity fluctuations of typical leading-order DMD modes for $m = 23$ at the approach condition. In this example, the second and fourth DMD modes take the frequencies corresponding to the first and second unstable regimes, respectively. In both DMD modes, we find large fluctuations near the outer wall, and the difference in the peak radial position is not as clear as that of the eigenfunctions between the first and second modes plotted in figure 14. Since coherent structures between the two unstable regimes seem to interact in the outer-wall boundary layer, separation between them on the same station might be difficult in the post-processing.

Figure 22.

Typical DMD modes on the mid-station plane at $61.7\,\%$ RPM ($m = 23$). (a) Second DMD mode ($\omega = 31.1$). (b) Fourth DMD mode ($\omega = 0.11$). (i) Axial-velocity fluctuation. (ii) Azimuthal-velocity fluctuation.

At the cut-back condition, the coherence becomes weaker than the approach condition; namely the turbulent kinetic energy contained in the leading-order POD modes becomes less, and this makes the extraction of coherent modes challenging. Figure 23 similarly plots the radial profiles of typical DMD modes for $m = 18$ at the cut-back condition. The second DMD mode takes a frequency in the second unstable regime, while the frequency of the third DMD mode is close to the first unstable regime. The latter mode shows relatively large fluctuations in the mid span, which is consistent with the linear stability analysis in figure 15, but also near the outer wall. In many fast co-rotating DMD modes (i.e. large $\omega > 0$), we find similar structures (not shown). In contrast, the former mode forms a peak clearly in the outer-wall boundary layer, as predicted by the linear stability analysis. The number of DMD modes found in the second unstable regime is relatively small at the cut-back condition, while the DMD modes with large fluctuations in the mid span are rather dominant. In fact, figure 18(b) supports this observation.

Figure 23.

Typical DMD modes on the mid-station plane at $87.5\,\%$ RPM ($m = 18$). (a) Second DMD mode ($\omega = 27.1$). (b) Third DMD mode ($\omega = 78.5$). (i) Axial-velocity fluctuation. (ii) Azimuthal-velocity fluctuation.

5.3. Pressure fluctuations on the outer wall

By azimuthally decomposing the wall-pressure fluctuations at the middle of the interstage, we can directly observe the rotor–stator and rotor-related interactions and the cut-on boundaries of the duct acoustic modes. Unlike the inlet/exhaust mode analysis, it is non-trivial to prescribe duct acoustic modes for the mode-extraction method (Ovenden & Rienstra Reference Ovenden and Rienstra2004) due to the swirl velocity, as explained in § 4.3. Thus, we only consider the azimuthal-mode decomposition along a single ring at $x/D \approx 0.18$, and figure 24 displays the decomposed sound pressure level (SPL) as a function of $m$ over the frequencies at the approach and cut-back conditions. The shaded area on the top left corner of figure 24(a) indicates the region in which aliased positive $m$ modes are contaminated due to a low sampling rate. Although non-acoustic components are included on the wall-pressure fluctuations, a V-shape cut-on region can be clearly seen for both conditions. Here, the skewed cut-on boundaries are very well predicted by (4.8) by setting $\chi = 0.35$. High SPL can be observed along these boundaries, especially on the co-rotating side (i.e. $m > 0$), which can also be found in the experiment (Premo & Joppa Reference Premo and Joppa2002). Some SPL can be discerned along the boundary on the counter-rotating side, but these modes are decaying once the cut-on boundaries become symmetric in the inlet or the exhaust.

Figure 24.

Azimuthal-mode decomposition of the wall SPL in the interstage. Horizontal white dotted lines denote the BPF and the harmonics. Skewed V-shape black dashed lines represent the cut-on boundaries predicted by (4.8). Inclined magenta dashed line denotes the rotor speed interacting with the threshold $m_0$; (a) $61.7\,\%$ RPM ($m_0 = -39$), (b) $87.5\,\%$ RPM ($m_0 = -57$).

Importantly, fluctuations are noticeable even outside the cut-on zone along multiple lines from the left bottom to the right top with the slope representing the rotor speed. When a disturbance of the azimuthal mode $n_B$ spinning at the rotor speed $\varOmega$ is interacting with a quasi-stationary object of the azimuthal mode $m_0$ moving at $\varOmega _0$ ($\approx 0$), it creates the azimuthal mode of $m = n_B + m_0$; hence, the mode line appears in figure 24 along

(5.1)\begin{equation} f\, ({\rm Hz}) = \frac{n_B \varOmega - m_0 \varOmega_0}{2 {\rm \pi}} = \frac{\varOmega}{2 {\rm \pi}} \left( m - m_0 \right) - \frac{\varOmega_0}{2 {\rm \pi}} m_0 \approx \frac{\varOmega}{2 {\rm \pi}} \left( m - m_0 \right) . \end{equation}

If multiple quasi-stationary structures of $m_0$ are moving with different speeds of $\varOmega _0$, the second term $\varOmega _0 m_0 / ( 2 {\rm \pi})$ becomes obscure, and lines only along the rotor speed remain to be noticeable. Namely, figure 24 suggests that multiple quasi-stationary structures are interacting with broadband motions spinning at the rotor speed. The threshold azimuthal-mode count $m_0$ drawn by the magenta dashed line approximately corresponds to the cut-off wavenumber of the VSTG ($| m_0 | \lesssim 39$ and $\lesssim 57$ at the approach and the cut-back conditions, respectively). When the $n-1$th harmonic of the fan-blade count $B$ (i.e. $n_B = nB$) interacts with the $k_n-1$th harmonic of the OGV count $V_n$ (i.e. $m_0 = k_n V_n$ and $\varOmega _0 = 0$), it yields the Tyler & Sofrin (Reference Tyler and Sofrin1962) relation, generating the $nB + k_n V_n$ azimuthal mode at the frequency $n B \varOmega / ( 2 {\rm \pi})$, where the BPF is given by $f_{{ BPF}} \equiv {B \varOmega }/{( 2 {\rm \pi})}$.

To extract coherent modes among the wall-pressure fluctuations, the snapshot POD method (Sirovich Reference Sirovich1987) is applied in the physical domain. On the $x - \phi$ domain in the stationary frame, snapshot POD modes on the outer wall are unwrapped, and a few representative modes are displayed from figures 25 to 28. The azimuthal-mode distribution averaged over $x$ is also plotted for each POD mode. Moreover, time histories of the POD-mode coefficients are calculated, and a frequency spectrum of each POD mode is given as well. Although the sampling rate was low at the approach condition, it is still meaningful to show it beyond the Nyquist frequency; thus, the spectrum including both sides (i.e. symmetric spectrum across $3f_{{ BPF}}$) being depicted in figures 25(c) and 26(c). At the approach condition, the first POD mode captures the rotor–stator interaction tones; namely, 2BPF with $m = -10$ and 4BPF with $m = 34$. The second POD mode (not shown) is very similar to the first POD mode. The fourth (and sixth, not shown) POD mode has a peak around $2.8 f_{{ BPF}}$ in the broadband with the azimuthal peak of $m = 31$. We can find the corresponding spot along a stripe of the broadband noise in figure 24(a), suggesting that the aforementioned stripes have strong coherency.

Figure 25.

First snapshot POD mode of the wall-pressure fluctuation at $61.7\,\%$ RPM. (a) Unwrapped pressure contour. (b) Azimuthal-mode distribution averaged over $x$. (c) Frequency spectrum of the POD-mode coefficient.

Figure 26.

Fourth snapshot POD mode of the wall-pressure fluctuation at $61.7\,\%$ RPM. Same layout as figure 25.

Figure 27.

Third snapshot POD mode of the wall-pressure fluctuation at $87.5\,\%$ RPM. (a) Unwrapped pressure contour. (b) Azimuthal-mode distribution averaged over $x$. (c) Frequency spectrum of the POD-mode coefficient.

Figure 28.

Fifth snapshot POD mode of the wall-pressure fluctuation at $87.5\,\%$ RPM. Same layout as figure 27.

Likewise, figures 27 and 28 show representative snapshot POD modes at the cut-back condition. Unlike the approach condition, the sampling rate was high enough so that each figure (c) plots one side of the spectrum. The third (and fourth, not shown) POD mode captures the 2BPF tone with $m = -10$ in figure 27 (although not shown, the first two POD modes take the BPF tone with $m = -32$). Compared with the approach condition, these tone intensities are clearly higher downstream, indicating the conventional rotor–stator interaction. In contrast, the fifth (and sixth, not shown) POD mode represents broadband noise around $4.2 f_{{BPF}}$ with $m \approx 66$ in figure 28. This peak appears as a strong SPL spot in figure 24(b) along a stripe. We should remark that when the POD mode is found close to the cut-on boundary, i.e. non-dimensional wavenumber of $k_x \ll 1$ so that the unwrapped pressure contour depicts almost horizontal stripes regardless of tones or broadband noise. In this sense, the wavenumber structures of these acoustic modes are somewhat similar to the aforementioned second unstable modes. Yet, because of the difference in the frequency range, these POD modes can be confirmed as the acoustic modes (note that the pressure field in the aforementioned quasi-cylindrical surface is dominated by these acoustic modes).

The results above suggest that the wall-pressure fluctuations in the interstage are largely influenced by the rotor–stator interaction tones satisfying the Tyler–Sofrin relation (Tyler & Sofrin Reference Tyler and Sofrin1962) as well as the interaction between disturbances moving with the rotor speed and those generated downstream of the rotor (likely spinning slowly). The latter disturbances may be governed by the instability in the second regime, which is found to dominate coherent hydrodynamic structures near the outer wall. Within the skewed cut-on boundaries of the interstage flow, strong coherent acoustic modes can be found in the broadband noise as a result of such an interaction. Although it has been revealed by Golubev & Atassi (Reference Golubev and Atassi1998) that acoustic modes possess a vortical component and hydrodynamic modes possess pressure fluctuations under swirl velocity, which has been confirmed in figures 10 and 12, such couplings do not seem to be the dominant interaction mechanism. These broadband noise generation seems to be rather caused by the nonlinear product of two different types of non-acoustic disturbances.

5.4. Contribution to fan broadband noise

Finally, we relate the stripe patterns with the sound properties radiating from the inlet and the exhaust. Figures 29 and 30 display the azimuthal-mode decomposition of the inlet/exhaust duct acoustic-mode PWL at the approach and cut-back conditions, respectively. Suzuki et al. (Reference Suzuki, Spalart, Shur, Strelets and Travin2019) described the procedure of the mode-extraction method (Ovenden & Rienstra Reference Ovenden and Rienstra2004) specifically tailored for the fan-noise simulations. In general, these IDDES results are consistent with the azimuthal-mode decomposition of the wall-pressure fluctuations measured by Premo & Joppa (Reference Premo and Joppa2002). A particular interest here is the stripe patterns drawn by the inclined lines: as explained in § 5.3, the slope of these lines represents the rotor speed, and the azimuthal-mode numbers at the root of the lines (i.e. at $0$ Hz) indicate the mode counts of the quasi-stationary structures. The magenta line corresponds to the threshold wavenumber of the VSTG, introduced in § 3 and figure 24 ($| m_0 | \le 39$ and $\le 57$ at the approach and cut-back conditions, respectively). The yellow lines represent the local azimuthal peaks of the inlet/exhaust PWL spectra identified in figure 31 below.

Figure 29.

Azimuthal-mode decomposition of the duct-mode PWL at $61.7\,\%$ RPM. Horizontal white lines and inclined magenta line are the same as figure 24. Yellow inclined lines denote the azimuthal peaks identified in figure 31. (a) Inlet. (b) Exhaust.

Figure 30.

Azimuthal-mode decomposition of the duct-mode PWL at $87.5\,\%$ RPM. (a) Inlet. (b) Exhaust.

Figure 31.

Integrated PWL along the inclined lines on the azimuthal PWL contours in the inlet and the exhaust (scaled by dB on the left axis). Integrated SPL of the azimuthal SPL contours in the interstage is also plotted (offset by $25$ dB); (a) $61.7\,\%$ RPM, (b) $87.5\,\%$ RPM.

By integrating the PWL along these inclined lines, the radiating sound power can be calculated as a function of the azimuthal-mode number associated with the potentially stationary structures (i.e. $m$ values at $0$ Hz). In figure 31, such PWLs are compared between the inlet and the exhaust together with the SPL on the interstage wall similarly integrated along the inclined lines in figure 24 (the interstage SPL is offset by $25$ dB). Here, similar azimuthal-mode patterns are depicted between the inlet and the exhaust as well as the interstage, particularly at the approach condition. At the cut-back condition, the PWL is rapidly decreased toward greater $-m$ in the inlet due to the transonic rotor speed, while the PWL is spread over a wide azimuthal-mode range in the exhaust and the interstage. More importantly, many common peaks, indicated by the vertical dotted lines, can be identified among the three regions, but they are smeared in the inlet and the exhaust at the cut-back condition. This suggests an existence of a broadband noise-generation mechanism at both approach and cut-back conditions, and it contributes commonly to the inlet and exhaust PWL spectra. This noise-generation mechanism appears to play a major role in the fan broadband noise at the approach condition but the relative importance becomes less at the cut-back condition, particularly in the inlet. By comparing figures 29(a) with 30(a), it is clear that the noise directly related to the transonically spinning rotor, which is distributed along the cut-on line on the co-rotating side, is significantly increased at the cut-back condition.

To crudely assess the magnitude of the PWL contained in the stripes of the IDDES results, the following steps are taken (for a detailed procedure, refer to Suzuki et al. Reference Suzuki, Shur, Strelets and Travin2022, particularly figure 10): first, the azimuthally decomposed PWL along the inclined lines in figure 29 or 30 is filtered using the five-point smoothing (Savitzky & Golay Reference Savitzky and Golay1964) in frequency. Second, troughs are identified in the azimuthal distribution of the PWL within the azimuthal bundle confined by the aforementioned threshold $m_0$ of the VSTG. By linearly connecting the troughs together with the PWL distribution outside the azimuthal range, the background PWL is then defined, and the logarithmic difference between the total and the background PWLs is eventually treated as the segregated PWL in the stripes. Thus, the extra contribution in the stripes is graphically extracted.

Figure 32 plots the breakdown of the spectra from the mode extraction at the approach condition by summing the background and segregated PWL over the azimuthal modes. Fluctuations in the PWL spectra are somewhat high even in the total spectra due to a narrow bandwidth to compare with the test spectra in the series of this study; however, these figures capture the trends of the PWL in the stripes. In the inlet, the segregated spectrum in the stripes appears to have comparable PWL to the rest of the background PWL, which should characterize the conventional rotor–stator interaction noise. In the exhaust, the PWL in the stripes steeply declines at higher frequencies, partially because of the greater total PWL in the exhaust at higher frequencies. As visually observed, the contribution from the stripes seems to be more relevant to the inlet broadband noise.

Figure 32.

Breakdown of the PWL spectra at $61.7\,\%$ RPM assessed by the steps given by Suzuki et al. (Reference Suzuki, Shur, Strelets and Travin2022). Line patterns: $-~-~-$, total; $\cdot \cdot \cdot \cdot \cdot$, background; —–, segregated. (a) Inlet. (b) Exhaust.

Likewise, figure 33 plots the breakdown at the cut-back condition. The contribution from the segregated spectrum relative to the total PWL spreads over wide frequencies, but it is typically lower by $6$ dB or more; hence the relative contribution of this noise-generation mechanism seems to be weaker at $87.5\,\%$ RPM than $61.7\,\%$ RPM, which can be observed by comparing figures 29 and 30. Although the procedure of this breakdown is nothing but a pattern recognition, these results suggest that such a noise-generation mechanism appreciably contributes to the fan broadband noise at the approach condition, but becomes less important at the cut-back condition.

Figure 33.

Breakdown of the PWL spectra at $87.5\,\%$ RPM assessed by the steps given by Suzuki et al. (Reference Suzuki, Shur, Strelets and Travin2022). Line patterns are the same as figure 32. (a) Inlet. (b) Exhaust.