Unilateral Route Method to Estimate Practical Mutual Inductance for Multi-Coil WPT System

Abstract

Multi-coil WPT systems require mutual inductance information between coils to increase the power transmission efficiency. However, in the high frequency (HF) bands such as 6.78 MHz and 13.56 MHz, the presence of surrounding coils changes the value of the mutual inductance between the two coils due to the parasitic element effect of the coils. These parasitic effects make it harder to estimate the mutual inductance among three or more coils. In contrast to ideal mutual inductance, which has a constant value regardless of frequency and surrounding coils, we define the practical mutual inductance as the mutual inductance varied by parasitic elements. In this paper, a new method is presented to estimate the practical mutual inductance between multiple coils in the HF band. The proposed method simply configures the expression of practical mutual inductance formula because only one of two bilateral dependent voltage sources generated by mutual inductance is considered. For several coils placed along the same axis, the practical mutual inductances between coils were measured with respect to the distance between them to validate the proposed method. The practical mutual inductance obtained from the proposed method was consistent with the simulated and measured values in HF band.

1. Introduction

Recently, multi-coil wireless power transfer (WPT) systems using three or more coils have been proposed to achieve a high power transfer efficiency (PTE) at greater distance than in the case of a conventional two-coil WPT system [1,2,3,4,5,6,7,8,9,10,11,12]. The multi-coil WPT systems are categorized into the domino structure and multiple transmitting (or receiving) coils structure. The domino structure means that one or more coils are placed between transmitting and receiving coils [1,2,3,4,5,6]. To achieve a long transmission distance, the structure adjusts the capacitance of each coil or gaps among coils. On the other hand, multiple transmitting (or receiving) coils correspond to the multiple input single output (MISO) or single input multiple output (SIMO) system in the wireless communication field [7,8,9,10,11,12]. The WPT system with multiple transmitting coils usually maximizes the power in the receiving coil through entering power of different magnitude or phase into transmitting coils. Sometimes, this technique is called magnetic beamforming [7,8].

The trouble in implementing these multi-coil WPT systems is that the mutual inductances (or coupling coefficients) among coils should be basically known. For the domino structure, the coupling coefficient information is used to estimate the target free-resonant frequencies, that is the resonant capacitances, or target distances between adjacent coils. For the multiple transmitting coils, power can be transferred into a not-aligned receiving coil by controlling the phase and magnitude of signal components input to each transmitting coil, where the mutual inductance information is used to calculate the phase and magnitude of input power [7,8]. This technique is referred to as magnetic beamforming. Similarly, for the multiple receiving coils, the output power can be maximized by adjusting the phase of the received power based on the mutual inductance information. Therefore, the accurate value of mutual inductance between coils is important to obtain high efficiency of multi-coil WPT systems such as domino, MISO, and SIMO structures in common. On the other hand, inaccurate mutual inductance causes interruption of maximum performance for the multi-coil systems.

The ideal mutual inductance between coils is completely separated from the operating frequency or presence of surrounding coils. In practice, when the mutual inductance between coils is measured, the measurement frequency and the presence of surrounding coils affect the measurement results. These phenomena can be inferred by the effect of the parasitic elements of the coil. We refer to the mutual inductance affected by the parasitic elements as the practical mutual inductance.

As is well known, the parasitic effect of a component is commonly frequency dependent, which becomes large as the frequency increases [13]. That is, the parasitic effects are not significant and can be ignored in low frequency (LF) bands such as 100–205 kHz. Therefore, most previous researches on the mutual inductance for the WPT system have focused on the estimation of the ideal mutual inductance of two coils with respect to the alignment status [14,15,16,17,18]. For WPT systems used in the HF band such as 6.78 MHz and 13.56 MHz, the practical mutual inductances between the transmitting and receiving coils are different from the theoretical value due to the effect of the parasitic component. S. Hackl et al. proposed an equivalent circuit model for the radio frequency identification (RFID) system with two coils operating at a very close distance considering not only the parasitic capacitors but also the capacitive coupling between coils [19]. However, the capacitive coupling between coils is generally taken account for 1 cm or less gap between coils, while the WPT system in the HF band is mainly used for a distance of tens or several cm. In addition, the method is difficult to expand to three or more coils because of complicated formulation.

