Multiphysics Investigation of an UltrathinVehicular Wireless Power Transfer Module for Electric Vehicles

: The functional and spatial integration of a wireless power transfer system (WPTS) into electric vehicles is a challenging task, due to complex multiphysical interactions and strict constraints such as installation space limitations or shielding requirements. This paper presents an electromagnetic– thermal investigation of a novel design approach for an ultrathin onboard receiver unit for a WPTS, comprising the spatial and functional integration of the receiver coil, ferromagnetic sheet and metal mesh wire into a vehicular underbody cover. To supplement the complex design process, two-way coupled electromagnetic–thermal simulation models were developed. This included the systematic and consecutive modelling, as well as experimental validation of the temperature- and frequency-dependent material properties at the component, module and system level. The proposed integral design combined with external power electronics resulted in a module height of only 15mm. The module achieved a power of up to 7.2 kW at a transmission frequency of f 0 = 85kHz with a maximum efficiency of 92% over a transmission distance of 110mm to 160mm. The proposed simulations showed very good consistency with the experimental validation on all levels. Thus, the performed studies provide a significant contribution to coupled electromagnetic and thermal design wireless power transfer systems.


Introduction
Wireless power transfer systems (WPTSs) are considered as a promising technology to increase the acceptance of electric vehicles (EVs) [1,2]. Compared to conductive charging, WPTSs offer vehicle charging without user intervention. This means that no cable handling and plugging are required. Thus, WPTSs provide a low-maintenance charging technology for EVs with high comfort and safety.
The working principle of WPTSs for electric vehicles is based on the electromagnetic interaction of an offboard ground pad module (GPM) and an onboard car pad module (CPM). The GPM typically comprises four main elements: a rectifier circuit, an inverter, a matching network and a transmitter coil. The CPM consists of a receiver coil, a matching network, a rectifier and, if necessary, a DC/DC converter. The architecture of an exemplary WPTS is illustrated in Figure 1.
This paper focuses on the onboard CPMs, which are usually mounted on the underside of EVs [3]. This conceptional position results in high mechanical requirements due to underfloor impact scenarios [3]. Therefore, conventional CPMs are robust and stiff constructions [3][4][5]. In detail, an underbody cover (UBC) typically encases the functional unit comprising a safety case, a CPM housing and a shielding to prevent interactions of the electromagnetic field with  This paper presents the novel design of an ultrathin CPM providing a reduction of the z-dimension by 62 % compared to the referenced conventional CPM presented in [3]. Figure 2 qualitatively compares the referenced CPM in the conventional design and the proposed ultrathin design. The novel design was based on the spatial and functional integration of the receiver coil, ferrites and shielding in the underbody cover. Unlike conventional designs, this module was designed as a flexible structure, which required the application of novel materials. For this purpose, ferromagnetic sheets as alternatives to the conventional rigid and brittle ferrites were selected. These sheets and metal wire meshes, used as the shielding, were integrated into a plastic composite for the first time. Furthermore, the power electronics were excluded from the CPM housing, which was advantageous in terms of reducing the z-dimension. Additionally, the power electronics can be freely positioned in EVs. However, this results in special efforts in terms of contacting, shielding, and assembly, which were described comprehensively in [6].
The thermal management of CPMs featuring a high volumetric power density is a major challenge [7][8][9]. The local proximity of the coil, ferrites and shielding, being the main generators of the power losses of CPMs, is very challenging [7,10,11]. In order to support the development and optimization process and reduce elaborate experimental testing, accurate electromagnetic-thermal simulation models from the component to system level are required. As a consequence, several studies have focused on the simulation and experimental investigation of the electromagnetic and thermal behaviour of WPTSs. The simulation approaches can be categorized into order-reduced and coupled electromagneticthermal models. The order-reduced models of WPTSs reduce the complexity by analytically calculating the power losses of each component and running the simulation only in the thermal domain. Accordingly, various studies have presented approaches for the analytical quantification of the power losses of the coil [12][13][14][15], ferrites [12,14,15] and shielding [14,16]. The calculated power loss densities serve as heat sources for simulating the temperature fields. Based on this approach, good forecasting qualities of the electromagnetic-thermal behaviour can be achieved [10,14,15]. However, the modelling of the temperature-and frequency-dependent material properties has proven to be insufficient. In particular, the losses of the metal shields due to eddy current induction typically have to be taken from a previous magnetic simulation [14,15] or experimental measurement [10]. In addition, the prediction of the thermal distribution under the influence of external magnetic fields is only possible to a limited extent using order-reduced models [10,14]. To overcome these deficits, several authors have proposed coupled electromagnetic-thermal models, which bidirectionally link the simulation of the electromagnetic and thermal fields. Kim et al. [17] presented a two-way electromagnetic-thermal simulation analysis of a DD-Layout, which provided a very good forecasting quality compared to the experimental validation. However, the modelling of power losses and the thermal distribution of the individual components was not given. Niu et al. [18] provided a bidirectional coupled electromagnetic-thermal simulation of a GPM considering misalignment, but an experimental validation was not conducted. Alsayegh et al. [12] compared order-reduced modelling with coupledfield simulation. The results showed a good agreement of coupled-field simulation and the experiment, whereas the deviation between the order-reduced model and respective experimental studies was very large. A component-wise investigation of the thermal losses was not given. In summary, the simulation models found in the literature are either not sufficiently accurate regarding the transient or steady-state behaviour, cannot accurately predict total losses at the component level from the current densities or magnetic field densities or have not been experimentally validated at all. This paper determines the electromagnetic-thermal behaviour and efficiency of the proposed CPM design. The presented work approaches the investigation in a bottom-up manner, from the component to the system level, through simulations and experiments. Figure 3 illustrates the methodological approach of this paper.
At each level, specific electromagnetic-thermal simulation models were set up and experimentally validated. At the component level, the transformer components, the coil, ferromagnetic sheet and shielding, were examined individually to verify the electromagnetic and thermal behaviour of each component including the component-specific modelling of the power losses, as well as the boundary conditions. On this basis, the componentspecific substitutional models developed at the component level were combined at the module level. To verify the modelling without the influence of external fields, only the CPM was examined, neglecting the interaction with the GPM. Finally, the electromagnetic performance of the entire magnetic circuit was investigated, taking into account the CPM and GPM.
The paper is structured according to the methodological approach illustrated in Figure 3. Sections 2 and 3 provide the basis for the investigation. In this context, the proposed WPTS is detailed in Section 2, and the proposed ultrathin CPM design including its electromagnetic and thermal specification is described in Section 3. On this basis, Section 4 details the proposed electromagnetic-thermal simulation models at the component, module and system level. Section 5 describes the experimental setups for the validation of the respective simulations. In Section 6, the simulation results are validated with experimental measurements at the component, module and system level. On this basis, the forecasting quality of the proposed modelling is discussed with respect to the electromagnetic-thermal behaviour. Finally, Section 7 provides conclusions and future prospects.

