Next Article in Journal
Advanced AI Hardware Designs Based on FPGAs
Previous Article in Journal
Speeding Up Velocity Consensus Control with Small World Communication Topology for Unmanned Aerial Vehicle Swarms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Method of Coupling Coefficient Estimation Based on the Bifurcation Phenomena in Inductive Power Transfer

1
Department of Electric Drives and Traction, Czech Technical University in Prague, 166 27 Prague, Czech Republic
2
Department of Electrical Engineering and Renewable Energy, Oregon Institute of Technology, Wilsonville, OR 97070, USA
*
Author to whom correspondence should be addressed.
Electronics 2021, 10(20), 2548; https://doi.org/10.3390/electronics10202548
Submission received: 20 September 2021 / Revised: 14 October 2021 / Accepted: 15 October 2021 / Published: 18 October 2021
(This article belongs to the Section Power Electronics)

Abstract

:
Inductive power transfer (IPT) applications, such as stationary charging of electric vehicles (EVs), at least moderate coupling between the coils to achieve high efficiency, but the coefficient k typically varies between of 0.1 to 0.4, depending on the displacement of the coils according to SAE J2954. Thus, accurate and reliable methods for estimation of k are required for positioning of the EV to achieve optimal alignment with the charging pad. Additionally, in IPT, numerous control strategies are available for regulating output power and optimizing system efficiency that require an accurate estimate of the mutual inductance or k . However, existing estimation methods tend to require detailed a-priori information of a large number of circuit parameters, or they need measurement of currents or voltages in both primary and secondary sides. This paper presents a preliminary evaluation of a novel, primary-side method to estimate k , which is based solely on the frequency response of the input phase while operating the system in bifurcation. The method does not require any additional measurements of the system parameters. The theoretical background of the method is presented together with the description of the measurement procedure. The method is experimentally verified and compared with two currently used estimation methods. According to the presented experimental evaluation, the proposed method estimates k with an error of 3.62% with respect to the reference over the evaluated range of 0.08 to 0.36. In addition, we demonstrate that the presented method is resilient to detuning.

1. Introduction

Inductive power transfer (IPT) is a developing technology for simplifying and streamlining the electric vehicle (EV) charging experience [1,2,3]. Based on the movement of the vehicle, this application can be divided in dynamic, semi-dynamic or stationary charging [4,5]. In dynamic charging, an EV is charged while it travels over an electrified road [6]. Multiple scenarios of the electrification track are possible: it can be formed by a single elongated primary coil (similar to materials handling systems), or by segmented track, where each segment powered or switched individually, or by a sequence of pads similar to those used for EV charging [6,7,8].
Semi-dynamic charging is derived from the dynamic charging, but instead of powering vehicles in-move, they are powered at the intersections, traffic lights or bus stops (electrified public transport vehicles), where the vehicles are temporary stationary [4].
Stationary IPT chargers replace conductive chargers and can be fully integrated in the parking space, so that no part of the charger is visible [9]. Several standards were developed for stationary EV charging, e.g., SAE J2954 [10] or IEC 61980 [11]. The modern IPT systems for EV charging achieve similar levels of efficiency (>90%) and transferred power (tens and hundreds of kW) as the conductive chargers [2,12,13,14].
Unfortunately, the power transfer capability and efficiency of IPT systems are dependent on the coil alignment [15]. Thus, to ensure sufficient efficiency, the coils need to be correctly positioned [10], which requires feedback of their position or, preferably, the coupling coefficient k . Additionally, several control algorithms have been proposed, such as [16,17], which aim to provide power regulation and optimize efficiency under coil displacement. However, they require the value of the coupling coefficient k or mutual inductance M for operation.
There are several methods for estimating the coupling coefficient k [18,19,20,21], but they usually require knowledge of a large number circuit parameters [18,20,21] or require voltage or current information from both the transmitter and receiver [19]. If these parameters are measured inaccurately or they change (e.g., due to the displacement or temperature), the estimate quality will decrease. Also, some of the methods are based on complicated procedures requiring increased computational power [19,20].
In this paper, we present a preliminary evaluation of a novel primary-sided method for estimating the coupling coefficient, which is based solely on the frequency response of the input phase while the IPT system is intentionally operated in bifurcation. The proposed method also requires only a minimal number of calculations.
Bifurcation is a phenomenon during which the single frequency at which the input current is in phase with voltage, splits into three such frequencies. These frequencies are called the zero-phase angle (ZPA) frequencies. The original ZPA frequency is the resonance frequency f 0 , and the additional ZPA frequencies are denoted f φ 1 and f φ 2 in this work. The ZPA frequency f φ 1 is lower than f 0 , and f φ 2 is higher than f 0 . The onset of bifurcation can be initiated for example by appropriately adjusting the load or coupling coefficient k . As bifurcation increases, the ZPA frequencies f φ 1 and f φ 2 approach their respective asymptotes f 1 a and f 2 a , which are functions solely of the coupling coefficient k and f 0 .
The proposed method is developed and applied for the series–series (SS) compensated IPT system. The most practical way to achieve deep bifurcation for the SS compensation is to short-circuit the output. Under such a condition, the frequency response of the input phase φ i n is measured to find the ZPA frequencies f φ 1 f 1 a and f φ 2 f 2 a , which are then used to calculate the coupling coefficient k .
Since the method requires shorting the load, the method would not be used while power is being transferred to the EV. Instead, this method is applicable for initial estimation of k during start-up (i.e., before the charging begins) or intermittent estimation of k (i.e., in between periods of charging). In practice, the load could be short circuited using an active rectifier (in which both of the bottom transistors are turned on), through the DC/DC converter between the rectifier and EV battery, or a dedicated switch at the DC side after the rectifier.
In comparison with the current estimation methods, the proposed method requires (coarse) measurement of the frequency response of the input phase. However, its main advantage is that it does not require a-priori information about the system (except for an approximate frequency range for the frequency response measurement). It also requires only the primary side information (i.e., no wireless feedback between the secondary and primary is needed). The preliminary measurement results presented in this paper also show that the method is resilient to the detuning due to change of inductances in displacement. The method’s estimate of k only differed from the reference value by an average absolute error of 3.62% over the range k = 0.08 to k = 0.36.
The second chapter of the paper describes the initial AC equivalent circuit model, of which the method is based. The third chapter follows by deriving of the asymptotes of the ZPA frequencies f φ 1 and f φ 2 . The fourth chapter outlines the presented method itself, which is then experimentally verified in the fifth chapter and compared to the selected currently used methods for coupling coefficient estimation. The practical aspects of the selected methods are further discussed in the sixth chapter.

2. AC Equivalent Circuit

