Next Article in Journal
A New Class of Edge Filter Based on a Cross-correlation-like Equation Derived from Fractional Calculus Principles
Next Article in Special Issue
Development of Dehydrogenation System for Liquid Organic Hydrogen Carrier with Enhanced Reaction Rate
Previous Article in Journal
Open DGML: Intrusion Detection Based on Open-Domain Generation Meta-Learning
Previous Article in Special Issue
Electrochemical Deposition of Bismuth on Graphite Felt Electrodes: Influence on Negative Half-Cell Reactions in Vanadium Redox Flow Batteries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Efficiency of Two-Degree-of-Freedom Galloping Energy Harvesters with Two Transducers

Institute of Applied Mechanics, Faculty of Mechanical Engineering, Poznan University of Technology, 60-965 Poznań, Poland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(13), 5427; https://doi.org/10.3390/app14135427
Submission received: 9 May 2024 / Revised: 17 June 2024 / Accepted: 18 June 2024 / Published: 22 June 2024

Abstract

:
This paper examines the energy efficiency of three variations of the two-degree-of-freedom transverse galloping energy harvester. These variants differ in the number and placement of electromechanical transducers. By utilizing the harmonic balance method, the limit cycles of mathematical models of the devices were determined. Analytical expressions derived from the models were then used to formulate the efficiency of the systems. It was demonstrated that efficiency depends on flow speed and can be comprehensively characterized by the following criteria parameters: peak efficiency, denoting the maximum efficiency of the system, and high-efficiency bandwidth, which describes the range of flow speeds within which the efficiency remains at no less than 90% of peak efficiency. The values of these parameters are heavily reliant on two other parameters: the speed at which the system achieves peak efficiency, referred to as the nominal speed, and also the flow speed at which the system undergoes Hopf bifurcation, namely the critical speed. Comparative analysis revealed that only the device equipped with two electromechanical transducers can potentially outperform a simple one-degree-of-freedom system. For selected parameters, this gain reached nearly 10%.

1. Introduction

The ongoing energy crisis and burgeoning concept of the Internet of Things (IoT) [1] are continually motivating factors for researchers from various scientific fields to develop newer technologies that enable cheaper and more efficient generation of usable energy. Among the most innovative of these technologies are those that allow energy recovery from vibrations induced by flow. A commercial example of this is the Vortex Bladeless power plants [2], whose vibrations are excited by von Karman vortices. According to the manufacturer, the cost of energy production using these turbines is only 70% of the cost associated with traditional wind turbines of comparable size [2]. Moreover, their sleek, bladeless design allows for closer placement than standard turbines, facilitating more efficient use of wind farm space.
The discussed advantages of vortex-induced vibration energy harvesters may be partially overshadowed by a significant limitation. Due to their reliance on resonance with the vortex frequency, these devices can generate energy only within a narrow flow speed range. An alternative type of power plant without this limitation is the galloping energy harvester (GEH). Transverse galloping [3] is a phenomenon where vibrations are induced due to the flow separation from a bluff body with an appropriate shape, leading to negative, speed-dependent aerodynamic (hydrodynamic) damping. When the flow speed exceeds a critical value, the negative aerodynamic (hydrodynamic) damping surpasses structural damping, resulting in the overall damping in the system becoming negative. In this situation, the damping will supply energy to the system rather than dissipate it, causing the system to become unstable.
The key advantage of using transverse galloping as an excitation mechanism for energy harvesters is that, beyond the critical speed, the amplitude of oscillation increases indefinitely with the flow speed and remains independent of the frequency or the presence of von Karman vortices.
Early studies [4] explore its potential for energy harvesting, emphasizing the crucial role of aerodynamic coefficients in its efficiency. Subsequent research [5] delves into analytical, numerical, and experimental analyses of GEHs, indicating the detrimental impact of structural damping on the performance of the device.
Due to the critical importance of the aerodynamics of the flowing body, many subsequent works have been devoted to examining the efficiency provided by GEHs of various shapes. Building on this trend, it has been demonstrated, for instance in [6], that rhomboidal shapes exhibit a significant propensity for galloping. The article [7] is dedicated to the study of the properties of a device equipped with a bluff body in the shape of various variants of an isosceles triangle. Furthermore, research [8,9] shows that adding a stream splitter favors power generation. The studies [10,11,12] have also demonstrated the benefits of employing non-typical shapes. According to the reports presented in [13,14,15,16], the efficiency of a GEH can also be increased by giving its surface an appropriate geometric metastructure.
Another extensively investigated area in the development of GEHs is the examination of the impact of nonlinear elasticity on its performance. Studies [17,18] propose a method for analyzing the effect of the nonlinear elasticity coefficient on the characteristics of Hopf bifurcation. The study [19] focuses on improving the energy harvesting performance of galloping systems through the use of magnetic coupling. Conversely, the work [20] presents modeling and experimental investigations of asymmetric distances using magnetic coupling based on a galloping piezoelectric energy harvester, aiming to enhance the stability of the generated power. In the articles [21,22], it has been demonstrated that incorporating nonlinear stiffness into the system can significantly increase its efficiency. Research [23] demonstrates that a GEH composed of two oscillators coupled by a nonlinear magnetic interaction can exhibit a lower critical speed. Furthermore, an appropriate selection of the distance between the magnetic segments can contribute to an increase in the system’s efficiency. A similar coupling, but in reference to a single-degree-of-freedom system, was discussed in [24]. This approach resulted in a tristable system, which exhibited an efficiency over 35% greater than that of a standard GEH.
According to the literature, enhancing the energy extraction efficiency of a GEH can be also achieved by extending its mechanical structure with an additional degree of freedom. In [25], two variants of GEHs with two degrees of freedom (GEH2Ds), differing in the location of the applied aerodynamic force, were compared to a single-degree-of-freedom device. Experimental comparisons of the electric voltage generated by the variants were presented in [26] and extended to numerical studies of systems with up to three additional masses in [27]. Comprehensive analytical studies of the voltage generated by a GEH2D, considering the possibility of two different vibration modes, were presented in [28]. A two-mass system, whose segments are magnetically coupled, was discussed in [29]. It was demonstrated that such a system can exhibit both a reduced critical speed and enhanced efficiency.
The literature review clearly indicates that extending the mechanical structure to include an additional degree of freedom is seen as a potential way to increase the efficiency of a GEH. However, the literature does not discuss the fact that within the family of two-mass systems, three device variants can be distinguished, differing in the placement of the electromechanical transducer: it can be located between the stationary base and the first mass, between the masses, or in both of these places.
The aim of this work is therefore to analytically examine the efficiency of all three GEH2D variants, compare them, and assess whether any of them can exhibit greater efficiency than the device of basic design.