Proposed Wireless Power Transfer System
The setup of the utilized WPTS is shown in Figure 4. The GPM used a reduced and simplified implementation of the power electronics without power factor correction (PFC) as the test system was operating on an isolated grid. The AC/DC stage was a silicon B6-rectifier, which transforms the 3-phase input voltage into DC feeding a DC link capacitor C DC,link . As the power inverter, a SiC-MOSFET half-bridge module was utilized, which switches up to 100 A RMS with up to a 100 kHz frequency. The serialserial-compensated resonance transformer consisted of the inductors of GPM L GPM and CPM L CPM , as well as the offboard and onboard compensation capacities C GPM and C CPM of the matching networks. As the onboard rectifier, a SiC-B4 bridge was used, which converts the AC from the resonance transformer into DC feeding a DC link capacity C filt . An actively controlled or passive DC/DC conversion was not implemented per the design as the output voltage of the CPM resonant tank was above the battery voltage. The power control strategy was considered to be a mix between frequency-control and interleave mode operation. However, in the experimental setup, the power was controlled manually via the output voltage of the isolated power supply.   Table 1 summarizes the main system specifications and design constraints of the exemplary WPTS. According to SAE J2954 [19], the onboard coil area was restricted to an area of 300 mm × 300 mm. Interoperability was required within lateral misalignments of x ± 75 mm and y ± 150 mm and air gaps between 100 mm and 210 mm. The DC link voltage of the CPM was fixed to a high-voltage (HV) battery voltage. The DC link voltage of the GPM was adjustable between 0 V and 540 V. Based on the system specifications and the design constraints, the geometric and material specifications of the transformer components of the CPM and GPM were defined. Figure 5 shows the setup of the exemplary WPTS. The transformer parameters are summarized in Table 2.  To maximize the given installation space while minimizing the z-dimension, the coil had a rectangular shape with 15 single-layer windings using a round stranded copper wire with 800 × 0.1 mm. As an alternative to the commonly used ceramic ferrite plates, a flexible 320 mm × 320 mm ferromagnetic sheet was used, providing a flexible and highly permeable foil. It was assembled from multiple, alternately stacked layers of ferromagnetic film, as well as nonmagnetic, adhesive and protective coatings made from polyethylene terephthalate (PET), as given in [20]. Overall, 30 layers were stacked to avoid saturation and high hysteresis losses, resulting in a total sheet thickness of 0.9 mm.
The lateral dimensions of the transmitting coil were 650 mm × 500 mm with a total number of 13 windings. The coil was underlaid with a layer of conventional BFM8 ferrite plates [21], which offer low loss characteristics in the required frequency range. The same 800 × 0.1 mm stranded wire was used on the primary side, although each winding was laid in double to ensure current carrying capacity (see Figure 5d).