The proposed k -parameter estimation method is based on the AC equivalent circuit model of the IPT system depicted in Figure 1. The arrangement represents a single transmitter and receiver with the series–series compensation topology. The primary and secondary coil inductances L p and L s , respectively, are compensated by capacitances C p and C s to form series resonant circuits with resonant frequencies f p and f s . The frequencies f p and f s are set to the same value of f 0 , which is the common resonant frequency f 0 . The capacitances C p and C s are then calculated as:
C p = 1 ω 0 p 2 L p
C s = 1 ω 0 s 2 L s
where ω 0 p = ω 0 s = ω 0 = 2 π f 0 for the primary tuned to the secondary.
The coils are coupled by the mutual inductance M . The circuit is loaded by a purely resistive load, which is represented in Figure 1 by the equivalent load resistance R L . The equivalent series resistances (ESRs) representing losses of each side are summarized in R p and R s . The circuit is powered from a source, which is considered sinusoidal (first harmonic approximation). Therefore, all the voltages and currents are represented by phasors, which are marked as X ^ . The symbol X without this marking refers to phasor amplitude. The input voltage V ^ i n is described by its amplitude V i n , its frequency f (the operating frequency of the circuit) and its phase φ i n with respect to the input current I ^ i n , which is taken as the reference. V i n is considered constant. The IPT system output is characterized by the output voltage V ^ o u t and current I ^ o u t .
The input impedance of the circuit in Figure 1 is described by:
Z ^ = R e { Z ^ } + j · I m { Z ^ } = ω 4 C s 2 M 2 R s L ( ω 2 C s L s 1 ) 2 + ω 2 C s 2 R s L 2 + j ( ω L 1 1 ω C p ω 3 C s M 2 ( ω 2 C s L s 1 ) ( ω 2 C s L s 1 ) 2 + ω 2 C s 2 R s L 2 )

3. Theoretical Background of the Estimation Method

As discussed in the introduction, the phase bifurcation is the splitting of a single ZPA frequency into three. Specifically, the ZPA frequencies are described by condition:
I m { Z ^ } = 0
The formulas describing the phase bifurcation frequencies can be obtained from (4). To simplify the calculation, the normalization process presented by Wang et al. in [22] can be advantageously used. In the normalization, the ESRs R p and R s are neglected. The IPT model is normalized as follows. First, the normalized operating frequency u is defined as the operating frequency f divided by f 0 :
u = f f 0
Second, the imaginary part of the input impedance I m { Z ^ } is divided by the real part of the reflected impedance at resonance R e { Z ^ r 0 } to obtain normalized I m { Z ^ n } :
I m { Z ^ n } = I m { Z ^ } R e { Z ^ r 0 }
R e { Z ^ r 0 } = ω 0 2 M 2 R L
A series of substitutions is then applied after which I m { Z ^ } is expressed solely in terms of u and the loaded quality factors Q L p and Q L s , which are defined as:
Q L p = R L L p ω 0 M 2 ,   Q L s = ω 0 L s R L
In the last step, Equation (4) is transformed into the phase bifurcation equation by substituting I m { Z ^ n } for I m { Z ^ } :
( u 2 1 ) ( ( Q L p Q L s 2 Q L s ) u 4 + ( Q L p 2 Q L p Q L s 2 ) u 2 + Q L p Q L s 2 ) u ( u 2 1 ) 2 Q L s 2 + u 3 = 0
The detailed derivation of the phase bifurcation equation is described in [22].
When solved for u , Equation (9) has six roots, three positive ones and three negative ones. The negative roots are neglected as u cannot be negative. The first root u φ 0 = 1 describes the resonance frequency f 0 , and second and third roots describe the ZPA frequencies u φ 1 and u φ 2 , respectively:
u φ 1 = 2 Q L p Q L s 2 Q L p ( 4 Q L s 3 4 Q L p Q L s 2 + Q L p ) Q L p 2 Q L s ( Q L p Q L s 1 )
u φ 2 = 2 Q L p Q L s 2 + Q L p ( 4 Q L s 3 4 Q L p Q L s 2 + Q L p ) Q L p 2 Q L s ( Q L p Q L s 1 )
and when the roots u φ 1 and u φ 2 become real, the phase bifurcation occurs.
Analyzing Equations (10) and (11) for the ZPA frequencies u φ 1 and u φ 2 shows that these frequencies are asymptotic and that the asymptotes depend solely on the coupling coefficient k . The asymptotic behavior can be demonstrated by decreasing the equivalent load resistance R L towards zero. Figure 2 for IPT system described in Section 5.1. shows the course of the ZPA frequencies u φ 0 , u φ 1 and u φ 2 with decreasing R L while the remaining equivalent circuit parameters remain constant. As it can be seen, the ZPA frequencies u φ 1 and u φ 2 begin close to the resonance frequency u 0 . As R L decreases further, u φ 1 and u φ 2 recede until they stabilize at the asymptotes. The deep bifurcation, i.e., very low R L , is required for the roots of the bifurcation equation approach the asymptotes closely. The easiest way to achieve this is to short-circuit the output of the IPT system.
To obtain the equation for the coupling coefficient k calculation, it is necessary to obtain the formulas for the asymptotes. In the first step, the loaded quality factors Q L p and Q L s given by (8) are substituted in the roots u φ 1 and u φ 2 of the bifurcation equation, which are given by (10) and (11), respectively. Then their asymptotic values are calculated in the limit as R L goes to zero. In the next step the normalized frequency u given by (5) is substituted together with the coupling coefficient k :
k = M L 1 L 2
The final results are the equations for the asymptotes f 1 a and f 2 a :
f 1 a = f 0 1 + k
f 2 a = f 0 1 k
This result is interesting, because the asymptotes of the frequencies describing the frequency splitting of the output power P o u t maxima can be approximated by a similar set of formulas [23,24]. Our examination of the frequencies of the input current amplitude I i n maxima showed that they also have the same asymptotes.
The coupling coefficient k can be calculated from the asymptotes f 1 a and f 2 a given by (13) and (14), respectively, as:
k = f 2 a 2 f 1 a 2 f 2 a 2 + f 1 a 2
Note that the issue of coupling coefficient extraction is also important in distributed-element filter design; and, as discussed in [25], an equation identical to (16) is used to characterize coupled RF/microwave resonators based on the splitting of the scattering parameter S 21 .

4. Proposed Estimation Method (Bifurcation Asymptotes Method)

