Realisation of Fractal Grid-Induced Turbulence Strength with PTFV: Effects of Grid Geometry

The unravelling of multilength-scale insert-generated turbulence, particularly, the induced-forcing plays critical role in the fundamental comprehension of energy formation and decay as a function of grid conformation. This study experimentally investigates the flow mechanical characteristics at ReDh = 4.1 × 104 for a regular-grid (RG), single-square-grid (SSG) and six 2D planar space-filling square-fractal-grids (SFG) of different fractal iterations (N), thickness ratios (tr) and blockage ratios (σ) via piezoelectric thin-film flapping velocimetry (PTFV). Thin-film’s tip-deflection (δrms) and voltage response (Vrms) analysis along the grids’ centreline reveals increasing flow fluctuation strength with increasing σ, tr and decreasing N, owing to higher shedding intensity of lower frequency, larger scale energy-containing vortices from thicker first iteration bar. However, higher: energy dissipation rate, centreline mean velocity decrement rate and local deceleration experienced in the turbulence decay region of larger tr grid, along with additional fractal scales lead to less potent flow-structure-interplay on thin-film undulation. More importantly, SSG-generated turbulence enables the generation of average (Vrms, δrms) and millinewton turbulence forcing Frms that are respectively, 9× and 5× larger than RG of similar σ, and 2× larger than the best performing N = 3 SFG. Our findings disclose the importance of grid geometrical management for effective utilisation of turbulence-generating grids in engineering applications.


