1. Introduction
A VEH has proven worthy of having the capacity to sustainably supply electrical power to wireless sensor nodes (WSNs) and body sensor networks (bodyNETs) [
1] by scavenging ambient vibrational energy and converting it into useful electrical energy. The scope of scavenging environmental energy has recently attracted interest because it opens a pathway to realizing the sustainable development goals of cutting down carbon footprints that are mostly produced as waste product during energy generation. Previous work concluded structural and electrical based optimization of the energy harvester [
2]. An innovative approach on the two-degree-of-freedom (2DOF) linear vibration energy harvester for train application and optimization was analyzed on Multiphysics using the result from numerical magnetic field simulation [
3]. The model prototype was fabricated and found consistent with the obtained theoretical results. A study on dynamic responses in 2DOF systems, based on different electrical coil connection and geometry to ascertain which connection mode gives an optimum performance, concluded that series connection is the best in terms of harvested voltage/power, operational bandwidth and normalized power [
4]. A path to achieving power maximization on a size optimized/miniature A-battery sized VEH design was reported [
5]. The author coupled a non-magnetic (tungsten) inertial mass alongside the axially oriented oscillating magnet to compensate for resonance due to miniaturization. Refs. [
6,
7] reported a ring-shaped, coil–magnet transducer architecture of VEH with Halbach configuration where a linear Halbach array concentrates its magnetic field in the inner space of the mechanism where a vertically centered single coil was located to increase the resonant mass within a fixed dimension of the transducer. In a separate attempt, the level of the flux density in the power line was measured using an FEMM software during fabrication of the indoor power line based magnetic field energy harvester [
8]. In a similar endeavor, measuring the flux density using FEMM on a Halbach array was mentioned in [
7].
Without a loss of generality, this paper focuses on realizing an approach to ensure an accurate prediction of the optimum overall size that will maximize the coupling coefficient and power output on the electromagnetic transducer of a VEH.
2. Harvesters Governing Equation
A vibration energy harvester is a device that scavenges and transforms ambient vibration into useable electrical energy that can power sensor nodes. The VEH comprises a coil placed in the field of a permanent magnet such that, during vibration, the coil that is fixed to the free end of a fixed-free mechanical structure will freely oscillate. As much as engineers have keen interest in realizing the above objectives, cost and size optimization remain a valuable pearl held in high esteem during fabrication/design. This work presents a finite element simulation approach to realize size optimization based on the level of the magnetic flux density/coupling in the iron–magnet–coil part of an electromagnetic vibration energy harvester. The focus in this work will be to optimize the iron–magnet–coil geometry with the view to realize more compact, lightweight and cost-effective iron–magnet–coil designs. The general geometry employed to fully characterize the transduction iron–magnet–coil, which will be modeled in the FEMM software, is shown in
Figure 1.
From
Figure 1,
,
and
are the respective heights of coil, pure iron plate which prevent flux leakages, separation distance between coil and magnets and outer magnet while
and
are the respective widths/thicknesses of the materials earlier mentioned. Equation (1) gives an expression for computing the terms
and
.
The sufficient clearance between the coil and the magnet,
, was fixed to 1.5 mm to prevent contact during excitation. The total width of the iron–magnet–coil obtained from
Figure 1 is shown in Equation (2).
When the geometry is visualized on a 3D plane, the model protrudes by a fixed length
into the page. During excitation, the magnetic flux of the permanent magnet couples into the freely oscillating coil, hence, voltages are induced in the coil according to the principle of electromagnetic induction.
is defined as the degree of coupling that was obtained [
9].
where
and
are effective turn, flux density, effective length and coil fill factor. If none of the magnetic flux couples into the coil, an approximately zero voltage is induced in the coil. Therefore, the coil materials are carefully selected (copper in this case) to ensure that the highest possible degree of coupling is realized. The selected copper wire material is wound into N turns on a nonconductive circular brace to achieve a total coil width/thickness
.
The Maxwell theory reported divergence and the curl of the flux density where
,
and
are the current density, displacement field, magnetic field, electric field and magnetizing field.
μ is the permeability of the magnetic material.
The physical meaning of Equations (4) and (5) asserts that, for any magnetic system/magnet, there are no isolated magnetic poles, and circulating magnetic fields are produced by changing electric currents. In the eventuality of using more than one magnet, Equation (4) sets an order for which the transduction magnet must be aligned to allow for continuous flux linkage between the several magnets in such a manner that no pole is isolated. In the region of no charge, and are all equal to zero in Equations (4) and (5). Therefore, the equations predicted a non-changing flux value since no external source of electric charge is in the system.
3. Iron–Magnet–Coil Simulation and Coupling Equations
Before the flux density was simulated on FEMM, an initial approach was taken to characterize the flux on a
magnet using a Gauss meter actualized to a value of 0.3322 T. However, the FEMM software predicted an average magnet flux density of 0.3421 T, which is also measured on the surface center in
Figure 2a, and the flux density across the red line is shown in
Figure 2b. This measured value diverges from empirical value to about 4.73 %; however, the level of divergence is considered sufficiently accurate while
Figure 2b shows that higher flux occurs at the edges of the magnets.
During FEMM simulation of the coil–magnet model, a total of eight (8) magnets of
were paired into four (4) groups as shown in
Figure 1 and according to Equation (2). FEMM, however, predicted average magnetic flux density of 0.04993 T and 0.1007 T at the center of the single and paired magnets, respectively, as shown in
Figure 2b, while
Figure 2a shows the FEMM pattern on the paired magnet.
Adequate flux/coupling prediction requires insight about the distribution of the flux fields in the coils (i.e., flux density per unit volume (
)). The predicted
B values were found consistent with the magnet dimension realizing a
value of
, 2.234 ×
and
on the
and
paired and
single magnet configurations, respectively. From the above, the magnetic flux density on any NdFeB N52 permanent magnet of known volume
was obtained as
where
is the magnetic flux density per unit volume. It was obtained as the ratio of the magnetic flux density associated with each magnet geometry on the FEMM to its volume.
Considering the transducer geometry, a need arose to normalize
and
associated with a single magnet to
and
to account for the coil area where the magnetic flux density was measured in the FEMM. Hence, the normalized equation for predicting magnetic flux density in any coil geometry with volume
was obtained as
Using Equation (7), we reformulate Equation (3) to an equation as shown in Equation (8).
where
was obtained as the product of the fill factor
and effective length
. A plot of variation of the flux measured on different coil geometry with thicknesses of different components of the model is shown in
Figure 3a. From
Figure 3b a line of fit and fit equation between
and the coil thickness measured in millimeters
is shown in Equation (8).
From Equations (3), (8) and (9) an empirical relation between the magnet flux density per unit volume of the transduction coil was obtained as
Equations (8) and (10) are sufficient to make a prediction of the flux density per volume of a coil and the coupling coefficient on any coil geometry, respectively.
Using
as the reference configuration, while keeping effective length
and packing density factor
approximately equal over different width size, implies that, to ensure that
remains as shown in Equation (1), the term
and
will approximately change with each configuration according to Equation (11).
where
,
and
are the ratio of the coil width, the total coil turn and width of the
coil while
and
are the total coil turns and width of the reference coil. An equation like Equation (11) in terms of the volume ratio
also exists.
4. Results
As earlier mentioned, and are fixed thicknesses of the magnet and air space separating the coil and the magnet, respectively. During FEMM simulation of the iron–magnet–coil part of the harvester, these two parameters are fixed while and were varied to realize different models of the coils while ensuring that the variable and are equal across different realized coil configurations and that was carefully chosen to ensure that an approximately zero flux leakage occurs on the iron.
The result and legends from the FEMM simulation are respectively shown in
Figure 3a. The location of the coil corresponding to each simulation is outlined in thick red lines. A summary of the flux density
and leakage sufficient iron cladding thickness over different coil width (
) is shown in
Table 1. Different coil width
, reported in
Table 1, was achieved by winding the copper wire having a total effective length
into a coil loop of total turns
while ensuring equal
and
[
10] for all geometries.
The dotted green line in
Figure 3b shows the level of flux at which the harvesters become size (thickness) optimized in term of the flux density (
) and degree of coupling (
). This is because the overall width of the magnet–coil is the minimum corresponding to the 4 mm coil thickness. Based on this estimation, the optimized flux, coupling coefficient, coil thickness and overall transducer thickness for the model herein described was predicted at a value of 0.4373 T,
Tmm, 4.00 mm and 18.4 mm, respectively, corresponding to the intersection of the flux density on the iron cladding and the transduction coil. This optimum point corresponds to
T
. Since
is linearly dependent on
and
is an inverse and direct relationship, in the last two columns of
Table 1, respectively, therefore an inverse relation will exist between
and
contrary to expectation since it is very basic to think that the flux and harvested power will become size optimized on highest possible coupling. However, this is not the case because higher coupling will always produce an undesirable larger damping, hence a reduced harvested power [
9].
During design, it is advised to concentrate the transducer mass in the non-magnetic coil brace to ensure accuracy of flux prediction while targeting expected resonance.