Next Article in Journal
Comparison of Torrefaction and Hydrothermal Treatment as Pretreatment Technologies for Rice Husks
Previous Article in Journal
Numerical Study on the Thermal Performance of a Single U-Tube Borehole Heat Exchanger Using Nano-Enhanced Phase Change Materials
Previous Article in Special Issue
Scaling-Factor and Design Guidelines for Shielded-Capacitive Power Transfer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Terminations for a Single-Input Multiple-Output Resonant Inductive WPT Link

by
Giuseppina Monti
1,*,†,
Mauro Mongiardo
2,†,
Ben Minnaert
3,
Alessandra Costanzo
4 and
Luciano Tarricone
1
1
Department of Engineering for Innovation, University of Salento, 73100 Lecce, Italy
2
Department of Engineering, University of Perugia, 06123 Perugia, Italy
3
Department of Industrial Science and Technology, Odisee University College of Applied Sciences, 9000 Ghent, Belgium
4
Department of Electrical, Electronic and Information Engineering Guglielmo Marconi, University of Bologna, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Energies 2020, 13(19), 5157; https://doi.org/10.3390/en13195157
Submission received: 28 May 2020 / Revised: 19 September 2020 / Accepted: 22 September 2020 / Published: 3 October 2020
(This article belongs to the Special Issue Modelling of Wireless Power Transfer)

Abstract

:
This paper analyzes a resonant inductive wireless power transfer link using a single transmitter and multiple receivers. The link is described as an ( N + 1 ) –port network and the problem of efficiency maximization is formulated as a generalized eigenvalue problem. It is shown that the desired solution can be derived through simple algebraic operations on the impedance matrix of the link. The analytical expressions of the loads and the generator impedances that maximize the efficiency are derived and discussed. It is demonstrated that the maximum realizable efficiency of the link does not depend on the coupling among the receivers that can be always compensated. Circuital simulation results validating the presented theory are reported and discussed.

1. Introduction

In recent years, several applications have been proposed for resonant inductive Wireless Power Transfer (WPT) [1,2,3,4]. In fact, resonant inductive WPT is an effective solution for wirelessly energizing electronic devices and several optimal design strategies have been investigated in the literature.
Usually, the goal is to recharge a single device and the focus is on maximizing either the power delivered to the load or the power transfer efficiency. In this regard, the most widely adopted scheme is that using a single transmitter, thus corresponding to a Single-Input Single Output (SISO) configuration. In a SISO configuration the link consists of just two magnetically coupled resonators: a transmitting resonator connected to the source and a receiving resonator connected to the load (i.e., the device to be recharged). SISO configurations have been widely investigated in the literature and it has been demonstrated that the link has to be terminated on its conjugate image impedances for maximizing both the power on the loads and the efficiency [5,6,7].
More recently, schemes using multiple transmitters and/or multiple receivers have been also investigated. The use of Multiple Input Single Output (MISO) schemes could be adopted to obtain an almost constant performance on a given area/volume this being useful if the position of the receiver is affected by small uncertainties (as in the case of embedded devices). In this regard, some interesting results are reported in [8] where it is demonstrated that a two-dimensional region of nearly constant power transfer efficiency can be obtained by using four transmitters. In [9] the use of a linear array of transmitters, activated two at time, is suggested for providing a constant output voltage to a load moving along a linear path. The problem of maximizing the efficiency and the power on the load in MISO schemes has been also analyzed and some interesting results have been reported in [10,11]. In particular, in [10] the solution for maximizing the efficiency has been formulated as a convex optimization problem. In [11] the optimal loads for both the maximum power and the maximum efficiency solutions have been presented for the case of a link using either two–transmitter and a single load or a single transmitter and two–load. In [12], a more abstract approach was used to maximize the efficiency by modeling the MISO-WPT system as a linear circuit whose input-output relationship is expressed in terms of a small number of unknown parameters that can be thought of as transimpedances and gains.
As per schemes using a Single Transmitter and Multiple Receivers (SIMO), they are adopted to recharge multiple devices with a single transmitter [13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]. In [20] the use of a multiple-output scheme is suggested for the recharge of electric vehicles. The problem of maximizing the power delivered to the loads has been solved in [21], where the expressions of the optimal loads have been derived by using the maximum power transfer theorem for an N–port.
As per the problem of efficiency maximization, in [23] the use of suitable matching networks is suggested. In [24], the specific case of a link using two receivers is analyzed and it is demonstrated that for some specific configurations of the receivers it is convenient to use a non-synchronous scheme with receivers resonating at a frequency different from that of the transmitter. In [22] a SIMO system with constant output voltage and operating at 6.78 M Hz is presented. The efficiency of the proposed WPT link is optimized by tuning the input voltage at the transmitter side.
In [25], the loads for maximizing the efficiency have been derived from the expression calculated for the case of a link using one receiver and that using two receivers. However, the analysis is performed assuming that the coupling among the receivers can be neglected, this representing a limitation for real applications. The presence of possible couplings among the receivers has been analyzed in [26,27]. It is demonstrated that for given loads a coupling among the receivers can be compensated by using suitable compensating reactances; however, in these papers it is assumed that the loads are given (i.e., they are not optimized).
Finally, for the problem of efficiency maximization, elegant and comprehensive analysis of all possible configurations (i.e., the SIMO, MISO and MIMO configurations) have been presented in [29,30]. A very elegant and general approach is presented in [29]; where, starting from the impedance or scattering matrix of a multiport the efficiency of a generic MIMO-WPT system is expressed by the Rayleigh quotient. However, the method is not applied on an inductive WPT system and the optimal loads are only expressed as function of the port currents and impedance matrix elements. In [30], the optimal loads are derived from the first-order necessary condition consisting of imposing the zeroing of the first-order partial derivatives of the efficiency with respect to the input and output currents. The optimal solution derived in this way is validated by checking the second order derivatives. The developed analysis is general and overcomes some limitations present in the previous literature. For instance, for the SIMO case a generic number of possibly coupled receivers is considered. Similarly, for the MISO case, the formulas are presented for a generic number of possibly coupled transmitters. However, the analysis developed in [30] is based on the assumption that all the couplings among the transmitters and the receivers are purely inductive; this assumption limits the applicability of the approach to practical applications where the conductivity of the propagation channel is negligibly small.
In this paper, referring to the SIMO configuration, similarly to [29], the problem of finding the optimal loads maximizing the efficiency is formulated as a generalized eigenvalue problem. The presented theory is valid for any strictly passive and reciprocal network in SIMO configuration and is applied in detail for the first time in this paper to the case of a resonant inductive WPT link. The application of the presented theory just requires the knowledge of the impedance matrix of the SIMO network that can be the result of measurements, simulations or theoretical derivation. The network must not satisfy any particular hypothesis except that of being passive and reciprocal; consequently, the proposed approach is also applicable in the case of non-purely inductive couplings (including the case of a propagation channel with non-negligible values of the conductivity).
The general theory is first presented for a generic ( N + 1 ) –port network in SIMO configuration and then applied to the specific case of a resonant inductive WPT link; the analytical expressions of the complex loads maximizing the efficiency are derived and discussed. Additionally, the importance of suitably selecting the generator impedance for maximizing the total output power corresponding to the maximum efficiency solution is discussed. The correctness of the derived expressions is validated by the results reported in [30] and by numerical data presented in this paper.
The paper is organized as follows. In Section 2 the problem of efficiency maximization is solved for a generic SIMO ( N + 1 ) –port network. In Section 3 the derived equations are specialized for the case of an inductive WPT link, the optimal expressions of the loads and the generator impedances are reported. In Section 4 theoretical formulas are validated through circuital and full-wave simulations. Finally, some conclusions are drawn in Section 5.

