Analysis and Optimization of Three-Resonator Wireless Power Transfer System for Predetermined-Goals Wireless Power Transmission

Abstract: Three-resonator wireless power transfer (WPT) systems have been proposed to improve the power transfer efficiency (PTE) and power delivered to the load (PDL) in recent years. However, analysis formulas of a three-resonator WPT system are complicated, and the parameters for clarifying the transfer characteristics of this system are difficult to extract. In this paper, concise formulas for analyzing PTE and PDL of the three-resonator system are derived by introducing three factors. Diagram discriminance based on the derived formulas is proposed to obtain the frequency splitting criterions of PTE and PDL in this system. Further, at the transfer distances, where the PTE and PDL are low at original frequency due to frequency splitting phenomenon, the two predetermined-goals of maximizing PTE and PDL are achieved by optimizing coupling strength between the three resonators. The third predetermined-goal of obtaining a constant amount of PDL transfer at maximum PTE is also implemented based on basic algorithms in numerical software. Finally, Simulation and measurement results verify the correctness of analyzing the transfer characteristics of three-resonator WPT system using the presented concise formulas and discriminance. Moreover, effectiveness of realizing the three predetermined-goals via the proposed optimization method is confirmed with experiments.


Introduction
Several means of wireless power transfer (WPT), including far-field, near-field, capacitive and magnetic coupling, have attracted much attention in recent years [1][2][3][4][5][6].In particular, the way of magnetic coupling WPT (MCWPT) has gained rapid development for its high efficiency and strong robustness.To extend the transmission distance of traditional two-resonator MCWPT system further, a method of inserting intermediate resonators (IX) between the transmitter (TX) and the receiver (RX) to form three-resonator system has been proposed recently [7][8][9][10][11][12][13][14][15][16].It has proved in [7] that three-resonator WPT system with IX close to TX or RX can achieve a relatively high power transfer efficiency (PTE) at original resonate frequency compared two-resonator counterpart.Ahn and Hong [8] pointed out that power transmission can be maximized at original resonant frequency via inserting a resonator in traditional two-resonator system.However, the conclusion that power transmission is always maximized at original resonant frequency in a three-resonator system is not entirely accurate, which will be proved in this paper.Although the analysis in [9] has proven the effectiveness of improving energy exchange rate in three-resonator system, the PTE and power delivered to the load (PDL) linking to the energy exchange rate was not quantitatively analyzed.Zhang et al. [10] has indicated that three-resonator WPT system can extend the power transfer distance, compared with traditional two-resonator WPT system.However, the acquisition of the optimal inserting position of IX for maximal PTE lacked theoretical analysis in this reference.Kim et al. [11] analyzed the efficiency characteristics of a WPT system with an IX which is coaxially and perpendicularly set to TX (or RX).In [12], due to strong coupling between TX and IX, a small current in TX can be obtained in a three-resonator system, which results in small loss in source resistance and improves PTE in such system.However, all the analyses in [10][11][12] were limited to original resonant frequency, which cannot reveal the relations between optimal transfer characteristics and system's working frequency.Although resonators with different resonant frequencies have been used in [13] to obtain a constant output current for various load resistor, further works which provide in-depth understanding of power transfer characteristics and guidelines for obtaining predetermined-goals WPT are needed in three-resonator system.
In order to maximize the PTE, optimizing the geometries of coils of TX and RX has been proposed [17][18][19].However, the optimized coils with the proper size, turns and conductor width for the maximizing PTE in these references were only effective for a certain transfer distance.A constant amount of PDL need to be transferred at maximum PTE for practical application in a WPT system [20][21][22].However, on-chip frequency-tracking oscillation class-D PA in [20] and additional compensation capacitors in [21] have increased the cost and volume of WPT system.
In this work, a detailed analysis of power transfer characteristics of three-resonator WPT versus working frequency and the mutual inductances between the three resonators is presented using circuit theory.Furthermore, a design method by adjusting rotation angles of TX and RX, corresponding to optimizing coupling strength between IX and TX (or RX), is used to achieve three predetermined-goals of maximizing PTE transfer, maximizing PDL transfer and a constant PDL transfer at maximum PTE at original frequency.Optimizing coupling coefficients between the resonators instead of optimizing coils' structure and of using additional compensation circuit to implement frequency tracking method can reduce the system cost and complexity to realize the three goals.The outline of this paper is as follows.Section 2 derives analytic formulas of the PTE and PDL for three-resonator system and analyzes the characteristics of PTE and PDL of a case based on practical resonators using diagram discriminance.Section 3 presents the optimization schemes for the three goals.Experimental verification of the analysis results and the effectiveness of proposed optimizing method are given in Section 4. Finally, Section 5 briefly summarizes the contributions of this work.

System Analysis
Figure 1a shows the simplified circuit of a magnetic coupling three-resonator WPT system.To concisely derive the analytical expressions of this WPT system, three parameters are introduced: inductance scaling factor, α 21 and α 31 (α 21 > 0 and α 31 > 0), expressing as α 2  21 " L 2 {L 1 and α 2 31 " L 3 {L 1 ; resistance scaling factor, β 21 and β 31 (β 21 > 0 and β 31 > 0), expressing as β 2  21 " R 2 {R 1 S and β 2 31 " R 1 L {R 1 S , where R 1 L " R 3 `RL and R 1 S " R 1 `RS ; generalized mistuning factor ω v " Q S0 pω{ω 0 ´ω0 {ωq, where Q S0 " ω 0 L 1 {R 1  S is loaded quality factor of lossy TX at the resonant frequency ω 0 .The currents I TX , I IX and I RX in TX, IX and RX can be solved by KVL via the following matrix equation: where Q S " ωL 1 {R 1 S is the loaded quality factor of lossy TX at working frequency ω, and V S is root mean square (RMS) amplitude of the sinusoidal voltage source.Since the cross coupling between TX Energies 2016, 9, 274 3 of 20 and RX are usually negligible under the conditions of k TR << k TI and k TR << k IR , the currents I TX , I IX and I RX are calculated from Equation (1) as follows: S rA 1 ´p1 `jω v q A 2 s (2a) S rA 1 ´p1 `jω v q A 2 s (2b) where: where:  

Transfer Characteristics of Power Transfer Efficiency
The PTE as a function of kTI, kIR and ω is defined as the ratio between the PDL and the total power supplied from the sources: where: and load resistance scaling factor βL

Transfer Characteristics of Power Transfer Efficiency
The PTE as a function of k TI , k IR and ω is defined as the ratio between the PDL and the total power supplied from the sources: where: and load resistance scaling factor S , ω 2 v and Q s are both functions of frequency ω.Obviously, to obtain analytical expression of extremal frequencies for the local extrema of F(ω, k TI , k IR ) (or PTE) is difficult.However, if the F(ω, k TI , k IR ) value calculated with a constant Q SCo approaches that with the variable Q S , the extremal frequencies expression can be expressed in terms of Q SCo .Using the resonator with parameters of R 1 = 5.1 Ω, L 1 =133 µH and resonant frequency f 0 = 1.96MHz as TX, IX and RX, Figure 2a shows the Q S versus normalized frequency ω/ω 0 .It varies from 144.5 to 176.6 with ω/ω 0 varying from 0.9 to 1.1, and Q S0 = 161 (C o=0 ) at ω = ω 0 .For comparing, Q S+ = 190 and Q S´= 140 are introduced to calculate F(ω, k TI , k IR ). Figure 2b shows function F(ω, k TI , k IR ) versus normalized frequency with variable Q S and the three constant Q S0 (C o=0 , + and ´) values at k TI = 0.0053, k IR = 0.0666, α , where x P p´8, 8q, which is used to investigate the local extrema of quadratic function with one variable in mathematics.f 1 (x) is the derivative of f(x) with respect to x: f 1 (x) = 2x + b e0 .x 0 = ´be0 /2 obtained by solving f 1 (x) = 0 is the point for minimum f(x).From Equation (3), it can be seen that PTE maximizes when f `ω2 v ˘minimizes.Properties of f(x) and frequency splitting situations of PTE are analyzed using diagram discriminance in Table 1.
) and ) at high (low) frequency range, respectively.The normalized frequencies of the local minimums of F(ω, kTI, kIR) for QS, QS0, QS+ and QS− are ) and ), respectively.It can be seen that NF  ) at high (low) frequency range respectively, therefore, we substitute QS0 for QS in Equation (3) to study the PTE characteristics.Frequency characteristics of PTE can be analyzed via ωv only.be and ce are rewritten as be0 and ce0 when QS is substituted by QS0.The denominator in Equation ( 3) is extracted to form a function which is used to investigate the local extrema of quadratic function with one variable in mathematics.fʹ(x) is the derivative of f(x) with respect to x: fʹ(x) = 2x + be0.x0 = −be0/2 obtained by solving fʹ(x) = 0 is the point for minimum f(x).From Equation (3), it can be seen that PTE maximizes when   Low Frquency Range Table 1.Properties and graph of f(x) = x 2 + be0x + ce0 for analyzing power transfer efficiency (PTE).
Graph and extrema for f(x), where Coupling region for PTE Non-frequency splitting region

