Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Dynamic Modeling and Vibration Characteristic Analysis of Fiber Woven Composite Shaft–Disk Rotor with Weight-Reducing Holes
Previous Article in Journal
Enhancing Ice Nucleation: The Role of Surface Roughness in Electrofreezing Using Laser Shock Processed Al6061 T6 Electrodes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Analysis and Closed-Form Solution for Overhead Line Magnetic Energy Harvester Behavior

School of Electrical and Computer Engineering, Ben-Gurion University of the Negev, P.O. Box 653, Beer-Sheva 84105, Israel
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(19), 9146; https://doi.org/10.3390/app14199146
Submission received: 18 August 2024 / Revised: 26 September 2024 / Accepted: 30 September 2024 / Published: 9 October 2024

Abstract

:

Featured Application

Magnetic energy harvesting of an overhead line.

Abstract

Recently, much attention has been given to the development of various energy harvesting technologies to power remote electronic sensors, data loggers, and communicators that can be installed on smart grid systems. Magnetic energy harvesting is, perhaps, the most straightforward way to capture a significant amount of power from a current-carrying overhead line. Since the harvester is expected to have a small size, the high currents of the distribution system easily saturate its magnetic core. As a result, the operation of the magnetic harvester is highly nonlinear and makes precise analytical modeling difficult. The operation of an overhead line magnetic energy harvester (OLMEH) generating significant DC power output into a constant voltage load was investigated in this paper. The analysis method was based on the Froelich equation to analytically model the nonlinearity of the core’s BH characteristic. The main findings of this piecewise nonlinear analysis include a closed-form solution that accounts for both the core and rectifiers’ nonlinearities and provides an accurate prediction of OLMEH transfer window length, output current, and harvested power. Continuous and discontinuous operational modes are identified and the mode transition boundary is obtained quantitatively. The theoretical investigation was concluded by comparison with a computer simulation and also verified by the experimental results of a laboratory prototype harvester. A good agreement was found.

1. Introduction

Advanced ultra-low power VLSI and wireless communications techniques spurred the development of various sensor networks and wearable electronics [1,2,3]. Such devices can be used to great advantage in various domestic, commercial, vehicular, military, and health applications.
A common concern of these technologies is their dependence on power sources. Batteries, biofuel cells, and supercapacitors have to be periodically replaced or recharged, posing certain application difficulties, whereas energy sources like solar cells, thermoelectric generators, piezoelectric and triboelectric generators, and various energy scavengers/harvesters can allow for a long service life while avoiding the need for maintenance and battery replacement. Energy harvesting is a particularly attractive approach to developing miniature self-powered systems designed to operate in remote locations or locations with restricted access [4].
The power distribution network crisscrossing the urban environment is a readily available, predictable, and reliable energy source. Therefore, recent studies have investigated the possibility of energy scavenging from power lines [5].
Thermal energy harvesting exploits the temperature difference developed by the current-carrying conductor and the ambient. DC electric current can be attained by applying thermo-electric generators clamped to the cable [6,7].
Energy can also be harvested from the electric field of the power line [8,9,10,11]. This suggests creating a capacitor-like structure by wrapping a conductive surface around the (isolated) line conductor and taking advantage of the displacement currents. Alternatively, the harvester can be designed to operate at a safe standoff distance from a high-voltage line. Since the line voltage is well stabilized, the electric field harvesters have the advantage of a nearly constant output current.
Magnetic energy harvesters (MEHs) rely on the magnetic field generated by a current-carrying conductor. Several types of MEHs were developed [4] that exploit different physical principles of interaction with the magnetic field such as variable reluctance [12,13], magnetostrictive [14,15], and ferrofluid [16,17,18] generators. The magnetic field generated by the line current can also be converted into the mechanical vibration of a miniature resonant cantilever. The output power can then be obtained by using a piezoelectric [19,20] or inductive [21] transducer mounted on the free end of the cantilever beam. Yet, the mentioned harvester types can develop only a tiny amount of power in the range of up to a few mW. This may be sufficient to operate miniature wireless sensors; however, in case a substantially higher power is needed, a different approach is needed. A current transformer (CT) type MEH is, perhaps, the most promising concept that can allow the harvesting of a significant amount of power from a power line. Two configurations are mainly considered in the current literature: the “stick on” or “freestanding” harvester, which is mounted in the close vicinity of a bus bar, and the “clamp on”-type harvester, which is configured to “embrace” the current-carrying conductor. Although the output power of MEHs varies with the line current amplitude, the average daily current profile can be measured and the system can be designed to its minimum expected value. Thus, MEHs can be considered as a quite dependable energy source.
Traditionally, CTs are used as sensors to accurately measure high AC currents. The most desirable feature of a CT is its high linearity in a wide range of line currents. Thus, the preferable operating conditions of a CT’s magnetic core are in the linear region of its BH characteristic. To help attain such a condition, the CT is usually operated under a nearly short-circuited secondary. The construction of an overhead line magnetic energy harvester (OLMEH) is quite similar to that of a CT; yet, since the task of an OLMEH is often charging a battery, the operating conditions of an OLMEH are quite different. Figure 1 shows an OLMEH as a part of a simple power processing scheme operating with a constant voltage load (CVL). The circuit employs a full-wave bridge rectifier to perform the AC to DC conversion. The rectifier’s action imposes a square wave voltage (whose amplitude equals the CVL voltage) across the OLMEH’s secondary winding. This may drive the magnetic core into the saturation region of the BH curve. While saturated, the OLMEH can provide no output current for a considerable portion of a line cycle. Thus, predicting the output power of OLMEHs is a challenging problem.
Furthermore, the increased power comes at the expense of increased volume and weight. To minimize the mentioned physical parameters while attaining the desired output power, OLMEHs’ core material, geometry, size, and number of turns have to be properly selected. Therefore, the aim of this article is to reexamine the workings of OLMEHs and obtain a complete analytical description of device operation, which can further serve as a theoretical foundation to develop a reliable set of engineering design rules. The applied methodology is as follows.