2. Derivation of the Solution for the General Case

The problem analyzed in this paper is a WPT link using a Single-Input Multiple-Output (SIMO) configuration: a single transmitter is wirelessly connected to N receivers. In this section, the general case is analyzed, no specific assumption is made on the coupling mechanism among the transmitter and the receivers, it is only assumed that the network is passive and reciprocal.
By using a network formalism, the link is modeled as an ( N + 1 ) –port network N , see Figure 1, described by its impedance matrix Z . The input port is connected to a sinusoidal source V G with internal impedance Z G and the output ports are connected to an N-port load N L with impedance matrix Z L . Generally, in practical cases, N L consists of a set of N uncoupled load impedances and, consequently, Z L is a diagonal matrix
Z L = diag ( Z L , n ) ,
with n = 1 , , N .
In real applications, the generator could be a complex network, comprising a DC-AC converter and other circuitry. Accordingly, in general, Z G is the input impedance of the network adopted for generating the power to be provided at the input port of the network N . The same consideration applies for each load. In fact, in real applications each load can be a more or less complicated network which in most cases includes a rectifier for converting the AC power at the output port of the network into a DC signal. Accordingly, the generic impedance Z L i is the input impedance of the network connected to the output port i of the link.
The vectors of voltage and current phasors at the network ports, V and I , and the matrix Z can be partitioned as
V i V o = Z ii Z io Z oi Z oo I i I o
where V i and I i represent voltage and current at the input port, while V o and I o are the N-vectors of voltages and currents at the output ports.
By replacing the load equation
V o = Z L I o
in (2), and by eliminating I o , the impedance seen at the input port of N can be derived as
Z in = V i I i = Z ii Z io Z oo + Z L 1 Z oi .
In a similar way, by combining (2) with the source equation
V i = V G Z G I i
and eliminating I i , the relation between voltages and currents at the output ports can be cast in the form
V o = V th + Z out I o
where
V th = Z oi V G Z ii + Z G
is a set of N Thévenin equivalent voltage sources and
Z out = Z oo Z oi Z io Z ii + Z G
is the equivalent impedance matrix of the network N with the input port closed on the impedance Z G . The network N can be thus represented by the equivalent circuit of Figure 2.
The maximum power transfer between the source and the input port of N can be achieved when the conjugate match condition
Z G = Z in * ,
where * denotes conjugation, is satisfied. In this case, the power delivered to N is equal to the generator available power
P AG = | V G | 2 8 Re Z G .
As far as the output side is concerned, it can be proved [31] that the power delivered by N is maximized when the output currents I o assume the values I oM given by
I oM = Z out + Z out 1 V th
where denotes conjugate transpose, and consequently the available power at the output ports of N is
P a = 1 4 V th Z out + Z out 1 V th .
It can be noted that for N > 1 , the optimal load is not univocally defined. In fact, the optimal currents can be obtained by any impedance matrix Z Lm such that
Z LM I oM = Z out I oM = V oM
where V oM are the voltages at output ports for I o = I oM . Equation (13) also shows that it is possible to realize Z LM as a set of N independent passive impedances provided that the possible zero elements of I oM corresponds to zero elements of V oM , and that the phase difference between any two corresponding elements of I oM and V oM is ≥ 90 in absolute value.
According to the previous discussion, also the problem of determining the impedances Z G and Z L which provide the simultaneous maximum power transfer at the input and output ports has not a unique solution.
To simplify the calculation of the optimal terminations, it is convenient to determine the corresponding optimal currents, which, on the contrary, are univocally defined.
Making use of (2), the total power delivered to the loads P o , i.e., the sum of the powers delivered to each load P o i
P o = n = 1 N P o i ,
can be expressed as a function of the port currents as
P o = 1 4 V o I o + I o V o = = 1 4 I o Z oo + Z oo I o + I o Z oi I i + I i * Z oi I o
and, similarly, the input power can be expressed as
P i = 1 4 V i * I i + I i * V i = = 1 4 I i * Z ii + Z ii * I i + I i * Z io I o + I o Z io I i .
The previous equations can be cast in the form
P o = 1 4 I A I P i = 1 4 I B I
where the matrices A and B are defined as
A = 0 Z oi Z oi Z oo + Z oo
B = Z ii + Z ii * Z io Z io 0 .
The power gain of N , defined as the ratio between the output and the input power, can thus be expressed as
G p = P o P i = I A I I B I .
In the context of WPT the quantity expressed in (20) is usually referred to as the efficiency of the link, in this paper it will be referred to as G p in analogy with the terminology adopted in the context of two-port networks.
The power gain is maximized when the maximum power transfer is realized at the output port. Since G P is a generalized Rayleigh quotient, the maximum of G P can be determined by solving a generalized eigenvalue problem.
As a matter of fact, using the quotient rule and taking into account the fact that A and B are Hermitian matrices, the differential of G P can be calculated as
δ G p = 2 ( δ I A I ) ( I B I ) ( I A I ) ( δ I B I ) I B I 2 .
Hence, requiring δ G p = 0 yields
A I I A I I B I B I = 0 ,
which can be rewritten as
A x = λ B x
and can be recognized as a generalized eigenvalue problem with λ = G p being the eigenvalue and x = I the corresponding eigenvector.
Since, by hypothesis, N is passive and the power supply is provided only at the input port, the maximum power gain, G M , and the corresponding currents (up to an arbitrary factor) can be determined by solving (23) with the constrains
λ 1 P i 0 P o 0 .
After determining the optimal currents
I M = I iM I oM
by (23), the corresponding voltages
V M = V iM V oM
can be obtained by (2). Hence the source impedance providing maximum power transfer at the input port can be calculated by letting
Z GM = V iM * I iM * ,
while the maximum power transfer at the output port is obtained with any load N L whose impedance matrix satisfies (13). In particular, if the previously enunciated conditions are satisfied, N L can be realized as a set of uncoupled loads with impedances
Z LM , n = V oM , n I oM , n
with n = 1 , , N .
It is worth observing that the theory presented in this section is completely general, it can be applied to any passive SIMO network; moreover, for its application it is sufficient to know the impedance matrix of the network. It is possible to derive the maximum achievable power gain and the optimal loads starting from the impedance matrix, which can be the results of measurements, theoretical calculation or a numerical analysis. Figure 3 summarizes how to apply the proposed approach for the determination of the load impedances maximizing the efficiency of a SIMO Resonant Inductive WPT Link.

