Next Article in Journal
Hybrid Neural Network Reduced Order Modelling for Turbulent Flows with Geometric Parameters
Next Article in Special Issue
Irreversibility Analysis for Eyring–Powell Nanoliquid Flow Past Magnetized Riga Device with Nonlinear Thermal Radiation
Previous Article in Journal
Fluid–Structure Interaction and Non-Fourier Effects in Coupled Electro-Thermo-Mechanical Models for Cardiac Ablation
Previous Article in Special Issue
Variable Energy Fluxes and Exact Relations in Magnetohydrodynamics Turbulence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of a Symmetrical Ferrofluid Sloshing Vibration Energy Harvester

Department of Mechanical & Aerospace Engineering, North Carolina State University, Raleigh, NC 27606, USA
*
Author to whom correspondence should be addressed.
Fluids 2021, 6(8), 295; https://doi.org/10.3390/fluids6080295
Submission received: 15 February 2021 / Revised: 12 August 2021 / Accepted: 13 August 2021 / Published: 22 August 2021
(This article belongs to the Special Issue Fluids in Magnetic/Electric Fields)

Abstract

:
Ferrofluid sloshing vibration energy harvesters use ferrofluid sloshing movement as a moving magnet between a fixed coil to induce current and, in turn, harvest energy from external excitations. A symmetric ferrofluid sloshing vibration energy harvester configuration is introduced in this study which utilizes four external, symmetrically placed, permanent magnets to magnetize a ferrofluid inside a tank. An external sinusoidal excitation of amplitude 1 m/s2 is imparted, and the whole system is studied numerically using a level-set method to track the sharp interface between ferrofluid and air. The system is studied for two significant length scales of 0.1 m and 0.05 m while varying the four external magnets’ polarity arrangements. All of the system configuration dimensions are parametrized with the length scale to keep the system configuration invariant with the length scale. Finally, a frequency sweep is performed, encompassing the structure’s first modal frequency and impedance matching to obtain the system’s energy harvesting characteristics.

1. Introduction

Energy harvesting from different ambient sources of energy has become a significant area of research. Energy harvesting entails converting ambient or unused energy into useful electrical energy through a miniature transduction mechanism. It is used mainly in places where large energy conversion devices, such as hydro/steam turbines, IC engines, etc., cannot convert the available ambient energy. Energy harvesting aims to provide a sustainable power source to stand-alone devices, such as sensors, actuators, etc., by utilizing ambient energy and without needing to store energy [1,2,3]. Energy harvesting devices have the potential to replace bulky batteries and deliver lightweight solutions for sensing and actuation and make them essentially wireless [4,5,6]. It is also helpful in increasing general energy conversion devices’ efficiency (for example, thermoelectric devices, which are used for waste heat recovery from heat engines).
Energy harvesting usually requires a transduction mechanism to convert ambient/waste energy into a useful electric energy form [7,8,9]. Figure 1 depicts the working of an energy harvesting device schematically. The energy harvesting device captures ambient or waste energy through an energy capturing mechanism. This captured energy is then used by a transduction mechanism that converts it into useful energy (electric power). The most general examples are thermoelectric devices, piezoelectric devices, ionic polymer metal composites, etc. [1,2,10]. However, recently, energy harvesters have used a liquid state transduction mechanism to obtain electric power through a ferrohydrodynamic vibration energy harvester [11,12]. Fluidic transduction can be applied in energy harvesting by both solid and liquid state transducers. For solid-state transducers, the fluid acts as the energy source for the forcing function required to deform part of the harvester. In liquid state transducers, the fluid acts as both the energy carrier and deformed harvester. This is shown as a schematic in Figure 2 Solid-state transduction mechanisms use any instability created in the flow field to excite the transduction element at a particular frequency (usually the natural frequency of the structure) and harvest energy via inherent electromechanical coupling of the harvester.
Sloshing refers to the movement of liquids inside partially filled containers under external excitation [13,14]. It is characterized as an internal free surface flow. Sloshing has typically been an area of interest for the shipping industry, the automobile industry, and intergalactic rocketry [13]. Most of these applications have tried to mitigate the effects of sloshing since, at large scales, a sloshing fluid can cause significant structural damage to the container [15]. Therefore, most research about sloshing has been conducted to mitigate the effects of sloshing using baffles, etc. The solid-state nature of the transduction mechanism in most current vibration energy harvesters usually limits their capabilities, especially in conformability to different shapes and very diverse sensitivity to broadband external excitations. Both of these problems can be addressed if the transduction mechanism is fluidic instead of a solid, and this has been achieved through a Ferrofluid sloshing vibration energy harvesting system [16].
Ferrofluids consist of stable ferrous nanoparticles in colloidal suspension, forming nanoscale magnetic dipoles [17,18] once subjected to an external magnetic field. Each particle is coated with a surfactant such that in the absence of an external magnetic field, the magnetic dipoles are randomly oriented in the carrier fluid and do not coagulate due to the fields generated by the ferromagnetic cores [17]. However, when an external magnetic field is applied, the average direction of the fluid magnetization becomes parallel to the magnetic field lines. This property is known as superparamagnetism and can be used in flow manipulation applications by affecting the ferrofluid through an external magnetic field [17,18,19,20]. One such application, which was briefly discussed above, is to use ferrofluids for energy harvesting. Here the ferrofluid, due to superparamagnetism, becomes a magnet when placed in an external magnetic field. Combining this with an external excitation and a pick-up coil, the ferrofluid sloshing motion can be utilized to harvest energy. This is known as ferrofluid sloshing vibration energy harvesting. Ferrofluid sloshing has been used for generating power from external oscillations. Performance characteristics of one such energy harvesting configuration were conducted, and a plethora of knowledge regarding design rules for such a configuration was proposed [21]. Table 1 shows a list of published Ferrofluid sloshing vibration energy harvesting configurations, which used only a liquid state transduction mechanism consisting of a ferrofluid. The highest power output has been around 80 mW for a 10 m/s2 acceleration.
In addition, different configurations have used Ferrofluids to augment a solid-state transduction mechanism. For instance, Ferrofluid is used as a bearing to reduce friction between a moving solid member and a stationary shell, and power is harvested through electromagnetic induction when external excitations are supplied to the system [22,23,24,25]. In these systems, the ferrofluid does not participate as a transduction mechanism but aids in the movement of solid magnetic material in a coil.
In this study, a unique configuration for a ferrofluid sloshing vibration energy harvester (VEH) is proposed, exploiting the harvester’s two symmetry planes by introducing four symmetrically placed powerful magnets to generate sixteen unique combinations of magnetic fields. The idea behind the symmetrical placement of magnets is that it becomes easy to manufacture if the structures are symmetric. Moreover, the sloshing tank proposed in this study has a rectangular elevation, with height being half that of length. The third dimension is assumed to be equal to the length of the tank. This helps in reducing the overall size of the system and increasing the harvested power per unit volume.
The sloshing phenomena is defined numerically by the Navier-Stokes equations, the Maxwell equations, and the Level-set method, and is simulated using COMSOL Multiphysics. The level-set method is utilized to capture the movement of the free surface. The level-set method has been used to perform numerical studies on ferrofluid droplets [28,29,30] and, in this work, similar formulations are adapted to depict a sharp interface between the ferrofluid and air. For convenience and efficiency, the flow simulation is performed as a 2-D flow instead of a 3-D flow using simplifying assumptions published in the literature. Section 2 defines the set-up of the system, and Section 3 comprehends the governing equations used to simulate the multiphysics flow.

2. System Set-Up