2. Efficiency of Reference Variant

While the primary focus of this study lies in analyzing the efficiency of systems with two degrees of freedom, understanding the characteristics of this class of devices can be facilitated by comparing them with a system possessing well-established properties—a linear device with one degree of freedom, hereafter referred to as the reference device (Figure 1). Therefore, it is justified to conduct an analysis of the efficiency of this system, particularly as we aim to highlight a certain intriguing feature of this variant, which has not been previously addressed in the literature. According to articles [3,4], the mathematical model of a GEH takes the form shown in Equation (1) with parameters detailed in Table 1.
m ^ z ¨ + c ^   z ˙ + k ^ z θ ξ = F L ^ u ^ = 1 2 ρ ^ u ^ 2 h ^ a 1 z ˙ u ^ + a 3 z ˙ u ^ 3 ,
C p ^ ξ ˙ + ξ R ^ + θ ^ x ˙ = 0 .
The efficiency η of such a system is defined as the ratio of the average power generated by the device P g ^ over a time period equal to the period T ^ of the voltage function ξ(τ) to the surface flow power density P f ^ :
η L = P g ^ P f ^ = 1 T ^ 0 T ^ ξ 2 τ R ^ d τ 1 2 ρ ^ h ^ u ^ 3 ,
By introducing the following dimensionless quantities:
y = z h ^ ,   v = ξ θ ^ m ^   h ^   ω ^ n 2 ,   ω ^ n = k ^ m ^ ,   c = c ^ m ^ ω ^ n ,   u = u ^ h ^ ω ^ n ,   ρ = ρ ^ h ^ 2 2 m ^ ,   r = C p ^ R ^ ω ^ n , κ = θ ^ 2 C p ^   m ^    ω ^ n 2 ,   t = τ ω ^ n ,
the device dimensionless mathematical model can be written as
y ¨ + c y ˙ + y v = ρ a 1 u y ˙ + a 3 y ˙ 3 u ,
v ˙ + v r + κ y ˙ = 0 .
Assuming that T c is the period of the function v ( t ) , efficiency in terms of dimensionless parameters can be expressed as
η L = 1 T c 0 T c v 2 t κ   r d t ρ u 3 .
To characterize the efficiency η L of the system, it is essential to derive the voltage function v ( t ) generated by it, which can be accomplished using the harmonic balance method. Consequently, it was assumed that the limit cycle of the system would be described by a set of solutions in the following form:
y = A y c o s ω t ,
v = A v cos ω t + φ ,
where the four unknown quantities A y , A v , ω = ω ^ ω ^ n , and φ represent the dimensionless vibration amplitude, dimensionless voltage amplitude, dimensionless vibration frequency, and the phase shift between the oscillator vibrations and voltage oscillations, respectively. The parameter ω ^ denotes the unknown dimensional frequency of the system’s vibrations. Substituting solutions of the form Equation (6) into the model Equation (1) leads to the transformation of the system of differential equations into a system of algebraic equations:
A y 1 ω 2 c o s ω t + A y 3 ρ ω 3 a 3 u s i n 3 ω t + A y ω u ρ a 1 c s i n ω t A v   c o s ω t + φ = 0 ,
A v cos ω t + φ r ω sin ω t + φ A y κ   ω s i n ω t = 0 .
This condition will be satisfied for every moment of time t if and only if the sum of the coefficients with corresponding time functions equals zero. From Equation (7b), the following was deduced:
A v c o s φ r ω   s i n φ r = 0 ,
A v   r ω c o s φ + s i n φ + A y κ r ω r = 0 .
Equation (8a) shows that
t g φ = 1 ω r ,    sin φ = 1 ω r 2 + 1 ,    cos φ = ω r ω r 2 + 1 .
Based on Equation (8b) and Equation (9), the expression describing the relationship between the vibration amplitude and the voltage amplitude can be derived:
A v = A y r κ ω 1 + r 2 ω 2 .
By balancing the harmonics of Equation (7a), the remaining two algebraic equations necessary to determine an approximate solution of the device model are obtained. After the previously derived relations Equations (9) and (10) are taken into account, they take the form
A y   k e L   ω A y ω 2 1 = 0 ,
A y u ρ ω a 1 + A y 3 3 ρ ω 3 a 3 4 u A y c ω A y e L ω = 0 ,
where k e L = κ r 2 ω 2 1 + r 2 ω 2 and e L = κ r 1 + r 2 ω 2 are the piezoelectric stiffness and electrical damping of the linear system. Based on the above system of equations, it can be shown that
ω 1 2 ω 2 2 = r 2 1 + κ 1 ± 4 r 2 + r 2 1 + κ 1 2 2 r 2
A y 2 = 4 u c + e L u ρ a 1 3 ρ ω 2 a 3 .
Equation (12a) shows that ω 2 2 < 0 , regardless of the system parameters. Therefore, in the following part of the work, the notation ω 1 = ω has been adopted. Returning now to the general definition of efficiency Equation (5) and substituting the expressions Equations (6b), (10), and (12) into it, we obtain
η L = 2 e L c + e L u ρ a 1 3 u 2 ρ 2 a 3 .
Figure 2 illustrates the efficiency characteristics of the system. The efficiency is depicted as the ratio η L / η p L , and this representation is maintained throughout the study. The figure also includes the analogous characteristic obtained numerically for initial conditions y ( 0 ) = 0.1 , y ˙ ( 0 ) = 0 , n ( 0 ) = 0 , utilizing the fourth-order Runge–Kutta method.
In this figure, it can be observed that the system undergoes Hopf bifurcation at a certain speed, from now on referred to as the critical flow speed u c r . Its value results directly from Equation (13):
η L = 2 e L c + e L u ρ a 1 3 u 2 ρ 2 a 3 = 0 u = u c r = c + e L ρ a 1 ,
Another noteworthy observation from the same set of graphs is that irrespective of the parameter set, there exists a specific nominal flow speed u = u p L for which the efficiency η L attains a maximum value, denoted as the peak efficiency η p L . These quantities are given by the expressions
δ η L δ u = 0 u = u p L = 2 c + e L ρ a 1 = 2 u c r ,
η p L = η L u p L = a 1 2 e L 6 a 3 c + e L .
The aforementioned identities were previously derived in [3,4]. However, the remainder of the article discusses entirely original content. An unexplored property of GEHs, which can be inferred from Figure 2 or deduced from the identity u p L = 2 u c r , is noteworthy. A system with a low critical speed u c r will experience a more pronounced decline in efficiency due to the deviation of the flow speed u at which it operates from the nominal speed u p L . Let the measure of this phenomenon be the flow speed bandwidth in which the system efficiency η L does not fall below 90% of the maximum efficiency η p L , hereafter referred to as the high-efficiency band B L . According to the definition, B L is given as
η L = 0.9 η p L u 2 u 1 = 20 ± 4 10 9 c + e L ρ a 1 ,
B L = u 2 u 1 = 4 10 9 c + e L ρ a 1 1.4 u c r 0.7 u p L .
The quantities critical speed u c r , nominal speed u p L , peak efficiency η p L , and high-efficiency bandwidth B L (Figure 2) will be further referred to as the criterion parameters.