3. The Case of an Inductive Resonant Coupling

In this section, the specific case of a WPT link consisting of ( N + 1 ) magnetically coupled resonators is considered (Figure 4). More specifically, the link consists of ( N + 1 ) magnetically coupled inductors, L i , each one loaded by a suitable compensating capacitor, C i , realizing the resonance condition at the operating angular frequency (i.e., ω 0 = 1 / L i C i ). The inductor losses are modeled by series resistors R i related to the quality factors of the coupled resonators:
Q n = ω L n R n .
The coupling between the inductors L m and L n is described by the coupling factor k m n related to the mutual inductance M m n
k m n = M m n L m L n .
Accordingly, the network is described by the following impedance matrix:
Z = R 0 j ω M 01 j ω M 02 j ω M 0 N j ω M 01 R 1 j ω M 12 j ω M 1 N j ω M 02 j ω M 12 R 2 j ω M 2 N j ω M 0 N j ω M 1 N j ω M 2 N R N .
By introducing the normalization matrix d :
d = diag 1 ω L n , n = 0 , , N ,
it is possible to obtain the following normalized expression for the impedance matrix of the network:
z = d Z d = 1 Q 0 j k 01 j k 02 j k 0 N j k 01 1 Q 1 j k 12 j k 1 N j k 02 j k 12 1 Q 2 j k 2 N j k 0 N j k 1 N j k 2 N 1 Q N .
Referring to Section 2 and to the Appendix A, for the specific analyzed case it is possible to derive:
z ˜ ii = 2 Q 0 z io = j k 01 j k 0 N z oi = z io T z ˜ oo = diag 2 Q n
and
c 0 = n = 1 N k 0 n 2 Q n c 1 = c 0 2 Q 0 c 2 = c 0 .
Accordingly, by introducing the parameter α :
α = 1 + n = 1 N k 0 n 2 Q 0 Q n ,
for the analyzed case, the solving equation is (see the Appendix A):
( α 2 1 ) λ 2 2 ( α 2 + 1 ) λ + ( α 2 1 ) = 0 .
Equation (37) has two eigenvalues:
G M = α 1 α + 1 , G M 1 = α + 1 α 1 .
It is evident that G M 1 > 1 ; as a consequence, only G M satisfies the first constrain expressed in (24).
By choosing to normalize the input current to 1
i iM = 1 ,
the following normalized eigenvectors can be obtained
i oM , n = j k 0 n Q n α + 1 .
The corresponding normalized voltages are:
v iM = α Q 0 ,
v oM , n = 1 α + 1 m = 1 m n N k 0 m k n m Q m + j k 0 n α .
Hence the optimal normalized source impedance is given by:
z GM = v iM * i iM * = α Q 0 ,
however, the optimal N-port load network N L can be realized as a set of uncoupled loads with normalized impedances:
z LM , n = v oM , n i oM , n = r LM , n j x LM , n
r LM , n = α Q n
x LM , n = 1 k 0 n Q n m = 1 m n N k 0 m k n m Q m
The corresponding unnormalized expressions are:
R LM , n = α R n
X LM , n = R n k 0 n m = 1 m n N k 0 m k n m Q m .

Discussion of the Results

According to the above reported formulas, the following considerations can be drawn.
  • From (38) it is evident that the maximum realizable efficiency of the link only depends on the quality factors of the resonators and on the couplings between the transmitting and the receiving resonators; however, it does not depend on the couplings among the receivers. This means that a possible coupling among the receivers can be always compensated.
  • In general, the optimal loads are complex quantities.
  • For identical resonators with the same quality factor the real part of the optimal loads is the same for all the loads.
  • The imaginary parts of the optimal loads are zero for uncoupled receiving resonators; this means that they play a role of compensation.
  • By comparing the expression of the reactive parts of the optimal loads with those of the optimal loads reported in [21] for the maximum power case, it can be easily verified that they are coincident. This means that the same compensating reactances are required for both the maximum efficiency case (i.e., for maximizing the power gain) and the maximum power case.
  • The proposed approach also provides the optimal value of the generator impedance, see (27); however, G p does not depend on the generator. The value provided for Z G in (27) is that maximizing the power entering the network, and then the power delivered to the loads, when the loads are those maximizing G p .