Figure 3 represents the configuration with four magnets symmetrically placed at the center planes of the harvesting tank.
The four magnets have length dimensions in terms of the significant length scale (i.e., the tank’s length). The magnets’ thickness is a tenth of the length of the tank, and their width is a fifth of the tank. Such a formulation helps in scaling the system with the actual length of the configuration. The magnets are placed at a distance of L / 4 from their adjacent tank wall. The height of the tank ( H ) is half of its length, and the fill level of the ferrofluid in the tank is h . The magnets are labeled from 1 through 4.
For convergence of the magnetic field, a large domain is created around the harvester. This is shown in Figure 4, where the domain is ten times larger than the harvester’s length scale. The harvester receives an external excitation, which is periodic and causes the ferrofluid to slosh inside.

3. Governing Equations

The flow inside the domain is governed by Navier-Stokes equations and Maxwell’s equations. The Navier-Stokes equations are as follows
Continuity:
ρ · V = 0
Momentum:
ρ V t + ρ ( V · ) V = p + μ 2 V + F m + ρ   g + F st  
where, V is the velocity vector, F m is the magnetic force, g is the external acceleration on the system, which consists of the external excitation and the acceleration due to gravity, and F st is the force of the surface tension.
Here the external excitation is given as:
x = X 0 Sin ( ω t )
where X 0 is the amplitude of external excitation and ω the angular frequency related to the excitation frequency f as, 2 π f . In this study, f is varied such that the set of frequencies contain the first modal frequency. This makes the acceleration:
x = X 0 Sin ( ω t )
Since the flow is a separated two-phase flow, with an interface separating the two, the properties of the fluids will vary sharply across the interface. This needs to be accounted for using an interface tracking method, which will keep the boundary between the fluids intact and allow the interface to advect with changing flow field. This requires a separate computation step and a separate mathematical treatment/modeling of the moving interface. Several methods are very effective in performing interface tracking computationally such as the volume of fluid method, the Level-set method, the phase-field method, etc. [31]. These methods usually work based on conservation of some physical or attributed numerical quantity, usually the mass or the area under the surface, the gradient of the distance function, etc., to characterize the free surface. This then is used to derive a partial differential equation form of the conserved quantity and is solved with the N-S equations to yield the moving interface profile. One such method mentioned above is the level-set method [32]. The interface is defined as a signed distance function in the level set method and is tracked as a higher dimensional field applied to the flow field, which in turn defines how the flow moves [33,34,35]. However, the level set method proposed initially [34] did not conserve mass or the area correctly, especially in fluid dynamics problems. This prompted further research into conservative level set methods that conserve the mass exactly. One such method is used in this work where the mass is conserved exactly, and the surface is expressed as a constant level set [36,37]. Although level set methods are inherently well suited for representing complex topologies, they become even more effective in tracking sharp interfaces with a conservative level-set formulation [38,39,40]. The mathematical formulation of this sharp interface, the conservative level-set method, is shown below.
Level-set equations for interface tracking:
ϕ t + V · ϕ = γ   · ( ϵ ls ϕ ϕ ( 1 ϕ ) ϕ | ϕ | )
where ϕ is the level set variable, γ is the re-initialization parameter and ϵ ls is the parameter which controls the thickness of the fluid-air interface. In this work, the value of the re-initialization parameter is set to X 0 ω , which is the characteristic velocity scale for the sloshing setup. The separating surface is given by ϕ = ϕ 0 = 0.5 , i.e., a set of all surfaces where the level set function ϕ = 0.5 . This means that the first fluid is from ϕ = 0   to   ϕ = 0.5 and the second fluid from ϕ = 0.5   to   ϕ = 1 . This level-set variable creates a higher dimensional field over the flow field and describes all the properties of the two fluids in its terms to represent the separation between the two fluids as a sharp interface. All fluid properties must be represented in terms of the level set variable to fully couple the flow field with the level-set field. Similarly, all interfacial forces must be expressed in terms of the level-set variable.
The surface tension force is expressed in terms of the level-set variable as follows:
F st = σ m δ κ n i + δ   s σ m
where, σ m is the surface tension coefficient, δ is a Dirac delta function located at the interface, κ is the curvature of the surface, n i is the normal vector of the liquid-gas interface, and s is the surface gradient operator.
The interface normal is given by:
n i = ϕ | ϕ |
The Dirac delta function is given by:
δ = 6   | ϕ |   | ϕ ( 1 ϕ ) |
The curvature of the surface is given by:
κ = · n i
The surface gradient operator is given by:
s = ( I n i n i T )
Additionally, the angle that the fluid makes with the wall while wetting is assumed to be 90°. This is a standard assumption for solving sloshing flows [27,41].
For the fluid domain inside the tank, the properties must be expressed in terms of the level set variable. Density is given as:
ρ = ρ 1 + ( ρ 2 ρ 1 ) ϕ  
Viscosity is given as:
μ = μ 1 + ( μ 2 μ 1 ) ϕ
The indices 1 and 2 denote fluid 1 and fluid 2, which is the ferrofluid and air, respectively. Similarly, to calculate the Kelvin body force, the magnetic properties of the fluids must be expressed in terms of the level set variable:
F m = M · H
where, H is the magnetic field, and the magnetization M is given by the constitutive relation:
M = χ m H
where, χ m is the magnetic susceptibility and is obtained from the level set function for the whole domain as:
χ m = χ m 1 + ( χ m 2 χ m 1 ) ϕ
Maxwell’s conservation laws for the system are [17]:
· B = 0
× H = 0
× E = B t
The magnetic field can also be written as:
H = V m
where, V m is the magnetic scalar potential.
The boundary conditions for the domain are the no flux boundary condition for the ends of the domain surrounding the harvester. The magnets are such that they have a uniform magnetization of 400 kA/m. The following are the constitutive relations for different regions:
B = μ 0 H ( for   air )   ; B = μ 0 ( H + M s )   ( for   the   magnets )   & B = ( 1 + χ m ) H   for   sloshing
The electromotive force from sloshing is picked up by an external coil with N turns (which are varied to optimize power harvested). The electromotive force is given as:
ϵ = N d Φ dt
where Φ is the magnetic flux, obtained by integrating the magnetic flux density over the area. The electromotive force due to the changing magnetic field inside the harvester is picked up by an external coil, which can be wound either along the harvester’s length and gets pitched in the y-direction, or it can be wound along the height of the tank while pitched along the x-direction. Equation (21) represents the expression of inductance of one such rectangle, where μ r is the relative permeability of the material enclosed in the coil.
L rect = N 2 μ 0 μ r π [ 2 ( W + H ) + 2 H 2 + W 2 H ln ( H + H 2 + W 2 W ) W ln ( W + H 2 + W 2 H ) + Hln ( 2 H d 2 ) + Wln ( 2 W d 2 ) ]
The current in the wound coil can be found by the series L-R circuit equation:
L c di dt + i   ( R L + R c ) = V  
where R L is the load resistance and R C the coil resistance, through which the power output of the harvester is measured, L c = L rect is the inductance of the wound coil, and V is the electromotive force obtained from sloshing. Finally, it is assumed that if a coil is wound along the length (i.e., horizontally wound), the number of turns would be half that of when it is wound along the height (i.e., vertically wound coil). Figure 5 shows the 3-D arrangement of the tank with a horizontally wound coil and Figure 6 shows the arrangement of the tank in 3-D for a vertically wound coil. The four magnets are also shown in the figures, with the depth dimension being equal to that of the height.
For vertical winding (along with the height), the magnetic flux’s x component would be used. For horizontal winding (along with the length), the y component of the magnetic flux would be used to calculate the electromotive force. This is formulated in Equation (23).
ϵ vertical = N vertical d Φ x dt   and   ϵ horizontal = N horizontal d Φ y dt

