UFR 3-33 Evaluation

From KBwiki
Jump to navigation Jump to search

Front Page


Test Case Studies


Best Practice Advice



Unsteady results

Based on the experimental and numerical unsteady data the flow field is characterized using the systematic classification map of the unsteady flow patterns given by Savory and Toy (1986). Seven regions are highlighted in Fig. 14, which shows the results of the instantaneous flow field at an arbitrary instant in time:

  • (1) Upstream of the hemisphere the horseshoe vortex system dominates. The hemispherical bluff body acts as an obstacle which leads to a positive pressure gradient so that the boundary layer separates from the ground forming the horseshoe vortex system.
  • (2) The stagnation area is located close to the lower front surface of the hemisphere where the stagnation point is found at the surface at an angle of about ≈ 166° (definition of in Fig. 14(a)).
  • (3) The flow is accelerated along the contour of the hemisphere. Therefore, this region is called the acceleration area. The acceleration leads to a high level of vorticity near the hemispherical surface.
  • (4) At an angle of ≈ 90° the flow detaches along a separation line.
  • (5) As a consequence of the flow separation a recirculation area appears. This region is separated from the outer field by a dividing streamline.
  • (6) Strong shear layer vorticity can be observed leading to the production of Kelvin-Helmholtz vortices which travel downstream in the flow field.
  • (7) The extent of the recirculation area behind the hemisphere is characterized by the reattachment point. In this region the flow impinges on the wall and a splatting effect occurs, redistributing momentum from the wall-normal direction to the streamwise and spanwise directions.

UFR3-33 flow regions.png Download movie

Fig. 14: Visualization of flow regions and characteristic flow features of the flow past the hemisphere: (1) horseshoe vortex system, (2) stagnation area, (3) acceleration of the flow, (4) separation point, (5) dividing streamline, (6) shear layer vorticity, (7) reattachment point.

The 3D geometry generates a 3D flow field illustrated in Fig. 15. The complex flow patterns are visualized using iso-surfaces of the pressure fluctuations (p′/(ρair U²) = -2.47 × 10-4) as recommended by Garcia-Villalba et al. (2009):

  • Just upstream of the bluff body the horseshoe vortex system dominates and leads to necklace-vortices that stretch out on both sides into the wake region.
  • The flow detaches from the surface of the hemisphere along the indicated separation line and the vortices roll up. Shedding type and frequency vary in time along the separation line.
  • In the wake these small roll-up vortices interact with the horseshoe vortices just after the hemisphere and in some cases they merge, generating big entangled vortical hairpin-structures. These patterns travel downstream forming a vortex chain. Schematic 3D sketches and explanations of the formation of the hairpin-structures around and behind the hemisphere can be found in the literature (Tamai et al., 1987, Acarlar and Smith, 1987).

UFR3-33 unsteady vortical structures.png

Fig. 15: Snapshot of unsteady vortical structures visualized by utilizing the iso-surfaces of the pressure fluctuations (p′/(ρair U²) = -2.47 × 10-4) colored by the spanwise instantaneous velocity based on LES prediction.

The shape of the instantaneous vortices described above depends on the shedding type. In order to determine the type and frequency of the shedding processes, the velocity spectra are given at the two monitoring points P1 and P2 in Fig. 16(a). Both points are located in the wake (see their position in Figs. 16(b) and (c)) and chosen based on the analysis of Manhart (1998). At these two locations the dominant shedding frequencies are clearly visible. In order to be sure to capture all frequencies of the wake flow and to get smooth velocity spectra, the data have to be collected with an adequate sampling rate during a long time period. Measurements are more appropriate than LES predictions for this purpose. The measurements were carried out with a sampling rate of 1 kHz and were collected over a period of 30 minutes with the hot-film probe described in Section Constant temperature anemometer.

UFR3-33 unsteady velocity spectra.png

Fig. 16: Velocity spectra at the monitoring points P1 and P2 in the wake regime of the hemisphere (experimental data).

Fig. 16(a) provides the power spectral density (PSD) at each location (in blue for P1 and in red for P2). With the help of the PSD maxima the shedding frequencies and the corresponding Strouhal numbers were determined:

  • At P1 on the symmetry plane the PSD is high between 7.9 Hz ≤ f1 ≤ 10.6 Hz (0.23 ≤ St1=f1 D / U ≤ 0.31). A maximum is reached at about f1 = 9.2 Hz (St1 ≈ 0.27).
  • At P2 a distinct frequency peak is found at f2 = 5.5 Hz corresponding to a Strouhal number of St2 ≈ 0.16.