Equation (15) can be advantageously used for estimation of the coupling coefficient k . Thus, we propose the following procedure based on the bifurcation asymptotes, which estimates k using the frequency response of the input phase φ i n .
The bifurcation asymptotes method has the following steps:
(1)
Short-circuit the output of the IPT system.
(2)
Find the approximate values of the frequencies f 1 a   and f 2 a by searching for either zero crossings of φ i n or input current amplitude I i n maxima with a large step of f .
(3)
Measure φ i n with increased precision (smaller step of f ) near the approximated values of f 1 a and f 2 a .
(4)
For each value of f 1 a and f 2 a select the frequency of φ i n < 0 and φ i n > 0 nearest to φ i n = 0 . Use a linear approximation to calculate the zero-crossing frequency of φ i n ( φ i n = 0 ) from the selected frequencies and phases. The obtained frequencies are approximately equal to f 1 a and f 2 a .
(5)
Apply (15) to calculate the coupling coefficient k from the approximated values of f 1 a and f 2 a .
The individual steps are examined in further detail in the following discussion. In order to use (15), it is necessary to achieve such values of Q L p and Q L s that f φ 1 and f φ 2 would approach f 1 a and f 2 a , respectively. While there are multiple ways to achieve such a change of Q L p and Q L s , e.g., by decreasing f 0 , minimizing of R L by short-circuiting the output is the most practical.
Because the method is based on measuring the frequency response, a knowledge of a frequency range to examine is necessary. However, this range depends on k based on (13) and (14), which is unknown. We find that it is most convenient is to start at f 0 , and then extend the search for f 1 a and f 2 a both above and below f 0 . Therefore, an approximate value of f 0 is the only information about the IPT system required for the method. An IPT system with f 0 of 100 kHz was used for the method verification.
Measurement of the full frequency response with a small step of frequency is time consuming. This could render the method unusable for more dynamic applications. However, because we are only looking for two frequency values (i.e., f 1 a and f 2 a ), it is also unnecessary. To reduce the number of steps, two measures were taken.
At first, a rough approximation of frequencies f 1 a and f 2 a is generated by employing a large measurement step—in testing of the method a step of 1 kHz was used. In this search, the positions of f 1 a and f 2 a in the frequency range can be evaluated based on multiple criteria. For example, if the evaluation is based on the input phase φ i n , the positions of f 1 a and f 2 a can be localized based on the change of sign of φ i n . Frequencies f 1 a and f 2 a can be also approximated by searching for the frequencies corresponding to maxima of the input current amplitude I i n , as these have the same asymptotes as the ZPA frequencies f φ 1 and f φ 2 . Alternatively, if the primary side DC link current is measured, the frequencies of its maxima can be used as well, as they match I i n maxima. In the verification of the presented method, the frequencies of I i n maxima were used to generate a first approximation of f 1 a and f 2 a .
The second measure is to use the linear approximation to find f 1 a and f 2 a with increased precision. Figure 3 shows an example of the frequency response of the input phase φ i n for deep bifurcation (i.e., f φ 1 f 1 a and f φ 2 f 1 a ) with marked f 0 , f 1 a and f 2 a . As the figure shows, the shape of the frequency response φ i n is close to linear near the ZPA frequencies f φ 1 and f φ 2 .
Thus, the values of φ i n are measured near the approximated values of f 1 a and f 2 a obtained in the first step with a smaller frequency step. In the verification of the method, we measured the ±1 kHz surroundings of the approximated frequencies f 1 a and f 2 a with a step of 0.5 kHz. From the resulting data series, the frequencies and phases of the points nearest the zero crossing of φ i n were selected for each asymptote: the point f i a 1 , φ i a 1 with phase below zero and the point f i a + 1 , φ i a + 1 with phase above zero. Both points are marked by hollow circles in Figure 3. From these two points, the values of asymptote frequencies are approximated by:
f i a = f i a + 1 f i a 1 f i a + 1 φ i a 1 φ i a + 1 φ i a + 1
The obtained values of f 1 a and f 2 a (full circles in Figure 3) are then used to calculate the coupling coefficient according to (15).
The number of measured frequency points depends on k . For example, if k = 0.3 is measured and f 0 = 100 kHz, then f 1 a = 87.70 kHz and f 2 a = 119.52 kHz according to (13) and (14), respectively. The first approximation of f 1 a and f 2 a with a step of 1 kHz then requires measurement of 34 points and the second refining for the linear approximation requires 10 steps if the ±1 kHz surroundings of frequencies f 1 a and f 2 a are measured with a step of 0.5 kHz. In total, 44 points are measured. However, if the method is used to estimate k increasing from a small value (e.g., parking of an EV over a charging station), the number of points can be further reduced by using f 1 a and f 2 a estimated from the previous step in the first approximation instead of f 0 .
Because the asymptotes f 1 a (13) and f 2 a (14) are the same for phase bifurcation, frequency splitting of the output power P o u t maxima and splitting of input current amplitude I i n maxima [23,24], alternative approaches to finding the zero crossings of φ i n can be used for asymptote approximation, e.g., finding the frequencies of the maxima of I i n or I o u t . However, the linear approximation cannot be used to estimate the positions of the maxima; thus, the step of the frequency response would need to be much smaller, resulting in increased measurement time.

5. Experimental Verification

The estimates of k provided by the bifurcation asymptotes method are evaluated with respect to a reference value measured with a calibrated LRC bridge and compared with the results of other active (based on the measurement of currents and voltages) estimation methods. A direct method based on Faraday’s law (i.e., measuring the induced electromotive force (EMF)) and the method introduced by Jiwariyavej et al. in [18] are used for comparison.
In the first part of the chapter, the measurement setup is described together with the measured sequence of coupling coefficients, which is used for evaluation of all three estimation methods. Then, the procedures of the Faraday’s law method and Jiwariyavej’s method are outlined, and the estimation results are compared with the reference for all three selected methods. In the last part of the chapter, the match to the reference of all three methods is evaluated.

5.1. Measurement Setup

The measured setup used for verification is depicted in Figure 4 (photographs) and in Figure 5 (schematics). The setup was designed to closely resemble the AC equivalent circuit in Figure 1.
The circuit is powered by a linear amplifier driven by the function generator BK Precision 4063 which generates a sinusoidal voltage. The linear amplifier is supplied by power sources Diametral P230R51D. Both the primary and the secondary use DD pads [15,26]. Each coil inductance is compensated by a capacitive decade to provide control over f 0 p and f 0 s . Nominal resonance frequency was selected as 100 kHz.
The output is loaded by a custom-built resistive decade, designed for high frequency operation. However, for lower values of R L it has a slight inductive character, but this is considered in L s and the calculation of C s , which compensates it. The currents are calculated from the measured voltages at current shunts. One of two available oscilloscopes Rohde & Schwarz RTH1004 measures the voltages with RT-ZI10 probes and the shunt voltages corresponding with currents with P2220 probes set to 1:1 attenuation. The voltages are measured at designated points at the switchboards, which connect together all components of each side. The function generator and oscilloscope are controlled by an acquisition MATLAB script running on the laptop. The setup was designed with a focus on the IPT transformer operation and the bifurcation phenomena; thus, the linear amplifier and the resistor decade are used in place of an inverter and output rectifier, respectively. The component parameter values were measured with the LRC bridge Hameg HM8118 with Kelvin probe. The system parameters for the aligned coils (position x = 0 mm) are listed Table 1, except R L , which is specific for each method.

5.2. Measured Sequence Used for Method Evaluation