It is worth observing that all the achieved results are in a perfect agreement with those reported in [30]. With respect to previously proposed approaches, the theory presented in this paper has the advantage of being completely general, it is valid for any strictly passive and reciprocal network. Additionally, the application of the presented theory just needs the impedance matrix of the link that can be the result of measurements or simulations or theoretical evaluation. In fact, the optimal loads are obtained by solving the eigenvalue problem expressed in (23). To solve the eigenvalue problem one just needs the matrices A and B that can be directly computed from the Z matrix, see (18) and (19).

4. Validation of the Results

To validate the theoretical data, full-wave and circuital simulations have been performed. The commercial tool CST Microwave Studio has been used for full-wave simulations, while the NI AWR Design Environment has been adopted for circuital simulations. Four different WPT links in SIMO configuration have been analyzed. The first three analyzed cases have identical resonators and a different number of receivers, as detailed in the following.
All the analyzed coils (transmitter and receivers) have the same dimensions; they are circular loops with a radius of 5 m m designed by using a copper wire with a radius of 0.3 m m . An operating frequency f 0 of 500 M Hz has been assumed. First, the single loop has been analyzed so to calculate the equivalent inductance. From full-wave simulations at f 0 each loop corresponds to an inductance of about 20.9 n H . Accordingly, a series capacitor of 4.84 p F has been added to each loop so to make them resonating at f 0 . The relative positions of the transmitting and the receiving coils assumed for the three analyzed cases are illustrated in Figure 5.
To calculate the impedance matrices, each link has been analyzed through full-wave simulations as a multiport network. The following impedance matrices have been obtained:
Z Case 1 = 0.118 1.41 j 1.41 j 1.41 j 0.118 2.62 j 1.41 j 2.62 j 0.118 ,
Z Case 2 = 0.118 2.66 j 2.66 j 1.17 j 2.66 j 0.118 1.84 j 1.86 j 2.66 j 1.84 j 0.118 1.86 j 1.17 j 1.86 j 1.86 j 0.118 ,
Z Case 3 = 0.118 0.81 j 0.81 j 1.37 j 1.37 j 0.81 j 0.118 2.60 j 0.56 j 2.58 j 0.81 j 2.60 j 0.118 2.58 j 0.56 j 1.37 j 0.56 j 2.58 j 0.118 2.60 j 1.37 j 2.58 j 0.56 j 2.60 j 0.118 .
By comparing the general expression of the impedance matrix of a resonant inductive WPT link given in (31) with the numerical values calculated through circuital simulations, the values reported in Table 1, Table 2 and Table 3 have been derived for the coupling coefficients. By using (47) and (48) it is possible to calculate the expressions of the optimal loads. The values calculated for the three analyzed examples are summarized in Table 1, Table 2 and Table 3. With all the resonators the same quality factor, the resistive part of the loads is the same for all the receivers.
As per the imaginary parts, for the analyzed cases all the calculated values of X LM , n are negative, thus corresponding to load impedances with an inductor L M , n in series configurations with the resistive part R L M , n .
The analytical values of the optimal loads have been validated through circuital simulations.
Two different sets of circuital simulations have been performed. A first set of simulations has been performed by modeling the links with lumped elements equivalent circuits with the parameters summarized in Table 1, Table 2 and Table 3; Figure 6 illustrates the circuit analyzed for case 1. The resistors R n that appear in Figure 6 are related to the quality factors of the resonators, Q n , through (29). A second set of simulations has been performed by modeling the analyzed links as ( N + 1 ) –port black–box networks described by the impedance matrices provided by full–wave simulations, being N the number of receivers (i.e., N = 2 for case 1, N = 3 for case 2, N = 4 for case 3). In more detail, referring to case 1, simulations have been performed by replacing the network in the dashed square of Figure 6 with a three–port black-box component described by the impedance matrix of the link calculated through circuital simulations.
First, the optimal values provided by the theory for the resistive parts of the loads have been validated. Simulations have been performed by varying the resistive part of the loads, the values calculated for G P are given in Figure 7a–c.
The results obtained for the case of purely resistive loads (i.e., Z L n = R L ) are compared with those obtained for the case of loads with the compensating inductances given in Table 3, i.e., for Z L n = R L + j ω L L M , n . In the figures, the triangles have been used for the results obtained by modeling the link with the impedance matrix provided by full-wave simulations; however, the solid lines have been used for the results obtained by modeling the links with the lumped elements equivalent circuit. It can be seen that the results obtained for the two representations of the links are coincident.
As per the optimal values of R L , in each figure the value of R L for which G P is maximized is highlighted by a dashed vertical line. It can be seen that the values calculated through circuital simulations confirm the theoretical values. Finally, with regard to the compensating inductances L L M , n , it seems that they play a more or less important role in maximizing G P depending on the analyzed case. For instance, according to the achieved results the compensating reactances play a marginal role in maximizing G P for case 2 while they seem to be more relevant for case 3. However, for all the three analyzed cases it is confirmed that their presence allows obtaining the maximum value of G P provided by the theory.
The behavior of the power delivered to the loads as function of the generator impedance has been also investigated. Simulations have been performed by using as source a voltage generator with a series impedance R G . The total output power has been calculated by varying R G when the loads assume the optimal values provided by the theory. The results are given in Figure 8a–c, data obtained for Z L n = R L M , n and Z L n = R L M , n + j ω L L M , n are compared. The dashed vertical lines highlight the values of R G maximizing P o for the case Z L n = R L M , n + j ω L L M , n . It can be verified that these values are in a perfect agreement with the optimal values provided by the theory. As it can be seen, the use of the generator impedance provided by (43) allows maximizing the power delivered to the loads when they are set to maximize G P . Additionally, it can be seen that the compensation reactances are crucial to maximize the power transferred to the loads.
It is worth observing that the output power illustrated in Figure 8a–c is the total output power delivered (i.e., the sum of the power delivered to the loads) when the network operates at maximum G P . From the figures it is evident that if R G is not optimized, although the network operates with efficiency values close to one, only a small portion of the power available from the generator is delivered to the load.
Finally, the case of a link consisting of three coils with different dimensions has been analyzed (case 4). The geometry analyzed through full-wave simulations is illustrated in Figure 9. All the coils have been designed by using a copper wire with a radius of 0.3 m m . The radius of the transmitting coil is 10 m m , those of the first and second receivers are 7.5 m m and 5 m m , respectively. Also, in this case an operating frequency f 0 of 500 M Hz has been assumed.
The impedance matrix as calculated from full-wave simulations is:
Z Case 4 = 0.229 10.08 j 8.66 j 10.08 j 0.145 2.82 j 8.66 j 2.82 j 0.118 .
The parameters derived for the equivalent circuit are summarized in Table 4.
The simulated results obtained for G P are illustrated in Figure 10. Circuital simulations have been performed by modeling the link as a three–port black–box component described by the impedance matrix calculated through full–wave simulations. In this case, the two receivers have slightly different values of the quality factors; accordingly, the theory predicts slightly different values for R L M , 1 and R L M , 2 . To verify the expected optimal values, simulations have been performed by terminating the receivers ports on the impedances Z L 1 = R L 1 + j ω L L M , 1 and Z L 2 = R L 2 + j ω L L M , 2 and by varying both R L , 1 and R L , 2 . From Figure 10 it can be seen that circuital simulations confirm the theory, a maximum of about 0.98 is obtained for G P when R L 1 = 11.07 Ω and R L 2 = 9.01 Ω ; however, from the figure it can also be seen that values of G P very close to its maximum (i.e., values greater than 0.97) are obtained for a wide range of values of R L 1 and R L 2 .
Finally, the behavior obtained for the total output power P o as function of the generator impedance is illustrated in Figure 11. In this case, simulations have been performed by terminating the ports of the receivers on the optimal load impedances, i.e., Z L 1 = R L M , 1 + j ω L L M , 1 and Z L 2 = R L M , 2 + j ω L L M , 2 . As for the previously analyzed cases, simulations confirm the importance of suitably selecting the generator impedance for maximizing P o when the loads are those maximizing G P .