Critical frequency splitting
Frequency splitting region Number of local maximum for PTE 1 -2 From Table 1, PTE maximizes at 2 v  = 0, namely ω = ω0, when x0 < 0, the corresponding coupling region of PTE is named non-frequency splitting region.Using 2 v 0 x   , the two splitting frequencies, ωe−H/L, for maximum PTE can be derived in frequency splitting region, where x0 > 0. Table 1.Properties and graph of f(x) = x 2 + b e0 x + c e0 for analyzing power transfer efficiency (PTE).
Graph and extrema for f(x), where x P p´8, 8q ) and ) at high (low) frequency range, respectively.The normalized frequencies of the local minimums of F(ω, kTI, kIR) for QS, QS0, QS+ and QS− are ) and ), respectively.It can be seen that NF  ) at high (low) frequency range respectively, therefore, we substitute QS0 for QS in Equation (3) to study the PTE characteristics.Frequency characteristics of PTE can be analyzed via ωv only.be and ce are rewritten as be0 and ce0 when QS is substituted by QS0.The denominator in Equation ( 3) is extracted to form a function which is used to investigate the local extrema of quadratic function with one variable in mathematics.fʹ(x) is the derivative of f(x) with respect to x: fʹ(x) = 2x + be0.x0 = −be0/2 obtained by solving fʹ(x) = 0 is the point for minimum f(x).From Equation (3), it can be seen that PTE maximizes when   Low Frquency Range Table 1.Properties and graph of f(x) = x 2 + be0x + ce0 for analyzing power transfer efficiency (PTE).
Graph and extrema for f(x), where Coupling region for PTE Non-frequency splitting region

Critical frequency splitting
Frequency splitting region Number of local maximum for PTE 1 -2 From Table 1, PTE maximizes at 2 v  = 0, namely ω = ω0, when x0 < 0, the corresponding coupling region of PTE is named non-frequency splitting region.Using 2 v 0 x   , the two splitting frequencies, ωe−H/L, for maximum PTE can be derived in frequency splitting region, where x0 > 0. From Table 1, PTE maximizes at ω 2 v = 0, namely ω = ω 0 , when x 0 < 0, the corresponding coupling region of PTE is named non-frequency splitting region.Using ω 2 v " x 0 , the two splitting frequencies, ω e´H/L , for maximum PTE can be derived in frequency splitting region, where x 0 > 0.

Locations of x
By substituting ω 2 v = 0 into Equation (3), The PTE at original frequency, PTE 0 , can be presented as: where: b p " β For the same reason as PTE in Equation ( 3), all Q S in Equation ( 6) can be substituted by Q S0 .b p and c p are rewritten as b p0 and c p0 when Q S is substituted by Q S0 .The expression under the square root sign is extracted to form a function g `ω2 v ˘in term of ω 2 v , g `ω2 , obviously ω 2 v ě 0. Substituting variable x for ω 2 v in g `ω2 v ˘, the expression is rearranged as a cubic function of x, g(x) = x 3 + b p0 x 2 + c p0 x + d p0 , where x P p´8, 8q, which is used to investigate the local extrema of cubic function with one variable in mathematics.g 1 (x) is the derivative of g(x) with respect to x: g 1 (x) = 3x 2 + 2b p0 x + c p0 and ∆ = (2b p0 ) 2 -4(3c p0 ), where factor ∆ is used to judge whether there are local extrema or not for g(x).
Two cases of ∆ > 0 and ∆ = 0 are listed in Table 2 for illuminating the characteristics of g(x) and frequency splitting situations of |I RX | (or PDL).x 1 and x 2 obtained by solving g'(x) = 0 are the stationary points for g(x) when ∆ ě 0: x 2 " ´´b p0 `bb 2 p0 ´3c p0 ¯{3 (7b) Graphs and extrema for g(x), where x P p´8, 8q Table 2. Properties and graphs of g(x) = x3 + bp0x2 + cp0x + dp0 for analyzing RX I (or power delivered to the load (PDL)).

Graphs and extrema for g(x)
, where

Frequency number of local maximum for
1 Critical coupling condition for judging the border between non-frequency splitting region and frequency splitting region with two splitting frequencies. 2Critical coupling condition for judging the border between frequency splitting region with two splitting frequencies and with three splitting frequencies. 3Critical coupling condition for judging the border between non-frequency splitting region and frequency splitting region with three splitting frequencies.

Graphs and extrema for g(x)
, where