2. The Applied Methodology

This paper is concerned with the analysis of an overhead line energy harvester (OLMEH) illustrated in Figure 1. The OLMEH in Figure 1 is of a clamped-on CT type, the concept for which can be found in [22,23,24,25,26]. An analysis of the OLMEH’s behavior, while modeling the nonlinear BH characteristic using the atan(x) function, was described in [23]. The choice of such an approximation function results in analytical complications and necessitates the application of numerical solutions to attain the final results. Alternatively, ref. [26] opted to model the BH curve by a simple piecewise linear function for the sake of deriving crude analytical results. In the low flux density region, the linear approximation allows easy analytical description of the device; however, in the high flux density region, where the core approaches the saturation knee, the linear approach is inadequate.
One can argue that, in our age, a researcher assisted with modern computational tools can crack any nonlinear problem. To attain the solution, it is sufficient to just apply simulation or to properly introduce basic laws and system descriptions to the computer and then obtain the desired numerical results and plots. Yet, there is a clear and undisputable advantage to the old school approach of deriving an analytical solution to the problem. The latter can show the explicit effect each variable has on the end result and so provide a deeper understanding of the physical process. The analytical solutions can also facilitate the derivation of straightforward designing rules. Hence, an extra effort was put in to revisit the earlier findings, particularly of [23]. This paper aims to develop a complete analytical model that faithfully describes OLMEH operation. The proposed modeling approach is based on the Froelich equation, which represents the nonlinearity of the BH curve by a first-order rational polynomial function. The choice of the Froelich equation as the approximation function is imperative. The advantage of the undertaken approach over the numerical [23] is that the Froelich approximation allows deriving closed-form analytical solutions that consider the complex behavior of the magnetic core, yet provides meaningful analytical results suited for engineering applications. In summary, the contribution of this paper over the earlier counterpart is the analytical solutions for the: (a) power transfer window; (b) charging output power, and (c) prediction of the operational modes boundary.
The rest of this paper is organized as follows: Section 3 reviews the modeling of the BH curve by the Froelich equation. Section 4 presents basic assumptions, preliminary simulation waveforms, and OLMEH models. Section 5 presents the analysis of the discontinuous current mode. The continuous current mode of OLMEHs is described in Section 6. Section 7 discusses the boundary conditions for continuous to discontinuous current mode change. Analytical solutions for output current and power are derived in Section 8. Verification and experimental results are given in Section 9. Finally, conclusions are given.

3. Review of the Froelich Equation

The Froelich equation is a first-order rational polynomial function that uses only two parameters, a and b, and provides a wide-range approximation to the non-hysteretic B H characteristic of a magnetic core material:
B H = a H b + H = a H b + H       H > 0 a H b H       H < 0 ,
here, a is given in [T] units, whereas b stands in [At/m]. The inverse function, H(B), can be derived as:
H B = b B a B = b B a B       B > 0 b B a + B       B < 0
In practice, the parameters a and b can be established by curve-fitting to the experimental data. Measured BH characteristics of a silicon steel core sample (EILOR MAGNETIC CORES) vs. the fitted BH curve, drawn by using the extracted Froelich parameters, are shown in Figure 2. Here, a = 2.05 [T] and b = 109.4 [A/m].
It can be seen that Froelich equation parameter a represents the saturation flux density of the core, B s . Yet, this value is unattainable in practice since it suggests imposing an infinite magnetic field intensity on the magnetic material:
B s = lim H a H b + H = a .
For the sake of a brief comparison, one possible approach to attaining a linearized BH curve is suggested here. Other approaches are also possible. Consider a piecewise linear approximation of the BH characteristic that includes the (B′, H′) operating point where the flux density at B = B s / 2 = a / 2 . Substituting in (2) yields the magnetic field intensity H = b . Thus, the equivalent magnetic permeability within the linear segment of the core can be approximated as the ratio:
μ e q = B H = a 2 b ,
The considerations above also suggest that the saturation magnetic field intensity of the piecewise-linear model, H s , is:
H s = B s μ e q = 2 b
Simulation is a helpful tool that can facilitate a detailed study of OLMEHs. Although PSIM requires a tedious trial-and-error procedure to adjust the magnetic core parameters to fit the experimental BH curve, once completed, this proved to be a well-spent effort. As shown in Figure 3a, the comparison of the simulated BH curve (obtained by PSIM v. 9.1) stands in good agreement with the fitted Froelich approximation (1).
A triple comparison of the approximated BH curve is shown in Figure 3b. The Froelich approximation, here with the mentioned parameters, vs. the B = 2 π B s a t a n ( H β ) approximation used by [23], here Bs = 2.05 and β = 126, and vs. the piecewise-linear approximation used by [26], is calculated by (3) and (5). As expected, the comparison reveals a significant error of the piecewise-linear approach in the vicinity of the saturation knee, above B s / 2 , whereas the Froelich and atan(x) approximations stand in excellent agreement with each other. As is shown in this paper, the advantage of the Froelich approximation is that it lends itself to analytical treatise and results in an analytical solution, whereas using the atan(x) approximation necessitates the application of numerical analysis [23].

4. The Basic Assumptions and the Equivalent Circuit