3. Efficiency of Two-Degree-of-Freedom System

According to the information provided in Section 1, the presence of additional mass in systems with two degrees of freedom implies the necessity to consider the placement of the electromechanical transducer. It may be positioned between the main mass and the stationary base (Figure 3a), between the masses (Figure 3b), or in both locations (Figure 3c).
To evaluate the influence of piezoelectric placement on efficiency, the efficiency characteristics of the variant with both transducers will be derived. Subsequently, the consequence of the absence of one of them will be examined.
The dynamics of devices within the discussed category can be described by a general dimensionless mathematical model of the form Equation (20). In addition to the identities in Equation (3), the following also holds:
y i = z i h ^ ,   v i = ξ θ i ^ m i ^   h ^   ω ^ n 2 ,   c i = c ^ i m ^ 1 ω ^ n ,   k 2 = k ^ 2 k ^ 1 ,   κ = θ i ^ 2 C p   i ^   m 1 ^    ω ^ n 2 , θ = θ ^ 1 θ ^ 2 ,   M = m ^ 2 m ^ 1 ,   i = 1,2
y ¨ 1 + y k 2 y 2 y 1 + c 1 y ˙ 1 c 2 y ˙ 2 y ˙ 1 v 1 + v 2 = ρ a 1 u y ˙ 1 + a 3 y ˙ 1 3 u ,
κ 1 v ˙ 1 + v 1 r 1 + y ˙ 1 ,
M y ¨ 2 + k 2 y 2 y 1 + c 2 y ˙ 2 y ˙ 1 v 2 = 0 ,
κ 2 v ˙ 2 + v 2 r 2 + y ˙ 2 y ˙ 1 = 0 .
To derive the solutions of the GEH2D mathematical model, the procedure outlined in [28] was adapted. It was assumed that the approximate solution of the model Equation (20) will have the following form:
y 1 = A 1 c o s ω 1 t + B 1 s i n ω 1 t + G 1 c o s ω 2 t + H 1 s i n ω 2 t ,
v 1 = n 1 c o s ω 1 t + n 2 s i n ω 1 t + n 3 c o s ω 2 t + n 4 s i n ω 2 t ,
y 2 = A 2 c o s ω 1 t + B 2 s i n ω 1 t + G 2 c o s ω 2 t + H 2 s i n ω 2 t ,
v 2 = ϑ 1 c o s ω 1 t + ϑ 2 s i n ω 1 t + ϑ 3 c o s ω 2 t + ϑ 4 s i n ω 2 t ,
Substituting the above identities into Equations (20b) and (20d) and then balancing the harmonics c o s ω 1 t , s i n ω 1 t , c o s ω 2 t , s i n ω 2 t allows for the derivation of the relationships between the amplitudes of voltage and the amplitudes of vibration.
n 1 = B 1 ω 1 ε D 11 A 1 δ D 11 ,       n 2 = A 1 ω 1 ε D 11 B 1 δ D 11 , n 3 = H 1 ω 2 ε D 12 G 1 δ D 12 ,       n 4 = G 1 ω 2 ε D 12 H 1 δ D 12 ,
ϑ 1 = θ ε D 21 θ B 2 B 1 ω 1 + δ D 21 A 1 θ A 2 , ϑ 2 = θ ε D 21 A 1 θ A 2 ω 1 + δ D 21 B 1 θ B 2 , ϑ 3 = θ ε D 22 θ H 2 H 1 ω 2 + δ D 22 G 1 θ G 2 , ϑ 4 = θ ε D 22 G 1 θ G 2 ω 2 + δ D 22 H 1 θ H 2 .
where δ D i j = κ i r 2 ω j 2 1 + r 2 ω j 2 and ε D i j = κ i r 1 + r 2 ω j 2 represent the piezoelectric stiffness and electric damping of the i -th piezoelectric for vibration with the j -th frequency.
Harmonic balancing of c o s ω 1 t and s i n ω 1 t of the algebraic equation obtained by substituting the solutions of Equation (21) into Equation (20c) and using the identities Equation (22) leads to the following system of equations:
A 2 = B 1 p 1 + A 1 q 1 ,
B 2 = B 1 q 1 A 1 p 1 ,
where
p 1 = ( k 2 + θ δ D 21 ) ( c 2 + θ 2 ε D 21 ) ω 1 ( c 2 + θ ε D 21 ) ( k 2 + θ 2 δ D 21 M ω 1 2 ) ω 1 ( c 2 + θ 2 ε D 21 ) 2 ω 1 2 + ( k 2 + θ 2 δ D 21 M ω 1 2 ) 2 ,
q 1 = ( c 2 + θ ε D 21 ) ( c 2 + θ 2 ε D 21 ) ω 1 2 + ( k 2 + θ δ D 21 ) ( k 2 + θ 2 δ T 21 M ω 1 2 ) ( c 2 + θ 2 ε D 21 ) 2 ω 1 2 + ( k 2 + θ 2 δ D 21 M ω 1 2 ) 2 .
By substituting the above identities into the harmonic balance equations for c o s ω 1 t and s i n ω 1 t of expression (20a), the following identities can be formulated:
A 1 α 1 + B 1 β 1 A 1 3 ρ a 3 ω 1 ( ( A 2 2 + B 2 2 ) ω 1 2 + 2 ( G 2 2 + H 2 2 ) ω 2 2 ) 4 u = 0 ,
A 1 β 1 B 1 α 1 + B 1 3 ρ a 3 ω 1 A 2 2 + B 2 2 ω 1 2 + 2 G 2 2 + H 2 2 ω 2 2 4 u = 0 ,
where
α 1 = ( p 1 ( k 2 + θ 2 δ D 21 ) + ( c 1 + ε T 11 + c 2 + ε D 21 ( 1 θ q 1 ) ) ω 1 u ρ ω 1 a 1 ) ,
β 1 = 1 + δ T 11 + δ D 21 θ + k 2 1 θ q 1 ω 1 p 1 c 2 + θ 2 ε D 21 + ω 1 .
Adapting the above procedure to balance the harmonics c o s ω 2 t and s i n ω 2 t of Equations (20c) and (20a) allows for the derivation of the following identities:
G 2 = H 1 p 2 + G 1 q 2 ,
H 2 = H 1 q 2 G 1 p 2 ,
G 1 α 2 + H 1 β 2 G 1 3 ρ a 3 ω 2 ( 2 ( A 1 2 + B 1 2 ) ω 1 2 + ( G 1 2 + H 1 2 ) ω 2 2 ) 4 u = 0
G 1 β 2 H 1 α 2 + H 1 3 ρ a 3 ω 2 ( 2 ( A 1 2 + B 1 2 ) ω 2 2 + ( G 1 2 + H 1 2 ) ω 2 2 ) 4 u = 0
where
p 2 = ( k 2 + θ δ D 22 ) ( c 2 + θ 2 ε D 22 ) ω 2 ( c 2 + θ ε D 22 ) ( k 2 + θ 2 δ D 22 M ω 2 2 ) ω 2 ( c 2 + θ 2 ε D 22 ) 2 ω 2 2 + ( k 2 + θ 2 δ D 22 M ω 2 2 ) 2 ,
q 2 = ( c 2 + θ ε T 22 ) ( c 2 + θ 2 ε D 22 ) ω 2 2 + ( k 2 + θ δ D 22 ) ( k 2 + θ 2 δ D 22 M ω 2 2 ) ( c 2 + θ 2 ε D 22 ) 2 ω 2 2 + ( k 2 + θ 2 δ D 22 M ω 2 2 ) 2 ,
α 2 = ( p 2 ( k 2 + θ 2 δ D 22 ) + ( c 1 + ε T 12 + c 2 + ε T 22 ( 1 θ q 2 ) ) ω 2 u ρ ω 2 a 1 ) ,
β 2 = 1 + δ D 12 + δ D 22 θ + k 2 1 θ q 2 ω 2 p 2 c 2 + θ 2 ε D 22 + ω 2 .
By adding Equation (24a) multiplied by B 1 to Equation (24b) multiplied by A 1 , an expression was obtained that allows the frequency ω 1 to be explicitly formulated:
1 + k 2 + δ D 11 k 2 q 1 + δ D 21 θ 1 θ q 1 p 1 c 2 + θ 2 ε D 21 ω 1 ω 1 2 = 0 .
Similarly, by adding Equation (26c) multiplied by H 1 to Equation (26d) multiplied by G 1 , an equation was obtained from which the frequency ω 2 can be derived:
1 + k 2 + δ D 12 k 2 q 2 + δ D 22 θ 1 θ q 2 p 2 c 2 + θ 2 ε D 22 ω 2 ω 2 2 = 0 .
Note that in each pair of parameters ( p 1 , p 2 ), ( q 1 , q 2 ), ( ε D 21 , ε D 22 ), ( δ D 11 ,   δ D 12 ), ( δ D 21 , δ D 22 ), the only difference is the frequency in their definition— ω 1 or ω 2 . Considering the similarity between Equations (28) and (29), one can conclude that the frequencies ω 1 and ω 2 must be equal, thereby excluding the possibility of polymodal vibrations in the system. The expression describing the vibration frequency ω 2 = ω 1 can therefore be derived by solving only one of the above equations.
The relationship between the vibration amplitude and the system frequency remains unknown. The first of the equations necessary to determine this relationship was obtained by subtracting Equation (24b) multiplied by B 1 from Equation (24a) multiplied by A 1 . The second one is the difference of Equation (26c) multiplied by G 1 and Equation (26d) multiplied by H 1 :
p 1 k 2 + θ 2 δ D 21 + c 1 + c 2 1 q 1 + ε D 11 + ε D 21 θ 1 θ q 1 ω 1 u ρ ω 1 a 1 3 ρ ω 1 A y 2 ω 1 2 + 2 G y 2 ω 2 2 a 3 4 u = 0 ,
p 2 k 2 + θ 2 δ D 22 + c 1 + c 2 1 q 1 + ε D 12 + ε D 22 θ 1 θ q 2 ω 2 u ρ ω 2 a 1 3 ρ ω 2 2 A y 2 ω 1 2 + G y 2 ω 2 2 a 3 4 u = 0 ,
where A y 2 = A 1 2 + B 1 2 and G y 2 = G 1 2 + H 1 2 . The set of Equation (30) has three non-trivial solutions in terms of A y 2 and G y 2 , which, after recalling the ω 2 = ω 1 identity, take the following form:
A y 2 = 4 u p 1 k 2 + θ 2 δ D 21 + c 1 + c 2 1 q 1 + ε D 11 + ε D 21 θ 1 θ q 1 u ρ a 1 3 ρ ω 1 3 a 3 , G y 2 = 0
A y 2 = 0 , A y 2 = 0 , G y 2 = 4 u p 1 k 2 + θ 2 δ D 21 ω 1 + c 1 + c 2 1 q 1 + ε D 11 + ε D 21 θ 1 θ q 1 u ρ a 1 3 ρ ω 1 3 a 3 ,
A y 2 = 4 u 3 p 1 k 2 + θ 2 δ D 21 ω 1 + c 1 + c 2 1 q 1 + ε D 11 + ε D 21 θ 1 θ q 1 u ρ a 1 9 ρ ω 1 2 a 3 , G y 2 = 4 u 3 p 1 k 2 + θ 2 δ D 21 ω 1 + c 1 + c 2 1 q 1 + ε T 11 + ε D 21 θ 1 θ q 1 u ρ a 1 9 ρ ω 1 2 a 3 .
These expressions, along with the previously derived identities, enable the explicit formulation of solutions for the system Equation (20) in the form Equation (21). It should be noted that the identities given by Equations (31a) and (31b) correspond to the same solution in the form of Equation (21); thus, only one of them, namely Equation (31a), will be further analyzed. Moreover, no set of system parameters and initial conditions has been found that would lead to the excitation of vibrations with amplitude Equation (31c). Therefore, this solution was considered unstable, and the efficiency has been derived based only on Equation (31a):
η D = 1 T c 0 T c v 1 2 t r 1 κ 1 + v 2 2 t r 2 κ 2 d t ρ u 3 = 2 λ ( μ u ρ a 1 ) 3 u 2 ρ 2 a 3 ,
where
λ = ( ε D 11 + ε D 12 θ 2 ( 1 2 θ q 1 + θ 2 ( p 1 2 + q 1 2 ) ) ) ,
μ = c 1 + c 2 + θ ε D 11 1 θ q 1 + p 1 k 2 + θ 2 δ T D 21 ω 1 .
The critical speed u c r D and nominal speed u p D of the GEH2D are given by the following expressions:
η D = γ 1 2 λ ( μ u ρ a 1 ) 3 u 2 ρ 2 a 3 = 0 u = u c r D = μ ρ a 1 ,
δ η D δ u = 0 u = u p D = 2 μ ρ a 1 = 2 u c r D .
Now it is possible to derive the peak efficiency of the system:
η p D = η D u p D = λ   a 1 2 6 γ 2 a 3 .
To fully determine the features of the analyzed variant, it is necessary to define its high-efficiency bandwidth B D . It is
η D = 0.9 η p D u 2 u 1 = 2 0.9 ρ a 1 μ ± 0.4 0.9 ρ a 1 μ ρ a 1 2 ,
B D = u 2 u 1 = 160 81 μ ρ a 1 1.4 u c r D 0.7 u p D .
Figure 4 depicts the relationship between efficiency η D and flow speed u , represented by the function Equation (32), compared with an analogous relationship obtained through numerical integration of the model Equation (20) for the following initial conditions: y 1 ( 0 ) = 0.1 , y ˙ 1 ( 0 ) = 0 , n 1 ( 0 ) = 0 , y 2 ( 0 ) = 0 , y ˙ 2 ( 0 ) = 0 , n 2 ( 0 ) = 0 . Moreover, the figure shows the values characterizing the efficiency of the variant— u c r D , u p D , η p D , and B I .

