3.1. Phase-averaged statistics

Vortex shedding and the oscillatory motion of the separated flow region is a distinguishing, general feature of the flow. Decomposing the flow in a periodic and random part can aid in the characterisation of this oscillatory motion. Phase-averaged statistics are obtained based on the temporal evolution of the lift coefficient, which exhibits regular oscillations with a frequency corresponding to the vortex shedding. Around each lift extremum, a tolerance window of duration $\pm 1 D/U_0$ was defined. Snapshots near lift minima were mirrored about the horizontal symmetry plane $y = -0.5$ and averaged together with those near lift maxima to produce statistically coherent phase-averaged fields. In total, $35$ and $21$ snapshots were used for the cases at ${\textit{Re}}=3000$ and ${\textit{Re}}=14\,000$ , respectively. The two phases analysed, $\phi _{\textit{max}}$ and $\phi _{\textit{min}}$ , correspond to the averaged solution obtained on the upper and lower sides of the body, respectively. The case at ${\textit{Re}}=8000$ is omitted for brevity, as it presents intermediate behaviour.

Figure 3.

Phase-averaged velocity streamlines and pressure field at (a) ${\textit{Re}}=3000$ and (b) ${\textit{Re}}=14\,000$ . The upper side of the rectangle corresponds to phase $\phi _{\textit{max}}$ , while the lower side corresponds to $\phi _{\textit{min}}$ . Cross symbols mark the mean reattachment points on the body sides.

Figure 4.

Phase-averaged distributions of (a) skin friction coefficient, (b) Nusselt number, (c) pressure coefficient and (d) standard deviation of pressure coefficient. Line styles denote results at ${\textit{Re}} = 3000$ (dotted line) and ${\textit{Re}} = 14\,000$ (solid line). Curve colours indicate averages over all samples (black), and at two phases $\phi _{\textit{max}}$ (red) and $\phi _{\textit{min}}$ (blue).

The phase-averaged flow structure is displayed in figure 3, showing velocity streamlines superimposed to the pressure field for ${\textit{Re}} = 3000$ and $14\,000$ . In phase with the streamwise oscillation of the low-pressure region, the separation bubbles along the rectangle sides periodically elongate and contract, while its height remains nearly constant. A positive peak in the lift coefficient corresponds to a longer separation bubble on the upper side. The displacement of the reattachment point increases with ${\textit{Re}}$ . As displayed in figure 4(a), showing the phase-averaged friction coefficient $c_f$ , the reattachment coordinate $x_r$ is in the range $4.04 \lt x_r \lt 4.29$ at ${\textit{Re}}=3000$ and within $4.07 \lt x_r \lt 4.55$ at ${\textit{Re}}=14\,000$ . When the separation bubble shortens (at $\phi _{\textit{min}}$ ), the reverse flow strengthens, as apparent from the reduced skin friction in the interval $2\lt x\lt 4$ . However, separation of the reverse boundary layer is observed at both phases. This indicates that strong reverse flow does not penetrate into the upstream portion of the backflow.

The unsteadiness of the separation bubble affects scalar transfer at the wall, as shown by the local Nusselt number in figure 4(b). A longer separation bubble at $\phi _{\textit{max}}$ leads to enhanced scalar near reattachment and after it, while the reduced bubble at $\phi _{\textit{min}}$ corresponds to diminished transfer with respect to the mean. This difference is attributed to vortices impinging on the wall, which enhance the wall-normal scalar gradient more effectively than the reattached turbulent boundary layer. These impingements increase wall-normal velocity fluctuations and promote entrainment of fluid at low scalar concentration from the free stream. The phase difference in scalar transfer is amplified as ${\textit{Re}}$ increases.

Figure 4(c) shows the phase-averaged pressure coefficient $c_p$ . The pressure difference between the two phases generates the fluctuating lift force on the body and grows with ${\textit{Re}}$ along the entire lateral face, consistent with the increase in lift amplitude reported by Cimarelli et al. (Reference Cimarelli, Corsini and Stalio2024) (see their figure 5d). While pressure minima maintain similar magnitudes across phases, their positions and pressure recovery rates change. At $\phi _{\textit{max}}$ (red curves), the pressure minimum shifts upstream with increasing ${\textit{Re}}$ , and the recovery becomes smoother, leading to lower pressure levels in the reattachment region. At $\phi _{\textit{min}}$ (blue curves), pressure is higher near the leading edge and recovers more steeply downstream, resulting in higher pressure in the reattached flow. The $c_p$ distribution motivates also the persistence of the secondary separation bubble, as the more energetic reverse flow is accompained by a stronger adverse pressure gradient.

Figure 4(d) shows the standard deviation of the phase-averaged pressure coefficient, $c_{p_{\textit{std}}}$ , which represents the amplitude of random pressure fluctuations at a fixed phase. According to phase-averaged decomposition (Cantwell & Coles Reference Cantwell and Coles1983), total fluctuations from the global mean include contributions from periodic large-scale motion and from local random turbulence: $ p = P + \widetilde {P} + p^{\prime \prime }$ , where $p^\prime = \widetilde {P} + p^{\prime \prime }$ , $\widetilde {P}$ is the periodic mean component and $p^{\prime \prime }$ is the random component. Fluctuations at $\phi _{\textit{max}}$ and $\phi _{\textit{min}}$ are similar and decrease in intensity with increasing ${\textit{Re}}$ . The difference between the global mean (black curves) and phase-locked results (coloured curves) arises from the contribution of large-scale periodic fluctuations. The growing gap with ${\textit{Re}}$ suggests that periodic, coherent unsteadiness becomes increasingly dominant relative to random turbulence.

4.1. Mean flow

The analysis begins with the mean momentum distribution around the 5 : 1 rectangular cylinder. Effects of a Reynolds number increase on the turbulent flow are considered. A variation in the Reynolds number changes the location of laminar to turbulent transition, affecting the pressure distribution and the momentum transfer along the shear layer and in the base region (Chapman, Kuehn & Larson Reference Chapman, Kuehn and Larson1958).

Figure 5.

Mean (a,c) streamwise and (b,d) vertical velocity components at (a,b) ${\textit{Re}}=3000$ and (c,d) $14\,000$ . Dotted and dashed lines, from top to bottom, designate the turbulent/non-turbulent interface $y_\varOmega$ , the shear-layer centreline $y_{sl}$ and the boundary of the backflow region $y_b$ . Cross symbols indicate the locations of maximum and minimum values.

Figure 5 shows the contours of the streamwise ( $U$ ) and vertical ( $V$ ) mean velocities at ${\textit{Re}}=3000$ and $14\,000$ , where the crosses indicate velocity extrema. To provide a support for the interpretation of statistical quantities above the rectangle, figure 5, as well as most of those following, include indications of the turbulent/non-turbulent interface of the leading-edge shear layer $y_\varOmega (x)$ , the shear-layer centreline $y_{sl}(x)$ , the line encompassing the backflow region $y_b(x)$ and the dividing streamline $y_\psi (x)$ . The turbulent/non-turbulent interface $y_\varOmega$ can be interpreted as the upper boundary of the shear layer, as it lies between the free flow and the sheared region. It is defined here as the set of outermost points of the rotational region with absolute values of vorticity greater than $0.01|\varOmega _{z}|_{\textit{max}}(x)$ , where $|\varOmega _{z}|_{\textit{max}}(x)$ is the maximum absolute value of the mean spanwise vorticity at each $x$ -location. The shear-layer centreline $y_{sl}$ is identified by the position of $|\varOmega _{z}|_{\textit{max}}(x)$ . Concerning the boundary of the backflow region $y_b$ , it separates regions of the flow with different characteristics beneath the shear layer, as will be shown later, and is representative of the recirculation region. This boundary corresponds to the line where the mean streamwise velocity is zero, $U=0$ . Finally, the dividing streamline $y_\psi$ is the mean streamline originating from the leading edge, and thus connects the separation point of the shear layer to its mean reattachment point.

Two distinct regions of negative streamwise velocity are observed, one above the rectangle and the other in its wake. These regions correspond to two mean separation bubbles induced by the leading and trailing edges of the rectangle, hereafter referred to as primary and wake separation bubbles. Minima of the streamwise velocity component are found in the backflow region inside the primary separation bubble, reaching $U_{\textit{min}}=-0.34$ and $-0.37$ for the cases at ${\textit{Re}}=3000$ and $14\,000$ , respectively. In the case at high ${\textit{Re}}$ , low values of streamwise velocity extend more upstream toward the leading edge, due to the influence of a more extended favourable pressure gradient along the wall (see figure 7c of Cimarelli et al. Reference Cimarelli, Corsini and Stalio2024).

The shear layer separating from the leading edge bounds the regions of flow reversal ( $U \lt 0$ ) on the lower side and high streamwise mean velocity ( $U \gt U_0$ ) on the upper side. The maximum value of the streamwise velocity component is attained close to the irrotational flow region, where $U_{\textit{max}}=1.32$ is reached for both the Reynolds numbers considered. While the location of the leading-edge shear layer proves to be almost independent from ${\textit{Re}}$ in the range analysed, its spreading primarily depends on the transition process (and therefore on ${\textit{Re}}$ ) and on the distance from the wall. In figure 5, $y_{sl}$ remains constant with respect to ${\textit{Re}}$ , while $y_\varOmega$ begins to spread at varying distances from the leading edge. The proximity of the wall to one side of the free shear layer limits the entrainment of fluid on that side, inducing a lower pressure in that region. This pressure difference generates a downward curvature in the flow, which further constrains the inflow of fluid into the turbulent shear layer. The process intensifies until the shear layer impinges on the rectangle wall. The mean reattachment point, $x_r$ , nearly coincides with the intersection of $y_b$ with the rectangle wall. These locations are found to move slightly downstream for ${\textit{Re}}$ ranging between $3000$ and $14\,000$ (see figure 10 of Cimarelli et al. Reference Cimarelli, Corsini and Stalio2024), following a behaviour similar to the reattaching flow behind a backward-facing step (Adams & Johnston Reference Adams and Johnston1988a ).