Frequency number of local maximum for
1 Critical coupling condition for judging the border between non-frequency splitting region and frequency splitting region with two splitting frequencies. 2Critical coupling condition for judging the border between frequency splitting region with two splitting frequencies and with three splitting frequencies. 3Critical coupling condition for judging the border between non-frequency splitting region and frequency splitting region with three splitting frequencies.From Equation ( 6), it can be seen that |I RX | maximizes when g `ω2 v ˘minimizes.If ∆ > 0 and x 1 < x 2 < 0, |I RX | maximizes only at ω 2 v " 0, the frequency number of local maximums of |I RX | is one and its coupling region names non-frequency splitting region; if ∆ > 0 and x 1 < 0 < x 2 , |I RX | maximizes only at ω 2 v " x 2 , the frequency number of local maximums of |I RX | is two and its coupling region names frequency splitting region; if ∆ > 0 and 0 < x 1 < x 2 , |I RX | maximizes at ω 2 v " x 2 and ω 2 v " 0, the frequency number of local maximums of |I RX | is three and its coupling region names frequency splitting region; if ∆ = 0 and 0 ě x 1 = x 2 , |I RX | maximizes only at ω 2 v " 0, the frequency number of local maximums of |I RX | is one and its coupling region names non-frequency splitting region; if ∆ < 0, g'(x) = 0 has no solution and |I RX | always maximizes at ω 2 v " 0, the frequency number of local maximums of |I RX | is one and its coupling region names non-frequency splitting region.There are three kinds of critical point to define the frequency splitting regions of |I RX |: one occurs when ∆ > 0 and x 1 < x 2 = 0, it is the critical point for differentiating between non-frequency splitting region and frequency splitting region with two splitting frequencies; one occurs when ∆ > 0 and 0 = x 1 < x 2 , it is the critical point for differentiating between frequency splitting region with two splitting frequencies and frequency splitting region with three splitting frequencies; one occurs when ∆ = 0 and 0 < x 1 = x 2 , it is the critical point for differentiating between non-frequency splitting region and frequency splitting region with three splitting frequencies.The two splitting frequencies, ω p´H/L , in frequency splitting region deviating from ω 0 , can be obtained by solving ω 2 v " x 2 .

Coupling regions for |I
By substituting ω 2 v = 0 into Equation ( 6), The |I RX | at original frequency, |I RX | 0 , can be presented as: As for multi-resonator system, the expressions of PTE and PDL can be derived by reflected impedance method [14][15][16].Compared with these references, the characteristics of PTE and PDL with respect to non-original frequency, related to ω v , are considered in our study.Based on reflected impedance method, the expressions of PTE and PDL for n-resonator system are n ´1 and n order function of ω 2 v , respectively.

Theoretical Illumination with a Case Based on Practical Resonators
In this section, a case is presented to clarify the analysis in Sections 2.1 and 2.2.To reveal the transfer characteristics in terms of the transmission distance, the coupling coefficients calculated in the aforementioned procedure need to be transformed into transmission distance first.The mutual inductance equation proposed in [23] is revamped in this paper for two single layer coils with inclined axes (Figure 3) as follow: where: p ; g " ´m1 , ..., 0, ..., m 1 ; p " ´m2 , ..., 0, ..., m 2 H i , r i and N i are the height, radius and turn number of coil i (i = 1,2).D and S are the longitudinal and transverse distances between coils' centers, respectively.θ is the angle between coil planes.m 1 and m 2 are the subdivision number for the height of Coils 1 and 2. K(k) and E(k) are complete elliptic integrals of the first and second kind, respectively.
Energies 2016, 9, 274 8 of 19 Large values of m1 and m2 increase the calculation accuracy of mutual inductance.In this study, we calculate with m1 = m2 = 5 and S = 0 m.We use three identical resonators, which consist of identical coils and series lumped capacitors, to serve as the TX, IX and RX respectively.The parameters of the real resonator are as follows: the multi-strand Litz wire is used to construct the coils, its structures is 0.1 mm 80 strands   and overall outer diameter is and VS = 10 V, the distance between TX and RX, DTR, is fixed at 70 cm, and the distance between IX and TX, DTI, changes from 10 cm to 60 cm. Figure 4a shows the judgment conditions of xj (j = 0, 1 and 2) and  Large values of m 1 and m 2 increase the calculation accuracy of mutual inductance.In this study, we calculate with m 1 = m 2 = 5 and S = 0 m.
We use three identical resonators, which consist of identical coils and series lumped capacitors, to serve as the TX, IX and RX respectively.The parameters of the real resonator are as follows: the multi-strand Litz wire is used to construct the coils, its structures is ∅ 0.1 mm ˆ80 strands and overall outer diameter is ∅ 1.18 mm; the radius, height and turn number of the coil are 15.75 cm, 1.5 cm and 11, respectively.The measured value of self-inductance is 133 µH, total parasitic resistance is 5.1 Ω, and original resonant frequency of the resonator is 1.96 MHz.The parameters for the sample case are: 1 Ω, R L = 45.9Ω) and V S = 10 V, the distance between TX and RX, D TR , is fixed at 70 cm, and the distance between IX and TX, D TI , changes from 10 cm to 60 cm.
Figure 4a shows the judgment conditions of x j (j = 0, 1 and 2) and ∆ for frequency splitting for PTE and |I RX | of the sample case.Combining x 0 in this figure with frequency splitting discriminance in Table 1, it can be found that PTE is in the non-frequency splitting region, critical frequency splitting position and frequency splitting region when D TI <, = and > 0.408 m, respectively.Ignoring the cross coupling between TX and RX, the values of PTE versus distance and frequency are shown in Figure 5a.The frequencies for local maximum (FLM) of PTE under different D TI , shown as black lines, are achieved by Equation ( 4).The frequencies for maximum (FM) of PTE are obtained by searching the Energies 2016, 9, 274 9 of 20 maximum in global frequency domain at each D TI , shown as white dotted lines.It can be seen that the critical frequency splitting distance is at 0.408 m, which validates the analysis result of judging the position for frequency splitting.The FM matches well with the FLM in non-frequency splitting region where D TI < 0.408 m, and is very close to the FLM at high splitting frequency in frequency splitting region where D TI > 0.408 m.In the greater degree of frequency splitting region, i.e., the region of splitting frequencies keeping farther away from each other, the FM slightly deviates from FLM, seeing the distribution of white dotted line and black line at larger D TI in Figure 5a.The reason is that the splitting frequencies are calculated under the constant Q S0 .The values of PTE versus distance and frequency with considering cross coupling are shown in Figure 5b.The FM is off FLM slightly to low frequency in non-frequency splitting region because of the effect of cross coupling.In frequency splitting region, the FM is close to the FLM at low splitting frequency under the effect of cross coupling.Comparing Figure 5a,b, the transfer distances for maximal PTE for with and without considering cross coupling are the same.All these show that the analytical solutions obtained by simplified model, namely model ignoring cross coupling, can well show the transfer characteristics of practical complete model, namely model considering cross coupling.Combining x 1 , x 2 and ∆ in Figure 4 and frequency splitting discriminance in Table 2 with the analysis in Section 2.2, we can see: |I RX | is in the non-frequency splitting region when 0.373785 < D TI < 0.373826 m (∆ > 0, x 1 < x 2 < 0), D TI = 0.373826 m (∆ = 0, x 1 = x 2 < 0) and 0.373826 < D TI < 0.536 m (∆ < 0); |I RX | is in the frequency splitting region with two splitting frequencies when 0.212 < D TI < 0.373785 m (∆ > 0, x 1 < 0 < x 2 ); |I RX | is in the frequency splitting region with three splitting frequencies when D TI < 0.212 m and D TI > 0.536 m (∆ > 0, 0 < x 1 < x 2 ).D TI = 0.373785 m is the critical coupling distance between non-frequency region and frequency region with two splitting frequencies (∆ > 0, x 1 < x 2 = 0); D TI = 0.212 m is the critical coupling distance between frequency region with two splitting frequencies and frequency region with three splitting frequencies (∆ > 0, 0 = x 1 < x 2 ); D TI = 0.536 m is the critical coupling distance between non-frequency region and frequency region with three splitting frequencies (∆ = 0, 0 < x 1 = x 2 ).Ignoring the cross coupling between TX and RX, the transfer characteristics of |I RX | versus distance and frequency are shown in Figure 5c.Plotting |I RX | instead of PDL in this figure is that the color gradation of the former is much clearer than the latter in two-dimensional color image.The distributions of |I RX | well verify the analysis of frequency splitting situation.At f 0 , The FM of |I RX | matches well with FLM which is obtained by Equation ( 8) in non-frequency splitting region (0.373785 < D TI < 0.536 m) and frequency splitting region with three splitting frequencies (D TI > 0.536 m), and it is very close to the FLM at high splitting frequency in frequency splitting region with two and three splitting frequencies (0.212 < D TI < 0.373785 m and D TI < 0.212 m). Figure 5d shows the values of |I RX | considering the cross coupling.One major difference from Figure 5c ignoring the cross coupling is that FM in Figure 5d is close to the FLM at low splitting frequency in the two kinds of frequency splitting regions (0.212 < D TI < 0.373785 m and D TI < 0.212 m).
Comparing the simulated values of simplified model in Figure 5c and complete model in Figure 5d, the maximal |I RX | of them and the transfer distance for the maximal |I RX | of them are basically the same.
By this case and the theoretical analysis in this section, we can found that PTE of three-resonator WPT system has at most two splitting frequencies and PDL (|I RX |) has at most three splitting frequencies.

Optimization for Predetermined-Goals Wireless Power Transmission
Via the analysis in Section 2, we can see that frequency splitting for PTE and PDL takes place when IX gets close to TX or RX, which cause the PTE and PDL decrease at original resonant frequency, f 0 .In this section, we adjust k TI and k IR to obtain three kinds of predetermined-goal WPT at f 0 for different inserted distance, D TI : maximizing PTE transfer; maximizing PDL transfer; obtaining a constant amount of PDL transfer at maximum PTE.MATLAB optimization toolbox is used to determine the k TI and k IR for the three kinds of optimization goals.The frequency tracking method, used in [24,25] to maximize the PTE and PDL, is also presented in this section for comparing.

Setting the Optimization Goals
As mentioned in Section 2, analytical Equations ( 3) and ( 6) of simplified model are competent to reveal the real transfer characteristics of complete model when cross coupling between TX and RX can be ignored (k TR << k TI and k TR << k IR ).Equations ( 5) and ( 9) as functions of tuning coupling coefficient vector k, k = [k TI , k IR ], are used to define the optimization functions for the three kinds of predetermined-goal at f 0 .There is a maximum value for k, i.e., k max = [k TI,max , k IR,max ], for a certain relative positions of TX, IX and RX, when the axes of them are arranged collinearly, namely Axe TX , Axe IX and Axe RX in Figure 1b are collinear, corresponding to θ TX = θ RX = 0 ˝.By using solver f mincon in MATLAB, the objective optimization function and constrains of the three predetermined-goals are denoted as following: Maximizing PTE transfer: Maximizing PDL transfer: Obtaining a constant amount of PDL transfer at maximum PTE: where PDL C is a constant amount of PDL, according to the requirement of charged devices, for the third optimal goal.

Numerical Results
The case, analyzed in Section 2.3, is used for the numerical calculation to clarify the feasibility of the optimization method for the first two predetermined-goals.The PTE curves of collinear axes arrangement of the three resonators at f 0 , maximizing PTE using optimizing k at f 0 , and maximizing PTE using frequency tracking method are shown in Figure 6a.The maximal coupling coefficient vector, k max = [k TI,max , k IR,max ], and optimal coupling coefficient vector, k = [k TI , k IR ], for maximizing PTE for each inserting distance, D TI , are plotted in Figure 6b.It can be seen that the values of PTE obtained by optimizing k are equal to those of collinear axes arrangement scenario and frequency tracking method when D TI < 0.37 m.This can be explained by the fact that the local maximal PTE only takes place at f 0 when IX gets close to TX (Figure 5a,b).PTE values of frequency tracking method are higher than those of collinear axes arrangement scenario when IX moves over 4.1 m, which is the frequency splitting distance for PTE, actually the accurate splitting distance shown in Figure 4 is 4.08 m.For the inserting distances of IX close to RX, optimizing k for maximal PTE method obtains higher PTE than frequency tracking method and collinear axes arrangement scenario, for instance the former can improve the PTE by up to 10% and 18% comparing with the latter two at D TI = 0.5 m, respectively.The PTE values of maximizing PDL via optimizing k, which are always lower than those of maximizing PTE at each D TI , are also compared in Figure 6a.The optimal k TI for maximizing PTE always maintains the maximum value, k TI,max , at whole D TI , while the optimal k IR gradually deviates from k IR,max to lower value starting at distance D TI = 0.37 m.The reason for this optimal k is that the strength of k TI,max is not enough high to cause frequency splitting when IX moves to TX in collinear axes arrangement scenario while the frequency splitting takes place with increasing k IR,max when IX moves to RX.
Energies 2016, 9, 274 11 of 19 method are higher than those of collinear axes arrangement scenario when IX moves over 4.1 m, which is the frequency splitting distance for PTE, actually the accurate splitting distance shown in Figure 4 is 4.08 m.For the inserting distances of IX close to RX, optimizing k for maximal PTE method obtains higher PTE than frequency tracking method and collinear axes arrangement scenario, for instance the former can improve the PTE by up to 10% and 18% comparing with the latter two at DTI = 0.5 m, respectively.The PTE values of maximizing PDL via optimizing k, which are always lower than those of maximizing PTE at each DTI, are also compared in Figure 6a.The optimal kTI for maximizing PTE always maintains the maximum value, kTI,max, at whole DTI, while the optimal kIR gradually deviates from kIR,max to lower value starting at distance DTI = 0.37 m.The reason for this optimal k is that the strength of kTI,max is not enough high to cause frequency splitting when IX moves to TX in collinear axes arrangement scenario while the frequency splitting takes place with increasing kIR,max when IX moves to RX.In Figure 7a, PDL of maximizing PDL goal by optimizing k at f0 is compared with PDLs of collinear axes arrangement scenario at f0, maximizing PTE goal by optimizing k at f0 and maximizing PDL goal by frequency tracking method.As in non-frequency splitting region, 0.373785 < DTI < 0.536 m (Figure 4a,b), and near RX frequency splitting region, 0.536 < DTI < 0.6 m (Figure 4a), the maximal PDL remains at f0 (Figure 5c,d), the PDLs obtained via maximizing PDL goal by frequency tracking In Figure 7a, PDL of maximizing PDL goal by optimizing k at f 0 is compared with PDLs of collinear axes arrangement scenario at f 0 , maximizing PTE goal by optimizing k at f 0 and maximizing PDL goal by frequency tracking method.As in non-frequency splitting region, 0.373785 < D TI < 0.536 m (Figure 4a,b), and near RX frequency splitting region, 0.536 < D TI < 0.6 m (Figure 4a), the maximal PDL remains at f 0 (Figure 5c,d), the PDLs obtained via maximizing PDL goal by frequency tracking method are the same with those of collinear axes arrangement scenario when 0.373785 < D TI < 0.6 m (Figure 7a).While the PDLs of maximizing PDL goal by frequency tracking method are higher than those of collinear axes arrangement scenario in near TX frequency splitting region, 0.1 < D TI < 0.373785 m (Figure 4a), due to the maximal PDLs occurring at low splitting frequency in this region (Figure 5d).Compared with collinear axes arrangement scenario and maximizing PDL goal by frequency tracking method, the maximizing PDL goal by optimizing k can improve PDLs both below and above the position of D TI = 0.41 m.The corresponding tuning k for maximizing PDL goal are compared with the k max for collinear axes arrangement scenario in Figure 7b.It can be found that decreasing k TI properly under D TI < 0.41 m and decreasing k IR properly under D TI > 0.41 m can improve the PDL at f 0 .The PDLs of maximizing PTE goal by optimizing k, which are comparatively plotted in Figure 7a, are always lower than those of maximizing PDL goal by optimizing k.The third predetermined-goal of obtaining a constant amount of PDL transfer at maximum PTE is instanced under transfer distances D TR = 0.5, 0.7 and 0.9 m respectively.The source and load resistance values and resonators used for this optimizing goal are the same with the case for the first two predetermined-goals; the constant PDL C is set as 0.5 W in this goal.According to the optimization function and constrains in Equation ( 15), PTEs of this goal by optimizing k are compared with those of maximizing PTE goal by optimizing k and collinear axes arrangement scenario in Figure 8a; and PDLs of this goal by optimizing k are compared with those of maximizing PDL goal by optimizing k and collinear axes arrangement scenario in Figure 8b.It is clear from Figure 8b that the constant PDL = 0.5 W can be obtained with 0.1 m ď D TI ď 0.4 m, 0.16 m ď D TI ď 0.6 m and 0.36 m ď D TI ď 0.72 m for D TR = 0.5, 0.7 and 0.9 m respectively.It is also worth pointing out that no matter how the k is tuned, PDL cannot achieve the power of 0.5 W when D TI < 0.16 m for D TR = 0.7 m and D TI < 0.36 m and D TI > 0.72 m for D TR = 0.9 m.The lower bound of D TI = 0.16 m for 0.5 W in D TR = 0.7 m case and lower (upper) bound of D TI = 0.36 (0.72) m for 0.5 W in D TR = 0.9 m case are determined by the maximizing PDL goal by optimizing k.For D TR = 0.5, 0.7 and 0.9 m cases, the PTEs are maximized at D TI = 0.21, 0.31 and 0.42 m, respectively, where the corresponding collinear axes arrangement scenarios achieve the PDL values of 0.5 W. For D TR = 0.5, 0.7 and 0.9 m cases, the optimal values of k for obtaining a constant PDL transfer at maximum PTE, comparing with those of collinear axes arrangement scenario, are shown in Figure 8c.

Experimental Verification
In this section, we discuss the experimental verification of the theoretical analysis and calculation presented in Section 2 and the numerical optimization results of the three predetermined-goals shown in Section 3.2.

Experimental Setup
The power transfer characteristics of the three-resonator system developed in this work are tested with the aid of the two-port VNA, as shown in Figure 9a.The part number of this VNA is E5071C, Agilent ENA series.The angle vector, θ = [θ TX , θ RX ] shown in Figures 1b and 9a, is tuned to achieve the optimal k for the three predetermined-goals WPT.According to the measured S-parameters based on port impedance z 0 = 50 Ω of the network analyzer, the PDL (or I RX ) and PTE of the WPT network with arbitrary R S and R L values can be determined by transforming calculation method [22,26], and the block diagram for the procedures is presented in Figure 9b.
with the aid of the two-port VNA, as shown in Figure 9a.The part number of this VNA is E5071C, Agilent ENA series.The angle vector, θ = [θTX, θRX] shown in Figures 1b and 9a, is tuned to achieve the optimal k for the three predetermined-goals WPT.According to the measured S-parameters based on port impedance z0 = 50 Ω of the network analyzer, the PDL (or IRX) and PTE of the WPT network with arbitrary RS and RL values can be determined by transforming calculation method [22,26], and the block diagram for the procedures is presented in Figure 9b.

Verification of the Analysis of Three-Resonator Wireless Power Transfer System
The case in Section 2.3 is measured to verify the analysis of the energy transfer characteristics in three-resonator WPT system.Theoretical results with and without considering the effects of the cross coupling, compared with measured counterparts at four representative inserted distances for PTE and PDL, are shown in Figure 10.From Figure 10a, the measured datum at DTI = 0.2, 0.34 and 0.4 m, which are located in the theoretical non-frequency splitting region (DTI < 0.41 m) of PTE, each of them has only one local maximum value of PTE; and at DTI = 0.5 m located in the theoretical frequency splitting region (DTI > 0.41 m), the measured curve has two local maximum values of PTE.Similarly in Figure 10b, the measured datum at DTI = 0.15 and 0.55 m located in the theoretical frequency splitting region with three splitting frequencies of PDL (DTI < 0.212 m and > 0.536 m) and at DTI = 0.3 m located in the theoretical frequency splitting region with two splitting frequencies of PDL (0.212 m < DTI < 0.373785 m) have three and two local maximum values of PDL, respectively.Transmission characteristic of PDL at DTI = 0.3 m indicates that the verdict of achieving maximum power transmission always at original resonant frequency in three-resonant system presented in [8] is inaccurate.At the distance of DTI = 0.41 m in non-frequency splitting region, only one maximal measured value of PDL appears.All of these demonstrate that two splitting frequencies arise at most for local maximum PTE in a three-resonator WPT system, and three or two splitting frequencies would appear for local maximum PDL when relative positions of TX, IX and RX are changed.

Verification of the Analysis of Three-Resonator Wireless Power Transfer System
The case in Section 2.3 is measured to verify the analysis of the energy transfer characteristics in three-resonator WPT system.Theoretical results with and without considering the effects of the cross coupling, compared with measured counterparts at four representative inserted distances for PTE and PDL, are shown in Figure 10.From Figure 10a, the measured datum at D TI = 0.2, 0.34 and 0.4 m, which are located in the theoretical non-frequency splitting region (D TI < 0.41 m) of PTE, each of them has only one local maximum value of PTE; and at D TI = 0.5 m located in the theoretical frequency splitting region (D TI > 0.41 m), the measured curve has two local maximum values of PTE.Similarly in Figure 10b, the measured datum at D TI = 0.15 and 0.55 m located in the theoretical frequency splitting region with three splitting frequencies of PDL (D TI < 0.212 m and > 0.536 m) and at D TI = 0.3 m located in the theoretical frequency splitting region with two splitting frequencies of PDL (0.212 m < D TI < 0.373785 m) have three and two local maximum values of PDL, respectively.Transmission characteristic of PDL at D TI = 0.3 m indicates that the verdict of achieving maximum power transmission always at original resonant frequency in three-resonant system presented in [8] is inaccurate.At the distance of D TI = 0.41 m in non-frequency splitting region, only one maximal measured value of PDL appears.All of these demonstrate that two splitting frequencies arise at most for local maximum PTE in a three-resonator WPT system, and three or two splitting frequencies would appear for local maximum PDL when relative positions of TX, IX and RX are changed.

Verification of the Feasibility for Predetermined-Goals Wireless Power Transfer in Three-Resonator System
To realize the three kinds of predetermined-goal WPT in practical measurement, the optimal values of k for these goals achieved in the aforementioned procedure are transformed into tuning angle θ = [θ TX , θ RX ].Setting S = 0 m in Figure 3, the coupling coefficient, k, between two coils versus rotation angle of Coil 2, θ, at different transfer distances, D, are calculated using Equations ( 10) and (11) and plotted in Figure 11.Using the case in Section 2.3, at DTI = 0.45 m, The measured and calculated values of PTE of collinear axes arrangement scenario and maximizing PTE goal by tuning angle vector to θ = [θTX, θRX] = [0°, 75.4°], which is the optimal angle for maximum PTE according to Figures 6b and 11, are shown in Figure 12a.It can be seen that measured results agree with the calculated ones of system model with considering the cross coupling between TX and RX, and maximal PTE at original resonant frequency 1.96 MHz of maximizing PTE goal by tuning rotation angle of TX and RX is higher than the counterpart at low splitting frequency of collinear axes arrangement scenario using frequency tracking method.Using the same case, similar comparison for PDL of collinear axes arrangement scenario and maximizing PDL goal by tuning angle vector to θ = [θTX, θRX] = [79.7°,0°] at DTI = 0.3 m are included in Figure 12b.We can also find that maximal PDL of maximizing PDL goal is higher than that of collinear axes arrangement scenario.Using the case in Section 2.3, at D TI = 0.45 m, The measured and calculated values of PTE of collinear axes arrangement scenario and maximizing PTE goal by tuning angle vector to θ = [θ TX , θ RX ] = [0 ˝, 75.4 ˝], which is the optimal angle for maximum PTE according to Figures 6b and  11, are shown in Figure 12a.It can be seen that measured results agree with the calculated ones of system model with considering the cross coupling between TX and RX, and maximal PTE at original resonant frequency 1.96 MHz of maximizing PTE goal by tuning rotation angle of TX and RX is higher than the counterpart at low splitting frequency of collinear axes arrangement scenario using frequency tracking method.Using the same case, similar comparison for PDL of collinear axes arrangement scenario and maximizing PDL goal by tuning angle vector to θ = [θ TX , θ RX ] = [79.7 ˝, 0 ˝] at D TI = 0.3 m are included in Figure 12b.We can also find that maximal PDL of maximizing PDL goal is higher than that of collinear axes arrangement scenario.Using the case in Section 2.3, at DTI = 0.45 m, The measured and calculated values of PTE of collinear axes arrangement scenario and maximizing PTE goal by tuning angle vector to θ = [θTX, θRX] = [0°, 75.4°], which is the optimal angle for maximum PTE according to Figures 6b and 11, are shown in Figure 12a.It can be seen that measured results agree with the calculated ones of system model with considering the cross coupling between TX and RX, and maximal PTE at original resonant frequency 1.96 MHz of maximizing PTE goal by tuning rotation angle of TX and RX is higher than the counterpart at low splitting frequency of collinear axes arrangement scenario using frequency tracking method.Using the same case, similar comparison for PDL of collinear axes arrangement scenario and maximizing PDL goal by tuning angle vector to θ = [θTX, θRX] = [79.7°,0°] at DTI = 0.3 m are included in Figure 12b.We can also find that maximal PDL of maximizing PDL goal is higher than that of collinear axes arrangement scenario.11), which results in inexact tuning angles and makes the measurement results inaccurate.The frequencies for the local maximum of model with cross coupling are consistent with the counterparts of model ignoring cross coupling at D TI = 0.2, 0.4 and 0.55 m.The reason is because small cross coupling between TX and RX is achieved by tuning rotation angles of TX and RX.By combining the analysis results of the D TR = 0.7 m case in Figure 8 with the measured results in Figure 13, it can be verified that PTE obtains its maximum value with preset PDL of 0.5 W at D TI = 0.31 m, where optimal angle vector θ = [0 ˝, 0 ˝].
Energies 2016, 9, 274 16 of 19 sensitive to θ (see curves of D = 0.2 or 0.15 m in Figure 11), which results in inexact tuning angles and makes the measurement results inaccurate.The frequencies for the local maximum of model with cross coupling are consistent with the counterparts of model ignoring cross coupling at DTI = 0.2, 0.4 and 0.55 m.The reason is because small cross coupling between TX and RX is achieved by tuning rotation angles of TX and RX.By combining the analysis results of the DTR = 0.7 m case in Figure 8 with the measured results in Figure 13, it can be verified that PTE obtains its maximum value with preset PDL of 0.5 W at DTI = 0.31 m, where optimal angle vector θ = [0°, 0°].According to Figure 8 and the analysis for this figure and the verification result in Figure 13, the PTE maximizes at an optimal inserted distance DTI when optimal angle vector θ = [0°, 0°] for a certain transfer distance DTR.At original resonant frequency of 1.96 MHz, calculated and measured values of PTE and PDL of the third optimizing goal with two preset constant PDLC of 0.5 W and 1 W, obtained with optimal DTI and θ = [0°, 0°], are shown in Figure 14b,c.In the case of PDLC = 0.5 W, the DTR is effective to obtain this predetermined PDL for up to a distance of 1.16 m; however, when PDLC = 1 W, the effective working distance is maximized at 0.99 m.This result indicates that the effective transmission distance DTR decreases as PDLC increases.The ratio of optimal DTI to DTR for the two constant PDLC are plotted in Figure 14a.As can be seen from this figures, if transfer distances between TX and RX are at closer ranges, such as DTR < 0.76 m and 1.07 m for PDLC = 0.5 W and 1 W respectively, the optimal inserted positions of IX should be set close to TX, i.e., DTI/DTR < 0.5; if not, the optimal inserted positions of IX are placed closed to RX.From Figure 14b,c, PTEs of model with cross coupling are nearly consistent with those of model ignoring cross coupling, and PDLs of the two models are increasingly differ with decreasing DTR because of the increasing effects of cross coupling between TX and RX.However, it is effective using the proposed optimal method to achieve the preset constant PDLs at maximum PTEs starting from DTR = 0.4 m and 0.5 m for PDLC = 0.5 W and 1 W respectively, for the maximum offsets below PDLC = 0.5 W at DTR = 0.4 m and PDLC = 1 W at DTR = 0.5 m are 12.9% and 9.21% respectively.The effects of cross coupling on energy transfer can be addressed through power flow analysis [22].The power flow from TX and IX to Rx can be calculated  According to Figure 8 and the analysis for this figure and the verification result in Figure 13, the PTE maximizes at an optimal inserted distance D TI when optimal angle vector θ = [0 ˝, 0 ˝] for a certain transfer distance D TR .At original resonant frequency of 1.96 MHz, calculated and measured values of PTE and PDL of the third optimizing goal with two preset constant PDL C of 0.5 W and 1 W, obtained with optimal D TI and θ = [0 ˝, 0 ˝], are shown in Figure 14b,c.In the case of PDL C = 0.5 W, the D TR is effective to obtain this predetermined PDL for up to a distance of 1.16 m; however, when PDL C = 1 W, the effective working distance is maximized at 0.99 m.This result indicates that the effective transmission distance D TR decreases as PDL C increases.The ratio of optimal D TI to D TR for the two constant PDL C are plotted in Figure 14a.As can be seen from this figures, if transfer distances between TX and RX are at closer ranges, such as D TR < 0.76 m and 1.07 m for PDL C = 0.5 W and 1 W respectively, the optimal inserted positions of IX should be set close to TX, i.e., D TI /D TR < 0.5; if not, the optimal inserted positions of IX are placed closed to RX.From Figure 14b,c, PTEs of model with cross coupling are nearly consistent with those of model ignoring cross coupling, and PDLs of the two models are increasingly differ with decreasing D TR because of the increasing effects of cross coupling between TX and RX.However, it is effective using the proposed optimal method to achieve the preset constant PDLs at maximum PTEs starting from D TR = 0.4 m and 0.5 m for PDL C = 0.5 W and 1 W respectively, for the maximum offsets below PDL C = 0.5 W at D TR = 0.4 m and PDL C = 1 W at D TR = 0.5 m are 12.9% and 9.21% respectively.The effects of cross coupling on energy transfer can be addressed through power flow analysis [22].The power flow from TX and IX to Rx can be calculated by P TR " Re `´jωM TR I TX I RX ˘and P IR " Re `´jωM IR I IX I RX ˘, respectively.Power loss in resonator of RX is P Loss_RX " ´|I RX | 2 R 3 .In Figure 15, the curves of P TR , P IR , P Loss_RX and PDL (PDL = P TR + P IR + P Loss_RX ) are plotted versus D TR when PDL C = 1 W and 0.5 W, respectively.It can be observed that the power delivered from IX to RX (P IR ), which is affected gradually by cross coupling with decreasing D TR , mainly determines the volume of PDL.