In HF band, the practical mutual inductance between two coils can be easily measured using a vector network analyzer (VNA). However, it cannot be used to measure the practical mutual inductance among three or more coils since public VNAs usually have only two ports. Moreover, it is obvious that the number of transmitting or receiving coils tends to increase steadily. As the number of coils increases, the parasitic effect becomes greater. Therefore, an effective method is highly required to estimate the practical mutual inductance among multiple coils with parasitic effects using a two-port VNA in HF band. In this paper, a new method is proposed to obtain the practical mutual inductance among multiple coils using a general two-port VNA under the assumption that the coils are unilateral instead of reciprocal characteristics. The proposed method is suitable for general HF band WPT systems used at several cm distance and the simple expression of the formula allows to determine the effect of the presence of surrounding coils through every possible route among coils. We also verify the proposed method by comparison with measured and simulated mutual inductance among multiple coils. 2. Theoretical Analysis This section presents theoretical analysis for the proposed unilateral route method. 2.1. Practical and Ideal Mutual Inductance between Two Coils

Figure 1a,b shows the equivalent circuit model of two ideal coils and two coils with parasitic element, respectively. If the two coils are considered as a two-port network, the two-port network can be represented by various two-port parameters such as impedance parameter and scattering parameter. For the two ideal coils of Figure 1a, the mutual inductance between the two coils is independent of the operating frequency. From the definition of the impedance parameter, the transfer impedance between two ideal coils is given by

_Z21=_V2I1|I2=0=jω_M21,

where V_i and I_i are the voltage and current of the ith port, and ω is the angular frequency. M₂₁ is the mutual inductance between the coils and the two-port network is reciprocal so that Z_ij = Z_ji and M_ij = M_ji. Equation (1) means that Z_ij can be found by exciting port j with the current I_j, open-circuiting all other ports, and measuring the open-circuit voltage at port i. From Equation (1), the mutual inductance (M₂₁) between the coils can be easily obtained by measuring the transfer impedance (Z₂₁). This method is commonly used to measure the mutual inductance between the coils.

On the other hand, as the operating frequency increases, the coil is not presented as a single inductor but combined with parasitic elements such as resistor and capacitor. However, the parasitic capacitor is only considered as shown in the equivalent model of Figure 1b. This is because the parasitic resistor of coils is usually much smaller than the impedance (ωL) of the inductor in HF band, R << ωL₁, and then the reactance is dominant in determining the current flowing on the inductor and series parasitic resistor. Moreover, the parasitic capacitor of real coil connected in parallel makes a closed-loop, which causes the current path. Therefore, while the real coil is open-terminated in HF band, the practical mutual inductance between two coils is changed by the current path. From the circuit theory, the transfer impedance between the two coils with parasitic capacitors is given by

_Z21=_V2I1|I2=0=_{i1I1⋅v2i1⋅V2v2|I2=0}=1jω_C11jω_C1+jω_L1+(^ω2 _M12 ²1jω_C2+jω_L2)⋅jω_M21⋅1jω_C21jω_C2+jω_L2

The transfer impedance Z₂₁ is presented by multiplying the ideal transfer impedance jωM and other terms. From Equation (2), causes of the parasitic effect are categorized into two. First, not all of the input current (I₁) of port 1 flows through L₁ due to the presence of C₁ and only part (i₁) of the current flows; then, the induced voltage (v₂) on L₂ decreases. In addition, the induced voltage (v₂) generates the current (i₂) flowing closed current path, which consists of L₂ and C₂, even when port 2 is opened. Then, the current (i₂) flowing through L₂ induces the voltage (v₁) on L₁ by the mutual inductance. Consequently, this lowers the amount of the current (i₁) flowing through L₁ and also lowers the induced voltage (v₂). This effect is expressed as the first term of the right side of Equation (2). Second, the voltage (v₂) induced in Coil 2 is divided into L₂ and C₂ so that only a part of the induced voltage (V₂) is delivered into port 2. This voltage distribution is expressed by the third term on the right side of Equation (2).

Similar to Equation (1), we define the transfer impedance using the practical mutual inductance (M_{21(practical)}) as

_Z21=jω_{M21(practical)}.

Comparing Equation (3) with Equation (2) yields that M_{21(practical)} is obviously different with M₁₂ and strongly dependent on the coils’ inductance, L₁ and L₂, and the parasitic capacitance, C₁ and C₂, as well as the angular frequency, ω.

The term in parentheses in the first term denominator of Equation (2) means that some of the voltage induced from Coil 1 to Coil 2 is induced again to Coil 1 through mutual inductance. Except for special cases where the coils are very close to each other or operate near the self-resonant frequency (SRF) of the coils, it is small enough to be negligible compared to other terms. Therefore, Equation (2) can be approximately expressed as follows in combination with Equation (3)