The fluid entrained from the backflow region into the shear layer balances the flow reversed by the pressure rise through the reattachment zone. The flow rate per unit span entrained from the backflow region, $Q_e$ , increases from $0.0689 \, U_0 D$ at ${\textit{Re}} = 3000$ to $0.0754 \, U_0 D$ at ${\textit{Re}} = 14\,000$ . This quantity is computed as the line integral of the mean velocity field $\boldsymbol{U}$ across the curve $\varGamma$ , which represents the upstream boundary of the backflow region:

(4.1) \begin{equation} Q_e = \int _{\varGamma } \boldsymbol{U} \boldsymbol{\cdot }\boldsymbol{\hat {n}} \, {\textrm{d}}s. \end{equation}

In the equation, $\hat {n}$ is the unit vector normal to the curve and ${\textrm{d}}s$ is the arc length element. The curve $\varGamma$ is parametrised as $\varGamma = \lbrace (x, y_b(x)) \, | \, x \in [0, x_0]\rbrace$ , where $y_b(x)$ is the height at which the mean streamwise velocity vanishes and $x_0$ is the streamwise location where the volumetric flux changes sign (indicated by an asterisk in figure 6). Figure 6 shows the volume flux across $y_b$ as a function of the streamwise direction. Positive values of the volume flux indicate fluid flowing outward from the backflow region, while negative values correspond to incoming flow. The slight increase in $Q_e$ with ${\textit{Re}}$ is consistent with the enlargement of the primary separation bubble, as its vertical thickness remains approximately constant over the analysed ${\textit{Re}}$ range. As the Reynolds number increases and the transition of the laminar shear layer shifts upstream, a higher fraction of fluid is entrained from the backflow region into the shear layer close to the leading edge.

Figure 6.

Volume flux across the line of zero streamwise velocity as a function of the streamwise distance normalised with the reattachment length: dotted line, ${\textit{Re}} = 3000$ ; dashed line, ${\textit{Re}} = 8000$ and solid line, ${\textit{Re}} = 14\,000$ . The asterisk symbol marks the zero-crossing point $x_0$ .

Maxima of the vertical mean velocity ( $V_{\textit{max}}=1.14$ and $1.20$ at ${\textit{Re}}=3000$ and $14\,000$ ) are observed at the leading edge, where the laminar boundary layer formed on the front face of the rectangle separates. The minimum value of the vertical velocity component is found above the reattachment region and corresponds to $V_{\textit{min}}=-0.172$ and $-0.156$ at ${\textit{Re}}=3000$ and $14\,000$ , respectively.

4.2. Reynolds stresses

Reynolds stresses reflect how turbulence redistributes momentum, affecting dominant phenomena such as shear-layer growth, reattachment and wake formation. In this section, we analyse their evolution across the Reynolds number range $3000 \lt {\textit{Re}} \lt 14\,000$ , also establishing connections to findings from prior works. Reynolds stresses in the flow around the rectangular 5 : 1 cylinder have been examined in previous numerical studies by Li et al. (Reference Li, Wang, Qiu, Zhou, Fu and Liu2024a ), Cimarelli et al. (Reference Cimarelli, Leonforte and Angeli2018a ) and Chiarini & Quadrio (Reference Chiarini and Quadrio2021), based on DNS data ranging from ${\textit{Re}}=1000$ to $3000$ , as well as in the experimental study by Kumahor & Tachie (Reference Kumahor and Tachie2022), based on PIV measurements at ${\textit{Re}}=16\,200$ .

Figure 7.

Distributions of the Reynolds normal stresses (a,c) $\langle u^{\prime 2} \rangle$ and (b,d) $\langle v^{\prime 2} \rangle$ at (a, b) ${\textit{Re}}=3000$ and (c,d) $14\,000$ . The cross symbols mark the locations of the maximum values.

Figure 8.

Distributions of the (a,c) Reynolds normal stress $\langle w^{\prime 2} \rangle$ and (b,d) shear stress $\langle u^{\prime } v^{\prime } \rangle$ at (a,b) ${\textit{Re}}=3000$ and (c,d) $14\,000$ . The cross symbols mark the locations of maximum value of $\langle w^{\prime 2} \rangle$ and the minimum value of $\langle u^{\prime } v^{\prime } \rangle$ .

In the statistically two-dimensional flow considered, the only non-zero components of the Reynolds-stress tensor are the normal stresses $\langle u^{\prime 2} \rangle$ , $\langle v^{\prime 2} \rangle$ and $\langle w^{\prime 2} \rangle$ , and the shear stress $\langle u^{\prime } v^{\prime } \rangle$ . Throughout this section, we use the term Reynolds stress to refer to quantities $\langle u_i^{\prime } u^{\prime}_{\!j} \rangle$ , excluding the negative sign that is introduced when these terms are transferred from the left-hand side to the right-hand side of the mean momentum equation to represent turbulent stresses. The spatial evolution of the four terms in the present flow configuration is shown in figures 7 and 8 for the cases at ${\textit{Re}} = 3000$ and $14\,000$ . For sharp-edged bluff bodies with uniform incoming flow, the regions of maximum Reynolds stresses are strictly related to the process of laminar to turbulent transition and occur along the turbulent free shear layers. Turbulent fluctuations are initiated by the Kelvin–Helmholtz-like roll-up of the vortex sheet separating from the leading edge and by the pairing of vortices. By increasing the Reynolds number, the instability mechanisms are triggered at shorter distances from the leading edge and turbulence arises more rapidly under the action of mean shears of growing intensity. To illustrate the deformation experienced by the fluid in the shear layer, the principal mean strain rate $S_{\lambda }$ , i.e. the largest eigenvalue of the mean strain-rate tensor $ S_{\textit{ij}} = ({1}/{2}) ( \partial U_i / \partial x_j + \partial U_j / \partial x_i )$ , is evaluated along the shear-layer centreline $y_{sl}$ and shown in figure 9. Note that $S_{\textit{ij}}$ is a second-order symmetric and traceless tensor, thus its eigenvalues are real and opposite, indicating stretching ( $S_{\lambda }$ ) in one direction and compression ( $-S_{\lambda }$ ) in the other. The data in figure 9 are presented on a logarithmic scale to highlight the power-law behaviour of $S_{\lambda }$ . In the initial, laminar stage of the shear layer, $S_{\lambda }$ exhibits larger values with ${\textit{Re}}$ , but converges to similar values at the trailing edge over the range of ${\textit{Re}}$ examined. Two distinct power-law decay regions, marked by the dash-dotted lines in figure 9, are observed. The first region exhibits a weaker decay, associated with transitional turbulence mechanisms, while the second shows a steeper decay, corresponding to the fully turbulent shear layer. This variation indicates a change in the deformation of fluid elements due to turbulence. The decay rates are found to be almost invariant with respect to ${\textit{Re}}$ and can be associated with the different growth stages of the shear layer thickness, as described by Cimarelli et al. (Reference Cimarelli, Corsini and Stalio2024).

Figure 9.

Largest eigenvalue of the mean strain-rate tensor along the leading-edge shear layer centreline: dotted line, ${\textit{Re}} = 3000$ ; dashed line, ${\textit{Re}} = 8000$ ; and solid line, ${\textit{Re}} = 14\,000$ . The dash-dotted lines mark the two power-law decay regions.

The streamwise Reynolds normal stress $\langle u^{\prime 2} \rangle$ is the dominant component in the leading-edge shear layer for all the Reynolds numbers considered. At ${\textit{Re}}=3000$ , high streamwise fluctuations arise in the rear part of the primary separation bubble, reaching a peak of $0.116$ at $(x, y) = (3.37, 0.45)$ . As shown by Cimarelli et al. (Reference Cimarelli, Corsini and Stalio2024), in the case at ${\textit{Re}} = 3000$ , the intense turbulent fluctuations originating in this region are due to the positive combination of the transition process and the shedding of large-scale vortices from the separated region, occurring in close proximity. As the Reynolds number increases, transition mechanisms shift upstream along the shear layer centreline and separate (at least in terms of distances) from the vortex shedding region. Hence, levels of $\langle u^{\prime 2} \rangle$ remain significant across the bulk of the shear layer and aligned with the streamwise direction. The maximum values of $\langle u^{\prime 2} \rangle$ are reached in the early turbulent region and increase with ${\textit{Re}}$ . In detail, the peak values of $0.114$ and $0.122$ are found at $x=1.38$ and $0.79$ for ${\textit{Re}} = 8000$ and $14\,000$ , respectively.