Conclusions
In this paper, we have introduced three parameters to derive concise formulations of the PTE and PDL based on coupled circuit model of three-resonator WPT system.Frequency splitting

Conclusions
In this paper, we have introduced three parameters to derive concise formulations of the PTE and PDL based on coupled circuit model of three-resonator WPT system.Frequency splitting

Conclusions
In this paper, we have introduced three parameters to derive concise formulations of the PTE and PDL based on coupled circuit model of three-resonator WPT system.Frequency splitting characteristics and calculating formulas for splitting frequencies of PTE and PDL are achieved by diagram discriminance.It has been proven via theoretical analysis and experimental measurement that there are non-frequency splitting region and frequency splitting region with two splitting frequencies for PTE, and non-frequency splitting region and frequency splitting region with two and three splitting frequencies for PDL.By transfer characteristic analysis of simplified circuit model and tuning rotation angles of TX and RX, an optimizing method using basic optimization algorithm is presented to achieve three predetermined-goals WPT of maximizing PTE transfer, maximizing PDL transfer and obtaining a constant amount of PDL transfer at maximum PTE at original resonant frequency.Comparing with frequency tracking method, the first two goals of maximizing PTE or PDL based on the proposed optimizing method can achieve higher PTE or PDL at different D TI .Experimental results have demonstrated the effectiveness of the proposed method.For the third optimizing goal, theoretical analysis and experimental results show that the maximum PTE always appears at tuning angles of θ = [θ TX , θ RX ] = [0 ˝, 0 ˝] and an optimal D TI with a certain D TR .The effect of cross coupling on power transfer is further analyzed based on power flow analysis.
For three-resonator WPT system, this investigation provides a deeper understanding of the WPT characteristics.The proposed optimizing method can effectively improve the PTE or PDL, especially for a fixed transmission distance case and can obtain a constant PDL transfer at maximum PTE for different transmission distance.