4. Comparison and Discussion

As suggested by Equations (34), (35), and (38), critical speeds, nominal speed, and high-efficiency bandwidths are interrelated in the same manner in both analyzed GEH variants. Consequently, comparing the efficiency of these systems can only be achieved by comparing their peak efficiencies. Since the efficiency of both devices depends on their nominal speeds, it is necessary to compare systems with the same nominal speeds. The electrical damping of the linear system e L will therefore be selected in such a way that this condition is satisfied. In the following derivation, it is assumed that the structural damping of the linear system c is equal to the structural damping c 1 of the damper connecting the lower mass of the GEH2D with the base, i.e., c = c 1 . The shape of the flowing body, represented by the coefficients a 1 and a 3 , and the density of the fluid remain the same for all variants.
u p I = u p L ,
2 μ + γ 2 ρ a 1 = 2 e L + c 1 ρ a 1 ,
e L = γ 2 + μ c 1 .
According to Equations (16) and (39c), the efficiency of a system with one degree of freedom with a nominal speed equal to the nominal speed of the GEH2D is equal to
η L u p D = μ c 1 a 1 2 6 μ a 3 .
The ratio of peak efficiencies of compared systems with the same nominal speed is therefore given by
η L u p D η D u p D = s 1 + s 2 s 1 + s 3 ,
where
s 1 = ε T 11 ( ( k 2 + θ 2 δ D 21 ) ( k 2 + θ 2 δ T 21 2 M ω 1 2 ) + ( c 2 + θ 2 ε D 21 ) 2 ω 1 2 + M 2 ω 1 4 ) ,
s 2 = M ω 1 2 ( k 2 ε D 21 c 2 δ D 21 ) ( θ 1 ) θ + M ω 1 2 ( c 2 + θ ε D 21 ) ,
s 3 = θ 2 ε D 21 M ω 1 2 ( 2 k 2 ( θ 1 ) + M ω 1 2 ) + ( c 2 2 ω 1 2 + k 2 2 ) ( θ 1 ) 2 .
Clearly, according to Equation (41), the discussed system with two degrees of freedom will have greater efficiency than the reference system if and only if s 3 > s 2 . For variant with only one lower transducer (Figure 3a), where ε D 21 = δ D 21 = 0 , this inequality takes the form
0 c 2 M 2 ω 1 4 ,
which implies that such a device cannot be more efficient than the reference system, regardless of its parameters. This conclusion contradicts the results presented in [25]. The reason for this discrepancy lies in the fact that in the cited work, the operating conditions of the compared systems were not standardized—the devices had different critical speeds and, consequently, different nominal speeds. Despite the unquestionable value of this article, the conclusion stated therein can be subject to questioning. In the case of the variant with only the upper transducer (Figure 3b), where ε T 11 = δ T 11 = 0 , the condition for efficiency improvement takes the following form:
ε D 21 ε D 21 + c 2 ,
which indicates that this variant does not offer efficiency enhancement. Similarly, an identical and impossible-to-satisfy requirement is associated with the special case of the third variant, where there are two identical transducers, meaning ε D 11 = ε T 21 , δ D 11 = δ D 21 , and θ = 1 . However, the inequality s 3 > s 2 can be satisfied for the most general variant—the one equipped with two different transducers. In this case, it can be reduced to the condition
θ 1 θ ε D 21 k 2 2 θ 1 θ + c 2 δ D 21 M ω 1 2 + c 2 2 ε D 21 ω 1 2 θ 1 θ + ε D 21 k 2 M ω 1 2 2 θ 1 > M 2 ω 1 4 ( c 2 ε D 21 ( θ 1 ) θ ) ,
which, after substituting the solution of Equation (27), can be solved numerically for a chosen parameter. For the parameters presented in Table 2, solving the above inequality for θ results in obtaining a threshold value of approximately θ 1.13 . An example of the efficiency characteristics of the GEH with e L given by Equation (39c) and GEH2D is shown in Figure 5, where condition Equation (45) has been met by adopting parameters according to Table 2, with θ = 1.3 > 1.13 . It should be noted that parametric analysis does not allow one to indicate the optimal gain that the discussed system can provide. Determining this fact can be considered as a potential direction for further research.