These results suggest the presence of two vortex shedding types in the wake:

  • At the top of the hemisphere the flow detachment generates a chain of arch-type-vortices observed in the symmetry plane at P1 (see Fig. 17) with a shedding frequency in the range 7.9 Hz ≤ f1 ≤ 10.6 Hz.
  • The second type is a von Karman-shedding process occurring at a shedding frequency of f2 = 5.5 Hz on the sides of the hemisphere captured at point P2.

UFR3-33 unsteady vortex shedding top.png Download movie

Fig. 17: Vortex shedding from the top of the hemisphere visualized by the pressure fluctuations of the LES in the symmetry plane (Click on the figure to see the animation).

The second shedding process involves two clearly distinguishable shedding types that switch in shape and time (see Fig. 18):

  • A quasi-symmetric process in which the vortices detach in a double-sided symmetric manner (visualized by the velocity magnitude near the wall in

Fig. 18(a) and schematically depicted in Fig. 18(c)) and form arch-type-vortices (Sakamoto and Arie, 1983) or symmetric-vortices (Okamoto and Sunabashiri, 1992);

  • A more classical quasi-periodic vortex shedding resulting in a single-sided alternating detachment pattern (visualized by the velocity magnitude near the wall in Fig. 18(b) and schematically depicted in Fig. 18(d)).

This alternating behavior was also noted by Manhart (1998). He assumed that the symmetric shedding type is mainly driven by small-scale, less energetic turbulent structures in the flow field. It nearly completely vanished in his predictions when performing a large-eddy simulation on a rather coarse grid, where the small-scale flow structures cannot be resolved appropriately.

UFR3-33 unsteady vortex shedding bottom.png Download movie

Fig. 18: Visualization of the two vortex shedding types present in the wake behind the hemisphere (LES data) (Click on the figure to see the animation).

Comparison of numerical and experimental time-averaged results

Mean velocities in the symmetry plane

This section presents the time-averaged results of the flow around the hemisphere. The instantaneous results were averaged over a time period of about 1370 units of the dimensionless time (normalized by U and D) for the laser-Doppler measurements which ensures a sufficiently large data sample. Due to the small time step, the LES results are averaged over a period of 86 dimensionless time units, following the advice given by Garcia-Villalba et al. (2009). To describe the main features of the flow, two-dimensional color plots of

  • the time-average velocity field and
  • the corresponding Reynolds stresses

are depicted identifying characteristic regions that are mentioned above. Here, the focus is on the symmetry plane and one specific spanwise plane in the wake.

Figure 19 presents the velocity field around the hemisphere focusing on the streamwise and the wall-normal components (LDA measurements (left), large-eddy simulation (right)). The streamwise velocity component u/U is shown in Figs. 19(a) and (b). First, the oncoming flow upstream of the hemisphere in the region -1.5 ≤ x/D ≤ -0.75 is investigated. The experiment shows that the thickness of the approaching boundary layer matches well the height of the hemisphere with z/D ≈ 0.5. A comparable velocity distribution can be seen in the large-eddy simulation. The development of a small recirculation area occurs close to the lower front of the hemisphere between -0.75 ≤ x/D ≤ -0.5 which is connected to the horseshoe vortex system as a result of the reorganization of the approaching boundary layer. The turbulent boundary layer detaches from the ground at = -0.97 due to the positive pressure gradient (stagnation area) located at the bottom front of the hemisphere at about = 166°in the measurements and at about = 161°in the simulation. The extent of the horseshoe vortex depends on the turbulence intensity of the approaching flow. Although the inflow conditions of the synthetic turbulence inflow generator were adjusted to the experimental boundary layer, the horseshoe vortex shows a slightly larger upstream extent in case of the numerical simulation.

UFR3-33 time-averaged results velocity contours.png

Fig. 19: Comparison of the experimental and numerical time-averaged velocity components and streamlines in the symmetry x-z-plane at y/D = 0.