_{M21(practical)}=_Z21jω≃1jω_C1jω_L1+1jω_C1⋅_M21⋅1jω_C2jω_L2+1jω_C2.

Equation (4) is equivalent to the circuit diagram of Figure 2 that shows the dependent voltage source on L₁ removed in Figure 1b. That is, it is assumed that the coils are unilateral. If the coils’ parameters, L₁, C₁, L₂ and C₂, are known, the ideal mutual inductance M₂₁ can be obtained by measuring the transfer impedance Z₂₁ and using Equation (4). Conversely, if the coil’s parameters and ideal mutual inductance M₂₁ are known, the practical mutual inductance M_{21(practical)} in the HF band can be estimated by using Equation (4).

2.2. Measurement Method for the Coil’s Parameters and Ideal Mutual Inductance between Two Coils

As mentioned above, the coil’s parameters and ideal mutual inductance should be known in order to estimate the practical mutual inductance between two coils. Figure 3 shows the proposed procedure to obtain the coil’s parameters and ideal mutual inductance. First, parameters of each single coil are measured without surrounding coils or the object. Since the parasitic effect rarely appears in the LF band, the input impedance Z₁₁ of a single coil is simply expressed as R₁ + jωL₁. Therefore, the self-inductance and resistance of each coil are directly obtained by measuring the input impedance of the coil using a VNA at the lowest possible frequency in the LF band. Next, the parasitic capacitance C₁ of the coil can be obtained from the SRF of the coil as given by

_C1=1_ω0 ² _L1=1^(2π_f0)2 _L1,

where f₀ is the SRF at which the reactance of the coil is zero, and it can be easily measured from the Smith Chart of S₁₁ or the input impedance plot on the VNA. Finally, the ideal mutual inductance M₂₁ between two coils can also be obtained from measuring the transfer impedance Z₂₁ between the coils using the VNA in the LF band. Therefore, the parameters of coils measured through the procedure of Figure 3 can be used to estimate the practical mutual inductance between two coils. Moreover, if the ideal mutual inductances among three or more coils are well defined, the measured parameters can also be used to estimate the practical mutual inductance.

2.3. Transfer Impedance of Multi-Coil Using the Proposed Circuit Model

For three coils shown in Figure 4a, we should estimate the transfer impedance Z₃₁ between Coils 1 and 3 to obtain the practical mutual inductance between the coils. From the definition of impedance parameter, the current applied to the first port can directly induce the voltage on the opened third port and also indirectly induce the voltage via the second port. Strictly, the induced voltage makes a current in the closed loop formed by the parasitic capacitor, and the current in the third inductor can affect other inductors. However, since the high order mutual effects are weak, under the assumption that the coils are unilateral, the transfer impedance between Coils 1 and 3 is obtained as the sum of two routes in Figure 4b. Thus, the practical mutual inductance between Coils 1 and 3 is given by

_{M31(practical)}=_Z31jω≃_ZR−1+_ZR−2jω,

where Z_R−1 and Z_R−2 are the transfer impedances of the direct and one-hop indirect routes between Coils 1 and 3, respectively, as shown in Figure 4b. Using Equation (4) the transfer impedance of the direct route is given by

_ZR−1=1jω_C1jω_L1+1jω_C1⋅jω_M31⋅1jω_C3jω_L3+1jω_C3.

From unilateral characteristics, the transfer impedance of the indirect one-hop route is also easily obtained by

_ZR−2=VI=_i1I⋅jω_M21⋅−1_Z2⋅jω_M32⋅V_v3=1jω_C1jω_L1+1jω_C1⋅jω_M21⋅−1jω_L2+1jω_C2⋅jω_M32⋅1jω_C3jω_L3+1jω_C3,

where M_21, M₃₂ and M₃₁ are ideal mutual inductances, and the minus sign is due to the dot notation between Coils 1 and 2.

For four coils shown in Figure 5, there are four possible routes from Coil 1 to Coil 4 under the unilateral assumption. The transfer impedance of the fourth route is given by

_ZR−4=VI=_i1I⋅jω_M21⋅−1_Z2⋅jω_M32⋅−1_Z3⋅jω_M43⋅V_v4.

The practical mutual inductance is summarized as

_{M41(practical)}=_Z41jω≈_ZR−1+_ZR−2+_ZR−3+_ZR−4jω.

If the number of coils is n, the transfer impedance is determined as