The presented method is evaluated together with selected existing methods by estimating the coupling coefficient k during the secondary pad displacement. The secondary pad is horizontally moved in the x-axis (see Figure 6) in 10 mm steps from the initial position of 0 to 200 mm. The step is decreased to 5 mm near the DD pads null position (see Figure 9c in [15]).
The coupling coefficient values used as a reference for method evaluation were measured by the following procedure using the calibrated LRC bridge.
(1)
Both pads are disconnected from the rest of the circuit.
(2)
The primary coil inductance L 1 a is measured while the secondary coil is open.
(3)
The secondary coil is short circuited. The primary inductance is measured again and labeled as L 1 b . L 1 b is related to L 1 a via:
L 1 b = L 1 a ( 1 k 2 )
(4)
The coupling coefficient k is then calculated from (17) as:
k = 1 L 1 b L 1 a
Figure 7 depicts the resulting k values (a) together with the primary and secondary inductances L 1 and L 2 (b). The figure clearly shows the null point of the DD pads near x = 75 mm. The resulting sequence of k covers the interval of 0.08 to 0.35, which is an aggregated range of k , which the IPT systems presented in SAE J2954 as examples for EV stationary charging gain under displacement [10]. As Figure 7b shows, the movement of the secondary coil changes the coil inductances, which in turn detunes the primary and secondary circuits. The impacts of detuning will be discussed for each method individually.
The methods for estimating k are compared with the reference measurement presented in Figure 7a. For each displacement position, the relative error with respect to the reference value is calculated according to:
Δ k r = k M k r e f k r e f [ % ]
where k M is the value estimated by the method and k r e f is the reference value.

5.3. Bifurcation Asymptotes Method

The proposed method for estimating k was implemented and tested according to the procedure outlined in Section 4. The method employed the setup shown in Figure 5, only the load was short-circuited at its connection point at the secondary switchboard. For each displacement position, the input current amplitude I i n was measured with a step size of 1 kHz to find the approximate values of f 1 a and f 2 a . These were then refined by measurement of the input phase φ i n measured with a step size of 0.5 kHz in interval ±1 kHz around the previously approximated frequencies f 1 a and f 2 a , and consequent linear interpolation. The obtained estimates of k “BA 0 mm” and “BA AT” are compared in Figure 8a with the reference “ref”. Label “BA 0 mm” marks the measurement with the system tuned to f 0 = 100 kHz in the initial position x = 0 mm and “BA AT” refers to the measurement when the system was manually tuned by adjusting the primary and secondary capacitances for each misalignment position to f 0 = 100 kHz (AT stands for always tuned).
For the system tuned at x = 0 mm, the method works well until approximately k = 0.17, as it has a relative error of approximately −2% (see Figure 8b), which is calculated according to (19). The negative sign means that the method underestimates k in comparison with the reference.
As k decreases less than k = 0.17, the relative error becomes positive, and its value increases up to 10%. For k smaller than approximately 0.08, the method fails altogether. However, for most practical systems [10], k is larger than 0.08 (typically, k > 0.2 is expected for high efficiency) so the method’s breakdown at low values of k is not expected to pose a limitation in practice.
The increase of the relative error for k below 0.17 and the method’s failure for k lower than 0.08 is primarily caused by detuning. Figure 9 shows that with increasing displacement the primary resonance frequency f 0 p increases above the secondary resonance frequency f 0 s . This detuning alters the frequency response of φ i n . Because f 0 p > f 0 s , the frequency response moves below zero in comparison with the case f 0 p = f 0 s —see Figure 10a for the tuned system and Figure 10b for the detuned system.
Under the displacement, the frequency responses change both with detuning and with changing k , as illustrated in Figure 11. In deep bifurcation, the frequency response of φ i n is relatively sharp near the zero crossings and it approaches its limiting values of −90° and +90° between the crossings if k is sufficiently high (see Figure 11a,b). Thus, the impacts of detuning are smaller for high k , as even with the shift the frequency response resembles a straight line at the zero crossing and can be approximated as a linear function there (see Figure 11b). However, as k decreases due to displacement, the frequency response curve becomes less pronounced, and it begins to cross zero near its peak, where it does not resemble a straight line (see Figure 11c). This increases the relative error. With a further decrease of k and increase of detuning, former zero crossings are shifted below zero and only a single zero crossing occurs—Figure 11d. Thus, for low k and increased detuning, the method fails.
The method’s performance for low k can be improved by maintaining the tuned system for specific inductance values under displacement by adjusting the primary and secondary compensation capacitances. The compensation adjustment can be done by e.g., switched capacity, or manually, as in our case. By removing the detuning, it was possible to obtain an approximate value of k in positions 65, 70, 180, and 190 mm—see “BA AT” in Figure 8a. However, even though as k decreases the measurement error grows up to 30%, as shown in Figure 8b. The high error is caused by the fact that the bifurcation is not deep enough such that frequencies f φ 1 and f φ 2 are not adequately meeting their asymptotes f 1 a and f 2 a , respectively. After k decreases under 0.04, not even tuning for a specific displacement position could help to achieve any meaningful results.
In the remaining points (75, 80, and 200 mm), only one zero crossing of the input phase was measured, even when the measurement step was decreased to 0.25 kHz. This may be caused either because the lack of bifurcation or the measurement step is still too coarse.
In applying the bifurcation asymptotes method, it is necessary to consider the impacts of the bifurcation phenomena on the input impedance Z ^ . When short-circuiting the output, the input impedance is relatively high at f 0 (e.g., Z = 99.3 Ω for x = 0 mm). However, near the ZPA frequencies f φ 1 and f φ 2 the impedance is very low (e.g., Z = 0.73 Ω), which is determined by the ESRs R p and R s . Thus, to limit the currents to 5 to 7 A in the primary side corresponding with the regular operation, the input voltage V i n was decreased to 5 V when applying the method.

5.4. Faraday’s Law Method

For comparison, we consider another common approach for extracting k by directly measuring the induced EMF in the secondary coil and solving for M using Faraday’s law of electromagnetic induction. The adjusted measurement setup is depicted in Figure 12. Both pads are placed in their positions. The secondary coil is open, and the primary coil is in series with a resistor and a capacitor, which provide control of the coil current I 1 . The primary coil current was limited between 5.8 and 6.6 A during the measurement process.
The primary coil is excited with current I 1 , which is measured with a precision current shunt. The I 1 value is controlled with a resistor and a capacitor, which reduces the inductive character. The value of I 1 is adjusted to equal the value of I i n during regular operation, to consider ferrite saturation and similar effects. The voltage V 2 is measured across the secondary coil terminals, and the mutual inductance M is calculated as:
M = V 2 ω I 1
The coupling coefficient k is calculated from M and the measured coil inductances L 1 , L 2 according to (12).
The estimated values of k obtained by Faraday’s law method are compared with the reference measurements in Figure 13a. The relative error depicted in Figure 13b is very low even for low k values. The distribution of error shows that the method generally estimates k approximately by 1.55% above the reference. Faraday’s law method is unaffected by detuning.

5.5. Method of Jiwariyavej