As simple as the OLMEH’s power processing circuit in Figure 1 may look, the nonlinear properties of the magnetic core and the rectifiers make an exact analysis prohibitive. Therefore, the following simplifying assumptions were adopted:
(1)
Since clamping the OLMEH on the line does not change the line current, and since the output power drawn by the OLMEH is negligible in comparison to the power carried by the high voltage transmission line, the OLMEH can be considered as fed by a sinusoidal current source.
(2)
The hysteresis phenomenon is neglected and the nonlinearity of the core is modeled by the Froelich Equations (1) and (2).
(3)
A simplified cantilever model is adopted to represent the OLMEH magnetic circuit with a nonlinear magnetizing inductance, Lm, placed on the secondary winding, whereas the leakage inductance and the winding resistance are considered negligible.
(4)
Also, ideal rectifiers with no voltage drop and no reverse recovery are assumed.
The assumptions above lead to the OLMEH model illustrated in Figure 4a. The model was simulated in PSIM v 9.1. The key simulated waveforms shown in Figure 5 reveal the intricate details of OLMEH operation.

5. Analysis of the Discontinuous Current Mode

An inspection of the OLMEH’s key simulated waveforms, presented in Figure 5, reveals that the rectifier’s input current, iin, vanishes for a substantial time interval. Henceforth, this operational mode is dubbed the discontinuous current mode (DCM).
While in DCM, each line half-cycle is comprised of two distinct time subintervals that emerge according to the conduction state of the rectifier; see Figure 5. These are designated as State 1 and State 2. The equivalent circuits of State 1 and State 2 are shown in Figure 4b and Figure 4c, respectively.
State 1: Commences at the instance the rectifiers are cut off at zero current; thus, the CVL receives no current from OLMEH, ib = |iin| = 0. This occurs because the magnetizing current steers all the secondary current away from the rectifier, i m t = i 2 t ; see Figure 5. Hence, for the low power line angle, ω t , values, see Figure 5, the expression for the magnetizing current can be approximated by:
i m t = i 2 t = I 2 m s i n ω t I 2 m ω t ,
where I1m is the peak line current, I 2 m = I 1 m / N is the peak secondary current, and N is the secondary number of turns.
During State 1, the magnetic field intensity in the core is dictated by the line current constrain:
H t I 1 m ω t l c = N I 2 m ω t l c < 0
where lc is the magnetic path length of the OLMEH’s core.
Therefore, according to (7) the line current generates the magnetic field, H, that, according to (1), forces the magnetic flux density, B, to decrease. Hence, the core is taken out of deep saturation towards a shallower saturation; see Figure 5. Flux density variation induces a voltage in the OLMEH’s secondary winding, which appears across the rectifier’s AC terminals as:
v i n t = N A c d B d t = N A c d d t a H b H = a b l c A c N 2 ω I 2 m b l c N I 2 m ω t 2
Here, Ac is the cross-section of the harvester’s core. Note, that for negative, H < 0, values, the Froelich Equation (1) with the (bH) term in the denominator is used in (8).
Equation (8) describes the concave segment seen in the Vin voltage waveform in Figure 5.
State 1 terminates at the instant, t1, when the voltage developed by the OLMEH’s winding equals the CVL voltage, thus turning the rectifiers on. The time instant, t 1 , can be found by applying the condition, v i n t 1 = V b to (8), whence
t 1 = b l c ω N I 2 m a b A c l c ω I 2 m V b
Note that (9) returns a negative result, t 1 < 0 , because State 1 terminates before the secondary current zero crossing, designated as time reference t = 0; see Figure 5.
Combining (1), (7), and (9) yields the value of the flux density, B 1 = B t 1 , at which State 1 is terminated:
B 1 = B t 1 a N I 2 m ω t 1 l c b N I 2 m ω t 1 l c
State 2 commences at t = t 1 + . Here, the rectifiers’ conduction imposes a constant CVL voltage, Vb, across the harvester’s secondary. As a result, the magnetic flux density in the core starts increasing linearly with time; see Figure 5:
B t = 1 N A c t 1 t V b d t + B 1 = V b N A c t t 1 + B 1
The flux density crosses zero at the instant:
t Z = t 1 N A c V b B 1
It is worthwhile noting that during State 2 the winding voltage constrains the flux density; hence, the magnetic field intensity follows according to (2). By Ampere’s law, the rising magnetic field intensity generates the magnetizing current:
i m t t > t Z = l c N H t = l c N b V b N A c t t Z a V b N A c t t Z
State 2 is terminated at t 2 , see Figure 5, when the rectifier’s input current is reduced to zero:
i i n t 2 = i 2 t 2 i m t 2 = 0
The Condition (14) can be rewritten using (13):
i m t 2 = l c N b V b N A c t 2 t Z a V b N A c t 2 t Z = I 2 m s i n ω t 2
The solution for the State 2 termination instant, t 2 , can be found by solving (15) numerically, whereas the analytical solution of (15) can be obtained by approximating the sine function on the right-hand side. Approximation is expected to represent a sinusoidal segment in the range of interest which, according to Figure 5, is located on the falling slope of the sine function. The first that comes to mind is to apply a truncated Taylor series. However, using the high-order polynomial makes the analytic solution of (15) an impossible task. Therefore, several options may be considered to approximate the sine function in the range 2 3 π < ω t < 5 6 π . The simple linear approximation of the type π ω t cannot provide sufficient accuracy within the desired range. A downgoing parabola of the type 4 π 2 ω t π ω t is another candidate that can provide acceptable accuracy in the full range of 0 < ω t < π . When substituted in (15), this will result in a cubic equation, which can be solved, however, with some difficulties. As a compromise between complexity, accuracy, and range, an approximation based on a “mirrored Froelich equation” is suggested here to approximate the falling segment of the sine function as:
s i n ω t K 1 π ω t K 2 ω t
with K1 = 3.232 and K2 = 6.003, (16) agrees with the s i n ω t function at 120°, 150°, and 180°, thus providing good accuracy within the range of interest. A comparison plot of the mentioned approximation functions is shown in Figure 6.
Substitution of (16) into (15) yields:
l c N b V b N A c t 2 t Z a V b N A c t 2 t Z = I 2 m K 1 π ω t 2 K 2 ω t 2
As shown in Appendix A, the solution for (17), i.e., State 2 termination instant, t 2 , is quite tedious. The result can be written using the normalized solution, t 2 N , (A21), as:
t 2 = T 2 t 2 N