Another interesting location is the separation point on the surface of the hemisphere. It marks an important characteristic for the validation of numerical simulations. Its position depends on various influence parameters (Reynolds number, turbulence intensity of the boundary layer and surface roughness). The flow detaches at an angle of ≈ 90°in case of the laser-Doppler measurements. A comparable angle of ≈ 92°is evaluated for the LES. The separated flow leads to the development of a free shear layer. This phenomenon can be observed as a strong velocity gradient between the outer flow field and the recirculation area in the wake regime. The size of the recirculation area stretches up to x/D ≈ 1.0. It is interrelated to the turbulence intensity of the approaching boundary layer. According to previous studies (Toy et al., 1983; Savory and Toy, 1988; Tavakol et al., 2010; Kharoua, 2013) the turbulence level of the oncoming flow influences the length of the recirculation area. With increasing turbulence intensity the location of the separation point is shifted further downstream on the hemisphere.

The wall-normal velocity component w/U is depicted in Figs. 19(c) and (d). The flow field close to the bottom wall is not resolved in the experimental investigation due to the restrictions of the chosen setup. A notable region is the area of increasing velocity at the front side of the hemisphere at -0.45 ≤ x/D ≤ -0.15 and 0.25 ≤ z/D ≤ 0.45 as a result of the acceleration of the fluid after passing the stagnation area. The size of this area and the velocity magnitude are almost identical between experiment and numerical simulation. Similar characteristics can be detected at about 0.5 ≤ x/D ≤ 1.5 and 0.40 ≤ z/D ≤ 0.85 above the recirculation area. Furthermore, a comparison of the streamline plots of the experiment and the numerical simulation is given in Figs. 19(e) and (f). In conclusion, the overall velocity distributions found in the experiment and the numerical simulation are very similar. A closer view using specific velocity profiles at certain positions within the flow field provides a more detailed insight into the quantitative data.

Figure 20 presents the velocity distribution at specific locations along the symmetry plane for the streamwise (Fig. 20(a)) and the wall-normal (Fig. 20(b)) component. The results of the large-eddy simulation are presented as blue solid lines superimposed by the discrete measurements of the LDA outlined as black squares. The profile at x/D = -0.6 in front of the hemisphere represents the position of the horseshoe vortex system with a strong backflow in the near-wall region that is well predicted by the large-eddy simulation. Another representative position of the flow field is located at x/D ≥ 0.25. The results show an excellent coincidence concerning the developing shear layer and the velocity distribution in the wake.

UFR3-33 time-averaged results velocity lines.png

Fig. 20: Comparison of the experimental (black symbols) and numerical (blue lines) time-averaged streamwise and wall-normal velocity in the symmetry x-z-plane at y/D = 0 and x/D = {-1.5, -1, -0.6, -0.25, 0, 0.25, 0.5, 1, 1.5} (only every second measurement point is displayed)

Reynolds stresses in the symmetry plane

Figures 21(a) and (b) refer to the normal Reynolds stress u′u′/U². The turbulence intensity in the approaching boundary layer is visible in case of the laser-Doppler measurement. In the large-eddy simulation the incoming turbulence intensity is not as strong even through an equal turbulence intensity level is imposed at the STIG window. A certain part of the generated turbulent fluctuations is damped by the numerical discretization scheme (non-equidistant, flux blending includes 5% of a first-order upwind scheme). The highest Reynolds stresses appear in the free shear layer which is connected to the rapid roll-up process of the vortical structures. This area extends into the upper recirculation region with high turbulent mixing rates. The splatting process arising at the reattachment point produces also streamwise fluctuations. However, these are not visible in the figure, because their associated magnitude is much lower than in the shear layer.

UFR3-33 time-averaged results reynolds stresses contours.png

Fig. 21: Comparison of the experimental and numerical time-averaged Reynolds stresses in the symmetry x-z-plane at y/D = 0.

The spanwise normal component v′v′/U² is depicted in Figs. 21(c) and (d). The experiment reveals high Reynolds stresses in the recirculation area as well as for the near-wall region. Note that the spanwise normal component v′v′/U² around the reattachment area is very high. Its values are comparable with the normal Reynolds stress u′u′/U² in the free shear layer. It is assumed that the spanwise velocity fluctuations are associated with the splatting process taking place in the reattachment region and with the detaching vortices at the sides of the hemisphere. Parts of the momentum are redistributed from the wall-normal component to the lateral component. The results of the large-eddy simulation support these observations delivering higher normal Reynolds stresses in the lower wake flow. Nevertheless, the experimental results reveal significantly higher v′v′/U² distributions in the upper part of the recirculation compared with the numerical simulation. The reason for this deviation is presently unclear.