Proposed Ultrathin CPM
This section describes the design and manufacturing process of the proposed ultrathin CPM and its electromagnetic and thermal specifications, as these parameter are key for the simulation of the electromagnetic-thermal behaviour.

Design and Manufacturing
In order to protect the active components of CPMs, mechanical designs are required to take the underfloor impact scenarios given in [3] into account. In contrast to conventional CPM designs with a high robustness and stiffness, this paper presents a novel ultrathin CPM design, which allows moderate deflections caused by mechanical bending loads. In this respect, the spatial and functional integration of the transformer components into the plastic composite of the underbody cover combined with the externalization of the power electronics enabled a significant reduction of the z-dimension [6]. The proposed design is schematically illustrated in Figure 6. The functionally integrated design of the proposed CPM was realized by utilizing the epoxy resin Araldite ® CW 2243-2L [22], providing a high thermal conductivity, which guarantees a sufficient structural connection and thermal coupling between the lossy integrated components. To ensure mechanical deflection, sufficient elongation of the integrated components was required, while ensuring the structural integrity of the plastic composite. Accordingly, the use of conventional aluminium shielding sheets acting as a separating layer and conventional ferrites with brittle fracture behaviour are not effective. Therefore, flexible ferromagnetic sheets substituted the conventional ferrites, and a metallic wire mesh was used as the shielding in the proposed CPM design. The demonstrator was scaled to the lateral dimensions of 900 mm × 600 mm with a z-dimension of 15 mm and a weight of 9.1 kg.
Utilizing the available tools, the manufacturing process of the proposed demonstrator was divided into 3 stages: (I) components, (II) module and (III) system. In the first stage, the components to be integrated, as well as the underbody cover as a long fibre-reinforced thermoplastic (LFT) product were produced. In the second stage, the coil, ferromagnetic sheet and epoxy resin were combined into a plastic-integrated functional unit using resin transfer moulding (RTM). Finally, in the third stage, the functional unit, the metal mesh and the underbody cover were combined with the epoxy resin using vacuum-assisted resin infusion (VARI). This procedure can be reduced to a 2-step process if an RTM tool of adequate size is available for the production of the CPM. Figure 7 illustrates the prototypical manufacturing process.

