2. Numerical methodology

We use the numerical database of the wake of a circular disk at $ {\textit {Re}} = 5 \times 10^4$ from CS2020. In particular, we analyse the datasets of stratified wakes at $ {\textit {Fr}} = 2$ and $10$ from their numerical database. CS2020 use high-resolution LES to numerically solve the filtered Navier–Stokes equations system along with density diffusion equation under the Boussinesq approximation.

These equations are as follows:

continuity,

(2.1)\begin{equation} \frac{\partial u_{i}}{\partial x_{i}} = 0; \end{equation}

momentum,

(2.2)\begin{equation} \frac{\partial u_{i}}{\partial t} + \frac{\partial (u_{i}u_{j})}{\partial x_{j}} ={-}\frac{\partial p}{\partial x_{i}} + \frac{1}{ {\textit{Re}}}\frac{\partial}{\partial x_{j}}\left[\left(1 + \frac{\nu_{s}}{\nu}\right)\frac{\partial u_{i}}{\partial x_{j}}\right] - \frac{1}{ {\textit{Fr}}^{2}}\rho'\delta_{i3}; \end{equation}

and density diffusion,

(2.3)\begin{equation} \frac{\partial \rho}{\partial t} + \frac{\partial (\rho u_{j})}{\partial x_{j}} = \frac{1}{ {\textit{Re}} {\textit{Pr}}}\frac{\partial}{\partial x_{j}}\left[\left(1 + \frac{\kappa_{s}}{\kappa}\right)\frac{\partial \rho}{\partial x_{j}}\right], \end{equation}

where $u_{i}$ corresponding to $i=1, 2$ and $3$ refer to velocity in the streamwise ($x_{1}$ or $x$), lateral ($x_{2}$ or $y$) and vertical ($x_{3}$ or $z$) directions, respectively. Gravity acts in the vertical direction (2.2). The density field is decomposed into a background profile, $\rho _b(z) = \rho _o + (\mathrm {d} \rho _b/\mathrm {d} z)z$ (where ${\rm d}\rho_b/{\rm d}z$ is a constant), and density deviation ($\rho '$). Thus, $\rho (x,y,z,t) = \rho _{b}(z) + \rho '(x,y,z,t)$. In (2.2), $\nu _{s}$ and $\nu$ refer to the subgrid kinematic viscosity obtained from LES and kinematic viscosity of the fluid, respectively. Likewise, $\kappa _{s}$ and $\kappa$ in (2.3) refer to the subgrid density diffusivity and density diffusivity of the fluid, respectively.

Equations (2.1)–(2.3) are non-dimensionalized using the following parameters: (i) free-stream velocity ($U_{\infty }$) for the velocity field, (ii) diameter of disk ($D$) for spatial locations $x_{i}$, (iii) dynamic pressure ($\rho _{o}U_{\infty }^{2}$) for pressure field, (iv) advection time scale ($D/U_{\infty }$) for time $t$, and (iv) $-(\mathrm {d} \rho _b/\mathrm {d} z)D$ for density deviation. There are three non-dimensional parameters of interest: (1) body-based Reynolds number ($ {\textit {Re}}$) defined as $U_{\infty }D/\nu$; (2) body-based Froude number ($ {\textit {Fr}}$) defined as $U_{\infty }/ND$, where $N$ is the buoyancy frequency, $N^{2} = -g/\rho _{o}(\mathrm {d} \rho _b/\mathrm {d} z)$; and (3) Prandtl number ($ {\textit {Pr}}$) defined as $\nu /\kappa$ which is set as 1 in CS2020 simulations. Here $\kappa _{s}$ is also set equal to $\nu _{s}$ for the LES simulations.

A cylindrical coordinate system is adopted and the disk is represented using the immersed boundary method of Balaras (Reference Balaras2004) and Yang & Balaras (Reference Yang and Balaras2006). Spatial derivatives are computed using second-order central finite differences and temporal marching is performed using a fractional step method which combines a low-storage Runge–Kutta scheme (RKW3) with the second-order Crank–Nicolson scheme. The kinematic subgrid viscosity ($\nu _{s}$) and density diffusivity ($\kappa _{s}$) are obtained using the dynamic eddy viscosity model of Germano et al. (Reference Germano, Piomelli, Moin and Cabot1991). At the inlet and outlet, Dirichlet inflow and Orlanski-type convective (Orlanski Reference Orlanski1976) boundary conditions are specified, respectively. The Neumann boundary condition is used at the radial boundary for the density and velocity fields. To prevent the spurious propagation of internal waves upon reflection from the boundaries, sponge regions with Rayleigh-damping are employed at radial, inlet and outlet boundaries.

The radial and streamwise domains span $0 \leqslant r/D \leqslant 80$ and $-30 \leqslant x/D \leqslant 125$, respectively. A large radial extent facilitates weakening of the IGWs before they hit the boundary and thereby also controls the amplitude of spurious reflected waves. The distribution of grid points are as follows: $N_{r} = 531$ in the radial direction, $N_{\theta } = 256$ in the azimuthal direction and $N_{x} = 4608$ in the streamwise direction, resulting in approximately 530 million elements. The grid resolution is excellent by LES standards in all three directions. Readers may refer to CS2020 for more details on the grid resolution and numerical scheme.

3. Spectral POD – theory and present application

In this work we employ SPOD to study the dynamics of coherent structures in stratified wakes, rather than the more commonly employed snapshot POD (Sirovich Reference Sirovich1987). Spectral POD enables the identification of dominant structures evolving coherently in both space and time by exploiting temporal correlation among flow snapshots. This approach is particularly well suited for flow configurations like turbulent wakes which are known to be dominated by mechanisms operating at specific frequencies, e.g. VS, pumping of recirculation bubble, shear layer breakdown, to name a few (Berger, Scholz & Schumm Reference Berger, Scholz and Schumm1990). On the contrary, snapshot POD assumes each snapshot of the flow to be an independent realization. As a result, the temporal coherence of POD modes is not guaranteed. Furthermore, it can be also shown that the coefficients dictating the temporal evolution of snapshot POD modes are broadband, i.e. containing contributions from a range of frequencies (Towne et al. Reference Towne, Schmidt and Colonius2018). Spectral POD requires a larger amount of time-resolved data compared with snapshot POD. Hence, snapshot POD has dominated the literature compared with SPOD.

3.1. Theory of SPOD for statistically stationary stratified flows

For the SPOD analysis of stratified wakes, the fluctuating density fields $(\rho '(\boldsymbol {x},t))$ and velocity fields $(\boldsymbol {u}'(\boldsymbol {x},t) = [{u}_{r}'(\boldsymbol {x},t),{u}_{\theta }'(\boldsymbol {x},t),{u}_{x}'(\boldsymbol {x},t)])$ are taken together as a single state-space field $\boldsymbol {\varLambda }(\boldsymbol {x},t) = [\boldsymbol {u}'(\boldsymbol {x},t), \rho '(\boldsymbol {x},t)]^{{\rm T}}$. Following Lumley (Reference Lumley1970), we seek POD modes $\boldsymbol {\varPsi }(\boldsymbol {x},t)$ that have maximum ensembled-average projection on $\boldsymbol {\varLambda }(\boldsymbol {x},t)$, expressed as

(3.1)\begin{equation} \max_{\boldsymbol{\varPsi}} \frac{\langle |\{\boldsymbol{\varLambda}(\boldsymbol{x},t), \boldsymbol{\varPsi}(\boldsymbol{x},t)\}|^{2}\rangle}{\{\boldsymbol{\varPsi}(\boldsymbol{x},t),\boldsymbol{\varPsi}(\boldsymbol{x},t)\}}, \end{equation}

where $\langle \cdot \rangle$ denotes the ensemble average. We define the inner product $\{\boldsymbol {\varLambda }^{(1)}(\boldsymbol {x},t), \boldsymbol {\varLambda }^{(2)}(\boldsymbol {x},t)\}$ as

(3.2)\begin{equation} \{\boldsymbol{\varLambda^{(1)}}(\boldsymbol{x},t), \boldsymbol{\varLambda^{(2)}}(\boldsymbol{x},t)\} = \int_{-\infty}^{\infty}\int_{\varOmega} \boldsymbol{\varLambda^{(2)*}}(\boldsymbol{x},t)\,\mathrm{diag}\left(1,1,1,\frac{g^{2}}{\rho^{2}_{o}N^{2}} \right)\boldsymbol{\varLambda^{(1)}}(\boldsymbol{x},t)\,\mathrm{d}\kern0.7pt\boldsymbol{x}\,\mathrm{d} t, \end{equation}

where $(.)^{*}$ denotes the Hermitian transpose. The so-defined inner-product norm ensures that the obtained POD modes are optimal in terms of capturing two times the overall sum of turbulent kinetic energy (TKE) and turbulent potential energy (TPE), where ${\rm TKE} = \langle u'_{i}u'_{i} \rangle /2$ and ${\rm TPE}={g^{2}}/{2\rho ^{2}_{o}N^{2}}\langle \rho '\rho '\rangle$.

Following Holmes et al. (Reference Holmes, Lumley, Berkooz and Rowley2012), (3.1) can be expressed as a Fredholm-type integral eigenvalue problem as

(3.3)\begin{equation} \int_{-\infty}^{\infty}\int_{\varOmega}R_{ij}(\boldsymbol{x},\boldsymbol{x}',t,t')\boldsymbol{W}(\boldsymbol{x}')\varPsi^{(n)}_{j}(\boldsymbol{x}',t')\,\mathrm{d}\boldsymbol{x}'\, \mathrm{d} t' = \lambda^{(n)} \varPsi^{(n)}_{i}(\boldsymbol{x},t), \end{equation}

where $\boldsymbol {W}(\boldsymbol {x})$ is a positive-definite Hermitian matrix accounting for the weights of each variable as defined in (3.2). In (3.3), $\lambda ^{(n)}$ and $\varPsi ^{(n)}_{i}(\boldsymbol {x},t)$ correspond to the $n{\rm th}$ eigenvalue and the $i{\rm th}$ component of the $n{{\rm th}}$ eigenmode. The kernel $R_{ij}(\boldsymbol {x}, \boldsymbol {x}',t, t')$ which is the two-point two-time correlation tensor, is defined as

(3.4)\begin{gather} R_{ij}(\boldsymbol{x}, \boldsymbol{x}',t, t') = \langle u'_{i}(\boldsymbol{x},t) u'_{j}(\boldsymbol{x}',t') \rangle, \quad i, j = 1, 2, 3, \end{gather}

(3.5)\begin{gather}R_{i4}(\boldsymbol{x}, \boldsymbol{x}',t, t') = \langle u'_{i}(\boldsymbol{x},t) \rho'(\boldsymbol{x}',t') \rangle, \quad i = 1, 2, 3, \end{gather}

(3.6)\begin{gather}R_{4j}(\boldsymbol{x}, \boldsymbol{x}',t, t') = \langle \rho'(\boldsymbol{x},t) u'_{j}(\boldsymbol{x}',t') \rangle, \quad j = 1, 2, 3, \end{gather}

(3.7)\begin{gather}R_{44}(\boldsymbol{x}, \boldsymbol{x}',t, t') = \langle \rho'(\boldsymbol{x},t) \rho'(\boldsymbol{x}',t') \rangle. \end{gather}

For statistically stationary flows, such as the turbulent stratified wake in the present case, the kernel $R_{ij}(\boldsymbol {x}, \boldsymbol {x}',t, t')$ is only a function of time difference $\tau = t-t'$, $\boldsymbol {x}$ and $\boldsymbol {x}'$. Furthermore, it can be Fourier transformed in the temporal direction as

(3.8)\begin{equation} R_{ij}(\boldsymbol{x}, \boldsymbol{x}',\tau) = \int_{-\infty}^{\infty} S_{ij}(\boldsymbol{x}, \boldsymbol{x'}, f){\rm e}^{{\rm i}2{\rm \pi} f\tau}\,\mathrm{d} f, \end{equation}

where $S_{ij}(\boldsymbol {x},\boldsymbol {x}', f)$ is the Fourier transform of the kernel $R_{ij}(\boldsymbol {x}, \boldsymbol {x}',\tau )$. Using (3.8), the Fredholm-type eigenvalue problem in (3.3) can be transformed into an equivalent eigenvalue problem which is solved at each frequency $f$, following Towne et al. (Reference Towne, Schmidt and Colonius2018),

(3.9)\begin{equation} \int_{\varOmega}S_{ij}(\boldsymbol{x},\boldsymbol{x}',f)\boldsymbol{W}(\boldsymbol{x}') \varPhi^{(n)}_{j}(\boldsymbol{x}',f)\,\mathrm{d}\boldsymbol{x}' = \lambda^{(n)}(f) \, \varPhi^{(n)}_{i}(\boldsymbol{x},f), \end{equation}

where $\lambda ^{(n)}(f)$ are the eigenvalues at $f$ and $\varPhi ^{(n)}_{i}(\boldsymbol {x},f) = \varPsi ^{(n)}_{i}(\boldsymbol {x},t)e^{-i2{\rm \pi} f t}$ are the modified eigenmodes. The eigenvalues are ordered such that $\lambda ^{(1)}(f) \geqslant \lambda ^{(2)}(f) \geqslant \dots \geqslant \lambda ^{(n)}(f)$. The sum over all the eigenvalues at frequency $f$ equates to two times the total fluctuation energy content, i.e. $\langle u'_iu'_i \rangle + {g^{2}}/{\rho ^{2}_{o}N^{2}}\langle \rho '\rho '\rangle$ at that frequency. The obtained eigenmodes in the frequency space are spatially orthogonal to each other such that

(3.10)\begin{equation} \int_{\varOmega} \boldsymbol{\varPhi}^{*(n)}(\boldsymbol{x},f)\boldsymbol{W}(\boldsymbol{x})\boldsymbol{\varPhi}^{(m)}(\boldsymbol{x},f)\,\mathrm{d} \boldsymbol{x} = \delta_{mn}, \end{equation}

where $\delta _{mn}$ is the Dirac-delta function.

3.2. Numerical implementation of SPOD for current work

In this work we primarily present results from SPOD on 2-D planes at various $x/D$ – ranging from $x/D = 10$ to $100$ – sampled at a spacing of approximately $5D$. The domain of $10\leqslant x/D\leqslant 100$ spans: (i) $5 \leqslant Nt_2 \leqslant 50$ for $ {\textit {Fr}} = 2$ and (ii) $1 \leqslant Nt_{10} \leqslant 10$ for $ {\textit {Fr}} = 10$ in terms of buoyancy time. In the radial direction the SPOD domain spans $0 \leqslant r/D \leqslant 10$, resulting in a total of $N^{SPOD}_r = 333$ points. In the azimuthal direction the number of grid points is $N_{\theta } = 256$.

For numerical implementation, the mean-subtracted data, consisting of $N$ temporal snapshots, is divided into $N_{blk}$ blocks with an overlap of $N_{ovlp}$ snapshots. Each block contains $N_{freq}$ entries: $\boldsymbol {Q} = [\boldsymbol {q}^{(1)}, \boldsymbol {q}^{(2)}, \boldsymbol {q}^{(3)}, \ldots \boldsymbol {q}^{(N_{freq})}]$. Here, $\boldsymbol {q}^{(i)} = [\boldsymbol {u'}^{(i)}, \rho '^{(i)}]^{{\rm T}}$, where $\boldsymbol {u'}$ and $\rho '$ are velocity and density fluctuations, respectively. Thereafter, discrete Fourier transform of each block is performed in the temporal direction and the ensemble of $N_{blk}$ Fourier realizations of any given frequency, let us say $f$, is collected as $\hat {\boldsymbol {Q}}_f = [\boldsymbol {q}^{(1)(f)}, \boldsymbol {q}^{(2)(f)}, \boldsymbol {q}^{(3)(f)}, \ldots \boldsymbol {q}^{(N_{blk})(f)}]$. Once $\hat {\boldsymbol {Q}}_{f}$ is obtained, SPOD eigenvalues and eigenvectors corresponding to $f$ are given by the following eigenvalue decomposition:

(3.11)\begin{equation} \hat{\boldsymbol{Q}}^{*}_{f}\boldsymbol{W}\hat{\boldsymbol{Q}}_{f} \boldsymbol{\varGamma}_{f} = \boldsymbol{\varGamma}_{f}\boldsymbol{\varLambda}_{f}. \end{equation}

Here $\boldsymbol {\varLambda }_{f} = {\rm diag}(\lambda ^{(1)}_{f},\lambda ^{(1)}_{f}, \ldots \lambda ^{(N_{blk})}_{f})$ is a diagonal matrix containing eigenvalues ranked in the decreasing order of energy content from $i = 1$ to $N_{blk}$. The corresponding spatial eigenmodes $\hat {\boldsymbol {\varPhi }}_{f}$ can be obtained as $\hat {\boldsymbol {\varPhi }}_{f} = \hat {\boldsymbol {Q}}_{f}\boldsymbol {\varGamma }_{f}\boldsymbol {\varLambda }_{f}^{-1/2}$. In (3.11), $\boldsymbol {W}$ is a diagonal matrix of size $4N^{SPOD}_{r}N_{\theta }$, containing the numerical quadrature weights multiplied by coefficients required to form the energy quantities given in (3.2).

The parameters for SPOD are set as follows: (i) total number of snapshots $N=7168$ with consecutive snapshots separated by $\Delta tD/U_{\infty } \approx 0.09$ and $0.104$ for $ {\textit {Fr}} = 2$ and $10$, respectively; (ii) number of frequencies $N_{freq} = 512$; and (iii) overlap between blocks $N_{ovlp} = 256$, resulting in total of $N_{blk} = ({N-N_{ovlp}})/({N_{freq} - N_{ovlp}}) = 27$ SPOD modes at each frequency. Interested readers are referred to Towne et al. (Reference Towne, Schmidt and Colonius2018) and Schmidt & Colonius (Reference Schmidt and Colonius2020) for more details on the theoretical aspects and numerical implementation of SPOD.

Most of the results are obtained from SPOD analyses at constant $x/D$ planes with modes maximizing the two times sum of TKE and TPE. However, for some results, we perform additional SPOD analyses. For example, to illustrate the streamwise variation of a certain leading-order SPOD mode in § 7, we perform SPOD analysis on fluctuating velocity and density fields at the centre-vertical plane ($y=0$ plane) with a reduced number of snapshots $N=5376$ and half-resolution in vertical and streamwise directions. Here $N_{freq}$ and $N_{ovlp}$ are kept the same as SPOD on fixed $x/D$ planes. The spatial resolution and $N$ are reduced to avoid memory limitations since large matrices with complex double precision have to be stored in the intermediate steps of SPOD. Also in § 6 we present results from SPOD analyses of the $ {\textit {Fr}} =2$ wake (at constant $x/D$ planes) with (i) density fluctuations replaced by pressure fluctuations, and (ii) norm defined such as to maximize the sum of $\langle p'p'\rangle$ and $\langle u'_iu'_i \rangle$. Here $N,N_{freq}$ and $N_{blk}$ are kept the same as in the previous paragraph. The motivation behind performing this additional set of SPOD analyses is explained in § 6.

4. Flow visualizations

Three-dimensional visualizations of the $Q$ criterion and planar views of the vorticity and velocity fields in this section provide a first look at the vortical and unsteady IGW structure of the simulated wakes. The structure of the steady (in a frame attached to the disk) lee wave field is not discussed in this paper. To emphasize the large-scale coherent structures, the instantaneous velocity fields have been filtered using a SciPy Gaussian low-pass filter (Gaussian_filter) in all three directions with standard deviation $\sigma = 5$ before calculating the $Q$ criterion and vorticity fields.

Figure 1 shows that, in both wakes, circular vortex rings appear immediately downstream of the disk. At $ {\textit {Fr}} = 2$, the buoyancy-induced anisotropy between horizontal and vertical directions commences in the near wake. The wake contracts in the vertical at $x/D \approx 5$ (visible in the side view given in figure 1b) owing to the oscillatory modulation by the lee wave. The top view (figure 1a) shows a distinct large-scale waviness in the intermediate wake, shown by the dashed black line. Its approximate wavelength is $\lambda /D \approx 1/ {\textit {St}}_{VS}$, where $ {\textit {St}}_{VS}$ is the VS frequency. Likewise, large-scale VS structures separated approximately by $\lambda /D \approx 1/ {\textit {St}}_{VS}$ can also be identified in the $ {\textit {Fr}} = 10$ wake (figure 1c). The value of $ {\textit {St}}_{VS}$ and the spatial behaviour of the VS mode will be made precise formally using SPOD in the subsequent sections.

Figure 1. Isosurfaces of instantaneous $Q$ criterion at $Q = 0.01$: (a,b) $ {\textit {Fr}} = 2$ and (c) $ {\textit {Fr}} = 10$. Streamwise domain is limited to $0 < x/D < 20$ for clarity.

Figure 2 shows the instantaneous vertical vorticity ($\omega _{z}D/U_{\infty }$) on the central horizontal plane ($z = 0$) for the $ {\textit {Fr}} = 2$ (a) and $ {\textit {Fr}} = 10$ (b) wakes. Similar to figure 1, $\omega _z$ is calculated using filtered velocity fields to emphasize large-scale features. In both wakes, the complex spatial distribution of vorticity in the immediate downstream of the disk gives way to a well-defined coherent distribution of opposite signed vortices in the intermediate to late wakes. For the $ {\textit {Fr}} = 2$ wake, spatial coherence is visible as early as $x/D \approx 20$. Beyond $x/D \approx 20$, the regions of opposite signed $\omega _{z}$ remain separated till the end of the domain. On closer inspection, a streamwise undulation of length $\lambda /D \approx 1/ {\textit {St}}_{VS}$ can be observed in figure 2(a). At this point, it is important to emphasize that the $ {\textit {Fr}} = 2$ wake remains actively turbulent throughout the computational domain as demonstrated by CS2020 through spectra and visualizations of the turbulent dissipation rate. From $x/D \approx 40$ onward, the $ {\textit {Fr}} =2$ wake resides in the SST regime. Different regimes of stratified turbulence are discussed briefly in § 1. The strong signature of coherence in the $ {\textit {Fr}} = 2$ wake is not a consequence of the transition into the weakly turbulent state of the Q2D regime noted in previous works, e.g. by Spedding (Reference Spedding1997).

Figure 2. Instantaneous snapshot of vertical vorticity at the central horizontal plane ($z = 0$): (a) $ {\textit {Fr}} = 2$ and (b) $ {\textit {Fr}}=10$.

The $ {\textit {Fr}} = 10$ wake also shows a distinct wavy motion with non-dimensional wavelength ${\approx }1/ {\textit {St}}_{VS}$. However, the separation between the regions with opposite signed vorticity is not as well defined as in the $ {\textit {Fr}} = 2$ wake. According to CS2020, the $ {\textit {Fr}} = 10$ wake stays in the weakly stratified regime (WST) from $x/D \approx 10$ to $50$ and thereafter stays in the intermediately stratified regime (IST) till the end of the domain.

To conclude this section, instantaneous snapshots of fluctuating spanwise velocity ($u'_{y}/U_{\infty }$) are shown in figure 3 at locations in the near, intermediate and far wake at $ {\textit {Fr}} = 2$ (a–c) and $ {\textit {Fr}} = 10$ (d–f). An ellipse with major and minor axes equal to $2L_{Hk}$ and $2L_{Vk}$, where $L_{Hk}$ and $L_{Vk}$ are the TKE-based wake widths in the horizontal and vertical directions, respectively, is also shown. Here $L_{Hk}$ is defined by $\mathrm {TKE}(x,y=L_{Hk},z=0) =\mathrm {TKE}(x,r=0)/2$ and $L_{Vk}$ by $\mathrm {TKE}(x,y=0,z=L_{Vk}) =\mathrm {TKE}(x,r=0)/2$. The disk centreline is denoted by $r=0$. It is worth noting that using the sum of TKE and TPE to define the wake widths (not shown here) result in values similar to $L_{Hk}$ and $L_{Vk}$ for both $ {\textit {Fr}} = 2$ and $ {\textit {Fr}} = 10$ wakes. Following CS2020, we use the TKE-based definitions in the rest of the results and discussions. This ellipse, based on $L_{Hk}$ and $L_{Vk}$, is used to approximately demarcate the wake core from the outer wake. In subsequent sections this definition of the wake core will prove to be useful for the interpretation of some SPOD results.

Figure 3. Instantaneous snapshots of the fluctuating spanwise velocity $u'_{y}/U_{\infty }$ shown for $ {\textit {Fr}} = 2$ (a–c) and $ {\textit {Fr}} = 10$ (d–f): (a,d) at $x/D = 10$, (b,e) at $x/D = 50$ and (c,f) at $x/D = 100$. Dashed close curve in white shows wake core.

At $ {\textit {Fr}} = 2$, an appreciable effect of buoyancy is already present in the near wake as shown in figure 3(a) for $x/D = 10$, which corresponds to $Nt_{2} = 5$ in buoyancy time units. At the same streamwise location, the $ {\textit {Fr}} = 10$ wake still has a circular cross-section with an imprint of the $m=1$ azimuthal mode which was found to be energetically important in the unstratified wake (Nidhan et al. Reference Nidhan, Chongsiripinyo, Schmidt and Sarkar2020). As both the wakes evolve downstream, buoyancy has a progressively increasing effect on the wake core as well as the surrounding outer wake. By $x/D = 50$, vertically flattened wake cores can be observed in figure 3(b,e) for both the wakes, more so at $ {\textit {Fr}} = 2$ than at $ {\textit {Fr}} = 10$. It is also worth noting that the wake core of $ {\textit {Fr}} = 2$ consists of distinct layers by $x/D = 50$. The $ {\textit {Fr}} =2$ wake also shows a significant amount of IGW activity in the outer region, i.e. outside the ellipse in figure 3(b). Farther downstream at $x/D = 100$, the $u'_{y}$ field of $ {\textit {Fr}} = 2$ (figure 3c) shows IGWs occupying a significant portion of the outer wake with the wake core being further flattened and comprising an increased number of horizontally oriented layers. The $ {\textit {Fr}} = 10$ wake core also starts showing appreciable IGW activity in the ambient by $x/D = 100$ ($Nt_{10}= 10$), as shown in figure 3(f).

5. Characteristics of SPOD eigenvalues and eigenspectra

We start the discussion of SPOD modes by evaluating their overall contribution to fluctuation energy and by their eigenspectra. There are significant effects of buoyancy as elaborated below.

5.1. Cumulative modal contribution to fluctuation energy

Figure 4 shows the variation of cumulative energy ($\xi (n)$) as a function of SPOD modal index ($n$) at four downstream locations: $x/D = 20$, 40, 80 and $100$. To calculate $\xi (n)$, the energy across all resolved frequencies $ {\textit {St}}$ at each modal index up to $n$ is summed and normalized by the total energy as

(5.1)\begin{equation} \xi(n;x/D) = \frac{\displaystyle\sum_{i=1}^{n}\sum_{St} \lambda^{(i)}(f;x/D)}{\displaystyle\sum_{i=1}^{N_{blk}}\sum_{St} \lambda^{(i)}(f;x/D)}, \end{equation}

where $N_{blk}$ is the total number of SPOD modes at a given $ {\textit {St}}$. Comparison among the various $x/D$ curves shows that the energy captured by leading SPOD modes in both wakes increases with downstream distance. This behaviour is in contrast to the unstratified wake where the relative importance of the dominant SPOD modes decreases with increasing $x/D$ (Nidhan et al. Reference Nidhan, Chongsiripinyo, Schmidt and Sarkar2020). Although both stratified wakes exhibit an increasing dominance of the leading modes as $x/D$ increases, there is a quantitative difference in that the jump of modal energy fraction from its $x/D = 20$ value is larger for the $ {\textit {Fr}} = 2$ wake relative to the $ {\textit {Fr}} = 10$ wake.

Figure 4. Variation of the cumulative fraction of energy, $\xi (n)$, as a function of modal index $n$ for (a) $ {\textit {Fr}}= 2$ and (b) $ {\textit {Fr}} = 10$ wakes, shown till $n = 15$ SPOD modes.

As discussed in the introduction, CS2020 found that the $ {\textit {Fr}} =2$ wake traversed the WST, IST and SST regimes during its streamwise evolution and the $ {\textit {Fr}} = 10$ wake accessed only the WST and IST regimes. Readers are referred to § 1 for an introduction to these regimes in stratified turbulence. These transitions also appear in the evolution of the modal energy content $\xi (n;x/D)$. For example, the $ {\textit {Fr}} = 2$ wake in figure 4(a) shows a transition at $x/D \approx 40$ whereby the $\xi (n)$ curves for $x/D \geqslant 40$ collapse onto a single profile. This result is consistent with CS2020 who find that $x/D \approx 40$ ($Nt_{2} \approx 20)$ is the location where the $ {\textit {Fr}} = 2$ wake transitions from IST to SST. The $ {\textit {Fr}} = 10$ wake was found by CS2020 to stay in the WST regime till $x/D \approx 50$ ($Nt \approx 5$) and thereafter transitioned to the IST regime. For the $ {\textit {Fr}} = 10$ wake in figure 4(b), the $\xi (n)$ curves collapse separately, i.e. there is one curve showing collapse between $x/D = 20$ and 40 which lies in the WST regime, and there is another showing collapse between $x/D = 80$ and $100$ which lies in the IST regime. Plots of $\xi (n)$ at other values of $x/D$ (not shown here) confirm that locations with $x/D \leqslant 50$ collapse on the $x/D = 20,40$ curve and locations with $x/D \geqslant 80$ collapse on the $x/D = 80,100$ curve.

The energy summed over frequencies instead of modal indices is now examined. Figure 5 shows the variation of $\xi ( {\textit {St}})$ calculated as

(5.2)\begin{equation} \xi( {\textit{St}}; x/D) = \frac{\displaystyle\sum_{f ={-}St}^{St}\sum_{i=1}^{N_{blk}} \lambda^{(i)}(f;x/D)}{\displaystyle\sum_{i=1}^{N_{blk}}\sum_{St} \lambda^{(i)}(f;x/D)}. \end{equation}

Figure 5 shows that $\xi ( {\textit {St}})$ increases for low-$ {\textit {St}}$ modes with increasing $x/D$ in both wakes, which is a trend also seen for $\xi (n)$. This is yet another indication of the increasing importance of the coherent modes as buoyancy effects come into play in these stratified wakes. Besides, for both wakes, $\xi ( {\textit {St}})$ increases steeply between $ {\textit {St}} = 0.1$ and $0.2$ at all downstream locations. The reason behind this sharp increase will be discussed shortly. Another observation of interest is that almost all the fluctuation energy at large $x/D$ is captured by the modes with $ {\textit {St}} < 1$ in both wakes.

Figure 5. Variation of the cumulative fraction of energy, $\xi ( {\textit {St}})$, as a function of $ {\textit {St}}$ for (a) $ {\textit {Fr}}= 2$ and (b) $ {\textit {Fr}} = 10$. The plots are shown for $0 \leqslant {\textit {St}} \leqslant 1$ for both cases. Inset plots show zoomed-in variation of $\xi ( {\textit {St}})$ for $0 \leqslant {\textit {St}} \leqslant 0.2$.

From $x/D = 20$ to $40$ there is a large jump in $\xi ( {\textit {St}})$ for the $ {\textit {Fr}} = 2$ wake in figure 5(a). As mentioned previously, $x/D = 40$ also marks the arrival of the $ {\textit {Fr}} = 2$ wake into the SST regime. Also, the $\xi ( {\textit {St}})$ curves collapse for locations $x/D = 80$ and $100$. On analysing other streamwise locations (not shown here), we find that the $\xi ( {\textit {St}})$ curves for $x/D \geqslant 70$ collapse together similar to the previously shown $\xi (n)$ curves of the $ {\textit {Fr}} = 2$ wake. One difference is that the collapse of $\xi (n)$ commences closer to the body at $x/D \approx 40$.

Contrary to the $ {\textit {Fr}} = 2$ wake where the change in $\xi ( {\textit {St}})$ from $x/D = 20$ to $x/D = 40$ was large, the corresponding change for the $ {\textit {Fr}} = 10$ wake (figure 5b) is small and consistent with an absence of regime change. However, the $ {\textit {Fr}} = 10$ wake exhibits a significant jump of $\xi ( {\textit {St}})$ between $x/D = 40$ and $80$, which lie in the WST and IST regime, respectively.

To summarize, figures 4 and 5 have the following implications. First, the relative importance of the leading SPOD modes increases with $x/D$ for the stratified wakes, which is in stark contrast to their behaviour in the unstratified wake (Nidhan et al. Reference Nidhan, Chongsiripinyo, Schmidt and Sarkar2020). Second, the trend of increasing dominance of leading SPOD modes is more pronounced for the strongly stratified wake of $ {\textit {Fr}} = 2$ as compared with $ {\textit {Fr}} = 10$. Third, transitions between WST, IST and SST regimes discussed by CS2020 for the turbulence statistics are also qualitatively reflected in the energetics of SPOD modes too.

5.2. Spectral POD eigenspectra of $ {\textit {Fr}} = 2$ and $10$ wakes

Figure 6 shows the SPOD eigenspectra of the $ {\textit {Fr}} = 2$ (a,c,e,g) and $ {\textit {Fr}} = 10$ (b,d,f,h) wakes at various downstream locations. The spectrum of the leading SPOD mode ($\lambda ^{(1)}$) shows a distinct spectral peak in the vicinity of $ {\textit {St}} = 0.13 - 0.15$ at all locations and for both wakes. This pronounced peak is the reason why there was a sharp increase of $\xi ( {\textit {St}})$ within $0.1 < {\textit {St}} < 0.2$ for both wakes in figure 5.

Figure 6. Spectral POD eigenspectra of the 25 most energetic modes, $\lambda ^{(1)}$ to $\lambda ^{(25)}$, for the $ {\textit {Fr}} = 2$ (a,c,e,g) and $ {\textit {Fr}} = 10$ (b,d,f,h) wakes at four streamwise locations: (a,b) $x/D = 20$, (c,d) $x/D = 40$, (e,f) $x/D =80$ and (g,h) $x/D = 100$. Dark to light shade corresponds to increasing model index $i$ in $\lambda ^{(i)}$.

In the $ {\textit {Fr}} = 10$ wake the $\lambda ^{(1)}$ eigenspectrum at all locations has a distinct peak at $ {\textit {St}} \approx 0.13$, which is very close to the VS frequency of the unstratified wake ($ {\textit {St}} = 0.135$) at the same Reynolds number (Nidhan et al. Reference Nidhan, Chongsiripinyo, Schmidt and Sarkar2020). Spectral POD eigenspectra at $x/D < 2$ (not presented here) show that this spectral peak has its origin near the wake generator and corresponds to VS in the $ {\textit {Fr}} = 10$ wake.

Unlike the $ {\textit {Fr}} = 10$ wake, the spectral peak in $\lambda ^{(1)}$ for the $ {\textit {Fr}} = 2$ wake shifts slightly from $ {\textit {St}} \approx 0.15$ at $x/D = 20$ to $ {\textit {St}} \approx 0.13$ at $x/D = 40$ and onward. At the far wake location of $x/D = 80$, the peak in the $\lambda ^{(1)}$ eigenspectrum broadens to reach $ {\textit {St}} \approx 0.11$. Near-body SPOD eigenspectra (not shown here) for the $ {\textit {Fr}} = 2$ wake show a prominent peak at $ {\textit {St}} \approx 0.15$ (slightly larger relative to the unstratified and $ {\textit {Fr}} = 10$ wakes) just downstream of the recirculation zone at $x/D \approx 2$. Furthermore, the pressure spectrum (not presented here) in the immediate proximity of the disk, at $x/D = 0.5$ and $r/D=0.5$, also peaks at $ {\textit {St}} \approx 0.15$, indicating that this frequency corresponds to the VS mechanism for the $ {\textit {Fr}} = 2$ wake. The shift in the spectral peak towards lower $ {\textit {St}}$ at later $x/D$ is consistent with the sphere-wake study of Spedding (Reference Spedding2002a) who report a gradual reduction in the dominant wake $ {\textit {St}}$ during $40 < Nt < 100$ (see figure 5 of their paper).

In the $ {\textit {Fr}} = 2$ wake there is a large gap (demarcated in red) between the $\lambda ^{(1)}$ and $\lambda ^{(2)}$ spectra for frequencies with $ {\textit {St}} < 0.2$. Beyond $ {\textit {St}} \approx 0.2$, values of all $\lambda ^{(i)}$ fall sharply. This large difference between $\lambda ^{(1)}$ and $\lambda ^{(2)}$ implies that the dynamics of the $ {\textit {Fr}}=2$ wake is low rank, i.e. it is dominated by the leading SPOD mode. The sharp drop-off in energy at higher $ {\textit {St}}$ points to the dominance of low-frequency energetic structures with $ {\textit {St}}$ in $[0, 0.2]$, specifically around the VS frequency.

In terms of low-rank behaviour, the $ {\textit {Fr}} = 10$ wake shows a peculiar difference from the $ {\textit {Fr}} = 2$ wake. Although the gap between $\lambda ^{(1)}$ and $\lambda ^{(2)}$ is significantly less compared with that for $ {\textit {Fr}} =2$, there is a significant gap between $\lambda ^{(2)}$ and $\lambda ^{(3)}$ around the VS frequency $ {\textit {St}} \approx 0.13$, shown in red in the right column of figure 6. Furthermore, the variation of $\lambda ^{(2)}$ with $ {\textit {St}}$ is very similar to that of $\lambda ^{(1)}$. On further investigation, we find that the SPOD eigenmodes of $\lambda ^{(1)}$ and $\lambda ^{(2)}$ at the VS frequency have a similar spatial structure, but with a rotation in their orientation. We hypothesize that $\lambda ^{(1)}$ and $\lambda ^{(2)}$ modes at the VS frequency are the manifestation of $m=1$ and $m=-1$ azimuthal modes in the weakly stratified $ {\textit {Fr}} = 10$ wake.

Figure 7 shows the fraction of energy in each SPOD mode as a function of $ {\textit {St}}$ for both wakes at a representative location of $x/D = 50$. In both wakes, the leading SPOD mode at the VS frequency capture at least $40\,\%$ of the total energy contained in the VS frequency. This also holds true for the near ($x/D = 10$) and far ($x/D = 100$) locations in both wakes (not discussed here for brevity). Also, in both wakes, less than five SPOD modes are required to capture $75\,\%$ of the total energy in the vicinity of the VS frequency, as indicated by the solid white line in figure 7.

Figure 7. Fraction of energy (at a given $ {\textit {St}}$) accounted by each SPOD mode as a function of frequency at $x/D = 50$ for (a) $ {\textit {Fr}} = 2$ and (b) $ {\textit {Fr}} = 10$ wakes. The solid and dashed white lines indicate the number of SPOD modes required to retain 75 % and 50 % of the total fluctuation energy respectively at each frequency.

The key takeaway from figures 6 and 7 is twofold: (i) the VS frequency is the leading contributor to the fluctuating energy content of both $ {\textit {Fr}} = 2$ and $10$ wakes and (ii) its dynamics is primarily governed by a few leading SPOD modes. Previous experimental studies of Chomaz et al. (Reference Chomaz, Bonneton and Hopfinger1993) and Lin et al. (Reference Lin, Boyer and Fernando1992a) have showed the existence of the VS mode in the near wake at moderate stratification using hot-wire measurements (at few select locations) and shadowgraph techniques. The present SPOD analysis enables us to establish the dominance of the VS mode in stratified wakes from near the body to 100 body diameters downstream by providing an ordered set of $\lambda ^{(i)}$ eigenvalues for different $ {\textit {St}}$.

6. The energetics of the VS mode

A comparison between SPOD eigenspectra of the stratified wakes (figure 6) and the unstratified wake (figure 8 in Nidhan et al. Reference Nidhan, Chongsiripinyo, Schmidt and Sarkar2020) reveals that both types of wakes are dominated by VS which gives rise to a distinct spectral peak in the vicinity of $ {\textit {St}} \approx 0.13$. For the unstratified wake, besides the VS structure, which appears in the azimuthal mode $m=1$, a double helix ($m=2$) mode with a peak at $ {\textit {St}} \rightarrow 0$ is also found to be energetically important (Johansson & George Reference Johansson and George2006; Nidhan et al. Reference Nidhan, Chongsiripinyo, Schmidt and Sarkar2020). In the stratified wake, as elaborated below, we find that the VS mode is persistent, is linked to unsteady IGWs, and is thereby responsible for the accumulation of fluctuation energy outside the wake core.

Stratification qualitatively affects the streamwise evolution of the energy in different frequencies. The evolution of the frequency-binned energy is shown for the stratified wakes in figure 8(a,b). For the unstratified case ($ {\textit {Fr}} = \infty$), the azimuthal modes $m=1$ and $m=2$ are shown in figures 8(c) and 8(d), respectively. For the stratified wakes in figure 8(a,b), the spectral peak in the vicinity of $ {\textit {St}} \approx 0.13$ remains prominent for significant downstream distances, especially for $ {\textit {Fr}} = 2$. A somewhat wide band ($0.1 \leqslant {\textit {St}} \leqslant 0.2$), centred around $ {\textit {St}} \approx 0.13$ of the VS mode, is excited for the stratified wakes. Furthermore, this band persists into the far wake. Even at $x/D = 100$ this band has a larger energy density than at other frequencies. Such persistence in the energetic dominance of the VS mode (and neighbouring frequencies) is absent in the unstratified $ {\textit {Fr}} = \infty$ case where the energy at the two peaks of (i) $ {\textit {St}} = 0.135$ in the $m=1$ mode (figure 8c) and (ii) $ {\textit {St}} = 0$ in the $m=2$ mode (figure 8d) declines sharply with increasing $x/D$.

Figure 8. Contour maps, $x/D- {\textit {St}}$, showing the variation of total energy in the leading 15 SPOD modes: (a) $ {\textit {Fr}} = 2$, (b) $ {\textit {Fr}} =10$, (c) $ {\textit {Fr}} = \infty$, $m=1$ (vortex shedding) mode, and (d) $ {\textit {Fr}} = \infty$, $m=2$ (double helix) mode. The contour level of each plot is set between zero and the maximum value of the energy over all $(x/D, {\textit {St}})$ pairs in that plot. Results for $ {\textit {Fr}} = \infty$ wake are taken from Nidhan et al. (Reference Nidhan, Chongsiripinyo, Schmidt and Sarkar2020).

Figure 9 shows the streamwise evolution of energy in the leading 15 SPOD modes and in a frequency band around the VS frequency. The energy in the $ {\textit {Fr}} = 2$ wake remains almost constant till $x/D = 60$ and starts decaying slowly thereafter. On the other hand, the $ {\textit {Fr}} = 10$ wake shows an initial decay in the VS mode energy which closely follows that of the $ {\textit {Fr}} = \infty$ wake till $x/D = 20$. Subsequently, buoyancy effects set in for the $ {\textit {Fr}} = 10$ wake to slow down the energy decay.

Figure 9. Evolution of the energy contained in the leading 15 SPOD modes at the VS frequency is shown for $ {\textit {Fr}} = 2, 10$ and $\infty$ wakes. The energy is normalized by its value at $x/D = 10$. Energy across $0.11 \leqslant {\textit {St}} \leqslant 0.15$ is summed for the $ {\textit {Fr}} = 2$ wake since it has a broader spectral peak (see figure 6).

To investigate the reason behind the downstream persistence of the VS spectral peak in stratified wakes, the total energy in the leading 15 SPOD modes is partitioned into two components: (i) energy of the wake core, $E_{core}$, and (ii) energy of the outer wake, $E_{outer}$. The energy in each of the regions is calculated as

(6.1)\begin{gather} E_{core}(x/D, {\textit{St}}) = \sum_{n=1}^{15} \,\int_{A \in \varOmega} \lambda^{(n)}(x/D, {\textit{St}})\varPhi^{*(n)}_{i}(x/D, {\textit{St}};y,z)\varPhi^{(n)}_{i}(x/D, {\textit{St}};y,z)\,\mathrm{d} A, \end{gather}

(6.2)\begin{gather}E_{outer}(x/D, {\textit{St}}) = \sum_{n=1}^{15} \, \int_{A \in \mathcal{H} - \varOmega} \lambda^{(n)}(x/D, {\textit{St}}) \varPhi^{*(n)}_{i}(x/D, {\textit{St}};y,z)\varPhi^{(n)}_{i}(x/D, {\textit{St}};y,z)\,\mathrm{d} A, \end{gather}

where $\varOmega$ denotes the wake core at a given $x/D$, as defined in § 4 and ${\rm d}A = {\rm d}y{\rm d}z$. Here $\mathcal {H}$ denotes the area of the circular cross-section bounded by $0 \leqslant r/D \leqslant 10$ at a given $x/D$, and $\varPhi ^{(n)}_{i}$ corresponds to the $n{\rm th}$ SPOD eigenmode for a given $x/D$ and $ {\textit {St}}$.

The energy in the wake core peaks around the VS mode frequency, $ {\textit {St}} \approx 0.12 - 0.13$, for both wakes (see figure 10a,c). With increasing $x/D$ (or $Nt$), the VS signature in the wake core decays for both wakes. The energetics of the outer wake is remarkably different; $E_{outer}$, which starts off with a small value across all $ {\textit {St}}$ at $x/D \approx 10$ in the $ {\textit {Fr}} = 2$ wake, develops a peak at $ {\textit {St}} \approx 0.15$ at $x/D \approx 20$. Note that this peak is the same as the peak in the SPOD eigenspectrum for the entire wake (figure 6a). Farther downstream, there is significant energy content in the outer wake for $x/D \approx 16 - 80$ ($Nt_{2} \approx 8 - 40$) with a spectral peak located at $ {\textit {St}} \approx 0.13 - 0.15$. The spectral peak is broad, i.e. nearby frequencies with $0.1 \leqslant {\textit {St}} \leqslant 0.2$ also have comparable energy levels. For the $ {\textit {Fr}} = 10$ wake, $E_{outer}$ picks up only beyond $x/D = 60$ ($Nt_{10} = 6$), and thereafter increases progressively in the vicinity of $ {\textit {St}} \approx 0.13$ till the end of the domain. We also find that the qualitative nature of the variation of energy in the outer wake and wake core found in figure 10 does not change when the number of modes over which energy is summed is decreased from 15 to 3 (not presented here for brevity).

Figure 10. Energy partition between wake core and outer wake: (a) wake core of $ {\textit {Fr}}=2$, (b) outer wake of $ {\textit {Fr}} =2$, (c) wake core of $ {\textit {Fr}} =10$ and (d) outer wake of $ {\textit {Fr}} = 10$. The first 15 SPOD modes are considered and the contour level of each plot is set between zero and the maximum value of the energy over all $(x/D, {\textit {St}})$ pairs in that plot.

In figure 11 we sum up the SPOD energies across $ {\textit {St}} \in [-0.4, 0.4]$ separately for the wake core and the outer wake and compute their percentage contribution to the entire area-integrated fluctuation energy as

(6.3)\begin{gather} \chi_{core}(x/D) = \frac{\displaystyle\sum_{|St|=[0, 0.4]}E_{core}( {\textit{St}},x/D)}{E_{k}^{T}(x/D) + E_{\rho}^{T}(x/D)} \times 100, \end{gather}

(6.4)\begin{gather}\chi_{outer}(x/D) = \frac{\displaystyle\sum_{|St|=[0, 0.4]}E_{outer}( {\textit{St}},x/D)}{E_{k}^{T}(x/D) + E_{\rho}^{T}(x/D)} \times 100, \end{gather}

where $E^{T}_{k}(x/D)$ and $E^{T}_{\rho }(x/D)$ are area-integrated TKE and TPE in the circular region of $0 \leqslant r/D \leqslant 10$ at the $x/D$ location under consideration. The streamwise evolution of $\chi _{core}$ and $\chi _{outer}$ are shown in figure 11.

Figure 11. Streamwise variation of $\chi _{core}$, $\chi _{outer}$ and $\chi _{core} + \chi _{outer}$: (a) $ {\textit {Fr}} = 2$, (b) $ {\textit {Fr}} =10$.

For the $ {\textit {Fr}} = 2$ wake (figure 11a), $\chi _{outer}$ increases monotonically until $x/D \approx 60$ followed by a slight decrease. At its peak, $\chi _{outer}$ constitutes up to $50\,\%$ of the total fluctuation energy, becoming even larger than $\chi _{core}$. In the $ {\textit {Fr}} = 10$ wake (figure 11b), $\chi _{outer}$ remains negligible till $x/D = 60$, followed by a monotonic increase. The increase in the value of $\chi _{outer}$ is accompanied by a decrease in the wake core contribution. The percentage of total energy captured by the leading 15 SPOD modes and $| {\textit {St}}| \in [0, 0.4]$, i.e. $\chi _{wake} + \chi _{outer}$ (shown in green), increases for both wakes from its initial value at $x/D = 10$. This reinforces a main finding of this work that stratified wakes display an increased coherence as they evolve downstream.

Figure 10 suggests that the unsteady IGWs in the outer wake radiate from the VS mode at intermediate to late $Nt$. Nevertheless, to further establish causation between the unsteady IGW emission and the VS mode, we perform additional SPOD analyses for the $ {\textit {Fr}} = 2$ wake. In these analyses, we replace the fluctuating density field $\rho '$ with the fluctuating pressure field $p'$. These SPOD analyses are performed at $x/D = 10, 20, \ldots, 90,100$. At all locations, the eigenspectra obtained from these modified SPOD analyses show a prominent peak at the VS frequency with a large gap between $\lambda ^{(1)}$ and $\lambda ^{(2)}$ for $ {\textit {St}} < 0.2$, qualitatively akin to the left column of figure 6.

Using $p'$ along with $u'_i$ enables us to reconstruct the pressure transport term in the radial direction, $\langle p'u'_r \rangle$, which accounts for the energy transferred radially from the wake core to the IGW dominated outer wake region through pressure work (de Stadler & Sarkar Reference de Stadler and Sarkar2012; Rowe et al. Reference Rowe, Diamessis and Zhou2020). We reconstruct $\langle p'u'_r \rangle$ contours using the leading 15 SPOD modes and frequencies in the range of (i) $ {\textit {St}} \in [0.1,0.2]$ and (ii) $ {\textit {St}} \in [0.1,0.3]$ as follows:

(6.5)\begin{equation} \langle p'u'_r \rangle(x;y,z) = \sum_{St }\sum_{n=1}^{15}\lambda^{(n)}(x; {\textit{St}}) \, \varPhi_{u_r}^{(n)}(x;y,z, {\textit{St}}) \, \varPhi_{p}^{(n)*}(x;y,z, {\textit{St}}), \end{equation}

Figure 12 shows the actual (a–c) and reconstructed (d–i) $\langle p'u'_r \rangle$ at three streamwise locations $x/D = 20, 40$ and $60$ for the $ {\textit {Fr}} = 2$ wake. The actual $\langle p'u'_r\rangle$ shows a strong signature of IGW flux in the outer wake region at all three downstream locations in figure 12. We found that the nonlinear transport term was negligible outside the wake core (not shown here). Hence, the primary source of the energy transfer to the outer wake is the pressure-work term due to the IGW radiation. The reconstructed $\langle p'u'_r\rangle$ using $ {\textit {St}} \in [0.1, 0.2]$ (middle row) shows qualitative agreement with the spatial distribution of actual $\langle p'u'_r\rangle$, both in the wake core as well as outer wake region, at all downstream locations. As more frequencies are included (bottom row), adjacent to the VS frequency, the accuracy of reconstruction increases. The key spatial characteristics of $\langle p'u'_r\rangle$ remain similar in both reconstructions, showing that frequencies in the vicinity of the VS frequency satisfactorily capture the key dynamics of unsteady IGW generation.

Figure 12. Contours of $\langle p'u'_r \rangle$ for the $ {\textit {Fr}} = 2$ wake obtained from (i) temporal averaging (a–c), reconstruction using the leading 15 SPOD modes and (ii) $ {\textit {St}} \in [0.1, 0.2]$ (d–f) and (iii) $ {\textit {St}} \in [0.1,0.3]$ (g–i). Three streamwise locations $x/D = 20$, 40 and $60$ (left to right column) are shown. Dashed closed curve in white shows wake core.

To further elucidate the causal link between the VS mode and the IGW generation, we plot the $x/D- {\textit {St}}$ variation of the integrated $\langle p'u'_r\rangle$ in the outer wake region (similar to (6.2)), reconstructed using 3 and 15 SPOD modes (figure 13). The evolution of integrated $\langle p'u'_r\rangle$ in the outer wake of the $ {\textit {Fr}} = 2$ wake is very similar to that of the energy in the outer wake region (figure 10b). It starts off at a small value at $x/D \approx 10$, develops a broad peak centred at $ {\textit {St}} \approx 0.15$ between $20 \leqslant x/D \leqslant 80$, and gradually starts declining beyond $x/D = 80$. Increasing the number of modes from $3$ to $15$ makes the active $ {\textit {St}}$ region broader while intensifying the reconstructed values.

Figure 13. Contour maps, $x/D- {\textit {St}}$, showing the variation of integrated $\langle p'u'_r\rangle$ in the outer wake region using (a) the leading 3 SPOD modes and (b) the leading 15 SPOD modes for the $ {\textit {Fr}} = 2$ wake.

Figures 10, 12 and 13 firmly establish that the VS mode energy radiates out of the wake core instead of being acted on by nonlinear interactions in the turbulent wake responsible for the usual energy cascade. Therefore, unlike their unstratified counterpart, the stratified wakes exhibit a persistent VS spectral peak when the energy in the full domain of influence (denoted by $\mathcal {H}$) of the wake is taken into account as in the SPOD results of figure 8(a,b).

7. Spatial structure of SPOD eigenmodes

7.1. Spatial structure of the VS eigenmode

The spatial structure of the dominant eigenmodes sheds further light on the manner in which buoyancy helps spread unsteady flow perturbation to well outside the turbulent core of the wake. Figure 14 shows the real part of the normalized (by $L_{\infty }$ norm) leading SPOD mode, $\varPhi ^{(1)}_{y}(y,z, {\textit {St}};x)/||\varPhi ^{(1)}_{y}(y,z, {\textit {St}};x)||_{\infty }$, of the lateral velocity $u_y$. The plotted modes correspond to the VS mechanism, which is at $ {\textit {St}}$ corresponding to the eigenspectrum peak and are shown for selected values of $x/D$. The ellipsoid wake core (dashed white curve) with dimensions $2L_{Hk}$ and $2L_{Vk}$ is also shown. At $ {\textit {Fr}} =2$ (upper row), the wake core exhibits flattening from $x/D = 10$ ($Nt =5$) onward and the eigenmodes in the core show horizontal layering for $x/D \geqslant 50$. The layering becomes visible in the eigenmodes at $x/D \approx 30$ (not shown here). The region outside the core has little activity at $x/D =10$ but shows IGW phase lines at $x/D = 50$ and 100. There is a clear and continuous transition of the eigenmode from its layered core to an IGW structure in the outer region at the far downstream locations. The flattening of the wake core and the IGW related spread of the eigenmode is delayed for the $ {\textit {Fr}} = 10$ wake (bottom row) relative to $ {\textit {Fr}} = 2$ since equivalent $Nt$ values occur farther downstream.

Figure 14. Shape of the leading SPOD mode (real part corresponding to $\lambda ^{(1)}$) for spanwise velocity, $\varPhi ^{(1)}_{y}(y,z, {\textit {St}};x/D)$: (a) $x/D = 10, {\textit {Fr}} = 2, {\textit {St}} = 0.15$, (b) $x/D = 50, {\textit {Fr}} = 2, {\textit {St}} = 0.13$, (c) $x/D = 100, {\textit {Fr}} = 2, {\textit {St}} = 0.13$, (d) $x/D = 10, {\textit {Fr}} = 10, {\textit {St}} = 0.13$, (e) $x/D = 50, {\textit {Fr}} = 10, {\textit {St}} = 0.13$ and (f) $x/D = 100, {\textit {Fr}} = 10, {\textit {St}} = 0.13$. At each $x/D$, the shown mode corresponds to the peak in the eigenspectrum of $\lambda ^{(1)}$. Real part of each mode is shown. Dashed closed curve in white shows wake core.

Comparing figure 14(b,c) with figure 3(b,c), there are striking similarities in the layered structure of the $ {\textit {Fr}} = 2$ wake core between the dominant eigenmodes and the instantaneous snapshots at the far downstream locations of $x/D = 50$ and $100$. Although SPOD only guarantees that the obtained modes optimally capture the prescribed energy norm of the flow (see § 3.1), these modes do generally contain the imprints of actual flow structures, as is the case here. The outer wake shows that distinct IGWs are associated with the wake core structure of dominant eigenmodes at late $Nt$ for both $ {\textit {Fr}} = 2$ and $10$ wakes. For the $ {\textit {Fr}} = 2$ wake, IGW activity in the outer region of the eigenmodes shown in figure 14 is negligible at $x/D = 10$ ($Nt_{2} = 5$) while it is readily noticeable at $x/D = 50$ ($Nt_{2} = 25$) and $x/D = 100$ ($Nt_{2} = 50$). The IGWs are found to be emitted within $30^\circ - 60^\circ$ with the $y$-axis. For $ {\textit {Fr}} = 10$, the IGWs found at $x/D = 100$ ($Nt_{10} = 10$) are emitted at ${\approx }45^\circ$ from the horizontal. A comparison between figures 3 and 14 reveals that the IGW in the dominant eigenmodes (figure 14) represent the IGWs in actual snapshots (figure 3) to a satisfactory extent, emphasizing that the VS mechanism is an important IGW generation mechanism in stratified wakes. The leading VS modes at different locations show asymmetry about the $y=0$ line, in both $ {\textit {Fr}} = 2$ (at $x/D = 10$) and $ {\textit {Fr}} = 10$ (at $x/D = 50$ and $100$) wakes. This could be a consequence of the presence of a very-low-frequency mode in the wake (Grandemange, Gohlke & Cadot Reference Grandemange, Gohlke and Cadot2013; Rigas et al. Reference Rigas, Oxlade, Morgans and Morrison2014).

To analyse the streamwise coherence of the leading SPOD eigenmode at the VS frequency, we conduct an additional SPOD analysis for the $ {\textit {Fr}} = 2$ wake, using the fluctuating density and velocity fields, at the centre-vertical plane ($y=0$). Specific details of this SPOD analysis are mentioned towards the end of § 3.2. The SPOD eigenspectrum (not shown here for brevity) shows a broad peak at $ {\textit {St}} \approx 0.13$ .

Figure 15 shows the spatial structure of the spanwise component of the leading SPOD mode at $ {\textit {St}} \approx 0.13$. The VS mode appears to be strongly coherent in the streamwise direction with a wavelength of $\lambda /D \approx 1/St_{VS}$. It has two distinctive features: (i) emergence of a well-defined IGW signature beyond $x/D \approx 20$ and (ii) gradual transition of the opposite signed lobes into V-shaped structures as the wake progresses downstream. These structures get progressively thinner and shallower (with respect to the $x$-axis) as $x/D$ increases.

Figure 15. Shape of the leading SPOD mode for spanwise velocity, $\varPhi ^{(1)}_{y}(x, z, {\textit {St}} \approx 0.13; y=0)$ in the centre-vertical plane for the $ {\textit {Fr}} = 2$ wake. Real part of the mode is shown in domain $z, x \in [-10, 10] \times [2, 100]$.

7.2. Spatial structure of high-$ {\textit {St}}$, high-$n$ eigenmodes

To contrast the structure of less energetic SPOD eigenmodes with the dominant SPOD eigenmodes, the $u_{y}$ eigenmode at $n=15$ and $ {\textit {St}} = 0.40$ is plotted in figure 16 at the same downstream locations of $x/D = 10, 50$ and $100$ considered previously. It should be noted that these SPOD modes have low energy, $O(10^{-2})$ that of the dominant SPOD modes. Visual inspection shows that the spatial coherence in the wake core, which is a characteristic of dominant SPOD modes, is lost for the high-$n$ and high-$ {\textit {St}}$ modes similar to the result in the snapshot POD study of Diamessis et al. (Reference Diamessis, Gurka and Liberzon2010). For both $ {\textit {Fr}} = 2$ and $10$ wakes, $\varPhi ^{(15)}_{y}(y,z, {\textit {St}} = 0.40)$ in the wake core is dominated by small-scale turbulence. For the $ {\textit {Fr}} = 2$ wake, the distinct layered structure found in the leading VS eigenmodes at $x/D = 50$ and $100$ is absent in the low-energy mode at the same locations. Nevertheless, buoyancy-induced anisotropy is evident at $Nt \geqslant 5$ in both wakes even in these low-energy modes with high $n$ and $ {\textit {St}}$. Moreover, the $\varPhi ^{(15)}_y(y,z, {\textit {St}} \approx 0.40)$ mode also shows IGWs in the outer wake at $x/D = 50$ and $100$ in the $ {\textit {Fr}} = 2$ wake (figure 16b,c), albeit with a smaller wavelength than for the VS mode. Contrary to the $ {\textit {Fr}} = 2$ wake, the $\varPhi ^{(15)}_{y}(y,z, {\textit {St}} = 0.4)$ mode for the $ {\textit {Fr}} = 10$ wake does not show any IGW in figure 16(e,f).

Figure 16. Shape of the fifteenth SPOD mode (real part corresponding to $\lambda ^{(15)}$) at $ {\textit {St}} = 0.40$ for spanwise velocity, $\varPhi ^{(15)}_{y}(y,z, {\textit {St}};x/D)$: (a) $x/D = 10, {\textit {Fr}} = 2$, (b) $x/D = 50, {\textit {Fr}} \!=\! 2$, (c) $x/D \!=\! 100, {\textit {Fr}} \!=\! 2$, (d) $x/D \!=\! 10, {\textit {Fr}} = 10$, (e) $x/D = 50, {\textit {Fr}} = 10$ and (f) $x/D = 100, {\textit {Fr}} = 10$. Dashed closed curve in white shows the wake core.

8. Reconstruction using SPOD modes

In this section we demonstrate the effectiveness of SPOD modes in reconstructing the following turbulence statistics: (i) TKE, $\langle u'_{i}u'_{i} \rangle /2$, (ii) lateral production $\mathcal {P}_{xy} = \langle -u'_{x}u'_{y}\rangle \partial \langle U \rangle /\partial y$ and (iii) buoyancy flux $\mathcal {B} = \langle -\rho 'u'_{z} \rangle / {\textit {Fr}}^{2}$. The reconstruction from SPOD modes is performed as follows:

(8.1)\begin{gather} \mathrm{TKE}(x;y,z) = \frac{1}{2}\sum_{n=1}^{\varLambda}\, \sum_{ {\textit{St}} ={-} {\textit{St}}_r}^{ {\textit{St}} = {\textit{St}}_r} \lambda^{(n)}(x; {\textit{St}}) \, \varPhi_{i}^{(n)}(x;y,z, {\textit{St}}) \, \varPhi_{i}^{(n)*}(x;y,z, {\textit{St}}), \end{gather}

(8.2)\begin{gather}\mathcal{P}_{xy}(x;y,z) = \sum_{n=1}^{\varLambda}\, \sum_{ {\textit{St}} ={-} {\textit{St}}_r}^{ {\textit{St}} = {\textit{St}}_r} -\lambda^{(n)}(x; {\textit{St}}) \, \varPhi_{x}^{(n)}(x;y,z, {\textit{St}}) \, \varPhi_{y}^{(n)*}(x;y,z, {\textit{St}}) \frac{\partial \langle U \rangle}{\partial y}, \end{gather}

(8.3)\begin{gather}\mathcal{B}(x;y,z) = \frac{1}{ {\textit{Fr}}^2}\sum_{n=1}^{\varLambda}\, \sum_{ {\textit{St}} ={-} {\textit{St}}_r}^{ {\textit{St}} = {\textit{St}}_r} -\lambda^{(n)}(x; {\textit{St}}) \, \varPhi_{\rho}^{(n)}(x;y,z, {\textit{St}}) \, \varPhi_{z}^{(n)*}(x;y,z, {\textit{St}}). \end{gather}

In (8.1)–(8.3) the values of $\varLambda$ and $ {\textit {St}}_r$ determine the set of modes used for reconstruction. The so-obtained turbulence statistics vary spatially in spanwise and vertical directions for different $x/D$.

Throughout this section, two sets of low-order truncation are used for reconstruction: (i) $n \leqslant 5$, $| {\textit {St}}| \leqslant 0.20$ (R1) and (ii) $n \leqslant 15$, $| {\textit {St}}| \leqslant 0.40$ (R2). While the R1 truncation primarily takes the VS mode into account for both wakes, R2 also accounts for some of the low-energy modes which reside at relatively higher $n$ and $ {\textit {St}}$. It should be noted that R1 and R2 set of modes account for approximately 0.7 % and 4.34 % of the total SPOD modes in both wakes.

Figure 17 compares the reconstructed TKE with its actual value for the $ {\textit {Fr}} =2$ wake at $x/D = 20$, $50$ and $100$. The actual TKE decays in magnitude, expands horizontally and narrows vertically with increasing $x/D$. At $x/D = 50$ and 100, the TKE contours display horizontal layering. At all three locations, reconstruction using the R1 set of modes (b,e,h) gives a fairly accurate estimate of the shape and spatial extent of the TKE contour. The layering at $x/D = 50$ and $100$ is also captured by the R1 reconstruction. These layers were also present in the reconstruction using only $n=1$ and $|St| \leqslant 0.2$ modes (not shown here), indicating the low-rank nature of layering in stratified wakes. On further increasing $[n, {\textit {St}}]$ as in the R2 reconstruction (c,f,i), the overall shape and structural features of the reconstructed TKE remain unchanged, while the magnitude increases, particularly at intense TKE locations, increasing the overall accuracy. It can also be ascertained visually that the accuracy of R1 and R2 increases with downstream distance pointing to the increasing coherence of the wake as it progresses downstream.

Figure 17. Contours of TKE for the $ {\textit {Fr}} = 2$ wake obtained from temporal averaging (a,d,g), reconstruction from R1 set of modes (b,e,h), and reconstruction from the R2 set of modes (c,f,i). Three streamwise locations $x/D = 20$, $50$ and $100$ are shown.

Figure 18 pertains to the reconstruction of the lateral production $\mathcal {P}_{xy}$ in the $ {\textit {Fr}} = 2$ wake. We limit ourselves to the lateral component since it dominates its vertical counterpart after the onset of buoyancy-induced suppression of vertical turbulent motions (Brucker & Sarkar Reference Brucker and Sarkar2010; de Stadler & Sarkar Reference de Stadler and Sarkar2012; Redford et al. Reference Redford, Lund and Coleman2015). In the IST and SST regimes of the disk wake, $\mathcal{P}_{xy}$ is the dominant component of turbulent production. The actual $\mathcal {P}_{xy}$ (a,d,g) shows two off-axis lobes of intense production primarily located near the horizontal centre plane ($z/D = 0$). With increasing $x/D$, these lobes flatten owing to buoyancy. With respect to the lateral production, the R1 and R2 set of modes capture the spatial distribution accurately for the $ {\textit {Fr}} = 2$ wake as shown in the middle and right columns of figure 18, respectively. Although SPOD modes are optimal for capturing the area-integrated sum of $\langle u'_iu'_i\rangle$ and $\langle \rho '\rho '\rangle$/$ {\textit {Fr}}^{2}$ by construction, we find that these modes provide an excellent low-order approximation for the production too.

Figure 18. Contours of $\mathcal {P}_{xy}$ for the $ {\textit {Fr}} = 2$ wake obtained from temporal averaging (a,d,g), reconstruction from R1 set of modes (b,e,h), and reconstruction from R2 set of modes (c,f,i). Three streamwise locations $x/D = 20$, 50 and $100$ are shown.

Finally, we explore the effectiveness of buoyancy flux ($\mathcal {B}$) reconstruction in figure 19. Unlike TKE and $\mathcal {P}_{xy}$, $\mathcal {B}$ is not a same signed quantity in the turbulent wake, as can be seen from figure 19(a,d,g). The R1 reconstruction of $\mathcal {B}$ (middle column) accurately captures the structural features of $\mathcal {B}$ at all locations: (i) layers of positively and negatively signed $\mathcal {B}$ at $x/D = 50$, $100$; and (ii) IGWs in the outer wake which carry significant $\mathcal {B}$ at $x/D = 20$ and $50$. On closer inspection, the R1 truncation is found to underpredict the strength of $\mathcal {B}$ in these outer regions with intense buoyancy flux. Including higher $ {\textit {St}}$ and $n$ modes for reconstruction, as done for R2, significantly improves the quality as shown in the right column.

Figure 19. Contours of $\mathcal {B}$ for the $ {\textit {Fr}} = 2$ wake obtained from temporal averaging (a,d,g), reconstruction from R1 set of modes (b,e,h), and reconstruction from R2 set of modes (c,f,i). Three streamwise locations $x/D = 20$, $50$ and $100$ are shown.

The reconstruction trends of these statistical quantities are also investigated for the $ {\textit {Fr}} = 10$ wake, but are not shown here for brevity. Qualitatively, the trends are similar to that of the $ {\textit {Fr}} = 2$ wake, wherein the R1 set captures the structural features of these quantities very satisfactorily. Further addition of high-$n$ and high-$ {\textit {St}}$ modes in the R2 truncation improves the quantitative prediction of these statistics, particularly in the region where they are found to be intense in the actual data.

To conclude this section, the streamwise variations of the wake-core-integrated TKE and $\mathcal {P}_{xy}$ for the $ {\textit {Fr}} = 2$ and $ {\textit {Fr}} = 10$ wakes are shown in figure 20. The corresponding variation for $\mathcal {B}$ is not shown here as it fluctuates between small positive and negative values, unlike TKE and $\mathcal {P}_{xy}$ which decay monotonically with $x/D$.

Figure 20. Streamwise variation of wake core TKE and $\mathcal {P}_{xy}$ reconstructed from R1 and R2 truncations: (a) TKE for $ {\textit {Fr}} = 2$, (b) TKE for $ {\textit {Fr}} = 10$, (c) $\mathcal {P}_{xy}$ for $ {\textit {Fr}} =2$ and (d) $\mathcal {P}_{xy}$ for $ {\textit {Fr}} = 10$. Here, $\int _{C} \ (.) \,\mathrm {d} C$ denotes the integration in the wake core.

For the $ {\textit {Fr}} = 2$ wake, wake core TKE shows two distinct decay rates: (i) $\mathrm {TKE} \propto x^{-1.11}$ in the IST regime spanning $10 \leqslant x/D \leqslant 40$ ($5 \leqslant Nt_2 \leqslant 20$) and (ii) $\mathrm {TKE} \propto x^{-0.30}$ in the SST regime spanning $60 \leqslant x/D \leqslant 100$ ($30 \leqslant Nt_{2} \leqslant 50$). The quality of TKE reconstruction improves monotonically from R1 to R2 at all downstream locations for both wakes, as is observed in figure 20(a,b). For $ {\textit {Fr}} = 2$, the TKE contained in the R1 set of modes stays approximately constant for $x/D \geqslant 40$ (figure 20a). It is only after high-$ {\textit {St}}$ and high-$n$ modes are added, as in R2, that the reconstructed TKE follows the decay rate of actual TKE. Reconstructed TKE from further lower-order truncations (not shown here), i.e. with lesser $n$ and $ {\textit {St}}$ than in R1, showed an increase in wake core TKE at large $x/D$, opposite to the decrease in the actual value. For the $ {\textit {Fr}} = 10$ wake, reconstruction from low-order truncations decay quite similar to the actual TKE (figure 20b).

The accuracy of R1 and R2 increase approximately three fold and two fold from $x/D = 10$ to $100$ for the TKE reconstruction in the $ {\textit {Fr}} = 2$ wake, suggesting development of low-rank dynamics in the $ {\textit {Fr}} = 2$ wake. By $x/D = 100$, R1 and R2 modes capture ${\approx } 63\,\%$ and ${\approx }82\,\%$, respectively, of the wake core TKE for $ {\textit {Fr}} = 2$. On the other hand, the reconstruction quality of the moderately stratified $ {\textit {Fr}} = 10$ wake changes only slightly from $x/D = 10$ to $100$ for both low-order truncations: (i) TKE in R1 modes changes from ${\approx } 29\,\%$ of total TKE at $x/D = 10$ to ${\approx } 38\,\%$ at $x/D = 100$ and (ii) TKE in R2 modes changes from ${\approx }45\,\%$ to ${\approx }55\,\%$ between $x/D = 10$ and $100$.

Figure 20(c,d) shows the reconstruction trends for the wake core $\mathcal {P}_{xy}$ term in the $ {\textit {Fr}} = 2$ and $10$ wakes, respectively, along with its actual variation obtained from temporal averaging (shown in red). The wake core $\mathcal {P}_{xy}$ for $ {\textit {Fr}} = 2$ decays as $x^{-0.85}$ throughout the spatial domain under consideration (figure 20c). Both R1 and R2 provide very good reconstruction of $\mathcal {P}_{xy}$ beyond $x/D \approx 30$ and exhibit better approximations relative to that for TKE. With its additional modes, R2 follows the behaviour of the actual value of $\mathcal {P}_{xy}$ very closely. The wake core $\mathcal {P}_{xy}$ for the $ {\textit {Fr}} = 10$ wake shows a faster decay rate of $x^{-2}$ in $10 \leqslant x/D \leqslant 70$ (figure 20d). Beyond $x/D \approx 80$, it decays at a slower rate of $x^{-0.52}$. Similar to the $ {\textit {Fr}} = 2$ wake, R2 reconstructs the actual wake core $\mathcal {P}_{xy}$ very well.

The visually good reconstruction of $\mathcal {P}_{xy}$ by the R2 set of modes can be quantified for both wakes. At $x/D = 10$ and $100$, R2 already accounts for $\approx 66\,\%$ and $\approx 92\,\%$ of the actual $\mathcal {P}_{xy}$, respectively, for the $ {\textit {Fr}} = 2$ wake. For the $ {\textit {Fr}} = 10$ wake, the R2 set of modes capture $\approx 80\,\%$ of the actual $\mathcal {P}_{xy}$ at both $x/D = 10$ and $100$. The SPOD modes provide a better low-order truncation for the lateral production as compared with the TKE for both wakes. This is similar to the trend observed by Nidhan et al. (Reference Nidhan, Chongsiripinyo, Schmidt and Sarkar2020) for the unstratified wake at the same $ {\textit {Re}}$.