4. Numerical Modeling and Validation

In this section, the numerical modeling details are discussed. The combined Navier-Stokes, level-set, and Maxwell’s equations are solved in COMSOL Multiphysics using its inbuilt multiphysics capacity. COMSOL is a finite element-based software, which solves partial differential equations using the weak form and shape functions. For computational fluid dynamics applications, especially where Navier-Stokes equations are involved, finite element-based codes need to account for stability errors in the crosswind and streamline direction due to numerical dispersion. This is especially significant when the same order of elements resolves velocity and pressure fields. These stabilizations are inbuilt in COMSOL and are based on Streamline upwinding methods for both pressure and velocity fields [42,43]. Further, as the problem is unsteady, the solution must be marched in time to obtain the system’s state at each time of interest during the simulation. For this, an implicit adaptive time-stepping scheme is used.
The solution is advanced by solving segregated steps at each time step when the tolerance level (which is set to be 1 × 10−4) is reached. This convergence criterion is for both solution and the residual. The following segregated steps of dependent variables are solved in the following order:
(a)
Level set variable
(b)
Velocity and pressure
(c)
Magnetic scalar potential.
To simplify and streamline the computations, the first-order discretization is used for the pressure and velocity fields. Similarly, for the level-set, variable linear elements are used. For the magnetic scalar potential, the element order is maintained at two everywhere in the domain.
Grid convergence studies and adaptive time step control studies were performed on two magnet configurations. It was found that the mesh inside the sloshing domain is both accurate and efficient for a grid size of 0.001 and for a maximum time step value of 2 milliseconds.
To validate COMSOL and the level set method used in the computations, a comparative study was run replicating the published experimental and computational results by Liu et al. [44]. The combination 1H is simulated, and the results are compared in Figure 7. The numerical scheme predictions adopted in the present work match very closely to the reported computational and experimental data in [44] as shown in Figure 7a,b, respectively. The peak voltage near the resonance frequency almost matches perfectly with the experimental data. However, since the spectrum of frequencies in the frequency sweep performed in the present work is not as closely spaced as reported by Liu et al., the internal energy pump, which occurs at around 1.5 Hz, is not as prominent in our simulation of their configuration as shown in Figure 7b.
Based on this validation, it is inferred that the level-set formulation and the chosen solvers can simulate the complexities involved in this flow. To advance the work on ferrofluid VEH, a new tank design is proposed in this work, as mentioned in previous sections. This configuration is then subjected to a frequency sweep, with various combinations between the four fixed magnets. The north and south poles of one or more magnets are flipped to generate a new magnetic field distribution. Table 2 lists all of the 16 cases which may arise with different combinations of North and South pole orientations of the magnets. The table lists South/North labels for the magnets, depicting which pole of the magnet faces the tank.
Two tanks are studied, one where L is 10 cm and the other where L is 5 cm. The fill level is a fraction of the length scale and is set to be L/5. For the two tanks, the frequency sweep is performed to capture the tank’s first resonance frequency within the sweep.
Finally, to maintain energy application parity across tank sizes and frequencies, the external acceleration (a) is kept as constant, at 1   m / s 2 .
i.e., a = X 0 ω 2 = 1
Hence when the frequency changes, the amplitude is adjusted to keep the acceleration constant.
To provide additional validation, a video (Video S1) illustrating the surface movement of the ferrofluid from the 2-D simulation for case 10 and frequency 2.05 Hz for the large tank is provided in the supplementary materials. To compare this 2-D simulation with a 3-D full-scale simulation, a 3-D simulation of the whole tank was performed, with the arrangement shown in Figure 5 and Figure 6. To compare this simulation with the 2-D simulation, a mid-plane video (Video S2) is taken from the 3-D simula-tion. The surface evolution of the ferrofluid in these two videos, over 4 time periods, is very similar. Hence, it is safe to assume that the 2-D simulation is a very good ap-proximation of the actual 3-D simulation.
Finally, the properties used for this study are of ferrofluid EFH3 manufactured by Ferrotec corporation [45] and are listed in Table 3.

5. Results and Discussion

In this section, simulation results and the performance of the two energy harvesters are discussed. An impedance matching scheme was followed to extract results for the highest extractable power from all these simulation runs for both tanks. Here, the external electrical circuit is considered a series R-L circuit, and the equations for this are formulated in the previous section. To determine the optimal impedance for the system, many pick-up coils were selected, and those pick-up coils were assumed to be connected to a load resistor, the resistance of which is varied from 1 ohm to 105 ohms to generate an impedance matching plot. This was performed for all the frequencies and both the tanks.
For example, Figure 8 and Figure 9 show impedance matching plots for the L = 10 cm tank for the magnet orientation case 7 for frequency of 2.05 Hz using horizontal and vertical pick-up coils of 1000 turns. For the horizontal coil, the optimal load resistor was found to be 103.5 ohms, and for the vertical coil 102.75 ohms. Appendix B shows the time evolution of voltage, current, and external acceleration for this 2.05 Hz case.
The following two subsections present simulation runs, and impedance matching results for the two tanks studied in this work.

5.1. Analysis of Output from L = 10 cm Tank