_Zn1≈∑k=1²ⁿ⁻²_ZR−k.

3. Experimental Results and Discussion

To verify the proposed method in the HF band, we analyzed the practical mutual inductance by changing the number of coils as shown in Figure 6. Two types of coils were used for the experiment, and parameters of the coils are summarized in Table 1. Coil_a has three turns and a single layer printed on an FR-4 substrate and dimensions 95.7 mm × 105.7 mm. Coil_b has six turns and two layers printed on an FR-4 substrate of the same dimensions as the Coil_a. All coils were arranged coaxially in order to maximize the parasitic effect and the effect of surrounding coils. The practical mutual inductance among multiple coils was measured with respect to the distance using the four-port VNA, Agilent’s E5071, as shown in Figure 7. Therefore, the [S] matrix for three coils that are directly connected to ports of the multi-port VNA was measured and then converted into the [Z] matrix using the equation in [20] on the Keysight’s ADS.

Additionally, we simulated the practical mutual inductance from the full circuit models with all parasitic components and dependent voltage sources using the ADS, and the measured practical mutual inductances are compared with the simulated results as well as those from the proposed method in several cases. 3.1. Two Coils Analysis

For two coils of Figure 6a, where the two coils used are identical to Coil_a in Table 1, the S parameters were measured without surrounding coils and the practical mutual inductances were estimated by substituting the Z parameters converted from S parameters into Equation (4). Figure 8 shows the measured and simulated practical mutual inductances with those of the proposed estimation method. As the distance between the coils decrease, the ideal mutual inductance increases and then the practical mutual inductance also increases. On the other hand, at 13.56 MHz, the practical mutual inductance has a higher level than 6.78 MHz. This result comes from the first and third terms of the right-hand side of Equation (4), being a function of frequency. As the frequency increases, the magnitude of the denominators of the terms dramatically increases more than those of the numerators. The mutual inductances estimated by the proposed method are in good agreement with the simulated and measured values over the distance. The closer the distance between the coils and the higher the frequency, the larger the practical mutual inductance becomes due to the parasitic effect. The proposed method accurately simulates these effects.

3.2. Three Coils Analysis

For three coils of Figure 6b, the experiments were conducted on the two cases of Table 2, where Coil 2 was located in the middle of Coils 1 and 3. The self-inductance of Coil 2 for Case#2 is larger than the value for Case#1. Figure 9 shows the simulated, measured, and estimated practical mutual inductances with respect to distance at 6.78 MHz and 13.56 MHz. It is noticed that the variation in the mutual inductance with Case#2 is greater than Case#1. This means that the larger the inductance of the surrounding coils, the greater the influence on the mutual inductance between other coils. In addition, the variation of the practical mutual inductance is greater at 13.56 MHz than at 6.78 MHz. That is, the parasitic effect is larger at a higher frequency. On the other hand, the proposed method exactly estimates the practical mutual inductance all over the distance, even if the parameters of the coil change.

To show the accuracy of the proposed method, we calculate the mean and normalized root-mean-square error (NRMSE) of difference between the results of measurement and the proposed method shown in Figure 8 and Figure 9, and summarize them in Table 3. While the mean and NRMSE increase as the number of coil and frequency increases, the mean and error are less than 20×10⁻¹⁰ (0.002 μH) and 3%, respectively, for all cases. This means that the proposed method accurately estimates the practical mutual inductance as measured by a multi-port VNA.

3.3. Analysis of Unilateral Method for the Number of Coils

Figure 10 shows the practical mutual inductance as a function of the frequency at several distances for two, three, and four coils as shown in Figure 6. The parameters of all arranged coils are identical with Coil_a in Table 1. The ideal mutual inductance does not depend on the frequency, but the practical mutual inductances are proportional to frequency. In addition, as the number of coils increases, the mutual inductance increases further due to the parasitic effect. In these cases, the results of the proposed method are almost identical to the simulated value.

Figure 11 shows the simulated transfer impedance of each route for four coils, which is normalized by the transfer impedance of the direct route, Z_R−1. The transfer impedances of one- and two-hop indirect routes are presented in Figure 11a,b, respectively. As the number of hops increase, the transfer impedance decreases dramatically. In other words, the normalized transfer impedance of one-hop routes have a value smaller than 0.1, while those of two-hop routes have a very small value of less than 0.005. Thus, the two-hop route is almost unaffected as an element to obtain the practical mutual inductance in a four-coil array. This means that among²ⁿ⁻²routes for n coils, only the direct and one-hop routes are practically dominant, so the transfer impedance between the primary and the nth coil can be approximately expressed as