6. Continuous–Discontinuous Mode Boundary

The discontinuous current mode (DCM) regime described above is the result of a relatively high volt per turn applied across the harvester’s winding that brings the core to saturation before the end of a half-cycle. In DCM, the instant the rapidly rising magnetizing current intercepts the secondary current the rectifiers are cut off at zero current at t 2 < T / 2 ; see Figure 5. The rectifiers resume conduction at the instant the voltage induced in the winding (due to the varying magnetic field generated by the line current) equals the CVL voltage; see (8). However, for lower values of the applied voltage, the build-up of the magnetizing current is slower and the interception point occurs at the instant, t 2 , at which the induced voltage is already of a sufficient magnitude to initiate the immediate rectifier’s conduction (yet, in the reversed polarity). Thus, the dead time in the rectifiers’ input current, iin, see Figure 5, vanishes. As the result, for lower CVL voltages the rectifier can conduct for the entire half-cycle. Therefore, this operational regime is dubbed the continuous current mode (CCM). Since in CCM the rectifiers’ conduction lasts for the entire half-cycle and the operation is symmetrical, the magnetic flux density in the core becomes a true triangular waveform. Simulated waveforms of an OLMEH at the CCM-DCM boundary are shown in Figure 7c. Here, the applied CVL voltage, V b * , brings the positive peak value of the magnetic flux density infinitesimally close to B 1 (recall that B1 is negative). Applying the CCM condition t 2 t 1 = T / 2 to (11) yields the DCM-CCM boundary condition:
V b * N A c T 2 + 2 B 1 = 0
Substitution of (9) and (10) into (19) and further manipulation gives:
V b * = 2 ω a N A c π 2 π + b l c I 1 m 2 π b l c I 1 m b l c I 1 m 2
However, under a constant voltage load condition the CCM-DCM mode transition is determined by the line current, which can be derived from (20):
I 1 m * = 4 ω a b N A c l c V b π V b 2 ω a N A c 2
The mode changes under varying CVL voltage and line current conditions are illustrated in Figure 7. In the given simulation example, the calculated DCM-CCM boundary occurs at V b * = 23.16   V for I 1 * = 100   A rms; see Figure 7c. DCM prevails in the low-line current range, I 1 < I 1 * , as in Figure 7a or for a high voltage range V b > V b * as in Figure 7b, whereas CCM prevails in the high line current range, I 1 > I 1 * , as in Figure 7d or for a low voltage range V b < V b * as in Figure 7e.

7. Output Current and Power Considerations

In DCM, CVL charging commences at t 1 and terminates at t 2 . Hence, the time interval t 1 , t 2 is defined as the power transfer window. The charging current at the rectifier’s output can be found as the per-half-cycle average of the difference between the secondary current, i 2 t , and the magnetizing currents, i m t , drawn by the OLMEH’s magnetizing inductance:
I b = 1 T / 2 t 1 t 2 i 2 t i m t d t
Calculation of the average secondary current component available for charging, which is the first term in (22), is straightforward:
I 2 c h = 1 T / 2 t 1 t 2 i 2 t d t = 1 T / 2 t 1 t 2 I 2 m sin ω t d t = I 2 m π cos ω t 1 cos ω t 2
However, finding the per-half-cycle average of the magnetizing current, which is the second term in (22), requires additional consideration. Since the flux density changes sign at t = t Z , see Figure 5, two subintervals t 1 , t Z and t Z , t 2 each with appropriate models should be considered to produce the correct result. Therefore:
I m a g = 1 T / 2 t 1 t 2 i m t d t = 1 T / 2 t 1 t Z i m t d t + 1 T / 2 t Z t 2 i m t d t = I m a g 1 + I m a g 2
The average of the negative (recycled) term of the magnetizing current, I m a g 1 , during the time interval t 1 , t Z , see Figure 5, should be calculated according to the negative flux density model, see (2), with B 1 as the initial condition:
I m a g 1 = 1 T / 2 t 1 t Z i m t d t = 1 T / 2 t 1 t Z l c N b V b N A c t t 1 + B 1 a + V b N A c t t 1 + B 1 d t
The (recharging) magnetizing current term, I m a g 2 , should be calculated within the time interval t Z ,   t 2 , see Figure 5, according to the positive flux density model, (2), with B = 0 as the initial condition:
I m a g 2 = 1 T / 2 t Z t 2 i m t d t = 1 T / 2 t Z t 2 l c N b V b N A c t t Z a V b N A c t t Z d t
Skipping the details of the derivation of (25) and (26), which is quite tedious, substitution into (24) yields the per half-cycle average magnetizing current:
I m a g = ω π b l c N 2 t Z t 2 t 1 ω π b l c N   a N A c V b l n 1 V b t 2 t Z a N A c 1 + V b t Z t 1 a N A c 1 + B 1 / a
Applying (23) and (27), the OLMEH’s per half-cycle average output power in the DCM mode can be found as:
P o = V b I 2 c h I m a g
While in CCM, the magnetizing current is symmetric and has a zero average during half a cycle; see Figure 5. This implies that the energy stored in the core is recycled to the CVL within the half-cycle (note that in DCM part of the energy is recycled to the system during State 1, not charging the CVL). Therefore, the average output power in the CCM mode is simply:
P o = V b I 2 c h 2 2 π N V b I 1 r m s
The output power in CCM mode (29) is a linearly increasing function of the CVL voltage, Vb, whereas in the DCM mode, the output power (28) appears as a bell-shaped function with a distinct extremum. This means that an increase of the CVL voltage does not necessarily bring an increase in the output power but depends on the operating point. Hence, the OLMEH designer should properly select the OLMEH parameters while aiming at the maximum power point. OLMEH design issues are out of scope of this paper and will be examined elsewhere.