Figure 3 .
Figure 3. Schematic of two coils for mutual inductance calculation.

Figure 3 .
Figure 3. Schematic of two coils for mutual inductance calculation.

Figure 6 .
Figure 6.(a) Comparison between maximizing PTE goal by optimizing coupling coefficient vector, k, and three other power transfer scenarios; (b) Comparison between maximal coupling coefficient vector, [kTI,max, kIR,max], of collinear axes arrangement scenario and optimal coupling coefficient vector, [kTI, kIR], for maximizing PTE goal.

Figure 6 .
Figure 6.(a) Comparison between maximizing PTE goal by optimizing coupling coefficient vector, k, and three other power transfer scenarios; (b) Comparison between maximal coupling coefficient vector, [k TI,max , k IR,max ], of collinear axes arrangement scenario and optimal coupling coefficient vector, [k TI , k IR ], for maximizing PTE goal.

Figure 7 .
Figure 7. (a) Comparison between maximizing PDL goal by optimizing coupling coefficient vector, k, and three other power transfer scenarios.(b) Comparison between maximal coupling coefficient vector, [kTI,max, kIR,max], of collinear axes arrangement scenario and optimal coupling coefficient vector, [kTI, kIR], for maximizing PDL goal.

kDFigure 8 .Figure 7 .
Figure 8. Predetermined-goal of obtaining a constant PDL of 0.5 W at maximum PTE by optimizing