The second method used for comparison was presented by Jiwariyavej et al. in [18]. The variant for a single receiver was used. This method does not directly measure k but the mutual inductance M . The calculation of M is based on the equation for the real part of input impedance:
R e { Z ^ } = R p + ( ω M ) 2 R L + R s ( R L + R s ) 2 + X s 2
where X s is the secondary side reactance:
X s = ω L s 1 ω C s
The mutual inductance M is then calculated as:
M = 1 ω ( R e { Z i n } R p ) ( ( R L + R s ) 2 + X s 2 ) R L + R s
The real part of input impedance R e { Z i n } is calculated from the input current and voltage amplitudes I i n and V i n , respectively, and the input phase:
R e { Z i n } = V i n I i n c o s φ i n
The coupling coefficient k is then calculated according to (12). Obtaining the value of k according to this method requires preliminary measurement of R L , R p , R s , L p , L s , and C s .
The method was tested using the measurement setup shown in Figure 5, with the resistive load set to 31.8 Ω. Due to the higher number of preliminary parameter measurements, the method was measured in reduced number of points. The primary and secondary capacitances were adjusted to keep the system tuned at resonance frequency f 0 = 100 kHz against variations in coil inductance; see the illustration in Figure 7b. The specific values of the circuit parameters measured for each displacement position were used to calculate M according to (23) and consequently k according to (12).
Values of k obtained using the method of Jiwariyavej are compared to the reference values in Figure 14a. In addition to the plotted values, the method was also applied at displacement of 75, 190, and 200 mm, but did not provide any useful results. Thus, despite tuning for each displacement position, the method of Jiwariyavej suffers a similar limitation for small values of k as the bifurcation asymptote method, which was tuned only for x = 0 mm.
The relative error of the method of Jiwariyavej is depicted in Figure 14b. For high k values (above approximately 0.23), this method exhibits good agreement with the reference. On the other hand, as k is decreased, the error grows significantly and reaches 19.3% at k = 0.08.

5.6. Measurement Comparison

A comparative evaluation of the methods for estimating k is summarized in Figure 15. Figure 15a shows how the methods compare directly to the reference and Figure 15b shows their relative errors.
The evaluation of the methods is divided in two intervals k 0.08 , 0.23 and k 0.23 , 0.36 based on the relative error results shown in Figure 15b. The bottom limit of k = 0.08 was selected as the final value, at which point all three methods yielded a result. Thus, only the bifurcation asymptote method with tuning at x = 0 mm was evaluated. For each interval, the method is evaluated based on the average error Δ k ¯ r and average absolute error | Δ k ¯ r | , which are calculated as:
Δ k ¯ r = 1 n M i = 1 n M Δ k r , i = 1 n M i = 1 n M k M , i k r e f , i k r e f , i [ % ]
Δ k ¯ r = 1 n M i = 1 n M | Δ k r , i | = 1 n M i = 1 n M | k M , i k r e f , i | k r e f , i [ % ]
while | Δ k ¯ r | quantifies the disagreement between the method’s results from the reference in the selected interval. Δ k ¯ r describes the shift of the estimate with respect to the reference, i.e., if Δ k ¯ r 0 then the values are distributed evenly around the reference. The results are summarized in Table 2.
All three estimation methods show a good match to the reference for k above 0.23. From these values it is evident that the Jiwariyavej’s method exhibits the best match to the reference ( | Δ k ¯ r | = 0.20 %), with the estimates which are distributed evenly around the reference ( Δ k ¯ r = 0.04 %). The estimates of the method based on Faraday’s law are shifted above the reference ( Δ k ¯ r = 1.91 %), while the estimates of the bifurcation asymptotes method are shifted below the reference ( Δ k ¯ r = 2.01 %).
However, with decreasing k all methods except Faraday’s law method decrease in precision (interval k 0.08 , 0.23 ). Faraday’s law method is the most accurate ( | Δ k ¯ r | = 1.39 %), but its estimates are again shifted above the reference ( Δ k ¯ r = 1.39 %). The second-best method is the bifurcation asymptotes method ( | Δ k ¯ r | = 4.36 %) which is now shifted above the reference ( Δ k ¯ r = 3.87 %). The Jiwariyavej’s method performs much worse ( | Δ k ¯ r | = 10.38 %) and its value begins to move further below the reference with decreasing k —see Figure 15b, Δ k ¯ r = 9.30 %.
Across the whole measured interval of k 0.08 , 0.36 , Faraday’s law method is the most accurate ( | Δ k ¯ r | = 1.55 %), followed by the bifurcation asymptote method ( | Δ k ¯ r | = 3.62 %) and then Jiwariyavej’s method ( | Δ k ¯ r | = 6.31 %). The estimates of the two first methods are on average shifted above the reference ( | Δ k ¯ r | = 1.55 % and | Δ k ¯ r | = 3.62 %, respectively), while the third is shifted below the reference ( Δ k ¯ r = 5.6 %).
The shifts of the estimates with respect to the reference could be partially explained by the different measurement currents used in each method (see Table 3) and the reference, which was measured with the LRC meter and thus with very small currents. The high currents could result in detuning caused by the dependance of the coil inductances on the currents through the coils and ferrite non-linearity [27], as the ferrites are used as the magnetic cores of the DD pads. This effect could especially impact the Jiwariyavej’s method, as the primary currents were significantly higher in the second interval.
Some of the differences could also be explained as a result of measurement errors due to limited instrument accuracy, loading effects, and environmental factors such as temperature, etc.

6. Discussion

As evaluated in the previous chapter, the bifurcation asymptote method provides relatively accurate estimates ( | Δ k ¯ r | = 3.62 %) across the interval of k 0.08 , 0.36 . However, the estimate precision is not the only consideration for evaluating the usability of the method. Table 4 summarizes some of the practical aspects of all three evaluated methods, which should also be considered.
While the method based on Faraday’s law provides very accurate estimates, it requires the ability to disconnect the secondary coil from the rest of the circuit and measure the coil voltage across its open terminals. Thus, this method is more suitable for characterizing the coils in the laboratory setting than performing real-time estimations of k for practical applications such as EV charging.
Jiwariyavej’s method for a single receiver also provides accurate estimates, but only for higher k values: k 0.08 , 0.36 . Another disadvantage is that it requires additional preliminary measurements including L s to calculate the value of M . Due to this, use of Jiwariyavej’s method during operation is limited to air-core coils in which the inductance does not vary with displacement. The method is not very useful for the coils with magnetic cores, as their inductances change with displacement.
The main downside of the bifurcation asymptotes method is that it requires short- circuiting of the load; therefore, it cannot provide estimation of k during the operation of an IPT system. However, the method can be used to estimate k intermittently during a charging cycle and at critical times such as at start-up. In addition, the method involves measuring the frequency response; thus, a variable frequency source is required, and a certain amount of time given by transients must be allotted for the extraction process to occur.
In return, however, the bifurcation asymptotes method provides significant advantages. Namely, the method requires no a-priori information of the system component values, it is based purely on primary-side measurements, it estimates k with good accuracy ( | Δ k ¯ r | = 3.62 % for k 0.08 , 0.36 ), and it is resilient to the detuning caused by the displacement of the coils with magnetic cores. Thanks to this, the bifurcation asymptotes method presents an interesting option for estimating k that can be readily employed by control strategies that require an accurate estimation of this variable. The method is especially suitable for implementation at the startup of stationary charging applications. For example, the method can be used to estimate k immediately after an EV is parked over a stationary wireless charger to determine if the EV is suitably positioned to allow for efficient wireless charging before high power is applied to the primary.