8. Verification and Experimental Results

Verification of the theoretical predictions was conducted first by simulation. The theoretical and simulated results for key parameters were found to be in excellent agreement; see Table 1 and Table 2.
An experimental laboratory prototype OLMEH shown in Figure 1 was constructed by stacking together three pairs of 10H10 C-cores made of silicon steel (EILOR MAGNETIC CORES). The magnetic path length was lc = 120 mm, and the total core cross-section was Ac = 1800 mm2. The total number of secondary turns was N = 40. The SB256 bridge rectifier was used.
To test the performance of the prototype OLMEH, an experimental test bed was constructed, see Figure 8, using two back-to-back connected isolation transformers capable of sustaining a low-side current of up to 200 A rms. The OLMEH line current was regulated using a VARIAC and a high-power resistive load. In Figure 8, the OLMEH is illustrated successfully charging a battery.
Magnetic properties of the prototype OLMEH were tested and Froelich coefficients were extracted from the experimental data as mentioned.
The prototype OLMEH was tested at several line currents and CVL voltage levels. During the experiments, electronic load was used to emulate the voltage sink load. Key waveforms of the OLMEH are shown in Figure 9. Here, the line current was I1 = 100 A rms. The cusp portion of the rectifiers’ input voltage waveform, Vin, manifests State 1, whereas the flat portion reveals State 2. Accordingly, the CVL current is zero during State 1, whereas CVL is provided with a pulse of current during State 2. This stands in agreement with the theoretical prediction above.
For Vb = 25 V, the measured output power was 49.4 W, see Figure 9a, whereas for Vb = 35 V the output power was 44.6 W, see Figure 9b. Also note that for Vb = 35 V the duration of State 1 is relatively long and is seen in Figure 9b. However, at the reduced CVL voltage, Vb = 25 V, the duration of State 1 is shrinking; see Figure 9a. This shows that OLMEH is approaching the CCM-DCM boundary.
The plot of the OLMEH’s output power vs. the CVL voltage at a line current of 100 A rms and 180 A rms is shown in Figure 10.
A comparison of the calculated vs. the measured output power as a function of the line current at fixed CVL voltage Vb = 30 V is presented in Figure 11. An excellent match was obtained.
A comparison of the calculated and the experimental results for line current I1 = 100 A rms is shown in Figure 12. CCM mode is observed at Vb < 20 V, whereas DCM settles in at a higher CVL voltage level.

9. Conclusions

This paper offers a model, analysis, and analytical solutions for a clamped-type overhead line energy harvester with a full-wave rectifier with natural commutation loaded by a constant voltage load.
This paper also reviews the Froelich equation and adopts it as a main tool to investigate harvester operation. This paper further suggests a model and qualitative description of the energy conversion process, and presents the simulation results followed by a full nonlinear analysis supported by experimental results. The contribution of this paper is the complete analytical solutions for charging and idle time intervals, the solution for the average output current, and the solution for the average output power.
The proposed theory was verified by comparison to simulation and experimental data collected from a laboratory prototype. An excellent match was observed. The prototype harvester was shown to provide a substantial charging power of up to 50 W in 25 V CVL at 100 A rms line current.
Compared to the earlier models [23,26], the choice of the Froelich equation as the approximation function to describe the core BH characteristic is critical, and proved to be well-suited to facilitate the derivation of complete analytical solutions to the key parameters of OLMEHs. The proposed approach offers both meaningful engineering insight and an accurate analytical description of the OLMEH’s performance. The undertaken approach also allowed finding the CCM-DCM regime boundary that has not yet been discussed in the literature.
Yet, the solutions are still an approximation. Firstly, because the Froelich equation does not consider the hysteresis feature and, secondly, since (16) is valid in a limited range of the power line angle the analysis result can be applied in cases where the width of the power transfer window exceeds 100°.
This experimental work revealed that handling the OLMEH prototype affects its magnetic properties. This is mainly attributed to the variations of the technological air gap remaining between the two C-core halves of a split magnetic core that depends on the mounting force, surface imperfections, etc. Therefore, in the authors’ opinion, the derivation of a more precise and, thus, more complex theoretical model (that includes the hysteresis effect) is of questionable practical value due to the mentioned, and unavoidable in practice, core parameter variations.

Author Contributions

Conceptualization, A.K.; methodology, A.A.; software, A.A.; validation, A.A. and M.S.; formal analysis, A.A.; investigation, A.A.; resources, A.K.; data curation, M.S.; writing—original draft preparation, A.A.; writing—review and editing, A.A.; visualization, A.A. and M.S.; supervision, A.K.; project administration, A.K.; funding acquisition, A.K. 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