Electromagnetic and Thermal Specifications
In order to investigate the coupled electromagnetic-thermal behaviour of the proposed CPM design, the electromagnetic and thermal material properties of the CPM components were required. Table 3 summarizes the relevant properties of the CPM components for the investigation of the thermal management of the proposed design. As the capacitive effects were neglectable, the permittivity r was considered to be 1 in all cases. To enable modelling the ferromagnetic sheet and the copper mesh as solid regions, being a common simplification, the macroscopic material properties of these composite regions were derived from the properties of the individual materials by using the mixing rule from [23], if either isotropic behaviour can be expected or, in the case of conductivities, the materials act in parallel (e.g., in-plane characteristics for stacked materials). For materials acting in series (e.g., out-of-plane characteristics), we used the mixing rule: where χ represents the corresponding material property, respectively, χ C is the corresponding material property of the composite material and η is the volume fraction of the corresponding component. Deviations from this rule due to anisotropy are explained in Section 4.

Electromagnetic-Thermal Simulation
The simulation evaluation of the thermal management of the WPTS requires an electrically and thermally coupled analysis [8,16]. Although a two-way coupled electromagnetic-thermal simulation is computational more expensive, as both solvers converge iteratively to a common solution, it enables accurately predicting the transient temperature distribution of CPMs [12,17]. Therefore, each simulation result given in this paper was achieved from two-way electromagnetic-thermal simulations using the FEM simulation software Comsol Multiphysics. This section describes the modelling of the power losses of each transformer component, as well as the derivation of the substitutional simulation models at the component, module and system level. The frequency-and temperature-dependent loss density of the litz wire p litz ( f , T) is calculated from the sum of the individual losses due to the skin effect p skin ( f , T) and the proximity effect p prox ( f , T).
The skin effect losses and the proximity effect losses of the litz wire are defined as the sum of the losses of the litz wire strands. Due to the computational complexity, not every single strand was modelled. Instead, the loss density of a single strand was applied to the whole litz wire and weighted with the volumetric fill factor η litz . Accordingly, the skin effect loss density can be calculated as: where σ Cu (T) is the linear approximation of the temperature-dependent conductivity of copper, G skin ( f , T) is a frequency-and temperature-dependent factor modelling the skin effect [13] andĴ strand is the amplitude of the current density within a single strand. Note that J strand was considered constant across the cross-section of the conductor, as we assumed a perfect litz wire. The fill factor η litz and the current densityĴ strand are calculated via: where n strand is number of strands per litz wire, A strand is the cross-sectional area of a single strand and A litz is the cross-sectional area of the litz wire. The proximity effect losses of the litz wire are calculated in similar way. However, the losses are proportional to the square of the external magnetic field strength amplitudeĤ, which results from the solution of the magnetic simulation.
where G prox ( f , T) is a temperature-and frequency-dependent factor modelling the proximity effect [13]. The specific material properties for the thermal simulation of the litz wire are the density ρ litz , heat capacity c p,litz and thermal conductivity λ litz , calculated according to the rule of mixtures in Equation (1).

Ferromagnetic Sheet
The magnetisation curve of the ferromagnetic film material was nonlinear. Therefore, the average energy method [26] was used in order to calculate an adequate magnetic field density for the magnetic loss calculations. For a harmonic magnetic field excitation H(t) =Ĥcos(ωt), the time-averaged magnetic field density amplitudeB AE is calculated as: where T P is the period of the time harmonic excitation and B(H) is the nonlinear BH-curve.
The loss density of the ferromagnetic sheet p ferr is calculated from the given core loss curves from the material datasheet. Analogous to the litz wire model, both the power density and the magnetic field density need to be scaled by a volumetric fill factor η ferr ≤ 1, as the ferromagnetic sheet is a composite material. That is, the magnetic field density has to be scaled up, as the field is mostly concentrated in the highly permeable ferromagnetic layers and the power density has to be scaled down, as the losses are only located in the ferromagnetic layers: as p(B) is nonlinear, η ferr does not cancel out.
Observations from experiments have shown a nonisotropic in-plane loss behaviour. The origin of this deviating behaviour is unknown to the authors and has not been described in the material datasheet. Thus, the nonisotropic loss behaviour was captured with an adjustment factor γ ferr = 0.9, which was derived from the component-level experiments. The adjustment factor was applied to the in-plane magnetic field components B x and B y in opposite ways.