Figure 10 shows how the ferrofluid moves inside the tank near the resonance frequency for the magnet combination case 1. The rest of the 15 cases are displayed in Appendix A.
The arrows in the above plot depict the magnetic flux lines. The flux lines in a magnet originate from the North Pole and terminate into the South Pole. In the above plot, magnet 1′s South Pole faces the harvester, while magnet 2′s north, magnet 3′s south, and magnet 4′s north, face the harvester. Now, for each of these 16 cases of magnet arrangement, a frequency sweep is performed. Each frequency sweep produces a voltage signal from the harvester. This is then transformed into RMS voltage to quantify the performance characteristics of the harvester. Now to study the maximum performance that the harvester is capable of, nine coils are modeled with different numbers of turns.
Additionally, the power harvested from the harvester cannot be the power that goes directly to the pick-up coil. A load resistor is attached to the circuit, and the power through the load resistance is calculated to measure the harvester’s performance. This forms a series L-R circuit, and Equations (20)–(22) are used to model this circuit. Secondly, to check the maximum power output in the load resistor for the same temporal variation in magnetic flux density, 25 load resistors are modeled. These resistors, combined with the nine coils, are used to estimate peak voltage (RMS) and peak power output (RMS). These two outputs are not necessarily linked to each other, as the peak power may or may not correspond to the peak voltage reported. This is performed to analyze the voltage and power output potential from the two tanks separately. Of course, the ultimate performance can only be analyzed through the power output from the system. However, that is also a function of the circuit one builds to harvest the system’s energy. In our case, it is a series L-R circuit represented by Equation (22). The calculation of the coil resistance is based on the number of turns the pickup coil has. The coil material is chosen to be copper. Table 4 summarizes the properties of the pick-up coils modeled for L = 10 cm tank.
The coil properties are calculated by assuming that the tank depth is 10 cm and both the coils envelop the length along the winding direction completely. This gives the coil diameter for the horizontal coil as H N and for the vertical coil as L N , where H is the height of the tank and L is the length of the tank. The tank depth is maintained to be equal to the tank’s largest dimension to discount 3-D effects on the 2-D simulation performed in this study. To calculate the coil’s inductance from Equation (21), the relative permeability is calculated by taking a weighted average of the ferrofluid and air’s relative permeabilities. This comes out to be 2, considering air as 1 and ferrofluid as 4.52 ( 1 + χ m ). Finally, to check the power output performance of each coil 25 load resistors ( R L ) are used, which is usually known as impedance matching. The load resistances are varied from 1 to 105 ohms, with an increment of 100.25 ohms. This covers the total range of coil resistances shown in Table 4.
Figure 11 shows the performance measured by max RMS voltage output when a horizontally wound coil is used for voltage pickup from the tank. These plots are for all the 16 magnet configurations, and almost all the cases have low voltages and have almost the same range. The max voltage noted is around 0.03 V, and it is from magnet case 6 and magnet case 11; both seem identical. It must be noted that the cases seem to be symmetric in their behavior, hence making a pair of eight cases, which must be valid for an oscillatory motion. The RMS voltages are calculated after the first two time periods have elapsed. At this point, it is assumed that the system has equilibrated from the initial jerk, which propels the system in an oscillatory motion. The simulations were run for five time periods, and hence the RMS quantities used are from the last three time periods.
For, vertical winding the RMS voltages are shown in Figure 12. They are an order of magnitude higher than the horizontal ones. Cases 7 and 10 achieved the highest RMS voltage output of around 0.15 V. Upon close observation, cases 7 and 10 look very similar to each other solely from the magnets’ placement point of view. This is similar to what was observed for max RMS voltages for the horizontal coil. This can be attributed to symmetric magnet placement in cases 7 & 10 and all the other pairs of cases. The symmetric magnet placement gives rise to magnetic field profiles inside the tank, which are mirror images of each other, which can be deduced from the direction of the magnetic flux lines for these two cases. Secondly, it is observed that there is a secondary peak near 1.5 Hz for both the horizontal and vertical windings. This is due to the internal coupling between modes, generating a non-linear internal energy pump [26,27]. The peak is smeared off since the frequencies in our sweep are equidistant from each other.
The RMS current is calculated by solving for instantaneous current from Equation (22) for all of the 9 coils and 25 load resistors to determine the power through the load resistor. Then, the RMS value of this current is taken, and power through the load resistance is calculated using I r m s 2 × R L using the value of RMS current pertaining to that particular load resistor and that particular coil. The maxima of the power output for the above 225 cases (i.e., 9 coils and 25 load resistors) is taken for each of the 16 cases of magnet arrangements and all of the 27 frequencies. Figure 13 shows the maximum power through the horizontal winding for all 16 cases, and all of the 27 frequencies are plotted. The maximum power obtained is 0.15 micro-watts ( μ W ) or 0.015 ( μ W g ) for case 6 and case 11. Here too the symmetry observed for the max voltages is evident. The maximum power value is obtained for a frequency of 2.3 Hz, which points towards shifting the peak from the first modal frequency of 2.0865 Hz. This can be attributed to the magnetic body force exerted on the fluid in the y-direction. Since the magnetic field in the y-direction is considered for electromotive force through the horizontal winding, the shift in the occurrence of the peak power output from the modal frequency is very evident. Hence, the frequency for peak power shifts to the right of the modal frequency for the horizontal winding.
For the vertical coil, the max RMS power for all of the 27 frequencies and 16 magnet cases is plotted in Figure 14. The maximum power is obtained at 2.05 Hz for case 7 and is almost near to 20 micro-watts ( μ W ) or 2 ( μ W g ) , which is two orders of magnitude higher than that from the horizontal coils. It must be noted that the frequency where the max power output is observed is almost near to the modal frequency of 2.0865 Hz. The shift in the peak occurs towards the left from the modal frequency and is very small compared to that in the horizontal winding case. This can be attributed to the fact that the vertical coil’s output voltage is a function of the magnetic field in the x-direction. The range of the max power extracted using different magnet configurations is very high. The max power out is as low as 0.5 ( μ W ) to the peak of 20 ( μ W ). This indicates that placing magnets around the tank in a certain way facilitates the power output. It must be restated that the magnets’ position relative to the tank does not change; only the polarity changes. With that change, almost two orders of magnitude of difference are observed in the max power harvested from the device’s vertical winding. This is also true for the horizontal winding case, but the gap between the lowest and highest max powers extracted is not as high as in the case of vertical winding. This provided additional design rules to harness this work into enhancing designs for ferrofluid sloshing vibration energy harvesting systems.

5.2. Analysis of Output from L = 5 cm Tank