The vertical and spanwise Reynolds normal stresses $\langle v^{\prime 2} \rangle$ and $\langle w^{\prime 2} \rangle$ are characterised by a lower intensity in the leading-edge shear layer with respect to the streamwise component. This difference becomes wider with the growth of the Reynolds number. In particular, while at ${\textit{Re}}=3000$ , the maximum values of $\langle v^{\prime 2} \rangle _{\textit{max}} = 0.090$ and $\langle w^{\prime 2} \rangle _{\textit{max}} = 0.090$ are observed respectively at $(2.96,0.47)$ and $(3.23,0.40)$ , and at ${\textit{Re}}=14\,000$ , $\langle v^{\prime 2} \rangle _{\textit{max}} = 0.067$ and $\langle w^{\prime 2} \rangle _{\textit{max}} = 0.080$ are found at $(0.71, 0.38)$ and $(0.89, 0.42)$ , respectively. Vertical fluctuations on the low-velocity side of the shear layer rapidly fall as the solid wall is approached. Due to the curvature of the mean streamline close to the reattachment region, turbulent motions are conveyed towards the wall and wall-normal fluctuations are suppressed because of the impermeability condition. On the other side, fluctuations parallel to the wall are generated close to the wall in the reattachment region and in the backflow region at decreasing distances from the leading edge with ${\textit{Re}}$ .

The streamwise Reynolds normal stress also prevails in the shear layer separating from the trailing edge of the rectangle. Here, turbulent motions developed in the reattached boundary layer, while transported downstream, are amplified by the mean velocity gradients in the wake region and acquire streamwise turbulent energy. The magnitude of $\langle u^{\prime 2} \rangle$ in this region increases with the Reynolds number. In the case at ${\textit{Re}} = 14\,000$ , ${\langle u^{\prime 2} \rangle }_{\textit{max}} = 0.084$ is obtained at $(6.03, 0.07)$ . Unlike the leading-edge shear layer, high streamwise fluctuations are not centred in the bulk of the trailing-edge shear layer, but are spread across the vertical direction. One of the reasons lies in the vertical oscillatory motion induced on the trailing-edge shear layer by the von Kármán vortex shedding motion occurring in the wake region. Other factors that contribute to the fast growth of the shear-layer thickness are the turbulent state of the boundary layer separating at the trailing edge and the supply of turbulent motions generated upstream in the reattached boundary layer.

The large scale oscillations in the vertical direction occurring in the wake are associated with high vertical fluctuations centred in the wake centreline, i.e. the symmetry plane $y=-0.5$ . More intense oscillations take place at greater Reynolds numbers, as shown by the significant rise of $\langle v^{\prime 2} \rangle$ for $x\gt 5$ in figure 7(b,d). The maximum value in the wake increase from $\langle v^{\prime 2} \rangle _{\textit{max}} = 0.070$ at ${\textit{Re}}=3000$ to $0.074$ and $0.105$ at ${\textit{Re}}=8000$ and $14\,000$ , respectively. Moreover, high levels of $\langle v^{\prime 2} \rangle$ last over longer distances in the wake at larger ${\textit{Re}}$ (not shown in figure).

The Reynolds shear stress $\langle u^{\prime } v^{\prime } \rangle$ provides a relevant contribution to momentum transfer across the flow, as well as to the production of turbulent kinetic energy. As shown in figure 8(b, d), $\langle u^{\prime } v^{\prime } \rangle$ is mostly negative around the rectangle and exhibits high magnitude levels in the core of the turbulent leading-edge shear layer and in the wake. The minimum value, located where streamwise and vertical fluctuations are best anti-correlated, reaches $-0.053$ and is located within the separation bubble at $(3.21, 0.47)$ for ${\textit{Re}}=3000$ , while equals $-0.045$ and is located in the wake region at $(6.89, 0.064)$ for ${\textit{Re}}=14\,000$ . Regions of positive values for the shear stress $\langle u^{\prime } v^{\prime } \rangle$ are limited to the most upstream part of the separation bubble below the shear layer and within the wake separation bubble.

To enable a quantitative comparison of the Reynolds stresses generated in each section $x$ at different Reynolds numbers, we compute the integrals of the Reynolds stress tensor components, $\langle u_{i}^{\prime } u_{j}^{\prime } \rangle$ , across the turbulent flow, from the rectangle wall to the outer edge of the shear layer $y_{\varOmega }$ ,

(4.2) \begin{equation} \gamma _{\textit{ij}}(x) = \int _{0}^{y_{\varOmega }} \big|\big\langle u_{i}^{\prime } u_{j}^{\prime } \big\rangle (x,y)\big| \, {\textrm{d}} y, \end{equation}

where the absolute value ensures that $\gamma _{12}$ is positive. The resulting integrals for the three ${\textit{Re}}$ values investigated are reported in figure 10. Globally, a larger amount of $\langle uu\rangle$ rises along the shear layer at higher Reynolds numbers. However, cross-stream Reynolds normal stresses $\langle v^{\prime 2} \rangle$ and $\langle w^{\prime 2} \rangle$ do not experience a similar increase. Although higher vertical and spanwise fluctuations occur near the leading edge by increasing the Reynolds number, lower values of $\gamma _{22}$ and $\gamma _{33}$ are observed in the downstream half of the rectangle. The same is true for the shear component $\gamma _{12}$ .

Figure 10.

Integrals of the Reynolds stresses across the separated flow above the rectangle, (a) $\gamma _{11}$ , (b) $\gamma _{22}$ , (c) $\gamma _{33}$ and (d) $\gamma _{12}$ : dotted line, ${\textit{Re}} = 3000$ ; dashed line, ${\textit{Re}} = 8000$ ; and solid line, ${\textit{Re}} = 14\,000$ .

4.3. Mean momentum balance

Turbulence significantly affects the mean flow (see § 3) by redistributing momentum through Reynolds stresses. This section quantifies this effect by analysing the mean momentum balance. We highlight how the pressure gradient and the divergence of the Reynolds stress tensor act to bend the shear layer, promote reattachment and sustain the backflow. To isolate these effects, we reformulate the momentum equations in a streamline-aligned coordinate system, which allows direct interpretation of turbulent momentum transport across the mean dividing streamline.

Previous sections compared flow features across three different Reynolds numbers. From this point forward, the analysis focuses on the case at ${\textit{Re}} = 14\,000$ , where the flow approaches an asymptotic high-Reynolds-number regime.

For the statistically two-dimensional and steady flow under consideration, the conservation equations for mean momentum in the streamwise and vertical directions reduce to

(4.3) \begin{equation} \displaystyle \displaystyle U_j \frac {\partial U_i}{\partial x_j} = \displaystyle - \frac {\partial P}{\partial x_i} - \frac {\partial \big\langle u_{i}^{\prime } u_{j}^{\prime } \big\rangle }{\partial x_j}. \end{equation}

The left-hand side represents the mean advection of momentum that is balanced by the gain or loss in mean momentum resulting from the action of the mean pressure field and the Reynolds stresses. The viscous diffusion term is omitted from the mean momentum balance, as it is negligible across most of the domain for the Reynolds numbers considered, except in the initial part of the separated shear layer, where turbulence is still developing and Reynolds stresses remain low. This simplification was verified by comparing its magnitude to those of the pressure gradient and Reynolds stress divergence.

Figure 11.

Field lines of (a) the mean pressure gradient $(-\partial P / \partial x, -\partial P / \partial y)$ and (b) the divergence of the Reynolds stress tensor $(-\partial \langle u^{\prime } u_j^{\prime }\rangle / \partial x_j, -\partial \langle v^{\prime } u_j^{\prime }\rangle / \partial x_j)$ , with $j=1,2$ . The colour represents the magnitude of the vector field. Results are evaluated at ${\textit{Re}} = 14\,000$ . The green line marks the dividing streamline originating from the leading edge $y_\psi$ .

The effect of the mean pressure gradient force on the flow is first examined. Figure 11(a) shows the vector field of the mean pressure gradient around the rectangle. The incoming flow about $y = -0.5$ loses streamwise momentum because of the increasing pressure levels near the stagnation point, $(x,y) = (0, -0.5)$ , and deviates vertically, gaining $y$ -momentum, under the influence of a favourable pressure gradient. Downstream the leading edge, curved mean streamlines which circumvent the front corner are accompanied by a transverse mean pressure gradient that is proportional to the centripetal acceleration. Peak pressure gradients are thus observed near the leading edge, $(x,y) = (0, 0)$ , in the region where high velocities occur together with small radii of curvature of the streamlines. In the shear layer, the pressure gradient field grows in intensity (indicated by a darker shade), starting from the initial stages of the turbulent shear layer. The pressure gradient peak is located at $x \approx 0.4$ at ${\textit{Re}} = 14\,000$ . Here, the pressure gradient force is directed towards the shear-layer centreline, drawing fluid from the outside region. Downstream the centre of the primary separation bubble, $(x,y) = (2.17, 0.33)$ , the mean pressure raises and slows down the streamwise component of the mean flow. The presence of the impermeable wall is felt by the mean flow moving downstream through a vertical adverse pressure gradient. This downwash is observed starting from $x \approx 2.2$ .

The divergence of the Reynolds stress tensor, shown in figure 11(b), represents the momentum changes due to turbulent transport within the flow. Turbulent motion significantly affects the momentum distribution along the separated shear layers and in the vicinity of the reattachment region. In the leading-edge shear layer, two peaks of Reynolds stress divergence are observed on the upper and lower sides. On the upper side, the divergence of Reynolds stresses removes streamwise ( $x$ ) momentum while supplying wall-normal ( $y$ ) momentum. Conversely, on the lower side, $x$ -momentum is gained while $y$ -momentum is drained. Due to the transport of mean momentum by turbulent motion, fluid is pushed from the shear layer centreline to the free stream and the backflow region. Turbulence is responsible for the spreading of the shear layer and for smoothing the longitudinal velocity gradient between its two sides, similar to what is observed in plane mixing layers (Castro & Haque Reference Castro and Haque1987; Rogers & Moser Reference Rogers and Moser1994). The $x$ -momentum lost in the upper side is released in the lower side. As the shear layer spreads and approaches the wall, turbulent transport extracts $y$ -momentum from the region near the wall and promotes flow reattachment.

