Design of a High Power, LCC-Compensated, Dynamic, Wireless Electric Vehicle Charging System with Improved Misalignment Tolerance

Dynamic wireless power transfer (DWPT) systems are becoming increasingly important for on-the-move electric vehicle (EV) charging solutions, to overcome range anxiety and compensate for the consumed energy while the EV is in motion. In this work, a DWPT EV charging system is proposed to be implemented on a straight road stretch such that it provides the moving EV with energy at a rate of 308 Wh/km. This rate is expected to compensate for the vehicle’s average energy consumption and allow for additional energy storage in the EV battery. The proposed charging system operates at an average power transfer efficiency that is higher than 90% and provides good lateral misalignment tolerance up to ±200 mm. Details of the proposed system’s design are presented in this paper, including EV specifications, inductive link and compensation network design and power electronic circuitry.


Introduction
The design of efficient wireless electric vehicle (EV) chargers is one of the most extensively addressed research topics in transportation electrification and intelligent transportation systems (ITS). This is because the use of wired chargers introduces a set of limitations that degrade the reliability of these environmentally-friendly vehicles and hinder their adoption by the general public. Plug-in electric vehicles (PEVs), i.e., wire-charged EVs, need to be physically connected to charging sockets to charge their batteries, which introduces high risks of electrocution due to bare wires, harsh weather conditions, etc. [1]. This has motivated research efforts towards wireless charging solutions that eliminate this physical connectivity and ease the charging process [2,3]. Nevertheless, wireless EV charging systems introduce a different set of challenges, which are mainly concerned with the impact of the charging magnetic fields on the human body for humans inside and surrounding the wirelessly-charged vehicle. The IEEE C95.1 Standard for Safety Levels with Respect to Human Exposure to Radio Frequency Electromagnetic Fields [4] and the International Commission on Non-Ionizing Radiation Protection (ICNIRP) 2010 Guidelines for Limiting Exposure to Time-Varying Electric and Magnetic Fields [5] both define a maximum permissible human exposure to a magnetic field density of 27 µT rms, assuming continuous charging fields. This safety exposure limit is acknowledged by the Society of Automotive Engineers (SAE) in their J2954 standard [6] in which the specifications of stationary wireless EV chargers are defined. The SAE J2954 standard proposes the use of resonant inductive power transfer (RIPT) to wirelessly couple the necessary power to charge the EV battery without exceeding the maximum magnetic field limit for safe human exposure. coupling factor variations. In both studies, an average mutual coupling value was utilized for LCC network tuning based on the desired coupling variation range.
In contrast to off-peak compensation network tuning, DWPT system designs in [7,[15][16][17][18][19][20] utilized closed-loop primary and secondary control to reduce fluctuations in output power at different coupling factors. In [7], power and current control loops were implemented on the primary side, offering improved efficiency and output power levels with up to 15 cm lateral misalignment. On the other hand, a secondary-side control system was proposed in [15,16]; it adjusts the secondary voltage levels to track the maximum power operating point of a S-S-compensated DWPT system at different coupling factor variations. The authors in [17] also proposed power control on the secondary side, aiming to simplify the primary-side circuitry by shifting the burden of power management to the secondary side. To leverage on the advantages of both primary and secondary side control loops and prevent primary side overload, a dual-side power control approach was adopted in [18] to avoid efficiency reduction at large coupling factor variations, also for a S-S-compensated WPT system. A similar approach was also adopted in [19] using LCC-S compensation, with the use of a low-latency communication link to guarantee real time information exchange between the primary and secondary sides. However, the utilization of an active communication link during the WPT process adds to the complexity of the system design and requires a detailed study of the different communication link parameters and the security of the information exchange process [2]. Accordingly, the authors in [20] implemented dual-side control loops independently, without the need for wireless communication between the primary and secondary sides.
Each of the aforementioned studies extensively addresses a single aspect of the DWPT system design to achieve output power stabilization, improved misalignment tolerance or maximum power transfer efficiency. Nevertheless, in order to simultaneously address all three objectives, this work integrates the different strategies addressed in those studies, including LCC-LCC network tuning and dual-side control, with a robust inductive link design, to leverage their corresponding advantages in improving the misalignment tolerance while addressing the maximum power and maximum efficiency objectives of dynamic wireless EV charging systems and acknowledging the EV energy demands. Hence, the main contributions of this paper can be summarized as follows: • Deriving a simplified relationship between the maximum power and maximum efficiency operating points of an LCC-LCC-compensated WPT system, to be used for tuning the compensation components for improved misalignment tolerance. • Providing a set of guidelines and design procedures for the design of a high power, high efficiency, misalignment tolerant EV DWC system, given a set of road and vehicle specifications, and a desired minimum received energy level to address the energy demands of a typical EV. This involves: -Tuning the LCC compensation components to ensure resonance operation, maximum power transfer and a sufficiently high power transfer efficiency at variable coupling conditions. -Improving the design of the inductive link to provide better coupling performance during misalignments. -Implementing closed-loop primary-side and secondary-side control to track the maximum power operating point in different coupling conditions.
The rest of this paper is organized as follows. Section 2 presents detailed mathematical analysis of the output power and the power transfer efficiency of the LCC-LCCcompensated DWPT system. The top-level design scenario is then presented in Section 3, followed by details of the design process and algorithm in Section 4. Details of the inductive link design and the finite element method (FEM) simulations are described in Section 5. Final circuit-level analysis and simulations are then reported and discussed in Section 6 before the paper is concluded in Section 7.

Modeling and Analysis
The block diagram of a typical inductive EV wireless charging system is shown in Figure 1, following the standardized stationary wireless EV chargers model presented in the SAE J2954 Standard [6].  This block diagram consists of several blocks that make up the wireless power transfer system. Starting from the primary side, electrical power from the mains supply is first rectified into a DC voltage to maximize its real power using an AC/DC rectifier and a power factor correction circuit (PFC). This DC signal is then input into a high-frequency inverter to be up-converted to the operating frequency of interest. Compensation networks are then required to help operate the inductive link in resonance conditions. On the secondary side, an AC-DC rectifier is utilized after the LCC compensation, in order to convert the coupled AC power to DC power that can charge the EV battery. Between the rectifier and the battery, a DC/DC converter may be used to aid in the output power control process.

Mains
In order to evaluate the grid-to-vehicle power transfer efficiency of this wireless EV charging system, the efficiency of each block in Figure 1 needs to be studied. However, the design of the primary-side rectifier and power factor correction (PFC) circuits is extensively addressed in the literature for a wide range of applications [21][22][23] with no specific EV-related details in its design. Hence, it is out of the scope of this paper. In addition, grid-related issues, including grid supply-demand analysis and load balancing, are also widely studied in the literature [24][25][26]. While they mainly impact the overall grid energy management and EV charging coordination and billing, these issues do not particularly affect the specific design of the different components of the DWPT EV charging system shown in Figure 1. Accordingly, in this work, grid power is assumed to be readily available and the study of grid-related issues is left out of the scope of this paper.
As shown in Figure 1, the inductive link consists of the primary and secondary coils, and their corresponding ferrite and shielding layers. The compensation networks are required to resonate the primary and secondary sides at the desired operating frequency, to maximize the efficiency of the power transfer from the primary side to the secondary one. A suitable compensation topology for dynamic charging systems is one in which the secondary side current is independent of the secondary voltage and only depends on the input AC voltage, i.e., constant current source operation. This is achieved by both series-series and LCC-LCC compensation networks [27]. Nevertheless, the presence of an additional series inductor and a parallel capacitor in LCC topologies makes the load current directly proportional to mutual inductance, hence preventing overshooting currents in cases of misalignment [28]. This can be observed by studying the schematic of an LCC-compensated inductive link shown in Figure 2.   The schematic in Figure 2 comprises N primary circuits connected in parallel to a single sinusoidal voltage source, representing a DC/AC inverter, and a single secondary circuit corresponding to a single EV. The maximum number of primary circuits that can be connected in parallel to a single inverter is limited by the maximum current rating and power handling capability of the inverter. The design of the inverter power ratings is part of a higher-level optimization problem in which the available grid supply, the number of vehicles requesting charge at a time and the corresponding charging energy demand are all analyzed at a specific DWC system location to determine the required charging power levels and the corresponding lengths of the charging lanes [29]. This is, however, beyond the scope of this work. Accordingly, for the remainder of this paper, it is assumed that a single inverter connected to N primary coils is used to supply power wirelessly to a single EV at a time, and the value of N is selected based on the primary coil structure, as detailed in Sections 4 and 5.
The input AC waveform, V invAC , is the fundamental harmonic of the output AC waveform of the inverter. It is connected to each primary coil, L p i , through an LCC compensation topology comprising a series compensation inductor, L cp i , with its equivalent series resistance (ESR); R cp i ; a parallel compensation capacitor, C cp i ; and a series compensation capacitor, C p i , where the subscript i = 1, 2, . . . , N is the index that maps to the corresponding primary coil. Each primary coil of inductance L p i has an ESR of R p i .
The AC-DC rectifier on the secondary side together with the battery and other potential secondary side circuits are all modeled as an AC load resistance, R LAC , connected to the secondary coil of inductance-L s -and ESR, R s , through the corresponding secondary LCC compensation network-L cs , C cs and C s -and the ESR of the compensation inductor, R cs . For better visibility, the dissipation resistors of the compensation capacitors are omitted from the schematic in Figure 2, but were included in the simulations conducted. The mutual coupling between the primary and secondary coils is denoted by M p i s and is used with the coils' respective self inductances to calculate the coupling factor, k p i s using the expression, k p i s = Without loss of generality and for the sake of simplicity, the following assumptions are made:

1.
Adjacent primary coils are placed at a sufficient separation distance, D, such that when the secondary coil is perfectly aligned on top of one of the primary coils, its mutual inductance with the previous and following coils is negligible [30]; i.e., M p i s = 0 and M p i+2 s = 0 when M p i+1 s = M ps max for i = 1, 2, . . . , N where M ps max is the mutual inductance at perfect alignment. This assumption is validated by the FEA simulation results reported later in this paper.

2.
The distance D is also large enough such that the mutual inductance between the two adjacent primary coils is very small and can be neglected. 3.
All the primary coils are identical with same geometries, equal self-inductances and equal ESRs.

4.
Each set of simultaneously active primary coils, hereafter referred to as a section of N primary coils, is connected in parallel to the same sinusoidal inverter voltage, V invAC . 5.
AC power from the mains grid is readily available for the dynamic EV wireless charging system. 6.
The EV is traveling at a constant speed, U, along a flat 0% grade road. This is an expected regulation on DWC lanes that helps maximize the net energy received by the traveling EVs by avoiding an increased energy consumption during the acceleration/deceleration of EVs on the charging lanes. 7.
The initial EV battery SoC, SoC i , for an EV demanding dynamic wireless charging falls between 20% and ∼ 90%. This is required to ensure that the battery is charged in the constant current (CC) charging mode for a typical Li-ion battery [31].
It should be noted that EVs with SoC i > 90% are more effectively charged with a constant voltage (CV) rather than a constant current. Nevertheless, since dynamic wireless charging (DWC) systems are mainly intended to address range anxiety of the EV drivers, EVs with SoC i > 90% are not expected to require DWC and are hence excluded from this study. Accordingly, only CC charging mode is assumed in this work.
Hence, acknowledging the aforementioned assumptions, the circuit in Figure 2 is analyzed as follows. By Faraday's law of induction, the electromotive force (EMF) induced across each inductor coil can be expressed in terms of the flux linkage λ as ε = dλ dt , where the flux linkage between mutually-coupled inductors is defined as λ p i = L p i I p i − M p i s I s , and λ s = M p i s i p i − L s i s . Accordingly, Kirchoff's voltage law (KVL) equations are derived for the schematic in Figure 2 in the frequency domain, while neglecting the ESRs of the coils for simplicity. The equations are as follows.
Hence, in order to maximize the power transferred from the primary side to the secondary side of the inductive link, the system needs to operate at a resonance frequency, f 0 , such that minimum power is lost during the power transfer [32,33]. To achieve resonance, the values of L cp i , L p i , C cp i and C p i need to be designed such that the input impedance seen at the primary side is purely real and its reactive component is equal to zero-i.e., the input impedance has zero phase angle (ZPA) at f 0 and the corresponding angular resonance frequency, ω 0 = 2π f 0 [34]. The expression for the input impedance is derived in [34,35] for an inductive link with a single primary coil as follows: For the multi-primary inductive link setup assumed in this work, as shown in Figure 2, and since all the primary circuits are assumed to be identical with identical component ratings, the mutual inductance variable M p i s is assumed to incorporate the joint coupling between the primary and secondary sides. This includes coupling during the perfect alignment of the secondary coil with one primary coil and its partial alignment between two primary coils. Accordingly, the variable M p i s is the mutual inductance profile of the overall system, incorporating both lateral misalignments of the secondary coil and its longitudinal displacement from one primary coil to another. This is reflected in the simulations described in the following sections of this paper.
Hence, by setting the imaginary component of Z in to zero, the values of the compensation components required to achieve resonance at the angular resonance frequency ω 0 need to satisfy the following expressions [19,34,36]: 1 1 It is observed in Equations (6)-(8) that in order to achieve resonance, the values of the compensation components are independent of M p i s and R LAC , which is required to ensure that the system remains in resonance despite variations in the coupling conditions during misalignment [27]. In addition, by substituting the resonance Equations (6)- (8) in (1)-(4), it is observed that the current through the primary coil, I p i , is independent of M p i s at resonance conditions, and only depends on the input AC voltage and the value of the compensation inductor, L cpi . This means that for a constant V invAC , the current through the primary coil is constant at a given value of the compensation inductor, which is important to ensure that the coil is operating at its rated current condition [28]. The current through the load, I cs , is also constant for a constant V invAC , and is independent of R LAC while only depending on the values of the compensation inductors. Nevertheless, it should be noted that the ZPA condition would be slightly relaxed during the design by phase, delaying the inverter current and making the circuit slightly inductive, to achieve the zero-voltage switching (ZVS) condition. This is required to minimize the switching losses of the MOSFET devices in the inverter [28,37].
The average power transferred to the load can be calculated using the expression P R LAC = V LAC I cs , which can be written in terms of inverter AC voltage, V invAC , as: As observed in (9), the output power is directly proportional to the mutual inductance M p i s , which highly improves the light load and no load conditions by decreasing the output power as the mutual inductance decreases. In addition, the dependence of P R LAC on L cp i and L cs makes the selection of these component values dependent on the maximum output power level required by the DWPT system. Accordingly, based on the aforementioned analysis, the compensation components need to be designed such that the ZPA condition occurs only at a single value of ω 0 , for all variations in the coupling conditions, to maximize the output power while eliminating the current harmonics otherwise associated with frequency bifurcation. Frequency bifurcation is the condition in which the frequency to achieve ZPA is not unique [38], and it typically occurs with very small load resistance values [39]. In S-S-compensated power transfer systems, a limitation is imposed by the authors in [40][41][42][43] on the maximum allowed coupling factor to avoid frequency bifurcation, which is referred to as the critical coupling factor. However, for LCC-LCC-compensated networks, the additional compensation inductors, L cp i and L cs , introduce an additional degree of freedom that allows the system to be designed to eliminate frequency bifurcation without limiting the maximum required coupling factor. Accordingly, the values of L cp i and L cs are calculated based on k ZPA ; the minimum coupling factor beyond which the ZPA condition is guaranteed to occur. Hence, by assuming that L p i = L s and L cp i = L cs , the values of the compensation inductors can be calculated using: In addition, since the LCC-compensated topology in this work is desired to operate as a constant current source, the values of the compensation inductors are selected based on the desired constant output load current, I cs , by replacing the fraction V LAC/P LACmax by I cs . Hence, the expression to calculate the value of the compensation inductor L c for constant current source design is given by: The selection of k ZPA for tuning the compensation inductors is addressed in [9,10,44] using coupling factor averaging techniques. However, in this work, the authors propose to design the compensation inductors using the coupling factor for the desired misalignment tolerance level, i.e., by setting k ZPA = k min . Accordingly, this guarantees that the ZPA condition occurs only at ω 0 for all k > k min , and hence the transferred power is maximized over this coupling range.
In addition to maximizing the power received at the load, the power transfer efficiency of the LCC-LCC-compensated system needs to be evaluated. This is done by acknowledging the effect of the coil ESRs, and assuming the secondary coil couples with one primary coil at a time. The coil quality factors are defined as Q p i = ω 0 L p i R p i and Q s = ω 0 L s R s , representing the ratios of the coils' self inductance to their equivalent ESRs at the system operating frequency. Accordingly, the power transfer efficiency can be expressed in terms of Q p i and Q s as follows: From (12), the expression kQ = k 2 p i s Q p i Q s can be used as a figure-of-merit (FoM) for the design of the primary and secondary coils. The higher the FoM, the higher the coupling and quality factors of the coils and the lower the losses in their ESRs.
Conventionally, the power transfer efficiency can be maximized by differentiating it with respect to R LAC to determine the optimum value of the load and the corresponding maximum achievable efficiency. By performing this differentiation operation, the expression of the optimum load, R LAC opt , is obtained as follows: However, the value of R LAC is one of the design constraints specified by the desired maximum output power and constant output current levels, using the expression R LAC des = 8 π 2 P max/I 2 DC , where I DC is the DC battery current during CC charging mode, P max is the maximum output power and 8 π 2 is the DC-AC impedance conversion factor assuming ideal diode rectifiers. Under this assumption, I DC can be related to I cs using I cs = pi 2 √ 2 I DC . Accordingly, Equation (13) can be re-arranged to obtain the optimal value of the coupling factor k that maximizes the power transfer efficiency in terms of R LAC des , by assuming that k 2 p i s Q p i Q s 1 and Q p i = Q s = Q. The optimal value of k is then obtained as follows: Acknowledging the aforementioned assumptions and by substituting L cs = L c| Pmax from (10) into (14), it can be observed that the condition for maximum power transfer, i.e., ZPA at constant current source operation, and that for maximum efficiency, can be related using Hence, using (11) and (15), the effect of the compensation inductor design on the maximum transferred power and the AC-AC efficiency at maximum lateral misalignment is highlighted. The relationship established in (15) identifies the maximum power and maximum efficiency operating conditions based on the input-output voltage ratio. This is utilized to simultaneously address the misalignment tolerance, maximum power and maximum efficiency objectives of dynamic EV charging systems. In this work, k ZPA is selected to maximize the transferred power at the required misalignment tolerance, while setting a lower bound on the power transfer efficiency throughout the power transfer operation and maintaining high received power at the perfect alignment condition. This is further detailed in the design process in Section 4.

Dynamic EV Charging Scenario
While the sole objective of a static wireless charging system is to increase the SoC of the EV battery, dynamic wireless charging systems aim to compensate for the energy consumed by the EVs while on the move to prevent excessive depletion of the EV batteries and hence prolong the battery's lifetime [45]. Accordingly, to maximize the gains from dynamic charging, charging lanes need to be placed in urban areas with medium-range driving speeds such that the dependence on the EV battery is minimized for trips within the city boundaries where charging lanes are implemented. Ultimately, this translates into ensuring that the SoC at the end of the charging lane exceeds that at the beginning of the lane. In CC charging mode, and based on the design of the charging controller [46], the SoC levels can be linearly related, within the assumed SoC range from 20% to 90%, to the energy consumed by the EV and the energy received from dynamic charging using the expression: where E Rx is the energy received from the DWPT system, E con is the amount of energy consumed by the EV as it traverses along the road and E max is the maximum EV battery capacity. Let T s be the time taken by the EV to traverse along a charging lane of length X chg km. The energy received by the vehicle from the WPT system, at a constant velocity U, is calculated using: where P LAC is the received power on the secondary side, obtained using (9), and x is the direction of the forward motion of the EV. Accordingly, in order to incorporate the grid power supply for EV charging as a parameter in the design process, a lower bound on the received energy is required, at a minimum acceptable efficiency, to deem the design feasible and effective. The received electric energy E Rx is not fully converted to mechanical energy due to the electrical losses of the EV motor. According to [47,48], these losses can be estimated to be around 10%-20%, which means that the received energy from the DWPT system needs to be at least 1.25 times the energy consumed by the EV. In practice, however, it is desired to further increase the battery's SoC to allow the EV to travel over non-charging lanes without depleting its battery. Hence, within a city's infrastructure, it is suggested to implement charging lanes at equal separation distances around the city, such that the average distance traveled by the vehicle without charging is approximately equal to the total length of the deployed charging lanes. This requires receiving at least twice the vehicle's consumed energy over each charging lane, while accounting for the efficiency of the power transfer process. Hence, the energy design objective in this work is defined as follows: This lower bound on the received energy addresses the ultimate objective of minimizing the utilization of the EV battery during urban city drives, and is used to determine the minimum required AC supply power to achieve this objective for any given set of lane and EV specifications, and inductive link parameters.

Design Methodology
Based on the aforementioned analysis, the design of the complete LCC-LCC-compensated wireless charging system involves the following key stages:

1.
Definition of the charging lane and EV specifications.

2.
Setting the energy design objective as per (18) using estimated E con .

3.
Design and FEM simulations of the inductive link to maximize FoM and improve misalignment tolerance.

4.
Circuit-level design and simulations of the DWPT system to achieve the power, efficiency and misalignment tolerance objectives.
Following the analysis presented in Section 2, an optimized DWPT system design maximizes the power transferred at different misalignment conditions while maintaining a satisfactory power transfer efficiency. For this work, the minimum allowed AC-AC power transfer efficiency was set to 90%. The choice of k ZPA for maximum output power was then initially set to the maximum allowed misalignment tolerance, and the circuit was simulated to verify that the received energy meets the objective in (18) while maintaining η AC−AC ≥ 90% over the entire coupling range for k > k ZPA . In case the received energy does not satisfy the energy design objective while high efficiency is maintained, k ZPA can be further reduced to allow higher power transfer at low coupling conditions, at the expense of a decrease in efficiency. Otherwise, if the efficiency condition is not met with the selected k ZPA , k ZPA needs to be increased together with the maximum power limit P max to provide the required energy to compensate for the EV consumption while achieving the minimum efficiency requirements. This is further clarified in the flow chart shown in Figure 3.

Charging Lane Design
The deployment of dynamic charging lanes in a city infrastructure depends on the city's arterial network, traffic flow and traffic density at different times throughout a day, amongst several other factors [49]. Nevertheless, to demonstrate the design process proposed in this work, a simplified single lane charging scenario is used for the analysis. This consists of a single charging lane of length X chg km, with an average free-flow vehicle traveling speed of U km/hr. The charging lane is divided into sections, each of which is powered from a dedicated inverter circuit. The design of each section depends on the structure of the primary coils, which can be a single elongated track per section or a segmented set of primary charging pads [3,50,51]. These are demonstrated in Figure 4, where a single road section of length S km can have N primary coils, each of length L and spaced by a distance D, or a single track primary coil of length NL + (N − 1)D. In this work, the spacing between the adjacent sections is set to L + D to allow for the on/off switching between the sections.  As shown in Figure 4, the track primary coil is an elongated structure implemented such that a constant magnetic field is presented to the electric vehicle as it traverses along its length [50,52]. The length of this primary track is several times larger than the secondary coil embedded beneath the chassis of the EV, which is intended to provide a constant mutual inductance profile and enhance longitudinal misalignment tolerance. However, this structure causes significant power losses and safety concerns, due to stray fields occurring when only a single vehicle is traversing over the long track at a time. For a segmented primary structure, on the other hand, the size of each segment is comparable to that of the secondary coil embedded in the vehicle's body, which means that the length of a single section is almost equal to the length of a typical EV. This is expected to improve the coupling performance as it focuses the coupling magnetic fields within the area enclosed between the primary and secondary coils of comparable sizes. This accordingly reduces stray fields, minimizes field leakage to the surroundings and improves the power transfer efficiency [3,51,53,54].

Vehicle Specifications
The EV model selected for this work was the Nissan Leaf 2019 model [55], whose dimensions, weight and other specifications are presented in Table 1 as reported in [55,56] based on the recommended vehicle operation conditions. The EV energy consumption for the EV model in this work was estimated from U.S. Department of Energy report for baseline on-road testing of the Nissan Leaf at 70 km/h, zero acceleration and 0% road grade [57]. This was used to calculate the design objective using (18), and the corresponding values are reported in Table 2.

Inductive Link Design and FEM Simulations
The first step in designing the inductive link is to determine the resonant operating frequency, f 0 , at which the inductive link shall resonate. The SAE J2954 standard [6] identifies 85 kHz as the nominal operating frequency for static wireless electric vehicle chargers. Since a single secondary coil is expected to be fitted at the bottom of the EV, it is required to operate for both static and dynamic charging modes. Hence, 85 kHz is proposed to be used as the operating frequency for dynamic wireless charging systems and is used for the design presented in this work. At f 0 = 85 kHz, a heuristic design approach was followed using FEM simulations to optimize the dimensions, geometry and construction of the inductive link to maximize the kQ FoM over a range of lateral misalignments. That is detailed as follows.

Ferrite Layer
The ferrite layer is a layer of ferrite material placed below the primary coil and on top of the secondary one, with a higher magnetic permeability than the copper used for coil construction. The main role of the ferrite layer is to intensify the magnetic fields between the two coils and hence increase the coupling factor, k, by increasing their mutual inductance. A ferrite material with high magnetic permeability, µ r , is required to minimize the dissipated energy and reduce core losses. In addition, the ferrite link geometry needs to be selected to provide enhanced coupling performance at different lateral misalignments. The ferrite layer specifications in this work follow the authors' earlier work in [58,59] in which I-core ferrite bars with µ r = 3000 provide satisfactory coupling performance over a wide range of lateral misalignments. These were used for all FEM simulations reported hereafter.

Shielding Layer
In addition to the field enhancement provided by the ferrite layer, effective electromagnetic shielding is required to ensure that the magnetic field density surrounding the inductive link does not exceed the maximum allowed limit of 27 µT, as dictated by IEEE C95.1 standard [4] and ICNIRP guidelines [5]. In this work, an aluminum sheet was used as a shielding layer-specifically, aluminum 1050A alloy that has a conductivity of σ = 33.9 MS/m [41]. In order to ensure effective shielding, the thickness of the aluminum shield needs to be at least twice the skin depth of the aluminum alloy at the link operating frequency, ω 0 [41]. This is because, after each skin depth of thickness, the field intensity is reduced by a factor of 1 e , and hence twice the skin depth shall result in a more than 70% reduction in field intensity. The skin depth of aluminum alloy is calculated using δ = 2 ω 0 µ 0 σ to be approximately equal to 0.3 mm. Nevertheless, to ensure effective shielding over a wide operating frequency range, the aluminum shield thickness is set to be six times the skin depth, i.e., 1.8 mm, based on earlier work by the authors in [60]. In addition, to reduce eddy current losses in the aluminum shield, a vertical gap of 5 mm is maintained between the ferrite layer and the aluminum shield.

Coil Design
Using I-core ferrite bars and a 1.8 mm thick aluminum shield, the coil structure was designed and simulated to maximize the kQ FoM over the widest range of lateral misalignments. Based on earlier studies [2,3], rectangular coils are proven to provide better magnetic flux distribution compared to circular coils and are less complex to design than double-D (DD) coils. In addition, planar coil structures provide better space efficiency in comparison with multi-layer designs. Hence, planar rectangular primary and secondary coils were selected for this work.
The size of the inductive link affects the self and mutual inductances of the coils, and according to [61], it needs to be proportional to the vertical separation distance, i.e., the air gap, between the coils. The area of the secondary coil in an EV charging system is restricted to the available space in the bottom chassis of the vehicle, and the vehicle-to-ground clearance is expected to be between 15 cm and 30 cm for most EVs. According to [62], a reasonable secondary coil occupancy is around 480,000 mm 2 for a typical sedan EV. The outer charging pad dimensions, on the other hand, need to be larger to ensure effective shielding. This was selected to be 1000 × 800 mm 2 . In addition, according to [63], the optimum value of the width of the coil is three times the height of air gap. Since the average air gap was estimated to be ∼ 200 mm, the outer coil width was set to 600 mm. Hence, for an area of 480, 000 mm 2 , the outer length of the coil was set to 800 mm.
Since the coils selected for the design of EV charging pads in this work were planar coils, the numbers of turns of the coils contributed to their self and mutual inductances, their ESRs and their coupling performance. A first-order approximation of Neumann's inductance calculation model is presented in [63], based on the geometric parameters of two air-core rectangular coils with no ferrite and aluminum shielding. Nevertheless, in order to account for the planar nature of the desired rectangular coils and the effects of the turn-to-turn spacing, ferrite and aluminum shielding layers, a set of FEM simulations was conducted on Maxwell simulator for different numbers of turns while setting N p = N s . The obtained results were fit onto the first order approximation model to obtain a corrected contour plot of coupling factor values at different primary and secondary turns. This is shown in Figure 5.
Accordingly, it was observed that for equal primary and secondary windings and dimensions, the rate of increase in coupling factor decreases as number of turns increases. This is shown in Figure 6.
Nevertheless, increasing the number of turns increases the amount of conductor wires used in coil construction and increases the ESR of the coil. The number of turns for both primary and secondary coils was then selected to be 16 based on earlier work by the authors in [60], as shown in Figure 6, and the coils' outer dimensions were set to 800 × 600 mm 2 . The proposed inductive link design is summarized in Table 3.    Table 3, the coupling profiles of the coils at different lateral misalignments were obtained as shown in Figure 7. A single road section was then constructed using the proposed inductive link design, following the structure in Figure 4b with three primary coil segments with coil lengths of 800 mm and a coil-to-coil spacing of 100 mm while sharing a common aluminum shielding pad of length 2800 mm. The corresponding track primary coil with the same outer area and same number of turns was also constructed and simulated. These are shown in Figure 8 and the corresponding coupling profiles are shown in Figure 9. As shown in Figure 9, the segmented primary configuration meets the two design conditions defined in the assumptions listed at the beginning of this paper: k p i s ≈ 0 and k p i+2 s ≈ 0 when k p i+1 s = k p i s max and k p i p j < for i = j. In addition, it was observed that the coupling factor obtained with an elongated track primary coil was relatively constant, yet significantly smaller than that obtained using a segmented primary coil. This, as highlighted in Section 3, was due to the smaller flux linkage between the secondary coil and the longer track coil in comparison to equal-sized primary coil segments. This reduced flux linkage is due to the fact that the entire track provides a single flux path that only partially couples with the smaller secondary coil. Accordingly, mutual inductance is much smaller in comparison with the large self inductance of the track primary coil, and the coupling factor is low, as shown in Figure 9. Since the transmitted AC power is directly proportional to mutual inductance as per (9) and the received EV energy is the integral of the AC power received along the length of the track as per (17), the ratio of the received energy using segmented vs. track primary structures can be related to the area under their coupling profiles. Hence, by observing the area under the two coupling profiles in Figure 9, it can be concluded that the segmented primary coils shall provide higher received energy in comparison with track coils. Hence, segmented primary coil structures were selected for the DWPT system design presented in this work. Coupling factor, k k p1s -Segment k p2s -Segment k p3s -Segment k ps -Track Figure 9. Coupling profile of the proposed 3-primary road segment: k p 1 s Segment , k p 2 s Segment and k p 3 s Segment , vs. a single track primary segment, k ps Track .

Circuit-Level Design and Simulation
Once the inductive link design was finalized, FEM simulation results were exported into a circuit-level simulation on MATLAB/SIMULINK platform to simulate the operation of the compensation networks, the inverter and rectifier circuits and the EV battery together with the inductive link model. The inverter, rectifier circuits and control loops were designed first, and the EV battery was modeled; then the complete system was integrated and simulated accordingly.
The implementation of feedback control loops in the design of DWPT power electronic circuits is required to compensate for the variation of mutual coupling between the primary and secondary coils during the vehicle's motion, while providing the necessary limitations on the input/output currents to protect the circuit components. While LCC-compensated circuits are inherently tolerant to misalignments up to a certain extent, as explained in Section 2, feedback control loops are required, in addition to the compensation components' design, to allow tracking of the maximum power operating point of the system by maintaining ZPA operation at different coupling and loading conditions during the vehicle's motion. In this work, primary side control is implemented using inverter power and current control loops, inspired by the designs in [7,8,[18][19][20], to achieve maximum power point tracking. On the secondary side, the controller is utilized to ensure that the output DC current does not exceed the constant charging current desired to charge the battery.

Inverter Design and Control
In the analysis presented so far, the input voltage was modeled as a sinusoidal waveform denoted by V invAC , which represents the fundamental harmonic of the square wave output voltage of a class D DC-AC inverter. The duty cycle of the inverter square wave output voltage needs to be controlled to track the maximum power point at different coupling conditions. This is achieved by controlling the phase shift, ϕ p , between the two inverter arms, which, in turn, controls the duty cycle θ and the magnitude, v m , of the fundamental sinusoidal signal, V invAC , and controls the output power of the inverter. The relationship between v m and θ is shown in the following expression [7,64,65]: where V inDC is the DC voltage input to the inverter, and the phase shift ϕ p and the duty cycle θ are related as θ = 180 − ϕ p [66]. The feedback loops implemented to perform this primary-side phase shift control are shown in Figure 10. According to Figure 10, the inverter output power is compared to the desired system power level, P max , to track the maximum power operating point, and the error signal, m p , is used to control ϕ p of the switching control signal S(t), and hence control θ and the magnitude of the inverter voltage. Accordingly, a pulse-width-modulated (PWM) inverter is operated whose output power is controlled based on the desired maximum output power. In addition, the amplitude of the inverter current, I inv , is also compared with the maximum allowed current I max , based on the maximum current handling capability of the inverter MOSFETs, to generate the error signal m I . m I modulates the amplitude of S l (t) and S r (t), the left and right switching control signals, with switching frequency f s , while the power error signal, m p , modulates the phase of S r (t) only to introduce the desired phase shift between the inverter arms. S l (t) and S r (t) are then compared with a triangular waveform at the same f s to generate the desired switching signals to the gates of the MOSFETs to provide the required power and current control.

Rectifier Design and Control
The conversion of the received AC waveform into a DC signal takes place through a rectifier-regulator pair, to provide a regulated DC output that can charge the EV battery. The input to the AC-DC conversion stage is a sinusoidal current, and a boost rectifier circuit is selected for this work to ensure operation in continuous current mode (CCM) [22]. In particular, a symmetrical bridgeless boost rectifier topology was selected, as it provides reduced control complexity and a lesser number of active semiconductor devices per conduction path [67][68][69]. Nevertheless, the conventional model utilizes an inductor in the AC signal path, which adds to the value of the series compensation inductor, L cs , and disturbs the resonance condition. In this work, the model was slightly modified by shifting the inductor to the DC side of the boost rectifier circuit, in order to perform the desired regulation function without disturbing the resonance of the system. This modification is demonstrated in Figure 11. As power control is implemented on the primary side of the proposed DWPT system, the key requirement for secondary side control is to ensure constant current is delivered to the EV battery to maintain CC charging mode with minimum current variations. This is implemented on the secondary side, and the control loop is implemented as shown in Figure 12. As shown in Figure 12, the output current of the rectifier is compared with the desired constant DC battery current I DC , and the current error signal, m Irec , is used to modulate the switching control signal of the rectifier MOSFETs, S rec (t). The switching control signal is set to provide a phase shift of ϕ r = π/2 to the modified symmetrical bridgeless boost rectifier circuit, to maintain the phase characteristics and provide effective rectification of the boost rectifier, without disturbing the resonance condition.

EV Battery Modeling
Li-ion batteries are utilized particularly often for EV energy storage systems, as they offer high energy densities, good life cycle performance, low self-discharge and long service lifetimes [70,71]. While electrochemical and mathematical battery models exist in the literature [72,73], electrical models are the most suitable for circuit-level simulations and analysis, as they use an equivalent circuit model of the battery and incorporate its transient charging/discharging behavior [74].
Accordingly, a first-order Thevenin equivalent electrical model of a Li-ion cell was used to evaluate the performance of the proposed DWPT system on EV battery charging. The model consists of a constant open-circuit voltage, V o , connected to an ohmic resistor, R o , in series with a parallel RC network consisting of R 1 and C 1 that models the transient response of the battery. For the Nissan Leaf 2019 EV model selected for this work, the 40 kWh battery is composed of two parallel modules, each consisting of 96 Li-ion cells with nominal voltage of 3.8 V [75,76]. Accordingly, the EV battery was modeled on SIMULINK using the circuit shown in Figure 13, where I bat is the charging current in each battery module and I DC = 2 * I bat for the two n-cell battery modules in the Nissan Leaf 2019 battery pack. Based on the experimental analysis conducted in [75,77], the internal resistances of a Li-ion battery remain relatively constant within the SoC range assumed in Section 2, i.e., for 20% ≤ SoC ≤ 90%. Hence, the values of R o and R 1 in Figure 13 are constants over the DWPT system operation range, and the value of C 1 is also constant for a particular battery time constant, τ = R 1 C 1 , which can be determined by experimental battery testing. The open circuit equivalent voltage V o , on the other hand, is dependent on the battery SoC, and hence the values are obtained from look-up tables and/or graphs based on tests conducted for each battery model at different charging/discharging currents [31,78,79].

SIMULINK Test Bench
With the primary and secondary control loops, the EV battery model and the inductive link design parameters, the complete DWPT system test bench was built on SIMULINK platform following the schematic in Figure 14. The desired P max was set to 40 kW with input DC voltage of 500 V. A 40 kW power supply for EV charging lanes was considered reasonable, as it corresponds to an on-demand charging process in which the required power is drawn only as the vehicle passes over the charging coils. In fact, ongoing research and development efforts into DWPT EV charging systems have considered a power supply range from 20 to 100 kW based on the size of the vehicle, the air gap between transmitter and receiver coils, the power transfer efficiency and other factors [7,8,80]. Accordingly, the DWPT system ratings are summarized in Table 4. Accordingly, based on the self inductances of the primary and secondary coils obtained from FEM simulations and the corresponding ESRs and coupling profiles, the values of the compensation components were calculated using an initial value for k ZPA and the desired P max , as shown in the DWPT design algorithm in Figure 3. The value of k ZPA was initially selected to be the coupling factor at a lateral misalignment of ±200 mm, i.e., kZPA ≈ 0.15, according to Figure 7. A summary of the components values is shown in Table 5. The frequency bifurcation phenomenon is first studied by plotting the phase of the input impedance evaluated using (5) for the compensation components' values calculated at k ZPA = 0.15. This is shown in Figure 15. By observing the phase of the input impedance in Figure 15, it is noted that the ZPA condition occurs only at ω 0 for all values of k > k ZPA , and hence the frequency bifurcation effect is eliminated for the DWPT system designed for a lateral misalignment tolerance of ±200 mm. The next step is then to determine the amount of energy received by this system and its AC-AC power transfer efficiency to satisfy the remaining design objectives.
To simulate a realistic system operation, the ESRs of the compensation components also need to be acknowledged. The corresponding ESRs of the compensation inductors are estimated by assuming their quality factors to be equal to 750. A dissipation factor of DF = 0.1% is also assumed for the capacitors, from which their ESRs are calculated as R c = DF ω 0 C . The parameters of all circuit components are then used in the SIMULINK test bench, and the output power at different secondary coil positions is plotted and mapped to the corresponding coupling factors, as shown in Figure 16. By observing the power profile in Figure 16, it is noted that the maximum power is transferred to the EV battery at k ZPA rather than at the maximum coupling point. This is essential to improving the misalignment tolerance of the DWPT by allowing higher output power levels at reasonably lower values of the coupling factor. The received power at k max is slightly lower, with a variation of only 3%. This variation is acceptable as long as the received energy satisfies the EV requirements in Table 2 and the average overall efficiency is maintained above 90%. The expected received energy in the EV battery is hence calculated by integrating the power profile in Figure 16, and the result is summarized in Table 6. The AC-AC power transfer efficiency of the proposed system design is shown in Figure 17.
Since the expected received energy per km exceeds the design objective in Table 2, and the average efficiency exceeds 90% for all regions of k > k ZPA , the proposed system satisfies the requirements of a misalignment tolerant, dynamic, wireless EV charging system that both compensates for the EV energy consumption and increases the stored energy in the EV battery. In fact, the difference between the received energy per km from Table 6 and the EV energy consumption per km in Table 2 is the surplus energy that is stored in the EV battery to enable EV motion over road sections with non-charging lanes. Hence, by substituting these energy figures in the SoC calculation equation in (16), the battery SoC can be estimated to increase at a rate of 0.5% per km of on-road wireless charging cells in addition to the energy utilized for the motion. Nevertheless, Equation (16) is a simplified SoC estimation method and does not account for the battery temperature, aging and other factors that affect the battery energy levels. More sophisticated SoC estimation methods are reported in the literature [81] and shall be used in future work to further enhance the battery model. In addition, future enhancements of the proposed system model shall also include the coordination of the energy exchange between the EV motor and the EV battery, and the potential utilization of supercapacitors [82,83] to ensure effective utilization of the energy received from the DWPT system. Nevertheless, in order to assess the performance of the proposed system, the results presented in this section are compared with the reported results for other DWPT system designs in the literature [9,10], as shown in Table 7. Table 7. Comparison between proposed system and other DWPT charging systems.

Reference
Ref. [9] Ref. [ As shown in Table 7, while the presented system selects k ZPA for the maximum power point to be equal to the coupling factor at maximum allowed misalignment, the studies presented in [9,10] both used an averaging equation to obtain their corresponding coupling factors to be used for tuning the compensation components. The objective of the authors in [9] was to minimize fluctuations in the voltage gain over the desired misalignment range, whereas the authors in [10] aimed to maximize the power transfer efficiency over that range. On the other hand, the work presented in this paper simultaneously addresses compensation network tuning, inductive link design and power electronic circuit analysis to achieve high misalignment tolerance, maximum power and maximum power transfer efficiency. The performance of the proposed system is hence superior to the other two systems in terms of the misalignment tolerance range, the maximum allowed power variations and the power transfer efficiencies at the maximum allowed misalignment point, while exhibiting comparable efficiencies at perfect alignment conditions. It is essential to highlight that the reported results from [9,10] in Table 7 are mainly the results of the simulations conducted by the authors in both works, which thereby present a fair comparison to the results presented in this work.

Conclusions
This work demonstrates the design process of an LCC-LCC-compensated DWPT system for charging electric vehicles while maximizing the received power at different misalignment conditions and maintaining a sufficiently high power transfer efficiency. By leveraging the inherent misalignment tolerance capabilities of LCC-LCC-compensated systems, effective control loops were designed to ensure maximum power point tracking with minimal degradation in efficiency. While the proposed design utilizes a particular inductive link structure, the detailed design process presented in this work shall be effectively utilized for different inductive link designs given they have sufficiently high kQ figures-of-merit. Accordingly, it is concluded that the optimal design of the DWPT system requires the design of the compensation components to satisfy the ZPA design objectives, while maximizing the kQ FoM of the inductive link and tuning the respective power electronic circuitry to optimize the joint operation of the different system blocks. This paper serves as a guideline for DWPT system developers aiming to design functional dynamic EV wireless charging systems based on sets of design specifications and objectives. Funding: This work is jointly supported by the American University of Sharjah through SCRI grant number SCRI 18-CEN-10, and by Sharjah Research Academy (SRA), Sharjah, United Arab Emirates.

Conflicts of Interest:
The authors declare no conflict of interest.