Abstract
An omnidirectional acoustic energy harvester (AEH) based on six Helmholtz resonators is proposed in this work. Compared with the previous structure, the insufficiency of the directionality and conversion efficiency of energy collection can be effectively improved due to the coupling of six resonators. Based on the distributed parameter model, the relationship of the electrical output, the input frequency with the structure size is obtained. The simulation results show that the maximum output voltage is 70.95 mV at the resonant frequency of 35 kHz. When the external load resistance is 14 kΩ, the maximum output power is 0.45 μW. Moreover, the energy conversion efficiency of this omnidirectional AEH can reach 23%, which is improved greatly compared with the traditional structure. Therefore, this AEH will have a wide range of application prospects in medical implantation equipment and other fields.
Introduction
Energy harvester is a development trend that is expected to replace traditional batteries in the future. The energy conversion methods are divided into electrostatic [1], electromagnetic [2] and piezoelectric [3], etc., and external excitation includes vibration energy [4], solar energy [5], magnetic energy [6], and acoustic energy [7] etc. Due to the excellent performance of piezoelectric materials and the significant advantages of acoustic as the excitation source, miniature piezoelectric energy harvesters that use acoustic as the excitation source have gradually become a research hotspot [8–11]. In response to the problem of insufficient external excitation, researchers at home and abroad mainly amplify sound waves with the cavity-shaped structure, and then use the piezoelectric effect to achieve the conversion from sound energy to electrical energy. In 2006, Horowitz et al. proposed an AEH with Helmholtz cavity. After the magnification of the cavity structure, its output power is 6 pW under the sound pressure level of 149 dB [12]. In 2008, F. Liu et al. presented the development of an AEH using an electromechanical Helmholtz resonator, and experimental results indicated that approximately 30 mW of output power was harvested for an incident sound pressure level of 160 dB with a flyback converter [13]. Their research provides a reliable basic reference for subsequent researchers. However, the traditional AEH based on Helmholtz resonant cavity has some performance defects, including: single energy harvesting direction, narrow resonance frequency band, and low output power. In order to solve the above problems, scholars mainly use the following three methods to improve output performance.
By designing multiple resonant cavities to improve output performance. In 2013, X. Peng et al. designed an acoustic energy harvester using dual Helmholtz resonators. Compared with the traditional Helmholtz resonator, the voltage generated by this structure is increased by 400%, and the output power was 16 times of the original [14]. By designing a novel cavity structure to improve output performance. In 2016, F. Khan et al. proposed an AEH with a tapered Helmholtz resonator, and they improved the output performance by changing the cavity structure. It could receive an output power of 90.61 μW with an excitation of 130 dB [15]. By designing a novel coupling structure. In 2020, P. Eghbalia et al. designed an AEH with an axenic latticed resonator backed by an acoustic rectangular tube. It could arrive at a large magnification factor of around 10.5 for a 100 dB sound pressure level at resonance [16].
However, none of these structures could solve the problems of insufficient external excitation, single energy collection direction, and narrow resonance frequency band at the same time. In our previous study, a quarter spherical acoustic energy harvester was designed. It could harvest multi-directional acoustic energy, had a widened resonance frequency band, and had an improved energy conversion efficiency. Its output voltage was about 28 mV, and its output power was about 0.05 μW [17].
When it comes to practical application, the higher standards is needed for some key performance indicators such as the output bandwidth, collection direction and energy conversion efficiency. On the basis of the quarter-spherical structure, an omnidirectional AEH based on six Helmholtz resonators is proposed in this work, which can collect the sound energy in the environment from all directions. As a result, some key performance indicators are significantly improved. In Section 2, the theoretical model of this omnidirectional AEH is discussed, and the relationship of the output voltage with the input sound pressure is obtained. In Section 3, the output performance of this omnidirectional AEH, including sound pressure level, output voltage, and output power, are studied. Finally, some conclusions are drawn in Section 4.
Principle and theory
As shown in the Fig. 1, the omnidirectional acoustic energy harvester contains six independent resonators and an external piezoelectric layer, which can collect the sound energy in the environment from all directions. Taking the rectangular coordinate system as a reference, a single set of energy harvesters A, B, and C can collect energy in the opposite directions of x, y, and z, and the energy harvesters D, E, and F can be used to collect energy in the positive direction of x, y, and z. All energy harvesters have the same size.
When excited by incident acoustic wave, the air in the neck moves downward from the static equilibrium position into the cavity. After that, the air in the cavity is compressed. Then, a downward pressure is formed and acts on the piezoelectric layer. Similarly, when the air in neck area moves down to the extreme position, the internal pressure is greater than external incentives, and it will push the air back to the upward direction. Due to inertia, the air in neck area returns to the equilibrium position and continues to move upward, and the pressure in cavity is lowered. As a result, the film produces a reverse strain. In this periodic, an equal amount of positive and negative bound charges appears on the upper and lower surface of the piezoelectric film. Thus, when the external charge has accumulated enough, an effective output voltage is produced.
The shell theoretical model is more suitable for this structure. However, it is difficult to obtain an analytical solution for the shell theoretical model. Therefore, a simplified beam theoretical model is used to study the one-way coupling mechanism of the omnidirectional acoustic energy harvester based on six Helmholtz resonators in this work. Only the output performance of the electrical energy generated by vibration is studied, and the coupling mechanism is not studied. The structural parameters that affect the performance can be obtained according to the simplified beam theoretical model, and the finite element software Comsol is used to simulate its output performance.

