1. Introduction
In addition to the regulations on wireless charging and its commercialization, it is necessary to consider several other important technical issues in the design of the system. Key factors of power electronic systems like energy efficiency and power density are leading to continual increase and have positive environmental, as well as economic impact. In connection to these facts, also the wireless power transfer (WPT) systems will face contradictory requirements imposed on their system components. Here we are talking about power electronic circuitry, coupling elements, compensation capacitors and thermal management, as well as EM (electro-magnetic) shielding system. The regulations of the power semiconductor systems determined by international organizations like, U.S Energy Star, 80 Plus, Climate Savers, German Blue Angel and so forth are completely different from the regulations defined for the WPT operation. It is therefore a complex task to design high-performance WPT system that will comply with all technical and safety regulations at once. Currently there are various studies describing design procedures of WPT coupling coil system [
1,
2,
3,
4,
5], studies related to control strategies for power delivery regulation [
6,
7,
8] and many studies related to power electronic system analysis and evaluation of the operational properties [
9,
10,
11,
12,
13]. Many of them are focusing on the evaluation of the qualitative indicators of WPT systems, whose are efficiency and power density at all but also the amount of delivered power in dependency on distance of the transferred power. When talking about battery charging in automotive industry, then the series-series compensation topology is the best choice [
14,
15,
16]. Typical equivalent circuit of this topology is being shown on
Figure 1, while secondary side is replaced by equivalent load, which considers the effect of secondary AC/DC converter together with battery charging module. As we focus more deeply on this topology, here it must be noticed that commercial use is also an alternative, while biocompatibility is very important issue. Here we are talking about safety standards [
17,
18,
19], which are defining limits for the safe interaction between EM field and living organisms. The ways for keeping safety limits within definitions are based on the EM shielding. The literature proposes several analyses focused on this topic but only a few of them are supported with the leakage fields’ measurement.
Considering previous facts, proper design of any WPT charging system represents complex and hard task in order to be successful with the design. If we consider that currently manufacturers of power electronic components offer circuit simulators enabling to provide highly accurate circuit simulations, the point is that there is still lack of simulation tools, which could simplify design procedure of WPT system. Power converters are currently available in the form of complete modules practically suitable for wide application spectrum. The coupling coil system therefore represents very important system component, which requires increased attention considering optimized design. For this purpose, an automated in-loop design procedure is beneficial, while it considers most important issues of the coil system design together with the possibility of optimization of system parameters like efficiency, EM performance, efficiency and so forth. Given approach can be realized using simulation tool. Therefore, this article describes the design and development of simulation model of the WPT system, verified by experimental measurements. As was mentioned, considering automotive industry, the series-series compensation topology (
Figure 1) is proper choice for charging, therefore design process considerers given topology. Proposed model can provide precise information related to magnetic fields distribution (electric and magnetic intensities related to regulations), faster and more accurate determination of operation conditions (i.e., efficiency, peak power, transfer distance, component stresses). The advantages of proposed procedure are related to the high validity results for direct WPT system development enabling further optimization of the mentioned systems.
2. Description of Geometry Calculations of Coupling Elements of Wireless Charging System and Its FEM (Finite Element Method) Model
The simplified model of coupling coils of WChS was developed in COMSOL environment using combination of two physic interfaces:
Magnetic fields (mf)
Electrical circuit (ec)
The model is created as an 2D axisymmetric model within air environment. This approach simplifies markedly final model generation and optimizes computation time required for identification of the electromagnetic fields. For the coils themselves and shielding system, the geometry of rectangle blocks is being used (
Figure 2).
The mentioned model (
Figure 2) is fully reconfigurable in following ways:
Coil Geometry/Material Parameters
Application of Magnetic shielding/Geometry/Material Parameters
Application of Electric shielding/Geometry/Material Parameters
System parameters (Power Transfer Distance)
The model itself is divided to four sub models, which are used for wireless charging system optimization. The first two models are used to determine values of primary and secondary coil inductance and the value of parasitic resistances. The third model is used to determine mutual inductance of non-compensated WChS (without resonant capacitors added within the circuit). The fourth model is used to identify the shape of the transfer characteristics and to calculate the electromagnetic field radiation of designed system of coupling elements, while these main results are valid under consideration of series-resonant compensation of the WChS.
(a) Calculation of coil´s geometrical properties considering WChS efficiency requirements
The procedure of the calculation of basic geometrical parameters is discussed here. The first input parameter is radius of the coil´s wire
RW. Its value is in relationship to required efficiency of wireless power transfer
η, expected nominal load
RLOAD, the value of transferred power
POUT, recommended current density
J (3 A/mm
2) [
1] and the number of parallel wires
NW as follows (1):
Here it is required to define the minimal required value of the input DC voltage, while it is defined by next Equation (2):
Consequently, from the value of specific resistance of the wire
ρ (copper material), from value of operating frequency
f and from the value of magnetic permeability
μ it is possible to determine the skin depth using (3):
where
f is the operational frequency and
µ is the permeability. Using (3) and (1), the value of maximal radius of the wire (
RW) can be determined (4). Because skin depth is considered here and to avoid “skineffect,” litz wire shall be used, while the number of the wires is defined using (5):
where
RLW is radius of one wire of litz wire,
NLW is number of wires used to design litz wire. Using the value of nominal load of the system and considering efficiency of the power transfer it is possible to determine maximal admissible parasitic resistance of the coil (
Rcoil) (6)
The important parameter, which influences the value of parasitic capacitance of the coil (
Ccoil) is radial displacement of individual turns. To reach minimal value of
Ccoil (7) is being used for determination of the value of radial displacement
CRAD.
Considering specific value of the resistance of used wire and considering estimated value of parasitic resistance of the coil, it is desirable to calculate total length (
lcoil) of the coil´s wire using (8).
Using (8) together with (7) and applying requirement on the outer diameter of designed coil it is possible to determine the number of coil turns (10) and its inner diameter (11).
where
Dout is coil´s outer diameter,
Ni is
i-th number of coil´s turn,
N is the number of turns,
Din is inner diameter of the coil.
The conditions listed above (9)–(11) determines, that the calculation of the coil´s geometry is optimized to reach defined requirements on the efficiency and on the amount of the transferred power.
(b) Physical boundary and subdomain settings of the FEM model
The previous mathematical model is used for determination of the coil geometry, while the FEM model is using following physical and domain settings (
Table 1) for its individual parts.
Standard equations for modelling of magnetic fields are being used within model (12)–(14):
where “
H” is Magnetic Field Density, “
J” is Current Density, “
B” is Magnetic Flux Density, “
A” is Magnetic Vector Potential, “
E” is Electric Field Intensity, “
σ” is Electrical Conductivity, “
v” is Velocity of moving geometry and “
Je” is Externally Generated Current.
For the domain defined by the “Ampère’s Law” condition the Equations (12) and (13) are used while the Equation (14) is used in simplified form (15). For the calculation of electrical intensity Equation (16) is being used.
where “
D” is Electric Flux Density and “ω” is Angular Frequency.
The current density in the Coil’s domains is calculated based on the equation (17).
where “
N” is Number of turns and “
Icir” is Total Current.
Above mentioned equations are representing basic set of mathematical expressions, which are being used within proposed modelling approach.
3. Basic Definitions for the in-Loop Coil Design and Its Validation
It was already discussed, that developed model uses four calculation stages for identification of the operational properties of WChS. These stages can be divided within several steps:
- -
Self-inductance and parasitic resistance determination of a Non-Shielded Coil
- -
Self-inductance and parasitic resistance determination of a Coil with magnetic shielding
- -
Mutual inductance determination of non-compensated WPT system with Non-Shielded Coils
- -
Transfer Characteristic identification of compensated WPT system—without shielding
- -
Transfer Characteristic identification of compensated WPT system—with magnetic shielding
- -
Electromagnetic radiation identification of WPT system in given distance (shielded/non-shielded)
Because the FEM model is integrated within in the loop procedure utilizing script/UI in MATLAB, the definition of the relationship of geometry parameters is required (18)–(20). In this way proposed model flexibly enables to identify main coils parameters that is, self-inductances and parasitic resistances.
where “
w” is a COMSOL parameter defining width of rectangle, “
h” is a COMSOL parameter defining the height of the rectangle, “
x” is a COMSOL parameter defining the distance of the rectangle from the revolution axis, “
N” is number of coil’s turns, “
RW” is wire radius, “
CRAD” is radial distance between turns of the coil and “
Din” is internal diameter of the coil.
Second stage is based on the mutual inductance determination; thus, it is also important to define relationship within mutual position (mutual distance, misalignment etc.) of the system of coils. The coils itself are defined by Equations (18)–(20) with one difference for secondary coil, which defines the transfer distance of the WChS (21).
where “
y” is a COMSOL parameter defining the distance from horizontal axis and “
Zshift” is distance between the coils.
The value of mutual inductance is consequently determined using (22).
where “
mf.VCoil_2” is the inducted voltage on secondary the secondary coil, “
mf.ICoil_1” is current at the primary coil and “
mf.omega” is angular frequency.
The accuracy of FEM model and in loop design procedure were verified by comparing the characteristics of the relevant coupling coil system suited for practical experimentations of WChS (
Figure 3). This system is composed of two coils located within frame that allows reconfiguration of the transfer distance of the system. The transfer distance is adjustable thanks to the guide lanes in corners of the frame and the helix in its center. Upper coil (secondary side) is static and the lower coil (primary side) is adjustable. The distance range is form 5 cm to 21 cm. In this way, it is possible to identify transfer characteristic of WChS very precisely as motion step is app. 2 mm and influence from holding frame is neglected due to use of ABS material. Here it must be noticed that practical design of coupling coils was subjected to the calculation procedure defined by (1)–(11).
The first results of model verification are seen in
Table 1, where the accuracy of the model for determination of the coil’s self-inductance is higher than 98%. The measurements have been realized with the use of high-precision LCR meter BK precision 891, whose measuring range is within 20 Hz–300 kHz and max. 1% error for selected operational range. Higher relative error (using Equation (23)) is seen for the case of parasitic resistance of the coil, while for primary and secondary side it achieves approximately 7%. These results are valid for non-shielded system. The
Table 2 also consists the geometrical parameters of the coupling coils.
Because modelling procedure offers implementation of the electromagnetic shielding, similar verification listed in
Table 1 was performed for fully shielded system, that is, both primary and secondary side uses configuration of shielding shown on
Figure 4a. At this point it must be discussed, that proposed in loop design offers possibility of the shielding reconfiguration according to example shown on
Figure 3b. For the proper shielding technic, the magnetic material needs to have the highest possible magnetic permeability (µ > 800 [−]—Powder Ferrite Cores) while the electrical conductivity needs to be as small as possible (
σ < 1 [S/m]—Powder Ferrite Cores). The electric shielding on the other hand must be designed from material with high electrical conductivity (
σ > 6 × 10
7 [S/m]—Copper). The two alternatives of magnetic shielding are seen at
Figure 3a and
Figure 4b respectively. The magnetic shielding is built from ferrite cores “N92” with dimensions 25 mm by 14 mm by 2.5 mm. The frame for the shielding is 3D printed from “PLA” so there is no interference with the transfer properties of the WPT system.
Verifications of the model accuracy have been done similarly to non-shielded system, so the measurements of the inductances and parasitic resistances are evaluated (
Table 3). Once again LCR BK 891 m was used. As for previous case, the accuracy of the model for determination of the coil’s self-inductance is higher than 98%, while relative error for parasitic resistance reached approximately 7%.
The evaluation of the relative error was performed similarly to evaluation of relative error of the self-inductance. However, for this situation, the value of the error of mutual inductance was evaluated within whole range of expected transferring distance. Almost for this variable the accuracy was more than 98% within investigated transferring range (60 mm–250 mm).
For the third part of verification of the coil system from
Figure 3 the value of the mutual inductance of the system was determined (coupling coefficient). The results from experimental measurement and simulation are compared and shown in
Figure 5.
Fourth part of the model verification was focused on the investigation of the transfer characteristic of proposed WChS. Evaluated variable for this case is the value of transferred power in dependency on operating frequency, while both non-shielded system and system with magnetic shielding are analyzed. The model is supplemented by electrical circuit physics. The settings of the added physics can be seen in
Table 4 (
Figure 6).
For further evaluation of the in the loop model´s accuracy, the measurement and simulation experiment was realized considering the power transfer within the distance of 206 mm. The amplitude of input voltage was set to 325 V and the value of the load was 11.9 Ω. These parameters settings are valid for both shielded and non-shielded WPT system. As can be seen from comparisons of the measurement and simulation (
Figure 7 and
Figure 8), the transfer characteristics relevant for the simulation model reflects the shape and values to those received from experimental measurement. If non-shielded system is under consideration, the average accuracy of the proposed model is more than 95%. In the case of the shielded system the accuracy of the proposed model is higher than 96% (
Figure 8.)
The last part of the model verification is focused on the evaluation of electro-magnetic field intensity at 3 different distances from the receiving coil. Receiving coil is located at the top part of proposed WChS, while measuring coil is placed as shown on
Figure 8 right. In order to identify magnitudes of the radiated field within z-axe, the measurement of the induced voltage was performed with the use of measuring coil (
Figure 9 left). Measuring coil is made of the same litz wire as in the case of the coils for the WChS, while the diameter of the coil is 100 [mm].
Once again, the comparisons have been provided analyzing the values of the inducted voltage in dependency on operating frequency of WChS, while detection point was selected as 20 mm from the top of receiving coil. Amplitude of input voltage for this verification is set to 325 V and the load is 11.9 Ω as in previous experiments. As can be seen at
Figure 10 (Non-Shielded situation) and
Figure 11 (Shielded situation), the relative error of the results between measurement and simulation model is for both cases less than 3%.
Previous analyses showed that individual calculation stages of the designed simulation models achieved accuracy within the range of 93%–97% as shown in
Table 5.
4. User Guided in Loop Design of Wireless Power Transfer System
Based on the previous evaluations it can be said that the proposed models are accurate and can be used for optimization of the WChS and its electromagnetic shielding. All the listed Equations (1)–(22) together with previously described and verified models are used within developed User guided in loop design, which is described below.
Figure 12 shows first initialization part of developed UI. It serves for determination of all required variables for determination of consequent parameters needed for calculations. These main input parameters are as follows: output power at secondary side of WChS, expected nominal load, requirement on outer diameter of the coil, transfer distance, efficiency of coil to coil power transfer and the value of operating frequency. Based on these input definitions, the geometrical properties of the coils are being calculated together with information about minimal input voltage of the system, which is required to reach nominal output power.
The calculation of the parameters given by (1)–(11) is provided in the loop, while during each calculation step, the length of the wire is incremented for given turn of the coil (8) with defined radial distance between the turns. In the case that inner diameter of the coil decreases below 10 mm and required length of the wire is not achieved, whole calculation process is repeated while two parallel strands are consequently considered. In this way, the diameter of the wire is reduced together with radial distance between turns. The maximum number of parallel strands is set to 5. Even for 5 parallel strands the value of required length of the coil´s wire is not determined, the UI will I inform the user to increase the value of outer diameter of the coil.
Within next in the loop design step (
Figure 13) it is possible to implement the results from previous calculation steps, while if required it is possible to define those parameters themselves. For both cases it is then necessary to “rewrite parameters” (
Figure 13) to store variables for next calculation purposes. Within this design step it is also possible to define the material properties of the coil´s wire, while pre-defined setting is copper material. Similar procedure is then required for the secondary side coil.
After the step, which is relevant for geometrical design of the coils, user receives results for the non-shielded WChS. To be able to optimize material and geometrical parameters of shielding system (magnetic as well as electric) next step is given for definition of requirements on the primary coil shielding components.
Figure 14 shows that it is possible to check which type of shielding is going to be used. After selection, it is possible to define required shape, that is, fully shielded system (full plate), inner wall or outer wall what is related to inner and outer parts of the coil or finally it is possible to use partially shielded system based on the consideration of ferrite rings. If the specification of the shielding shape is done, then it is possible to define required parameters of the shielding system, that is, Magnetic Shielding Thickness, Coil Shielding Gap, External Diameter, Ring Width, Number of Rings or Gap Between Rings. In this step the range of values can be given for thickness and the distance from coil for both parts of the shielding. Together with these specifications, material properties are required like relative permeability, relative permittivity and conductivity of shielding material. Consequently, similar process is required for secondary side coil. Here the possibility to select similar configuration is available “Primary Coil Equivalent.” Next step of proposed design procedure is simulation, while during this step each of previously described FEM models are created together with their solution. Final step is postprocessing of the results, whereby considering distribution of electromagnetic field within WChS surrounding both non-shielded and shielded systems are displayed. In case of the user defined range of values, the displayed results are for the last combination of the values combination, while all other combinations are stored within working folder. The user is then able to move the simulation range in the desired direction based on the set of results.
All models created for the calculation of required parameters are stored in the application folder and can be used outside the UI. Thus, the created UI can also serve as optimization tool for design of a reconfigurable model. With the use of such approach it is subsequently possible to investigate the influence of a larger number of parameters than this interface allows.
5. Evaluation of Proposed Modeling Approach with WChS Suited for Automotive Charging
At this point the use of designed user guide will be verified through the experimentation on the physical prototype of wireless charger for automotive application operated at nominal conditions. The main input to output specifications are as follows:
- -
Input Voltage = rectified 3-phase (3 × 230 Vac) voltage, VINDC = app. 570 V
- -
Operating frequency = 85 kHz
- -
Peak output power = 9 kW, nominal power = 7 kW
Main aim of experimentation is to optimize EM field distribution for proposed WChS using presented in loop design user guide. The specifications related to the coils design have been considering above mentioned input/output parameters. Initially, the simulation analysis of the magnetic field distribution was done for the non-shielded system (
Figure 15). From results it is seen, where the effective value of magnetic induction reaches limiting field level defined by ICNIRP 2010 for public exposure [
20,
21,
22]. It is highlighted by the white dashed lines for vertical and horizontal plane.
The achieved results present enormous magnetic intensity exposure for relatively far distance if practical use and long-term charging is considered. Therefore, the application of the shielding system was defined within user guide environment (
Figure 14). For transmitting (bottom) coil the partially shielded (rings) system was used, while for receiving coil fully shielded system was utilized. The material properties and geometric arrangement of used ferrite bricks are made of N92 ferrite material with dimensions of 25 mm × 14 mm × 2.5 mm.
The results from simulation analysis of the magnetic field distribution for shielded WChS is shown on
Figure 16. Here it is seen that significant reduction of the effective value of magnetic intensity is reached within vertical direction. Compared to non-shielded system, the limiting field level was reduced from 66 cm to 25 cm. From practical point of view, it is acceptable results if other construction properties of charged system are considered (location of receiving coil within car chassis). Horizontal reduction is also evident, while limiting area was reduced by 10 cm. So here it can be said that optimization purposes of the EM radiation of WChS can be easily done with the use of proposed in loop design methodology.
Achieved simulation results have been verified by the similar experimental measurement of the EM field distribution. Graphical interpretation of the experimental measurement set-up is being shown on
Figure 17. Input voltage is supplied by DC input source (Magna-Power TS 600). It is directly connected to primary side SiC full-bridge inverter. In order to meet good operational properties and to reduce the complexity of control, the output rectifier is based on SiC bridge technology supplying programable electronic load (EA-EL 9400–100). If target application is considered, it is then expected to use secondary side converter in order to regulate voltage and current required by the application.
Experimental measurement using laboratory test stand was focused on the determination of the electromagnetic field distribution around the WChS. The leakage magnetic field distribution was analyzed within the defined plane (as for simulation), which was array in the cutting plane with 180 cm × 180 cm dimensions. This approach is based on the regulations given by international standards ICNIRP. Received results from experimental measurement are represented on
Figure 18.
From the experimental results is seen (
Figure 8) that the field distribution is like the one received from simulation analysis. If focus is given on the limiting areas, the experiment shows that effective value of 27 µT is located at around 28 cm from receiving coil (distance is measured at the top of the coil excluding magnetic shielding). Within horizontal plane, limiting threshold is achieved at the 36 cm from the side of the coil. Both results confirm very high accuracy and reliability of the proposed in loop design methodology and simulation model as well.