In this subsection, the tank, which is half the previously discussed tank’s length, is studied. For this tank, the frequencies are varied from 2.2 Hz to 3.5 Hz, which passes through the first modal frequency of 2.96 Hz. The 16 cases of magnet configuration as mentioned in orientation remain the same. Since the length scale L parametrizes all the system dimensions, the whole system’s length becomes half of what it was for the tank studied earlier. This way, the present system becomes identical to the last tank but with a change in the significant length scale. Figure 15 illustrates the surface evolution of the ferrofluid for case 1, for frequency nearest to the resonant frequency. For case 1, the north face of magnet 2 and 4 face the tank while the south face of magnet 1 and 3 face the tank. The plots for the other 15 cases are placed in Appendix A.
Now, for each of these 16 cases of magnet arrangement, a frequency sweep is performed. Each frequency sweep produces a voltage signal from the harvester. This is then transformed into RMS voltage to quantify the performance characteristics of the harvester. Now to study the maximum performance that the harvester is capable of, nine coils are modeled with a different number of turns. Also, the power harvested from the harvester cannot be the power that goes directly to the pick-up coil. A load resistor is attached to the circuit, and the power through the load resistance is calculated to measure the harvester’s performance. This forms a series L-R circuit, and Equations (20)–(22) are used to model this circuit. Secondly, to check the maximum power output in the load resistor for the same temporal variation in magnetic flux density, 25 load resistors are modeled. These resistors combined with nine coils are used to estimate peak voltage (RMS) and peak power output (RMS). These two outputs are not necessarily linked to each other, as in, the peak power may or may not correspond to the peak voltage reported. This is done to analyze the voltage and power output potential from the two tanks separately. Of course, the ultimate performance can only be analyzed through the system’s power output; however, that is also a function of the kind of circuit one builds to harvest the energy from the system. In our case, it is a series L-R circuit represented by Equation (22).
The calculation of the coil resistance is based on the number of turns the pickup coil has. The coil material is chosen to be Copper. Table 5 summarizes the properties of the pick-up coils modeled for L = 5 cm tank. Since the tank’s dimension is halved, the resistance values and the nine coils’ inductance values will also change. The only constants here are the number of turns in the coils and the properties of the magnets used. The coil properties are calculated by assuming that the tank depth is 5 cm, i.e., equal to the length (same as for L = 10 cm tank), and both the coils (Horizontal and Vertical) envelop the length along the winding direction completely. This gives the coil diameter for the horizontal coil as H N and for the vertical coil as L N , where H is the height of the tank and L is the length of the tank. To calculate the coil’s inductance from Equation (21), the relative permeability is calculated by taking a weighted average of ferrofluid and air’s relative permeabilities. This comes out to be 2, considering air as 1 and ferrofluid as 4.52 ( 1 + χ m ). Finally, to check the power output performance of each coil 25 load resistors ( R L ) are used, which is usually known as impedance matching. The load resistances are varied from 1 to 105 ohms, with an increment of 100.25 ohms.
Figure 16 the performance measured by max RMS voltage output is displayed when a horizontally wound coil is used for voltage pickup from the tank. These plots are for all the 16 magnet configurations, and almost all the cases have low voltages and have almost the same range. The max voltage noted is around 0.008 V, and it is from magnet case 6 and magnet case 11, like in the case of L = 10 cm tank. However, the harvested voltage is an order of magnitude smaller as compared to the last tank. It must be noted that the cases seem to be symmetric in their behavior, hence making a pair of eight cases, as was observed in the previous tank. The RMS voltages are calculated after the first two time periods have elapsed, at which point it is assumed that the system has equilibrated from the initial jerk, which propels the system in an oscillatory motion. The simulations were run for five time periods, and hence the RMS quantities used are from the last three time periods.
For, vertical winding the RMS voltages are shown in Figure 17. They are an order of magnitude higher than the max RMS voltage for horizontal coil, as was observed for L = 10 cm tank. Cases 7 and 10 achieved the highest RMS voltage output of around 0.038 V. The symmetric magnet placement in cases 7 and 10 causes the max RMS voltage observed to be identical. As mentioned in the previous sub-section, the symmetric magnet placement gives rise to magnetic field profiles inside the tank, which are mirror images of each other, which can be deduced from the direction of the magnetic flux lines for these two cases. Similarly, this applies to the other 14 configurations. It can be observed from the voltages and power plots that eight pairs of cases show identical behavior. This would be true for any such system subjected to oscillatory external excitation.
Secondly, it is observed that there is a secondary peak that was observed for the L = 10 cm tank near 1.5 Hz for both the horizontal and vertical windings, which does not seem to appear for the L = 5 cm tank. This might be due to the inexistence of any non-linear coupling between the sloshing modes of the system, perhaps due to the reduction in size and the increase in the modal frequency from 2 Hz (for L = 10 cm tank) to 3 Hz (for L = 5 cm tank).
The RMS current is calculated by solving for instantaneous current from Equation (22) for all the 9 coils and 25 load resistors to determine the power through the load resistor. Then the RMS value of this current is taken, and power through the load resistance is calculated using I r m s 2 × R L , using the value of RMS current pertaining to that particular load resistor and that particular coil. Then the maxima of the power output for the above 225 cases (i.e., 9 coils and 25 load resistors) is taken for each of the 16 magnet arrangements and all the 27 frequencies. Figure 18 shows the maximum power through the horizontal winding for all 16 cases, and all of the 27 frequencies are plotted.
The maximum power obtained is 0.005 micro-watts ( μ W ) or 0.5 ( n W g ) for case 6 and case 11. Here too, the symmetry observed for the max voltages is evident. The maximum power value is obtained for a frequency of 2.5 Hz, pointing towards shifting the peak from the first modal frequency of 2.96 Hz. This again can be attributed to the magnetic body force exerted on the fluid in the y-direction. Here, the frequency for peak power shifts to the left of the modal frequency for the horizontal winding, unlike the L = 10 cm tank case where the frequency shifted towards the right.
For the vertical coil, the max RMS power for all the 27 frequencies and 16 magnet cases is plotted in Figure 19. The maximum power is obtained at 2.8 Hz for case 7 (and case 10) and is almost near to 0.7 micro-watts ( μ W ) or 70 ( n W g ) , which is two orders of magnitude higher than that from the horizontal coils. It must be noted that the frequency where the max power output is observed is nearer to the modal frequency of 2.95 Hz than the horizontal coils. The shift in the peak occurs towards the left from the modal frequency and is very small compared to that in the horizontal winding case. This can be attributed to the fact that the vertical coil’s output voltage is a function of the magnetic field in the x-direction. The range of the max power extracted using different magnet configurations is very high for the vertical winding, as was observed for L = 10 cm case. With that change, almost two orders of magnitude of difference are observed in the max power harvested from the device’s vertical winding. This can also be attributed to the fact that the system’s excitation is in the x-direction. This, coupled with the gravity in the y-direction, makes the harvesting of energy from changes in the magnetic field in the y-direction unattractive. This also makes the power harvested from horizontal windings from all of the 16 cases of magnet arrangements to have a minor variance. This can be useful while designing systems with both x and y excitations or perhaps angular excitation or combining all of the three.

6. Conclusions

A completely new design for a ferrofluid sloshing vibration energy harvester is introduced in this work. The design is inherently symmetric since its purpose is to exploit symmetry that can be easy to manufacture. Two tanks, one of length 10 cm and one of length 5 cm, are observed by providing a constant external acceleration of 1 m/s2. The tanks are then subjected to a frequency sweep for 16 different combinations of the four magnets by changing their polarity.
The most important results observed are as follows:
(a)
The maximum RMS voltages harvested reached almost 0.15 V
(b)
The maximum power after impedance matching was obtained from 10 cm tank, which is around a value of 20 µW or 2 ( μ W g ) .
(c)
For a 5cm tank, the highest power output was 0.7 µW, which is higher than what was obtained from the horizontal winding for L = 10 cm tank.
These observations will have a significant impact on the future design of Ferrofluid sloshing VEH’s. Usually, sloshing motion decreases with the decrease in the significant length scale. At low scales, the flow converts into a capillary flow. Hence, miniaturizing these sloshing energy harvesting systems after a particular extent would involve huge trade-offs. However, from the third observation, it can be concluded that miniaturization can still lead to good performance if conducted within limits and if the magnets are oriented optimally, to yield the highest result. Finally, these observations pave the way for future research into miniaturizing the systems. Perhaps a regime like the ‘high-frequency, low amplitude’ might be suitable for harvesting energy from miniaturized Ferrofluid VEH’s. Finally, it is also to be noted that for maintaining the viability of the 2-D simulation as an exact representation of a 3-D tank, the depth of the tank was assumed to be equal to the longest edge of the tank. This increased the wire’s length in winding around the tank as a coil and reduced the harvested power. The power harvested can be increased significantly when the coil is optimized or when a different ‘pick-up’ mechanism is used, which is more efficient than an inductive circuit. Hence, optimizing such systems’ electrical conditioning efficiency would be necessary for future studies to arrive at the most optimal designs for these harvesters.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/fluids6080295/s1.

Author Contributions

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

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Most of the relevant data has been presented in this manuscript. Raw data can only be provided at the discretion of the authors due to proprietary constraints. Requests for raw data can be made to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Surface Evolution for L = 10 cm and L = 5 cm Tank

