Hostname: page-component-6766d58669-88psn Total loading time: 0 Render date: 2026-05-18T02:07:50.914Z Has data issue: false hasContentIssue false

The unsteady wake transition behind a wall-mounted large-depth-ratio prism

Published online by Cambridge University Press:  23 November 2022

Arash Zargar
Affiliation:
Department of Mechanical Engineering, University of Alberta, Edmonton, AB T6G 2R3, Canada
Ali Tarokh
Affiliation:
Department of Mechanical Engineering, University of Alberta, Edmonton, AB T6G 2R3, Canada Department of Mechanical Engineering, Lakehead University, Thunder Bay, ON P7B 5E1, Canada
Arman Hemmati*
Affiliation:
Department of Mechanical Engineering, University of Alberta, Edmonton, AB T6G 2R3, Canada
*
Email address for correspondence: arman.hemmati@ualberta.ca

Abstract

The wake of a long rectangular wall-mounted prism is investigated at Reynolds numbers of $Re=250\unicode{x2013}1200$ by directly solving the Navier–Stokes equations. The aim of this study was to examine the unsteady transition mechanism in the wake of a large-depth-ratio (streamwise length to width) prism as well as to characterize the unsteady wake evolution at low Reynolds numbers. The results highlighted that increasing Reynolds number significantly altered the dominance of upwash and downwash flows in the time-averaged flow and changed the characteristics of coherent structures, including their size, dominant frequency and interaction with other structures in the flow. The wake is, therefore, categorized into three regimes within the transition process: steady regime at $Re \leq 625$, regular unsteady regime at $625 < Re < 750$ and irregular unsteady regime at $Re \geq 750$. Particularly, the wake started to exhibit unsteady features at $Re=625\unicode{x2013}650$, which transitioned to an early irregular unsteady wake topology at $Re=750$. At $Re \geq 675$, horseshoe vortices transformed to vortex loops. There were hairpin-like structures formed on the upper and side faces of the long prism. The results further indicated that the transition to unsteadiness is attributed to separated leading-edge shear-layer instabilities and interactions of the shear layer with the horseshoe structures. The wake was more complex due to the interactions of multiple coherent structures in the flow, which resulted in a multiple-periodicity wake system. A skeleton model is proposed for a large-depth-ratio prism, to incorporate details of the unsteady flow features and specify various flow coherent structures at low Reynolds numbers.

Information

Type
JFM Papers
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press.
Figure 0

Figure 1. Wake models for wall-mounted square prisms: (a) the original model developed by Wang et al. (2004); (b) the modified model by Wang & Zhou (2009).

Figure 1

Figure 2. Schematic of the computational domain (not to scale).

Figure 2

Table 1. Boundary layer thickness at the prism leading edge ($x/h = -4.17$ and $y/h = -7.5$) at different Reynolds numbers.

Figure 3

Figure 3. The spatial grid distribution for the long rectangular prism at $Re=1000$. Top view at $z/h=0.5$ (a) and side view at $y/h=0$ (b).

Figure 4

Figure 4. Comparing the profile of time-averaged streamwise velocity ($\bar {u}/U_0$) in the wake of the long rectangular prism at $Re=1000, x/h=2$ and $z/h=0.5$ using (a) different domain sizes and (b) different grid refinements.

Figure 5

Figure 5. Grid refinement analysis for the wake of the long rectangular prism at $Re = 250$ using (a) time-averaged coefficient of pressure distribution ($\overline {c_p} = {\bar {p}}/{0.5\rho U_0^2}$) on the surface of the prism at its middle height ($z/h = 0.5$) and (b) time-averaged streamwise velocity ($\bar {u}/U_0$) at $x/h=2$ and $z/h=0.5$.

Figure 6

Figure 6. The ratio of grid size, defined as $\varDelta _{max}=\max \{ \delta x , \delta y , \delta z \}$, to Kolmogorov scale at (a) $z/h = 0.5$ and (b) $y/h = 0$.