5. Conclusions

The aim of the study was to investigate the efficiency of a galloping energy harvester with two degrees of freedom, considering its three variants, which differ in the number and location of electromechanical transducers. The realization of this objective commenced with the analysis of the reference variant with one degree of freedom. Utilizing the harmonic balance method, an approximate solution of the mathematical model of the system was derived, followed by the formulation of an expression describing the efficiency of the variant. Based on this, key criteria parameters were defined, providing comprehensive information about the variant’s efficiency: peak efficiency, high-efficiency bandwidth, critical speed, and nominal speed.
Subsequent sections of the study delineated the multitude of configurations that a system with two degrees of freedom can adopt. Three different subvariants of the device were characterized, differing from each other by the location and number of the electromechanical transducers. Parameters characterizing the efficiency of all subvariants were derived, demonstrating that the critical speed, nominal speed, and high-efficiency bandwidth are related to each other in the same manner as in the case of the reference system. It was then demonstrated that among the three indicated subvariants of the two-degree-of-freedom system, only the one with two transducers can be more efficient than the reference system.

Author Contributions

Conceptualization, F.S.; methodology, F.S.; software, F.S.; validation, R.S.; formal analysis, F.S. and R.S.; investigation, F.S.; resources, F.S.; data curation, F.S.; writing—original draft preparation, F.S.; writing—review and editing, R.S.; visualization, F.S.; supervision, R.S.; project administration, F.S.; funding acquisition, R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Polish Ministry of Higher Education, grant number 0612/SBAD/3628.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumar, S.; Tiwari, P.; Zymbler, M. Internet of Things is a revolutionary approach for future technology enhancement: A review. J. Big Data 2019, 6, 111. [Google Scholar] [CrossRef]
  2. Bladeless Vortex Official Site. Available online: https://vortexbladeless.com (accessed on 1 May 2024).
  3. Den Hartog, J.P. Mechanical Vibrations, 4th ed.; McGraw-Hill: New York, NY, USA, 1956. [Google Scholar]
  4. Barrero-Gil, A.; Alonso, G.; Sanz-Andrés, A. Energy harvesting from transverse galloping. J. Sound Vib. 2010, 329, 2873–2883. [Google Scholar] [CrossRef]
  5. Ludlam, D.V. Optimal Energy Harvesting Form Vortex Induced Vibration and Transverse Galloping Vibrations. Ph.D. Thesis, Universidad Politécnica de Madrid, Madrid, Spain, 2017. [Google Scholar]
  6. Ibarra, D.; Sorribes, F.; Alonso, G.; Meseguer, J. Transverse galloping of two-dimensional bodies having a rhombic cross-section. J. Sound Vib. 2014, 333, 2855–2865. [Google Scholar] [CrossRef]
  7. Wang, J.; Tang, L.; Zhao, L.; Zhang, Z. Efficiency investigation on energy harvesting from airflows in HVAC system based on galloping of isosceles triangle sectioned bluff bodies. Energy 2019, 172, 1066–1078. [Google Scholar] [CrossRef]
  8. Noel, J.; Yadav, R.; Li, G.; Daqaq, M. Improving the performance of galloping micro-power generators by passively manipulating the trailing edge. Appl. Phys. Lett. 2018, 12, 083503. [Google Scholar] [CrossRef]
  9. Zhao, D.; Hu, X.; Tan, T.; Yan, Z.; Zhang, W. Piezoelectric galloping energy harvesting enhanced by topological equivalent aerodynamic design. Energy Convers. Manag. 2020, 222, 113260. [Google Scholar] [CrossRef]
  10. Zhao, K.; Zhang, Q.; Wang, W. Optimization of Galloping Piezoelectric Energy Harvester with V-Shaped Groove in Low Wind Speed. Energies 2019, 12, 4619. [Google Scholar] [CrossRef]
  11. Zheng, J.; Li, Z.; Zhang, H. Low-Wind-Speed Galloping Wind Energy Harvester Based on a W-Shaped Bluff Body. Energies 2024, 17, 958. [Google Scholar] [CrossRef]
  12. Wang, J.; Ren, H.; Qin, W. High-performance piezoelectric wind energy harvester with Y-shaped attachments. J. Appl. Phys. 2019, 181, 645–652. [Google Scholar] [CrossRef]
  13. Yuan, Y.; Wang, H.; Yang, C.; Sun, H.; Tang, Y.; Zhang, Z. Exploring the Potential of Flow-Induced Vibration Energy Harvesting Using a Corrugated Hyperstructure Bluff Body. Micromachines 2023, 14, 1125. [Google Scholar] [CrossRef]
  14. Ewere, F.; Wang, G.; Frendi, A. Galloping Piezoelectric Energy Harvester with Bio-inspired Square Bluff Body. In Proceedings of the 23rd AIAA/AHS Adaptive Structures Conference Materials, Kissimmee, FL, USA, 5–9 January 2015. [Google Scholar]
  15. Wang, J.; Sun, S.; Hu, G.; Yang, Y.; Tang, L.; Li, P.; Zhang, G. Exploring the potential benefits of using metasurface for galloping energy harvesting. Energy Convers. Manag. 2021, 243, 114414. [Google Scholar] [CrossRef]
  16. Hai, W.; Sun, H.; Zhang, Z.; Tang, Y.; Yang, C.; Li, H. Improving the performance of galloping energy harvester with striped bluff body. AIP Adv. 2023, 13, 075019. [Google Scholar]
  17. Abdelkefi, A.; Yan, Z.; Hajj, M.R. Nonlinear dynamics of galloping-based piezoaeroelastic energy harvesters. Eur. Phys. J. Spec. Top. 2013, 222, 1483–1501. [Google Scholar] [CrossRef]
  18. Abdelkefi, A.; Yan, Z.; Hajj, M.R. Modeling and nonlinear analysis of piezoelectric energy harvesting from transverse galloping. Smart Mater. Struct. 2013, 22, 025016. [Google Scholar] [CrossRef]
  19. Li, H.T.; Ren, H.; Cao, F.; Qin, W.Y. Improving the galloping energy harvesting performance with magnetic coupling. Int. J. Mech. Sci. 2023, 237, 107785. [Google Scholar] [CrossRef]
  20. Zhang, H.; Zhang, L.; Wang, Y.; Yang, X.; Song, R.; Sui, W. Modeling and experimental investigation of asymmetric distance with magnetic coupling based on galloping piezoelectric energy harvester. Smart Mater. Struct. 2022, 31, 065007. [Google Scholar] [CrossRef]
  21. Sarbinowski, F.; Starosta, R. Assessment of the role of structural nonlinearity in galloping energy harvesters. Vib. Phys. Syst. 2021, 32, 2021209. [Google Scholar]
  22. Alhadidi, A.; Abderrahmane, H.; Daqaq, M. Exploiting stiffness nonlinearities to improve flow energy capture from the wake of a bluff body. Phys. D Nonlinear Phenom. 2016, 337, 30–42. [Google Scholar] [CrossRef]
  23. Yang, K.; Yang, K.; Wang, J.; Wang, J.; Yurchenko, D. A double-beam piezo-magneto-elastic wind energy harvester for improving the galloping-based energy harvesting. Appl. Phys. Lett. 2019, 115, 193901. [Google Scholar] [CrossRef]
  24. Wang, J.; Geng, L.; Zhou, S.; Zhang, Z.; Lai, Z. Design, modeling and experiments of broadband tristable galloping piezoelectric energy harvester. Acta Mech. Sin. 2020, 36, 592–605. [Google Scholar] [CrossRef]
  25. Lan, C.; Tang, L.; Hu, G.; Qin, W. Dynamics and performance of a two degree-of-freedom galloping-based piezoelectric energy harvester. Smart Mater. Struct. 2019, 28, 045018. [Google Scholar] [CrossRef]
  26. Hu, G.; Wang, J.; Qiao, H.; Zhao, L.; Li, Z.; Tang, L. An Experimental Study of a Two-Degree-of-Freedom Galloping Energy Harvester. Int. J. Energy Res. 2020, 45, 3365–3374. [Google Scholar] [CrossRef]
  27. Hu, G.; Wang, J.; Tang, L. A Comb-Like Beam based Piezoelectric System for Galloping Energy Harvesting. Mech. Syst. Signal Process. 2021, 150, 107301. [Google Scholar] [CrossRef]
  28. Hu, G.; Liang, J.; Tang, L.; Wang, J. Improved Theoretical Analysis and Design Guidelines of a Two-Degree-of-Freedom Galloping Piezoelectric Energy Harvester. J. Intell. Mater. Syst. Struct. 2021, 33, 210–230. [Google Scholar] [CrossRef]
  29. Zhao, L.; Tang, L.; Yang, Y. Enhanced piezoelectric galloping energy harvesting using 2 degree-of-freedom cut-out cantilever with magnetic interaction. Jpn. J. Appl. Phys. 2014, 53, 060302. [Google Scholar] [CrossRef]
