3.1. Load resistance analysis
In order to ensure that the PEHs in this paper can achieve the best energy harvest effect, we first explore the influence of load resistance on PEHs. Considering the impedance of vibrating piezoelectric ceramic, we analyzed the energy harvest effect of BARE-UPEH by changing the load resistance value. The experimental wind speed is 7.5m/s. In terms of spacing ratio, we chose three spacing ratios L/D: 1.5, 3.0 and 5.0, to explore the maximum voltage and harvested power of PEHs. It can be seen from Fig. 5 that the root mean square voltage (RMS Voltage) generated by BARE-UPEH at three spacing increases with the increase of load resistance, while the power reaches the peak at 1.8×106MΩ. At the spacing ratio of 1.5, 3.0 and 5.0, the maximum power generated by BARE-UPEH is 0.74629mW, 0.42681mW and 0.52156mW respectively. It can be concluded that the optimal load resistance of PEHs in this paper is 1.8×106 MΩ, both UPEH and DPEH circuits in the following experiments are conducted at this load resistance value.
3.2. Effect of spacing ratio on amplitude of UPEH and DPEH
As we all know, spacing ratio (L/D) is critical to the flow pattern and wake structure of two tandem PEHs in a cross flow. In other words, the aerodynamic coefficient depends on the spacing ratio [40]. As shown in the Fig. 6 (a) (b) (c), whether to install D-type bionic fins and change the installation position of D-type bionic fins, the maximum amplitude response of UPEH occurs at the spacing ratio of 1.5. This indicates that the D-type bionic fin may enhance the flow acceleration of fluid in the gap area between the two PEHs, thus generating the gap driving force and increasing the fluid force on UPEH [41]. As the spacing ratio increases, the gap driving force and flow acceleration in the gap area decrease, and the amplitude of the corresponding UPEH decreases. It shows that adding D-type bionic fins in tandem PEHs does not change this rule. It is worth mentioning that the maximum amplitude response of BARE-UPEH and U-UPEH occurs at 10 m/s, while D-UPEH occurs at about 5 m/s, indicating that a “lock-in zone” similar to vortex-induced vibration may occur near 5 m/s in D-DPEH, which will be analyzed in the voltage response.
As shown in the Fig. 6 (d) (e) (f), BARE-DPEH maximum amplitude is generated at 1.5 spacing ratio, which is consistent with the results of literature [42]. After installing D-type bionic fins upstream, the coupling effect of U-UPEH on downstream U-DPEH is enhanced, and the maximum amplitude produced by three-spacing U-DPEH under the experimental wind speed is close. For D-DPEH, at the spacing ratio of 5.0, due to the large spacing ratio, the vortex generated by D-UPEH has little impact on D-DPEH. D-DPEH is similar to the state of a single PEH. The amplitude response of D-DPEH is largely determined by the promotion of D-type bionic fin on vortex shedding. The maximum amplitude generation at the spacing ratio of 5.0 indicates that the enhanced intensity of D-type bionic fin on vortex shedding can break through the restriction of spacing ratio to amplitude response. An interesting phenomenon is found in this paper, from the perspective of the degree of influence of installing D-type bionic fins, on the amplitude of DPEH, the installation of D-type bionic fins directly downstream has a greater impact on the corresponding amplitude of DPEH than the installation of D-type fins upstream. To sum up, D-type bionic fin can make the maximum amplitude of DPEH change from small spacing to large spacing, and expand the application range of PEHs.
3.3. Output voltage analysis
In this paper, the voltage is divided into three regions: low velocity region 0 ~ 2.5m/s, medium velocity region 2.5 ~ 7.5m/s, and high velocity region 7.5 ~ 10m/s. As described in reference, the voltage trend is consistent with the amplitude trend, and the phenomenon is consistent [40]. This paper analyzes the three vibration modes and change trend based on voltage.
As shown in the Fig. 7 (a) (b), when the L/D = 1.5 (mode 1), BARE-UPEH and BARE-DPEH vibrate at the wind speed of 3 m/s and 1.5 m/s respectively, and then the voltage response slowly increases with the increase of wind speed, showing a typical wake induced galloping (WIG) phenomenon [43]. After adding D-type fins to the upstream PEH, the voltage of U-UPEH starts to respond at U = 3.5m/s, the voltage response increases first and then decreases, and the voltage response basically remains unchanged between 3.5 ~ 4.5m/s wind speed, similar to the situation of traditional VIV "locking zone". However, the difference is that the cut-in speed Ucut is different from the traditional vortex-induced vibration [44]. Traditional VIV occurs at a very low Ucut. Tang et al. [45] also observed a similar phenomenon and defined it as "coupled vortex-induced vibration". Interestingly, the voltage response of U-UPEH increases again when the wind speed is 7m /s, and the subsequent voltage response increases with the wind speed, showing the phenomenon of WIG. The RMS voltage generated at 10m/s reaches 84.15V, which is much higher than BARE-UPEH and D-UPEH. This may be caused by the D-type bionic fin increasing the velocity of surface fluid on both sides of U-UPEH. The fifth of Section 3 will further analyze the flow field to illustrate this phenomenon. The above phenomenon shows that adding D-type fins upstream of PEH can achieve complete coupling of the two vibration types of “coupled vortex-induced vibration” and WIG by controlling the wind speed, thus increasing the voltage harvested by UPEH. U-DPEH achieves better voltage response than BARE-DPEH in the "locked zone" of "coupled vortex-induced vibration", which shows the possibility that D-type bionic fin can cause greater voltage response with lower wind speed. The voltage of D-UPEH starts to respond at the speed of 3.5m/s, and also produces the phenomenon of "coupled vortex-induced vibration". The voltage response first increases and then decreases with the wind speed. At the same time, D-DPEH shows the phenomenon of WIG, and the voltage slowly rises with the wind speed.
As shown in the Fig. 7 (c) (d), when the L/D = 3.0 (mode 2), BARE-UPEH has almost no voltage response, indicating that BARE-UPEH does not vibrate at this time, which is consistent with the amplitude response mentioned above. The voltage response of U-UPEH and D-UPEH increases sharply in the high wind speed region, and the overall response is still at a low value, reaching the maximum value when the wind speed is 9m/s, and then the voltage response decreases with the increase of wind speed. Ucut of the U-DPEH is the largest, the voltage response rises slowly with the increase of wind speed and increases sharply when the wind speed is 6m/s. The overall phenomenon is WIG. The installation of D-type bionic fin makes the UPEH voltage response change from nothing to something in the high wind speed region. In the high wind speed region, the voltage response of U-DPEH is better than BARE-DPEH and D-DPEH, and the maximum RMS voltage can reach 35.10V. D-DPEH has the lowest Ucut and shows better voltage harvest capability in the medium and low wind speed areas. The voltage of D-DPEH increases first and then decreases in the high wind speed areas. The voltage response mode is chaotic. In this mode, the vortex shedding from the D-UPEH cannot be closed, and the fluid flow in this interval is the most complex [46].
As shown in the Fig. 7 (e) (f), when L/D = 5.0 (mode 3), the voltage response of BARE-UPEH and D-UPEH is relatively small within the experimental wind speed, the Ucut of U-UPEH is 8 m/s, and then WIG occurs. At the wind speed of 10 m/s, the RMS voltage is 6.169V. BARE-DPEH and D-DPEH show WIG phenomenon, but Ucut of D-DPEH is 4.5m/s, compared with BARE-DPEH and U-DPEH, the Ucut is lower, and the voltage response of D-DPEH is greater in the whole experimental wind speed range. This may be due to the installation of D-type bionic fins, which leads to the enhancement of D-DPEH vortex shedding [47]. U-DPEH also has the phenomenon of “coupled vortex-induced vibration”. The maximum RMS voltage generated in the “locking zone” is 35.77V, which is higher than BARE-DPEH at the same wind speed. Interestingly, installing D-type bionic fins will increase the natural frequency of the PEHs as it is shown in Fig. 4. Installing D-type bionic fins in the downstream energy harvester will theoretically increase the Ucut of D-DPEH, but D-DPEH has voltage response before D-UPEH and BARE-DPEH at the spacing ratio of 3.0 and 5.0, which indicates that installing D-type bionic fins in the middle and large spacing ratio may promote the vortex at the tail of D-DPEH to fall off ahead of time and reduce the Ucut. So as to expand the bandwidth of the voltage harvested by the PEH.
3.4. Evaluation of the output power enhancement of the proposed D-type bionic fin
Through various studies on the installation position and three modes of D-type bionic fins so far, we can determine that adding D-type bionic fins on upstream and downstream under the experimental wind speed can undoubtedly bring the maximum power increase, and can determine the single PEH with the best performance within the entire wind speed range considered in this study: U-UPEH, mode 1 (L/D = 1.5), wind speed 10m/s. As shown in Fig. 8 (a), the maximum power of U-UPEH is increased by 392.38%, and the maximum power of D-UPEH is increased by 116.19% compared with BARE-UPEH. Compared with mode 1 (L/D = 1.5), when mode 2 (L/D = 3.0) and mode 3 (L/D = 5.0), the power generated by UPEH is relatively small, and the output power of BARE-UPEH is almost zero. However, installing D-type bionic fins at upstream and downstream can still increase the power of U-UPEH.
As shown in Fig. 8 (b), the maximum power of U-DPEH is less than BARE-DPEH when L/D = 1.5 within the experimental wind speed range. However, considering that the excitation wind speed Ucut=7.5m/s of the maximum power generated by U-DPEH is lower than the excitation wind speed Ucut=10m/s of BARE-DPEH, and this maximum power is generated under "coupled vortex-induced vibration", the above voltage analysis also confirms this view, so U-DPEH also shows the potential of harvesting energy at lower wind speed. The maximum power of U-DPEH and D-DPEH at mode 2 (L/D = 3.0) and mode3 (L/D = 5.0) is also improved compared with BARE-DPEH, which is 13% and 8.8%. If we consider a single PEH, we recommend the U-UPEH in Mode 1 and if two PEHs harvest energy at the same time, we recommend the U-UPEH and U-DPEH in Mode 2.
3.5. Flow analysis
An interesting phenomenon observed in this paper is shown in Fig. 7(a), (b), in the case of mode 1 (L/D = 1.5), when the wind speed is U = 10m/s, U-UPEH shows a much higher energy harvest capability than BARE-UPEH, while the energy harvest capability of D-DPEH is lower than BARE-DPEH. Therefore, we carried out two-dimensional unsteady numerical simulation at this wind speed to further understand the physical differences between UPEH and DPEH in the process of flow-induced vibration when installing D-type bionic fins upstream under mode 1 (L/D = 1.5).
The length and width of the simulated flow field are 60D and 40D respectively. The PEHs are placed at 15D away from the entrance and middle position in the width direction. Unstructured triangular mesh is introduced to divide the calculation area. In order to optimize the calculation, the calculation area is divided into three areas according to the distance from the PEHs, and the grid is further refined near the wall of the PEHs to improve the simulation accuracy. The maximum grid size of these three areas is 0.12D, 0.08D and 0.04D respectively. The total number of meshes is about 600000, and the calculation time step is 0.001s. The inlet and outlet boundary conditions are 10m/s uniform wind speed and pressure outlet flow. The horizontal boundary of the fluid field and the boundary of the two PEHs are assumed to be sliding walls and non-sliding walls, respectively. Since the Reynolds number Re = 3.6×105, there is a turbulent wake behind the DPEH and the flow area behind the DPEH is an intermediate subcritical region with periodic shedding vortex, although the Reynolds average Navier-Stokes (RANS) method can greatly reduce the computational cost, it ignores the turbulence pulsation information, which makes the computational data in the turbulent wake distorted. Compared with the RANS method, large eddy simulation (LES) can capture the periodic shedding of the entire eddy by directly solving the large-scale eddy. After modeling, it can describe the complex small-scale turbulent eddy at different locations in time and space in detail. Considering the calculation accuracy and cost, LES is introduced to calculate the turbulent fluid flow around the proposed PEHs by CFD.
As shown in Fig. 9, when the incoming flow meets PEHs, a relatively low speed zone will appear in front of the UPEH, and then the air flow will flow along both sides of the UPEH and form a high speed zone on both sides of the UPEH and continue to flow between the UPEH and the DPEH to form a low speed zone, so then the high speed zone will continue to form on both sides of the DPEH, and the vortex shedding is caused by the speed difference, so then the vortex will be generated alternately in the S shape at the tail of the DPEH.
As shown in Fig. 9 and Fig. 10, the U-UPEH has higher flow velocity and greater pressure difference on both sides compared with BARE-UPEH, which indicates that the lateral excitation force generated by the pressure difference on both sides of U-UPEH is greater, so it is a good explanation for the amplitude response ability and voltage harvest ability of U-UPEH shown in the experiment. This also verifies the effect of manufacturing a microstructure similar to the bionic non-smooth surface structure mentioned in the design of D-type bionic fins to enhance the fluid vortex near the body surface, thus reducing the fluid resistance and increasing the flow rate. Compared with U-DPEH, the intensity of vortex generated by the wake of BARE-DPEH is larger, and the frequency of vortex shedding is faster. After BARE-DPEH, it can maintain a longer vortex area and a higher vortex intensity. At the same time, there is a greater velocity difference and pressure difference on both sides of BARE-DPEH, which fully verifies that Fig. 7 (b) BARE-DPEH harvests voltage greater than U-DPEH at the wind speed of 10m/s.
By integrating the pressure distribution around the UPEH and DPEH cross sections, the lift coefficient that changes with time is calculated. As shown in Fig. 11, U-UPEH and BARE-DPEH show higher lift coefficients than BARE-UPEH and U-DPEH, which also verifies that the amplitude response and voltage response discussed above are consistent with the experimental data. Since the vibration of the PEHs is determined by the lift coefficient and the drag coefficient, this paper only discusses how to improve the performance of the lift coefficient by adding D-type bionic fins.