(a) Overall viewing (Separate structure for easy observation). (b) Single group viewing of omnidirectional AEH.
The sound path equation in the cavity can be obtained by electric acoustic analogy [18,19]:
The highest sound pressure of the cavity is generated at the acoustic resonant frequency, which is determined by the cavity volume and neck size of the cavity. The acoustic resonant frequency f
a
and the maximum sound pressure amplification coefficient G
max are expressed as [20–22]:
The force F
0 on the piezoelectric layer can be obtained as
A distributed parameter model is used to analyze piezoelectric conversion, and its advantage lies on its ability to predict high-order modal vibration, considering the dynamic mode shape and the stress distribution on the cantilever beam during the vibration process. In the distributed parameter model, one end of the piezoelectric vibrator is clamped, and the other end is free. The upper and lower electrodes are connected to an external load. The external vertical excitation received by the piezoelectric vibrator is g(t), and the small rotation excitation at the clamping end is h(t). The damping of the system mainly comes from two aspects, which are aeroelastic damping 𝜉
a
and the strain rate damping 𝜉
s
. In the vibration process, both of them satisfy the proportional damping criterion. The Euler–Bernoulli equation can be used to obtain the dynamic equation of the piezoelectric vibrator [23]:
The model is solved by decomposing the vibration displacement. The absolute vibration displacement of the clamping end can be decomposed into a vertical vibration excitation and a rotational vibration excitation. With the mode superposition method, the relative vibration displacement can be expressed as the series of the mode function, which is composed of the standardized mode function and the generalized coordinates of the vibration mode [24,25].
Based on the boundary conditions, the solution of the mode can be calculated as:
With the orthogonal conditions and boundary conditions, the modal solution can be calculated as [26]:
When the tiny rotational motion is not considered, that is, h(t) = 0; and when in the first resonance mode, that is, r = 1, the formula ((16)) can be simplified as:
The resonator was made of copper, and the piezoelectric layer was made of PZT-5H. The highest output power can be obtained only when the mechanical resonance and the acoustic resonance are coupled. It is a simpler and more convenient method to change the acoustic resonant frequency parameters to make it close to the mechanical resonant frequency. According to the equation ((2)), the acoustic resonant frequency is related to the neck size and cavity volume. Therefore, by adjusting the size of the neck of the resonator, the resonant frequency of the sound pressure can be changed to achieve the coupling of the two resonance frequencies.
According to the results of the previous theoretical analysis, when the radius of the neck is 0.15 mm and the length of the neck is 0.115 mm, the mechanical resonant frequency is about 35.1 kHz, and the acoustic resonant frequency is about 35.109 kHz. The two are basically in resonance and are in the best working condition.
The material parameters and the structural dimensions of this omnidirectional AEH are shown in Tables 1 and 2, respectively.
Material parameters
Material parameters
Structural dimensions
In the simulation, Comsol 5.3a was used for finite element simulation to analyze the model, including four physical fields: solid mechanics, static electricity, pressure acoustics and electric circuit. The relationship of the output characteristics with the incident sound frequency of this system is studied. The sound field space uses a three-dimensional hexahedron with a side length of 4 mm, which is filled with ideal medium (air). The front, bottom and left sides of the cube are set as the incident pressure field. The back, top and right sides are set as hard sound field boundaries. The sound source is set as a far-field plane wave with a sound pressure of 110 pa, and the transmission direction is e = (−1, −1, 1). The energy harvester is located in the middle of the cube, and the boundary conditions of piezoelectric film satisfies mechanical clamping and electrical short-circuit. The positive electrode and the negative electrode are respectively provided on the upper and lower surface of the piezoelectric film, which are connected with an external resistor. The energy harvester periodically strains through the amplified incident sound wave. Due to the piezoelectric effect, the strain is converted into a periodic voltage output, which is collected by an external resistor. Finally, the external air grid is set as normal, and the internal energy harvester grid is set as fine, with a minimum unit of 0.02 mm to obtain more accurate results. On this basis, inspections of different mesh thicknesses are also carried out, and the results are basically unchanged.
When simulating, the input acoustic frequency is set as independent variable. According to the equation ((3)) and the structural size of the model, it can be calculated that the maximum multiple of the external sound pressure is about 51. According to the equation ((18)), the amplified sound pressure level is 34 dB [27].
When the external sound pressure of the system is 110 pa, the equivalent sound pressure level is 135 dB. Since the sound pressure is amplified by 51 times, the theoretically calculated sound pressure level inside the cavity is about 169 dB. As shown in Fig. 2, the maximum sound pressure level of the simulation result is 170 dB at 35 kHz, which is almost the same as the theoretical result. It reflects the sound pressure amplification of the resonant cavity structure. Therefore, according to the calculation and the simulation results, it can be found intuitively that the Helmholtz resonant cavity can amplify the input by 51 times compared with the cavityless structure, and the output voltage of the energy harvester will also differ by several tens of times. The advantage of the cavity structure meets the need of practical applications.