Figure A1. Case 2.
Figure A1. Case 2.
Fluids 06 00295 g0a1
Figure A2. Case 3.
Figure A2. Case 3.
Fluids 06 00295 g0a2
Figure A3. Case 4.
Figure A3. Case 4.
Fluids 06 00295 g0a3
Figure A4. Case 5.
Figure A4. Case 5.
Fluids 06 00295 g0a4
Figure A5. Case 6.
Figure A5. Case 6.
Fluids 06 00295 g0a5
Figure A6. Case 7.
Figure A6. Case 7.
Fluids 06 00295 g0a6
Figure A7. Case 8.
Figure A7. Case 8.
Fluids 06 00295 g0a7
Figure A8. Case 9.
Figure A8. Case 9.
Fluids 06 00295 g0a8
Figure A9. Case 10.
Figure A9. Case 10.
Fluids 06 00295 g0a9
Figure A10. Case 11.
Figure A10. Case 11.
Fluids 06 00295 g0a10
Figure A11. Case 12.
Figure A11. Case 12.
Fluids 06 00295 g0a11
Figure A12. Case 13.
Figure A12. Case 13.
Fluids 06 00295 g0a12
Figure A13. Case 14.
Figure A13. Case 14.
Fluids 06 00295 g0a13
Figure A14. Case 15.
Figure A14. Case 15.
Fluids 06 00295 g0a14
Figure A15. Case 16.
Figure A15. Case 16.
Fluids 06 00295 g0a15
Figure A16. Case 2.
Figure A16. Case 2.
Fluids 06 00295 g0a16
Figure A17. Case 3.
Figure A17. Case 3.
Fluids 06 00295 g0a17
Figure A18. Case 4.
Figure A18. Case 4.
Fluids 06 00295 g0a18
Figure A19. Case 5.
Figure A19. Case 5.
Fluids 06 00295 g0a19
Figure A20. Case 6.
Figure A20. Case 6.
Fluids 06 00295 g0a20
Figure A21. Case 7.
Figure A21. Case 7.
Fluids 06 00295 g0a21
Figure A22. Case 8.
Figure A22. Case 8.
Fluids 06 00295 g0a22
Figure A23. Case 9.
Figure A23. Case 9.
Fluids 06 00295 g0a23
Figure A24. Case 10.
Figure A24. Case 10.
Fluids 06 00295 g0a24
Figure A25. Case 11.
Figure A25. Case 11.
Fluids 06 00295 g0a25
Figure A26. Case 12.
Figure A26. Case 12.
Fluids 06 00295 g0a26
Figure A27. Case 13.
Figure A27. Case 13.
Fluids 06 00295 g0a27
Figure A28. Case 14.
Figure A28. Case 14.
Fluids 06 00295 g0a28
Figure A29. Case 15.
Figure A29. Case 15.
Fluids 06 00295 g0a29
Figure A30. Case 16.
Figure A30. Case 16.
Fluids 06 00295 g0a30

Appendix B. Waveform Plots for Acceleration, Voltage & Current

Figure A31. External Excitation.
Figure A31. External Excitation.
Fluids 06 00295 g0a31
Figure A32. Voltage evolution for horizontal coil with 1000 turns for case 7 of magnet arrangement and 2.05 Hz frequency of excitation.
Figure A32. Voltage evolution for horizontal coil with 1000 turns for case 7 of magnet arrangement and 2.05 Hz frequency of excitation.
Fluids 06 00295 g0a32
Figure A33. Current evolution for horizontal coil with 1000 turns for case 7 of magnet arrangement and 2.05 Hz frequency of excita-tion.
Figure A33. Current evolution for horizontal coil with 1000 turns for case 7 of magnet arrangement and 2.05 Hz frequency of excita-tion.
Fluids 06 00295 g0a33
Figure A34. Voltage evolution for vertical coil with 1000 turns for case 7 of magnet arrangement and 2.05 Hz frequency of excitation.
Figure A34. Voltage evolution for vertical coil with 1000 turns for case 7 of magnet arrangement and 2.05 Hz frequency of excitation.
Fluids 06 00295 g0a34
Figure A35. Current evolution for vertical coil with 1000 turns for case 7 of magnet arrangement and 2.05 Hz frequency of excitation.
Figure A35. Current evolution for vertical coil with 1000 turns for case 7 of magnet arrangement and 2.05 Hz frequency of excitation.
Fluids 06 00295 g0a35