Figure 7 .
Figure 7. (a) Comparison between maximizing PDL goal by optimizing coupling coefficient vector, k, and three other power transfer scenarios.(b) Comparison between maximal coupling coefficient vector, [kTI,max, kIR,max], of collinear axes arrangement scenario and optimal coupling coefficient vector, [kTI, kIR], for maximizing PDL goal.

Figure 8 .
Figure 8. Predetermined-goal of obtaining a constant PDL of 0.5 W at maximum PTE by optimizing coupling coefficient vector, k, in three-resonator WPT: (a) PTEs of collinear axes arrangement scenario, maximizing PTE goal by optimizing k and obtaining a constant PDL at maximum PTE by optimizing k; (b) PDLs of collinear axes arrangement scenario, maximizing PDL goal by optimizing k and obtaining a constant PDL at maximum PTE by optimizing k; and (c) values of k for collinear axes arrangement scenario and obtaining a constant PDL at maximum PTE.

Figure 8 .
Figure 8. Predetermined-goal of obtaining a constant PDL of 0.5 W at maximum PTE by optimizing coupling coefficient vector, k, in three-resonator WPT: (a) PTEs of collinear axes arrangement scenario, maximizing PTE goal by optimizing k and obtaining a constant PDL at maximum PTE by optimizing k; (b) PDLs of collinear axes arrangement scenario, maximizing PDL goal by optimizing k and obtaining a constant PDL at maximum PTE by optimizing k; and (c) values of k for collinear axes arrangement scenario and obtaining a constant PDL at maximum PTE.