Figure 7

Figure 7. Contour of the ratio of grid size to Kolmogorov length scale, defined as $\varDelta =(\delta x\times \delta y \times \delta z)^{1/3}/\eta$, for (a) the time-averaged field at $Re=1000$ and (b) instantaneous field at $Re=1200$ along the $xz$ plane at (top) $z/h=0.5$ and (bottom) $y/h=0$.

Figure 8

Figure 8. Isosurface of $Q^*=0.18$ in the wake at an arbitrary time for (a) $Re = 250$, (b) $Re = 350$, (c) $Re = 500$, (d) $Re = 575$, (e) $Re = 600$ and ( f) $Re = 625$. Here, separation bubbles are denoted by ‘SB’ and horseshoe vortex and vortex loops are respectively highlighted as ‘HV’ and ‘VL’.

Figure 9

Figure 9. Isosurface of $Q^*=0.18$ at an arbitrary time at (a) $Re = 650$, (b) $Re = 675$, (c) $Re = 700$, (d) $Re = 725$, (e) $Re = 750$, ( f) $Re = 850$, (g) $Re = 1000$ and (h) $Re = 1200$. Here, hairpin vortices on the upper and side faces are respectively denoted by ‘${\rm HP}_U$’ and ‘${\rm HP}_S$’, separation bubbles are denoted by ‘SB’ and horseshoe vortex and vortex loops are respectively highlighted as ‘HV’ and ‘VL’.

Figure 10

Figure 10. Long prism lift coefficient versus normalized time ($t^* = (t-t_0)U_0/h$) at (a) $Re = 500$ and (b) $Re = 1000$.

Figure 11

Figure 11. (ad) Phase portraits based on the instantaneous lift and drag forces and (eh) PSD of the lift force signal at different Reynolds numbers.

Figure 12

Figure 12. Instantaneous velocity profiles at $x/h = 0.5$ and $x/h = 5$ at an arbitrary time. The black dashed lines show the velocity signal at $-4 \leq y/h \leq 0$ and solid red lines display the profile at $0 \leq y/h \leq 4$.

Figure 13

Figure 13. Instantaneous vorticity distribution on the prism's mid-height plane ($z/h = 0.5$) at different Reynolds numbers: (a) $Re = 625$, (b) $Re = 650$, (c) $Re = 675$, (d) $Re = 700$, (e) $Re = 725$, ( f) $Re = 750$, (g) $Re = 850$, (h) $Re = 1000$ and (i) $Re = 1200$.

Figure 14

Figure 14. Isosurface of $Q^*=0.18$ in the wake of the long rectangular prism at six consecutive time instants at $Re=1000$: (a) $t= t_0$, (b) $t= t_0 + \tau$, (c) $t= t_0 + 2\tau$, (d) $t= t_0 + 3\tau$, (e) $t= t_0 + 4\tau$ and ( f) $t= t_0 + 5\tau$. Hairpin vortices on the upper and side faces are respectively denoted by ‘${\rm HP}_U$’ and ‘${\rm HP}_S$’, and horseshoe vortex and vortex loops are respectively highlighted as ‘${\rm HV}$’ and ‘${\rm VL}$’. Here, $\tau = 1.0417 h/U_0$.

Figure 15

Figure 15. Isosurface of $Q^*=0.18$ from (a,c) side view and (b,d) top view in the wake of the long rectangular (large-depth-ratio) prism at the same arbitrary time for the case of (a,b) $Re=650$ and (c,d) $Re=1000$.

Figure 16

Figure 16. The isosurface of $Q^*=0.18$ overlaid by the time-averaged streamwise velocity and instantaneous vorticity contours on five different planes at $Re = 1000$.

Figure 17

Figure 17. Instantaneous streamwise velocity profiles in the side and top wake, normal to the prism surfaces at $Re = 1000$. The distances are normalized using the prism length ($l = 5w = 4.17h$).

Figure 18