Figure 1. Rodel of reference variant.
Figure 1. Rodel of reference variant.
Applsci 14 05427 g001
Figure 2. Efficiency characteristics of GEH for κ = 1.3 , r = 3 , c = 0.1 , ρ = 0.02 , a 1 = 2.3 , a 3 = 18 .
Figure 2. Efficiency characteristics of GEH for κ = 1.3 , r = 3 , c = 0.1 , ρ = 0.02 , a 1 = 2.3 , a 3 = 18 .
Applsci 14 05427 g002
Figure 3. Subvariants of two-degree-of-freedom system: (a) transducer between base and lower mass, (b) transducer between masses, (c) two transducer.
Figure 3. Subvariants of two-degree-of-freedom system: (a) transducer between base and lower mass, (b) transducer between masses, (c) two transducer.
Applsci 14 05427 g003
Figure 4. Efficiency characteristic of GEH2D for parameter values presented in Table 2.
Figure 4. Efficiency characteristic of GEH2D for parameter values presented in Table 2.
Applsci 14 05427 g004
Figure 5. Comparison of GEH and GEH2D having equal nominal speeds.
Figure 5. Comparison of GEH and GEH2D having equal nominal speeds.
Applsci 14 05427 g005
Table 1. List of GEH parameters.
Table 1. List of GEH parameters.
Symbol P a r a m e t e r U n i t
m ^ Mass k g
k ^ Stiffness coefficient N / m
c ^ Damping coefficient k g / s
u ^ Flow speed m / s
ρ ^ Planar fluid density k g / m 2
h ^ Characteristic dimension m
θ ^ Piezoelectric coefficient N / V
C p ^ Equivalent capacity F
R ^ Circuit electrical resistance Ω
τ Time s
z = z τ Vibration vs. time function m
ξ = ξ τ Voltage vs. time function V
Table 2. GEH2D system parameter values.
Table 2. GEH2D system parameter values.
ParameterValue
M 0.5
k 2 0.3
c 1 0.1
c 2 0.25
κ 1 0.9
κ 2 1.8
r 1 3
r 2 2
θ 1.3
a 1 2.3
a 3 −18
ρ 0.02
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sarbinowski, F.; Starosta, R. On Efficiency of Two-Degree-of-Freedom Galloping Energy Harvesters with Two Transducers. Appl. Sci. 2024, 14, 5427. https://doi.org/10.3390/app14135427

AMA Style

Sarbinowski F, Starosta R. On Efficiency of Two-Degree-of-Freedom Galloping Energy Harvesters with Two Transducers. Applied Sciences. 2024; 14(13):5427. https://doi.org/10.3390/app14135427

Chicago/Turabian Style

Sarbinowski, Filip, and Roman Starosta. 2024. "On Efficiency of Two-Degree-of-Freedom Galloping Energy Harvesters with Two Transducers" Applied Sciences 14, no. 13: 5427. https://doi.org/10.3390/app14135427

APA Style

Sarbinowski, F., & Starosta, R. (2024). On Efficiency of Two-Degree-of-Freedom Galloping Energy Harvesters with Two Transducers. Applied Sciences, 14(13), 5427. https://doi.org/10.3390/app14135427

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop