UFR 3-33 Test Case
- 1 Test case study
- 1.1 Description of the problem
- 1.2 Description of the wind tunnel and of the test section
- 1.3 Measuring techniques
- 1.4 Generation of artificial turbulent boundary layer in the wind tunnel
- 1.5 Numerical simulation methodology
- 2 Data files for download
Test case study
Description of the problem
The investigated hemisphere (diameter D) is rigid and mounted on a smooth wall as depicted in Fig. 3. The hemispherical surface is also considered to be smooth. The structure is placed in a thick turbulent boundary layer whose velocity profile can be described by a 1/7 power law as reviewed by Counihan (1975). At a distance of 1.5 diameters upstream of the bluff body the thickness of the boundary layer δ corresponds to the height of the hemisphere, i.e., δ = D/2. The Reynolds number of the air flow (ρair = 1.225 kg/m³, μair = 18.27 × 10-6 kg/(m s) at ϑ = 20°C) was set to Re = ρair D U∞ / μair ≈ 50,000. The undisturbed free-stream mean velocity in x-direction outside the boundary layer at standard atmospheric conditions is U∞ = 5.14 m/s. The Mach number is low (Ma ≤ 0.3). Therefore, the air flow can be assumed to be incompressible. Moreover, the fluid is considered to be isothermal.
As shown in Fig. 3, the origin of the frame of reference is taken at the center of the base area of the hemisphere, where x denotes the streamwise, y the spanwise and z the vertical direction (wall-normal).
Fig. 3: Geometrical configuration of the wall-mounted hemisphere.
Description of the wind tunnel and of the test section
A Göttingen-type subsonic wind tunnel with an open test section was used for the experimental investigations. Its size and specifications are summarized in Fig. 4 and in the table below.
Fig. 4: Wind tunnel used in the experimental investigations.
In the empty test section the free-stream streamwise turbulence intensity Tu0,u = (u′)rms/U∞ is less than 0.2%. This value is based on high resolution single-wire constant temperature anemometry (CTA) measurements (see Section Constant temperature anemometer) that were conducted to specify the overall quality of the wind tunnel. According to LDA measurements (see Section Laser-Doppler anemometer) the free-stream turbulence intensities in both cross-stream directions are one order of magnitude smaller, which leads to a free-stream total turbulence level of:
A schematic illustration of the symmetry plane of the section with the model is given in Fig. 5(a) (The dimensions are relative to the diameter of the hemisphere D). The hemispherical model made of an aluminum alloy (D = 150 mm, average roughness Ra < 0.8 μm) was mounted on a flat plate with the following specifications:
- To correctly transfer the boundary layer from the nozzle, the leading edge of the plate is in alignment with the bottom of the rectangular nozzle of the wind tunnel.
- Moreover, to ensure a smooth transition of the near-wall flow from the nozzle, the flat plate covers the complete spanwise extension of the cross-section of the wind tunnel.
- A gap remains between the trailing edge of the flat plate and the receiver.
Fig. 5: Dimensions and position of the hemisphere in the test section.
The blocking ratio of the hemispherical structure in the wind tunnel is approximately 4.7% based on the projected area Ahemi = (π/8) D² of the hemisphere and the area of the cross-section of the nozzle. More details on the geometrical setup can be found in Wood et al. (2016).
To measure the flow field in a non-invasive way, a standard Laser-Doppler anemometer (LDA) system was used. It consists of the following components:
- The laser beam is generated by a water-cooled Coherent Innova 70C argon-ion laser.
- The beam is guided to a Dantec FiberFlow transmitter box, where it is split up into two different wave lengths (514.5 nm and 488 nm) which are used for each velocity component.
- The beams are sent from the transmitter box to a 2D optical probe.
- The aerosol generator TSI Six-Jet Atomizer 9306 is utilized to generate small droplets of Di-Ethyl-Hexyl-Sebacat (DEHS). The spherical droplets have an average size between 0.2 and 0.3 μm and feature a long-life cycle, which leads to excellent optical properties for LDA and other reflection based measurement techniques. They are atomized and seeded at the receiver of the wind tunnel.
- The Dantec BSA Flow Software is used for the measurements and for the positioning of the traverse.
In order to measure all three velocity components with a two-component probe the LDA system is configured in two ways depicted in Fig. 6:
- Configuration 1: The probe is facing the flat plate so that the optical axis of the laser beam points vertically at the surface of the plate (see Fig. 6(a)). This setup is primarily used to record the spanwise component v. Strong reflections occur close to the surface region and measurements are not valid in this region.
- Configuration 2: The probe is on the side of the test section so that the optical axis of the laser beam points horizontally (see Fig. 6(b)). This setup acquires the streamwise and wall-normal flow components u and w. With this horizontal configuration unwanted laser reflections of the surface are avoided leading to a better spatial resolution of the streamwise component. However, note that the measurement of the wall-normal component w is geometrically restricted (see Wood et al., 2016 for details).
Fig. 6: LDA configuration and measurement grid resolution of the symmetry x-z cross-section.
In Fig. 6(a) the dashed area in blue indicates the x-z-measurement plane. The associated experimental grid used for the LDA measurements is shown in Fig. 6(c). It contains 1239 measurement points. Their distribution was chosen based on the characteristics of the predicted flow field. In the near-wall region and the recirculation area close to the hemisphere more measurement points are placed to ensure a good resolution of the gradients.
Constant temperature anemometer
A standard constant temperature anemometer (TSI 1750 CTA) system was used to complete and expand the measurements. The CTA technique was applied to record the velocity spectra in the wake of the hemisphere (see Section Unsteady results). The spectra provide information about the vortex shedding processes and their frequencies as well as on the decay in the inertial subrange. The main reason to use CTA for this specific analysis instead of LDA is that raw data samples from LDA measurements show an irregular time signature since the measured droplets pass the measurement volume randomly. Furthermore, the CTA measurements were used to validate the LDA data of the artificial turbulent boundary layer (see Section Generation of artificial turbulent boundary layer in the wind tunnel).
Generation of artificial turbulent boundary layer in the wind tunnel
To adjust the inlet conditions described in Section Description of the problem in the experiment, the oncoming experimental boundary layer has to be artificially adjusted. The test section of a wind tunnel is not long enough to generate naturally developing turbulent boundary layers of a desired thickness (δ/(D/2) = 1 at a distance of x/D = -1.5 upstream of the hemisphere). Therefore, in order to rapidly thicken the boundary layer, turbulence generators were installed at the bottom wall inside the nozzle (see Fig. 7). The combination of the specific turbulence generators is the result of a preliminary study detailed in Wood et al. (2016).
The size, the form and the position of the turbulence generators are depicted for the symmetry plane in Fig. 8. Near the outlet the first vortex generator is a fence with rods shown in green. It is followed by a castellated barrier marked in blue. Both devices are used to generate large disturbances in the near-wall flow. The produced vortical structures travel to the third vortex generator illustrated in red, which is a fence with triangular spikes. It is used to break up the large structures to make turbulence more isotropic in the boundary layer.
Fig. 7: Generation of a thick turbulent boundary layer with turbulence generators mounted onto the bottom wall of the wind tunnel's nozzle.
Fig. 8: Close view on the position of the vortex generators inside the nozzle.
In order to ensure that this artificially thickened turbulent boundary layer corresponds to the predefined conditions, a combination of LDA and CTA measurements of the inlet boundary layer profile were carried out and are presented in Fig. 9. All data were recorded just after the beginning of the test section at x/D = -1.5 in the symmetry plane without placing the hemisphere into the test section.
The measured mean velocity distributions (LDA and CTA) u/U∞ are in close agreement with the reference (1/7 power law) exhibiting minor deviations in the region z/D = 0.25 (see Fig. 9(a)). The free-stream velocity u/U∞ ≈ 1 is reached at the edge of the boundary layer which has the desired thickness. Figure 9(b) presents the dimensionless velocity u+ plotted against the dimensionless wall-normal distance z+. The first measured point is located at a distance of Δz = 0.25 mm above the flat plate. The velocity distribution is nearly linear in the region 4 ≤ z+ ≤ 10. Therefore, the first two points close to the wall are still inside the viscous sublayer. This is in good agreement with the literature where it is often stated that the border between the viscous sublayer and the buffer layer is at about z+≈ 5. The friction velocity is estimated to uτ = = 0.225 m/s (uτ / U∞ = 4.38 × 10-2), where τw is approximated by μair Δu / Δz. Both classical laws of the wall (viscous sublayer: u+ = z+ and log-layer: u+ = 1/0.41 (z+) + 5.2) are correctly reproduced. Some discrepancies in the velocity distribution are observed in the log area. Indeed, the measured boundary layer does not exactly follow the 1/7 power law. Note that the mean velocity varies only little across the spanwise extent of the test section (not shown here for the sake of brevity).
Additionally, the normalized displacement thickness δ1 / δ and the momentum thickness δ2 / δ are evaluated from the experimental data to 1/8 and 7/72, respectively. This results in a shape factor of H = δ1/δ2 = 1.286 which confirms a classical property of a turbulent boundary layer. The Reynolds number based on δ2 is estimated as Reδ2 = 2503.
Figures 9(c) and 9(d) show the Reynolds stresses of the produced turbulent boundary layer. The streamwise fluctuations (u′)rms/U∞ measured by CTA and LDA are very similar (see Fig. 9(c)). Therefore, it is assumed that the other components which are solely measured by LDA, are also reliable. This is supported by the fact that Counihan (1969) and Schlatter et al. (2009) published similar distributions of turbulence intensities of a turbulent boundary layer flow on a flat plate. In the present case, the streamwise turbulent fluctuations (u′)rms/U∞ gradually increase from the free-stream (1.2%) to the near-wall region with a peak value of about 12.1% at z/D = 0.01. Contrary to the case of an empty test section without any turbulence generator presented in Section Description of the wind tunnel and of the test section, the free-stream streamwise turbulence intensity for the generated thick boundary layer (TBL) is significantly higher, i.e., TuTBL,u = (u′)rms/U∞ = 1.2%. Consequently, the total free-stream turbulence level also increases to:
Figure 9(d) presents the Reynolds shear stress u′w′/U∞², which is the only one of interest for the flow physics of the current test case. The other cross-components theoretically vanish due to the homogeneity of the spanwise direction. As expected for a positive mean velocity gradient at the wall, the Reynolds shear stress is negative.
Fig. 9: Inflow properties of the turbulent boundary layer at the inlet of the test section.
In conclusion, the predefined inlet conditions of a fully developed turbulent boundary layer with a desired thickness of δ/D ≈ 0.5 have been achieved by using customized vortex generators in the nozzle of the wind tunnel. The obtained data are verified by comparison with the 1/7 power law as the reference for the velocity distribution. Furthermore, an overall reasonable distribution of the turbulent fluctuations similar to a natural boundary layer commonly presented in the literature is achieved.
These experimental inflow profiles are available for download in the Section Experimental data/Inflow profiles.
Numerical simulation methodology
To predict the turbulent flow around the hemisphere based on the large-eddy simulation technique, the 3D finite-volume fluid solver FASTEST-3D is used. This in-house code is an enhanced version of the original one (Durst and Schäfer, 1996, Durst et al. 1996). To solve the filtered Navier-Stokes equations for LES, the solver relies on a predictor-corrector scheme (projection method) of second-order accuracy in space and time (Breuer et al., 2012). The discretization is based on a curvilinear, block-structured body-fitted grid with a collocated variable arrangement. The surface and volume integrals are calculated applying the midpoint rule. Most flow variables are linearly interpolated to the cell faces leading to a second-order accurate central scheme. The convective fluxes are approximated by the technique of flux blending (Khosla and Rubin, 1974, Ferziger and Peric, 2002) to stabilize the simulation.For the current case the flux blending includes 5 % of a first-order accurate upwind scheme and 95 % of a second-order accurate central scheme. A preliminary study shows that these settings are a good compromise between accuracy and stability. The momentum interpolation technique of Rhie and Chow (1983) is applied to couple the pressure and the velocity fields on non-staggered grids.
FASTEST-3D is efficiently parallelized based on the domain decomposition technique relying on the Message-Passing-Interface (MPI). Non-blocking MPI communications are used and offer a non negligible speed-up compared to blocking MPI communications (Scheit et al. 2014).
Subgrid-scale modeling in LES
Since LES is used, the large scales of the turbulent flow field are resolved directly, whereas the non-resolvable small scales have to be taken into account by a subgrid-scale (SGS) model. Different SGS models based on the eddy-viscosity concept are available in FASTEST-3D: The well-known and most often used Smagorinsky model (Smagorinsky, 1963), the dynamic Smagorinsky model according to Germano et al. (1991) and Lilly (1992), and the WALE model (Nicoud and Ducros, 1999). Owing to the moderate Reynolds number considered and the fine grid applied, the SGS model is expected to have a limited influence on the results. Nevertheless, in order to investigate and verify this issue, simulations of the flow around the hemisphere are carried out applying the above mentioned SGS models. The results are presented and analyzed in Wood et al. (2016). This SGS investigation shows that the Smagorinsky model with 0.065 ≤ Cs ≤ 0.1 or the dynamic Smagorinsky model basically lead to the same results. The WALE model with CW = 0.33 (value corresponding to the classical Smagorinsky model with Cs = 0.1 (Nicoud and Ducros, 1999) produces a nearly identical flow except for the region upstream to the hemisphere. Therefore, as the best compromise between accurate results and fast computations, the standard Smagorinsky model with Van Driest damping near walls and the constant set to Cs = 0.1 is used for the present case.
Computational domain and grid
To simulate the problem using a block-structured mesh, the chosen computational domain is a large hemispherical expansion with its origin at the center of the hemisphere (see Fig. 10(a)). This domain is originally divided into 5 geometrical blocks, so that nearly orthogonal angles are obtained on the surface of the hemisphere (see Fig. 10(b)) and in the entire volume. To prescribe the inlet and outlet boundary conditions described in the next paragraph, the upper, left and right blocks are divided along the x/D = 0 plane leading to 8 geometrical blocks (see Fig. 10(a)). Figure 10(c) shows the x-y cross-section of the grid at the bottom wall and Fig. 10(d) depicts the x-z cross-section in the symmetry plane. For the sake of visualization only every fourth grid line of the mesh is shown. The 8 geometrical blocks are later split into 80 parallel blocks for the distribution of the computation on a parallel computer. The outer domain has a radius of 10D. 240 grid points are distributed non-equidistantly based on a geometrical stretching in the expansion direction. 640 points are used at the circumference of the bottom of the hemisphere. The final grid contains 30.72 × 106 control volumes (CVs). In order to fully resolve the viscous sublayer, the first cell center is located at a distance of Δz/D ≈ 5 × 10-5 from the wall, which leads to averaged z+ values below 0.25 (see Figs. 10(e) and (f) and more than 50 points in the boundary layer on the hemisphere upstream to the separation. The geometrical stretching ratios are kept below 1.05. The aspect ratio of the cells on the hemispherical body are low, i.e., in the range between 1 and 10. This yields a dimensionless cell size in the two tangential directions below 29, which fits to the recommendation of Piomelli and Chasnov (1996) for a wall-resolved LES. Note that the resolution of the grid is chosen based on extensive preliminary tests not presented here. For this fine grid a small time step of Δt* = Δt U∞ / D = 3.084 × 10-5 is required ensuring a CFL-number below unity.
Fig. 10: Grid used for the LES predictions including the resolution at specific locations.
The boundary conditions used in the simulation are listed below and depicted in color in Fig. 11: Black for the walls, blue for the inlet and red for the outlet.
- At the bottom of the domain and on the hemisphere a no-slip wall condition is applied justified by the fine near-wall resolution mentioned above.
- A 1/7 power law with δ/D = 0.5 and without any perturbation is applied as inlet condition on the external surface of the domain for x ≤ 0. Moreover, this power law is applied for all CVs with x/D ≤ -2 (see the area with hatched lines on Fig. 11). This region (x/D ≤ -2) does not need to be solved for the problem. However, it could not be simply cut from the mesh because of the hemispherical form of the block-structured grid. Therefore, for all CVs with x/D ≤ -2 the flow field is not predicted, so that the mean velocity profile at x/D = -2 remains constant in time and perfectly fits the experiment. In order to approximate the turbulent boundary layer depicted in Fig. 9, perturbations produced by a turbulence inflow generator (described in Section Generation of artificial turbulent boundary layer in the wind tunnel) are injected in a 2D × D window at x/D = -1.5 (see Fig. 11 (b)).
- A zero velocity gradient boundary condition is defined for the outlet on the external surface of the domain for the geometrical blocks 5, 6 and 7 as defined in Fig. 10 (a). At the outlet of block 8 where the large-scale flow structures leave the computational domain, a convective boundary condition is applied with a convective velocity set according to the 1/7 power law. The fact that the simulation does not use symmetry boundary conditions or slip walls at the top or at the lateral sides, is in agreement with the free flow situation in the experiments. Indeed, the test section is open on the top and on the lateral sides.
Fig. 11: Boundary conditions.
Synthetic turbulence inflow generator (STIG)
As shown by several research groups (see, e.g., Druault et al., 2004, Schmidt and Breuer, 2015 or Wood et al., 2016) the application of appropriate temporally and spatially correlated velocity distributions as inlet boundary conditions is essential to predict realistic flow fields based on numerical simulations. In the present study the digital filter concept of Klein et al. (2003) is applied to provide instantaneous three-dimensional velocity distributions based on the definition of an integral time scale and two integral length scales.
Considering the typical coarse resolution of numerical grids at the inflow region, the application is limited due to possible damping effects of the numerical scheme on coarse grids. In order to avoid this loss of information, a shift towards finer resolved areas within the computational domain is meaningful. This methodology requires a source term formulation based on the artificial turbulence. In the study of Wood et al. (2016) the following formulation of the source term SΦ,syn is applied: The first expression of the source term denoted as the original formulation inserts the synthetic velocity fluctuations as a temporal derivative as follows:
Here, ρ describes the constant density, (Φ′)syn (= (u′)syn, (v′)syn and (w′)syn) the synthetically generated fluctuations and V the volume of the control volume.
To reduce the required development length of the flow and to avoid discontinuities in the conservation equations, the source terms are superimposed in a spatially spread influence area. Within this influence area the amplitude of the source terms follows the Gaussian shape used within the generation process of the synthetic turbulence. For the present case, it was found that two control volumes in both upstream and downstream direction provide a suitable influence area for the superposition of the fluctuations.
In the current investigation the required information of the integral time and length scales are assumed with the help of the 1/7 power law and the analytic expression of Prandtl's mixing length at z+ = 100 leading to L/D = 2.06 × 10-2. Based on the Taylor hypothesis and the characteristic velocity of the 1/7 power law at z+ = 100, the integral length scale is transformed into the integral time scale T = L U∞/ (D u(z+ = 100)) = 2.79 × 10-2. The integral scale in wall-normal direction and spanwise direction is assumed to be equal to the integral length scale predicted by the formulation by Prandtl (Ly/D = Lz/D = 2.06 × 10-2). As depicted above in Section Boundary conditions the synthetic fluctuations are introduced at x/D = -1.5. Based on the employed dimension of the STIG-window (-1 ≤ y/D ≤ 1 and 0 ≤ z/D ≤ 1), the applied equidistant grid resolution (Δy/D = 9.713 × 10-3 and Δz/D = 2.191 × 10-3) and the the normalized time step of the simulation (Δt U∞/ D = 3.084 × 10-5) the resulting normalized scales and support of the filter for the generation of the synthetic turbulence are: nt × ny × nz = 906 × 2 × 9 and Nt × Ny × Nz = 1812 × 4 × 18. The required reference velocity distributions are given by the experimentally measured velocity profile of the u-component (see Fig. 9(a) and Section Experimental data/Inflow profiles). Furthermore, the measured normal and shear components of the Reynolds stress tensor (see Fig. 9(c) and 9(d) and Section Experimental data/Inflow profiles) are used as reference input data. The missing parts of the Reynolds stress components close to the wall are filled up (see Fig. 12) with the help of the DNS data by Schlatter et al. (2009). In order to describe a two-dimensional turbulent boundary layer, the values of the other shear components u′v′ and v′w′ are set to zero. The inflow profiles used as reference for the generation of the STIG data can be downloaded below (see Section Numerical data/Inflow profiles).
Fig. 12: Distributions of the Reynolds stresses derived from the experiment as input parameters for the application of the synthetic turbulence inflow generator.
160,000 time steps are generated to provide a sufficiently long time series of inflow data. Figure 13 gives an example of contours of the instantaneous streamwise velocity at one arbitrarily chosen time step.
Fig. 13: Example of the instantaneous streamwise velocity distribution generated for the given STIG window.
Data files for download
The dimensionless experimental and numerical data of the present test case are provided here for download.
The experimental data files below contain the time-averaged artificial turbulent boundary layer measured by LDA and CTA at the inflow plane (x/D = -1.5). The velocity u/U∞ is normalized by the free-stream velocity, whereas the Reynolds stresses (u′u′/U∞², v′v′/U∞², w′w′/U∞², u′w′/U∞²) are normalized by the square of the free-stream velocity.
The experimental data files below contain the flow results obtained by the LDA setup presented before.
- The file "UFR3-33_LDA_results_yD0_symmetry_plane.dat" contains the dimensionless results in the symmetry x-z-plane at y/D = 0 for -1.5 ≤ x/D ≤ 2 and 0 ≤ z/D ≤ 1. It has 8 columns: The two first ones contain the x- and z-position of each monitoring point. The next two columns contain the streamwise and wall-normal mean velocities (u/U∞, w/U∞) at the point. The last four columns contain different Reynolds stresses (u′u′/U∞², v′v′/U∞², w′w′/U∞², u′w′/U∞²). All details can be found in the headers of the files.
- The file "UFR3-33_LDA_results_xD05_spanwise_plane.dat" contains the dimensionless results in the spanwise y-z-plane at x/D = 0.5 for -1 ≤ y/D ≤ 1 and 0 ≤ z/D ≤ 1. It has 6 columns: The two first ones contain the y- and z-positions of each monitoring point. The next two columns contain the streamwise and wall-normal mean velocities (u/U∞, v/U∞) at the point. The last two columns contain different Reynolds stresses (u′u′/U∞², v′v′/U∞²). All details can be found in the headers of the files.
The numerical data file below contains the inflow profiles used as reference for the generation of the STIG data. As for the experimental inflow profiles the velocity u/U∞ is normalized by the free-stream velocity, whereas the Reynolds stresses (u′u′/U∞², v′v′/U∞², w′w′/U∞², u′w′/U∞²) are normalized by the square of the free-stream velocity.
The numerical data files below contain the flow results obtained by LES.
- The file "UFR3-33_LES_results_yD0_symmetry_plane.dat (zipped version since it is rather large)" contains the LES raw dimensionless results of the symmetry x-z-plane at y/D = 0 for -1.5 ≤ x/D ≤ 2 and 0 ≤ z/D ≤ 1. It is a Tecplot file and has 8 columns: The two first ones contain the x- and z-positions of each monitoring point. The next two columns contain the mean velocities (u/U∞, w/U∞) at the point. The last four columns contain different Reynolds stresses (u′u′/U∞², v′v′/U∞², w′w′/U∞², u′w′/U∞²). All details can be found in the headers of the file.
- The file "UFR3-33_interpolated_LES_results_yD0_symmetry_plane.dat (zipped version)" contains the dimensionless LES results of the symmetry x-z-plane at y/D = 0 interpolated on a Cartesian grid (400 x 115) with -1.5 ≤ x/D ≤ 2 and 0 ≤ z/D ≤ 1. It has 8 columns: The two first ones contain the x- and z-positions of each monitoring point. The next two columns contain the mean velocities (u/U∞, w/U∞) at the point. The last four columns contain different Reynolds stresses (u′u′/U∞², v′v′/U∞², w′w′/U∞², u′w′/U∞²). Again all details can be found in the headers of the file.
- The file "UFR3-33_interpolated_LES_results_xD05_spanwise_plane.dat (zipped version)" contain the dimensionless KES results in the spanwise y-z-plane at x/D = 0.5 interpolated on a Cartesian grid (230 x 115) with -1 ≤ y/D ≤ 1 and 0 ≤ z/D ≤ 1. It has 6 columns: The two first ones contain the y- and z-positions of each monitoring point. The next two columns contain the mean velocities (u/U∞, v/U∞) at the point. The last two columns contain different Reynolds stresses (u′u′/U∞², v′v′/U∞²). Again all details can be found in the headers of the file.
Contributed by: Jens Nikolas Wood, Guillaume De Nayer, Stephan Schmidt, Michael Breuer — Helmut-Schmidt Universität Hamburg
© copyright ERCOFTAC 2021