Shielding
The electromagnetic shielding was a mesh, which was woven from small copper strands embedded in the epoxy resin Aradur ® HY 2966 [22]. On a macroscopic level, a full-geometry simulation would require a huge computational effort. Thus, an equivalent substitute material model for the electromagnetic and thermal simulation was derived from a full-geometry model. Thermal and electric conductivity both describe a heat flow q or current flow J based on the gradient of the temperature T or the electric potential Φ.
In order to derive η mesh,cond , a section of the copper mesh was modelled in detail together with a block of identical dimensions, but with varying conductivity coefficients, as shown in Figure 8. Although the mesh was integrated in resin, its thermal λ Res or electric conductivity σ res can be neglected, since λ Res λ Cu and σ Res σ Cu . An external heat flux was applied on one edge, while on the opposite edge, the temperature was fixed. This enforced a temperature gradient according to the effective thermal conductance. The inplane conductivity coefficient η mesh,cond,x,y of the equivalent block model was then tuned, such that the temperature difference of the block matched the temperature difference of the mesh. In the out-of-plane direction (z), the conductivity coefficient is calculated as: where A mesh denotes the covered area of the copper mesh and A cu denotes the area of the mesh that is optically covered by copper (i.e., without air gaps). However, as the mesh was thin with respect to the height of the CPM, the error of this approximation was expected to be low.

Module Level
Using the material parameters given in Table 3 and the models derived and validated at the component level, a two-way transient electromagnetic-thermal simulation of the CPM was performed. Due to the integral design of the proposed CPM, the dimensions of the integrated components were derived from a roentgenogram, as shown in Figure 9. In order to reduce the model complexity, the simulation model took advantage of the axial symmetry by only simulating a half model, in which only important geometric features from the physical model, such as the curved shielding, were applied. At the module level, the coil was excited with a current of I litz = 25 A RMS and a frequency of f = 80 kHz. The heat transfer coefficient α was assumed to be 10 W m −1 K −2 on the top side and 5 W m −1 K −2 on the bottom side. Heat transfer by radiation was neglected according to [12,16]. The loss densities were calculated according to the material models and were applied to their specific volumetric domains.

System Level
In order to model the effects of an external magnetic field, as it is required in WPTSs through the transmitting coil of the GPM, a second coil on a layer of BFM8 [21] ferrite was added to the model according to the parameters in Table 2. To reduce the complexity, the GPM coil was modelled as a block with an equivalent coil model. The enforced coil current was I GPM = 34 A RMS and phase shifted by π/2 with respect to the CPM current, as this is the usual operating point. The air gap between the CPM and GPM was s = 150 mm.

Experimental Validation
In order to validate the simulation results, experimental evidence was required. Accordingly, the virtual experiments were re-enacted physically at the component, module and system level. All tests were executed on a hardware-in-the-loop (HIL) test bench, which was completely made out of plastics and wood to avoid measurement errors due to the influence of magnetic or conductive surrounding materials. The corresponding experimental setups are described in the following subsections.

Component Level
To validate the component-specific power loss modelling, the electromagnetic-thermal behaviour of each transformer component-coil, ferromagnet sheet and metal wire meshwas analysed separately. For this purpose, first the coil embedded in the resin was analysed alone, that is the coil was excited with a harmonic current (I 2 = 25 A RMS, f 0 = 80.5 kHz), while the heating process was captured with the thermographic camera VarioCAM HD head 600 [27]. In alignment with the simulation, the bottom side of the coil was isolated thermally with foam to enforce natural convection on the top side only. Moreover, this coil was utilized as a magnetic field excitation in the metal wire mesh test (I 2 = 25 A RMS, f 0 = 90.5 kHz) and the ferromagnetic sheet test (I 2 = 25 A RMS, f 0 = 88 kHz). Again, foam was used as an insulator to thermally decouple the components from the excitation coil. Power measurements were taken for every component test to further check for plausibility. These measurements were each based on the DC input power of the power electronics. The influence of static losses such as the power electronics or supply lines could be factored out due to the comparable operating frequencies and the known coil resistance. The experimental and measurement setups are shown in Figures 10 and 11. Due to the different thermal mass and the loss mechanisms, the measuring times varied greatly from 5 min for the metal wire mesh to 90 min for the coil until the steady-state temperature was reached. To increase the measurement accuracy of the infrared camera, tape was applied to the edges and the centre of the metal wire mesh to set the local material emissivity.
Tb2 Tb3 T b1 , T b2 , T b3 Figure 11. Experimental setup of the thermal investigation and position of the fibre-optic therma sensors.