Figure 18. Contour of vorticity along with circulation of main structures on (a) $xy$ plane at $z/h=0.5$ and (b) $xz$ plane at $y/h=0$ for the case of $Re=1000$.

Figure 19

Figure 19. The PSD of velocity fluctuations at three locations in the wake of a long rectangular prism at (a) $Re=750$ and (b) $Re=1000$. Here $St_{hp1}, St_{hp2}, St_{H}$ and $St_{a}$ are captured at $w(-0.4, 0.5, 0.6), u(-0.4, 0, 1.5), u(-2, 1.0, 1.5)$ and $v(5, 0, 0.6)$, respectively.

Figure 20

Figure 20. The PSD of velocity fluctuations at three locations in the wake of a long rectangular prism at (a) $Re=650$ and (b) $Re=700$. Here $St_{hp1}$, $St_{hp2}$, $St_{H}$ and $St_{a}$ are captured at $w(-0.4, 0.5, 0.6)$, $u(-0.4, 0, 1.5)$, $u(-2, 1.0, 1.5)$ and $v(5, 0, 0.6)$, respectively.

Figure 21

Figure 21. Contours of instantaneous vorticity in spanwise ($y$) and normal ($z$) directions at mid-planes ($z/h = 0.5$ and $y/h = 0$) of the long prism at $Re = 1000$ during five consecutive instants: (a) $t = t_0$, (b) $t = t_0 + \tau$, (c) $t = t_0 + 2\tau$, (d) $t = t_0 + 3\tau$ and (e) $t = t_0 + 4\tau$. Here, $\tau = 1.0417 h/U_0$.

Figure 22

Figure 22. Cross-correlation analysis based on the dataset captured at $Re = 700$. Here $f_a$ is the frequency associated with arc-type wake structure with Strouhal number of $St_a = 0.2$.

Figure 23

Figure 23. Time-averaged two-dimensional streamwise streamlines in the near-wake region at (ad) $Re=750$ and (eh) $Re=1000$.

Figure 24

Figure 24. Time-averaged two-dimensional streamwise streamlines in the far-wake region at (ac) $Re=750$ and (df) $Re=1000$.

Figure 25

Figure 25. Time-averaged distribution of normal ($z$) velocity component on $yz$ planes in the near wake at (ad) $Re=750$ and (eh) $Re=1000$.

Figure 26

Figure 26. Time-averaged normal ($z$) component of velocity vector at $y/h=0$: (a) $Re=350$, (b) $Re=500$, (c) $Re=750$ and (d) $Re=1000$.

Figure 27

Figure 27. (b) Time-averaged coefficient of surface pressure variations with Reynolds number on the middle prism height ($z/h = 0.5$) $t$ axis shown schematically in (a).

Figure 28

Figure 28. Variation of the normalized time-averaged circulation with Reynolds number for the two-dimensional vortex structure (a) on the upper face ($y/h=0$) and (b) the side face ($z/h=0.5$). These results are calculated within (a) $-4.16 \leq x/h \leq -2$ and $0.416 \leq y/h \leq 0.7$ and (b) $-4.16 \leq x/h \leq -2$ and $1 \leq z/h \leq 1.35$.

Figure 29

Figure 29. Variations of the time-averaged distance of (a) saddle point on the upper prism face and (b) nodal point on the ground plane with Reynolds number.

Figure 30

Figure 30. (a) Time-averaged $Q$ criterion and time-averaged vorticity distribution plotted on the $x/h = -1$ plane at $Re=1000$. (b) Effect of changing Reynolds number on the strength of the horseshoe vortex parts.

Figure 31

Figure 31. Effect of Reynolds number on the location of horseshoe vortex components in (a) the streamwise direction from the prism leading edge and (b) the normal direction from the ground plane.

Figure 32

Figure 32. The (a,b) spanwise and (c,d) normal velocity profiles at streamwise mid-length of the prism ($x/h = -2.08$).

Figure 33

Figure 33. New model for the vortex structures of a long rectangular prism.