To clarify the mechanisms regulating shear layer bending and reattachment, we analyse the dominant terms in the mean momentum balance along the dividing streamline, $y_\psi$ , which originates from the leading edge (see the green solid line in figure 11). For this purpose, we examine the mean momentum balance in the streamline coordinate system $(\tau , \eta , z)$ (Finnigan Reference Finnigan1983), where $\tau$ is locally aligned with the mean velocity vector, $z$ corresponds to the Cartesian $z$ -axis, and $\eta$ is orthogonal to both $\tau$ and $z$ by the right-hand rule.

Reformulating the mean momentum equations in the streamline coordinate system offers several advantages. First, the terms in the transformed equations distinctly describe transport along or normal to the streamlines, removing any ambiguity caused by misaligned coordinates relative to the mean flow. In the second instance, it facilitates the direct interpretation of the Reynolds shear stress as the turbulent momentum transport across streamlines. Furthermore, the transformation emphasises quantities related to mean flow acceleration and streamline curvature, which are relevant characteristics of distorted shear flows. Finnigan et al. (Reference Finnigan, Raupach, Bradley and Aldis1990) were the first to apply the streamline flow equations developed by Finnigan (Reference Finnigan1983) to investigate flow over a hill. Recent studies employing streamline coordinate analysis to capture the properties of external flows include those by Morse & Mahesh (Reference Morse and Mahesh2021), Plasseraud, Kumar & Mahesh (Reference Plasseraud, Kumar and Mahesh2023) and Prakash et al. (Reference Prakash, Balin, Evans and Jansen2024).

The mean momentum equations are expressed in the streamline coordinate system as follows:

(4.4) \begin{equation} \left . \begin{array}{rr} \displaystyle U_{\tau } \frac {\partial U_{\tau }}{\partial \tau } &= \displaystyle - \frac {\partial P}{\partial \tau } - \frac {\partial \big\langle u_{\tau }^{\prime }u_{\tau }^{\prime } \big\rangle }{\partial \tau } - \frac {\partial \big\langle u_{\tau }^{\prime }u_{\eta }^{\prime } \big\rangle }{\partial \eta } + \frac {\big\langle u_{\tau }u_{\tau } \big\rangle - \big\langle u_{\eta }^{\prime }u_{\eta }^{\prime } \big\rangle }{L_{a}} + 2 \frac {\big\langle u_{\tau }^{\prime }u_{\eta }^{\prime } \big\rangle }{R},\\[10pt] \displaystyle \frac {U^{2}_{\tau }}{R} &= \displaystyle - \frac {\partial P}{\partial \eta } - \frac {\partial \big\langle u_{\tau }^{\prime }u_{\eta }^{\prime } \big\rangle }{\partial \tau } - \frac {\partial \big\langle u_{\eta }^{\prime }u_{\eta }^{\prime } \big\rangle }{\partial \eta } + 2 \frac {\big\langle u_{\tau }^{\prime }u_{\eta }^{\prime } \big\rangle }{L_{a}} + \frac {\big\langle u_{\tau }^{\prime }u_{\tau }^{\prime } \big\rangle - \big\langle u_{\eta }^{\prime }u_{\eta }^{\prime } \big\rangle }{R}, \end{array}\right \} \end{equation}

where viscous terms have been again neglected. In (4.4), the subscripts denote the directional components of the velocity vector (e.g. $U_{\tau }$ is the mean velocity component in the tangential $\tau$ direction). Here, $L_{a}$ and $R$ are the local radii of curvature of the $\eta$ and $\tau$ coordinate lines, respectively, and are defined by the following relations:

(4.5) \begin{align} \frac {1}{L_{a}} &= \frac {1}{U_{\tau }} \frac {\partial U_{\tau }}{\partial \tau }, \end{align}

(4.6) \begin{align} \frac {1}{R} &= \frac {1}{U_{\tau }}\left ( \varOmega + \frac {\partial U_{\tau }}{\partial \eta } \right )\!. \end{align}

In the streamline and streamline-normal momentum equations (4.4), the left-hand side represents advection along a streamline and the mean centripetal acceleration, respectively. On the right-hand side, additional terms appear in the expression for Reynolds stress divergence due to the curvilinear nature of the coordinate system adopted. For a detailed derivation of the streamline-coordinate momentum equations, the reader is referred to Finnigan (Reference Finnigan1983).

Figure 12.

Mean momentum balance evaluated along the dividing streamline $y_\psi$ . The dotted vertical line indicates the location of mean reattachment $x_r$ . (a) Terms in the streamline momentum equation. (b) Terms in the streamline-normal momentum equation. (c) Vectors of the mean pressure gradient (blue arrows) and the divergence of the Reynolds stress tensor (red arrows) along the dividing streamline $y_\psi$ , outlined by the dashed line.

Figure 12 shows the terms from the streamline and streamline-normal momentum equations along the dividing streamline $y_\psi$ . The mean advection term is substantially driven by the pressure gradient in both the tangential ( $\tau$ ) and normal ( $\eta$ ) directions. Turbulence effects are confined to the early turbulent and reattachment regions, where stress divergence terms gain relevance due to fluctuations production process and wall impermeability. The following analysis focuses on the range $0 \lt x \lt x_r$ .

In the streamline balance (see figure 12 a) the tangential pressure gradient supplies $\tau$ -momentum in the early turbulent region ( $x\lt 1$ ), compensating the negative peak of transport due to Reynolds tangential stress, $-\partial \langle u_{\tau }^{\prime }u_{\tau }^{\prime } \rangle / \partial \tau$ , and inducing a net acceleration along the streamline. At greater streamwise distances, $\tau$ -momentum is removed ( $U_{\tau } \partial U_{\tau } / \partial \tau \lt 0$ ) under the influence of the adverse tangential pressure gradient. When the flow approaches the wall ( $x\gt 3$ ), turbulent transport terms deliver momentum, balancing the pressure gradient at the mean reattachment location. Here, the Reynolds stress gradient $-\partial \langle u_{\eta }^{\prime }u_{\eta }^{\prime } \rangle / \partial \eta$ is the leading term, peaking at $x \simeq 4$ . Explicit effects of the curvature are represented by $(\langle u_{\tau }^{\prime }u_{\tau }^{\prime } \rangle - \langle u_{\eta }^{\prime }u_{\eta }^{\prime } \rangle ) / L_{a}$ , which extracts tangential momentum in the reattachment region.

As the trajectory of the dividing streamline indicates, fluid at the leading edge initially experiences a negative centripetal acceleration, i.e. the local centre of curvature lies in the direction of decreasing $\eta$ . In examining the streamline-normal balance (figure 12 b), the contributions to the centripetal acceleration are dominated by two main terms, the normal pressure gradient and the transport due to normal Reynolds stress, $-\partial \langle u_{\eta }^{\prime }u_{\eta }^{\prime } \rangle / \partial \eta$ . The magnitudes of both terms peak in the early turbulent region, with turbulent transport providing a centrifugal contribution, $-\partial \langle u_{\eta }^{\prime }u_{\eta }^{\prime } \rangle / \partial \eta \gt 0$ , and the pressure gradient contributing centripetally, $-\partial P / \partial \eta \lt 0$ . As the turbulent shear layer evolves ( $x\gt 1$ ), the turbulent transport term becomes negligible and the centripetal acceleration along the curved streamline is dictated solely by the normal pressure gradient. The pressure gradient term decreases in magnitude with the streamwise distance and changes its sign to positive while approaching the wall, inducing a change in streamline curvature ( $U^{2}_{\tau } / R \gt 2$ ) at $x \simeq 4$ . From reattachment onwards, the dividing streamline flattens against the wall, exhibiting zero curvature $1/R \approx 0$ . Here, the normal pressure gradient balances the normal Reynolds stress gradient.

Figure 12(c) provides an overview of the acceleration terms from pressure and turbulent stresses in the mean momentum balance along the dividing streamline. Notice that these two contributions are invariant in Cartesian and streamline coordinate systems, provided all terms emerging in the expressions are combined into a single vector. In general, the Reynolds stress divergence counteracts the pressure gradient, working to increase the thickness of the separation bubble in the early turbulent region of the shear layer and promoting the flow reattachment in the downstream region of the rectangle.

4.4. Reynolds-stress balances

This section investigates the production, redistribution, transport and dissipation of the individual components of the Reynolds-stress tensor. While the previous study by Chiarini & Quadrio (Reference Chiarini and Quadrio2021) has provided a detailed description of the dominant mechanisms at ${\textit{Re}}=3000$ , particularly the role of production in the separated shear layer and intercomponent redistribution via the pressure–strain tensor, here we extend the analysis to the highest Reynolds number case. Attention is given to the scaling and spatial displacement of peak values, as well as changes in the relative weight between production, dissipation and pressure–strain terms as the Reynolds number increases up to ${\textit{Re}} = 14\,000$ .

The budgets for the Reynolds stresses can be written as

(4.7) \begin{equation} \frac {\partial T_{k,{\textit{ij}}}}{\partial x_{k}} = P_{\textit{ij}} - \epsilon _{\textit{ij}} + \varPi _{\textit{ij}}, \end{equation}

where the left-hand side term represents the mean, pressure, turbulent and diffusive transport terms, all of which are combined and expressed as the divergence of the Reynolds-stress flux $T_{k,ij}$ . The Reynolds-stress flux is defined as