_Zn1≃∑k=1²ⁿ⁻²_ZR−k≃∑k=1n−1_ZR−k.

Finally, for four coils, the practical mutual inductances by all routes and the routes except two more hopping routes with respect to the frequency are as shown in Figure 12. Even if the two-hop routes are not taken into account, the calculated mutual inductances are in good agreement with the results of all the routes as well as the simulated results.

4. Conclusions The multi-coil WPT systems commonly use three or more coils to improve efficiency in HF bands as well as LF bands. These systems are basically required to estimate or measure the practical mutual inductance between coils. In this paper, the unilateral approach method is proposed to estimate the mutual inductance among multiple coils. With the proposed method, it is not only simple to configure the expression of formulas but also possible to estimate the practical mutual inductance through several routes for multiple coils since the assumption that only one voltage generated with respect to the coil pair was considered. To verify the validity of the proposed method, we coaxially placed several coils along the distance and confirmed the consistent results with simulation and measurement. In reality, it is difficult to analyze the mutual inductance in multiple coils using a public VNA that has only two ports. However, the proposed method can provide the mutual inductance using a public two-port VNA. Furthermore, it will be quite useful for the multi-coil WPT systems with a regular array of identical coils.

Generate image description

Figure 1. The equivalent circuit model: (a) Ideal two coils; and (b) real coils with parasitic capacitors.

Generate image description

Figure 2. The equivalent circuit model of two coils with Route 1.

Generate image description

Figure 3. The procedures to obtain the coil's parameter and the ideal mutual inductance between coils.

Generate image description

Figure 4. The equivalent circuit model: (a) Three coils with ideal mutual inductances; and (b) with two routes.

Generate image description

Figure 5. The equivalent circuit model: (a) Four coils with ideal mutual inductances; and (b) with four routes.

Generate image description

Figure 6. The experimental array of the coils: (a) Two coils; (b) three coils; and (c) four coils.

Generate image description

Figure 7. The experimental setup with a four-port vector network analyzer (VNA).

Generate image description

Figure 8. Practical mutual inductance between two coils with distance.

Generate image description

Figure 9. Practical mutual inductance between Coils 1 and 3 with distance: (a) 6.78 MHz; and (b) 13.56 MHz.

Generate image description

Figure 10. Practical mutual inductances between side coils.

Generate image description

Figure 11. Normalized transfer impedance in the array of four coils: (a) One-hop; and (b) two-hop routes.

Generate image description

Figure 12. Simulated and calculated results for the array of four coils.

Parameter	Coila	Coilb
Inductance (μH)	1.13	3.83
Parasitic cap. (pF)	8.22	15.17

Two Cases	Coil 1	Coil 2	Coil 3
Case#1	Coila	Coila	Coila
Case#2	Coila	Coilb	Coila

	6.78 MHz		13.56 MHz
	Mean (μH)	NRMSE (%)	Mean (μH)	NRMSE (%)
Two coils	−2.7 × 10−4	0.68	16.9 × 10−4	1.67
Case#1	−5.8 × 10−4	1.03	20.0 × 10−4	1.92
Case#2	−6.5 × 10−4	0.84	7.3 × 10−4	2.81

Author Contributions

Conceptualization, S.-J.J. and S.-H.L.; methodology, S.-J.J.; software, S.-J.J.; validation, S.-J.J. and S.-H.L.; formal analysis, S.-J.J. and D.-W.S.; investigation, S.-J.J. and S.-H.L.; resources, S.-H.L.; data curation, S.-J.J.; writing-original draft preparation, S.-J.J.; writing-review and editing, D.-W.S.; visualization, S.-J.J.; supervision, D.-W.S.; project administration, D.-W.S.; funding acquisition, D.-W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2018R1C1B6003854).

Conflicts of Interest

The authors declare no conflicts of interest.

References

1. Zhong, W.X.; Lee, C.K.; Hui, S.Y.R. General analysis on the use of tesla's resonators in domino forms for wireless power transfer. IEEE Trans. Ind. Electron. 2013, 60, 261-270.

2. Kim, J.H.; Park, B.C.; Lee, J.H. Optimum design of WPT relay system by controlling capacitance. Microw. Opt. Technol. Lett. 2014, 56, 1658-1661.