The wall-normal Reynolds stress w′w′/U² is presented in Figs. 21(e) and (f). High Reynolds stresses are present in the free shear layer and the recirculation region at 1 ≤ x/D ≤ 1.5. The Reynolds shear stress u′w′/U² is shown in Figs. 21(g) and (f). Both the measurement and the simulation show that the largest values are expected in the free shear layer.

The profiles of the Reynolds stresses are presented in Fig. 22. The complete upper flow field until x/D = 0 shows only minor differences between the laser-Doppler measurements and the large-eddy simulation for all Reynolds stress components. The streamwise Reynolds stress u′u′/U² is well predicted past the separation point. The results of the Reynolds shear stress u′w′/U² show close agreement between the experiment and the simulation.

UFR3-33 time-averaged results reynolds stresses lines.png

Fig. 22: Comparison of the experimental (black symbols) and numerical (blue lines) time-averaged Reynolds stresses in the symmetry x-z-plane at y/D=0 and x/D = { -1.5, -1, -0.6, -0.25, 0, 0.25, 0.5, 1, 1.5} (only every second measurement point is displayed).

Mean velocity in a spanwise y-z-plane

Figures 23(a) and (b) depict the streamwise velocity component along a chosen position (x/D = 0.5) in the spanwise y-z-plane. An almost symmetric velocity distribution of the streamwise flow component with regard to the symmetry plane at y/D = 0 is found as expected. The shear layer forms an arch-type structure that is related to the roll-up process of the detaching vortices. The recirculation region expands from -0.3 ≤ y/D ≤ 0.3 and 0 ≤ z/D ≤ 0.3. In the near-wall region the centers of the trailing necklace vortices are observed at the position y/D = 0.7. Figures 23(c) and (d) refer to the spanwise velocity component. The velocity distribution in the lower region closely behind the hemisphere is dominated by two counter-rotating vortices that are located symmetrically to the plane y/D = 0. The alternating direction of the velocity component across the spanwise direction indicates the rotation of the vortices in the opposite direction.

UFR3-33 time-averaged results velocity x-plane contours.png

Fig. 23: Comparison of the experimental and numerical time-averaged velocity components in the y-z-plane at x/D = 0.5.

Reynolds stresses in a spanwise y-z-plane

Referring to the velocity distributions for the spanwise plane discussed above, the corresponding Reynolds stresses are depicted in Fig. 24. Figures 24(a) and (b) show the streamwise Reynolds stress component u′u′/U². The minor differences in its size between the experiment and the large-eddy simulation is related to the applied grid resolution. The very fine mesh used in the large-eddy simulation leads to a better resolution of the gradients in the flow field. This can easily be perceived by the Reynolds stress distribution in the shear layer that reveals an overall thinner arch. The near-wall data of the experiment between 0 ≤ z/D ≤ 0.02 are erroneous due to optical reflections of the flat plate that occur in the utilized LDA setup (configuration 1) and are therefore not usable for further flow interpretation.

Finally, a view of the spanwise Reynolds stress distribution is given in Figs. 24(c) and (d) which confirms all significant effects already mentioned for the streamwise case. Additionally, this component has noticeably higher Reynolds stresses located in the region -0.15 ≤ y/D ≤ 0.15 compared with the streamwise Reynolds stresses. This seems to be connected to the two large vortices that form in this region which are connected to a strong spanwise movement. As observed in the case of the symmetry plane, the Reynolds stresses in the experiment are more pronounced. The discrepancies in the spanwise Reynolds stresses between the numerical simulation and the experiment should be examined in further studies to clarify which side of the investigation is causing this deviation.

UFR3-33 time-averaged results reynolds stresses x-plane contours.png

Fig. 24: Comparison of the experimental and numerical time-averaged Reynolds stresses in the y-z-plane at x/D = 0.5.

3D visualization of the time-averaged flow

The major benefit of the large-eddy simulation lies in its high spatio-temporal resolution. This leads to a large amount of flow field information that can be utilized to analyze even smallest flow structures in characteristic regions such as corner eddies. Besides this, large structures, like the horseshoe vortex system, can be explored in detail. A few chosen numerical results of the three-dimensional time-averaged flow field were used to provide a deeper insight into the characteristics of the flow field around the hemisphere.