Module Level
As shown in [10], module testing provides an appropriate substitution test to evaluate the thermal management of various lightweight designs for CPMs. Accordingly, the module testing determined the transient temperature distribution on the top and bottom of the proposed CPM design, taking into account only the CPM. Figure 11 illustrates the experimental setup of the module testing.
The CPM was excited with a harmonic current I coil = 25 A RMS, which was measured and monitored by a Rogowski coil of type PEM CWTMini HF06B [28]. Fibre-optic thermal sensors [29] were used for the thermal measurement on the bottom side, since their functionality is not affected by magnetic fields. To measure the temperatures as a result of the component-specific power losses, the sensors were positioned at the centre winding of the coil, at the edge of the ferromagnetic sheet and at the edge of the metal wire mesh. Additionally, the thermographic camera measured the surface temperature on the top side of the CPM. Although the system was designed to operate under forced air cooling, in this module-level test, natural convection was used, in order to reduce the complexity of both the experimental setup and the simulation model.

System Level
To validate the system simulation model and the module-level testing, the systemlevel testing determined the interaction of the CPM and GPM. The system was operated as a hardware-in-the-loop, which is shown in Figure 12. To ensure operation and testing under realistic electric conditions as given in automotive applications, an exemplary 400 V battery system was used. The power electronics was supplied by an isolated grid inverter, which provided a 3-phase grid with adjustable output voltage and frequency to enable the emulation of different grid types. The inverter itself was powered by the battery system, which closed the loop. The test bench was controlled via a measurement and control computer. This setup enabled operating either at a constant state of charge (i.e., charging mode) for efficiency tests or performing complete charging cycles (i.e., islanded mode). According to SAE J2954 [19], the CPM was fixed while the GPM was manually adjustable to examine the system behaviour during three-dimensional positioning. Analogous to the module testing, the temperature distributions were measured with the thermographic camera VarioCAM HD head 600 [27] on the top side and fibre-optic sensors [29] on the bottom side. The setup of the system testing is shown in Figure 13. In additional to the experimental investigation of the electromagnetic-thermal behaviour, efficiency measurements were performed "DC-to-DC", i.e., from the DC link voltage of the power electronics on the transmitting side to the battery voltage on the receiving side. Consequently, the input rectifier or the power factor correction of the primary power electronics was not accounted for. Misalignment and air gap requirements were considered according to SAE J2954 [19].

Results and Discussion
In the subsequent sections, the results of the simulations and the experiments performed are presented and compared at the component, module and system level. Finally, the results are summarized and discussed.  Table 4.  The experimental data of the ferromagnetic sheet showed the already mentioned anisotropic in-plane temperature distribution. As datasheet information on this anisotropy was not provided, a compensation parameter was used to approximate the observed behaviour. As shown in the simulation, the effect was captured qualitatively well; however, the transient time behaviour and the steady-state temperatures showed some differences. This might be related to bad measurement conditions resulting in a lower measured temperature compared to the simulation. In parallel, the modelling of the unexpected temperature anisotropy may also add to this observation. Further investigations and consultation with the manufacturer are necessary to mitigate this issue or to provide more accurate material data, respectively. However, the measured and the simulated power losses within the ferromagnetic material showed a good consistency with less than 10% lower losses in the simulation (see Table 4).

Component Level
Regarding the metal wire mesh, the simulated transient behaviour of the maximum temperature showed very good consistency with the measured infrared data. Furthermore, the observed and simulated heating pattern were consistent. Due to the applied tape (see Section 5.1), a cross-shaped heating pattern emerged in the measured infrared image.
From the experimental and simulation data, it can be concluded that the general validity of the proposed electromagnetic-thermal component-level models was successfully proven.