Figure 9 .
Figure 9. Experimental setup: (a) Measurement platform; and (b) Block diagram of a two-port network is terminated by arbitrary RS and RL for calculating the PTE and PDL of the measuring system.

Figure 9 .
Figure 9. Experimental setup: (a) Measurement platform; and (b) Block diagram of a two-port network is terminated by arbitrary R S and R L for calculating the PTE and PDL of the measuring system.

Figure 10 .
Figure 10.Measured and theoretical results of three-resonator WPT system under collinear axes arrangement scenario at different inserted distance, DTI, when transfer distance DTR = 0.7 m: (a) PTE vs. frequency; and (b) RX I vs. frequency.

Figure 10 .Figure 11 .
Figure 10.Measured and theoretical results of three-resonator WPT system under collinear axes arrangement scenario at different inserted distance, D TI , when transfer distance D TR = 0.7 m: (a) PTE vs. frequency; and (b) |I RX | vs. frequency.

Figure 12 .Figure 11 .
Figure 12.(a) Measured and theoretical PTE of collinear axes arrangement scenario and maximizing PTE by optimizing θ at DTI = 0.45 m when transfer distance DTR = 0.7 m; (b) Measured and theoretical PDL of collinear axes arrangement scenario and maximizing PDL by optimizing θ at DTI = 0.3 m when transfer distance DTR = 0.7 m.