References

  1. Briand, D.; Yeatman, E.; Roundy, S. Micro Energy Harvesting; Wiley-VCH: Hoboken, NJ, USA, 2015. [Google Scholar]
  2. Niell, A.; Elvin, E. Advances in Energy Harvesting Methods; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  3. Erturk, A.; Inman, D.J. Piezoelectric Energy Harvesting; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  4. Mitcheson, P.D.; Yeatman, E.M.; Rao, G.K.; Holmes, A.S.; Green, T.C. Human and machine motion for wireless electronic devices. Proc. IEEE 2008, 96, 1457–1486. [Google Scholar] [CrossRef] [Green Version]
  5. Wang, S.; Lin, L.; Lin, Z. Triboelectric nanogenerators as self-powered active sensors. Nano Energy 2015, 11, 436–462. [Google Scholar] [CrossRef] [Green Version]
  6. Roundy, S.; Wright, P.K. A piezoelectric vibration based generator for wireless electronics. Smart Mater. Struct. 2004, 13, 1131–1142. [Google Scholar] [CrossRef] [Green Version]
  7. Challa, V.R.; Prasad, M.G.; Shi, Y.; Fisher, F.T. A vibration energy harvesting device with bidirectional resonance frequency tunability. Smart Mater. Struct. 2008, 17, 15035. [Google Scholar] [CrossRef]
  8. Thambi, N.; Sastry, A.M. Powering MEMS portable devices—A review of non-regenerative and regenerative power supply systems with special emphasis on piezoelectric energy. Smart Mater. Struct. 2008, 17, 043001. [Google Scholar] [CrossRef] [Green Version]
  9. Cottone, F.; Vocca, H.; Gammaitoni, L. Nonlinear energy harvesting. Phys. Rev. Lett. 2009, 102, 080601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Erturk, A.; Vieira, W.G.R.; de Marqui, C.; Inman, D.J. On the energy harvesting potential of piezoaeroelastic systems. Appl. Phys. Lett. 2010, 96, 184103. [Google Scholar] [CrossRef]
  11. Bibo, A.; Masana, R.; King, A.; Li, G.; Daqaq, M.F. Electromagnetic ferrofluid-based energy harvester. Phys. Lett. A 2012, 376, 2163–2166. [Google Scholar] [CrossRef]
  12. Seol, M.; Jeon, S.; Han, J.; Choi, Y. Ferrofluid-based triboelectric-electromagnetic hybrid generator for sensitive and sustainable vibration energy harvesting. Nano Energy 2017, 31, 233–238. [Google Scholar] [CrossRef]
  13. Ibrahim, R.A. Liquid Sloshing Dynamics; University Press: Cambridge, MA, USA, 2005. [Google Scholar]
  14. Price, W.G.; Chen, Y.G. A simulation of free surface waves for incompressible two-phase ows using a curvilinear level set formulation. Int. J. Numer. Methods Fluids 2006, 51, 305–330. [Google Scholar] [CrossRef]
  15. Bautista-Jacobo, J.L.; Rodríguez-Morales, E.; Montes-Rodríguez, J.J.; Gámez-Cuatzín, H. Effect of baffles on the sloshing in road tankers carrying LPG: A comparative numerical study. Math. Probl. Eng. 2015, 2015, 359470. [Google Scholar] [CrossRef]
  16. Wang, D.; Chang, K. Electromagnetic energy harvesting from flow induced vibration. Microelectron. J. 2010, 41, 356–364. [Google Scholar] [CrossRef]
  17. Rosensweig, R.E. Ferrohydrodynamics; Dover Publications: Garden City, NY, USA, 1970. [Google Scholar]
  18. Hosking, R.J.; Dewar, R.L. Fundamental Fluid Mechanics and Magnetohydrodynamics; Springer: Singapore, 2016. [Google Scholar]
  19. Anand, N. Application of Ferrohydrodynamic Flows for Heat Transfer Enhancement and Energy Harvesting. Ph.D. Thesis, North Carolina State University, Raleigh, NC, USA, 2020. [Google Scholar]
  20. Anand, N.; Gould, R. Study of Enhanced Self Mixing in Ferrofluid Flow in an Elbow Channel Under the Influence of Non-Uniform Magnetic Field. In Proceedings of the 8th Joint Fluids Engineering Conference, ASME-JSME-KSME 2019, San Francisco, CA, USA, 28 July–1 August 2019. [Google Scholar] [CrossRef]
  21. Liu, Q.; Daqaq, M.F.; Li, G. Performance analysis of a ferrofluid-based electromagnetic energy harvester. IEEE Trans. Magn. 2018, 54, 1–14. [Google Scholar] [CrossRef]
  22. Wang, Y.; Zhang, Q.; Zhao, L.; Tang, Y.; Shkel, A.; Kim, E.S. Vibration energy harvester with low resonant frequency based on flexible coil and liquid spring. Appl. Phys. Lett. 2016, 109, 203901. [Google Scholar] [CrossRef]
  23. Wang, Y.; Zhang, Q.; Zhao, L.; Kim, E.S. Non-resonant electromagnetic broad-band vibration-energy harvester based on self-assembled ferrofluid liquid bearing. J. Microelectromechanical Syst. 2017, 26, 809–819. [Google Scholar] [CrossRef]
  24. Wu, S.; Luk, P.C.K.; Li, C.; Zhao, X.; Jiao, Z. Investigation of an electromagnetic wearable resonance kinetic energy harvester with ferrofluid. IEEE Trans. Magn. 2017, 53, 9. [Google Scholar] [CrossRef]
  25. Li, C.; Wu, S.; Luk, P.C.K.; Gu, M.; Jiao, Z. Enhanced bandwidth nonlinear resonance electromagnetic human motion energy harvester using magnetic springs and ferrofluid. IEEE/ASME Trans. Mechatron. 2019, 24, 710–717. [Google Scholar] [CrossRef] [Green Version]
  26. Alazemi, S.F.; Bibo, A.; Daqaq, M.F. A ferrofluid-based energy harvester: An experimental investigation involving internally-resonant sloshing modes. Eur. Phys. J. Spec. Top. 2015, 224, 2993–3004. [Google Scholar] [CrossRef]
  27. Alazmi, S.; Xu, Y.; Daqaq, M.F. Harvesting energy from the sloshing motion of ferrofluids in an externally excited container: Analytical modeling and experimental validation. Phys. Fluids 2016, 28, 077101. [Google Scholar] [CrossRef] [Green Version]
  28. Ki, H. Level set method for two-phase incompressible flows under magnetic fields. Comput. Phys. Commun. 2010, 181, 999–1007. [Google Scholar] [CrossRef]
  29. Capobianchi, P.; Lappa, M.; Oliveira, M.S.N. Deformation of a ferrofluid droplet in a simple shear flow under the effect of a constant magnetic field. Comput. Fluids 2018, 173, 313–323. [Google Scholar] [CrossRef] [Green Version]
  30. Anand, N.; Gould, R. Modeling Ferrofluid Sloshing Vibration Energy Harvesting Using Level-Set Method. 2020. Available online: https://www.comsol.co.in/paper/download/854041/COMSOLCOnference_2020_NAdishANand.pdf (accessed on 1 June 2021).
  31. Elgeti, S.; Sauerland, H. Deforming fluid domains within the finite element method. Arch. Comput. Methods Eng. 2014, 23, 323–361. [Google Scholar] [CrossRef] [Green Version]
  32. Sethian, J.A. Level Set Techniques for Tracking Interfaces: Fast Algorithms, Multiple Regions, Grid Generation and Shape/Character Recognition. 1994. Available online: http://www.ann.jussieu.fr/~frey/papers/levelsets/ (accessed on 5 July 2021).
  33. Tryggvason, G.; Scardovelli, R.; Zaleski, S. Fluid mechanics with interfaces. In Direct Numerical Simulations of Gas–Liquid Multiphase Flows; Cambridge University Press: Cambridge, UK, 2011; pp. 21–49. [Google Scholar]
  34. Sethian, J.A. Theory, algorithms, and applications of level set methods for propagating interfaces. Acta Numer. 1996, 5, 309–395. [Google Scholar] [CrossRef] [Green Version]
  35. Sethian, J.A.; Smereka, P. Level set methods for fluid interfaces. Annu. Rev. Fluid Mech. 2003, 35, 341–372. [Google Scholar] [CrossRef] [Green Version]
  36. Olsson, E.; Kreiss, G.; Zahedi, S. A conservative level set method for two phase flow II. J. Comput. Phys. 2007, 225, 785–807. [Google Scholar] [CrossRef]
  37. Walker, C.; Müller, B. A conservative level set method for sharp interface multiphase flow simulation. In Proceedings of the V European Conference on Computational Fluid Dynamics, ECCOMAS CFD 2010, Lisbon, Portugal, 14–17 June 2010; pp. 14–17. [Google Scholar]
  38. Gibou, F.; Chen, L.; Nguyen, D.; Banerjee, S. A level set based sharp interface method for the multiphase incompressible Navier-Stokes equations with phase change. J. Comput. Phys. 2007, 222, 536–555. [Google Scholar] [CrossRef]
  39. Fattori, M. A Level Set And Sharp Interface Approach for Simulating Incompressible Two-Phase Flow. Master’s Thesis, University of Waterloo, Waterloo, ON, Canada, 2014. [Google Scholar]
  40. Sharma, A. Level set method for computational multi-fluid dynamics: A review on developments, applications and analysis. Sadhana Acad. Proc. Eng. Sci. 2015, 40, 627–652. [Google Scholar] [CrossRef] [Green Version]
  41. Jasper, W.J.; Anand, N. A generalized variational approach for predicting contact angles of sessile nano-droplets on both flat and curved surfaces. J. Mol. Liq. 2019, 281, 196–203. [Google Scholar] [CrossRef]
  42. Lu, L.; Jiang, S.; Zhao, M.; Tang, G. Two-dimensional viscous numerical simulation of liquid sloshing in rectangular tank with/without baf fl es and comparison with potential fl ow solutions. Ocean Eng. 2015, 108, 662–677. [Google Scholar] [CrossRef]
  43. Battaglia, L.; Cruchaga, M.; Storti, M.; D’Elía, J.; Aedo, J.N.; Reinoso, R. Numerical modelling of 3D sloshing experiments in rectangular tanks. Appl. Math. Model. 2018, 59, 357–378. [Google Scholar] [CrossRef] [Green Version]
  44. Liu, Q.; Alazemi, S.F.; Daqaq, M.F.; Li, G. A ferrofluid based energy harvester: Computational modeling, analysis, and experimental validation. J. Magn. Magn. Mater. 2018, 449, 105–118. [Google Scholar] [CrossRef]
  45. FerroTec. EFH Series Educational Ferrofluid Type: EFH3; Ferrotec (USA) Corporation: Santa Clara, CA, USA, 2019. [Google Scholar]
