Characteristics of low Reynolds number shear-free turbulence at an impermeable base.
The interaction of turbulence with an impermeable surface is important in various geophysical phenomena, notably in the turbulence-sediment interaction including incipient sediment motion and sediment entrainment. An in-depth understanding of the physics of near-wall turbulence is essential to better understand the phenomena. Considerable research has been conducted to investigate the characteristics of turbulence at near impermeable surface, where the turbulence is generated using an oscillating grid in a water tank [1-3]. The grid-generated flow (without the presence of an impermeable boundary) is a statistically stationary, zero-mean turbulence that is horizontally homogeneous at depths beyond 2 to 3 mesh widths (M) below the grid's mid position. At this distance, the decay of the time-averaged r.m.s velocity components (m, v, w), with depth z below the grid's mid position, can be represented by the expression
u = v = [C.sub.1][f.sub.g][S.sup.3/2][M.sup.1/2] [z.sup.-1], (1a)
w = [C.sub.2][f.sub.g][S.sup.3/2][M.sup.1/2][z.sup.-1], (1b)
where [f.sub.g] is the frequency of oscillation, S is the stroke (twice the amplitude), and [C.sub.1] and [C.sub.2] are the horizontal and vertical constants, usually taken as 0.25 and 0.27, respectively . The derivation of (1a) and (1b) requires the assumption that the horizontal integral length scale of the turbulence l satisfies
l = [beta]z, (2)
where [beta] is constant and usually takes the value 0.1 .
Thus, the level of turbulence intensity can be described in terms of Reynolds number [Re.sub.l] = 2ul/v, where v is the water kinematic viscosity. Note that within the quasi-isotropic, quasihomogeneous region at z [??] 2.5M, due to the decreasing u and increasing l, respectively, the Reynolds number [Re.sub.l] remains constant.
When an initially quasi-isotropic turbulence field is made to interact with a solid boundary, the flow at the near impermeable surface region responds in two ways. First, the decreasing normal velocity component due to the kinematic boundary condition and the intercomponent energy transfer results in an amplification of the tangential velocity components. This wall blocking effect occurs within the boundary layer known as a blockage layer, where the thickness depends on the integral length scale l. Second, the no slip condition cancels all velocity components within the thinner viscous layer adjacent to the solid wall. These profile changes has been theoretically described by Hunt and Graham  who utilised the methods of rapid distortion theory. Their model describes the inviscid process of an instantaneous response of a turbulent shear-free boundary layer when an impermeable surface is inserted. The theory predicts that the tangential turbulent intensities increase at the expense of the normal components at near surface. At high [Re.sub.l], this standard profile of the fluctuating field in the bounded shear-free flow has been experimentally and numerically validated [3, 7, 8].
Utilising grid-generated turbulence has been proven to be a successful alternative to examine the sole effect of turbulent fluctuations in various geophysical phenomena such as incipient sediment motion, turbulence-particle effect in suspension, and sedimentation of particles in a water column [9-11]. The statistically stationary grid-generated flow is easy to measure and offers a readily repeatable flow condition, allowing the method to be exploited. Investigating the effect of turbulence on such phenomena, the level of turbulence intensity needed spans from low to high [Re.sub.l]. At low [Re.sub.l], however, Perot and Moin , in their simulation work, conjectured that there is no amplification of the tangential components. In such flow, the properties of statistical turbulence at near impermeable surface deviate from the predicted profile. As such, any geophysical phenomena investigated under (very) low [Re.sub.l] require a different interpretation of the flow behaviour. Thus, the minimum [Re.sub.l] needs to be established when turbulent flow used to investigate the geophysical phenomena is generated by an oscillating grid.
Experimental analyses by McDougall  and Fernando et al.  suggest that the u amplification within the near solid region does not occur for flow within the range 40 < [Re.sub.l] < 150, and this is served as an initial basis. This range, however, is a preconception and cannot as yet be asserted with confidence. The experiments reported here address this preconception by investigating the turbulence properties including the evolution of r.m.s fluid velocity, integral length scales, isotropy, and frequency spectrum; in particular the impermeable surface induced inhomogeneity for flow [Re.sub.l] < 150.
This paper is set out as follows. In Section 2 of this paper the experimental method is described, and the statistical properties of base-induced turbulence at high [Re.sub.l] (to verify the experimental setup adopted) are also presented. Section 3 is devoted to presenting turbulence measurements at low [Re.sub.l] and comparing them with both existing theoretical equation and turbulence profiles at high [Re.sub.l]. A summary of the results is given in Section 4.
2.1. Apparatus. A set of experimental data was performed to investigate the turbulence profile within the bulk flow, focusing on the region at near base. A sketch of the apparatus is shown in Figure 1. The experiments were performed in an acrylic tank, with cross-sectional dimensions of 35.4 x 35.4 [cm.sup.2] and a height of 50 cm. The tank was fixed within a rigid inner steel frame and filled with water to a depth of 40 cm. The oscillating-grid mechanism, consisting of a linear bearing, was positioned vertically through the central axis of the tank and connected to a rotating spindle. A motor was used to rotate the spindle and thus continuously oscillate the grid. The plan view of the square (d = 1 cm) uniform grid shown in Figure 1(b) consisted of a 7 x 7 array of aluminium bars with a length of 35 cm. The uniform mesh size M was 5 cm giving the grid a solidity of 36%, within the acceptable solidity limit [4,14]. The end condition of the grid was designed to ensure that the wall acted as a plane of reflection-symmetry, which has been shown to lessen the Reynolds-stress gradients within the fluid, thereby inhibiting the presence of undesirable secondary circulations (see Fernando and de Silva  for a detailed discussion). The opening between the (side) wall and the grid is approximately 2 mm, which was too small to cause any significant secondary circulation . For consistency, each experiment was performed with the stroke fixed at S = 8 cm with the grid lowered down so that the distance between the midplane of grid and the base tank is z = 3.2M =16 cm. As the stroke was fixed throughout, the turbulence intensity was varied by adjusting the frequency of oscillation within the range 0.8 < f [less than or equal to] 3.1 Hz.
For each experiment, the velocity measurements (between the grid and the tank base) were obtained using two-dimensional planar particle image velocimetry (PIV), applied in the vertical (x, z)-plane through the centre of the tank (y = 0). The reflective, neutrally buoyant tracer particles were illuminated using a narrow, vertical light sheet (with a thickness of 2 mm) that was directed through the midplane of the tank (Section A-A as shown in Figure 1(b)). The motion of the tracer particles was recorded using a high-speed Dantec NanoSense MkIII digital camera (with 1280 x 1024 pixel resolution) and the frame rate was set at 200 Hz. The camera was positioned so as to have a horizontal view through the tank sidewall and be normal to the vertical light sheet. The camera was connected to a computer, which was used to acquire and store the images. In each case, the instantaneous velocity fields u = ([??], [??]) were calculated from pairs of images with a cross-correlation algorithm using the Digiflow PIV software . The interrogation area was set at 20 x 20 pixels (equivalent to 3.8 x 3.8 [mm.sup.2]), with a 50% overlap between the analysis spots. For each experimental setup, the velocity fields were acquired for a period T = 120 seconds, which provided the spatial and the temporal data of the velocity motions in the bulk fluid. Post processing of the images showed that mean and r.m.s velocity components were essentially converged after time averaging over a period of at least 80 seconds. The measurement of velocity using PIV allowed a planar flow distribution to be obtained, from the bottom of the grid (i.e., z = S/2 + d/2) to the tank base, giving a physical region of approximately 160 x 200 [mm.sup.2]. Thus, the horizontal mean U(x, z) and r.m.s u(x, z) velocity components were obtained by temporally averaging over the acquired period. Note that similar procedure was applied for the mean and r.m.s vertical velocity components.
2.2. The Statistical Properties of Base-Induced Inhomogeneous Turbulence. The inhomogeneity behaviour within the blockage layer has been studied extensively using experimental apparatus similar to those described above and by numerical modelling [1, 3, 13]. This section presents the results of a series of experiments, that is, flow with high Re;, to investigate and compare the basic statistical properties of the turbulence generated within the previous experimental studies and the theoretical prediction of Hunt and Graham . However, prior to further analysis, it is practical to examine the significance of the mean flow occurrence in the generated grid turbulence.
The weak, secondary mean flow is an inherent characteristic of an oscillating-grid turbulence experiment . Therefore, although several design restrictions have been met (i.e., grid solidity less than 40%, oscillation frequency less than 6 Hz, and very small end openings), it is important to assess the turbulence intensities, in particular at near-base to ensure that the effect of the mean flow is insignificant. Using the operator <*> to denote spatial averaging, the r.m.s horizontal velocity components were spatially averaged across the horizontal direction x, defined as
<u(z)> = 1/2L [[integral].sup.L.sub.-L] u(x, z) dx, (3)
where L is the horizontal width from the central midplane of the tank, with typical value of 2L = 200 mm. Figure 2 shows that the presence of base only results in a weak secondary flow, evident by the nonzero values of the horizontal (<U>) and vertical (<W>) mean velocity components in the near-base region. However, the ratios <U>/<u> and <W>/<w> are typically between 0.1 and 0.4, indicating that the fluctuating components are always dominant compared to the mean components, and so the secondary flow has a minimal impact on the analysis.
Figure 3 plots both dimensionless r.m.s horizontal (u) and vertical <w> velocity against [xi], which henceforth denotes height above and normal to the base surface. The plots show a typical set of velocity profiles within the blockage layer, that is, at [xi]/[l.sub.0] = 0-2. The data is made dimensionless with [u.sub.0], [w.sub.0], and [l.sub.0], which are the values of <u>, <w>, and l measured at depth z = 2.5M (i.e., [xi] = 0.9M), respectively. The measured values of [u.sub.0], [w.sub.0], and [l.sub.0] for each set of experiment, which were used in the plots presented in this paper, are presented in Table 1. These data also explicitly show, as expected, that both [u.sub.0] and [w.sub.0] decreased as [Re.sub.l] decreased.
Figure 3(a) suggests that within the blockage layer the (u) component experiences an amplification between 0.3 [less than or equal to] [xi] [less than or equal to] 1, attaining a peak value at [xi]/[l.sub.0] [approximately equal to] 0.3. The amplification of the horizontal component is the result of an inter-component energy transfer produced when the turbulent eddies propagate towards and impinge on the tank base . The energy obtained is transferred from normal velocity components where, as shown in Figure 3(b), the (w) components are progressively damped as the distance to the base decreases [xi] [right arrow] 0. The strong self-damping of the normal turbulent intensity within the blockage layer is due to the kinematic boundary condition . The impingement of large-scale turbulent eddies eddy stretches the tangential vorticity components, amplifying the small-scale horizontal turbulence, increasing the horizontal velocity components, denoted here as <[u.sub.p]> . The distance [xi] [approximately equal to] 0.3 where the near base peak <[u.sub.p]> is observed lies within the range obtained by previous studies 0.2 [less than or equal to] [xi] [less than or equal to] 0.4 [7,19]. At this distance, the maximum amplification observed is about 20% higher than [u.sub.0] when [Re.sub.l] > 120, although for [Re.sub.l] < 120, the <u> amplified is <[u.sub.p]>/[u.sub.0] < 1. This result is in line with Hunt and Graham  in the fact that the intensity of amplification depends on the strength of incoming turbulence, where increasing the Reynolds number by a factor of 10 would result in an amplification of 1.3[u.sub.0]. For all the data, the peak <u> observed was within the range of 0.8[u.sub.0] to 1.2[u.sub.0].
The r.m.s vertical velocity w at near impermeable surface was seen to monotonically decrease to zero at the boundary. Hunt and Carlotti  reported that although there might be an amplification of w component, it is not significant due to the wall blocking effect. It is clear from Figure 3 that the profiles for both <u> and <w> components agree reasonably well with the theoretical prediction, although the magnitude for the velocity components was much lower as [Re.sub.l] decreased. Note that the measured <u> and <w> were lower up to a factor of 2 from the predicted values when [Re.sub.l] = 117. Note that for the range of [Re.sub.l] > 117 used in this experiment, consistently similar trends to the data presented in Figures 3(a) and 3(b) were observed.
The same data of [Re.sub.l] = 570 in Figure 3 was also compared with the results from previous experimental work and the theoretical prediction of Hunt and Graham , as shown in Figure 4, to assess the similarity profiles obtained with the previous work. The experimental data from previous work was obtained from Brumley and Jirka  and Fernando et al. , who analysed the interaction of oscillating grid turbulence with a free-surface and a solid plate, respectively. Although the dynamic boundary condition for a free surface (i.e., the work of Brumley and Jirka ) is different (where only w = 0), it appears that the anisotropy of the fields is insensitive to the precise nature of the dynamic boundary condition across the blockage layer . Thus, a comparison of data with the velocity profile obtained by Brumley and Jirka  is allowed. Figure 4 shows that the previous data and measured r.m.s horizontal velocity components produced similar amplification profiles. From <u>/[u.sub.0] = 1 at [xi] = 1, the <u>/[u.sub.0] amplified in the blockage layer reaching a maximum peak at [xi] = 0.2-0.4. On the other hand, the <w>/[w.sub.0] experienced a monotonous decrease towards the surface.
Next, the measurement of the integral length scale, for both horizontal l and vertical length [l.sub.z] scales, was computed. Note that the symbol for the horizontal integral length scale is l, to be in line with other research and, for the vertical integral length scale, the subscript z was added. These measures give both the quantitative turbulent eddy sizes available and how the base modifies the turbulent structure. The horizontal integral length scales were calculated by computing the autocorrelation function based on the r.m.s horizontal (w) velocity components, expressed as
[R.sub.uu] (r,x) = [bar.u(x)u(x + r)]/[bar.[u.sup.2]], (4)
and by obtaining the integral of [R.sub.uu](r) to the first zero crossing. Similar procedure was applied to obtain the vertical integral length scale [l.sub.z] using the autocorrelation function based on the r.m.s vertical (w) velocity components.
Figure 5 shows the results of both horizontal and vertical integral length scales for flows [Re.sub.l] > 117. The normalising factor [l.sub.z0] denotes the value of [l.sub.z] measured at depth z = 2.5M (similar approach when l0 was defined). The horizontal integral length scales l were seen to increase from [xi]/[l.sub.0] [approximately equal to] 1.5 and continued towards the base until [xi]/[l.sub.0] [approximately equal to] 0.5. The energetic eddies approach towards the boundary where the the eddies were stretched along the horizontal plane, increasing the diameter of the eddies. This is a rather similar behaviour observed for splat events, that is, when turbulent eddies or blobs of fluid moving towards the base are flattened and move parallel to the base surface, upon their impingement [3,12]. The increasing size of l is seen to be dependent on the Rei, where flows with [R.sub.el] > 325 can have an increase up to l/[l.sub.0] [approximately equal to] 1.5, whereas the maximum achieved for lower [Re.sub.l] is around 1.3. The vertical integral length scales [l.sub.z], as shown in Figure 5(b), were monotonically decreased in particular from [xi]/[l.sub.0] [approximately equal to] 1.5 towards the base. An increase of the [l.sub.z] was observed at [xi]/[l.sub.0] [approximately equal to] 0.3, indicating an antisplat event, that is, the collision of two splats, resulting in an ejection of fluid (in the normal direction) from the base.
Previous analysis shows that for flows [Re.sub.l] [greater than or equal to] 117, the turbulence generated using the described experimental setup changes according to the standard profile. Here, the experimental data obtained from flows generated at [Re.sub.l] < 117 is discussed.
3.1. The Evolution of r.m.s Fluid Velocity. The r.m.s velocity measurements of <u> and <w> within the blockage layer for flow conditions [Re.sub.l] < 117 are plotted in Figure 6. For a consistent comparison, here the <u> and <w> were obtained at similar distances (as described in Section 2), that is, at [xi]/[l.sub.0] = 0-2 with the predicted lines also shown. The plots demonstrate contrasting profiles to the one obtained in flows with higher Ret (as illustrated in Figure 3), suggesting appreciable changes in the turbulence structure within the blockage layer. The r.m.s horizontal velocity components experienced no amplification at 0 < [xi]/[l.sub.0] [less than or equal to] 0.5 and were consistently constant from 2 [greater than or equal to] [xi]/[l.sub.0] [greater than or equal to] 0.3. This result is coherent with the plot of the normal velocity components (in Figure 6(b)), where they do not monotonically decrease and remain relatively constant from 2 [greater than or equal to] [xi]/[l.sub.0] [greater than or equal to] 1. Thus no (significant) intercomponent energy transfer between the two velocity components occurred. Both r.m.s velocity profiles show a sharp decrease to zero (due to a viscous effect) from [xi]/[l.sub.0] [approximately equal to] 0.3, satisfying the no-slip condition at [xi] = 0. The base does not possess any suppression effect for either horizontal or vertical fluid motions unlike the cases for high [Re.sub.l].
The profiles shown in Figure 6 support the claim of Hunt and Graham , with regard to the importance of high [Re.sub.l] for the occurrence of amplified <[u.sub.p]>. It is believed that, at high [Re.sub.l], the viscous effect is negligible within the blockage layer and the energetic eddies caused a local disruption on the viscous sublayer . At low [Re.sub.l], however, the dominant viscous stress dampens the energy and deviates the velocity profile from the standard curve. The viscous effect is evident as the thickness of the viscous sublayer ([[delta].sub.v]) increases when [Re.sub.l] decreases . Using the relationship [[delta].sub.v] = 2[l.sub.0][Re.sup.-1.sub.1/2], the relative values of [[delta].sub.v]/[l.sub.0] = for [Re.sub.l] < 117 were calculated. The relative thickness of [[delta].sub.v] (for all data) spanned from 0 to 0.18 to 0.23[l.sub.0], shown here in Figure 6 as (two) dashed lines. The thickness of the viscous sublayer attached to the base was (rather) a large fraction of the integral scale [l.sub.0] and contributed to a sharp decrease of the tangential velocity components as the base is approached .
The influence of the impermeability condition can also be traced through the evolution of isotropy I (defined as the ratio of spatially averaged <u> and <w>). The plots of isotropy I = <w>/<u> for high Ret and low [Re.sub.l] at near-base were plotted in Figures 7(a) and 7(b), respectively, to give a direct comparison. The figures were plotted up to [xi]/[l.sub.0] = 3.5 to also include the quasi-isotropic homogeneous turbulence region. Figure 7(a) shows that I is (consistently) constant between the range 1.2 < I < 1.4 from [xi]/[l.sub.0] = 2-3.5. Note that it is difficult (if not impossible) to achieve the ideal value I = 1 due to the normal forcing direction of the grid . This is supported by the I obtained by Hopfinger and Toly , where I = [C.sub.2]/[C.sub.1] = 1.08 (recall (1a) and (1b)). The isotropy plot (shown here in Figure 7) is not only to illustrate the isotropy behaviour within the flow but also to be used to identify the start of the blockage layer, that is, when I starts to decrease . Data shows that, at [xi]/[l.sub.0] = 1.2, the previously constant I (at 1.2 < [xi]/[l.sub.0] < 3.5) experienced a decreasing value and monotonically decreased to zero at the base. Hence, the thickness of the blockage layer (using the experimental setup described) is [approximately equal to] 1.210 above the tank base.
For low [Re.sub.l] (as shown in Figure 7(b)), however, I values were not constant and exhibit different and inconsistent profiles, in particular at near base. Isotropy remains relatively constant until [xi]/[l.sub.0] [approximately equal to] 1.2, before suppression or amplification of isotropy at [xi]/[l.sub.0] [approximately equal to] 0.3. At low [Re.sub.l], the turbulence impacted on the base, which is perceived to be in a decaying state. In such conditions, the intercomponent energy transfer is crucially affected, where the intercomponent energy transfer near a flat (solid) surface depends crucially on the nature (steady vs decaying, homogeneous vs inhomogeneous) of the outer turbulence . For such flow, using I to determine the blockage layer is not effective.
3.2. Integral Length Scales. Figure 8 shows the measurements of both horizontal l and vertical [l.sub.z] integral length scales, made dimensionless with their respective values at z = 2.5M. In the case of a solid boundary with w [right arrow] 0 at the wall and at a distance [xi] < [l.sub.0], the eddies are flattened due to the impingement on the base. It can be said here that although it is not distinctive as the trends shown in Figure 5, the splat-antisplat events occur, even at low [Re.sub.l]. An increase of l towards the base was evident for some flows (for example, [Re.sub.l] = 113 ([??]) shown in Figure 8) and accompanied by a consistent decreasing [l.sub.z] until [xi]/[l.sub.0] [approximately equal to] 0.3 before experiencing a slight increase at near base.
Particular attention is paid for the lowest flow [Re.sub.l] = 74 (symbol [DELTA]), where at near base, [l.sub.z] relatively experienced no changes within the blockage layer, in contrast with the horizontal integral length scales where l was increased up to l/[l.sub.0] [approximately equal to] 1.6. The profiles indicate that the splat events occur due to low [Re.sub.l], where the pressure strain term decreases and subsequent increase of the effect from viscosity reduces the possibility of antisplat events . Changes of the horizontal and vertical integral length scales within the blockage layer at [Re.sub.l] > 74 suggest that the presence of base induced splat-antisplat events. The occurrence of the antisplat event is dependent on the [Re.sub.l], where at low [Re.sub.l] (for example, here at [Re.sub.l] = 74) the incoming turbulence at near base was in a decaying state, such that the energy contained within the eddies is weak and the strong viscous effects near the base dampen the energy.
3.3. Spectra. The fluid velocity measurement using PIV allows the analysis of simultaneous time series velocity components at various distances. Measurement of one-dimensional energy spectra at different distances above the base, at a distances of [xi] = 0.05, 0.5,1, and 2 cm, respectively, was performed. The energy spectra were calculated from the time series of horizontal and vertical velocity components, for each distance. Let f([t.sub.n]) be an input signal (a sequence) of either horizontal or vertical velocity time series and let X(/) be a frequency spectrum; the result of the discrete Fourier transformation of a signal f([t.sub.n]) is expressed as
X(f) = [N-1.summation over (n=0)]([t.sub.n])[e.sup.-(2[pi]i/N)nf], (5)
where n = 0,1, ..., N - 1, i is imaginary unit and [e.sup.-2[pi]i/N] is a primitive Nth root of unity.
By performing a fast fourier transformation to compute the relations in (5), the energy spectral density is calculated as
[phi](f) = X(f)X*(f) / 2[pi], (6)
where X(f) and X*(f) are the real and imaginary components, respectively. The horizontal and vertical energy spectra which are denoted by [[phi].sub.u](f) and [[phi].sub.w](f), respectively, were obtained using (5) and (6), shown in Figures 9 and 10 for flows of [Re.sub.l] [greater than or equal to] 117 and [Re.sub.l] < 117, respectively.
The plot of high [Re.sub.l] [greater than or equal to] 117 will be discussed first, where the one-dimensional spectra for velocity components, taken at near base up to the homogeneous region (i.e., at [xi]/[l.sub.0] = 0.032, 0.317, 0.634, and 1.27), were plotted together. In Figure 9 the frequency spectrum, which represents a homogeneous region, is plotted as dashed-dotted lines for both [[phi].sub.u](f) and [[phi].sub.w](f). As the Reynolds number Ret here is relatively small compared to the Reynolds number in naturally occurring turbulent flow, the [f.sup.-5/3] slope (showing inertial subrange) is not guaranteed. Restricting the Reynolds number restricts the frequency range; that is, the range may not be wide enough to exhibit an inertial subrange. However, the data in Figure 9 suggests that the inertial subrange exists within the frequency range of 0.078 [less than or equal to] f [less than or equal to] 1 for the horizontal component and at a frequency range of 0.4 [less than or equal to] f [less than or equal to] 1.5 for the vertical component. As the base is approached, the inertial subrange region loses its characteristics and very near to the base (at [xi]/[l.sub.0] = 0.032) an inertial subrange barely exists.
The plot shows that as the base was approached, the horizontal components received more kinetic energy and were seen to amplify at [xi]/[l.sub.0] = 0.317. Very near to the base [xi]/[l.sub.0] = 0.032, the energy decreases satisfying the no slip boundary condition. The amplification of [[phi].sub.u](f) is due to the transferred energy from the vertical components, where [[phi].sub.w](f) must decrease to satisfy the kinematic boundary condition, as supported in the analysis shown in Figure 9(b). The presence of tank base has no effect on the high frequency end of the energy spectrum, implying that the net transfer of turbulent kinetic energy from the vertical to the horizontal components of velocity is confined to small frequencies or large eddy scales. Note that similar profiles were obtained for all data and Figure 9 is a representative profile for all data.
The energy spectrum for [Re.sub.l] < 117, shown here in Figure 10, evidently demonstrates low energy for both [[phi].sub.u](f) and [[phi].sub.w](f). Note that the distance [xi] taken was consistent with the distance shown in Figure 9 (i.e., at 0.05, 0.5,1, and 2 cm) above the tank base. The energy for horizontal and vertical components steadily decreased as they approached the solid surface, notably at the low-frequency end. No amplification of kinetic energy in the horizontal direction was observed. At low [Re.sub.l], the flow is believed to have reached the state of decaying turbulence as it approached the tank base. At this state, the viscous effects near the base inhibit the amplification of the horizontal component . The [[phi].sub.u](f) were found to decrease as [xi] [right arrow] 0, where the kinetic energy was progressively dampened by the viscous forces for both large and small scale eddies. Similarly, the plot of [[phi].sub.w](f) (as shown in Figure 10(b)) evidently shows a comparable profile with [[phi].sub.u](f). The kinetic energy for both low and high frequency ends steadily decreased towards the tankbase, suggesting that the range of sizes of active eddies contributing to the energy spectrum decreased. At low flows of [Re.sub.l], the energy transfer from the large-scale eddies to the viscous dissipating (small-scale) eddies is much faster, as the external energy (from the movement of oscillating-grid) at high frequency end is low. Both energy spectrums eventually approached zero in a very thin viscous layer where measurement was made very close to the tank base (i.e., [xi]/[l.sub.0] = 0.032), shown as solid lines in Figure 10, showing a low energy of [approximately equal to][10.sup.-5]. Note that the noise peak at [approximately equal to]23 Hz in Figures 10(a) and 10(b) is speculated to be caused by mechanical vibration, probably due to the grid oscillation.
The experimental work presented here described the characteristics of experimentally oscillating-grid generated turbulence at impermeable surface, for which the base of an acrylic tank was used. The motivation for the study was to find the lower limit of [Re.sub.l], where the turbulence profile at near base follows the rapid distortion theory based on the calculation of Hunt and Graham . Data was obtained from a series of experiments with the turbulent Reynolds number ranging between 74 [less than or equal to] [Re.sub.l] [less than or equal to] 570. Although there are apparent discrepancies between the assumptions of the theory and the experimental conditions, there was clear qualitative agreement between our data trend and that predicted by Hunt and Graham  for flows of [Re.sub.l] [greater than or equal to] 117. The r.m.s horizontal velocity components amplified at [xi]/[l.sub.0] [approximately equal to] 0.3, where the energy is transferred from the decreasing vertical velocity components within the blockage layer. The changes of velocity components mimic the splat-antisplat events, where the eddies flatten, increase in diameter, and eject back the outer fluid. Additionally, the occurrence of splat events was evidenced by a notable increment of horizontal integral length scales and associated monotonous decrease of vertical integral length scales within the blockage layer. Subsequent antisplat events were observed where the vertical length scales [l.sub.z] were seen to increase at [xi]/[l.sub.0] [approximately equal to] 0.3.
For flows [Re.sub.l] < 117, however, the r.m.s horizontal velocity components experienced no amplification. Both (u) and (w) were relatively constant towards the tank base and sharply decreased to zero at [xi] [approximately equal to] 0.3[l.sub.0] due to the viscous effect. The increasing and simultaneously decreasing of horizontal and vertical integral lengthscales, respectively, suggest that, even at low [Re.sub.l] > 74, the splat-antisplat events occurred. However, at the lowest [Re.sub.l] = 74 (available in this study), only splat event was evident. The strong viscous effect at near base dampens the available tangential energy to drive the antisplat events.
The incoming flow [Re.sub.l] < 117 impinging on the tank base showed the characteristics of anisotropic and alters the turbulence profile at near tank base. For flows of [Re.sub.l] < 117, the tank base poses quite an impact on both large and small-scale eddies, where their kinetic energy monotonously decreases when approaching the base. A direct quantitative comparison of the velocity components for flows [Re.sub.l] < 117 with the Hunt and Graham  model shows a departure from the theoretical prediction. For low flows Ret <117, the turbulence at near base was damped, largely contributed by the significant presence of the viscous sublayer.
When using an oscillating grid to investigate the sole effect of turbulence on geophysical phenomena, the lower limit of [Re.sub.l] is within the range 114 [less than or equal to] [Re.sub.l] [less than or equal to] 116. However, we recognised that the precise range may differ with each experimental setup. This research narrowed down the previous possible range [Re.sub.l] limit (i.e., 40 < [Re.sub.l] < 150) that does not follow the theoretical profiles of Hunt and Graham  to a much smaller range of 114 [less than or equal to] [Re.sub.l] [less than or equal to] 116, where this can act as definitive limit. That is, for any experiments conducted below this value, a different physical interpretation of the turbulence at near base is acquired.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
 T. J. McDougall, "Measurements of turbulence in a zero-mean-shear mixed layer," Journal of Fluid Mechanics, vol. 94, no. 3, pp. 409-431, 1979.
 C. Liu, A. Huhe, and L. Tao, "Sediment incipience in turbulence generated in a square tank by a vertically oscillating grid," Journal of Coastal Research, vol. 39, pp. 465-468, 2006.
 J. Bodart, J.-B. Cazalbou, and L. Joly, "Direct numerical simulation of unsheared turbulence diffusing towards a free-slip or no-slip surface," Journal of Turbulence, vol. 11, article N48, pp. 1-18, 2010.
 E. J. Hopfinger and J. A. Toly, "Spatialy decaying turbulence and its relation to mixing across density interfaces," Journal of Fluid Mechanics, vol. 77, no. 3, pp. 155-175, 1976.
 N.-S. Cheng and A. W.-K. Law, "Measurements of turbulence generated by oscillating grid," Journal of Hydraulic Engineering, vol. 127, no. 3, pp. 201-208, 2001.
 J. Hunt and J. Graham, "Free-stream turbulence near plane boundaries," Journal of Fluid Mechanics, vol. 84, no. 2, pp. 209-235, 1978.
 L. Lofdahl, B. Johansson, and C. Andersson, "The influence of temperature on the measurements of Reynolds stresses in shear free turbulence near a wall," Experiments in Fluids, vol. 25, no. 2, pp. 160-164, 1998.
 G. Campagne, J.B. Cazalbou, L. Joly, and P. Chassaing, "Direct numerical simulation of the interaction between unsheared turbulence and a free-slip surface," in European Conference on Computation Fluid Dynamics (ECCOMAS CFD '06), 2006.
 J. J. Orlins and J. S. Gulliver, "Turbulence quantification and sediment resuspension in an oscillating grid chamber," Experiments in Fluids, vol. 34, no. 6, pp. 662-677, 2003.
 M. Belinsky, H. Rubin, Y. Agnon, E. Kit, and J. F. Atkinson, "Characteristics of resuspension, settling and diffusion of particulate matter in a water column," Environmental Fluid Mechanics, vol. 5, no. 5, pp. 415-441, 2005.
 M. Holzner, A. Liberzon, M. Guala, A. Tsinober, and W. Kinzelbach, "Generalized detection of a turbulent front generated by an oscillating grid," Experiments in Fluids, vol. 41, no. 5, pp. 711-719, 2006.
 B. Perot and P. Moin, "Shear-free turbulent boundary layers. I. Physical insights into near-wall turbulence," Journal of Fluid Mechanics, vol. 295, pp. 199-227, 1995.
 H. J. S. Fernando, E. J. List, and I. A. Hannoun, "Turbulence structure near a sharp density interface," Journal of Fluid Mechanics, vol. 189, pp. 189-209, 1988.
 H. J. S. Fernando and I. P. D. de Silva, "Note on secondary flows in oscillating-grid, mixing-box experiments," Physics of Fluids A, vol. 5, no. 7, pp. 1849-1851, 1992.
 S. I. Voropayev and H. J. S. Fernando, "Propagation of grid turbulence in homogeneous fluids," Physics of Fluids, vol. 8, no. 9, pp. 2435-2440, 1996.
 Digiflow, Dalziel Research Partners, Cambridge, UK, 2014.
 S. P. McKenna and W. R. McGillis, "Observations of flow repeatability and secondary circulation in an oscillating grid-stirred tank," Physics of Fluids, vol. 16, no. 9, pp. 3499-3502, 2004.
 J. C. R. Hunt and P. Carlotti, "Statistical structure at the wall of the high reynolds number turbulent boundary layer," Flow, Turbulence and Combustion, vol. 66, no. 4, pp. 453-475, 2001.
 B. H. Brumley and G. H. Jirka, "Near-surface turbulence in a grid-stirred tank," Journal of Fluid Mechanics, vol. 183, pp. 235-263, 1987.
 W. H. M. Wan Mohtar and R. J. Munro, "Threshold criteria for incipient sediment motion on an inclined bedform in the presence of oscillating-grid turbulence," Physics of Fluids, vol. 25, no. 1, Article ID 015103, 2013.
 J. Magnaudet, "High-Reynolds-number turbulence in a shear-free boundary layer: revisiting the Hunt-Graham theory," Journal of Fluid Mechanics, vol. 484, pp. 167-196, 2003.
 W. H. M. WanMohtar, The interaction between oscillating-grid turbulence and a sediment layer [Ph.D. thesis], School of Civil Engineering, University of Nottingham, Nottingham, UK, 2011.
W. H. M. Wan Mohtar and A. ElShafie Department of Civil & Structural Engineering, Faculty of Engineering and Built Environment, Universiti Kebangsaan Malaysia, 43600 Bangi, Malaysia
Correspondence should be addressed to W. H. M. Wan Mohtar; firstname.lastname@example.org
Received 11 April 2014; Revised 11 July 2014; Accepted 28 July 2014; Published 27 August 2014
Academic Editor: Jianhong Xu
TABLE 1: A summary of the [u.sub.0], [w.sub.0], and [l.sub.0] for each ow studied. [Re.sub.l] [u.sub.0] [w.sub.0] [l.sub.0] (cm/s) (cm/s) (cm) Flows with high [Re.sub.l] 570 1.606 2.154 1.78 474 1.502 2.094 1.58 325 1.336 1.742 1.22 209 0.604 0.871 1.73 117 0.570 0.814 1.35 Flows with low [Re.sub.l] 113 0.482 0.556 1.18 95 0.341 0.449 1.40 76 0.297 0.378 1.28 74 0.302 0.361 1.22
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Mohtar, W.H.M. Wan; ElShafie, A.|
|Publication:||The Scientific World Journal|
|Date:||Jan 1, 2014|
|Previous Article:||Diesel fuel from used frying oil.|
|Next Article:||Estimating tree height-diameter models with the Bayesian method.|