Comparison of the Efﬁciency and Load Power in Periodic Wireless Power Transfer Systems with Circular and Square Planar Coils

: In the article, a wireless charging system with the use of periodically arranged planar coils is presented. The efﬁciency of two wireless power transfer (WPT) systems with different types of inductors, i.e., circular and square planar coils is compared, and two models are proposed: analytical and numerical. With the appropriate selection of a load resistance, it is possible to obtain either the maximum efﬁciency or the maximum power of a receiver. Therefore, the system is analyzed at two optimum modes of operation: with the maximum possible efﬁciency and with the highest power transmitted to the load. The analysis of many variants of the proposed wireless power transfer solution was performed. The aim was to check the inﬂuence of the geometry of the coils and their type (circular or square) on the efﬁciency of the system. Changes in the number of turns, the distance between the coils (transmit and receive) as well as frequency are also taken into account. The results obtained from analytical and numerical analysis were consistent; thus, the correctness of the adopted circuit and numerical model (with periodic boundary conditions) was conﬁrmed. The proposed circuit model and the presented numerical approach allow for a quick estimate of the electrical parameters of the wireless power transmission system. The proposed system can be used to charge many receivers, e.g., electrical cars on a parking or several electronic devices. Based on the results, it was found that the square coils provide lower load power and efﬁciency than compared to circular coils in the entire frequency range and regardless of the analyzed geometry variants. The results and discussion of the multivariate analysis allow for a better understanding of the inﬂuence of the coil geometry on the charging effectiveness. They can also be valuable knowledge when designing this type of system. imitate an inﬁnite dielectric background. The model was made in the Comsol Multiphysics software. We have performed an analysis in the frequency domain, utilizing Magnetic ﬁelds physics in conjugation with the fragmentary Equivalent circuit . The were modeled by using a built-in current sheet approximation of the planar inductors ( Multi-turn coil boundary condition), while the lumped voltage source ( U t ), capacitors ( C ) and load ( Z l ) were attached to the coils by internal coupling with fragmentary Equivalent circuit . NDOF varied for different cases from 200.000 to 400.000 degrees of freedom.


Introduction
Power in the form of electricity generated in power plants is expensive, and the overall efficiency of the transmission from plants to end users is low (about 30%), in particular, due to energy losses in parasitic resistances. The power plant itself needs many resources to produce electricity and has many troubles related with environmental issues (e.g., pollution). Additionally, insufficient reactive power compensation, voltage asymmetry and current distortion reduce the efficiency and quality of electricity transferred to consumers [1]. Indicated problems were one of many, which have resulted in alternate solutions connected with power transmission and distribution. The other problems include cables and wires for supplying electrical devices, e.g., at home, workplace and school, since these are lossy elements that are prone to failure and inconvenient to install. Furthermore, devices using batteries can be disturbed and may stop working during battery replacement.
Considering the above, the high-efficiency wireless charging system would be beneficial in both economic and even social approach. Many entities are trying to explore the use of electric vehicles in order to reduce fuel and energy consumption [2,3]. Moreover, the wireless power transfer (WPT) system could result in huge benefits in the medical field in terms of charging or supplying medical devices (e.g., pacemaker and endoscopy) [4,5]. The focus of researchers working on WPT is to develop and research wireless powering equipment for different industries, e.g., medical, factory and private customer.
Wireless power transmission (WPT) has become one of the most important research focuses in this century. Portability is the main motivation for WPT as the number of portable devices is enormously increasing, and wired chargers will limit their portability. Several review works [6][7][8] described the WPT theory and its applications. Still, developments in the resonant coupled system [9][10][11][12], in the influence of couplings and frequency on efficiency and in human exhibition aspects [13,14] have not been widely considered.
The near-field inductive power transfer (IPT) [15,16] has gained popularity for charging mobile consumer devices. The IPT system can effectively transmit power from a source to the device by using the rule of electromagnetic induction. The advantage of this system is the lack of radiation. Inductive devices (e.g., charging pad for cell phone) work according to the same principles.
After the first experiment of the glowing bulb by using resonant coils conducted by MIT (Massachusetts Institute of Technology, Cambridge, MA, USA) in 2007, progress in resonant coupled system has occurred [17]. In 2008, Intel studied the resonant coupled by using surface coils, which are simpler to montage into mobile devices than the helix coils [18]. In [19], a progressive contactless solution for synchronous powering relative to a few receiving devices was discussed. The research studies on IPT have indicated that, for practical applications, the frequency range from 20 kHz to 200 MHz has been considered [20]. Generally speaking, the bandwidth of WPT should have been limited to industrial, scientific and medical (ISM) frequencies, e.g., 13.56 MHz [6]. However, these frequencies are higher than 6.78 MHz, which will cause additional skin and proximity effect losses in copper conductors wider than 0.05 mm (based on a round wire skin depth). Therefore, IPT systems were also designed to operate at lower frequencies, such as wireless electric vehicles power supply at 20.15 kHz [2] and charging of implantable medical sensors (from 300 kHz to 10 MHz) [4]. Hence, it can be acceptable to use WPT at the kHz range if the low-power system will not exceed electromagnetic safety limits. An example of such system may be also periodic WPT, proposed in this article, which can efficiently operate below 1 MHz since multi-turn coils with relatively high inductances are used.
In opposition to the traditional 2-coil approach [21], the 4-coil approach is designed by setting two indirect multi-turn coils between two loop coils. Every loop coil is a mechanism in which impedance matching takes place and works as a non-resonator to exchange energy between circuits and indirect coils [18]. Further works [22,23] related to multi-coil resonant coupling created interest in the 4-coil system due to increased efficiency. The disadvantage of this approach is it requires more space than any other transmission system.
WPT for many receivers using a single source coil has been described and analyzed in [24]. The downside to this approach is that the resonant frequency of the coils is divided when two receivers are close to each other, which reduces the efficiency. In the article [9], the systems of the coils as domino form resonators are presented. The authors of [10] considered linear resonator systems. These systems included several additional resonators located in the middle (between transmitter and receiver). Their task was to support energy transfer in the system. The analysis was conducted for many resonator configurations. However, parallel-series systems containing planar coils are still not completely developed.
The power transfer efficiency falls at any coupling larger or smaller than its critical value. For this reason, it is necessary to use an additional system to maintain the high efficiency. In this approach, the mentioned system must work without movable resonant frequency at coupling length changes. Different approaches (e.g., adjusting resonant parameters) were presented in [25][26][27][28][29][30][31]. Their goal was to transfer power with better efficiency in the resonant coupled system. The influence of axial and angular misalignment of the coils was presented in [32]. This is the main problem in WPT realization for mobile devices.
The power supply of many devices situated close to one another can be realized by using a net of periodically located coils. This approach raises the density of transferred power and also provides possibilities of synchronous energy supply for many devices. Proposed solutions can be used to power either one or multiple independent loads and, in some cases, replace conventional IPT systems. The developed periodic WPT system allows for the simultaneous supply/charging of many low-power receivers, such as mobile devices or sensors repeatedly distributed over hard-to-reach areas.
The article presents the wireless charging system with periodically arranged planar coils. The proposed analysis of the unit cell with periodic boundary conditions does not demand the creation of a 3D model with many coils [33]. This model has many degrees of freedom. The simplified model of equivalent circuit is an alternative for large matrix formulations [9,10,34]. The main goal of this work is to present and analyze WPT systems, including circular or square coils. The proposed system may be used for analyzing power transfer conditions. Two analysis approaches have been proposed: analytical and numerical. An analytical model, based on equivalent circuit of the coupled resonators, was developed to solve the periodic system. This approach resulted in the quick determination of the crucial parameters of the system (e.g., transmitter power, load power and efficiency), without the need to create large and complex numerical models. In addition, the numerical model of the proposed periodic WPT system was presented. This approach takes into account the periodic distribution of the coils and their diversified structure.
The results were consistent, which confirmed the correctness of the assumptions made. These methods can be used as alternatives to experimental prototypes and simplified analytical-empirical models currently used for analyzing the electrical and magnetic properties of WPT systems. Both approaches reduce the complexity of the commonly used numerical and analytical models. By the proposed appropriate selection of the load resistance, it was possible to determine either the maximum efficiency of the WPT system or the maximum power of the receiver. Therefore, the system was considered at two optimum modes of operation: with the maximum possible efficiency and with the highest power transmitted to the load. The calculations of exemplary periodic WPT systems were performed over a frequency range from 0.1 MHz to 1 MHz. The analysis was multi-variant because the influence of geometric parameters (type of the coil, coil radius, number of turns and distance between the coils) on the power transfer efficiency was analyzed.