Figure 11 .
Figure 11.Angle, θ, vs. coupling coefficient, k, at different distances, D, without lateral misalignments of the two coils in Figure 3.

Figure 12 .Figure 12 .
Figure 12.(a) Measured and theoretical PTE of collinear axes arrangement scenario and maximizing PTE by optimizing θ at DTI = 0.45 m when transfer distance DTR = 0.7 m; (b) Measured and theoretical PDL of collinear axes arrangement scenario and maximizing PDL by optimizing θ at DTI = 0.3 m when transfer distance DTR = 0.7 m.

Figure 13 .
Figure 13.Measured and theoretical results of the predetermined-goal of obtaining a constant PDL at f0 = 1.96MHz when transfer distance DTR = 0.7 m.
Power loss in resonator of

Figure 13 .
Figure 13.Measured and theoretical results of the predetermined-goal of obtaining a constant PDL at f 0 = 1.96MHz when transfer distance D TR = 0.7 m.

Figure 14 .
Figure 14.Energy transfer of obtaining predetermined constant PDL at maximum PTE at 1.96 MHz when PDLC = 1 W and 0.5 W: (a) Ratio of optimal DTI to DTR vs. DTR; (b) Comparison between the calculated and measured values of PTE; and (c) Comparison between the calculated and measured values of PDL.

Figure 15 .
Figure 15.Power delivered to RX and power loss in it at 1.96 MHz under considering cross coupling model when PDLC = 1 W and 0.5 W.

Figure 14 .Figure 14 .Figure 15 .
Figure 14.Energy transfer of obtaining predetermined constant PDL at maximum PTE at 1.96 MHz when PDL C = 1 W and 0.5 W: (a) Ratio of optimal D TI to D TR vs. D TR ; (b) Comparison between the calculated and measured values of PTE; and (c) Comparison between the calculated and measured values of PDL.

Figure 15 .
Figure 15.Power delivered to RX and power loss in it at 1.96 MHz under considering cross coupling model when PDL C = 1 W and 0.5 W.
TI , k IR ), corresponding to local maximums of PTE, occur at both high and low frequency ranges.The locally minimal F(ω, k TI , k IR ) values for each case are symbolized as F H´Q S (F L´Q S ), F H´Q S0 (F L´Q S0 ), F H´Q S`( F L´Q S`) and F H´Q S´( F L´Q S´) at high (low) frequency range, respectively.The normalized frequencies of the local minimums of F(ω, k TI , k IR ) for Q S , Q S0 , Q S+ and Q S´a re NF H´Q S (NF L´Q S ), NF H´Q S0 (NF L´Q S0 ), NF H´Q S`( NF L´Q S`) and NF H´Q S´( NF L´Q S´) , respectively.It can be seen that F H´Q S0 (F L´Q S0 ) and NF H´Q S0 (NF L´Q S0 ) are close to F H´Q S (F L´Q S ) and NF H´Q S (NF L´Q S ) at high (low) frequency range respectively, therefore, we substitute Q S0 for Q S in Equation (3) to study the PTE characteristics.Frequency characteristics of PTE can be analyzed via ω v only.b e and c e are rewritten as b e0 and c e0 when Q S is substituted by Q S0 .The denominator in Equation (3) is extracted to form a function 2 21 = 1, α 2 31 = 1, β 2 21 = 0.5 and β 2 31 = 5.For each quality factor, two local minimums of F(ω, k |I RX |, can present the characteristics of PDL because of PDL = R L |I RX | 2 .|I RX | is compiled according to Equation (2c) as follow:

Table 2 .
Properties and graphs of g(x) = x 3 + b p0 x 2 + c p0 x + d p0 for analyzing |I RX | (or power delivered to the load (PDL)).
for frequency splitting for PTE and RX I of the sample case.Combining x0 in this figure with frequency splitting discriminance in Table 1, it can be found that PTE is in the non-frequency splitting region, critical frequency splitting position and frequency splitting region when DTI <, = and > 0.408 m, respectively.Ignoring the cross coupling between TX and RX, the values of PTE versus distance and frequency are shown in Figure5a.The frequencies for local maximum (FLM) of PTE under different DTI, shown as black lines, are achieved by Equation (4).The frequencies for maximum (FM) of PTE are obtained by searching the maximum in global frequency domain at each DTI, shown as white dotted lines.It can be seen that the critical frequency splitting distance is at 0.408 m, which validates the analysis result of judging the position for frequency splitting.The FM matches well with the FLM in non-frequency splitting region where DTI < 0.408 m, and is very close to the FLM at high splitting frequency in frequency splitting region where DTI > 0.408 m.In the greater degree of frequency splitting region, i.e., the region of splitting frequencies keeping farther away from each other, the FM slightly deviates from FLM, seeing the distribution of white dotted line and black line at larger DTI in Figure5a.The reason is that the splitting frequencies are calculated under the constant QS0.The values of PTE versus maximal coupling coefficient vector, [kTI,max, kIR,max], of collinear axes arrangement scenario and optimal coupling coefficient vector, [kTI, kIR], for maximizing PDL goal.