Figure 25 presents the bottom wall streamlines based on the time-averaged velocity in the x-y-plane including the surface of the hemisphere.

UFR3-33 time-averaged streamlines.png

Fig. 25: Time-averaged streamlines near the bottom wall and on the surface of the hemisphere.

This view is used to examine the separation and reattachment behavior of the flow field:

  • Far upstream of the hemisphere the flow is divided by the separation streamline that wriggles widely around the obstacle and is connected to the separation of the boundary layer from the ground. This phenomenon is also observed by Martinuzzi and Tropea (1993) for the turbulent flow past a wall-mounted cube at Re = 4.3 × 105.
  • The upstream region close to the hemisphere is dominated by the horseshoe vortex system. At certain positions along the symmetry plane an alternating series of saddle and nodal points indicates either a separation or a reattachment of the flow and helps to separate single vortices. The points can be easily detected since the streamlines bundle up at these specific spots. A comparable formation of vortices is noticed for the wall-mounted cube by Martinuzzi and Tropea (1993).
  • In front of the hemisphere after the stagnation point, the flow field accelerates along the surface up to the separation line. This separation line stretches out along the circumference nearly down to the bottom wall. This is a significant difference to the turbulent flow past the axisymmetric bump (Simpson et al., 2002, Byun and Simpson, 2006, Byun and Simpson, 2010, Garcia et al., 2009) at Re = 1.3 × 105 (based on the hill height), where the separation line is shifted to the backside of the 3D hill.
  • Behind the obstacle, a classical recirculation area forms with a reattachment point located in the symmetry plane (see also Fig. 27). In the recirculation area two symmetric spiral flow pattern are present on the ground, which represent the footprint of the arch-type vortical structure. This pattern is also observable for other wall-mounted bluff obstacles such as the cube (Martinuzzi and Tropea, 1993) and the finite-height circular cylinder (Pattenden et al., 2005). However, such a structure is not mentioned for the axisymmetric 3D hill in Simpson et al. (2002), Byun and Simpson (2006), Byun and Simpson (2010) and Garcia et al. (2009).

One of the characteristic regions is the horseshoe vortex system. Figure 26 depicts the streamlines in the symmetry plane just in front of the lower base of the hemisphere to highlight the vortices forming the horseshoe vortex system. These vortices are paired. The vortices of a pair rotate in opposite directions. Each vortex of the system is classified by applying the indices chosen by Baker (1980) who, among others, presented the mechanism of the formation of the complex vortex system. Vortex 0 is related to the separation of the boundary layer just in front of the lower face of the hemisphere. The counter-rotating vortex 0' is also resolved in the LES but too small to be visualized in Fig. 26. It is followed by the primary vortex 1 which is the largest structure of the horseshoe vortex system. It is caused by the separation of the boundary layer due to the presence of vortex 0. A secondary counter-rotating vortex 1' is generated by the detaching boundary layer beneath the primary vortex. This complex separation process leads to an overall number of four visible vortices.

UFR3-33 streamstraces zoom horseshoe vortices stig.png

Fig. 26: Horseshoe vortex system: Time-averaged streamlines in the symmetry plane upstream of the hemisphere.

The second major flow region is the recirculation area visualized in Fig. 27. The illustrated 3D-filaments in the near wake are coiled up in an arc-shaped structure that characterizes the size of the recirculation area. Corresponding wall streamlines of the rear side of the hemisphere are given in Fig. 28. The streamlines depict the large backflow area on the surface of the hemispherical body. Two other symmetric spiral flow patterns appear on each side relating to the lateral separation points.

UFR3-33 visu recirculation rod stig.png

Fig. 27: Recirculation area: 3D-filaments based on the time-averaged velocity and colored by the mean wall-normal velocity.

UFR3-33 streamstraces behind hemisphere stig.png

Fig. 28: Time-averaged streamlines on the surface behind the hemisphere.

More details of the time-averaged flow can be found in Wood et al. (2016).

Contributed by: Jens Nikolas Wood, Guillaume De Nayer, Stephan Schmidt, Michael Breuer — Helmut-Schmidt Universität Hamburg

Front Page


Test Case Studies


Best Practice Advice


© copyright ERCOFTAC 2022