Description of the Analyzed Models of Wireless Power Transfer System
The analyzed WPT system is composed of periodically arranged pairs of transmitting (Tx) and receiving (Rx) coils. Each pair is called the WPT cell with outer dimensions d × d (Figure 1). Both planar coils in WPT cell are identical and possess the same radius (r) and number of turns (n) wounded around a dielectric carcass. Each coil has an additional compensating capacitor. The distribution of WPT cells results in a transmitting surface consisting of many Tx and sinusoidal voltage source (U t ), while the Rx coils of a receiving surface are connected to loads (Z l ). Two models of WPT system were analyzed. The first one was made of circular coils with smaller (r = 10 mm) and larger (r = 25 mm) radius ( Figure 2a). The second system had identical radii but consisted of square planar coils (Figure 2b). The compensating capacitors were connected in series with the coils. In the two proposed WPT models, it is possible to select the power conditions in each cell depending on the imposed requirements, e.g., to supply multiple independent receivers simultaneously. The periodic WPT structure consists of two main surfaces: transmitting and receiving, presented in Figure 3. Each surface (composed of Tx or Rx coils) includes a group of planar coils with the same winding direction. The considered cell Ax+j,y+k is an element of an array with identical inductors, where j is the number of column and k is the number of row in a grid (j, k∈N and N is the group of integers). Adjacent coils of the element Ax,y (e.g., Ax+1,y or Ax+1,y+1) are separated by the distance d where d ≈ 2r. Two models of WPT system were analyzed. The first one was made of circular coils with smaller (r = 10 mm) and larger (r = 25 mm) radius ( Figure 2a). The second system had identical radii but consisted of square planar coils (Figure 2b). The compensating capacitors were connected in series with the coils. Two models of WPT system were analyzed. The first one was made of circular coils with smaller (r = 10 mm) and larger (r = 25 mm) radius ( Figure 2a). The second system had identical radii but consisted of square planar coils (Figure 2b). The compensating capacitors were connected in series with the coils. In the two proposed WPT models, it is possible to select the power conditions in each cell depending on the imposed requirements, e.g., to supply multiple independent receivers simultaneously. The periodic WPT structure consists of two main surfaces: transmitting and receiving, presented in Figure 3. Each surface (composed of Tx or Rx coils) includes a group of planar coils with the same winding direction. The considered cell Ax+j,y+k is an element of an array with identical inductors, where j is the number of column and k is the number of row in a grid (j, k∈N and N is the group of integers). Adjacent coils of the element Ax,y (e.g., Ax+1,y or Ax+1,y+1) are separated by the distance d where d ≈ 2r. In the two proposed WPT models, it is possible to select the power conditions in each cell depending on the imposed requirements, e.g., to supply multiple independent receivers simultaneously. The periodic WPT structure consists of two main surfaces: transmitting and receiving, presented in Figure 3. Each surface (composed of Tx or Rx coils) includes a group of planar coils with the same winding direction. The considered cell A x+j,y+k is an element of an array with identical inductors, where j is the number of column and k is the number of row in a grid (j, k ∈ N and N is the group of integers). Adjacent coils of the element A x,y (e.g., A x+1,y or A x+1,y+1 ) are separated by the distance d where d ≈ 2r. Energies 2021, 14, x FOR PEER REVIEW 5 of 24 (a) (b) Figure 3. A two-dimensional fragment of the transmitting/receiving surface group of WPT system for the following: (a) circular coils (Ax,y-the WPT cell) and (b) square coils.
Every transmitter Tx is connected in parallel with a sinusoidal voltage source (Ut). The coils creating the receiving surface are linked to a separate load (Zl). The discussed models of the WPT system guarantee an increase in density of transmitted power in the area between the receiving and transmitting surfaces. This approach also allows the selection of power conditions depending on the requirements.

Numerical Approach to the Analysis of the Periodic WPT System
The proposed WPT systems were analyzed by using two approaches: 3D numerical model and equivalent circuit. The usage of simpler circuit model renders it possible to easily determine power flow at the design stage or initial analysis of the WPT system. Numerical methods (e.g., FEM, FDTD and FDFD) [35][36][37] allow for the creation of a complex model and determination of the magnetic field distribution. In this typical approach, it is necessary to prepare the 3D model and to impose suitable boundary conditions. However, the efficiency and accuracy of the solution depends on the model size (number of degrees of freedom, NDOF). Increasing NDOF will result in greater accuracy of the solution but also longer computation time. For an open boundary problem, the compromise between numerical methods and circuit approaches is the partial equivalent electrical circuit method (PEEC) [38,39]. It permits modeling geometry and analyzing electric and magnetic phenomena, however, NDOF is significantly reduced due to partitioning of the 3D structure by using larger and, thus, fewer finite elements, which become the coupled nodes in a peculiar kind of electrical circuit [39]. Then, lumped parameters, e.g., self-inductance, mutual inductance and resistance are calculated by an application of integral form of Maxwell's equations instead of their differential form typical for FDTD or FEM.
The proposed numerical analysis of the system composed of many WPT cells requires taking into account the following: coil geometry (size and shape), number of turns (n), number of WPT cells and elements of the electrical circuit connected to each coil. The wounds of coils can be made of an ultra-thin wire with diameter (w) and covered with insulation of a thickness (i). A compensating capacitor can be modeled as an element with lumped capacity (C) in series with an inductor. The voltage source with given amplitude (Ut) and specified frequency (f) is connected with Tx coil, while the load (Zl) is connected with Rx coil.
In order to reduce the model to the single cell Ax,y filled with air and containing only a pair of Tx and Rx coils (Figure 4), an infinite array of resonators was modeled by using periodic boundary conditions (PBC) [36,37] imposed on the four lateral surfaces of the Every transmitter Tx is connected in parallel with a sinusoidal voltage source (U t ). The coils creating the receiving surface are linked to a separate load (Z l ). The discussed models of the WPT system guarantee an increase in density of transmitted power in the area between the receiving and transmitting surfaces. This approach also allows the selection of power conditions depending on the requirements.

Numerical Approach to the Analysis of the Periodic WPT System
The proposed WPT systems were analyzed by using two approaches: 3D numerical model and equivalent circuit. The usage of simpler circuit model renders it possible to easily determine power flow at the design stage or initial analysis of the WPT system. Numerical methods (e.g., FEM, FDTD and FDFD) [35][36][37] allow for the creation of a complex model and determination of the magnetic field distribution. In this typical approach, it is necessary to prepare the 3D model and to impose suitable boundary conditions. However, the efficiency and accuracy of the solution depends on the model size (number of degrees of freedom, NDOF). Increasing NDOF will result in greater accuracy of the solution but also longer computation time. For an open boundary problem, the compromise between numerical methods and circuit approaches is the partial equivalent electrical circuit method (PEEC) [38,39]. It permits modeling geometry and analyzing electric and magnetic phenomena, however, NDOF is significantly reduced due to partitioning of the 3D structure by using larger and, thus, fewer finite elements, which become the coupled nodes in a peculiar kind of electrical circuit [39]. Then, lumped parameters, e.g., selfinductance, mutual inductance and resistance are calculated by an application of integral form of Maxwell's equations instead of their differential form typical for FDTD or FEM.
The proposed numerical analysis of the system composed of many WPT cells requires taking into account the following: coil geometry (size and shape), number of turns (n), number of WPT cells and elements of the electrical circuit connected to each coil. The wounds of coils can be made of an ultra-thin wire with diameter (w) and covered with insulation of a thickness (i). A compensating capacitor can be modeled as an element with lumped capacity (C) in series with an inductor. The voltage source with given amplitude (U t ) and specified frequency (f ) is connected with Tx coil, while the load (Z l ) is connected with Rx coil.
In order to reduce the model to the single cell A x,y filled with air and containing only a pair of Tx and Rx coils (Figure 4), an infinite array of resonators was modeled by using periodic boundary conditions (PBC) [36,37] imposed on the four lateral surfaces of the cell. A perfectly matched layer (PML) was placed at the top and bottom of the model to imitate an infinite dielectric background. The model was made in the Comsol Multiphysics software. We have performed an analysis in the frequency domain, utilizing Magnetic fields physics in conjugation with the fragmentary Equivalent circuit. The coils were modeled by using a built-in current sheet approximation of the planar inductors (Multi-turn coil boundary condition), while the lumped voltage source (U t ), capacitors (C) and load (Z l ) were attached to the coils by internal coupling with fragmentary Equivalent circuit. NDOF varied for different cases from 200.000 to 400.000 degrees of freedom.  The energy transport in the presented model and the magnetic field phenomena in the frequency domain can be solved by using the Helmholtz equation: where A denotes magnetic vector potential (Wb/m), ω denotes angular frequency (rad/s), μ0 denotes permeability of an air (H/m), σ denotes electrical conductivity (S/m), and Jext denotes external current density (A/m 2 ). Periodic boundary conditions on four external surfaces were defined as magnetic insulation: where nsurf is the surface normal vector. The external current density (Jext) results from the value of the voltage supply (Ut) connected to the transmitting coil. Equation (1) with Equation (2) and PML conditions can be solved by using the Finite Element Method (FEM) [37]. The aim is to find the unknown spatial distribution of the magnetic vector potential A in 3D Cartesian space at specified frequency f.

Analytical Approach to the Analysis of WPT System
Despite the availability of the computational units, it is always a more difficult task to make and solve the numerical model than the equivalent circuit in which there is no need to use boundary conditions. Therefore, in many cases, a simpler model is more desirable if it still provides a similar scope of analysis but less complexity and much faster modeling and calculation processes. Based on the proposed structure of WPT system, two grids (transmitting and receiving surface) with an infinite number of identical resonators have to be considered. In the electrical circuit, the analysis of the infinite periodic network can be reduced to the single cell Ax,y. The energy transport in the presented model and the magnetic field phenomena in the frequency domain can be solved by using the Helmholtz equation: where A denotes magnetic vector potential (Wb/m), ω denotes angular frequency (rad/s), µ 0 denotes permeability of an air (H/m), σ denotes electrical conductivity (S/m), and J ext denotes external current density (A/m 2 ). Periodic boundary conditions on four external surfaces were defined as magnetic insulation: where n surf is the surface normal vector. The external current density (J ext ) results from the value of the voltage supply (U t ) connected to the transmitting coil. Equation (1) with Equation (2) and PML conditions can be solved by using the Finite Element Method (FEM) [37]. The aim is to find the unknown spatial distribution of the magnetic vector potential A in 3D Cartesian space at specified frequency f.

Analytical Approach to the Analysis of WPT System
Despite the availability of the computational units, it is always a more difficult task to make and solve the numerical model than the equivalent circuit in which there is no need to use boundary conditions. Therefore, in many cases, a simpler model is more desirable if it still provides a similar scope of analysis but less complexity and much faster modeling and calculation processes. Based on the proposed structure of WPT system, two grids (transmitting and receiving surface) with an infinite number of identical resonators have to be considered. In the electrical circuit, the analysis of the infinite periodic network can be reduced to the single cell A x,y .
The system included identical coils stacked coaxially at the distance h, and they are made of thin wire with electrical conductivity σ. The length of all windings in the circular coil is described by Equation (3), while of the square coil is described by Equation (4) If Tx and Rx coils are identical, then their resistances are equal; thus, R t = R r = R c . Taking into account Equation (3), the formula for the resistance R c of circular coil is as follows.
Based on Equation (4), the R c of the square coil is as follows.
In an infinite periodic grid where each coil has identical electrical parameters and magnetic coupling occurs with its neighbors, it is possible to reduce an analysis to the single WPT cell ( Figure 5). Mutual inductances between coil in the cell A x,y and coils in cells A x+i,y+j affect the internal inductance of the coil in A x,y , which can be expressed for both types of the coils as follows.
The following components of Equation (7) are as follows: L c -effective inductance in (H); M x+i,y+j -mutual inductance in (H) between coils adjacent in horizontal plane, for i = 0 and j = 0; L self -self-inductance of the planar coil in (H). The system included identical coils stacked coaxially at the distance h, and they are made of thin wire with electrical conductivity σ. The length of all windings in the circular coil is described by Equation (3), while of the square coil is described by Equation (4) ( )( ) If Tx and Rx coils are identical, then their resistances are equal; thus, Rt = Rr = Rc. Taking into account Equation (3), the formula for the resistance Rc of circular coil is as follows.
( )( ) Based on Equation (4), the Rc of the square coil is as follows.
In an infinite periodic grid where each coil has identical electrical parameters and magnetic coupling occurs with its neighbors, it is possible to reduce an analysis to the single WPT cell ( Figure 5). Mutual inductances between coil in the cell Ax,y and coils in cells Ax+i,y+j affect the internal inductance of the coil in Ax,y, which can be expressed for both types of the coils as follows.
( ) The following components of Equation (7) are as follows: Lc-effective inductance in (H); Mx+i,y+j-mutual inductance in (H) between coils adjacent in horizontal plane, for i ≠ 0 and j ≠ 0; Lself-self-inductance of the planar coil in (H). Self-inductance is calculated using equation [40]: where dm is mean diameter and ν is a fill coefficient. Self-inductance is calculated using equation [40]: where d m is mean diameter and ν is a fill coefficient.
In Equation (8) coefficients c 1 , c 2 , c 3 and c 4 are dependent on the shape of the coil [40]. The values of the coefficients for circular and square inductors are presented in Table 1. Table 1. Geometrical coefficients for circular and square coils [40]. From Equation (7), an inductance of the considered coil in a segment A x,y takes the following form:

Type of Coil
where M pe is sum of mutual inductances in the periodic grid. Since these mutual inductances reduce the effective inductance of the coil, the inductance M pe is written with a minus. For the case when loads Z l = ∞ and there is no capacitor in series with Tx coils, at the arbitrary frequency f, one may find M pe as follows [11]: where I t,∞ =|I t,∞ |e jψ is the source current in (A); |I t,∞ | is the RMS value of the source current in (A); ψ is the phase angle between the source voltage and current in (rad). The underlined quantities, e.g., the current I t,∞ in Equation (12), denotes the complex number. According to the equivalent circuit shown in Figure 5b, instead of calculating inductances M tr,x+i,y+j , the mutual inductance M tr between Tx and Rx can be found as [11].
where U r,∞ =|U r,∞ |e jθ is the voltage induced in the receiving coil in (V), |U r,∞ | is the RMS value of the induced voltage in (V) and θ is the phase angle between the source voltage and induced voltage in (rad). After calculations of self-inductance (L self ), mutual inductance of periodically adjacent coils (M pe ) and mutual inductance between Tx and Rx (M tr ), it is possible to find the compensating capacity (C) at a specified frequency.

Analyzed Models
By the utilization of the proposed analytical and numerical models, many variants of the periodic WPT system were analyzed ( Table 2). The shapes of the coils (circular or square), radius (r), number of turns (n) and distance between Tx and Rx coils (h) were changed.  Based on the results obtained for several exemplary periodic WPT systems, the correctness of the proposed analytical model was verified by comparing the active power of the receiver: where I r is the current flowing through receiving coil in (A). The power of the transmitter is represented by the following: where I t is the current flowing through transmitting coil in (A). By using Equations (15) and (16), the power transfer efficiency was calculated as follows.
The comparison concerned the determination of the powers P o and P z and also efficiency η. Additionally, the analysis was divided into two cases: -Operation with maximum efficiency, where the optimal load impedance was described as follows.
-Operation with maximum transferred power, where the optimal load impedance was described as follows.
The discussed unit cell A x,y was subjected to analysis where it was assumed that the system was made of the infinite number of WPT cells. Each WPT cell consisted of pair of identical coaxial coils (Tx and Rx). The parameters regarding windings (wires) of the coils, which were used in models, are presented in Table 3. The analysis was performed in the frequency domain for the frequencies from f min = 100 kHz to f max = 1000 kHz. The amplitude of the voltage source was set to U t = 1 V. Table 3. Parameters of the wire used to form the coils.

Parameter Symbol Value
Diameter of the wire w 200 µm Thickness of the wire insulation i 5 µm Conductivity of the wire (copper) σ 5.6 × 10 7 S/m The numerical model was created in the Comsol Multiphysics software and was solved by using FEM. We utilized built-in, multi-turn coils approximation (current sheet approximation) and the electrical circuit connected with the 3D model. The lumped parameters of the electrical circuit from Figure 5b (Tables 4-7) for circular and square coils were found using preliminary numerical calculations and Equations (8), (13), (14), (18) and (19). Table 4. Self and mutual inductances as well as compensating capacities and load impedance (at f max ) for circular coils and r = 10 mm.  Table 6. Self and mutual inductances as well as compensating capacities and load impedance (at f max ) for square coils and r = 10 mm.

System Operating with Maximum Efficiency
The calculation results of the proposed WPT system, received by the numerical and analytical method, were compared. The structure of the WPT system was taken into account (radius and shape of the coil, number of turns and distance between coils), and the transmitter and receiver powers were calculated based on both analytical and numerical models.
The first operating mode was related with maximum efficiency. Impedances of the load (Z e ) were calculated based on Equation (18) for different number of turns and distance between Tx and Rx and presented in Figure 6 for the coil with r = 10 mm and in Figure 7 for the coil with r = 25 mm. The analysis showed that the values of impedances (Z e ) were always higher for the circular coil. An increase in the number of turns causes a nonlinear increase in the impedance. Doubling the distance h results in an over four times decrease in the impedance for both coil sizes and types of the coils.

System Operating with Maximum Efficiency
The calculation results of the proposed WPT system, received by the numerical and analytical method, were compared. The structure of the WPT system was taken into account (radius and shape of the coil, number of turns and distance between coils), and the transmitter and receiver powers were calculated based on both analytical and numerical models.
The first operating mode was related with maximum efficiency. Impedances of the load (Ze) were calculated based on Equation (18) for different number of turns and distance between Tx and Rx and presented in Figure 6 for the coil with r = 10 mm and in Figure 7 for the coil with r = 25 mm. The analysis showed that the values of impedances (Ze) were always higher for the circular coil. An increase in the number of turns causes a nonlinear increase in the impedance. Doubling the distance h results in an over four times decrease in the impedance for both coil sizes and types of the coils. In Table 8, the maximum power transfer efficiency values depending on the type of the coils, the number of turns and the distance between Tx and Rx were compared. In each case, higher power transfer efficiency values were obtained for the circular coil than for the square one.
In order to compare the results obtained with the analytical and numerical methods, the calculated absolute error is shown in Table 9. The results are consistent and do not exceed 0.4%.

System Operating with Maximum Efficiency
The calculation results of the proposed WPT system, received by the numerical and analytical method, were compared. The structure of the WPT system was taken into account (radius and shape of the coil, number of turns and distance between coils), and the transmitter and receiver powers were calculated based on both analytical and numerical models.
The first operating mode was related with maximum efficiency. Impedances of the load (Ze) were calculated based on Equation (18) for different number of turns and distance between Tx and Rx and presented in Figure 6 for the coil with r = 10 mm and in Figure 7 for the coil with r = 25 mm. The analysis showed that the values of impedances (Ze) were always higher for the circular coil. An increase in the number of turns causes a nonlinear increase in the impedance. Doubling the distance h results in an over four times decrease in the impedance for both coil sizes and types of the coils. In Table 8, the maximum power transfer efficiency values depending on the type of the coils, the number of turns and the distance between Tx and Rx were compared. In each case, higher power transfer efficiency values were obtained for the circular coil than for the square one.
In order to compare the results obtained with the analytical and numerical methods, the calculated absolute error is shown in Table 9. The results are consistent and do not exceed 0.4%.
n [-] h=25 mm (sqr) h=25 mm (cir) Figure 7. Load impedance at f max for obtaining a maximum efficiency for coil with r = 25 mm at the following distances: In Table 8, the maximum power transfer efficiency values depending on the type of the coils, the number of turns and the distance between Tx and Rx were compared. In each case, higher power transfer efficiency values were obtained for the circular coil than for the square one.
In order to compare the results obtained with the analytical and numerical methods, the calculated absolute error is shown in Table 9. The results are consistent and do not exceed 0.4%. The exact characteristics of the source power, receiver power and power transfer efficiency are presented on the graph in the following subsections: small cell (r = 10 mm) in Section 3.2.1 and large cell (r = 25 mm) in Section 3.2.2. The results of the numerical analysis are shown as dashed lines for circular coils and solid lines for square coils. The results of the analytical analysis, based on the previously presented formulas (Section 2.3), are marked with dots for each type of the coil and the distance between Tx and Rx.
3.2.1. Maximum Power Transfer Efficiency for Small Coil (r = 10 mm)  show the characteristics of the following: power (transmitter and receiver) and efficiency at the distance h = r/2 = 5 mm and h = r = 10 mm between coils for circular and square planar spiral coils.  The exact characteristics of the source power, receiver power and power transfer efficiency are presented on the graph in the following subsections: small cell (r = 10 mm) in Section 3.2.1 and large cell (r = 25 mm) in Section 3.2.2. The results of the numerical analysis are shown as dashed lines for circular coils and solid lines for square coils. The results of the analytical analysis, based on the previously presented formulas (Section 2.3), are marked with dots for each type of the coil and the distance between Tx and Rx.    The transmitter power (Pz) decreases over the entire frequency range, regardless of the type and size of the coil and the number of turns (Figures 8a, 9a, 10a and 11a). The transmitter power is almost the same for the circular and square coil at the distance h = r/2. For the distance h = r, Pz is greater for the circular coil in the initial frequency range (from 100 kHz to about 450 kHz), and Pz values are greater for the square coil with increasing frequency.
For all characteristics (Figures 8b, 9b, 10b and 11b), it can be noted that when the efficiency reaches approximately 50%, the receiver power (Po) reaches its maximum value. Then, the receiver power decreases as efficiency exceeded 50%. Despite the nonlinear frequency dependency of power and efficiency, the number of turns also has a significant impact. By increasing the number of turns, the receiver power will fall. The shape of the receiver power characteristics for both types of the coils is very similar. Maximum receiver power is achieved for h = 5 mm in the range 100 kHz ÷ 200 kHz and for h = 10 mm in the range 450 kHz ÷ 800 kHz.
The greatest difference in efficiency between circular and square coils occurs at n = 15 and reaches almost 20% (at f = 1MHz, h = 10 mm) and almost 15% (at f = 300 kHz, h = 5 mm) (Figure 8c). As the number of turns increases, this difference decreases and at n = 45 it does not exceed 8% regardless of the distance between the coils h ( Figure 11c)  The transmitter power (Pz) decreases over the entire frequency range, regardless of the type and size of the coil and the number of turns (Figures 8a, 9a, 10a and 11a). The transmitter power is almost the same for the circular and square coil at the distance h = r/2. For the distance h = r, Pz is greater for the circular coil in the initial frequency range (from 100 kHz to about 450 kHz), and Pz values are greater for the square coil with increasing frequency.
For all characteristics (Figures 8b, 9b, 10b and 11b), it can be noted that when the efficiency reaches approximately 50%, the receiver power (Po) reaches its maximum value. Then, the receiver power decreases as efficiency exceeded 50%. Despite the nonlinear frequency dependency of power and efficiency, the number of turns also has a significant impact. By increasing the number of turns, the receiver power will fall. The shape of the receiver power characteristics for both types of the coils is very similar. Maximum receiver power is achieved for h = 5 mm in the range 100 kHz ÷ 200 kHz and for h = 10 mm in the range 450 kHz ÷ 800 kHz.
The greatest difference in efficiency between circular and square coils occurs at n = 15 and reaches almost 20% (at f = 1MHz, h = 10 mm) and almost 15% (at f = 300 kHz, h = 5 mm) (Figure 8c). As the number of turns increases, this difference decreases and at n = 45 it does not exceed 8% regardless of the distance between the coils h (Figure 11c)  The transmitter power (Pz) decreases over the entire frequency range, regardless of the type and size of the coil and the number of turns (Figures 8a, 9a, 10a and 11a). The transmitter power is almost the same for the circular and square coil at the distance h = r/2. For the distance h = r, Pz is greater for the circular coil in the initial frequency range (from 100 kHz to about 450 kHz), and Pz values are greater for the square coil with increasing frequency.
For all characteristics (Figures 8b, 9b, 10b and 11b), it can be noted that when the efficiency reaches approximately 50%, the receiver power (Po) reaches its maximum value. Then, the receiver power decreases as efficiency exceeded 50%. Despite the nonlinear frequency dependency of power and efficiency, the number of turns also has a significant impact. By increasing the number of turns, the receiver power will fall. The shape of the receiver power characteristics for both types of the coils is very similar. Maximum receiver power is achieved for h = 5 mm in the range 100 kHz ÷ 200 kHz and for h = 10 mm in the range 450 kHz ÷ 800 kHz.
The greatest difference in efficiency between circular and square coils occurs at n = 15 and reaches almost 20% (at f = 1MHz, h = 10 mm) and almost 15% (at f = 300 kHz, h = 5 mm) (Figure 8c). As the number of turns increases, this difference decreases and at n = 45 it does not exceed 8% regardless of the distance between the coils h (Figure 11c)  The transmitter power (P z ) decreases over the entire frequency range, regardless of the type and size of the coil and the number of turns (Figures 8a, 9a, 10a and 11a). The transmitter power is almost the same for the circular and square coil at the distance h = r/2. For the distance h = r, P z is greater for the circular coil in the initial frequency range (from 100 kHz to about 450 kHz), and P z values are greater for the square coil with increasing frequency.
For all characteristics (Figures 8b, 9b, 10b and 11b), it can be noted that when the efficiency reaches approximately 50%, the receiver power (P o ) reaches its maximum value. Then, the receiver power decreases as efficiency exceeded 50%. Despite the nonlinear frequency dependency of power and efficiency, the number of turns also has a significant impact. By increasing the number of turns, the receiver power will fall. The shape of the receiver power characteristics for both types of the coils is very similar. Maximum receiver power is achieved for h = 5 mm in the range 100 kHz ÷ 200 kHz and for h = 10 mm in the range 450 kHz ÷ 800 kHz.
The greatest difference in efficiency between circular and square coils occurs at n = 15 and reaches almost 20% (at f = 1MHz, h = 10 mm) and almost 15% (at f = 300 kHz, h = 5 mm) ( Figure 8c). As the number of turns increases, this difference decreases and at n = 45 it does not exceed 8% regardless of the distance between the coils h (Figure 11c). For h = 5 mm and n = 35, the number of turns n = 45 had minor impact on the efficiency, which reached almost 90% at 1 MHz for both types of the coils (Figures 10c and 11c). For h = 10 mm, efficiency reached 67% for circular coils and 60% for square coils (Figure 11c). Doubling the distance between the coils caused the efficiency to decrease by 20-40%.  (Figures 10c and 11c). For h = 10 mm, efficiency reached 67% for circular coils and 60% for square coils (Figure 11c). Doubling the distance between the coils caused the efficiency to decrease by 20-40%.    (Figures 10c and 11c). For h = 10 mm, efficiency reached 67% for circular coils and 60% for square coils (Figure 11c). Doubling the distance between the coils caused the efficiency to decrease by 20-40%.    (Figures 10c and 11c). For h = 10 mm, efficiency reached 67% for circular coils and 60% for square coils (Figure 11c). Doubling the distance between the coils caused the efficiency to decrease by 20-40%.    Regardless of the type of the coil and the number of turns, the transmitter power (Pz) decreases over the entire frequency range (Figures 12a, 13a, 14a, 15a and 16a). For the large coil (r = 25 mm) and the h = 12.5 mm distance between the coils, the transmitter power is almost the same for the circular and square coil. However, for the distance h = 25 mm, Pz is greater for the circular coil in the initial frequency range, and it is greater for the square coil with increasing frequency. For the distance h = r, Pz is greater for the circular coil in the initial frequency range 100 kHz ÷ 300 kHz, and the Pz values became greater for the square coil with increasing frequency.

Maximum Power
On all characteristics, it can be noted that when the efficiency reaches approximately 50%, the receiver power Po reaches its maximum value. Then, the receiver power decreases as efficiency exceeded 50%. Maximum receiver power is achieved for h = 12.5 mm at f = 100 kHz and for h = 25 mm in the range 200 kHz ÷ 250 kHz. For larger distances between coils h = r = 25 mm, the values of the receiver power for the circular coils are greater than for the square coils up to the frequency of about 700 kHz, and the relationship is inverse.
Regardless of the size of the coils, the number of turns and the distance between the coils, the efficiency is always greater for circular coils than for square ones (Figures 12c,  13c, 14c, 15c and 16c). The efficiency of the analyzed system increases with increasing frequency and the number of turns. Regardless of the number of turns, doubling the distance between the coils caused the efficiency to decrease by 20-50% in lower frequency range and by 10-40% in higher frequency range. The greatest differences in efficiency between the circular and square coils occurred at n = 40 and reached almost 20% (h = 25 mm, f = 400 kHz) and almost 15% (h = 12.5 mm, f = 100 kHz) (Figure 12c). For large distances (h = 25 mm), as the number of turns increases, the difference in efficiency between the circular and square coils becomes smaller.

System Operating with Maximum Load Power
The second operating mode was related with the maximum load power. In order to determine the maximum power transferred to the receiver, the values of load impedance Regardless of the type of the coil and the number of turns, the transmitter power (Pz) decreases over the entire frequency range (Figures 12a, 13a, 14a, 15a and 16a). For the large coil (r = 25 mm) and the h = 12.5 mm distance between the coils, the transmitter power is almost the same for the circular and square coil. However, for the distance h = 25 mm, Pz is greater for the circular coil in the initial frequency range, and it is greater for the square coil with increasing frequency. For the distance h = r, Pz is greater for the circular coil in the initial frequency range 100 kHz ÷ 300 kHz, and the Pz values became greater for the square coil with increasing frequency.
On all characteristics, it can be noted that when the efficiency reaches approximately 50%, the receiver power Po reaches its maximum value. Then, the receiver power decreases as efficiency exceeded 50%. Maximum receiver power is achieved for h = 12.5 mm at f = 100 kHz and for h = 25 mm in the range 200 kHz ÷ 250 kHz. For larger distances between coils h = r = 25 mm, the values of the receiver power for the circular coils are greater than for the square coils up to the frequency of about 700 kHz, and the relationship is inverse.
Regardless of the size of the coils, the number of turns and the distance between the coils, the efficiency is always greater for circular coils than for square ones (Figures 12c,  13c, 14c, 15c and 16c). The efficiency of the analyzed system increases with increasing frequency and the number of turns. Regardless of the number of turns, doubling the distance between the coils caused the efficiency to decrease by 20-50% in lower frequency range and by 10-40% in higher frequency range. The greatest differences in efficiency between the circular and square coils occurred at n = 40 and reached almost 20% (h = 25 mm, f = 400 kHz) and almost 15% (h = 12.5 mm, f = 100 kHz) (Figure 12c). For large distances (h = 25 mm), as the number of turns increases, the difference in efficiency between the circular and square coils becomes smaller.

System Operating with Maximum Load Power
The second operating mode was related with the maximum load power. In order to determine the maximum power transferred to the receiver, the values of load impedance Regardless of the type of the coil and the number of turns, the transmitter power (P z ) decreases over the entire frequency range (Figures 12a, 13a, 14a, 15a and 16a). For the large coil (r = 25 mm) and the h = 12.5 mm distance between the coils, the transmitter power is almost the same for the circular and square coil. However, for the distance h = 25 mm, P z is greater for the circular coil in the initial frequency range, and it is greater for the square coil with increasing frequency. For the distance h = r, P z is greater for the circular coil in the initial frequency range 100 kHz ÷ 300 kHz, and the P z values became greater for the square coil with increasing frequency.
On all characteristics, it can be noted that when the efficiency reaches approximately 50%, the receiver power P o reaches its maximum value. Then, the receiver power decreases as efficiency exceeded 50%. Maximum receiver power is achieved for h = 12.5 mm at f = 100 kHz and for h = 25 mm in the range 200 kHz ÷ 250 kHz. For larger distances between coils h = r = 25 mm, the values of the receiver power for the circular coils are greater than for the square coils up to the frequency of about 700 kHz, and the relationship is inverse.
Regardless of the size of the coils, the number of turns and the distance between the coils, the efficiency is always greater for circular coils than for square ones (Figures 12c, 13c,  14c, 15c and 16c). The efficiency of the analyzed system increases with increasing frequency and the number of turns. Regardless of the number of turns, doubling the distance between the coils caused the efficiency to decrease by 20-50% in lower frequency range and by 10-40% in higher frequency range. The greatest differences in efficiency between the circular and square coils occurred at n = 40 and reached almost 20% (h = 25 mm, f = 400 kHz) and almost 15% (h = 12.5 mm, f = 100 kHz) (Figure 12c). For large distances (h = 25 mm), as the number of turns increases, the difference in efficiency between the circular and square coils becomes smaller.

System Operating with Maximum Load Power
The second operating mode was related with the maximum load power. In order to determine the maximum power transferred to the receiver, the values of load impedance (Z p ) were calculated by taking into account the number of turns, distance between the coils and mutual inductance for circular and square planar spiral coils for the following: r = 10 mm ( Figure 17) and r = 25 mm ( Figure 18). As with the previous approach (Section 3.2), it can be observed that impedance values are greater for circular coils than square coils. In both types of the coils, doubling the distance h results in an over seventeen-times decrease in the impedance (for r = 10 mm) but more than a twenty-times decrease in Z p for square coils (for r = 25 mm). The maximum load power values depending on the type of the coils, the number of turns and the distance (h) are presented in Table 10. In each case, higher power transfer efficiency values were obtained for the circular coil than for the square one. The maximum load power at the receiver is about 20% greater for the circular coil than for the square one. The maximum load power values depending on the type of the coils, the number of turns and the distance (h) are presented in Table 10. In each case, higher power transfer efficiency values were obtained for the circular coil than for the square one. The maximum load power at the receiver is about 20% greater for the circular coil than for the square one. The maximum load power values depending on the type of the coils, the number of turns and the distance (h) are presented in Table 10. In each case, higher power transfer efficiency values were obtained for the circular coil than for the square one. The maximum load power at the receiver is about 20% greater for the circular coil than for the square one. The compared values obtained with proposed circuit and FEM methods were performed by calculation of the absolute error for all frequencies and all cases (Table 11). The results are consistent and do not exceed 2.1 mW.  The compared values obtained with proposed circuit and FEM methods were performed by calculation of the absolute error for all frequencies and all cases (Table 11). The results are consistent and do not exceed 2.1 mW.    The compared values obtained with proposed circuit and FEM methods were performed by calculation of the absolute error for all frequencies and all cases (Table 11). The results are consistent and do not exceed 2.1 mW.             The transmitter power values (Pz) are always greater for circular coils than square coils, regardless of the number of turns and the size of the coils. Pz also decreases with the increasing number of turns (Figures 19a, 20a, 21a and 22a). For the distance between the coils h = r, the transmitter power decreases in the entire frequency range. For the distance between the coils h = 0.5 r, the transmitter power decreases in the initial frequency range 100 kHz ÷ ~500 kHz and then stabilizes.
The receiver power (Po) for h = 0.5 r (Figures 19b, 20b, 21b and 22b) increases with increasing frequency in the range 100 kHz ÷ 500 kHz, and this range decreases with increasing n, e.g., 100 kHz ÷ 300 kHz for n = 45. Then, the receiver power values stabilize when reaching the maximum value. However, for h = r, the receiver power increases in the whole range of analyzed frequencies. Regardless of the number of turns and the distance between the coils, the receiver power is always greater for circular coils than for square ones, e.g., the difference is almost 25% for n = 15 and f = 1 MHz (Figure 19b). The receiver power stabilizes faster for distance h = 0.5 r.  The transmitter power values (Pz) are always greater for circular coils than square coils, regardless of the number of turns and the size of the coils. Pz also decreases with the increasing number of turns (Figures 19a, 20a, 21a and 22a). For the distance between the coils h = r, the transmitter power decreases in the entire frequency range. For the distance between the coils h = 0.5 r, the transmitter power decreases in the initial frequency range 100 kHz ÷ ~500 kHz and then stabilizes.
The receiver power (Po) for h = 0.5 r (Figures 19b, 20b, 21b and 22b) increases with increasing frequency in the range 100 kHz ÷ 500 kHz, and this range decreases with increasing n, e.g., 100 kHz ÷ 300 kHz for n = 45. Then, the receiver power values stabilize when reaching the maximum value. However, for h = r, the receiver power increases in the whole range of analyzed frequencies. Regardless of the number of turns and the distance between the coils, the receiver power is always greater for circular coils than for square ones, e.g., the difference is almost 25% for n = 15 and f = 1 MHz (Figure 19b). The receiver power stabilizes faster for distance h = 0.5 r.  The transmitter power values (Pz) are always greater for circular coils than square coils, regardless of the number of turns and the size of the coils. Pz also decreases with the increasing number of turns (Figures 19a, 20a, 21a and 22a). For the distance between the coils h = r, the transmitter power decreases in the entire frequency range. For the distance between the coils h = 0.5 r, the transmitter power decreases in the initial frequency range 100 kHz ÷ ~500 kHz and then stabilizes.
The receiver power (Po) for h = 0.5 r (Figures 19b, 20b, 21b and 22b) increases with increasing frequency in the range 100 kHz ÷ 500 kHz, and this range decreases with increasing n, e.g., 100 kHz ÷ 300 kHz for n = 45. Then, the receiver power values stabilize when reaching the maximum value. However, for h = r, the receiver power increases in the whole range of analyzed frequencies. Regardless of the number of turns and the distance between the coils, the receiver power is always greater for circular coils than for square ones, e.g., the difference is almost 25% for n = 15 and f = 1 MHz (Figure 19b). The receiver power stabilizes faster for distance h = 0.5 r.  The transmitter power values (P z ) are always greater for circular coils than square coils, regardless of the number of turns and the size of the coils. P z also decreases with the increasing number of turns (Figures 19a, 20a, 21a and 22a). For the distance between the coils h = r, the transmitter power decreases in the entire frequency range. For the distance between the coils h = 0.5 r, the transmitter power decreases in the initial frequency range 100 kHz ÷~500 kHz and then stabilizes.
The receiver power (P o ) for h = 0.5 r (Figures 19b, 20b, 21b and 22b) increases with increasing frequency in the range 100 kHz ÷ 500 kHz, and this range decreases with increasing n, e.g., 100 kHz ÷ 300 kHz for n = 45. Then, the receiver power values stabilize when reaching the maximum value. However, for h = r, the receiver power increases in the whole range of analyzed frequencies. Regardless of the number of turns and the distance between the coils, the receiver power is always greater for circular coils than for square ones, e.g., the difference is almost 25% for n = 15 and f = 1 MHz (Figure 19b). The receiver power stabilizes faster for distance h = 0.5 r.
At the considered operation mode, the efficiency of the system tended to maximize at 50%; then, the maximum power is emitted on the receiver. For the circular coils, the system achieves 50% efficiency at lower frequency compared to the square coils. Regardless of the number of turns and the distance between the coils, the efficiency is always greater for circular coils than for square ones (Figures 19c, 20c, 21c and 22c). For distance h = 5 mm, along with an increase in the number of turns, the efficiency values are comparable for both types of coils, whereas the efficiency is always higher for the circular coils for h = 10 mm.

Results of the Maximum Load Power for Large Coil (r = 25 mm)
The transmitter power (P z ) is always greater for circular coils than for square coils, regardless of the number of turns and the size of the coils. The trend of all P z characteristics is similar to the graphs for smaller coils (r = 10 mm), but the values are lower. P z also decreases when increasing the number of turns (Figures 19a, 20a, 21a and 22a). For the distance between the coils h = r, the transmitter power decreases in the entire frequency range. For the distance between the coils h = 0.5 r, the transmitter power decreases in the initial frequency range and then stabilizes with increasing frequency. As with in the case of smaller coils (r = 10 mm), as the frequency increases, the P z values are similar for the smaller (h = 0.5 r) and larger (h = r) distance between Tx and Rx regardless of the type of coils. The same relationship also occurs for the receiver power (P o ) (Figures 23b, 24b, 25b,  26b and 27b).
For the large coil r = 25 mm (Figures 23b, 24b, 25b, 26b and 27b), the power of the receiver P o increases with increasing frequency, and then it stabilizes when reaching the maximum value, regardless of the distance between the coils. The P o values stabilize faster than for smaller distance between the coils (h = 12.5 mm). As the number of turns increases, the maximum power at the receiver is already achieved at lower frequencies.
Receiver power decreases as the number of turns increases, regardless of the type of the coils. Similarly, as in the case of smaller coils (r = 10 mm), regardless of the number of turns and the distance between the coils, the receiver power is always greater for circular coils than for square ones.
In this approach, the efficiency of the system tended to maximize at 50%; then, the maximum power is emitted on the receiver. Similarly, as in the case of smaller coils (r = 10 mm), the efficiency is greater for the circular coils than for the square ones, regardless of the n and h (Figures 23c, 24c, 25c, 26c and 27c). Moreover, for both the larger coil and circular coil, the system achieves 50% efficiency at a lower frequency compared to the square coils. For both types of the coils at smaller distance (h = 12.5 mm), the system reaches maximum efficiency faster than at h = 25 mm (in the frequency range 100 kHz ÷ 500 kHz).

Conclusions
Most of the modern electronic devices still use cables and plugs for charging. When the battery runs out, these devices become useless and must be recharged. The reality of wireless powering is to have power wherever it is needed and used. If all these efforts bring progress, we can expect wireless power technology to become widespread and have an impact on our daily lives. Now it is possible to use batteryless devices that are powered all the time by the WPT system within their defined area of work. However, research into different and better solutions is still ongoing. Currently, many electronic companies are developing the essential components of WPT in order to accelerate the execution of WPT technology into more cost-limited applications. If all these efforts bring progress, we hope that wireless power technology will become widespread and will have an impact on our daily lives.
In the article, a wireless charging system with the use of periodically arranged planar coils was presented. The discussed system contains periodically arranged planar circular or square coils. The maximum efficiency and maximum load power in the proposed periodic WPT system were solved based on exemplary structures with many magnetic couplings between element inductors. The models used for solving proposed structures, arranged by many circular or square planar coils, were developed. The main goal was to quickly estimate the parameters (e.g., power and efficiency), both with 3D numerical model and by using its equivalent circuit. The use of the electrical circuit represented by single cell of the WPT system eliminated the need to make very complex models, which are solved using numerical methods. In the article, the construction of the 3D numerical model and conditions that allow reducing the complexity of the numerical model was presented.
The presented solutions, which were adopted in the analytical and numerical models, permits the study of the influence of the type of coils, their geometry and the distance between the transmitter and receiver upon power transmission. Calculations were performed over a wide frequency range. Taking into account an appropriate selection of the load resistance, the system was considered at two following optimum modes of operation: with the maximum possible efficiency and with the highest power transmitted to the load. By easy adjustment of the number of turns and increasing the frequency of the current, it was possible to obtain high efficiency in power transmission (about 95% for both types of the coils) for the loads supplied by using the proposed system, without the use of intermediate coils or iron cores. Even at distances equal to the radius of the coil, the system was still able to reach near 85% efficiency for circular coils and near 80% for square ones. In addition, the results have shown that by using the periodically arranged planar coils system, better efficiency can be obtained with the circular coils than with the square ones: (1) For small coil (r = 10 mm):  The results from the equivalent circuit were comparable with results from the numerical model. The maximum absolute error did not exceed 0.4%. The obtained characteristics confirmed the correctness of the assumptions and proved that the analysis of an extensive grid of periodic resonators can be simplified to the single WPT cell and modelled as a two-port circuit with lumped parameters. Both proposed methods of the analysis are suitable for evaluating the properties of the proposed WPT system.