The data is presented in this manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

Throughout this paper, the following definitions of key parameters were adopted:
Time (continuous):t
State 1 termination instantt1
State 2 termination instantt2
Flux density zero crossing instanttZ
Power line angular frequency ω
Flux density: B t
Saturation flux densityBs
Field intensity: H t
Saturation field intensityHs
Froelich coefficients:a, b
Number of turnsN
Magnetic path lengthlc
Magnetic core areaAc
Line (primary) current instantaneousi1(t)
Line (primary) peak currentI1m
Secondary current instantaneousi2(t)
Secondary peak currentI2m
Magnetizing current, (instantaneous)im(t)
Magnetizing current, per-half-cycle averageImag
Rectifiers’ input (AC side) voltage v i n t
Rectifiers’ input (AC side) current i i n t
CVL voltageVb
CVL current i b t
Average output powerPo
Definitions of additional auxiliary parameters can be found in the text.

Appendix A. Calculation of the Power Transfer Window

In the following, to keep the solutions manageable, the variables are normalized with respect to the following base quantities:
Base time
T b a s e = T 2 = π ω
here, as usual, T is the line period and ω is the line angular frequency.
Base flux density:
B b a s e = a
Base field intensity:
H b a s e = b
Base current, referred to secondary:
I b a s e = H b a s e l c N = b l c N
Base voltage
V b a s e = N A c B b a s e T b a s e = ω a N A c π
Normalized time
t N = t T b a s e = ω t π
Normalized flux density:
B N = B B b a s e = B a
Normalized field intensity:
H N = H H b a s e = H b
Normalized secondary current
I N = I 2 m I b a s e = N I 2 m b l c = I 1 m b l c
Normalized CVL voltage
V N = V b V b a s e = π V b ω a N A c
Applying the condition, v i n t 1 = V b to (8) and normalizing:
V N = π I N 1 π I N t N 2
Whence, the normalized time instant, t1N, is:
t 1 N = 1 π I N 1 π I N V N
The normalized flux density at the rectifier’s cut-in instant is:
B 1 N = B t 1 a N I 2 m ω t 1 b l c 1 N I 2 m ω t 1 b l c = π I N ω t 1 π 1 π I N ω t 1 π = π I N t 1 N 1 π I N t 1 N
The normalized instant when flux density crosses zero is:
t Z N = t 1 N B 1 N V N
Equation (17) can be first transformed into
V b t 2 t Z a N A c V b t 2 t Z = N I 2 m b l c K 1 π ω t 2 K 2 ω t 2
then normalized
V N t 2 N t Z N 1 V N t 2 N t Z N = I N K 1 π 1 t 2 N K 2 π t 2 N
and after additional manipulation put into the form of a quadratic equation:
α t 2 N 2 + β t 2 N + γ = 0
where
α = π K 1 V N I N + π V N
β = π K 1 I N + π K 1 1 + t Z N V N I N + π t Z N + K 2 V N
γ = π K 1 I N + π K 1 t Z N V N I N + K 2 t Z N V N
Only a positive solution for t2N has a physical meaning:
t 2 N = β + β 2 4 α γ 2 α