The simulation results of sound pressure level.
When the two resonant frequencies of the energy harvester are coupled, the system reaches the best resonance state. The theoretical calculation and the simulation results are compared. Since the vibration of the system is mainly based on the first-order mode, the theoretical result can be obtained according to the equation ((17)). The maximum output voltage of the energy harvester is 67.61 mV at 35.1 kHz. For simulation results, the incident sound frequency is set as the independent variable, which ranges from 30 kHz to 40 kHz, with a 200 Hz interval. The resonant frequency is 35 kHz, and the maximum output voltage result is 70.95 mV. As shown in Fig. 3, the frequency interval of the simulation and calculation is very small, and the two curves have similar changing trends, so the results of voltage are relatively accurate. Compared with the theoretical result, the curve on the right side of the simulation result drops faster. The reason is that the damping value in the simulation is larger than that in the equation ((17)), which causes the peak voltage to decrease faster. When it comes to the resonant frequency band, 0.707 times the peak voltage is also regarded as the effective working voltage. So that the effective bandwidth of the simulation result is about 700 Hz. The effective bandwidth of the theoretical result is about 950 Hz, which is also caused by the difference of the damping value. With the electromechanical acoustic analogy method, the acoustics, mechanics, and kinematics in this model are studied using Euler–Bernoulli equation. The aeroelastic damping and the strain rate damping are considered in the theoretical model, and the finite element analysis of this model is carried out to incorporate the aeroelastic damping and the strain rate damping. In addition, the isotropic loss factor of Cu is studied to show better prediction of the performance of the proposed method.

Comparison of theoretical and simulated results.

Structure of the quarter spherical AEH.
The quarter-spherical energy harvester has been shown in Fig. 4. Comparing the voltage results of this structure with that of a quarter-spherical energy harvester in the same size, since all the simulation conditions are the same, the improvement in output performance can be visually demonstrated. It can be seen that the voltage performance of this structure has been significantly improved, the amplitude is about 153.4% as shown in Fig. 5.

Voltage comparison results of omnidirectional and quarter spherical AEH.
According to the displacement results obtained by the simulation, the superiority of the omnidirectional structure can be demonstrated. As shown in Fig. 6, the vibration displacement of the energy harvester at the resonant frequency of 35 kHz reaches the maximum, which is about 23.1 nm. At 34 kHz and 36 kHz near the resonant frequency, the vibration displacement is relatively small, respectively 17.6 nm and 4.15 nm, which also reflects the role of the resonance peak. It can be clearly found that whether it is at or near the resonance frequency, the vibration displacement of the overall resonator is not much different among the six cavities, less than 10%. It is found that the vibration is more balanced. The reason is that the six Helmholtz resonators of this AEH correspond to all directions in the environment, and the utilization rate of acoustic energy is also higher. Because the external piezoelectric film is connected as a whole, six Helmholtz resonators can drive the piezoelectric film together, so it can also generate a higher output voltage.

Vibration displacement of omnidirectional energy harvester at (a) 34 kHz (b) 35 kHz (c) 36 kHz.
The output power of the energy harvester is related to the output voltage and the external resistance. According to formula (19), the expression of output power can be expressed as [28].

The relationship between output power and (a) external resistance and (b) sound frequency.
As shown in Fig. 7(a), at the resonant frequency, when the resistance of the external resistor is 14 kΩ, the maximum output power result is 0.45 μW, which also reflects that the internal equivalent impedance of the energy harvester is 14 kΩ. When putting it into practical application, it should also be used as a reference standard to design external circuits. At the same time, after determining the best external resistance, the acoustic frequency is used as an independent variable to verify the result of the output power. At this time, the external resistance is 14 kΩ. It can be found from Fig. 7(b) that the maximum output power at the resonant frequency of 35 kHz is also 0.45 μW, which is the same as the result in the Fig. 7(a).
The input power is produced by the incident sound wave, which can be expressed as:
The energy conversion efficiency is:
Comparison of acoustic energy harvesters based on resonator structures
An omnidirectional Helmholtz resonant cavity acoustic energy harvester is proposed in this work. By setting six Helmholtz resonant cavities, the acoustic energy in all directions can be collected, thus greatly improving the energy conversion and output performance. When simulating, the input acoustic frequency is set as independent variable, the sound pressure level, output voltage and output power are set as dependent variables. Results show that the sound pressure amplification factor of the energy harvester is about 51, the maximum output voltage at the resonant frequency of 35 kHz is 70.95 mV, and the maximum output power is 0.45 μW. Moreover, the effective bandwidth is about 700 Hz, and the energy conversion efficiency reaches to 23%. The electrical output is more efficient and meets the need of practical application.
Footnotes
Acknowledgements
This work was supported by the National Natural Science Foundation of China (61904085, 61904089), Province Natural Science Foundation of Jiangsu (BK20190731), the China Postdoctoral Science Foundation (2017M621692), and Jiangsu Postdoctoral Foundation (1701131B).