(4.8) \begin{equation} T_{k,ij} = U_{k} \big\langle u_{i}^{\prime } u_{j}^{\prime } \big\rangle + \big [\big\langle p^{\prime } u_{i}^{\prime }\big\rangle \delta _{jk} + \big \langle p^{\prime } u_{j}^{\prime }\big\rangle \delta _{ik} \big ] + \big\langle u_{k}^{\prime }u_{i}^{\prime }u_{j}^{\prime }\big\rangle - \frac {1}{{\textit{Re}}} \frac {\partial \big\langle u_{i}^{\prime } u_{j}^{\prime } \big\rangle }{\partial x_{k}}, \end{equation}

where $\delta _{\textit{ij}}$ is Kronecker’s delta. The terms on the right-hand side of (4.7) correspond to the production tensor $P_{\textit{ij}}$ , the dissipation tensor $\epsilon _{\textit{ij}}$ and the pressure–strain tensor $\varPi _{\textit{ij}}$ . These are defined as follows:

(4.9) \begin{align} P_{\textit{ij}} &= -\left\langle u_i^{\prime } u_k^{\prime } \right\rangle \frac {\partial U_{j}}{\partial x_{k}} -\big\langle u_j^{\prime } u_k^{\prime } \big\rangle \frac {\partial U_{i}}{\partial x_{k}}, \end{align}

(4.10) \begin{align} \epsilon _{\textit{ij}} &= \frac {2}{{\textit{Re}}} \left\langle \frac {\partial u_i^{\prime }}{\partial x_k} \frac {\partial u_j^{\prime }}{\partial x_k} \right\rangle\!. \end{align}

(4.11) \begin{align} \varPi _{\textit{ij}} &= \left \langle p^{\prime }\, \frac {\partial u_i^{\prime }}{\partial x_j} \right \rangle + \left \langle p^{\prime }\, \frac {\partial u_j^{\prime }}{\partial x_i} \right \rangle\! . \end{align}

These terms can be combined to form the net source term $\varPsi _{\textit{ij}} = P_{\textit{ij}} - \epsilon _{\textit{ij}} + \varPi _{\textit{ij}}$ , representing the overall source or sink of Reynolds stresses. In the region where $\langle u_{i}^{\prime } u_{j}^{\prime } \rangle \gt 0$ , a negative value of $\varPsi _{\textit{ij}}$ ( $\varPsi _{\textit{ij}} \lt 0$ ) indicates a loss, whereas a positive value ( $\varPsi _{\textit{ij}} \gt 0$ ) indicates a gain. For cases where $\langle u_{i}^{\prime } u_{j}^{\prime } \rangle \lt 0$ , such as for $\langle u^{\prime } v^{\prime } \rangle$ , the opposite is true, i.e. $\varPsi _{\textit{ij}} \lt 0$ denotes a source of the shear stress, corresponding to an increase in the magnitude of the shear stress, and $\varPsi _{\textit{ij}} \gt 0$ indicates a sink. Each term in (4.7) is indexed according to the specific Reynolds stress component under consideration.

The non-zero components of the production tensor $P_{\textit{ij}}$ are shown in figure 13. Note that $P_{33}$ is zero, as the mean spanwise velocity and its spatial derivatives are zero. The production terms attain their highest values within the leading and trailing shear layers, where the work of the Reynolds stress against the mean velocity gradient is enhanced by intense turbulent fluctuations and high mean strain rates. In the leading-edge shear layer, $P_{11}$ acts as a source for $\langle u^{\prime 2} \rangle$ , drawing energy from the mean flow, while $P_{22}$ serves as a sink for $\langle v^{\prime 2} \rangle$ , transferring energy from the fluctuating motion to sustain the mean vertical flow. This reversed energy exchange, from turbulence to mean flow, intensifies with increasing Reynolds number. At ${\textit{Re}} = 3000$ , the negative peak of $P_{22}$ is approximately $9\,\%$ of the maximum of $P_{11}$ . By ${\textit{Re}} = 14\,000$ , the negative production level reaches $34\,\%$ . The increased energy transfer from turbulence to the mean vertical flow likely contributes to the weak Reynolds number sensitivity observed of the thickness of the main separation bubble. As shown by Cimarelli et al. (Reference Cimarelli, Corsini and Stalio2024), the separated region preserves its vertical size across the analysed ${\textit{Re}}$ range. The large-scale oscillations in the lee of the body correspond to a region where the production rate $P_{22}$ is positive, driven by the vertical deceleration of the mean flow ( $\partial V / \partial y \lt 0$ ) and significant vertical fluctuations $\langle v^{\prime 2} \rangle$ . The shear-stress production term $P_{12}$ is generally negative around the rectangle, except for small positive values in the backflow region, where the adverse pressure gradient induces the separation of the backward boundary layer. The negative values of $P_{12}$ , combined with $\langle u^{\prime } v^{\prime } \rangle \lt 0$ (see figure 8 d), indicate areas of positive production of $\langle u^{\prime } v^{\prime } \rangle$ , and thus enhance the shear stress magnitude. Table 1 shows the location and magnitude of peaks of the main components of the production tensor. As the Reynolds number increases, all components exhibit a marked intensification, with $P_{11}$ growing from approximately 0.5 at ${\textit{Re}} = 3000$ to nearly 2.5 at ${\textit{Re}} = 14\,000$ , and the peaks shifting upstream and closer to the wall. The component $P_{22}$ becomes significantly more negative, while $P_{12}$ also increases in magnitude, indicating intensified shear production in the separated region.

Table 1.

Peak values of production tensor terms at increasing Reynolds numbers.

Figure 13.

Production tensor terms (a) $P_{11}$ , (b) $P_{22}$ , and (c) $P_{12}$ at ${\textit{Re}}=14\,000$ . The cross symbols mark the locations of the maximum or minimum values of $P_{\textit{ij}}$ .

Figure 14 displays the non-zero components of the dissipation tensor, $\epsilon _{\textit{ij}}$ , highlighting regions of significant viscous dissipation within the flow. The dissipation $\epsilon _{12}$ is not shown as negligible. High dissipation levels are concentrated within the leading-edge shear layer, reaching peaks in the same locations of turbulence production. As the turbulent shear layer develops downstream and approaches the wall near reattachment, dissipation rates decrease. In comparison to the production and pressure-rate-of-strain tensors, the dissipation terms exhibit relatively lower magnitudes. The peak values of dissipation components, summarised in table 2, increase with the Reynolds number, albeit less steeply than the production terms. The dissipation peaks for $\epsilon _{11}$ and $\epsilon _{33}$ nearly triple from ${\textit{Re}} = 3000$ to ${\textit{Re}} = 14\,000$ , and their locations move upstream and downward, following the earlier onset of transition in the shear layer.

Table 2.

Peak values of dissipation tensor terms at increasing Reynolds numbers.

Figure 14.

Dissipation tensor terms (a) $\epsilon _{11}$ , (b) $\epsilon _{22}$ and (c) $\epsilon _{33}$ at ${\textit{Re}}=14\,000$ . The cross symbols mark the locations of the maximum values of $\epsilon _{\textit{ij}}$ .

Figure 15.

Pressure–strain tensor terms (a) $\varPi _{11}$ , (b) $\varPi _{22}$ , (c) $\varPi _{33}$ and (d) $\varPi _{12}$ at ${\textit{Re}}=14\,000$ . The cross symbols mark the locations of the maximum or minimum values of $\varPi _{\textit{ij}}$ .

Figure 16.

Pressure–strain normal terms $\varPi _{ii}$ evaluated along (a) the leading-edge shear layer centreline and (b) in the near wall region, at $y=0.024$ : solid line, $\varPi _{11}$ ; dashed line, $\varPi _{22}$ ; and dotted line, $\varPi _{33}$ .

In the Reynolds-stress equations (4.7), the pressure–strain tensor accounts for the intercomponent energy transfer and the tendency towards isotropy of the fluctuating motion. Figure 15 shows the pressure–strain tensor terms $\varPi _{\textit{ij}}$ . In the separated shear layer, the pressure–strain correlation redistributes energy from the streamwise component $\langle u^{\prime 2} \rangle$ ( $\varPi _{11}\lt 0$ ) into the cross-flow component $\langle v^{\prime 2} \rangle$ and $\langle w^{\prime 2} \rangle$ ( $\varPi _{22}\gt 0$ and $\varPi _{33}\gt 0$ ). This redistribution reaches a maximum in the early turbulent region of the shear layer. As displayed in figure 16(a), showing the normal terms $\varPi _{ii}$ along the centreline of the shear layer $y_{sl}$ , the energy transferred from $\langle u^{\prime 2} \rangle$ to $\langle v^{\prime 2} \rangle$ is larger than that transferred to $\langle w^{\prime 2} \rangle$ . Along $y_{sl}$ , the peak values $\varPi _{22,\,{max}}$ and $\varPi _{33,\,{max}}$ amount to $82\,\%$ and $26\,\%$ of $|\varPi _{11}|_{\textit{max}}$ , respectively. In the energy budget for $\langle w^{\prime 2} \rangle$ , there is no production term, making $\varPi _{33}$ the sole source of energy for the spanwise fluctuations. Thus, energy from the mean flow is primarily transferred to $\langle u^{\prime 2} \rangle$ , while $\langle v^{\prime 2} \rangle$ and $\langle w^{\prime 2} \rangle$ gain energy from the streamwise fluctuations through the action of pressure fluctuations. Near the wall, the impermeability condition causes the pressure–strain normal terms to dampen the vertical velocity fluctuations and transfer energy from $\langle v^{\prime 2} \rangle$ to the wall-parallel components, $\langle u^{\prime 2} \rangle$ and $\langle w^{\prime 2} \rangle$ . This process is associated with the strong vertical deceleration of vortical structures in the turbulent shear layer as they impinge upon the wall, driven by the reattaching mean flow. Figure 16(b) shows the profiles of $\varPi_{\textit{ii}}$ close to the wall at $y=0.024$ , the distance at which $\overline {y^+} = 10$ . Here, the overline denotes the streamwise average and the superscript $+$ indicates the normalisation with friction units. In this region, the pressure–strain correlations are significant in the reattachment area, though their magnitude is lower than in the separated shear layer. As the Reynolds number increases, the rate of energy drained from $\langle v^{\prime 2} \rangle$ near the wall and redistributed decreases. At ${\textit{Re}}=3000$ (not shown in the figure), a larger portion of energy is transferred to $\langle w^{\prime 2} \rangle$ , but this imbalance decreases with increasing Reynolds number, and at ${\textit{Re}} = 14\,000$ , the energy exchange between $\langle u^{\prime 2} \rangle$ and $\langle w^{\prime 2} \rangle$ becomes nearly equal. For the off-diagonal term $\varPi _{12}$ , it is commonly understood that pressure fluctuations introduce a term in the $\langle u^{\prime } v^{\prime } \rangle$ budget proportional to $-\langle u^{\prime } v^{\prime } \rangle$ , reducing the absolute value of the Reynolds shear stress (Monin & Yaglom Reference Monin and Yaglom1975). This reflects the natural tendency of turbulence to decay towards isotropy. Figure 15(d) shows $\varPi _{12}$ as positive within the core of the leading-edge shear layer, but exhibiting negative values on the upper side for $2\lt x\lt 4$ . Changes in pressure–strain peaks, reported in table 3, indicate a strengthened intercomponent energy redistribution with increasing Reynolds number. In particular, $\varPi _{11}$ exhibits increasingly intense negative values, while $\varPi _{22}$ and $\varPi _{33}$ show stronger positive peaks. For ${\textit{Re}} \lt 14\,000$ , approximately $75\,\%$ of the energy removed from the streamwise component is redirected to the vertical component, suggesting that the relative partitioning of turbulent kinetic energy within the evolving shear layer is preserved as ${\textit{Re}}$ increases.

Table 3.

Peak values of pressure-strain tensor terms at increasing Reynolds numbers.

Figure 17.

Distribution of the net source and sink terms, $\varPsi _{\textit{ij}} = P_{\textit{ij}} - \epsilon _{\textit{ij}} + \varPi _{\textit{ij}}$ overlaid with flux lines of $\boldsymbol{T_{\textit{ij}}}$ : (a) $\varPsi _{11}$ and $\boldsymbol{T_{11}}$ ; (b) $\varPsi _{22}$ and $\boldsymbol{T_{22}}$ ; (c) $\varPsi _{33}$ and $\boldsymbol{T_{33}}$ ; and (d) $\varPsi _{12}$ and $\boldsymbol{T_{12}}$ .

The combined sink and source terms for each of the Reynolds stress budgets (4.7),  $\varPsi _{\textit{ij}}$ , are analysed along with the lines of the Reynolds-stress flux vector $\boldsymbol{T_{\textit{ij}}}$ in figure 17. This representation summarises locations of net source, sink and paths of Reynolds stress transfer. In the $\langle u^{\prime 2} \rangle$ budget, energy is provided to the streamwise fluctuations directly from the mean flow along the free shear layers, and it is removed by $\epsilon _{11}$ and $\varPi _{11}$ , with $\varPi _{11}$ being the dominant sink away from the wall. Fluctuations $\langle u^{\prime 2} \rangle$ generated in the core of the leading-edge shear layer is primarily transported downstream towards the wake while also spreading towards the lower and upper sides of the shear layer. On the upper side, the upward transfer of $\langle u^{\prime 2} \rangle$ is mainly driven by the turbulent transport term, $\langle u^{\prime }u^{\prime }v^{\prime } \rangle$ , balancing the energy loss due to $\varPi _{11}$ in the outer region. On the lower side, the downward transfer of $\langle u^{\prime 2} \rangle$ is induced by both the turbulent transport and advection, ${\langle u^{\prime 2} \rangle } V$ . This downward transfer splits at $x \approx 2.8$ , feeding both upstream recirculation and the downstream vortex shedding regions. In the separation bubble, a singularity point where the $\langle u^{\prime 2} \rangle$ flux lines converge is identified at $(1.60, 0.33)$ , just along the boundary of the backflow region $y_b$ . A second region where $\varPsi _{11}\gt 0$ is observed in the flow impingement area ( $2.5\lt x\lt 4$ ) within the backflow region, resulting from the positive contribution of $\varPi _{11}$ . Along the reattached boundary layer ( $x \gt x_r$ ), $\langle u^{\prime 2} \rangle$ is fed back to the mean flow, $P_{11}\lt 0$ , or dissipated.

In the $\langle v^{\prime 2} \rangle$ budget, $\varPi _{22}$ acts as the primary source within the leading-edge shear layer, resulting in a net increase of $\langle v^{\prime 2} \rangle$ at the expense of $\langle u^{\prime 2} \rangle$ . In the early stages of the turbulent shear layer, negative production $P_{22}$ is outweighed by the positive contribution of $\varPi _{22}$ , except on the lower side of the shear layer, where $\varPsi _{22} \lt 0$ for $0.1\lt x\lt 0.6$ and $y_b \lt y \lt y_{sl}$ . In this area, $\langle v^{\prime 2} \rangle$ transported from the upper side of the shear layer is either returned to the mean flow kinetic energy $VV$ or dissipated by viscosity, $\epsilon _{22}$ . In the separated flow, flux lines curve downward to transfer $\langle v^{\prime 2} \rangle$ towards the wall, where $\varPsi _{22} \lt 0$ . Here, $\varPi _{22}$ is the dominant sink, affected by the wall impermeability constraint. A net rate of gain ( $\varPsi _{22}\gt 0$ ) is obtained along the trailing-edge shear layer and in the wake region where large-scale vertical oscillations occur. This excess of $\langle v^{\prime 2} \rangle$ is carried downstream by advection, $U \langle v^{\prime 2} \rangle$ , and across the wake by the turbulent transport term  $\langle v^{\prime }v^{\prime }v^{\prime } \rangle$ .

The lines of $\langle w^{\prime 2} \rangle$ flux closely follow the mean velocity streamlines, indicating that advection dominates the energy transfer process. The only other significant contributions to the transport of $\langle w^{\prime 2} \rangle$ come from the turbulent fluxes, $\langle w^{\prime }w^{\prime }u^{\prime }\rangle$ and $\langle w^{\prime }w^{\prime }v^{\prime }\rangle$ . The excess of $\langle w^{\prime 2} \rangle$ ( $\varPsi _{33}\gt 0$ ), generated in the core of the leading-edge shear layer through energy exchange with $\langle u^{\prime 2} \rangle$ ( $\varPi _{33}\gt 0$ ), is transported downstream and downward to the lower side of the shear layer, where it dissipates by viscous action, $\epsilon _{33}$ . Near the wall, the net energy gain provided by $\varPi _{33}$ in the impingement region is partly transported downstream along the reattached boundary layer ( $x\gt x_r$ ) and partly flows through the backflow region ( $x\lt x_r$ ), reaching both the leading edge and the upstream part of the shear layer.

In the $\langle u^{\prime } v^{\prime } \rangle$ budget, negative $\varPsi _{12}$ corresponds to a gain where $\langle u^{\prime } v^{\prime } \rangle \lt 0$ , thereby increasing the magnitude of the shear stress. Hence, lines of $\langle u^{\prime } v^{\prime } \rangle$ flux are viewed as carrying $\langle u^{\prime } v^{\prime } \rangle$ from regions where $\varPsi _{12} \lt 0$ to regions where $\varPsi _{12} \gt 0$ , which is opposite to the flux directions shown in figure 17(d). In the core of the leading-edge shear layer, production $P_{12}\lt 0$ exceeds the pressure–strain $\varPi _{12}\gt 0$ , resulting in a net production of $\langle u^{\prime } v^{\prime } \rangle$ . This produced $\langle u^{\prime } v^{\prime } \rangle$ is partially carried downstream by advection and partially outward by pressure and turbulent fluxes. The transport of $\langle u^{\prime } v^{\prime } \rangle$ in the upper side of the shear layer and at the wall is balanced by $\varPi _{12}$ , which serves as a sink. In the recirculating region, $\langle u^{\prime } v^{\prime } \rangle$ flux lines converges to a singularity point slightly downstream of the centre of rotation of the main separation bubble, located at $(1.95, 0.29)$ along $y_b$ .

5.1. Turbulent scalar fluxes

At high Reynolds numbers, and thus at high Péclet numbers ( ${\textit{Pe}} = {\textit{Re}} {\textit{Sc}}$ ), the turbulent motion becomes more effective than molecular diffusion in mixing and transporting the scalar. In these conditions, turbulent scalar fluxes, $\langle u_i^{\prime } \theta ^{\prime } \rangle$ , emerge as the primary mechanism for scalar transport.

Figure 18.

Scalar fluxes (a,c) $\langle u^{\prime } \theta ^{\prime } \rangle$ and (b,d) $\langle v^{\prime } \theta ^{\prime } \rangle$ at (a,b) ${\textit{Re}}=3000$ and (c,d) $14\,000$ . The cross symbols indicate the locations of maximum and minimum values.

Figure 18 shows the streamwise, $\langle u^{\prime } \theta ^{\prime } \rangle$ , and vertical, $\langle v^{\prime } \theta ^{\prime } \rangle$ , scalar fluxes around the rectangle at ${\textit{Re}}=3000$ and ${\textit{Re}}=14\,000$ . Negative streamwise scalar flux, $\langle u^{\prime } \theta ^{\prime } \rangle$ , is observed in the leading-edge shear layer and near thewall in the reattached boundary layer ( $x\gt x_r$ ). For both ${\textit{Re}}$ cases, ${\langle u^{\prime } \theta ^{\prime } \rangle }_{\textit{min}}$ occurs near the trailing-edge flow separation rather than in the leading-edge shear layer. As highlighted by Cimarelli et al. (Reference Cimarelli, Corsini and Stalio2024) in the analysis of scalar variance $\langle \theta ^{\prime } \theta ^{\prime } \rangle$ , while streamwise velocity fluctuations reach their highest values downstream of the early turbulent region of the shear layer, maximum scalar fluctuations arise close to the wall. This is driven by large-scale vortices in the separated flow, which transport low-concentration scalar fluid and impinge on the wall, where $\theta = 1$ . Positive values of $\langle u^{\prime } \theta ^{\prime } \rangle$ appear in the backflow region and the wake separation bubble, with maximum values located near the separation point of the backward boundary layer.

Vertical scalar flux is positive ( ${\langle v^{\prime } \theta ^{\prime } \rangle } \gt 0$ ) around the rectangle, except within the wake separation bubble. High-intensity regions of $\langle v^{\prime } \theta ^{\prime } \rangle$ are found along the leading-edge shear layer, in the reattachment region near the wall, and along the trailing-edge shear layer. Due to the wall impermeability, vertical velocity fluctuations are dampened near the wall, leading ${\langle u^{\prime } \theta ^{\prime } \rangle }_{\textit{max}}$ to occur in the early turbulent region of the shear layer. As Reynolds number increases, the magnitude of scalar fluxes decreases and the gradients near the wall become less steep, indicating reduced turbulence-scalar interaction in this region.

5.2. Mean scalar balance

The role of turbulence in shaping the mean scalar distribution is encapsulated in the mean scalar balance. The conservation equation for the mean passive scalar, $\varTheta$ , in statistically steady and spanwise homogeneous flows is given by

(5.1) \begin{equation} \displaystyle U_i \frac {\partial \varTheta }{\partial x_i} = \frac {\partial }{\partial x_i} \left ( - \left\langle u_i^{\prime } \theta ^{\prime }\right \rangle + \frac {1}{{\textit{Pe}}} \frac {\partial \varTheta }{\partial x_i} \right )\!. \end{equation}

The term on the left-hand side represents the rate of change of $\varTheta$ due to advection, while terms on the right-hand side describe the transport of $\varTheta$ through turbulent and molecular diffusion. In this analysis, the Péclet number ( ${\textit{Pe}}=9940$ ) is sufficiently high to make molecular diffusion negligible in regions away from the wall and outside the laminar section of the separated shear layer. The rate of change in $\varTheta$ due to advection is balanced by turbulent transport.

Figure 19.

Scalar transport term $-(\partial {\langle u^{\prime } \theta ^{\prime } \rangle } / \partial x +\partial {\langle v^{\prime } \theta ^{\prime } \rangle } / \partial y)$ with flux lines of $\langle \boldsymbol{u}^{\prime }\theta ^{\prime } \rangle$ at ${\textit{Re}} = 14\,000$ .

Figure 19 shows the turbulent transport term $-\partial \langle u_i^{\prime } \theta ^{\prime } \rangle / \partial x_i$ from (5.1), where positive values indicate a net addition of the mean scalar, while negative values denote scalar loss. The superimposed lines of the scalar flux vector, $\langle \boldsymbol{u}^{\prime }\theta ^{\prime } \rangle$ , display the spatial redistribution of $\varTheta$ , with scalar flux lines originating in regions of scalar removal and terminating in regions of scalar supply.

The distribution of the turbulent transport term varies significantly across different regions of the turbulent flow. In the leading-edge shear layer, turbulent motion carries  $\varTheta$ from the lower side to the upper side, with a more intense transfer observed in the early stages of turbulence development ( $x\lt 1$ ). Scalar removed from the wall region is redistributed into surrounding areas, specifically within the backflow region and along the reattaching boundary layer. In the reattachment region ( $3\lt x\lt 5$ ), the scalar flux is strong enough to carry $\varTheta$ toward the outer edge of the turbulent shear layer. In the upstream portion of the backflow region ( $0\lt x\lt 1$ ), turbulent transport drives $\varTheta$ upward through the turbulent shear layer. Downstream of the rectangle, the scalar flux draws $\varTheta$ from the trailing-edge shear layer and the core of the wake, redistributing it into the outer turbulent region of the flow. This analysis highlights the primary role of the scalar fluxes $\langle u_i^{\prime } \theta ^{\prime } \rangle$ in enhancing scalar mixing across the separated flow and the wake, promoting more uniform scalar distribution throughout these flow regions.

5.3. Scalar flux balances

To better understand how turbulent scalar fluxes are generated, redistributed, and removed, we examine the transport equations for $\langle u_i^{\prime } \theta ^{\prime } \rangle$ . These balances provide a detailed view of the mechanisms governing scalar-velocity correlations, which are responsible for turbulent scalar transport. This section analyses the individual budget terms to identify the dominant processes in separated and reattaching regions.

The balance equations of scalar flux $\langle u_i^{\prime } \theta ^{\prime } \rangle$ for steady conditions are expressed as

(5.2) \begin{equation} \frac {\partial T_{k,i\theta }}{\partial x_{k}} = P_{i\theta } - \epsilon _{i\theta } + \varPi _{i\theta }. \end{equation}

In analogy to (4.7), the left-hand side collects the spatial transport terms arising from advection, pressure, turbulence, and molecular diffusion, represented as the divergence of the flux $T_{k,i\theta }$ ,

(5.3) \begin{equation} T_{k,i\theta } = U_{k} \left\langle u_i^{\prime } \theta ^{\prime } \right\rangle + \langle p^{\prime } \theta ^{\prime } \rangle \delta _{ik} + \left \langle u_{k}^{\prime }u_{i}^{\prime }\theta ^{\prime }\right \rangle - \frac {1}{{\textit{Re}}{\textit{Sc}}} \left \langle u_i^{\prime } \frac {\partial \theta ^{\prime }}{\partial x_k} \right \rangle - \frac {1}{{\textit{Re}}} \left \langle \theta ^{\prime } \frac {\partial u_i^{\prime }}{\partial x_k} \right \rangle \!. \end{equation}

On the right-hand side of (5.2), $P_{i\theta }$ denotes the production term, representing the rate of generation of $\langle u_i^{\prime }\theta ^{\prime } \rangle$ through turbulent interactions with mean scalar and mean velocity gradients. The term $\epsilon _{i\theta }$ represents destruction, which quantifies the mean rate at which $\langle u_i^{\prime }\theta ^{\prime } \rangle$ is reduced by molecular diffusivity of momentum and scalar. This term is zero in isotropic turbulence and is usually neglected in non-isotropic turbulence at high Reynolds numbers or under the assumption of local isotropy (Launder Reference Launder1976). Finally, $\varPi _{i\theta }$ represents the pressure-scalar gradient correlation and it may be viewed as the counterpart of the pressure-strain term in the $\langle u_{i}^{\prime } u_{j}^{\prime } \rangle$ budget equations. However, since the sum of the $\varPi _{i\theta }$ terms may not be zero, the pressure-scalar gradient correlation does not necessarily redistribute scalar flux among components. These terms are defined as follows:

(5.4) \begin{align} P_{i\theta } &= -\left\langle u_i^{\prime } u_k^{\prime } \right\rangle \frac {\partial \varTheta }{\partial x_{k}} -\left\langle u_k^{\prime } \theta ^{\prime } \right\rangle \frac {\partial U_{i}}{\partial x_{k}}, \end{align}

(5.5) \begin{align} \epsilon _{i\theta } &= \left ( \frac {1}{{\textit{Re}}} + \frac {1}{{\textit{Re}}{\textit{Pr}}} \right ) \left \langle \frac {\partial u_i^{\prime }}{\partial x_k} \frac {\partial \theta ^{\prime }}{\partial x_k} \right \rangle\! , \end{align}

(5.6) \begin{align} \varPi _{i\theta } &= \left \langle p^{\prime }\, \frac {\partial \theta ^{\prime }}{\partial x_i} \right \rangle\! . \end{align}

Together, these terms form the net source term $\varPsi _{i\theta } = P_{i\theta }-\epsilon _{i\theta }+\varPi _{i\theta }$ representing the net scalar flux that is lost or gained. Similar to $\varPsi _{\textit{ij}}$ , the role of $\varPsi _{i\theta }$ as a source or sink depends on the sign of $\langle u_i^{\prime }\theta ^{\prime } \rangle$ . Each term in (5.2) is indexed to indicate the direction of the scalar flux under consideration ( $i=1$ for $\langle u^{\prime } \theta ^{\prime } \rangle$ or $i=2$ for $\langle v^{\prime } \theta ^{\prime } \rangle$ ).

Figure 20.

Production terms of scalar flux (a) $P_{1\theta }$ and (b) $P_{2\theta }$ at ${\textit{Re}}=14\,000$ . The cross symbols mark the locations of the minimum value of $P_{1\theta }$ and maximum value of $P_{2\theta }$ . The green line denotes the dividing streamline originating from the leading edge, $y_\psi$ .

Figure 20 shows regions of positive and negative production rates of scalar flux, $P_{i\theta }$ , for the flow at ${\textit{Re}}=14\,000$ . The scalar flux $\langle u^{\prime } \theta ^{\prime } \rangle$ is primarily generated within the leading-edge shear layer. In this area, negative values of $P_{1\theta }$ act as a source since ${\langle u^{\prime } \theta ^{\prime } \rangle } \lt 0$ . Near the wall, in the backflow region, positive values of $P_{1\theta }$ correspond to zones where ${\langle u^{\prime } \theta ^{\prime } \rangle } \gt 0$ and hence are sources of $\langle u^{\prime } \theta ^{\prime } \rangle$ . The minimum value $P_{1\theta ,\,{min}}=-0.67$ is found at $(x, y) = (0.40, 0.30)$ . The term $P_{2\theta }$ acts as a source for $\langle v^{\prime } \theta ^{\prime } \rangle$ above the rectangle, with peak intensity along the shear-layer centreline, with $P_{2\theta ,\,{max}}=0.24$ at $(x, y) = (0.49, 0.34)$ , and in the reattachment region ( $2.2\lt x\lt 5$ ).

Figure 21.

Production terms of (a) $\langle u^{\prime } \theta ^{\prime } \rangle$ and (b) $\langle v^{\prime } \theta ^{\prime } \rangle$ evaluated along the reattaching streamline $y_\psi$ : solid line, $P_{i\theta }$ ; dashed line, $P_{i\theta ,a}$ ; and dotted line, $P_{i\theta ,b}$ .

To elucidate the production mechanism of scalar fluxes, the individual contributions to $P_{1\theta }$ and $P_{2\theta }$ are examined. Each production term is expressed as the sum of two components, representing the effects of mean scalar gradients and mean velocity gradients interacting with turbulent motions, as follows:

(5.7) \begin{equation} P_{i\theta } = \underbrace {-\left\langle u_i^{\prime } u_k^{\prime } \right\rangle \frac {\partial \varTheta }{\partial x_{k}}}_{P_{i\theta ,a}} \underbrace {-\left\langle \theta ^{\prime } u_k^{\prime } \right\rangle \frac {\partial U_{i}}{\partial x_{k}}}_{P_{i\theta ,b}}. \end{equation}

Figure 21 illustrates the two terms $P_{i\theta ,a}$ and $P_{i\theta ,b}$ along with the total production rate $P_{i\theta }$ , evaluated along the dividing streamline $y_\psi$ . The production term $P_{1\theta }$ is evenly affected by the mean scalar and velocity gradients ( $P_{1\theta ,a} \approx P_{2\theta ,b}$ ) along the free shear layer ( $x\lt 3$ ) and in the reattached flow region ( $x\gt 4$ ). In contrast, the dominant contribution to $P_{2\theta }$ arises from $P_{2\theta ,a}$ , driven by the interaction of the wall-normal scalar gradient with vertical Reynolds stress, $-\langle v^{\prime 2} \rangle \partial \varTheta /\partial y$ . The contribution due to the mean velocity gradient term, $P_{2\theta ,b}$ , results in a loss of $\langle v^{\prime } \theta ^{\prime } \rangle$ within the shear layer and becomes negligible in the reattached region. This opposing effect reduces the magnitude of $P_{2\theta }$ .

Figure 22.

Destruction terms of scalar flux (a) $\epsilon _{1\theta }$ and (b) $\epsilon _{2\theta }$ at ${\textit{Re}}=14\,000$ . The cross symbols mark the locations of the minimum values of $\epsilon _{i\theta }$ . (c) Detailed view of the distribution of $\epsilon _{2\theta }$ .

Assuming the Reynolds number is high enough, the small scales can be considered locally isotropic, making the destruction terms $\epsilon _{i\theta }$ negligible compared with $P_{i\theta }$ and  $\varPi _{i\theta }$ . However, this assumption does not fully account for the energy transfer to higher wavenumbers, a process that requires some level of destruction (Younis, Speziale & Clark Reference Younis, Speziale and Clark1995). Furthermore, DNS results from Rogers, Mansour & Reynolds (Reference Rogers, Mansour and Reynolds1989), obtained for a fully developed turbulent channel flow, indicate that $\epsilon _{i\theta }$ can act as a production term under certain conditions. Figure 22 shows the destruction terms $\epsilon _{i\theta }$ for the case under consideration. These terms are approximately an order of magnitude smaller than $P_{i\theta }$ and $\varPi _{i\theta }$ , and the colourbar limits have been adjusted accordingly for clarity. In the streamwise balance, $\epsilon _{1\theta }$ is a sink for $\langle u^{\prime } \theta ^{\prime } \rangle$ , since $-\epsilon _{1\theta }\gt 0$ (or $-\epsilon _{1\theta }\lt 0$ ) where ${\langle u^{\prime } \theta ^{\prime } \rangle }\lt 0$ (or ${\langle u^{\prime } \theta ^{\prime } \rangle }\gt 0$ ). The minimum value, $\epsilon _{1\theta ,\,{min}} = -0.015$ , occurs along the shear-layer centreline $y_{sl}$ at $(0.51, 0.35)$ . For the vertical balance, $\epsilon _{2\theta }$ predominantly leads to a loss of $\langle v^{\prime } \theta ^{\prime } \rangle$ around the rectangle, except for a localised zone along $y_{sl}$ ( $0.2\lt x\lt 0.6$ ) where $\epsilon _{2\theta }\lt 0$ is observed, see figure 22(c) for the detailed view. This negative region, with $\epsilon _{2\theta ,\,{\textit{min}}} = -0.0087$ at $(0.37, 0.28)$ , corresponds to positive values of $\langle v^{\prime } \theta ^{\prime } \rangle$ , indicating a production of scalar flux due to viscous effects. It is worth noting that this localised gain in scalar flux contributes significantly less than other sink and source terms in the same regions. For instance, at the location of $\epsilon _{2\theta ,\,{\textit{min}}}$ , $P_{2\theta }=0.103$ .

With the viscous destruction negligible, the pressure-scalar gradient correlations $\varPi _{i\theta }$ provide the mechanism that limits the growth of the scalar fluxes by counteracting the rate of production of $\langle u_i^{\prime }\theta ^{\prime } \rangle$ . Figure 23 shows the distribution of $\varPi _{i\theta }$ . The streamwise term $\varPi _{1\theta }$ is found to oppose to $P_{1\theta }$ discussed previously, serving as a sink of $\langle u^{\prime } \theta ^{\prime } \rangle$ along the free shear layer and in the reattachment zone. The positive peak $P_{1\theta ,\,{max}}$ is found at $(0.43, 0.31)$ . Regarding $\varPi _{2\theta }$ , the centreline $y_{sl}$ approximately divides a gain region ( $\varPi _{2\theta }\gt 0$ ) on the upper side of the shear layer from a loss region ( $\varPi _{2\theta }\lt 0$ ) on the lower side. In the latter, the negative peak occurs at $(0.43, 0.31)$ , with $P_{2\theta ,\,{\textit{min}}} = -0.34$ .

Figure 23.

Pressure-scalar gradient correlation of scalar flux (a) $\varPi _{1\theta }$ and (b) $\varPi _{2\theta }$ at ${\textit{Re}}=14\,000$ . The cross symbols mark the locations of the maximum value of $\varPi _{1\theta }$ and minimum value of $\varPi _{2\theta }$ .

Figure 24.

Distribution of the net source terms of the scalar flux, $\varPsi _{i\theta } = P_{i\theta } + \varPi _{i\theta } - \epsilon _{i\theta }$ overlaid with flux lines of $\boldsymbol{T_{i\theta }}$ : (a) $\varPsi _{1\theta }$ and (b) $\varPsi _{2\theta }$ .

The net sources of the streamwise and vertical scalar fluxes, $\varPsi _{1\theta }$ and $\varPsi _{2\theta }$ , are analysed along with lines of the spatial flux vector $\boldsymbol{T_{i\theta }}$ shown in figure 24. Regions where $\varPsi _{i\theta }\gt 0$ ( $\varPsi _{i\theta }\lt 0$ ) indicate sources for the scalar flux if $\langle u_i^{\prime } \theta ^{\prime } \rangle \gt 0$ ( $\langle u_i^{\prime } \theta ^{\prime } \rangle \lt 0$ ), corresponding to areas of net increase in scalar flux magnitude. In the shear layer and wake, where $\langle u^{\prime } \theta ^{\prime } \rangle$ is predominantly negative, the transport lines of scalar flux in figure 24(a) indicate flux $\langle u^{\prime } \theta ^{\prime } \rangle$ being carried from source regions ( $\varPsi _{1\theta }\lt 0$ ) on the upper side of the shear layer to the lower side and into the wake, moving opposite to the direction indicated by the arrows in the figure. Within the backflow region, the excess of $\langle u^{\prime } \theta ^{\prime } \rangle$ generated near the wall is transported upstream towards the leading edge of the rectangle and upward, beginning at $x \approx 1.5$ .

For the vertical flux budget, the flux $\langle v^{\prime } \theta ^{\prime } \rangle$ gained in the upper side of the shear layer $\varPsi _{2\theta } \gt 0$ is redistributed downward towards the the lower side and into the wake, where it is ultimately removed. In the near-wall region, in particular within $2\lt x\lt 5$ , a strong source of $\langle v^{\prime } \theta ^{\prime } \rangle$ is evident, associated with the impingement of turbulent vortices formed in the shear layer. The flux $\langle v^{\prime } \theta ^{\prime } \rangle$ generated near the wall is subsequently transported to the lower side of the shear layer, where it is drained ( $\varPsi _{2\theta } \lt 0$ ).