References

  1. Cenkeramaddi, L.R.; Bhatia, J.; Jha, A.; Vishkarma, S.K.; Soumya, J. A Survey on Sensors for Autonomous Systems. In Proceedings of the 2020 15th IEEE Conference on Industrial Electronics and Applications (ICIEA), Kristiansand, Norway, 9–13 November 2020; pp. 1182–1187. [Google Scholar] [CrossRef]
  2. Sharma, S.; Kaur, A. Survey on Wireless Sensor Network, Its Applications and Issues. J. Phys. Conf. Ser. 2021, 1969, 012042. [Google Scholar] [CrossRef]
  3. De-la-Fuente-Robles, Y.-M.; Ricoy-Cano, A.-J.; Albín-Rodríguez, A.-P.; López-Ruiz, J.L.; Espinilla-Estévez, M. Past, Present and Future of Research on Wearable Technologies for Healthcare: A Bibliometric Analysis Using Scopus. Sensors 2022, 22, 8599. [Google Scholar] [CrossRef] [PubMed]
  4. Enayati, J.; Asef, P. Review and Analysis of Magnetic Energy Harvesters: A Case Study for Vehicular Applications. IEEE Access 2022, 10, 79444–79457. [Google Scholar] [CrossRef]
  5. Riba, J.-R.; Moreno-Eguilaz, M.; Bogarra, S. Energy Harvesting Methods for Transmission Lines: A Comprehensive Review. Appl. Sci. 2022, 12, 10699. [Google Scholar] [CrossRef]
  6. Guo, F.; Hayat, H.; Wang, J. Energy harvesting devices for high voltage transmission line monitoring. In Proceedings of the 2011 IEEE Power and Energy Society General Meeting, Detroit, MI, USA, 24–28 July 2011; pp. 1–8. [Google Scholar]
  7. Becker, T.; Kluge, M.; Schalk, J.; Otterpohl, T.; Hilleringmann, U. Power Management for Thermal Energy Harvesting in Aircrafts. In Proceedings of the IEEE Sensors Conference, Lecce, Italy, 26–29 October 2008; pp. 681–684. [Google Scholar]
  8. Kang, S.; Yang, S.; Kim, H. Non-intrusive voltage measurement of ac power lines for smart grid system based on electric field energy harvesting. Electron. Lett. 2017, 53, 181–183. [Google Scholar] [CrossRef]
  9. Moser, M.J.; Bretterklieber, T.; Zangl, H.; Brasseur, G. Strong and Weak Electric Field Interfering: Capacitive Icing Detection and Capacitive Energy Harvesting on a 220-kV High-Voltage Overhead Power Line. IEEE Trans. Ind. Electron. 2011, 58, 2597–2604. [Google Scholar] [CrossRef]
  10. Li, M.; Luo, A.; Luo, W.; Wang, F. Recent Progress on Mechanical Optimization of MEMS Electret-Based Electrostatic Vibration Energy Harvesters. J. Microelectromechanical Syst. 2022, 31, 726–740. [Google Scholar] [CrossRef]
  11. Khan, F.U.; Qadir, M.U. State-of-the-art in vibration-based electrostatic energy harvesting. J. Micromech. Microeng. 2016, 26, 103001. [Google Scholar] [CrossRef]
  12. Olszewski, O.Z.; Houlihan, R.; Blake, A.; Mathewson, A.; Jackson, N. Evaluation of vibrational PiezoMEMS harvester that scavenges energy from a magnetic field surrounding an AC current-carrying wire. J. Microelectromech. Syst. 2017, 26, 1298–1305. [Google Scholar] [CrossRef]
  13. Zorlu, Ö.; Topal, E.T.; Külah, H. A vibration-based electromagnetic energy harvester using mechanical frequency up-conversion method. IEEE Sens. J. 2011, 11, 481–488. [Google Scholar] [CrossRef]
  14. Zhangxian Deng and Marcelo J Dapino “Review of magnetostrictive vibration energy harvesters. Smart Mater. Struct. 2017, 26, 103001. [CrossRef]
  15. Liu, H.; Zhao, X.; Liu, H.; Yang, J. Magnetostrictive biomechanical energy harvester with a hybrid force amplifier. Int. J. Mech. Sci. 2022, 233, 1. [Google Scholar] [CrossRef]
  16. 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]
  17. 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]
  18. 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, 1–6. [Google Scholar] [CrossRef]
  19. Kanno, I. Piezoelectric MEMS for energy harvesting. J. Phys. Conf. Ser. 2015, 660, 012001. [Google Scholar] [CrossRef]
  20. Wang, Q.; Kim, K.-B.; Woo, S.; Sung, T. Magnetic Field Energy Harvesting with a Lead-Free Piezoelectric High Energy Conversion Material. Energies 2021, 14, 1346. [Google Scholar] [CrossRef]
  21. Xiao, H.; Peng, H.; Yuan, J. A Coil Connection Switching Strategy for Maximum Power Delivery in Electromagnetic Vibration Energy Harvesting System. In Proceedings of the 2021 IEEE Energy Conversion Congress and Exposition (ECCE), Vancouver, BC, Canada, 10–14 October 2021; pp. 5867–5872. [Google Scholar]
  22. Cheang, J.; Cheng, W.; Gavrilov, D.; Schiller, B.; Smagin, V.; Gouzman, M. High efficiency powering system for wireless sensor for ac monitoring in smart grid applications. In Proceedings of the 2014 11th International Conference & Expo on Emerging Technologies for a Smarter World (CEWIT), Melville, NY, USA, 29–30 October 2014; pp. 1–5. [Google Scholar]
  23. Moon, J.; Leeb, S.B. Analysis model for magnetic energy harvesters. IEEE Trans. Power Electron. 2015, 30, 4302–4311. [Google Scholar] [CrossRef]
  24. Vos, M.J. A magnetic core permeance model for inductive power harvesting. IEEE Trans. Power Electron. 2020, 35, 3627–3635. [Google Scholar] [CrossRef]
  25. Liu, Z.; Li, Y.; Yang, H.; Na, D.; He, Z. An accurate model of magnetic energy harvester in the saturated region for harvesting maximum power: Analysis, design and experimental verification. IEEE Trans. Ind. Electron. 2023, 70, 276–285. [Google Scholar] [CrossRef]
  26. Li, Y.; Duan, N.; Liu, Z.; Hu, J.; He, Z. Impedance-Matching-Based Maximum Power Tracking for Magnetic Field Energy Harvesters Using Active Rectifiers. IEEE Trans. Ind. Electron. 2023, 70, 10730–10739. [Google Scholar] [CrossRef]