Figure 1. Schematic representation of an energy harvesting device.
Figure 1. Schematic representation of an energy harvesting device.
Fluids 06 00295 g001
Figure 2. Classification of fluidic energy harvesting methods.
Figure 2. Classification of fluidic energy harvesting methods.
Fluids 06 00295 g002
Figure 3. Schematic representation of the proposed harvester configuration.
Figure 3. Schematic representation of the proposed harvester configuration.
Fluids 06 00295 g003
Figure 4. Harvester and the large domain surrounding it.
Figure 4. Harvester and the large domain surrounding it.
Fluids 06 00295 g004
Figure 5. Arrangement of tank with horizontally wound coil.
Figure 5. Arrangement of tank with horizontally wound coil.
Fluids 06 00295 g005
Figure 6. Arrangement of tank with vertically wound coil.
Figure 6. Arrangement of tank with vertically wound coil.
Fluids 06 00295 g006
Figure 7. Comparison between (a) 1H configuration of Liu [44] and (b) present work.
Figure 7. Comparison between (a) 1H configuration of Liu [44] and (b) present work.
Fluids 06 00295 g007
Figure 8. Impedance matching for horizontal coil of 1000 turns on L = 10 cm tank for frequency = 2.05 Hz.
Figure 8. Impedance matching for horizontal coil of 1000 turns on L = 10 cm tank for frequency = 2.05 Hz.
Fluids 06 00295 g008
Figure 9. Impedance matching for Vertical coil of 1000 turns on L = 10 cm tank for frequency = 2.05 Hz.
Figure 9. Impedance matching for Vertical coil of 1000 turns on L = 10 cm tank for frequency = 2.05 Hz.
Fluids 06 00295 g009
Figure 10. Surface evolution and magnetic flux density arrow surface for case 1, near resonance frequency for L = 10 cm tank.
Figure 10. Surface evolution and magnetic flux density arrow surface for case 1, near resonance frequency for L = 10 cm tank.
Fluids 06 00295 g010
Figure 11. RMS voltage plots for Cases 1–16 while considering horizontally wound coil for L = 10 cm tank.
Figure 11. RMS voltage plots for Cases 1–16 while considering horizontally wound coil for L = 10 cm tank.
Fluids 06 00295 g011
Figure 12. RMS voltage plots for Cases 1–16 while considering vertically wound coil for L = 10 cm tank.
Figure 12. RMS voltage plots for Cases 1–16 while considering vertically wound coil for L = 10 cm tank.
Fluids 06 00295 g012
Figure 13. Maximum power output from the horizontal coil for L = 10 cm tank.
Figure 13. Maximum power output from the horizontal coil for L = 10 cm tank.
Fluids 06 00295 g013
Figure 14. Maximum power output from the vertical coil for L = 10 cm tank.
Figure 14. Maximum power output from the vertical coil for L = 10 cm tank.
Fluids 06 00295 g014
Figure 15. Surface evolution and magnetic flux density arrow surface for case 1, near resonance frequency for L = 5 cm tank.
Figure 15. Surface evolution and magnetic flux density arrow surface for case 1, near resonance frequency for L = 5 cm tank.
Fluids 06 00295 g015
Figure 16. RMS voltage plots for Cases 1–16 while considering horizontally wound coil for L = 5 cm tank.
Figure 16. RMS voltage plots for Cases 1–16 while considering horizontally wound coil for L = 5 cm tank.
Fluids 06 00295 g016
Figure 17. RMS voltage plots for Cases 1–16 while considering vertically wound coil for L = 5 cm tank.
Figure 17. RMS voltage plots for Cases 1–16 while considering vertically wound coil for L = 5 cm tank.
Fluids 06 00295 g017
Figure 18. Maximum power output from horizontal winding for L = 5 cm tank.
Figure 18. Maximum power output from horizontal winding for L = 5 cm tank.
Fluids 06 00295 g018
Figure 19. Maximum power output from Vertical winding for L = 5 cm tank.
Figure 19. Maximum power output from Vertical winding for L = 5 cm tank.
Fluids 06 00295 g019
Table 1. Ferrofluid sloshing vibration energy harvesting configurations published in literature.
Table 1. Ferrofluid sloshing vibration energy harvesting configurations published in literature.
ReferenceResonant Frequency (Hz)Volume (cc)Power (µW)Power Density (µW/cc)Acceleration (m/s2)
[11]9 Hz44.230.850.023
[26]2.2 Hz104980,00076.2610
[27]2.1 Hz1000800.080.3
Table 2. Cases for Magnet orientation.
Table 2. Cases for Magnet orientation.
CaseMagnet 1Magnet 2Magnet 3Magnet 4
1South (k1 = −1)North (k2 = −1)South (k3 = −1)North (k4 = −1)
2South (k1 = −1)North (k2 = −1)South (k3 = −1)South (k4 = 1)
3South (k1 = −1)North (k2 = −1)North (k3 = 1)North (k4 = −1)
4South (k1 = −1)North (k2 = −1)North (k3 = 1)South (k4 = 1)
5South (k1 = −1)South (k2 = 1)South (k3 = −1)North (k4 = −1)
6South (k1 = −1)South (k2 = 1)South (k3 = −1)South (k4 = 1)
7South (k1 = −1)South (k2 = 1)North (k3 = 1)North (k4 = −1)
8South (k1 = −1)South (k2 = 1)North (k3 = 1)South (k4 = 1)
9North (k1 = 1)North (k2 = −1)South (k3 = −1)North (k4 = −1)
10North (k1 = 1)North (k2 = −1)South (k3 = −1)South (k4 = 1)
11North (k1 = 1)North (k2 = −1)North (k3 = 1)North (k4 = −1)
12North (k1 = 1)North (k2 = −1)North (k3 = 1)South (k4 = 1)
13North (k1 = 1)South (k2 = 1)South (k3 = −1)North (k4 = −1)
14North (k1 = 1)South (k2 = 1)South (k3 = −1)South (k4 = 1)
15North (k1 = 1)South (k2 = 1)North (k3 = 1)North (k4 = −1)
16North (k1 = 1)South (k2 = 1)North (k3 = 1)South (k4 = 1)
Table 3. Ferrofluid property data.
Table 3. Ferrofluid property data.
PropertyValue
Viscosity12 [mPa·s]
Density1420 [kg/m3]
Magnetic susceptibility3.52
Saturation magnetization65 [mT]
Table 4. Summary of properties of coils for L = 10 cm tank.
Table 4. Summary of properties of coils for L = 10 cm tank.
Coil NumberNumber of TurnsResistance of Horizontal Coil (ohm)Resistance of Vertical Coil (ohm)Inductance of Horizontal Coil (H)Inductance of Vertical Coil (H)
11003.4220.640.00830.0050
220027.385.130.0380.023
3400219.0441.070.170.10
4500427.8180.210.270.17
58001752.30328.550.750.47
610003422.47641.711.200.77
712005914.021108.881.771.14
8150011,550.832165.782.851.84
9200027,379.745133.705.253.42
Table 5. Summary of properties of coils for L = 5 cm tank.
Table 5. Summary of properties of coils for L = 5 cm tank.
Coil NumberNumber of TurnsResistance of Horizontal Coil (ohm)Resistance of Vertical Coil (ohm)Inductance of Horizontal Coil (H)Inductance of Vertical Coil (H)
11006.841.280.00420.0025
220054.7610.270.0190.012
3400438.0782.140.0840.053
4500855.62160.430.130.086
58003504.61657.110.370.24
610006844.931283.430.600.39
7120011828.052217.760.890.57
8150023101.664331.561.430.92
9200054759.4810267.402.631.71
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Anand, N.; Gould, R. Analysis of a Symmetrical Ferrofluid Sloshing Vibration Energy Harvester. Fluids 2021, 6, 295. https://doi.org/10.3390/fluids6080295

AMA Style

Anand N, Gould R. Analysis of a Symmetrical Ferrofluid Sloshing Vibration Energy Harvester. Fluids. 2021; 6(8):295. https://doi.org/10.3390/fluids6080295

Chicago/Turabian Style

Anand, Nadish, and Richard Gould. 2021. "Analysis of a Symmetrical Ferrofluid Sloshing Vibration Energy Harvester" Fluids 6, no. 8: 295. https://doi.org/10.3390/fluids6080295

APA Style

Anand, N., & Gould, R. (2021). Analysis of a Symmetrical Ferrofluid Sloshing Vibration Energy Harvester. Fluids, 6(8), 295. https://doi.org/10.3390/fluids6080295

Article Metrics

Back to TopTop