Introduction
Corrugated surfaces (Tinney et al. 2020;Guilarte Herrero et al. 2021), flaps (Li et al. 2002;Rosti et al. 2018) and turbulence-generating grids (Batchelor and Townsend 1947;Comte-Bellot and Corrsin 1966;Geipel et al. 2010;Fragner et al. 2015) are examples of passive flow control techniques widely adopted to generate fluid flow perturbations to alter the flow field to a desired state.The turbulence generated by planar fractal grids, which comprise self-repeating patterns superimposed at multiple scales, have been extensively studied both experimentally (Hurst and Vassilicos 2007;Valente and Vassilicos 2011;Nagata et al. 2013;Goh et al. 2021) and numerically (Laizet and Vassilicos 2011;Suzuki et al. 2013;Zhou et al. 2014;Mok et al. 2023) over the past decade and a half.Pioneering work by Hurst and Vassilicos (2007) has discovered the existence of a protracted turbulence generation region in the lee of square fractal grid (SFG), where the centreline turbulence intensity increases progressively to a peak, followed by a decay region.It is acknowledged that regular grid (RG) generates higher centreline turbulence intensity peak over a short region very close to the grid than SFG of the same blockage ratio σ and effective mesh size M eff 1 (Laizet and Vassilicos 2011;Melina et al. 2016).Nevertheless, the higher vorticity and turbulence intensity flow produced by SFG in the protracted region (Laizet and Vassilicos 2011) makes SFG highly attractive for thermal mixing (Teh et al. 2015;Skanthan et al. 2018), convective heat transfer (Hoi et al. 2019a, b) and energy harvesting (Ferko et al. 2019a, b) applications.
In general, grid turbulence can be divided into the near-field and far-field regions (Simmons and Salter 1934;Hearst and Lavoie 2016).The non-equilibrium near-field is identified as the region that experiences strong cross-stream inhomogeneities in which its turbulence evolution is directly influenced by the initial flow generating conditions (Hearst and Lavoie 2016).As one goes sufficiently far downstream, the homogeneous far-field turbulence generated by SFG behaves as one might expect for an RG undergoing powerlaw decay (Hearst andLavoie 2014, 2016).While rapid decay also exists close to the RG (Krogstad and Davidson 2012;Isaza et al. 2014), when the mesh length of RG is made identical to the largest grid bar length L 0 of SFG with both grids having equivalent σ, the extended region of rapid decay in the lee of SFG relative to RG results in a less energetic far-field (Hearst and Lavoie 2016).Therefore, apart from producing stronger turbulent diffusion and higher scalar mixing (Nedić and Tavoularis 2016), RG also generates higher turbulence intensity and Taylor microscale Reynolds number Re λ than the SFG in the farfield region (Hearst and Lavoie 2016).This contradicts the assertion that the scalar diffusion behind SFG is more efficient than RG, since the mesh sizes of the RGs used by Suzuki et al. (2010) and Laizet and Vassilicos (2015) were much smaller than the L 0 of SFG.
Hitherto, various research has been undertaken to examine the interesting interplay between the SFG geometric parameters and the unique multilength-scale flow generated.While the magnitude of the turbulence intensity peak is reported to be proportional to the ratio of the width to the length of largest grid bar t 0 /L 0 (Gomes-Fernandes et al. 2012), its peak distance x peak can be approximated from the wake-interaction length scale (Mazellier and Vassilicos 2010;Gomes-Fernandes et al. 2012), which represents the location on the centreline where the wakes from the largest grid bars encounter each other.It was discovered that an increase in thickness ratio t r (the ratio of largest to smallest bar width) shortens the x peak but at the same time increases the turbulence intensity and Reynolds number of the generated flow (Hurst and Vassilicos 2007;Mazellier and Vassilicos 2010).Weitemeyer et al. (2013) experimentally showed that the turbulence intensity heightens with increasing σ, whilst x peak is pushed further downstream with the addition of a fourth iteration to Hurst and Vassilicos 2007), where T 2 is the tunnel's square cross section and P the grid perimeter length.the SFG design.This is further supported by the experiments conducted by Omilion et al. (2018) where the high turbulence intensity leeward of SFG displays an elongated trend when additional fractal scales are added to the grid.Moreover, Nagata et al. (2017) demonstrated that the smallest grid elements of higher iteration SFG plays a significant role for the nearly constant Taylor microscale and integral length scale detected in the turbulence decay region.On the other hand, the presence of smaller geometrical iterations in SFG with fractal iteration N = 4 could lead to a less intense vortex shedding effects along the centreline as compared to a single square grid (SSG) (Melina et al. 2016).
All the above-mentioned experimental studies employ either costly hot wire anemometry (HWA) or particle image velocimetry (PIV) to quantitatively detect turbulence.The qualitative characterisation of fractal grid-induced turbulence based upon the direct interaction of the flow with a piezoelectric thin-film has only been recently initiated (Ferko et al. 2018;Ferko et al. 2019a, b;Lee et al. 2021).In the previous experimental work by Lee et al. (2021), mechanical characteristic of the centreline streamwise and cross-sectional turbulence generated from a space-filling SFG and RG with comparable σ and M eff were experimentally unveiled via piezoelectric thin-film flapping velocimetry (PTFV).Analysis of the thin-film's physical response revealed the superior advantages of SFG-generated turbulence over RG in terms of flow fluctuation strength and the corresponding turbulence forcing (turbulence-induced excitation force acting on the thin-film).
In view of the significant impact of grid geometry on the flow generated by an SFG and RG, a parametric study on the effects of SFG design towards the resulting turbulence mechanical characteristics is of utmost importance for effective utilisation of fractal grids in a wide variety of engineering applications.To date, the interplay between the grid geometric parameters and the mechanical characteristics of SFG-generated flow, particularly the turbulence forcing has yet to be explored.There also appears to be limited work on the individual influence of different σ and N towards the resulting flow characteristics, as the σ is not held constant in most studies (Nagata et al. 2017;Omilion et al. 2018) involving SFG of different fractal scales.Furthermore, the SFGs tested in the majority of the research (Hurst and Vassilicos 2007;Laizet and Vassilicos 2011;Valente and Vassilicos 2011;Nagata et al. 2013;Hearst and Lavoie 2016;Goh et al. 2021) are of low to medium σ with σ < 0.4.Essentially, the turbulent wakes and jets generated by high σ grid would become unstable, and the resulting turbulence is less likely to be homogeneous besides taking longer to reach uniformity (Corrsin 1963).However, through careful fine-tuning of the fractal-fin geometrical parameters, the intense turbulence level and vigorous flow fluctuation characteristics of the flow generated by SFG, partially-covered SFG and partiallycovered RG with σ = 0.49 have been recently proven to be highly critical and preferable in enhancing the thermal dissipation of plate-fin heat sink (Hoi et al. 2019a(Hoi et al. , 2022;;Chew et al. 2022).In addition, the potent flow fluctuation strength and turbulence forcing induced by SFG of σ = 0.49 allows high piezoelectric voltage output to be secured, thus presenting promising prospect for fluidic energy harvesting from high-blockage SFG-induced turbulence (Lee et al. 2021).
Hence, the present work seeks to investigate the fluid-structure interaction between a flexible piezoelectric thin-film and the flow generated by different geometric parameters high-blockage passive grids via in-house PTFV (Lee et al. 2021).Since the piezoelectric thin-film immersed in insert-generated flow domain flaps to generate voltage as a function of fluid flow undulation, the tip deflection and voltage response of a polyvinylidene fluoride (PVDF) piezoelectric cantilever beam are measured at various centreline streamwise distances leeward of an RG, an SSG and six SFGs of different N, t r and σ to disclose the respective flow fluctuation strength.The local equivalent turbulence statistics determined from the film undulation are also analysed to provide insights into the insert-induced centreline streamwise turbulent flow characteristics and induced forcing.The in-house developed PTFV would serve as a potential candidate for a low cost, in situ qualitative evaluation of the fundamental mechanical characteristics of the turbulent flow structures generated by SFGs of different geometric parameters, such that the SFG-generated flow can be better manipulated by restructuring the flow through proper turbulence management to best match the intended application.At the same time, the voltage data acquired from experiments could be particularly useful for future energy harvesting studies.For instance, piezoelectric thin-film could be placed in the lee of the turbulator that induces greatest piezoelectric voltage generation to maximise the electrical energy harvested for self-powering low-power sensors to monitor, control, and manage heating, ventilation and air conditioning (HVAC) systems.

Turbulence-Generating Grids
A total of eight turbulence-generating grids [ Grid t r N ( ) ] made from 10 mm thick transparent acrylic sheet are examined in this study: a regular grid [ RG(0.492) ], a single square grid [ SSG(0.493) ] and six different 2D planar space-filling square fractal grids [ SFG 6.5  2 (0.493) , SFG 6.5  3 (0.493) , SFG 6.5 3 (0.382) , SFG 6.5 3 (0.270) , SFG 4.0 3 (0.493) , SFG 9.8 3 (0.493) ].The sche- matic of the grids along with their geometrical details are presented in Fig. 1 and Table 1, respectively.Note that A o denotes the grid's central opening area and x * the wake-interac- tion length scale.x * was first introduced by Mazellier and Vassilicos (2010) as, where L 0 and t 0 are the length and lateral thickness (width) of the largest grid bar, respectively.Equation 1 has been widely used in literature (Laizet and Vassilicos 2011;Laizet et al. 2015;Hearst and Lavoie 2016;Melina et al. 2016;Nedić and Tavoularis 2016;Nagata et al. 2017) to compute x * with the intention of providing an estimation of the centreline turbulence intensity peak location and also to normalise the results obtained for fractal grids of different geometric characteristics.By taking into consideration the effect of inlet free-stream turbulence on the wake development in a water tunnel, Gomes-Fernandes et al. (2012) proposed a new wake-interaction length scale x peak * , where C d is the drag coefficient, while α refers to the parameter that takes into account the free-stream turbulence intensity and length scale.Considering that the free-stream turbulence is negligible for wind tunnel experiments (Gomes-Fernandes et al. 2012), x * is calcu- lated using Eq. 1 for the present study.
The RG considered here (Fig. 1h) has a blockage ratio (ratio of lateral plane area covered by grid to tunnel cross-sectional area) σ = 0.492, and its geometry is the same as that used by Lee et al. (2021).Apart from SSG (Fig. 1a) and SFG 2 (Fig. 1b) where both inserts are supported at the corners by four struts with 3 mm lateral thickness, the (1) rest of the inserts have cross-sectional area equivalent to the test section.All SFGs and SSG are designed to have almost similar first, second and third iterative bar lengths (L 0 ≈ 91 mm, L 1 ≈ 46 mm, L 2 ≈ 20 mm) except for SFG 9.8 3 (0.493) , of which, its geometry is made identical to the fractal insert numerically optimised by Hoi et al. (2019a) with L 0 = 77 mm.On the other hand, the first to the third iterative lateral bar thicknesses (t 0 , t 1 and t 2 ) are altered to individually vary the N, t r and σ of the SFGs while keeping the other two parameters constant.The three cases investigated in the present study include: (i) different N with fixed t r and σ (N = 1, 2, 3; t r = 6.5; σ = 0.493), (ii) different t r with fixed N and σ (t r = 4.0, 6.5, 9.8; N = 3; σ = 0.493) and (iii) different σ with fixed N and t r (σ = 0.270, 0.382, 0.493; N = 3; t r = 6.5).

Experimental Setup and Procedure
Experiments are performed in the same transparent acrylic wind tunnel with the same setup as that in Lee et al. (2021) study.Figure 2 shows the schematic of the experimental setup.The tunnel is 4.41 m in length and has a cross-section of 0.160 m × 0.160 m.Flow straighteners are placed at the inlet right after the bellmouth and in front of the axial drawthrough fan (Kruger MTD200, SG) at the outlet.The centreline air inlet velocity is fixed at v in = 4 m/s with hydraulic diameter-based flow Reynolds number Re Dh = 4.1 × 10 4 .Each turbulence-generating grid is inserted at 2.02 m from the inlet.
A unimorph piezoelectric thin-film (TE Connectivity LDT1-028K, CH), which consists of a PVDF film laminated to a Mylar sheet is cantilevered and placed parallel to the flow.Table 2 tabulates the properties of the piezoelectric cantilever beam (Measurement Specialties 1999, 2015;Danesh-Yazdi et al. 2016;Lee et al. 2021).With the piezoelectric output connected to a 10 MΩ external electrical load, the voltage signal is sampled and digitised for a period of 60 s at a sampling rate of 10 3 Hz using a data acquisition device (LabJack U3-HV, US) with a voltage resolution of 4.88 mV/bit.Under the illumination of an LED light source, the lateral fluctuations of the thin-film, which is fluorescently highlighted at the edges are simultaneously captured for 10 s using an industrial high-speed camera (fps4000-720, UK).The images are recorded at a rate of 10 3 fps with a pixel resolution of 1280 × 720.A simple LED circuit (Fig. 2b) (Lee et al. 2021) is also employed to ensure synchronised measurement of the voltage signal and piezoelectric undulation.For each of the grid investigated, the thin-film is positioned along the grid's centreline at different grid-film distances: 20 ≤ x ≤ 400 mm corresponding to 0.13 ≤ x/D h ≤ 2.50, where D h = 160 mm is the hydraulic diameter of the wind tunnel.The experiment is repeated six times for each x/D h to ensure experimental repeatability.Finally, the local centreline mean velocities in x, y and z-directions at x peak and x/D h = 0.13 and 2.50 are measured using a hotwire anemometer (Testo 405i, DE) of accuracy ± 0.3 m/s.The hand-held probe is first inserted vertically into the wind tunnel along the y-axis with its sensor head facing the streamwise air flow to measure the instantaneous longitudinal velocity (U x,ins ).After the probe is rotated 90° for instantaneous spanwise velocity (U z,ins ) measurement, it is inserted horizontally along the z-axis to measure the instantaneous vertical velocity (U y,ins ).For each velocity component, three repeated measurements are performed, and for each measurement the data is logged for a period of 60 s to obtain its time average velocity.

Table 1 Geometrical parameters of the grids
Parameter SSG(0.493)SFG 6.5 2 (0.493) SFG 6.5  3 (0.493) SFG 6.5 3 (0.382) SFG 6.5 3 (0.270) SFG 4.0 3 (0.493) SFG   7.69 [c]   Natural frequency (Hz) The piezoelectric thin-film tip deflection δ is determined from the experimentally recorded images using the previously developed Thin-Film Fluctuating Image-Detection Algorithm (Lee et al. 2021).The δ detected has an average uncertainty of 0.02 mm.The signal-to-noise ratio of δ is 4.5 for small amplitude fluctuation, which is 7.4× lower than the voltage signal-to-electrical noise ratio.The local equivalent lateral flow velocity fluctuation v' is then calculated from the time derivative of δ, i.e., d dt prior to evaluating the corresponding turbulence statistics and induced forcing in the following section.It is worthwhile to mention that when v in = 5.3 m/s (viz.the maximum v in in the absence of a turbulence-generating grid), the piezoelectric thin-film positioned at x/D h = 0.56 on the centreline of the test section without grids has RMS voltage V rms = 0.18 V and RMS tip deflection δ rms = 0.09 mm.Such inconsiderable thin-film response signifies that the wind tunnel has a rather low background turbulence.

Effects of Fractal Iteration N on Centreline Streamwise Flow Characteristics
The experimentally recorded RMS voltage V rms and RMS tip deflection δ rms of the piezoelectric thin-film, as well as their corresponding dominant frequencies f V and f δ at different streamwise distances x/D h along the centreline of four high σ (0.493) inserts: SSG, SFG 2 (N = 2) and SFG 3 (N = 3) with t r = 6.5, and RG are presented in Fig. 3a and b.The centreline streamwise evolution of the equivalent lateral integral length scale L v normalised against the tunnel width T is plotted in Fig. 3d to determine the size of the large-scale energy-containing eddies generated by the inserts.L v is evaluated from Eq. 3 via integrating the autocorrelation function of v' with respect to time-lapse ∆ up to the first zero-crossing time τ zc of the autocorrelation plot.
As seen in Fig. 3a, the V rms and δ rms induced by SSG and SFG-generated turbulence exhibit a similar profile as described in literature (Hurst and Vassilicos 2007;Mazellier and Vassilicos 2010;Melina et al. 2016).The V rms and δ rms in the turbulence generation region of SSG, SFG 2 and SFG 3 increase, respectively, to a peak at x peak /D h = 1.0, 1.2 and 1.0 before decreasing exponentially further downstream in the decay region.As for the RG, a power-law decay is observed throughout the entire streamwise distances.Once the SFG-generated flow is fully recovered from the inhomogeneity of the near-field at region adequately far from the grid, it is anticipated that the SFG-induced V rms and δ rms will also undergo a power-law decay if the experiments were to be conducted further downstream in the decay region (x/D h > 2.5) as supported by the studies conducted by Hearst andLavoie (2014, 2016).In view of the fact that the magnitude of the turbulence intensity peak is proportional to the t 0 /L 0 ratio of SFG (Gomes-Fernandes et al. 2012), the V rms and δ rms peaks are the most prominent for SSG, with values of 7.6 V and 3.1 mm, which are about 2 × , 3 × and 9 × larger than SFG 2 , SFG 3 and RG, respectively.SSG has an earlier peak point than SFG 2 as a result of its smaller wake-interaction length scale x * , which represents the location on the centreline where the wakes from the largest grid bars interact, bounding the x peak (Mazellier and Vassilicos 2010).However, due to the large t 0 /L 0 ratio of SSG (0.38) and SFG 2 (0.31), their x peak are larger than expected despite having a smaller x * than SFG 3 . (3) The reason is that a thick first iteration grid bar not only leads to smaller A o , which intensifies the contraction effect to further accelerate the flow, but also enlarges the size of the recirculation regions immediately in the lee of the largest bar, and consequently diverges the streamlines of the flow towards the tunnel walls (Laizet and Vassilicos 2011).Hence, the wakes from the largest bar tend to meet further downstream and cause a delay in wakeinteraction.In addition, an increase in x * creates an extended turbulence generation region, giving rise to an increase in the full width at half maximum (FWHM) of the peaks with increasing N. Nevertheless, the FWHM of the peaks for SFG 2 is slightly higher than anticipated due to the aforementioned wake-interaction-delay.
As shown in Fig. 3b, there is a general increase in the f V and f δ trends with increasing x/D h for all inserts, of which, the L v /T profiles in Fig. 3d are virtually a reverse of the thin-film response dominant frequencies f profiles.Without the addition of smaller scale bars, the SSG-generated flow, particularly in the domain proximate to grid (x/D h ≤ 0.6) is mainly dominated by the considerably large-scale (14 mm < L v < 16 mm), low-frequency (35 Hz < f < 40 Hz) eddies generated from the thick square bar, which evolve into smaller ones as proceeding downstream.With constant t r (for the case of SFG) and σ, t 0 increases as a function of decreasing N. The overall f is therefore the lowest for SSG followed by SFG 2 , SFG 3 and RG owing to the slowly rotating, larger L v /T vortices shed from the thicker first iteration grid bar, which act as the primary flow components in inducing thin-film flapping.When additional fractal scales are added to the grid, energy are re-distributed to a broader range of scales in the flow (Nedić et al. 2015).The presence of higher frequency, small-sized eddies could lead to less potent flow structures driving the thin-film fluctuations, thus an overall lower V rms and δ rms detected for higher N grid (Fig. 3a).
Figure 3c compares the centreline RMS turbulence-induced excitation force F rms acting on the thin-film (turbulence forcing) in the lee of RG and SFGs of different N. The time-dependent normal flow excitation force F y is calculated from the following equation (Sodano et al. 2004;Elvin et al. 2006), where is the electromechanical coupling coefficient, m, c and k are the mass, damping and stiffness of the piezoelectric beam.These parameters are determined using the same approach as our previous study (Lee et al. 2021) and are listed in Table 2.The F rms calculated for all inserts display a similar streamwise evolution as the V rms and δ rms in Fig. 3a with only a minute shift in the x peak of SFG 3 .As high as 28 mN is secured from the SSGgenerated flow at the x peak , followed by the flow generated from SFG 2 (16 mN), SFG 3 (9 mN) and RG (7 mN).Interestingly, in spite of the lower V rms and δ rms induced by RG-generated turbulence in comparison to SSG at x/D h = 0.13, SFG 2 and SFG 3 at x/D h = 0.25, the F rms for RG is larger than the three grids at the respective locations.This is the consequence of the additional forcing provided by the small-scale turbulence as supported by the energy spectra E v in Fig. 6, of which, the minor peak E v (as indicated by the red arrow) of RG is an order of magnitude higher than SSG at x/D h = 0.13 and (SFG 2 , SFG 3 ) at x/D h = 0.25.This result is in line with Lee et al. (2021) finding where the F rms calculated takes into account the broad-band random forcing across various eddy sizes including those smallscale eddies of higher frequencies, whereas V rms and δ rms merely reflect the vortex shedding intensity of the dominant large eddies.

Effects of Fractal Thickness Ratio t r on Centreline Streamwise Flow Characteristics
Figure 4a-d show the V rms and δ rms , f, F rms and L v /T against x/D h along the centreline of three SFGs (N = 3, σ = 0.493) of different t r .All inserts demonstrate an increase in V rms , δ rms and F rms with x/D h followed by an exponential decay beyond the peak (Fig. 4a, c).The V rms , δ rms and F rms of SFG 9.8 at x peak /D h = 0.8 are 1.5× to 1.6× higher than SFG 6.5 at x peak /D h = 1.0, and 1.6× to 1.7× higher in comparison with SFG 4.0 at x peak /D h = 1.2.This is consonant with the results reported in literature (Hurst and Vassilicos 2007;Mazellier and Vassilicos 2010) where turbulence intensity heightens and x peak shortens when increasing t r .The decrease in x * with t r brings about a decrease in both x peak /D h and FWHM of the peaks.By keeping the σ and N of the grid constant, first iteration bar thickens while second and third iteration bars thinned when increasing t r , giving rise to a wider range of eddy's scale in the flow.As depicted in Fig. 4b, d, the vortical structures generated from SFG 9.8 with large t 0 /L 0 ratio (0.32) are relatively low in frequency (f ≈ 40 Hz) and large in size (L v ≈ 15 mm) in the region very close to the grid.For all three inserts, we once again notice a general increase in f and decrease in L v /T with increasing x/D h .The V rms , δ rms and F rms in Fig. 4a, c are greater for higher t r grid in the turbulence generation region as a consequence of the larger size energy-containing eddies generated from the thicker first iteration bar.However, the V rms , δ rms and F rms for all three inserts coincide to almost similar values in the decay region, whereby at x/D h = 2.5, the values are the largest for SFG 4.0 and the lowest for SFG 9.8 .This could be explained by the higher energy dissipation rate experienced by larger t r grid as evidenced in Fig. 6, where both (4) m δ + c δ + k − V = F y major and minor peaks' intensities in the decay region of SFG 9.8 attenuate 3× and 4× faster than SFG 6.5 and SFG 4.0 , respectively.Considering that smallest eddies decay the fastest, the finer small-scale structures present in higher Reynolds number flow (Davidson 2004) generated along the centreline of smaller A o grid, as well as the smaller size vortices shed from the thinner second and third iteration bars, may have prompted the energy dissipation rate to increase with t r in the decay region.

Effects of Blockage Ratio σ on Centreline Streamwise Flow Characteristics
Figure 5a-d illustrate the downstream profile of the centreline V rms and δ rms , f, F rms and L v /T leeward of three different σ SFGs with N = 3 and t r = 6.5.The V rms , δ rms and F rms for all the three inserts in Fig. 5a, c demonstrate turbulence generation and decay albeit with a slightly irregular trend detected for SFG(0.270)'sF rms distribution.SFG(0.493) has peak V rms , δ rms and F rms that are approximately 2× larger than SFG(0.382)and 3 to 5× larger than SFG(0.270).The peak location is further and the FWHM is wider when reducing σ.This again is the result of larger x * owing to the smaller t 0 of a lower-blockage SFG and there- fore, a flatter peak is attained for SFG(0.270).
In view of the increase in all three iterations lateral bar thickness with σ, the multiscale jets produced from the smaller openings are accelerated.Moreover, there is an overall decrease in f and increase in size (L v /T) of the multilength-scale turbulent eddies generated from higher σ grid at each x/D h (Fig. 5b, d), giving rise to more effective (energy-containing) flow structures exciting the thin-film flapping.Unlike the variation in f and L v /T with respect to x/D h observed for the previous two cases, the f and L v /T for SFG(0.382) and SFG(0.270) are rather constant throughout the entire streamwise distances investigated.This suggests that turbulence generated from low-and medium-blockage SFGs have a higher degree of homogeneity compared to high-blockage SFGs.

Energy Spectra
Figure 6a presents the centreline streamwise 3D energy spectra E v of all eight inserts determined from the Fourier transformations of the v' autocovariance function while Fig. 6b shows the top view of the respective spectra.For each x/D h in Fig. 6a, the major peak appearing at frequency f close to the first resonant mode of piezoelectric beam, which was also observed by Laizet et al. (2015) in the longitudinal energy spectra for SFG, Isaza et al. (2014) for RG and Melina et al. (2016) for SSG, is the result of vortex shedding from the grid bars (Gomes-Fernandes et al. 2015;Melina et al. 2016;Nagata et al. 2017).These equivalent vortex shedding frequencies are in excellent agreement with the f V and f δ in Figs.3b, 4b and 5b.The high vortex shedding intensity observed in the turbulence generation region intensifies with t r and σ but weakens with N, which is identical to the results in Fig. 7a.The shedding effect is the most pronounced for SSG since it has the largest t 0 .The major peak E v (E v,max ) at x peak /D h leeward of SSG is approximately 2× and 40× higher than SFG 9.8 3 (0.493) and RG, respectively.One can also infer that the more intense the vortex shedding, the higher the equivalent vortex shedding energy as represented by the broader f range of high energy fluctuations with E v > 10 -4 (Fig. 6b).Consequently, more energetic flow fluctuations inducing the thin-film vibration, yielding larger δ rms and V rms .
On the contrary, minor peaks with relatively smaller E v tend to occur at higher f (200 Hz < f < 350 Hz) associated to the second vibration mode of the piezoelectric beam.These minor peaks are the most prominent for RG especially in the region near to grid, 3 (0.270) , SFG 6.5 3 (0.382) and SFG 6.5 3 (0.493) 1 3 signifying that higher frequency small-scale eddies possess more significant contribution to RG-induced turbulence when compared to SFG and SSG.Moreover, it is noteworthy that the slight shift in SFG 6.5 3 (0.382) F rms peak location from x/D h = 1.2 to x/D h = 1.4 (Fig. 5c) could be due to the additional forcing provided by the small-scale eddies, given that the minor peak E v at x/D h = 1.4 is 1.3× higher than at x/D h = 1.2 albeit the former has a major peak E v that is 1.1× lower than the latter.Last but not least, the streamwise evolution of the E v,max for all inserts are consistent with the streamwise profiles previously secured for V rms and δ rms in Figs.3a, 4a and 5a.Both the major and minor peaks in the decay region diminish as one proceeds downstream owing to energy dissipation.

Interplays Between N, t r and σ on the Overall Flow Mechanical Characteristics
Figure 7a-e plot the average V rms (V avg ), f V (f Vavg ), L v /T (L v /T avg ), F rms (F avg ) and damping ratio ζ (ζ avg ) downstream of the grids' centreline against N, t r and σ.It is immediately apparent from Fig. 7a-c that the V avg and L v /T avg increase with t r and σ but decrease with N, in contrast to f Vavg .This is attributable to the higher shedding intensity of lower frequency, larger scale energy-containing vortices from the thicker first iteration bar as evidenced in Figs. 6, 7b and c.The F avg calculated in Fig. 7d ranges from 2.0 to 18.4 mN and exhibits a similar trend as V avg with respect to N, t r and σ, denoting that increasing turbulenceinduced excitation force acting on the thin-film is generally accompanied by an increase in δ and V.One can also deduce that larger V signifies larger L v but lower frequency eddies present in the flow.In the context of energy harvesting for vortex-induced vibration-based devices where potent flow fluctuation strength and induced forcing are highly desirable, SFG 9.8 3 (0.493) with considerable V avg and F avg turns out to be the best performing grid among all the N = 3 SFGs, followed by SFG 6.5 3 (0.493) .Thus, despite the fact that the tur- bulent flow generated by high-blockage grid would become unstable and is less likely to be homogeneous (Corrsin 1963), the turbulence induced by SFG with higher σ is proven in our current work to be favourable for piezoelectric energy harvesting application.The V avg (4.6 V), L vavg (13.9 mm) and F avg (18.4 mN) are the highest while f Vavg (42.5 Hz) is the lowest in SSG-induced turbulence, of which, the V avg and F avg for SSG are twice of SFG 9.8 3 (0.493) , and (9×, 5×) higher than the RG, respectively.All the SFGs investigated demonstrate superior flow fluctuation strength and induced forcing than the RG, which is consistent with the earlier studies (Laizet and Vassilicos 2011;Ferko et al. 2019b;Lee et al. 2021), except for SFG 6.5  3 (0.270)where its V avg and F avg are 1.2 × and 1.7 × lower than RG owing to its 1.8× lower σ than the latter.
In spite of the f Vavg variation with respect to grid geometry, the piezoelectric beam positioned in the lee of all eight inserts vibrates with variations of less than 30% from its own natural frequency f n , denoting that f V and f δ are associated to the f n of piezoelectric beam amalgamate with the characteristic frequencies of vortices shed from the multi-scale grid bars.It is observed in Fig. 7c that the L v /T avg gradient between N = 2 and 3 is less steep than that of N = 1 and 2. Notwithstanding the fact that the increase in L v /T avg when N decreases from 3 to 2 is the result of larger t 0 as well as the removal of the third iteration bar, the thinner second iteration bar of SFG 6.5 2 (0.493) compared to SFG 6.5 3 (0.493) would offset the L v /T avg increment to a certain extent.This implies that the L v calculated in the present study (10.9 mm < L v < 15.8 mm) is essentially affected by t 0 , yet the size of the eddies generated from higher iteration bars also impose some minor effect on the L v .Given that the whiskers in the line plots represent the standard deviation of the values over the entire streamwise distances investigated, the extremely large whiskers observed for SSG and SFG 9.8 3 (0.493) could be an indication of the highly inhomogeneous flow generated.
Based on Lee et al. (2021) finding, the possible causes for the high ζ attained along the centreline leeward of SFG are (i) the high velocity jet flowing through the central opening, which may possibly hinder the up and down movement of the thin-film, and (ii) the multilength-scale flow structures generated, which makes the fluid domain around the thin-film to become rather crowded.Hence, the flow generated from RG of uniform geometry has the lowest ζ avg among the grids (Fig. 7e).In addition, one can also infer that the smaller the central opening (larger t 0 ), the stronger the central jet as a result of contraction effect, thereby higher ζ avg .Although at first glance, the ζ avg seems to have the same trend as V avg , L v /T avg and F avg , the ζ avg for SFG 9.8  3 (0.493) is rather low despite having larger t 0 .This could be due to its smaller L 0 and its asymmetrical fractal geometry in comparison to other SFGs, where the absence of vortices shed from the second iteration bar at the central opening (Fig. 1g) may have prompted the local multilength-scale flow domain to be less crowded, and subsequently lower ζ avg .For the case of different σ, the ζ avg for σ = 0.270 and σ = 0.382 are slightly above expectation compared to σ = 0.493.Unequivocally, the central jet is weaker for lower σ insert with large A o , however the thinner second and third iteration bars of SFG 6.5  3 (0.270) and SFG 6.5 3 (0.382) create large number of small-scale disordered eddies in the flow, inducing thin-film damping.It is also important to note that greater t 0 brings about larger sized eddies impinging the thin-film surface, which yields substantial film vibration to overcome the thin-film damping instigated by the stronger central jet.Therefore, we still obtain considerably high F avg and V avg albeit with high ζ avg .
The local centreline mean velocities U at x/D h = 0.13, 2.50 and x peak of each grid are normalised against v in (U/v in ) and plotted in Fig. 8a, where with U x , U y and U z denoting respectively, the local centreline mean velocities in x, y and z-directions measured using hotwire anemometer.The corresponding ratio of the U gradient in the turbulence decay region to that in the generation region (m d /m g ) are depicted in Fig. 8b.
As shown in Fig. 8a, the centreline jet behind the central opening of all grids at x/D h = 0.13 have 1.2 < U/v in < 2.6, where the U decreases with x/D h towards v in further downstream.At x/D h = 2.5, the U of all SFGs, excluding SFG 9.8 3 (0.493) remain higher than v in by about 18% to 34% when compared to SSG and RG, which is consonant with the earlier study by Melina et al. (2016).Apart from SFG 9.8 3 (0.493) and RG, the U at x/D h = 0.13 increases with t 0 in view of the fact that the jet accelerates more rapidly when flowing through smaller sized central opening.In close proximity to the grid, SSG and SFGs with t r ≥ 6.5 have U that are higher than RG of equivalent σ.The overall U decrement rate is remarkably high for SSG and SFG 9.8 3 (0.493) with the former having the fastest decrement, SFG 6.5 3 (0.270) the slowest, and RG somewhere in-between.Inter- estingly, the overall U decrement rate appears to be interrelated with the FWHM of the V rms and δ rms peaks, namely the slower the U decrement, the more protracted the turbulence generation region, viz.larger FWHM.As such, the U decrement rate displays the same trend as those of V rms , δ rms and F rms , i.e., increase for larger values of t r [in good agreement with the results reported by Hurst and Vassilicos (2007) and Mazellier and Vassilicos (2010)] and σ but smaller N.These further shed light on the lower V rms , δ rms and F rms detected in the decay region when increasing t r (Fig. 4a, c).
It is clear from Fig. 8b that the rate where U declines in the turbulence decay region is comparable to that in the generation region (m d /m g ≈1) for SFG 4.0 3 (0.493) , SFG 9.8 3 (0.493) and SFG 6.5 2 (0.493) .For the remaining grids, m d /m g increases with N and when σ reduces.This could again be explained by the extended turbulence generation region as a result of larger x * , which prompts U to decline slower in the generation region.Nonetheless, it is worthwhile to note that the last measurement location (x/D h = 2.5) may not be adequately far away from the insert to fully capture the U decrement in the decay region.This is especially true for SFG 6.5 3 (0.270)where the decrease in U has yet to reach a pla- teau because of its large x peak , leading to an exceptionally high m d /m g .

Probability and Frequency Distributions of Piezoelectric Voltage Response
The cumulative probability of centreline V 2 at the x peak are plotted in Fig. 9a to compare the probability distributions of the random piezoelectric voltage response towards the turbulent flows generated by the best performing grids (grid that generates greatest flow fluctuation strength and turbulence forcing) for different N (SSG), t r [ SFG 9.8  3 (0.493) ] and σ [ SFG 6.5 3 (0.493) ] against RG.It is apparent that the grid with larger V rms and δ rms has a higher probability in obtaining large V 2 .SSG has the highest V 2 among the four inserts at x peak , of which, top 5 percentile have V 2 > 225 V 2 while top 1 percentile has V 2 > 403 V 2 .This is followed by SFG 9.8 3 (0.493)where 95% of V 2 are below 66 V 2 while top 1 per- centile has V 2 > 115 V 2 .For SFG 6.5 3 (0.493) , 99% of V 2 are not more than 47 V 2 and this is narrowed down to 5 V 2 for RG.
In Fig. 9b, we present the centreline V histograms for the corresponding grids at x peak .Clearly, the V induced by all the four inserts are normally distributed.One can see that the voltage response is most often found near zero values and that an increase in the voltage magnitude is generally accompanied by a decrease in the frequency of Fig. 9 a Cumulative probability of centreline V 2 at x peak of SSG(0.493) , SFG 9.8 3 (0.493) , SFG 6.5 3 (0.493) and RG(0.492) , as well as b the corresponding V histograms occurrence.All histograms are in good agreement with the cumulative probability distributions in Fig. 9a.The frequency in attaining high V is the rarest for RG with V ranging between − 4 and 4 V only.The V histogram becomes more spread out for similar σ SFG as supported by the wider V range of SFG 6.5  3 (0.493) (− 12 to 12 V) and SFG 9.8 3 (0.493) (− 18 to 18 V) of larger t r .Remarkably, the removal of smaller scale bars from the insert causes the histogram for SSG to have the widest V range (− 31 to 30 V).All inserts have almost the same percentage of outlier in their respective voltage distributions, ranging from 0.7% (RG) to 1.1% (SSG).
Both the cumulative probability distribution and histogram plots in Fig. 9 reveal the variation of voltage response with respect to grid conformation.The thin-film flapping induced by large t r SFG presents higher centreline voltage output over the RG of the same σ, with SSG outperforming the rest.Hence, SSG is deemed a promising turbulator in producing potent turbulence strength flow that is highly desirable for piezoelectric energy harvesting from turbulence-induced vibration.

Conclusions
The effects of fractal iteration N, thickness ratio t r and blockage ratio σ on the mechanical characteristics of fractal grid-generated turbulence were experimentally explored by means of PTFV.Measurements of the piezoelectric thin-film tip deflection δ and voltage response V at different grid-film distances x/D h along the centreline leeward of a SSG and various SFGs revealed increasing flow fluctuation strength in the turbulence generation region followed by a decay further downstream of the grids.Results showed that the distance from the grid where V rms and δ rms peak is a decreasing function of σ and t r .An increase in N along with a decrease in t r and σ leads to an elongated turbulence generation region as a result of larger wake-interaction length scale x * and slower local centreline mean velocity U decrement.
In contrast to the average dominant frequency of voltage signal f V (f Vavg ), the heighten in the average V rms (V avg ), equivalent lateral integral length scale L v /T (L v /T avg ) and millinewton turbulence forcing F rms (F avg ) downstream of the grids' centreline with increasing t r and σ but decreasing N is the consequence of the higher shedding intensity of lower frequency, larger scale energy-containing vortices from the thicker first iteration bar, which acts as the primary flow components in inducing thin-film flapping.For higher σ SFG, the accelerated multi-scale jets produced from the smaller grid openings, in addition to the more energetic vortex shedding from all three iterations bars of greater thickness further promotes the thin-film undulation.However, the higher: energy dissipation rate, U decrement rate and local deceleration experienced in the decay region of larger t r grid give rise to lower V rms , δ rms and F rms , particularly in the decay regime of SFG 9.8 3 (0.493) .Interestingly, the additional fractal scales of larger N grid results in less potent, faster rotating, smaller size vortices driving the thin-film vibration.
The flow generated by medium to high-blockage SFGs exhibit greater fluctuation strength and turbulence forcing than the RG, with the (V avg , F avg ) of SSG being twice of the best performing N = 3 SFG [i.e., SFG 9.8 3 (0.493) ], and (9×, 5×) higher than the RG of equiv- alent σ.Although the high turbulence intensity multilength-scale turbulence generated from the SFG is favourable for turbulent mixing and heat transfer enhancement (Teh et al. 2015;Skanthan et al. 2018;Chew et al. 2022;Hoi et al. 2022), the superior fluctuation strength (V rms = 7.6 V, δ rms = 3.1 mm) and turbulence forcing (F rms = 28 mN) demonstrated by the high-blockage SSG-generated flow at x peak enables SSG to serve as a promising turbulator for fluidic energy harvesting application.

Fig. 2
Fig. 2 Schematic of a side-view and b top-view of the experimental setup

Fig. 3 a
Fig. 3 a RMS of experimentally recorded piezoelectric output voltage V rms and tip deflection δ rms [top right shows the full width at half maximum (FWHM) of the respective peaks], as well as b the dominant frequency of voltage signal f V and thin-film flapping f δ at different streamwise distances x/D h along the centreline of SSG(0.493) , SFG 6.5 2 (0.493) , SFG 6.5 3 (0.493) and RG(0.492) [Symbol: piezoelectric measurement; Line: camera detection]; with the corresponding grid-induced c RMS excitation force acting on the piezoelectric thin-film F rms , and d equivalent lateral integral length scale L v /T against x/D h

Fig. 6 a
Fig. 6 a Downstream evolution of the centreline energy spectra for all grids; b top view of the respective spectra

Fig. 7 a
Fig. 7 a Average V rms , b f V , c L v /T, d F rms and e damping ratio ζ downstream of the grids centreline as functions of N, t r and σ

Fig. 8 a
Fig. 8 a Normalised local centreline mean velocity U/v in at x/D h = 0.13, 2.50 and the peak location x peak of all grids; b the corresponding ratio m d /m g of U gradient in the turbulence decay region (x peak → x/D h = 2.50) to that in the generation region (x/D h = 0.13 → x peak )

Table 2
Properties of the piezoelectric cantilever beam (Measurement Specialties 1999, 2015; Danesh-Yazdi et al. 2016; Lee et al. 2021) Subscript p represents PVDF film; subscript s the Mylar substrate.Superscript [m] denotes experimentally measured and [c] the calculated values