Figure 1. Overhead line energy harvester under constant voltage load.
Figure 1. Overhead line energy harvester under constant voltage load.
Applsci 14 09146 g001
Figure 2. Comparison of the measured BH curve of the silicon steel core sample (EILOR MAGNETIC CORES) to its Froelich approximation.
Figure 2. Comparison of the measured BH curve of the silicon steel core sample (EILOR MAGNETIC CORES) to its Froelich approximation.
Applsci 14 09146 g002
Figure 3. PSIM—generated comparison plot of the simulated BH curve vs. the Froelich approximation (1) (a); MATHCAD—generated Froelich Equation (1) vs. atan(*) approximation, vs. the piecewise-linear approximation of the BH curve (b).
Figure 3. PSIM—generated comparison plot of the simulated BH curve vs. the Froelich approximation (1) (a); MATHCAD—generated Froelich Equation (1) vs. atan(*) approximation, vs. the piecewise-linear approximation of the BH curve (b).
Applsci 14 09146 g003
Figure 4. Simplified OLMEH model (a) and its DCM equivalent circuits during the positive half-cycle with: rectifier OFF (State 1) (b); rectifier ON (State 2) (c).
Figure 4. Simplified OLMEH model (a) and its DCM equivalent circuits during the positive half-cycle with: rectifier OFF (State 1) (b); rectifier ON (State 2) (c).
Applsci 14 09146 g004
Figure 5. OLMEH’s key simulated waveforms in the discontinuous mode.
Figure 5. OLMEH’s key simulated waveforms in the discontinuous mode.
Applsci 14 09146 g005
Figure 6. Comparison plot of the sin(x) function vs. the linear approximation, vs. the parabolic approximation, and vs. the “mirrored Froelich equation” (15) (for K1 = 3.232, K2 = 2.861).
Figure 6. Comparison plot of the sin(x) function vs. the linear approximation, vs. the parabolic approximation, and vs. the “mirrored Froelich equation” (15) (for K1 = 3.232, K2 = 2.861).
Applsci 14 09146 g006
Figure 7. OLMEH’s simulated waveforms in: DCM due to low current I1 = 50 A rms, Vb = 23 V (a); DCM due to high voltage I1 = 100 A rms, Vb = 25 V (b); DCM-CCM boundary I1 = 100 A rms, Vb = 22 V (c); CCM due to high current I1 = 125 A rms, Vb = 23 V (d); CCM due to low voltage I1 = 100 A rms, Vb = 22 V (e).
Figure 7. OLMEH’s simulated waveforms in: DCM due to low current I1 = 50 A rms, Vb = 23 V (a); DCM due to high voltage I1 = 100 A rms, Vb = 25 V (b); DCM-CCM boundary I1 = 100 A rms, Vb = 22 V (c); CCM due to high current I1 = 125 A rms, Vb = 23 V (d); CCM due to low voltage I1 = 100 A rms, Vb = 22 V (e).
Applsci 14 09146 g007
Figure 8. View of the experimental prototype OLMEH (a) and its test bench (b).
Figure 8. View of the experimental prototype OLMEH (a) and its test bench (b).
Applsci 14 09146 g008
Figure 9. Typical waveforms of the experimental overhead line energy harvester at: I1 = 100 A rms and Vb = 25 V (a); I1 = 100 A and Vb = 35 V (b); I1 = 75 A rms and Vb = 23 V (c); I1 = 150 A rms and Vb = 23 V (d). Vert. scale: I1—200 A/div; Vin—20 V/div; V1—2 V/div; Iin—5 A/div; hor. scale 5 ms/div.
Figure 9. Typical waveforms of the experimental overhead line energy harvester at: I1 = 100 A rms and Vb = 25 V (a); I1 = 100 A and Vb = 35 V (b); I1 = 75 A rms and Vb = 23 V (c); I1 = 150 A rms and Vb = 23 V (d). Vert. scale: I1—200 A/div; Vin—20 V/div; V1—2 V/div; Iin—5 A/div; hor. scale 5 ms/div.
Applsci 14 09146 g009
Figure 10. Measured OLMEH output power vs. the CVL voltage at a line current of 100 A rms and 180 A rms.
Figure 10. Measured OLMEH output power vs. the CVL voltage at a line current of 100 A rms and 180 A rms.
Applsci 14 09146 g010
Figure 11. Comparison of the calculated vs. the measured output power as a function of the line current at fixed CVL voltage Vb = 30 V.
Figure 11. Comparison of the calculated vs. the measured output power as a function of the line current at fixed CVL voltage Vb = 30 V.
Applsci 14 09146 g011
Figure 12. Comparison of the calculated power, Pocalc, vs. the measured power, Pomeas, output power (DCM) as a function of the CVL voltage at a fixed line current (100 A rms).
Figure 12. Comparison of the calculated power, Pocalc, vs. the measured power, Pomeas, output power (DCM) as a function of the CVL voltage at a fixed line current (100 A rms).
Applsci 14 09146 g012
Table 1. Comparison of the simulated vs. the calculated key parameters at different peak line current levels and at a load voltage of Vb = 24 V.
Table 1. Comparison of the simulated vs. the calculated key parameters at different peak line current levels and at a load voltage of Vb = 24 V.
I1pk
[A]
t1Sim
[mS]
t1Calc
[mS]
t2Sim
[mS]
t2Calc
[mS]
B1Sim
[T]
B1Calc
[T]
POSim
[W]
POCalc
[W]
50−1.35−1.437.057.07−1.25−1.2912.8713.08
100−1.15−1.198.238.27−1.51−1.5233.0533.21
200−0.92−0.939.059.11−1.66−1.6773.2573.34
Table 2. Comparison of the simulated vs. the calculated borderline voltage Vb* for different turns number N, at line current I1 = 100 A rms.
Table 2. Comparison of the simulated vs. the calculated borderline voltage Vb* for different turns number N, at line current I1 = 100 A rms.
N40506070
CalcSimCalcSimCalcSimCalcSim
23.1623.2528.9529.534.7535.240.5441.14
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abramovitz, A.; Shwartsas, M.; Kuperman, A. Nonlinear Analysis and Closed-Form Solution for Overhead Line Magnetic Energy Harvester Behavior. Appl. Sci. 2024, 14, 9146. https://doi.org/10.3390/app14199146

AMA Style

Abramovitz A, Shwartsas M, Kuperman A. Nonlinear Analysis and Closed-Form Solution for Overhead Line Magnetic Energy Harvester Behavior. Applied Sciences. 2024; 14(19):9146. https://doi.org/10.3390/app14199146

Chicago/Turabian Style

Abramovitz, Alexander, Moshe Shwartsas, and Alon Kuperman. 2024. "Nonlinear Analysis and Closed-Form Solution for Overhead Line Magnetic Energy Harvester Behavior" Applied Sciences 14, no. 19: 9146. https://doi.org/10.3390/app14199146

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

Article Metrics

Back to TopTop