1. Introduction
Underwater mooring platforms (UMPs) are a class of underwater devices which are tethered to the seabed using mooring cables. UMPs have broad applications in both civil and military missions, such as underwater monitoring (oceanographic sensors), communication (acoustic communication nodes), and defense (mooring mines), with expected performance durations typically ranging from months to years. Nowadays increasing emphasis has been put on the improvement of UMPs. The greatest concern for the UMPs is the duration of their underwater operation, which is mainly determined by the total energy contained in the on-board batteries and the uninterrupted consumption of energy by the electronic devices. Although most UMPs are designed with low-power configurations now, the operation durations are still not satisfying due to the limited battery energy. Extending the operational life of UMPs can significantly reduce the cost for missions where a sustained presence is required, because of the high costs associated with retrieving, repowering, and redeploying remote systems.
UMPs stay in the ocean environment once deployed. The ocean has proved to be a promising source of renewable energies, including ocean surface solar energy [
1], wave energy [
2,
3], current energy [
4], thermal energy [
5], and osmotic energy [
6]. If these types of energy can be extracted to recharge the batteries of the UMPs, the underwater operation duration can be greatly extended. Although some of the above types of energy have been studied to power ocean devices, such as solar-powered autonomous underwater vehicles (AUVs) [
7], solar-powered ships [
8], and thermal-powered underwater gliders [
9,
10], they are not suitable for powering UMPs. UMPs are usually moored tens to hundreds of meters beneath the ocean surface, where solar energy and wave energy attenuate quickly and could hardly be utilized. As for the thermal energy, the UMPs, which are expected to perform fixed-point monitoring, cannot move vertically to pass through a thermocline to gain the thermal energy like the underwater gliders [
10].
UMPs are often deployed where ocean currents are consistently available. The kinetic energy stored in ocean currents provides an ideal energy source to recharge the UMPs. In previous studies, specially-designed ocean current turbines (OCTs) have been proposed for UMPs to generate electricity from the currents [
11,
12]. The use of the OCTs could extend the operation life of the UMPs, however, the OCTs have complex structures and show poor feasibility [
12]. Additionally, the OCT applies extra thrust and torque on the mooring system when rotating, making it harder for the UMP to reach motion stability [
13].
The flow-induced motion (FIM) conversion technology provides another way to utilize the ocean current energy. It is known that flow passing blunt bodies experiences unsteady separation and forms a wake of repeating swirling vortices, known as the Karman vortex street. The FIM generates kinetic energy from the surface-shed vortices and converts the kinetic energy into electricity. The FIM can be roughly divided into two groups, vortex-induced vibration (VIV) and galloping [
14]. VIV is caused by periodically-shed vortices which generate an asymmetric pressure distribution around the body and provide periodic forces that lead to a vibration in the body. Galloping is known as a dynamic instability that is induced in an elastic structure due to internal turbulence of the fluid or any reason which provides initial disturbance. A detailed review of the FIM devices can be found in [
15]. The idea of the VIV is first proposed by Bernitsas and Raghavan [
16], and further developed by the Marine Renewable Energy Laboratory at the University of Michigan [
17,
18,
19]. They built an energy converter named VIVACE (Vortex-Induced Vibration Aquatic Clean Energy) that uses a passive circular cylinder with upward-downward motion induced by vortex shedding. Over the past few years, comprehensive research has been done to reveal the mechanism of the flow-structure interaction and to enhance the power efficiency of the VIV device, including the influences of the cross-sections [
20], the spring stiffness [
21], damping ratio [
22], mass ratio [
23], passive turbulence control structure [
24], and multiple cylinders [
25].
UMPs usually have blunt sections, so they will possibly experience vortex shedding on the surface. Inspired from the above research, the shed vortices from the UMP contain a certain kinetic energy and could possibly be recovered to provide power for the UMP. A vortex-induced vibration energy converter (FIVEC) for UMPs is then designed.
Figure 1a shows the schematic diagram of the UMP and the FIVEC. The UMP is tethered to the seabed by a main cable, which is connected to a link of the UMP. The function of the link is shifting the anchored point on the UMP and increasing the motion stability. A thin plate is used to collect kinetic energy from the vortices. The working principle of the FIVEC is shown in detail in
Figure 1b. At a certain Reynolds number, repeating vortices will shed from both sides of the cylinder [
26]. The counter-rotating vortices will then change the load distribution on the plate periodically and drive the plate to vibration. This vibration motion will be transmitted to the inside of the FIVEC for power generation. A detailed description of the device will be given in
Section 2.
The FIVEC is more suitable for UMPs than OCTs due to the much simpler structure. However, its performance has never been studied before. This paper will focus on evaluating the characteristics of motion, force, power, and flow structure of the FIVEC at different design parameters with a computational fluid dynamic (CFD) method. It should be noted that the working principle of the FIVEC design is similar to the thrust foils in turbulent wakes [
27,
28]. Shao and Pan numerically modeled a thrust foil in the wake of a D-cylinder [
27], they found that the thrust of the foil was increased in the wake and explained the interactions between the vortices and the foil position. Similarly, Beal et al. investigated the hydrodynamics of a dead fish and a thin plate in the turbulent wake of a D-cylinder [
28]; they found that the dead fish and the thin plate could generate thrust, lift and, most interestingly, mechanical energy from the vortices in the wake of the cylinder. These studies prove the feasibility of energy extraction from the vortices and give inspiration to the current research.
5. Results and Discussion
The output power performance of the FIVEC is significantly influenced by the design parameters
c,
ks, and
Lp, among which
c determines the portion of the total power transformed into kinetic energy,
ks affects the natural frequency of the plate, and
Lp affects the rotational inertial and hydrodynamic forces. A previous study [
21] has found that the VIV reaches the highest power harnessing efficiency when the oscillation frequency (
fosc) is close to the device natural frequency (
fn) and that the frequency ratio between the
fosc and the
fn is not affected significantly by the total damping. Therefore, the effects of
k on the performance of the device are not considered and the specific value of
ks is set so that
fn is equal to the vortex shedding frequency according to: .
Since the appropriate values of c and Lp for the best power performance are unknown for this new design, this research takes the following steps to investigate the effects of the two parameters:
- (1)
Perform simulations on fixed plates (ks = ∞, c = ∞) with different Lp, and find the rough range of Lp where stronger hydrodynamic forces occur. In addition, obtain the approximate frequency of the varying forces (f0);
- (2)
Select a constant Lp in the rough range, calculate the rotational inertial according to Equations (5) and (6). Then define a constant ks to make fn, obtained by Equation (17), match with f0. Perform simulations for the specific Lp and ks and study the response of the system with different c; and
- (3)
Select an optimal c based on the results of Step (2), evaluate the system performance when Lp varies. In this simulation group, ks varies for each case so that the natural frequency fn is unchanged and matches with f0.
5.1. Results for the Fixed Plates
5.1.1. Instantaneous Forces
A series of simulations are carried out with
Lp varying from 0 (without a plate) to 2.4
D. The case of
Lp = 1.5
D is selected for the illustration of the instantaneous forces on both the UMP and the plate. The time history of the coefficients of forces, including the coefficients of torque, lift, and drag on the plate (termed as
Cm,
Clp, and
Cdp, respectively) and the coefficients of lift and drag on the UMP (termed as
Clu and
Cdu, respectively), are shown in
Figure 7a. To give a clear view, the curves of
Clp,
Cdp, and
Clu are magnified by two, four, and 100 times, respectively. It can be seen that all the curves change irregularly at the beginning of the simulation and reach oscillatory convergence after 30 s. It should be noted that all the simulations in this research last for at least 50 s for the coefficients to reach convergence and the data in the last three cycles are used to calculate the averaged coefficients. The curves of torque and lift oscillate around zero and have an averaged zero value, while those of drag are biased. The
Cdp curve offsets to the negative direction which means that the plate generates thrust rather than drag in the wake of the UMP. The plate is located behind the UMP and is influenced by the vortices shed from the UMP surface. The velocity direction of the vortices which contact with the plate surfaces points left, therefore, the plate generates negative drag, this phenomenon was also observed in an experiment of dead fish in vortices [
28]. The
Cdu remains at a high level of about 0.60 and shows a very slightly oscillation, as shown in the zoomed-in plot.
The fast Fourier transform (FFT) spectral analysis for the force coefficients are shown in
Figure 7b. It is observed that the dominant frequencies are 0.4272 Hz for
Cm,
Clp, and
Cdp, and is doubled (0.8545 Hz) for
Cdp and
Cdu. This indicates that the drag experiences two cycles in a whole plate oscillation cycle. Another finding is that
Cdu has a second dominant frequency, 0.4272 Hz, which is the equal to the dominant frequencies of
Cm,
Clp, and
Cdp. The reason for the two dominant frequencies of
Cdu can be explained as follows: The first dominant frequency, 0.8545 Hz, is the result of the periodically shed vortices. The UMP in turn shed two contra-rotating vortices in a variation period (Figure 13, the flow structures will be discussed in depth in
Section 5.2 and
Section 5.3), each vortex has the same contribution to the drag, therefore, the frequency is twice that of the lift. The second dominant frequency, 0.4272 Hz, is the result of the interacting forces between the plate and the UMP and has the same frequency with that of the lift and torque.
5.1.2. Plate Length Effects
This section investigates the effect of length of the fixed plates on the force characteristics. To clearly characterize the forces on the UMP and the plate, the root mean square (RMS) of the lift and torque, and the averaged drag, are chosen in the following contents. The variation of the forces and their dominant frequency are shown in
Figure 8.
The
Clp increases with
Lp/D from
Lp/D = 0 to
Lp/D = 1.0 and remains at a high level in the range of 1.0 <
Lp/D < 2.0. As
Lp/D further increases, a abrupt drop in
Clp occurs. The
Clu is greatly reduced in the region of 0 <
Lp/D < 0.4, then increases slightly between 0.4 <
Lp/D < 1.0. After that
Clu decreases with the increased
Lp/D and is almost zero after
Lp/D = 2.2. Taking the total lift into consideration, the curve of the lift coefficient can be divided into three regions: a low suppression region (0 <
Lp/D ≤ 1.0), a amplification region (1.0 <
Lp/D ≤ 1.8), and a high suppression region (
Lp/D > 1.8). In the two suppression regions, the total lift of the system is reduced, which means that the UMP vibration motion could be possibly suppressed. While in the amplification region, the total lift is higher than that of the bare UMP (
Lp/D = 0), which means that a possible high output power of the blade could be obtained in this region. Therefore, a
Lp/D of 1.0 is selected for the damping analysis in
Section 5.2.
Since the lift and the torque has the same dominant frequency (
Figure 7b), FFT frequencies of the
Clp, alone, are used. The dominant frequency in
Figure 8b is 0.5249 Hz for the bare UMP, and remains unchanged with a value of 0.4587 Hz between 0.1 ≤
Lp/D ≤ 1.4. A nearly linear decrease exists in the region of 1.5 <
Lp/D ≤ 2.2. Although there must be a difference between the dominant frequencies of the fixed plate and the oscillating plate, the frequency of the fixed plate provides a guide for the design of the natural frequency of the oscillating plate. Since the dominant frequencies are very close when
Lp/D varies from 0 to 1.4 and that the frequency may change slightly in the cases of oscillating plates due to the reaction of the plate on the UMP, a constant natural frequency of 0.5000 Hz is chosen for all the simulations in this study.
The
Cdp shown in
Figure 8c is rather smaller compared with the
Cdu. This is because the plate is mainly influenced by the vortices rather than the freestream.
Figure 9 presents the pressure contours and streamlines for three different fixed plates. It can be seen that the magnitude of the velocity around the plate is much lower than the freestream. In addition, the direction of the vortices on the plate surface is almost parallel to the plate, meaning that the pressure drag can be ignored. The
Cdu decreases nearly linearly with increase of
Lp/D. This is because the UMP with a plate has a more streamlined shape. The drag coefficients of streamlined shapes are smaller than that of the UMP, which has a circular section.
The
Cm shown in
Figure 8d rises with
Lp/D from
Lp/D = 0 to
Lp/D = 1.7. It is easy to understand that the enlarged plate has a larger surface area and can generate a greater hydrodynamic torque. However, as
Lp/D further increases, the
Cm drops abruptly and remains at a low level after
Lp/D = 2.2. This can be explained by the pressure contours in
Figure 9. In the pressure contours of
Lp/D = 2.2, the plate is too long and it is not easy for the reattached vortices to shed from the tip. As a result, both sides of the plate are covered by vortices and the pressure distributions are almost symmetrical. The torque on the plate is then greatly reduced.
5.2. Effects of Damping on the System Performance
Simulations are carried out with constant Lp (Lp = 1.0D) and ks (1020 N/rad, resulting with a natural frequency of 0.5000 Hz). The effects of c are evaluated over a range varying from 0 to 1000 Ns/rad.
5.2.1. Oscillation Amplitude and Power
Figure 10 shows the variation of the oscillation amplitude,
θa, and power coefficient,
Cp, with respect to damping.
θa decreases with the increase of
c, due to the increased damping force which hinders the rotation of the plate. The maximum
θa occurs at zero damping, with a value of 10.38°. Although the amplitude is higher for a smaller damping, the plate power is also smaller because it is proportional to
c, according to Equation (7). As for
p, it rises with
c firstly to
c = 150 Ns/rad, and then decreases with further increased
c. The
Cp has a peak value of 0.0482. This indicates that the suitable damping for high power output lies near
c = 150 Ns/rad. Therefore,
c = 150 Ns/rad is chosen for the plate length analysis in
Section 5.3.
5.2.2. Performance of Forces
The RMS of the lift and torque coefficients, together with the averaged drag coefficients, are shown in
Figure 11.
Figure 11a suggests that the
Clp increases with
c for the whole examined range, but the growth rate decreases. It is known that the lift drives the plate while the damping force hinders the motion of the plate, so the plate with a higher damping is more capable to remain its position in flow. Therefore, the plate can generate a higher lift at a higher damping. Since the lift is the main force that generates the torque, the
Cm shows a similar trend with the
Clp, both increasing with
c but decreasing in growth rate (
Figure 11c).
The
Clu decreases with
c at first and reaches its minimum at
c = 300 Ns/rad, then it increases slowly with
c. An interesting finding is that the
Clu is greatly reduced at
c = 300 Ns/rad, causing the total lift coefficient to drop to 0.3549, which is 4.83% lower than that of a bare UMP (0.3729). To explain the variation of
Clu, the instantaneous lifts of both the UMP and the plate are shown in
Figure 12. There is a phase difference between the
Clu and the
Clp which makes the
Clp lag behind the
Clu. It is clearly shown that the phase difference increases with the damping. Since
Clp is the main force that generates torque and drives the rotation of the plate, the phase difference leads to an asynchronization in the UMP lift and the plate motion.
Figure 12 also shows the pressure contours around the UMP at
t = 40.96, which corresponds to the maximum
Clu. It is clear that the position of the plate differs at different damping. The rotation of the plate changes the flow near the suction (downstream) side of the UMP and affects the pressure on the UMP, as marked in the dotted circles. At
c = 300 Ns/rad, the pressure on the UMP is almost symmetry about the centerline, so the
Clu is greatly reduced. While for the other cases, the pressure on the top caused by the vortex shedding is lower than that on the bottom, so the
Clu is relatively higher.
The Cdp is negative and near zero for all damping, this means that the plate could generate a small thrust in the wake. The total drag coefficient decreases with damping slightly at the beginning and then remains nearly constant after c = 400 Ns/rad. Generally, damping has only a slight effect on the drags of the UMP and the plate, because the amplitude of the plate is relatively small.
5.2.3. Wake Structures
In order to show the influences of system damping on the wake of the UMP, the time-varying coefficients of torque,
Cm, and the angular position,
θ, are presented in
Figure 13. Three different
c, including 0 (free motion), 150 Ns/rad (peak
Cp) and 300 Ns/rad (high damping), are selected. It can be seen that the
Cm curves show stronger oscillations at high damping, while the
θ curves shows the opposite trend. Another finding is that a phase difference exists between
Cm and
θ, making
Cm lags behind
θ. This phase difference becomes more obvious at higher damping. Contours of vorticity are plotted at four typical points in a cycle of the
Cm curve, as shown in the curves. The wake structure are similar for different damping. A regular 2S vortex pattern (two single vortices shed per cycle) is observed for all simulation cases. Two vortices are shed per cycle of oscillation, the clockwise rotating vortex by the top shear layer and the counter-clockwise rotating one by the bottom shear layer. Both the vortex shedding frequency and the distance between two adjacent vortices are equal. Vortices shed from both sides of the UMP and then reattaches to the surface of the plate forming a thick shear layer flow. The shear layer extends to the tip of the plate and then separates and forms contra-rotating vortices. This separation delay phenomenon is also observed in a previous study [
34].
5.3. Effects of Plate Length on the System Performance
Ten sets of simulations are carried out at constant system damping (150 Ns/rad), natural frequency (0.5000 Hz), and inflow velocity (1 m/s). The effects of Lp are evaluated over a range varying from 0.2D to 1.8D.
5.3.1. Oscillation Amplitude and Power
Figure 14 shows the variation of the oscillation amplitude,
θa, and the power coefficient,
Cp, with respect to
Lp/D. θa increases with
Lp/D quickly from
Lp/D = 0.2 and reaches the peak at
Lp/D = 0.9 with a value of
θa = 7.82°. After that
θa drops approximately linearly with
Lp/D. The
Cp shows a similar trend with
θa. The peak
Cp is 0.0520 at
Lp/D = 0.7, which means that the plate can generate an averaged power of 13.86 W from the flow. Although the peak
Cp is smaller than conventional VIV devices (0.37 by Ding et al. [
14]) and ocean current turbines, this result is still satisfactory. Considering that the standby power of common underwater mooring buoys are on the order of tens to hundreds of milliwatts, the averaged power produced by the FIVEC could compensate for that.
5.3.2. Performance of Forces
In order to clearly show the performance of forces at different plate lengths, the RMS of the lift and torque coefficients, together with the averaged drag coefficients, are shown in
Figure 15. Although the
Clp increases with
Lp on the whole, there is a relatively flatter region (0.4 <
Lp/D < 1.2) where
Clp shows small fluctuation. The
Clu decreases with
Lp at first and reaches its minimum at
Lp/
D = 0.7, then it increases slowly. The
Cdp is negative and decreases with the increase of
Lp for all cases. The
Cdu is increased slightly when
Lp/D < 0.7 and drops nearly linearly with further increased
Lp. While for the
Cm, it shows a similar trend with the total
Cl, both have little change within the region of 0.4 <
Lp/D < 1.2.
In order to show the influences of plate length on the performance of the device, the time-varying coefficients of lift and amplitude are presented in
Figure 16. Typically, four different lengths of plates, including 0 (bare UMP), 0.2
D, 0.7
D (corresponding to the peak
Cp), and 1.4
D, are selected and compared. Despite for the amplitude, the shapes of the curves of
Clp and
Clu are also different for different plate lengths. For
Lp/D = 0 and 0.2, the curves are like rough sinusoidal curves, while for the other two cases, the curves are more complicated. Additionally, a phase difference exists between
Clp and
Clu, which is larger for longer plate lengths. Contours of vorticity are plotted at four typical points in a cycle of the
Cm curve, as shown in
Figure 14. The wake shows a regular 2S vortex pattern for all cases. This indicates that the plate have little effect on the far wake of the UMP. However, the near wake of the UMP is greatly influenced. The vortices shed from the UMP reattaches to the plate. The shear layer is then extended to the tip of the plate. The shifts in the shear layer are the reason for the variation of the forces.
To illustrate how
Lp influences the forces on the UMP and the plate clearly, the velocity contours with streamlines and the pressure contours near the FIVEC are shown in
Figure 17. The plate has three main effects on the near wake flow of the UMP. Firstly, the plate works as a barrier and prevents the vortices from rolling from one side to the other side of the UMP. It is clear to see that a vortex shed from the top of the UMP and rolls to the bottom, covering most area of the suction (downstream) surface at
Lp/
D = 0. This vortex significantly influences the flow near the bottom of the UMP and the pressure is recovered. The pressure difference between the top and the bottom side is then increased and, finally, the lift is increased. For the other cases with plates, the top vortices cannot easily move from the top to the bottom and the pressure on the bottom side of the UMP is hardly changed. This helps reduce the lift of the UMP. Secondly, the plate helps extend the length of the shear layer and delaying the shedding of vortices. The vortices separate from the tip of the plate and help extend the low-pressure area where separation occurs. Therefore, the pressure on the suction side of the UMP is increased and the drag is reduced. Thirdly, the push-and-pull movement of the plate causes a local flow near the suction side of the UMP, this flow is stronger at larger lengths of plate. Taking the case of
Lp/D = 1.4 for example, the plate is rotating counterclockwise, pushing the flow on the top and drawing the flow at the bottom. Therefore, the pressure on the bottom side of the UMP is much smaller than the other cases. The effects of local flow is the reason for why
Clu is increased at higher
Lp/D.