7. Conclusions

This paper presented a novel bifurcation asymptotes method for accurately estimating the coupling coefficient k for SS compensated IPT systems. The method involves temporarily operating the system in bifurcation by short-circuiting the load and measuring the zero phase angle frequencies at the primary. The method can be used by control algorithms that require the value of k , and is especially suitable for implementation at the startup of certain applications such as EV stationary charging. This paper presents a preliminary evaluation of the method, which contains the theoretical background and description of the method, experimental verification, and a comparison with existing k -estimation methods. The proposed method provides good estimates (average absolute error of 3.62 %) for a wide range of coupling coefficients ( k 0.08 , 0.36 ).
For future work, we intend to evaluate the method performance in a non-idealized setup by implementing it in an inverter. We also intend to explore how this method can be extended to other compensation topologies such as series–parallel, LCC–series, etc. Finally, we intend to perform more comparisons with current k -estimation methods (e.g., [18,19] cited in the Introduction).

Author Contributions

Conceptualization, M.K.; methodology, M.K. and A.D.S.; software, M.K.; validation, M.K. and A.D.S.; formal analysis, M.K.; investigation, M.K.; resources, M.K. and J.L.; writing—original draft preparation, M.K. and A.D.S.; writing—review and editing, M.K., A.D.S. and J.L.; visualization, M.K.; supervision, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank Pavel Skarolek for many constructive discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, S.; Mi, C. Wireless Power Transfer for Electric Vehicle Applications. IEEE J. Emerg. Sel. Top. Power Electron. 2015, 3, 4–17. [Google Scholar] [CrossRef]
  2. Feng, H.; Tavakoli, R.; Onar, O.C.; Pantic, Z. Advances in High-Power Wireless Charging Systems: Overview and Design Considerations. IEEE Trans. Transp. Electrif. 2020, 6, 886–919. [Google Scholar] [CrossRef]
  3. Covic, G.A.; Boys, J.T. Inductive Power Transfer. Proc. IEEE 2013, 101, 1276–1289. [Google Scholar] [CrossRef]
  4. Lukic, S.; Pantic, Z. Cutting the Cord: Static and Dynamic Inductive Wireless Charging of Electric Vehicles. IEEE Electrif. Mag. 2013, 1, 57–64. [Google Scholar] [CrossRef]
  5. Patil, D.; McDonough, M.K.; Miller, J.M.; Fahimi, B.; Balsara, P.T. Wireless Power Transfer for Vehicular Applications: Overview and Challenges. IEEE Trans. Transp. Electrif. 2018, 4, 3–37. [Google Scholar] [CrossRef]
  6. Mi, C.C.; Buja, G.; Choi, S.Y.; Rim, C.T. Modern Advances in Wireless Power Transfer Systems for Roadway Powered Electric Vehicles. IEEE Trans. Ind. Electron. 2016, 63, 6533–6545. [Google Scholar] [CrossRef]
  7. Shin, J.; Shin, S.; Kim, Y.; Ahn, S.; Lee, S.; Jung, G.; Jeon, S.J.; Cho, D.H. Design and Implementation of Shaped Magnetic-Resonance-Based Wireless Power Transfer System for Roadway-Powered Moving Electric Vehicles. IEEE Trans. Ind. Electron. 2014, 61, 1179–1192. [Google Scholar] [CrossRef]
  8. Miller, J.M.; Jones, P.T.; Li, J.M.; Onar, O.C. ORNL Experience and Challenges Facing Dynamic Wireless Power Charging of EV’s. IEEE Circuits Syst. Mag. 2015, 15, 40–53. [Google Scholar] [CrossRef]
  9. Boys, J.T.; Covic, G.A. The Inductive Power Transfer Story at the University of Auckland. IEEE Circuits Syst. Mag. 2015, 15, 6–27. [Google Scholar] [CrossRef]
  10. Wireless Power Transfer for Light-Duty Plug-In/Electric Vehicles and Alignment Methodology; SAE TIR J2954 NOV2017; SAE International: Warrendale, PA, USA, 2017.
  11. Electric Vehicle Wireless Power Transfer (WPT) Systems—Part 1: General Requirements; IEC 61980-1:2020; International Electrotechnical Commission: London, UK, 2020.
  12. Schumann, P.; Diekhans, T.; Blum, O.; Brenner, U.; Henkel, A. Compact 7 kW inductive charging system with circular coil design. In Proceedings of the International Electric Drives Production Conference (EDPC), Nuremberg, Germany, 15–16 September 2015; pp. 1–5. [Google Scholar] [CrossRef]
  13. Galigekere, P.V.; Pries, J.; Onar, C.O.; Su, J.G.; Anwar, S.; Wiles, R.; Seiber, L.; Wilkins, J. Design and Implementation of an Optimized 100 kW Stationary Wireless Charging System for EV Battery Recharging. In Proceedings of the 2018 IEEE Energy Conversion Congress and Exposition (ECCE), Portland, OR, USA, 23–27 September 2018; pp. 3587–3592. [Google Scholar] [CrossRef]
  14. Pries, J.; Galigekere, V.P.N.; Onar, O.C.; Su, G.J. A 50-kW Three-Phase Wireless Power Transfer System Using Bipolar Windings and Series Resonant Networks for Rotating Magnetic Fields. IEEE Trans. Power Electron. 2020, 35, 4500–4517. [Google Scholar] [CrossRef]
  15. Budhia, M.; Boys, J.T.; Covic, G.A.; Huang, C.Y. Development of a single-sided flux magnetic coupler for electric vehicle IPT charging systems. IEEE Trans. Ind. Electron. 2013, 60, 318–328. [Google Scholar] [CrossRef]
  16. Zakerian, A.; Vaez-Zadeh, S.; Babaki, A.; Moghaddam, M.F. Efficiency Optimization of a Dynamic Wireless EV Charging System Using Coupling Coefficient Estimation. In Proceedings of the 2019 10th International Power Electronics, Drive Systems and Technologies Conference (PEDSTC), Shiraz, Iran, 12–14 February 2019; pp. 629–634. [Google Scholar] [CrossRef]
  17. Liu, F.; Chen, K.; Zhao, Z.; Li, K.; Yuan, L. Transmitter-Side Control of Both the CC and CV Modes for the Wireless EV Charging System with the Weak Communication. IEEE J. Emerg. Sel. Top. Power Electron. 2018, 6, 955–965. [Google Scholar] [CrossRef]
  18. Jiwariyavej, V.; Imura, T.; Hori, Y. Coupling Coefficients Estimation of Wireless Power Transfer System via Magnetic Resonance Coupling Using Information From Either Side of the System. IEEE J. Emerg. Sel. Top. Power Electron. 2015, 3, 191–200. [Google Scholar] [CrossRef]
  19. Kobayashi, D.; Imura, T.; Hori, Y. Real-time coupling coefficient estimation and maximum efficiency control on dynamic wireless power transfer for electric vehicles. In Proceedings of the 2015 IEEE PELS Workshop on Emerging Technologies: Wireless Power (2015 WoW), Daejeon, Korea, 5–6 June 2015; pp. 1–6. [Google Scholar] [CrossRef]
  20. Chow, J.P.-W.; Chung, H.S.-H.; Cheng, C.-S. Use of Transmitter-Side Electrical Information to Estimate Mutual Inductance and Regulate Receiver-Side Power in Wireless Inductive Link. IEEE Trans. Power. Electron. 2016, 31, 6079–6091. [Google Scholar] [CrossRef]
  21. Nawada, H.; Takahashi, Y.; Hata, K.; Imura, T.; Fujimoto, H.; Hori, Y.; Yabumoto, T. Coupling Coefficient Estimation for Wireless Power Transfer System at Constant Input Power Operation. In Proceedings of the 2019 IEEE PELS Workshop on Emerging Technologies: Wireless Power Transfer (WoW), London, UK, 17–21 June 2019; pp. 288–291. [Google Scholar] [CrossRef]
  22. Wang, C.-S.; Covic, G.A.; Stielau, O.H. Power transfer capability and bifurcation phenomena of loosely coupled inductive power transfer systems. IEEE Trans. Ind. Electron. 2004, 51, 148–157. [Google Scholar] [CrossRef]
  23. Niu, W.Q.; Chu, J.X.; Gu, W.; Shen, A.D. Exact Analysis of Frequency Splitting Phenomena of Contactless Power Transfer Systems. IEEE Trans. Circuits. Syst. I. Regul. Pap. 2013, 60, 1670–1677. [Google Scholar] [CrossRef]
  24. Terman, F.E. Radio Engineers’ Handbook; McGraw-Hill: New York, NY, USA, 1943. [Google Scholar]
  25. Hong, J.-S.; Lancaster, M.J. Microstrip Filters for RF/Microwave Applications; John Wiley & Sons, Inc.: New York, NY, USA, 2001. [Google Scholar]
  26. Mohamed, A.A.S.; Shaier, A.A.; Metwally, H.; Selem, S.I. A comprehensive overview of inductive pad in electric vehicles stationary charging. Appl. Energy 2020, 262, 114584. [Google Scholar] [CrossRef]
  27. Bartoli, M.; Reatti, A.; Kazimierczuk, M.K. High-frequency models of ferrite core inductors. In Proceedings of the IECON’94-20th Annual Conference of IEEE Industrial Electronics, Bologna, Italy, 5–9 September 1994; pp. 1670–1675. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Equivalent AC circuit schematics for the series–series compensated IPT system.