5. Conclusions

In this paper, the case of a WPT link using a Single-Input Multiple-Output (SIMO) configuration is analyzed. The solution for maximizing the efficiency is derived from a generalized eigenvalue problem. The main advantage of the presented theory is its complete generality, no specific assumptions are made about the link. In fact, the proposed approach is valid for any strictly passive and reciprocal ( N + 1 ) –port network. Additionally, the desired solution can be derived directly from the impedance matrix with simple algebraic operations. Theoretical formulas for the optimal loads maximizing the efficiency are derived for the case of a resonant inductive link with a generic number N of possibly coupled receivers.
The results obtained this way for the optimal loads are coincident with those reported in the previous literature were the first-order condition on the partial derivatives of the efficiency has been exploited in order to find the desired solution. As a further validation of the derived formulas, several numerical examples have been analyzed through full–wave and circuital simulations. In particular, four different links have been considered. The first three analyzed links use identical resonators for the transmitter and the receivers and differ for the number of receivers. The last analyzed case is a link using two receivers; in this case, the three resonators have been designed so to have different values of the equivalent inductance. More specifically, the transmitter has been designed so to have a larger inductance with respect to the receivers, so to obtain higher values of the couplings with respect to the previously analyzed cases. For all the investigated cases, the impedance matrix has been calculated through full-wave simulations and the presented theory applied so to determine the maximum realizable efficiency and the optimal terminating impedances. The correctness of the analytical data has been verified through simulations performed for evaluating the efficiency of the links. According to the theoretical data, simulations confirm that the maximum realizable efficiency of a resonant inductive link in SIMO configuration does not depend on the coupling among the receiving resonators. In fact, it is demonstrated that a possible coupling among the receivers can be always compensated by using suitable complex loads.
The possibility of maximizing the power delivered to the loads when they are set to maximize the efficiency has been also discussed. In fact, the presented theory provides both:
  • the expressions of the loads that maximize the efficiency,
  • the expression of the generator impedance that allows maximizing the power entering the network when the loads are those maximizing the efficiency.
Maximizing the power entering the network for a given efficiency corresponds maximizing the power delivered to the loads. The reported results highlight the importance of also optimizing the generator impedance to avoid that. Despite the high efficiency values, only a small portion of the power available from the generator is transferred to the loads.
As future developments of the presented research, experimental tests will be performed to verify the application of the proposed theory to a real application. Furthermore, in a future work the analysis presented in this paper will be extended to a generic MIMO (Multiple Input Multiple-Output) system.

Author Contributions

Conceptualization, methodology, validation, G.M., M.M.; writing—original draft preparation, G.M., M.M.; writing—review and editing, B.M., A.C. and L.T. 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 remember the colleague Franco Mastri who suddenly passed away on 3 April 2020. He was a great colleague and a profound scientist. Fundamental discussions and studies on the theoretical modeling of near-field WPT systems were of great inspiration also for the results presented in this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. How to Solve the Generalized Eigenvalue Problem

The following generalized eigenvalue problem is considered:
A x = λ B x .
In this Appendix the procedure for calculating the eigenvalues λ and the corresponding eigenvectors x = I will be illustrated. By introducing the variable η :
η = 1 λ
it is possible to write:
B η A x = 0 .
By using (18) and (19), it is possible to obtain:
B η A = Z ˜ ii Z io + η Z oi Z io + η Z oi η Z ˜ oo ,
where the following definitions have been introduced:
Z ˜ ii = Z ii + Z ii * Z ˜ oo = Z oo + Z oo .
The eigenvalues can be obtained from:
det ( B η A ) = 0 .
Considering that for a matrix M partitioned in 4 submatrix M 11 , M 12 , M 21 , M 22 :
M = M 11 M 12 M 21 M 22 ,
the determinant is given by:
det ( M ) = det ( M 11 M 12 M 22 1 M 21 ) det ( M 22 ) ,
it is possible to write
det ( B η A ) = η N 1 ( c 0 η 2 + c 1 η + c 2 ) det ( Z ˜ oo ) .
where the coefficients c 0 , c 1 and c 2 are given by:
c 0 = Z oi Z ˜ oo 1 Z oi , c 1 = Z ˜ ii + Z io Z ˜ oo 1 Z oi + Z oi Z ˜ oo 1 Z io , c 2 = Z io Z ˜ oo 1 Z oi .
Accordingly, the non-trivial solutions are provided by the equation:
c 2 λ 2 + c 1 λ + c 0 = 0 .

References

  1. Karalis, A.; Joannopoulos, J.D.; Soljačić, M. Efficient wireless non-radiative mid-range energy transfer. Ann. Phys. 2008, 323, 34–48. [Google Scholar] [CrossRef] [Green Version]
  2. Monti, G.; Paolis, M.V.D.; Corchia, L.; Tarricone, L. Wireless Resonant Energy link for Pulse Generators Implanted in the Chest. IET Microw. Antennas Propag. 2017, 11, 2201–2210. [Google Scholar] [CrossRef]
  3. Carvalho, N.B.; Georgiadis, A.; Costanzo, A.; Stevens, N.; Kracek, J.; Pessoa, L.; Rogier, H. Europe and the future for WPT. IEEE Microw. Mag. 2017, 18, 56–87. [Google Scholar] [CrossRef] [Green Version]
  4. Rim, C.T.; Mi, C. Wireless Power Transfer for Electric Vehicles and Mobile Devices; Wiley-IEEE Press: Hoboken, NJ, USA, 2017. [Google Scholar]
  5. Inagaki, N. Theory of Image Impedance Matching for Inductively Coupled Power Transfer Systems. IEEE Trans. Microw. Theory Tech. 2014, 62, 901–908. [Google Scholar] [CrossRef]
  6. Mastri, F.; Mongiardo, M.; Monti, G.; Dionigi, M.; Tarricone, L. Gain expressions for resonant inductive wireless power transfer links with one relay element. Wirel. Power Transf. 2018, 5, 27–41. [Google Scholar] [CrossRef]
  7. Monti, G.; Costanzo, A.; Mastri, F.; Mongiardo, M. Optimal design of a wireless power transfer link using parallel and series resonators. Wirel. Power Transf. 2016, 3, 105–116. [Google Scholar] [CrossRef]
  8. Yoon, I.; Ling, H. Investigation of near-field wireless power transfer under multiple transmitters. IEEE Antenna Wirel. Propag. Lett. 2011, 10, 662–665. [Google Scholar] [CrossRef]
  9. Pacini, A.; Costanzo, A.; Aldhaher, S.; Mitcheson, P.D. Load- and Position-Independent Moving MHz WPT System Based on GaN-Distributed Current Sources. IEEE Trans. Microw. Theory Tech. 2017, 65, 5367–5376. [Google Scholar] [CrossRef] [Green Version]
  10. Lang, H.D.; Ludwig, A.; Sarris, C.D. Convex Optimization of Wireless Power Transfer Systems with Multiple Transmitters. IEEE Trans. Antennas Propag. 2014, 62, 4623–4636. [Google Scholar] [CrossRef]
  11. Monti, G.; Che, W.; Wang, Q.; Costanzo, A.; Dionigi, M.; Mastri, F.; Mongiardo, M.; Perfetti, R.; Tarricone, L.; Chang, Y. Wireless Power Transfer with Three-Ports Networks: Optimal Analytical Solutions. IEEE Trans. Circuits Syst. 2017, 62, 494–503. [Google Scholar] [CrossRef]
  12. Arakawa, T.; Goguri, S.; Krogmeier, J.V.; Kruger, A.; Love, D.J.; Mudumbai, R.; Swabey, M.A. Optimizing wireless power transfer from multiple transmit coils. IEEE Access 2018, 6, 23828–23838. [Google Scholar] [CrossRef]
  13. Kim, G.; Boo, S.; Kim, S.; Lee, B. Control of Power Distribution for Multiple Receivers in SIMO Wireless Power Transfer System. J. Electromagn. Eng. Sci. 2018, 18, 221–230. [Google Scholar] [CrossRef]
  14. Cai, W.; Ma, D.; Tang, H.; Lai, X.; Liu, X.; Sun, L. Highly Efficient Target Power Control for Two-Receiver Wireless Power Transfer Systems. Energies 2018, 11, 2726. [Google Scholar] [CrossRef] [Green Version]
  15. Liu, M.; Fu, M.; Wang, Y.; Ma, C. Battery Cell Equalization via Megahertz Multiple-Receiver Wireless Power Transfer. IEEE Trans. Power Electron. 2018, 33, 4135–4144. [Google Scholar] [CrossRef]
  16. Cai, W.; Lai, X.; Ma, D.; Tang, H.; Hashmi, K.; Xu, J. Management of Multiple-Transmitter Multiple-Receiver Wireless Power Transfer Systems Using Improved Current Distribution Control Strategy. Electronics 2019, 8, 1160. [Google Scholar] [CrossRef] [Green Version]
  17. Luo, C.; Qiu, D.; Lin, M.; Zhang, B. Circuit Model and Analysis of Multi-Load Wireless Power Transfer System Based on Parity-Time Symmetry. Energies 2020, 13, 3260. [Google Scholar] [CrossRef]
  18. Wagih, M.; Komolafe, A.; Zaghari, B. Dual-Receiver Wearable 6.78 MHz Resonant Inductive Wireless Power Transfer Glove Using Embroidered Textile Coils. IEEE Access 2020, 8, 24630–24642. [Google Scholar] [CrossRef]
  19. Ahn, D.; Kim, S.; Kim, S.; Moon, J.; Cho, I. Wireless Power Transfer Receiver with Adjustable Coil Output Voltage for Multiple Receivers Application. IEEE Trans. Ind. Electron. 2019, 66, 4003–4012. [Google Scholar] [CrossRef]
  20. Vu, V.; Phan, V.; Dahidah, M.; Pickert, V. Multiple Output Inductive Charger for Electric Vehicles. IEEE Trans. Power Electron. 2019, 34, 7350–7368. [Google Scholar] [CrossRef] [Green Version]
  21. Monti, G.; Dionigi, M.; Mongiardo, M.; Perfetti, R. Optimal Design of Wireless Energy Transfer to Multiple Receivers: Power Maximization. IEEE Trans. Microw. Theory Tech. 2017, 65, 260–269. [Google Scholar] [CrossRef]
  22. Fu, M.; Yin, H.; Liu, M.; Wang, Y.; Ma, C. A 6.78 MHz Multiple-Receiver Wireless Power Transfer System with Constant Output Voltage and Optimum Efficiency. IEEE Trans. Power Electron. 2018, 33, 5330–5340. [Google Scholar] [CrossRef]
  23. Lin, H.; Li, L. Efficiency Analysis and Optimization for Multiple-Receiver Magnetic Coupling Resonant Wireless Power Transfer System. In Proceedings of the 2020 5th International Conference on Computer and Communication Systems (ICCCS), Shanghai, China, 22–24 February 2020; pp. 742–747. [Google Scholar]
  24. Lee, K.; Chae, S.H. Comparative Analysis of Frequency-Selective Wireless Power Transfer for Multiple-Rx Systems. IEEE Trans. Power Electron. 2020, 35, 5122–5131. [Google Scholar] [CrossRef]
  25. 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, 3463–3477. [Google Scholar] [CrossRef]
  26. Fu, M.; Zhang, T.; Zhu, X.; Luk, P.C.; Ma, C. Compensation of Cross Coupling in Multiple-Receiver Wireless Power Transfer Systems. IEEE Trans. Ind. Inform. 2016, 12, 474–482. [Google Scholar] [CrossRef]
  27. Ishihara, M.; Fujiki, K.; Umetani, K.; Hiraki, E. Automatic Active Compensation Method of Cross-Coupling in Multiple-receiver Resonant Inductive Coupling Wireless Power Transfer Systems. In Proceedings of the 2019 IEEE Energy Conversion Congress and Exposition (ECCE), Baltimore, MD, USA, 29 September–3 October 2019; pp. 4584–4591. [Google Scholar]
  28. Duong, Q.T.; Okada, M. Maximum efficiency formulation for inductive power transfer with multiple receivers. IEICE Electron. Exp. 2016, 22, 20160915. [Google Scholar] [CrossRef] [Green Version]
  29. Yuan, Q.; Aoki, T. Practical applications of universal approach for calculating maximum transfer efficiency of MIMO-WPT system. Wirel. Power Transf. 2020, 7, 86–94. [Google Scholar] [CrossRef]
  30. Duong, Q.; Okada, M. Maximum Efficiency Formulation for Multiple-Input Multiple-Output Inductive Power Transfer Systems. IEEE Trans. Microw. Theory Tech. 2018, 66, 3463–3477. [Google Scholar] [CrossRef]
  31. Desoer, C. The maximum power transfer theorem for n-ports. IEEE Trans. Circuit Theory 1973, 20, 328–330. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of a SIMO WPT link.
Figure 1. Schematic representation of a SIMO WPT link.
Energies 13 05157 g001
Figure 2. Equivalent Thévenin representation of the circuit of Figure 1.
Figure 2. Equivalent Thévenin representation of the circuit of Figure 1.
Energies 13 05157 g002
Figure 3. Block diagram of the proposed approach to determine the optimal terminations for efficiency maximization of a SIMO Resonant Inductive WPT Link.
Figure 3. Block diagram of the proposed approach to determine the optimal terminations for efficiency maximization of a SIMO Resonant Inductive WPT Link.
Energies 13 05157 g003
Figure 4. Equivalent circuit of a WPT link with a single transmitter and N receivers, determined by its impedance matrix Z .
Figure 4. Equivalent circuit of a WPT link with a single transmitter and N receivers, determined by its impedance matrix Z .
Energies 13 05157 g004
Figure 5. WPT links analyzed through full-wave simulations. (a) Case 1: single transmitter and two-receiver link; (b) Case 2: single transmitter and three–receiver link; (c) Case 3: single transmitter and four-receiver link. In all cases the transmitter is the loop with the center in the point O 0 .
Figure 5. WPT links analyzed through full-wave simulations. (a) Case 1: single transmitter and two-receiver link; (b) Case 2: single transmitter and three–receiver link; (c) Case 3: single transmitter and four-receiver link. In all cases the transmitter is the loop with the center in the point O 0 .
Energies 13 05157 g005
Figure 6. Equivalent circuit analyzed for Case 1.
Figure 6. Equivalent circuit analyzed for Case 1.
Energies 13 05157 g006
Figure 7. Power gain calculated through circuital simulations by varying the resistive part of the loads R L n . (a) Case 1, the link has two receivers with R L 1 = R L 2 = R L ; (b) Case 2, the link has three receivers with R L 1 = R L 2 = R L 3 = R L ; (c) case 3, the link has four receivers with R L 1 = R L 2 = R L 3 = R L 4 = R L . The figure compares full-wave and circuital simulation results obtained for the case of purely resistive loads and for the case of loads with the compensating inductances given in Table 1, Table 2 and Table 3.
Figure 7. Power gain calculated through circuital simulations by varying the resistive part of the loads R L n . (a) Case 1, the link has two receivers with R L 1 = R L 2 = R L ; (b) Case 2, the link has three receivers with R L 1 = R L 2 = R L 3 = R L ; (c) case 3, the link has four receivers with R L 1 = R L 2 = R L 3 = R L 4 = R L . The figure compares full-wave and circuital simulation results obtained for the case of purely resistive loads and for the case of loads with the compensating inductances given in Table 1, Table 2 and Table 3.
Energies 13 05157 g007
Figure 8. Total output power as function of the generator impedance R G corresponding to the optimal loads Z L M , n = R L M , n + j X L M , n ; results obtained with and without the compensating reactances X L M , n . (a) Case 1; (b) Case 2; (c) Case 3.
Figure 8. Total output power as function of the generator impedance R G corresponding to the optimal loads Z L M , n = R L M , n + j X L M , n ; results obtained with and without the compensating reactances X L M , n . (a) Case 1; (b) Case 2; (c) Case 3.
Energies 13 05157 g008
Figure 9. Case 4: WPT link analyzed through full-wave simulations. The link has a single transmitter and two receivers.
Figure 9. Case 4: WPT link analyzed through full-wave simulations. The link has a single transmitter and two receivers.
Energies 13 05157 g009
Figure 10. Power gain calculated through circuital simulations for case 4. Results obtained by varying the resistive part of the loads when the reactive parts are set according to the optimal values provided by the theory (see Table 4).
Figure 10. Power gain calculated through circuital simulations for case 4. Results obtained by varying the resistive part of the loads when the reactive parts are set according to the optimal values provided by the theory (see Table 4).
Energies 13 05157 g010
Figure 11. Total output power calculated for case 4 by varying the generator impedance when the loads are those maximizing G P .
Figure 11. Total output power calculated for case 4 by varying the generator impedance when the loads are those maximizing G P .
Energies 13 05157 g011
Table 1. Parameters of the equivalent circuit and optimal loads of the WPT link illustrated in Figure 5a (Case 1).
Table 1. Parameters of the equivalent circuit and optimal loads of the WPT link illustrated in Figure 5a (Case 1).
L n , ( n = 0 , 1 , 2 ) C n , ( n = 0 , 1 , 2 ) Q f 0
( n H )( p F ) ( M Hz )
20.914.84557500
Coupling coefficients
k 01 k 02 k 12
0.02150.0215−0.0399
Optimal loads
α G M R G R L M , n , ( n = 1 , 2 ) L L M , n , ( n = 1 , 2 )
( Ω )( Ω )( n H )
16.940.891.9991.9990.835
Table 2. Parameters of the equivalent circuit and optimal loads of the WPT link illustrated in Figure 5b (Case 2).
Table 2. Parameters of the equivalent circuit and optimal loads of the WPT link illustrated in Figure 5b (Case 2).
L n , ( n = 0 , 1 , 2 , 3 ) C n , ( n = 0 , 1 , 2 , 3 ) Q f 0
( n H )( p F ) ( M Hz )
20.914.84557500
Coupling coefficients
k 01 k 02 k 03 k 12 k 13 k 23
0.04060.04060.0178−0.0281−0.0283−0.0283
Optimal loads
α G M R G R L M , n , ( n = 1 , 2 , 3 ) L L M , n , ( n = 1 , 2 ) L L M , 3
( Ω )( Ω )( n H )( n H )
33.4020.9423.9413.9410.8452.695
Table 3. Parameters of the equivalent circuit and optimal loads of the WPT link illustrated in Figure 5c (Case 3).
Table 3. Parameters of the equivalent circuit and optimal loads of the WPT link illustrated in Figure 5c (Case 3).
L n , ( n = 0 , 1 , 2 , 3 , 4 ) C n , ( n = 0 , 1 , 2 , 3 , 4 ) Q f 0
( n H )( p F ) ( M Hz )
20.914.84557500
Coupling coefficients
k 01 k 02 k 03 k 04 k 12 k 13 k 14 k 23 k 24 k 34
0.01230.0123−0.0209−0.0209−0.03950.00840.03930.03930.0084−0.0395
Optimal loads
α G M R G R L M , n , ( n = 1 , 2 , 3 , 4 ) L L M , n , ( n = 1 , 2 ) L L M , j , ( j = 3 , 4 )
( Ω )( Ω )( n H )( n H )
19.1220.9012.2562.2562.5221.414
Table 4. Parameters of the equivalent circuit and optimal loads of the WPT link illustrated in Figure 9 (Case 4).
Table 4. Parameters of the equivalent circuit and optimal loads of the WPT link illustrated in Figure 9 (Case 4).
L 0 L 1 L 2 C 0 C 1 C 2 Q 0 Q 1 Q 2 f 0
( n H )( n H )( n H )( p F )( p F )( p F ) ( M Hz )
45.3828.6720.912.233.534.84622621557500
Coupling coefficients
k 01 k 02 k 12
−0.0890.08940.0367
Optimal loads
α G M R G R L M , 1 R L M , 2 L L M , 1 L L M , 2
( Ω )( Ω )( Ω )( n H )( n H )
76.380.97417.4911.079.010.9480.852

Share and Cite

MDPI and ACS Style

Monti, G.; Mongiardo, M.; Minnaert, B.; Costanzo, A.; Tarricone, L. Optimal Terminations for a Single-Input Multiple-Output Resonant Inductive WPT Link. Energies 2020, 13, 5157. https://doi.org/10.3390/en13195157

AMA Style

Monti G, Mongiardo M, Minnaert B, Costanzo A, Tarricone L. Optimal Terminations for a Single-Input Multiple-Output Resonant Inductive WPT Link. Energies. 2020; 13(19):5157. https://doi.org/10.3390/en13195157

Chicago/Turabian Style

Monti, Giuseppina, Mauro Mongiardo, Ben Minnaert, Alessandra Costanzo, and Luciano Tarricone. 2020. "Optimal Terminations for a Single-Input Multiple-Output Resonant Inductive WPT Link" Energies 13, no. 19: 5157. https://doi.org/10.3390/en13195157

APA Style

Monti, G., Mongiardo, M., Minnaert, B., Costanzo, A., & Tarricone, L. (2020). Optimal Terminations for a Single-Input Multiple-Output Resonant Inductive WPT Link. Energies, 13(19), 5157. https://doi.org/10.3390/en13195157

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