Module Level
The results of the module-level simulation and measurements are shown in Figure 15. The infrared image and the associated simulation were taken after 90 min of operation. Both temperature distributions on the surface, aside from geometric simplifications in the model, showed a good consistency. Small deviations can be observed at the edges of the component, most notably higher measured temperatures at the rear, marked as (×) and on the left and right side. The transient behaviour, however, both at the measurement points of top and bottom side, was captured accurately in the simulation. At most, the difference between the simulation and measurement was less than 5 K over the full 90 min operation time. On the basis of the module-level experiments and simulations performed, it can be concluded that the adoption of the prevalidated models from the component level to the module level was valid and gave accurate predictions.

System Level
The results for the system-level simulation and validation are shown in Figure 16. In order to reach a sufficient steady-state approximation, the CPM was operated in the system for 120 min. As expected, in comparison to the module-level test, the average measured temperature level increased due to the additional magnetic field of the GPM. Again, both the simulation and measurement temperature distributions showed a good consistency with the exception of the rear edge, marked as (×). This hotspot was induced by the high magnetic field densities of the combined magnetic field of the CPM and GPM, which were concentrated locally due to the shielding effect of the copper wire mesh and the field guidance of the ferrite sheet. This effect could be observed in the simulation as well; however, the magnitude could not be reproduced. Since these losses occurred in the shielding, it can be concluded that the proposed copper mesh model is, as a result, not yet readily applicable to the system-level simulation. Thus, further research is required to specifically model the edge effects on metal wire meshes. It is suspected that locally occurring hotspots are strongly dependent on the local geometric details and sensitive to the effective local conductivity. Besides the deviations on the hotspot, the transient and steady temperature behaviour was well predicted by the simulation and consistent with the measurements. In total, after 120 min, the CPM reached a maximum temperature of 80°C with a maximum hotspot temperature of 115°C on the top side and 90°C on the bottom side at a power level of 7.2 kW. The results of the efficiency measurements are shown in Figure 17. The proposed CPM design achieved its highest efficiency of 92% at transmission distances of 115 mm and 130 mm. Both efficiency maps are very similar, with slightly decreased efficiency for the lower distance. This might be related to an out-of-tune resonant tank on the primary side, which increased the GPM current and thus losses. For the highest transmission distance of 155 mm, the efficiency dropped to 89% in the centre and to 75% at the maximum offsets, which was related to a very low coupling between the CPM and GPM. The observed efficiency asymmetry regarding the y-axis can be explained with a tilted GPM imitating an inclined parking position (rotation along the x-axis). Efficiency η max (%) Figure 17. Results of the efficiency analysis depending onthe air gap s.
In conclusion, despite the local hotspot and its impact on the surrounding region, the system-level simulation was consistent with the measurements, as the remaining top and bottom measurements showed. Thus, the proposed electromagnetic-thermal simulation models allowed an accurate prediction of the transient temperature distribution and the losses within the CPM components under the aforementioned restrictions. This provides the basis for a reliable WPTS design and optimization procedures in the future.

Conclusions
This paper presented a systematic investigation of the electromagnetic-thermal behaviour of a novel ultrathin onboard receiver module for vehicular wireless power transfer systems. The studied module was designed as an integral design, in which a receiver coil, a ferromagnetic sheet as a alternative to conventional ferrites and a metal wire mesh as the shielding were spatially and functionally integrated. The module is capable of receiving a power of up to to 7.2 kW at a transmission frequency of f 0 = 85 kHz with a maximum efficiency of 92% over a distance of 110 mm to 160 mm. Special emphasis was given to the validation of a novel two-way electromagnetic-thermal simulation model, which was developed systematically considering the component, module and system level. At each level, the simulation results were consistent with the experimental validation. However, at the component level, the transient time behaviour and the steady-state temperature of the ferromagnetic sheet must be improved quantitatively. However, the module simulation accurately reproduced the electromagnetic-thermal behaviour. At the system level, the simulation and measurement temperature distribution showed a good qualitative and quantitative consistency. Thus, the proposed two-way electromagnetic-thermal simulation model provided a significant contribution to a more sophisticated dimensioning of receiver modules, especially in regard to their integration in vehicle structures. However, future work is required to improve the quantitative prediction of local hotspots in the peripheral area of metal wire mesh.

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