3. Sun, Y.; Liao, Z.J.; Ye, Z.H.; Tang, C.S.; Wang, P.Y. Determining the maximum power transfer points for MC-WPT systems with arbitrary number of coils. IEEE Trans Power Electron. 2018, 33, 9734-9743.

4. Seo, D.W. Design method of three-coil WPT system based on critical coupling condition. IEEE J. Emerg. Sel. Top. Power Electron. 2019.

5. Sun, T.; Xie, X.; Li, G.; Gu, Y.; Wang, Z. Indoor wireless power transfer using asymmetric directly-strong-coupling mechanism. Microw. Opt. Technol. Lett. 2013, 55, 250-253.

6. Xu, B.; Li, Y. Investigation of surface wave propagation along a multi-coil wireless power transfer system. Microw. Opt. Technol. Lett. 2016, 58, 2261-2265.

7. Jadidian, J.; Katabi, D. Magnetic MIMO: How to charge your phone in your pocket. In Proceedings of the 20th Annual International Conference on Mobile Computing and Networking (MOBICOM), Maui, HI, USA, 7-11 September 2014; pp. 495-506.

8. Shi, L.; Kabelac, Z.; Katabi, D.; Perreault, D. Wireless power hotspot that charges all of your devices. In Proceedings of the 21th Annual International Conference on Mobile Computing and Networking (MOBICOM), Paris, France, 7-11 September 2015; pp. 2-13.

9. Sritongon, C.; Wisestherrakul, P.; Hansupho, N.; Nutwong, S.; Sangswang, A.; Naetiladdanon, S.; Mujjalinvimut, E. Novel IPT multi-transmitter coils with increase misalignment tolerance and system Efficiency. In Proceedings of the 2018 IEEE International Symposium on Circuits and Systems (ISCAS), Florence, Italy, 27-30 May 2018; pp. 1-5.

10. Nguyen, M.Q.; Ta, K.; Dubey, S.; Chiao, J.C. Frequency modes in a MIMO wireless power transfer system. In Proceedings of the 2017 IEEE Asia Pacific Microwave Conference (APMC), Kuala Lumpar, Malaysia, 13-16 November 2017; pp. 146-149.

11. Kim, D.H.; Ahn, D. Maximum efficiency point tracking for multiple-transmitters wireless power transfer. IEEE Trans. Power Electron. 2019.

12. Fu, M.; Zhang, T.; Ma, C.; Zhu, X. Efficiency and optimal loads analysis for multiple-receiver wireless power transfer systems. IEEE Trans. Microw. Theory Tech. 2015, 63, 801-812.

13. Alley, C.L.; Atwood, K.W. Electronic Engineering; Wiley: Hoboken, NJ, USA, 1973; p. 199.

14. Piri, M.; Jaros, V.; Fricaldsky, M. Verification of a mutual inductance calculation between two helical coils. In Proceedings of the 2015 16th International Scientific Conference on Electric Power Engineering (EPE), Kouty nad Desnou, Czech Republic, 20-22 May 2015; pp. 712-717.

15. Wu, L.; Lu, K.; Xia, Y. Mutual inductance calculation of two coaxial solenoid coils with iron core. In Proceedings of the 2018 21st International Conference on Electrical Machines and System (ICEMS), Jeju, Korea, 7-10 October 2018; pp. 1804-1808.

16. Ariffin, A.; Inamori, M. Influence of mutual inductance measurement for high efficiency wireless power transmission. In Proceedings of the 2018 7th International Conference on Renewable Energy Research and Applications (ICRERA), Paris, France, 14-17 October 2018; pp. 949-953.

17. Zhang, X.; Meng, H.; Wei, B.; Wang, S.; Yang, Q. Mutual inductance calculation for coils with misalignment in wireless power transfer. IET J. Eng. 2019, 16, 1041-1044.

18. Rakhymbay, A.; Khamitov, A.; Bagheri, M.; Alimkhanuly, B.; Lu, M.; Phung, T. Precise analysis on mutual inductance variation in dynamic wireless charging of electric vehicle. Energies 2018, 11, 624.

19. Hackl, S.; Lanschutzer, C.; Raggam, P.; Randeu, W.L. A Novel Method for Determining the Mutual Inductance for 13.56MHz RFID Systems. In Proceedings of the 2008 6th International Symposium on Communication Systems, Networks and Digital Signal Processing, Graz, Austria, 25-25 July 2008; pp. 297-300.

20. Pozar, D.M. Microwave Engineering, 4th ed.; Wiley: Hoboken, NJ, USA, 2011; p. 181.