Figure 1. Equivalent AC circuit schematics for the series–series compensated IPT system.
Electronics 10 02548 g001
Figure 2. The course of bifurcation with decreasing load resistance R L . The normalized frequency u 0 corresponds to the resonance frequency f 0 , while u φ 1 and u φ 2 with higher ZPA frequencies f φ 1 and f φ 2 , respectively. Their asymptotes are u 1 a and u 2 a , respectively.
Figure 2. The course of bifurcation with decreasing load resistance R L . The normalized frequency u 0 corresponds to the resonance frequency f 0 , while u φ 1 and u φ 2 with higher ZPA frequencies f φ 1 and f φ 2 , respectively. Their asymptotes are u 1 a and u 2 a , respectively.
Electronics 10 02548 g002
Figure 3. Frequency response of the input phase φ i n for the IPT system in deep bifurcation, with marked measured points, linear approximation and interpolated zero crossings. Based on the circuit parameters from Table 1.
Figure 3. Frequency response of the input phase φ i n for the IPT system in deep bifurcation, with marked measured points, linear approximation and interpolated zero crossings. Based on the circuit parameters from Table 1.
Electronics 10 02548 g003
Figure 4. Measurement setup. (a) Test stand for coil positioning. 1—Test stand. 2—Pads in the aligned position. (b) Setup without the test stand. 3—Laptop. 4—Function generator. 5—Linear amplifier. 6—Laboratory power sources. 7—Total stop. 8—Primary switchboard. 9—Capacitive decade used as the primary compensation. 10—Cable connecting the primary pad. 11—Cable connecting the secondary pad. 12—Secondary switchboard. 13—Capacitive decade used as the secondary compensation. 14—Resistive decade used as the load. 15—Oscilloscopes. 16—LRC bridge.
Figure 4. Measurement setup. (a) Test stand for coil positioning. 1—Test stand. 2—Pads in the aligned position. (b) Setup without the test stand. 3—Laptop. 4—Function generator. 5—Linear amplifier. 6—Laboratory power sources. 7—Total stop. 8—Primary switchboard. 9—Capacitive decade used as the primary compensation. 10—Cable connecting the primary pad. 11—Cable connecting the secondary pad. 12—Secondary switchboard. 13—Capacitive decade used as the secondary compensation. 14—Resistive decade used as the load. 15—Oscilloscopes. 16—LRC bridge.
Electronics 10 02548 g004aElectronics 10 02548 g004b
Figure 5. Measurement setup block diagram.
Figure 5. Measurement setup block diagram.
Electronics 10 02548 g005
Figure 6. Secondary pad displacement in x axis.
Figure 6. Secondary pad displacement in x axis.
Electronics 10 02548 g006
Figure 7. Reference values for x-axis displacement. (a) Coupling coefficient k. (b) The primary and secondary inductances L1 and L2.
Figure 7. Reference values for x-axis displacement. (a) Coupling coefficient k. (b) The primary and secondary inductances L1 and L2.
Electronics 10 02548 g007
Figure 8. Bifurcation asymptotes method. (a) Comparison with the reference measurement. (b) Relative error with respect to the reference over k. “BA 0 mm” refers to system tuned for x = 0 mm. “BA AT” refers to system tuned for each misalignment position.
Figure 8. Bifurcation asymptotes method. (a) Comparison with the reference measurement. (b) Relative error with respect to the reference over k. “BA 0 mm” refers to system tuned for x = 0 mm. “BA AT” refers to system tuned for each misalignment position.
Electronics 10 02548 g008
Figure 9. Detuning of the primary and secondary side.
Figure 9. Detuning of the primary and secondary side.
Electronics 10 02548 g009
Figure 10. Frequency response shift due to detuning, illustrated on the displacement position of x = 60 mm. (a) Both sides are tuned to f0 = 100 kHz. (b) The side resonance frequencies are given by Figure 9, i.e., f0p = 103.27 kHz and f0s = 100.13 kHz.
Figure 10. Frequency response shift due to detuning, illustrated on the displacement position of x = 60 mm. (a) Both sides are tuned to f0 = 100 kHz. (b) The side resonance frequencies are given by Figure 9, i.e., f0p = 103.27 kHz and f0s = 100.13 kHz.
Electronics 10 02548 g010
Figure 11. Influence of detuning on the bifurcation asymptote method. (a) Tuned system (x = 0 mm). (b) Detuned system with medium coupling coefficient k (x = 130 mm). (c) Detuned system with small k (x = 60 mm). (d) Detuned system with very small k (x = 180 mm). The discrepancies between the measured data points and calculated frequency response (especially in (d)) are caused by high sensitivity of frequency responses for low k and the influence of the ferrites.
Figure 11. Influence of detuning on the bifurcation asymptote method. (a) Tuned system (x = 0 mm). (b) Detuned system with medium coupling coefficient k (x = 130 mm). (c) Detuned system with small k (x = 60 mm). (d) Detuned system with very small k (x = 180 mm). The discrepancies between the measured data points and calculated frequency response (especially in (d)) are caused by high sensitivity of frequency responses for low k and the influence of the ferrites.
Electronics 10 02548 g011
Figure 12. Faraday’s law method for measurement setup.
Figure 12. Faraday’s law method for measurement setup.
Electronics 10 02548 g012
Figure 13. Method based on Faraday’s law. (a) Comparison with the reference measurement. (b) Relative error with respect to the reference over k.
Figure 13. Method based on Faraday’s law. (a) Comparison with the reference measurement. (b) Relative error with respect to the reference over k.
Electronics 10 02548 g013
Figure 14. Method presented by Jiwariyavej in [18]. (a) Comparison with the reference measurement. (b) Relative error with respect to the reference over k.
Figure 14. Method presented by Jiwariyavej in [18]. (a) Comparison with the reference measurement. (b) Relative error with respect to the reference over k.
Electronics 10 02548 g014
Figure 15. (a) Comparison of all tested estimation methods with the reference. (b) Relative error of the estimation methods with respect to the reference over k.
Figure 15. (a) Comparison of all tested estimation methods with the reference. (b) Relative error of the estimation methods with respect to the reference over k.
Electronics 10 02548 g015
Table 1. System parameters in the position x = 0 mm.
Table 1. System parameters in the position x = 0 mm.
ParameterSymbolValue
Input voltage V i n 5 V
Primary coil inductance L p 28.8 µH
Primary compensation capacitance C p 87.4 nF
Primary ESR R p 0.23 Ω
Secondary coil inductance L s 16.7 µH
Secondary compensation capacitance C s 151.2 nF
Secondary ESR R s 0.15 Ω
Resonance frequency f 0 100 kHz
Table 2. Overview of the average relative errors Δ k ¯ r and average absolute errors | Δ k ¯ r | of the evaluated methods.
Table 2. Overview of the average relative errors Δ k ¯ r and average absolute errors | Δ k ¯ r | of the evaluated methods.
Method
k 0.08 , 0.23  
k 0.23 , 0.36
k 0.08 , 0.36
Bifurcation asymptote method
Δ k ¯ r = 3.87 % | Δ k ¯ r | = 4.36 %
Δ k ¯ r = 2.01 % | Δ k ¯ r | = 2.01 %
Δ k ¯ r = 2.04 % | Δ k ¯ r | = 3.62 %
Faraday’s law method
Δ k ¯ r = 1.39 % | Δ k ¯ r | = 1.39 %
Δ k ¯ r = 1.91 % | Δ k ¯ r | = 1.91 %
Δ k ¯ r = 1.55 % | Δ k ¯ r | = 1.55 %
Jiwariyavej’s method
Δ k ¯ r = 9.30 % | Δ k ¯ r | = 10.38 %
Δ k ¯ r = 0.04 % | Δ k ¯ r | = 0.20 %
Δ k ¯ r = 5.6 % | Δ k ¯ r | = 6.31 %
Table 3. Currents used for testing of the estimation methods.
Table 3. Currents used for testing of the estimation methods.
Method I i n I o u t
Bifurcation asymptote method5 to 7 A5 to 10 A
Faraday’s law method5 to 7 AOpen
Jiwariyavej’s method 4.5   to   6   A   ( k 0.23 , 0.36 )
10.5   A   ( k 0.08 , 0.23 )
0.5   to   1   A   ( k 0.23 , 0.36 )
0.5   A   ( k 0.08 , 0.23 )
Table 4. Comparison of the practical aspects of the bifurcation asymptotes method with the selected estimation methods.
Table 4. Comparison of the practical aspects of the bifurcation asymptotes method with the selected estimation methods.
MethodBifurcation Asymptotes Method Faraday’s Law MethodJiwariyavej’s Method [18]
Measured parameter k M M
Use during operationIntermittentNoYes, but only for air-core coils
Required active measurements f ,   φ i n ,   I i n (optional) f ,   I 1   ( primary   coil   current   amplitude ) ,   V 2 (secondary coil voltage amplitude) f ,   I i n ,   V i n ,   φ i n
Number of frequency points required Varies   with   k 11
Required preliminary measurementsNone
( general   knowledge   of   f 0 to specify measurement range)
None   to   obtain   M ,
L p ,   L s   to   obtain   k
R L ,   R p ,   R s ,   L s ,   and   C s   to   obtain   M ,
plus   L p   to   obtain   k
Other requirementsAbility to short-circuit the load.
Source with variable frequency
Decreased source voltage
Open secondary coil Air - core   coils ,   or   L p   and   L s must be remeasured after displacement
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Košík, M.; Scher, A.D.; Lettl, J. Novel Method of Coupling Coefficient Estimation Based on the Bifurcation Phenomena in Inductive Power Transfer. Electronics 2021, 10, 2548. https://doi.org/10.3390/electronics10202548

AMA Style

Košík M, Scher AD, Lettl J. Novel Method of Coupling Coefficient Estimation Based on the Bifurcation Phenomena in Inductive Power Transfer. Electronics. 2021; 10(20):2548. https://doi.org/10.3390/electronics10202548

Chicago/Turabian Style

Košík, Michal, Aaron D. Scher, and Jiří Lettl. 2021. "Novel Method of Coupling Coefficient Estimation Based on the Bifurcation Phenomena in Inductive Power Transfer" Electronics 10, no. 20: 2548. https://doi.org/10.3390/electronics10202548

APA Style

Košík, M., Scher, A. D., & Lettl, J. (2021). Novel Method of Coupling Coefficient Estimation Based on the Bifurcation Phenomena in Inductive Power Transfer. Electronics, 10(20), 2548. https://doi.org/10.3390/electronics10202548

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