Next Article in Journal
Bending and Torsional Stress Factors in Hypotrochoidal H-Profiled Shafts Standardised According to DIN 3689-1
Previous Article in Journal
Formalising Autonomous Construction Sites with the Help of Abstract Mathematics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of One-Way and Two-Way Coupled Simulation for Thermal Investigation of Vehicular Wireless Power Transfer Modules †

1
Mercedes-Benz AG, 059/L423, 71059 Sindelfingen, Germany
2
Institut für Leichtbau und Kunststofftechnik, Technische Universität Dresden, 01307 Dresden, Germany
*
Author to whom correspondence should be addressed.
This is an extended article from a published conference paper.
These authors contributed equally to this work.
Eng 2023, 4(1), 816-828; https://doi.org/10.3390/eng4010049
Submission received: 31 January 2023 / Revised: 18 February 2023 / Accepted: 25 February 2023 / Published: 5 March 2023
(This article belongs to the Section Electrical and Electronic Engineering)

Abstract

:
The thermal management of the vehicular module is key to the design of efficient wireless power transfer systems. In order to predict the thermal behavior by simulation, the mutual interaction of the electromagnetic and thermal fields must be taken into account. This multiphysical coupling leads to extensive computational effort. One approach to reduce the complexity by limiting the interdependencies between the domains is one-way coupling. This paper examined the applicability of one-way and two-way coupling for the prediction of the thermal management of an exemplary vehicular wireless power transfer module. The electromagnetic–thermal behavior of the proposed module was systematically studied by experiments and simulations on the component and module level. The performed studies showed that both simulation approaches accurately capture the transient thermal behavior of the coil and ferrites on the component level, whereas the one-way coupled simulation underpredicts the power losses caused by shielding by more than 20%, leading to a steady-state temperature difference of 15 K. As a result, on the module level, the two-way coupled simulation model provides a more accurate representation of the electromagnetic–thermal behavior of the proposed car pad module. Therefore, the authors recommend using a two-way coupling approach for the thermal dimensioning of wireless power transfer modules for electric vehicles.

1. Introduction

A key aspect for the adoption of electromobility is the improvement of the convenience in the charging of electric vehicles (EVs). In this regard, wireless power transfer systems (WPTSs) provide a very user-friendly charging technology by wirelessly transferring power from an offboard ground pad module (GPM) to an onboard car pad module (CPM) without the need for manual intervention. The setup of an exemplary WPTS is illustrated in Figure 1.
This paper focused on the onboard CPM, which is usually mounted on the underside of EVs due to the operation principle. In particular, the limited installation space in the vehicle underbody requires CPM designs to provide high volumetric and gravimetric power densities combined with low vertical dimensions. For this reason, the authors already presented a sandwich design in [1].
Thermal management of high-volumetric-power-density CPMs is a major challenge due to the local proximity of the major transformer components—coil, ferrites, and shielding—which cause the largest power losses. Bosshard and Kolar [2] conducted a multi-objective optimization study on a 50 kW/85 kHz wireless power transfer system (WPTS) by analyzing the power loss distribution of the transformer components. Their findings indicated that the total losses attributed to the coil, ferrites, and shielding exceeded 200 W. Consequently, the prototype was designed with active air cooling, which resulted in a maximum steady-state temperature of 40 C in the coil winding. Moghaddami and Sarwat [3] performed a multi-physics analysis of a 22 kW WPTS based on a 2D simulation model. The results indicated a steady-state temperature over 70 C under forced air cooling conditions. The thermal management of different CPM designs, including the proposed sandwich design, was investigated by the authors of this paper under natural convection conditions in [1]. Depending on the CPM design, the results showed that the maximum temperatures surpassed the maximum operating temperature limit of 130 C for some plastic components used.
In order to develop efficient CPMs without having to perform a large number of experimental tests, accurate electromagnetic–thermal simulation models are required. For this reason, the authors of this paper investigated the prediction quality of one-way and two-way coupled simulations using the sandwich design as an example. However, this, study given in [4], was limited to the module level only considering the CPM without the influence of the external magnetic field caused by the GPM. In detail, the substitutional modeling of the transformer components was verified by experiments at the module level. However, detailed component-level simulations to derive this substitution model were not given therein.
This paper examined the applicability of one-way and two-way coupled simulation models to predict the electromagnetic–thermal behavior of WPTSs for electric vehicles. The paper aimed to provide a quantitative comparison of the respective prediction accuracies based on the analysis of an exemplary CPM design. The presented work approached the investigation in a two-step manner from the component to the module level. First, at the component level, the main transformer components—coil, ferrites, and shielding—were investigated individually to verify the electromagnetic and thermal behavior of each component. This included the component-specific modeling of the power losses, as well as the thermal boundary conditions. The second step was the combination of the component-level models at the module level. At this level, only the proposed CPM design was examined. The interaction with the GPM was neglected in order to verify the modeling without the influence of external fields.
This paper is structured into seven sections based on the study of related works provided in Section 2. Section 3 details the proposed CPM design in terms of the construction, electromagnetic specifications, and properties of the used materials. Section 4 describes the proposed one- and two-way coupled simulation approaches including modeling the power losses of the main transformer components. Section 5 presents the experimental setups for the validation of the simulation models. Section 6 presents the results of the one- and two-way coupled simulations, as well as the experimental measurements at the component and the module level. On this basis, the accuracy of the proposed models is discussed with respect to the electromagnetic–thermal behavior. Finally, Section 7 provides conclusions and future prospects.

2. Related Works

Coupled simulation approaches are key to predict electromagnetic–thermal behavior. In particular, two-way coupled simulations, linking the electromagnetic and thermal domain bidirectionally, provide very good consistency between simulation and experiment. However, this approach results in high computational costs and efforts due to the complexity of the models. One approach to improve computational efficiency is to use a one-way coupling limiting the exchange of information between the separately calculated electromagnetic and thermal fields. Thus, only the information of the electromagnetic field is considered for the calculation of the thermal field. Figure 2 shows a schematic comparison of one-way and two-way coupling for calculating electromagnetic–thermal behavior.
In order to predict the thermal management of WPTSs, only a few authors have examined one-way coupled approaches. Table 1 summarizes the relevant research. Kim et al. [5] performed a coupled electromagnetic–thermal simulation of a GPM with a double-D coil layout and achieved a very good agreement of the electromagnetic–thermal behavior compared to the experimental validation. However, the modeling of power losses and the thermal distribution of the individual components was not given. Zhang et al. [6] presented a comprehensive investigation of the thermal performance of a 30 kW WPTS. This study was performed on the system level under the consideration of CPM and GPM including misalignment scenarios and achieved a good consistency of the simulation and experimental validation. However, the temperature of the metal shielding was not given. The authors of this paper already studied the thermal management of the proposed CPM design on the module level. The study showed that power losses in metal parts can cause a significant discrepancy between the experiment and one-way coupled simulation.
Regarding two-way coupling, several authors have performed detailed studies for predicting the electromagnetic–thermal behavior of WPTSs. Amirpour et al. [8] found good agreement between the numerical and experimental results analyzing an in-road GPM on the module level. Bertoluzzo et al. [10] presented a magneto-thermal analysis of a WPTS under the consideration of a car frame. However, the analysis and experimental validation were not given. Moghaddami and Sarwat [3] performed a multiphysics analysis of 22 k W WPTS coupling time-harmonic electromagnetic, time-dependent, or steady-state thermal and fluid models. However, an experimental validation was not conducted in these studies. Rasekh et al. [9] conducted a temperature analysis for three substantial pad designs on the system level including misalignments. The results showed a very good agreement between the simulation and experiment. The authors of this paper provided a two-way coupled electromagnetic–thermal simulation model to investigate the electromagnetic–thermal behavior of a functionally integrated CPM with an innovative magnetic core and shielding material [7]. The study included the systematic modeling and experimental validation of the temperature- and frequency-dependent material properties at the component, module, and system levels. The results showed a very good consistency of the proposed simulations with the experimental validation on all levels.
In summary, both one-way and two-way coupling can be efficient approaches to predict the electromagnetic–thermal behavior of WPTSs. However, only the authors of this paper gave the transient temperature rise of the CPM components being crucial to examine the prediction quality of the simulation approach. To the knowledge of the authors, a quantitative comparison of the one-way and two-way coupled simulation strategies for calculating the electromagnetic–thermal behavior of automotive WPTSs was only given in [4]. However, this work was limited to the module level, without analyzing the prediction quality on the component level.

3. Description of Proposed CPM Design

Solid and stiff CPM designs provide benefits due to potential impact scenarios associated with the installation position in the vehicle underbody. Accordingly, a novel CPM design with an integrated 11 kW WPTS was already proposed as a sandwich design in [1]. For this reason, only a brief description of the demonstrator is given herein, including the electromagnetic specifications and material properties of the components used. However, the proposed demonstrator is intended for experimental testing under laboratory conditions. Therefore, a simplified generic prototype was used scaled to the dimensions of 390 mm × 390 mm length and width and a height of 30 mm.

3.1. Sandwich Design

Figure 3 shows the proposed CPM design and the dimensions of the magnetic active components. The design was based on a sandwich core using Rohacell 71 IG F [11], which serves as carrier of the active components. The coil area was dimensioned to a rectangular shape of 300 m m width × 300 m m length. The coil consists of 15 single-layer windings of a round stranded copper wire made of 840 strands with a diameter of 0.1 mm per strand. The ferrite layout consists of BFM8 [12] plates (length (l) × width (w) × thickness (t): 49 m m × 49 m m × 3 m m ). The ferrites completely cover the winding area of the coil. Additionally, 3F36 ferrites [13] are aligned and stacked along the entire inner and outer perimeter. To improve thermal coupling, the coil and the ferrites are embedded into the sandwich core using the thermal conductive resin CW 2243-2L [14]. Furthermore, thermal gap filler pads [15] (l × w × t: 145 m m × 145 m m × 0.5   m m ) improve the thermal transfer between ferrites and shielding, whereas the shielding is glued on top of the sandwich core using Betaforce 2850 [16]. Additionally, aluminum heat sinks (l × w × t: 145 m m × 20 m m × 10 m m ) are glued to the shielding using ACC Silicones SILCOTHERM [17]. On the bottom side, an organic sheet [18] is glued to the backside of the sandwich core using Betaforce 2850 [16].

3.2. Material Properties

In order to investigate the electromagnetic–thermal behavior of the proposed CPM designs, the electromagnetic and thermal material properties of the CPM components are required. Table 2 summarizes the relevant properties of the CPM components, which are necessary 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 the modeling of the coil component as a solid region, the macroscopic material properties of this composite region were derived from the properties of the individual materials by using the mixing rule from [19]:
χ C = η χ 1 + 1 η χ 2 ,
where χ 1 and χ 2 represent the component’s material properties, respectively, χ C is the corresponding material property of the composite material, and η is the volume fraction of the corresponding component. The resulting material properties of the CPM components are also provided in Table 2.

4. Simulation Model

The electromagnetic–thermal simulation, both one-way and two-way coupled, was calculated in Comsol Multiphysics. First, one- and two-way coupled simulations were performed on the component level. This means that the models for the litz wire, ferrites, and shielding were implemented, simulated, and validated before all components were combined into a final electromagnetic–thermal model.
The following sections describe the modeling of the magnetic and thermal behavior, with particular emphasis on the coupling between the temperature, the electromagnetic properties, and the power dissipation mechanisms.

4.1. Litz Wire

It was shown in [7] that the frequency- and temperature-dependent loss density of litz wire p litz ( f , T ) can be accurately calculated from the sum of skin effect losses p skin ( f , T ) and proximity effect losses p prox ( f , T ) .
p litz ( f , T ) = p skin ( f , T ) + p prox ( f , T )
The litz wire’s skin and proximity effect loss densities were calculated from the loss densities of a single strand, which were each multiplied by the volumetric litz wire fill factor η litz ; see Equations (3) and (4).
p skin ( f , T ) = η litz σ Cu ( T ) G skin ( f , T ) J ^ strand 2 ,
p prox ( f , T ) = η litz σ Cu ( T ) G prox ( f , T ) H ^ 2 A strand 2 ,
G skin ( f , T ) and G prox ( f , T ) are temperature- and frequency-dependent, unit-less factors modeling the skin and proximity effect, and σ Cu ( T ) is the linear approximation of the temperature-dependent electrical conductivity of copper. J ^ strand represents the amplitude value of a single strand’s current density, while H ^ is the amplitude value of the magnetic field strength H = ( H x , H y , H z ) to which a single strand is exposed. The volumetric litz wire fill factor η litz and the current density J ^ strand are calculated following Equations (5) and (6):
η litz = n strand A strand A litz = 0.41 ,
J ^ strand = 2 J strand = 2 I litz A strand n strand ,
where n strand denotes the number of strands per litz wire, A strand is the cross-sectional area of a single strand, A litz is the cross-sectional area of the litz wire, and I litz is the RMS value of the litz wire current.

4.2. Ferrites

The initial magnetic permeability μ i of the ferrite materials is strongly temperature dependent in the expected operating temperature range, albeit at an already high permeability level. In the previous work [4], the temperature-dependent permeability was expected to have a small effect; hence, the average energy method [22] was applied to account for a high flux density nonlinearity. However, the investigations carried out revealed that the flux densities are in the linear permeability range. Therefore, in this paper, the permeability was modeled as an interpolated look-up function μ i ( T ) . This was performed according to the curves provided in the datasheet. The loss density of the ferrites p core ( T ) is calculated from the frequency f and the amplitude value of the magnetic field strength B ^ according to the temperature-dependent Steinmetz equation [23] following Equation (7):
p core ( T ) = C m · f x · B ^ y · C t 2 T 2 C t 1 T + C t 0 ,
where C m , x , y , C t 2 , C t 1 , C t 0 are coefficients, which either are determined by the manufacturer or need to be approximated from the provided manufacturer datasheet. In the case of 3F36, the Steinmetz coefficients were supplied by the manufacturer, whereas the coefficients for BFM8 were determined from the material datasheet supplemented by the experimental data. The coefficients are listed in Table 3.

4.3. Shielding

The power loss density of the shielding material p sh was calculated from the current density J and the material’s electrical conductivity σ :
p sh ( T ) = 1 σ ( T ) J 2
The electrical conductivity σ ( T ) was modeled to be temperature dependent. This was realized by using linear approximation:
σ ( T ) = σ 0 1 + α 0 ( T T 0 ) ,
where α 0 is the temperature coefficient of the specific material, T 0 is the reference temperature, and σ 0 is the electrical conductivity at the reference temperature.

5. Experimental Setups

The experimental validation was divided into two steps: component and module testing. In the first step, at the component level, each active component of the CPM—coil, ferrites, and shielding—was investigated separately. In the second step, at the module level, the electromagnetic–thermal behavior of the proposed CPM was studied, neglecting the interaction with the GPM. The tests were executed on a hardware-in-the-loop (HiL) test bench for WPTS at TU Dresden. Fiber-optic temperature sensors [24] were used to measure the local temperatures, since their functionality is not affected by magnetic fields.

5.1. Component Testing

The experimental setups are shown schematically in Figure 4. For the validation of the component-specific power loss modeling, the electromagnetic–thermal behavior of each transformer component—coil, ferrites, and metallic shielding—was analyzed separately; see Figure 4a. For this purpose, the coil embedded in the resin was analyzed as a standalone component. The coil was excited with a harmonic current (current: I = 25 RMS, frequency: f = 80.5   kHz ), while the heating process was captured by a thermographic camera, the VarioCAM HD head 600 [25]. In order to force cooling by natural convection only on the top side, the bottom side of the coil was thermally insulated with foam. In addition, this coil was used as a magnetic field excitation for the investigation of the BFM8 ferrites ( I = 25 RMS, f = 88   kHz ) and the aluminum shield with the heat sinks ( I = 45 RMS, f = 85   kHz ), as shown in Figure 4b,c. In the shielding setup, foam was used as the thermal insulator to thermally decouple the components from the excitation coil. The harmonic current excitation of the coil was achieved by a series connection with a bank of capacitors. As the coil inductance varied during the component-level testing, different resonance frequencies occurred, explaining the three different excitation frequencies for the ferrites, shielding, and coil testing.

5.2. Module Testing

Figure 5 shows the experimental setup for module-level testing. By using module testing, the transient temperature distribution of the active components of the proposed CPM design was determined, while neglecting the influence of external fields. As shown in [1], this substitute testing is reasonable for investigating the thermal performance of CPMs. Assuming minimum battery voltage implies the most-challenging case for the thermal management of CPMs, since this results in the maximum CPM current I CPM . Considering this scenario, the proposed CPM was excited with a harmonic current I CPM = 33   A RMS ( f = 81 kHz ) measured and monitored by a Rogowski coil of type PEM CWTMini HF06B [26]. To measure the transient temperature rise of the magnetic active components, fiber-optic thermal sensors [24] were attached to the coil ( T coil ), ferrites ( T core ), and shielding ( T sh ). Furthermore, the thermographic camera, the VarioCAM HD head 600 [25], was used to measure the surface temperature T surf on the top side of the CPM. The system was considered without forced air cooling assuming natural convection.

6. Results and Discussion

6.1. Component Level

Figure 6 and Figure 7 show the results of the component-level simulations compared to the measurements. Due to the differences in the thermal mass and loss mechanisms, there was a wide variance in the measurement times to the steady-state (Figure 7). With respect to the coil and the ferrites, it can be observed that both the one-way and the two-way coupled simulation models resulted in an almost identical temperature distribution at the steady-state, which was in good agreement with the recorded thermographic image from the experiment (Figure 6). Since there was no feedback from the temperature field to the magnetic field, the exact alignment of the one-way and two-way coupled coil simulations was expected. In contrast, a strong temperature dependence of the initial permeability μ i was implemented for the ferrite material, resulting in only a small deviation of the observed temperature field. The transient temperature behavior of both models, calculated by averaging the surface temperature, matched the experimental data (Figure 7). Consequently, the transient curves of the calculated power losses P coil and P core of both simulation approaches agreed qualitatively and quantitatively.
For the shielding component, the two-way coupled simulation model predicted the transient behavior of the mean surface temperature very well compared to the measurements (Figure 7). As a result, the steady-state temperature distribution matched the experimentally recorded thermographic image (Figure 6). In contrast, the simulation results showed that the one-way coupled simulation model produced temperatures that were too low for the majority of the experimental time. Thus, both the transient behavior and the steady-state distribution of the top surface deviated by up to 20 K compared to the experimental data. The quantitative comparison of the time-dependent power losses P sh shown in (Figure 7) indicated significant differences between the simulation models. In both cases, a steady-state value was reached within the operating time. However, over the course of the experimental time, the calculated shielding losses increased using two-way coupling and decreased using one-way coupling. This was expected due to the positive temperature coefficient of the shielding material, which caused the electrical resistivity of the material to decrease as the temperature increased. This resulted in an increase in the skin depth, changing the eddy current distribution in the two-way coupled simulation. As a consequence, the total losses increased as well. On the contrary, the one-way coupled simulation had a fixed magnetic field solution, resulting in a constant induced electric field. Thus, an increase in material resistivity resulted in a decrease in induced eddy-currents and losses, respectively. However, the assumed heat transfer coefficient h of 10 W m 2 K 1 on the top surface was considered to be valid due to consistency of the two-way coupled simulation model compared to the measurements.

6.2. Module Level

Figure 8 shows the results of the component-level simulations in comparison with the measurements. In contrast to both simulation results, the transient temperature behavior of the measurement results indicated that the steady-state was not reached after the maximum operating time of 120 min. In most cases, a steady-state comparison was representative for the application of the CPM. However, the authors contend that the use of transient temperature characteristics for model validation is preferable, despite the increased computational cost. This is due to the fact that modeling errors are more easily detected using transient simulations. Furthermore, a transient validated model provides a more accurate representation of the system’s dynamic behavior, enabling the evaluation of different operating scenarios and conditions. Additionally, transient simulations are useful for investigating thermal stress scenarios that may not be apparent in steady-state simulations, while still providing a reasonable prediction of the steady-state temperature under complex heat transfer coefficient conditions. Therefore, only the experimental data resulting from an operating time of 120 min were collected for validation.
Comparing the simulation results, it can be observed that the two-way coupled simulation model provided a better representation of the experimental measurements than the one-way coupled model. Despite not reaching the measured temperatures after 120 min, the two-way coupled model demonstrated a closer agreement with the experimental data. In contrast, the one-way coupled model consistently predicted lower temperatures for the majority of the experimental duration, although the difference between both models was never higher than 5 K. A quantitative comparison indicated that there was a notable effect of additional shielding losses in the two-way simulation, as the highest deviation in temperature of almost 5 K was observable in the shielding after a simulation time of 120 min. The coil and ferrite exhibited similar behavior, but with a reduced deviation. This observation was consistent with the component-level analysis. In the case of the coil, this can be attributed to the lack of thermal–electromagnetic feedback in the one-way coupled model, and although thermal–electromagnetic feedback was implemented for the ferrites, the effect on the temperature or dissipated power was small, as was already observed at the component level. This can also be seen in Figure 9, in which the time-dependent power losses are compared componentwise. While the coil and ferrites only showed a slight difference in losses over the simulated time, the shield losses diverged right from the beginning. Furthermore, the differences in the coil and ferrite losses can be explained by the generally diverging temperature of the exemplary CPM caused by the shielding losses.
The computational time for the one-way coupled simulation was 1 h and 34 min, whereas the two-way coupled simulation required 8 h and 13 min. Both simulations were performed on an AMD Ryzen 3700X 8-core computer with 46 GB of RAM and using equivalent physics and meshing.
Overall, it was shown that the two-way coupled simulation yielded more accurate results than the one-way coupled simulation, which can mostly be attributed to the diverging power losses of the shielding. However, in combination with other heat sources at the module level, the overall temperature error may be reduced. Furthermore, despite the utilization of a two-way coupled simulation model, the authors acknowledge that the model does not perfectly capture the transient electromagnetic–thermal behavior of the examined CPM yet. The deviations observed were attributed to either inaccuracies in the identification of the material properties or an incorrect estimation of the material quantity within the CPM. Both factors have a significant impact on the total heat capacity and, therefore, on the transient thermal behavior. Additionally, the assumption of a constant heat transfer coefficient is a simplification, as it is a nonlinear function of temperature in reality. Furthermore, the components were validated at lower steady-state temperatures than the module, and higher temperatures may reveal other effects not captured in the transfer coefficient.

7. Conclusions

This paper presented a quantitative comparison of a one-way and a two-way coupling for predicting the electromagnetic–thermal behavior of an exemplary CPM. The investigated CPM was designed as a sandwich and capable of receiving a transmission power of 11 kW at a transmission frequency of 85 kHz. The one-way and the two-way simulation models were developed systematically at the component and module level. At the component level, the results showed that both simulation approaches accurately captured the transient thermal behavior of the coil and ferrites, but the one-way coupled simulation underpredicted the temperature of the shielding. This discrepancy was attributed to the strong impact of the shielding material’s positive temperature coefficient on the magnetic field simulation. The two-way simulation model reproduced this effect better, which was also evident at the module level. As a result, the one-way coupled simulation model produced lower temperatures due to the divergent calculation of shielding losses compared to the two-way coupled simulation. In conclusion, this study found that, in order to correctly predict the thermal behavior of applications with conductive components in alternating magnetic fields, a two-way coupled simulation is advantageous. Further work will include a review of the deviations in transient temperature prediction and a system-level investigation, including the operation with a ground pad module.

Author Contributions

Conceptualization, S.Z.; software, M.H.; validation, M.H.; writing—original draft preparation, S.Z. and M.H.; writing—review and editing, A.W.; visualization, S.Z.; supervision, N.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zimmer, S.; Helwig, M.; Lucas, P.; Winkler, A.; Modler, N. Investigation of Thermal Effects in Different Lightweight Constructions for Vehicular Wireless Power Transfer Modules. World Electr. Veh. J. 2020, 11, 67. [Google Scholar] [CrossRef]
  2. Bosshard, R.; Kolar, J.W. Multi-Objective Optimization of 50 kW/85 kHz IPT System for Public Transport. IEEE J. Emerg. Sel. Top. Power Electron. 2016, 4, 1370–1382. [Google Scholar] [CrossRef]
  3. Moghaddami, M.; Sarwat, A. Time-Dependent Multi-Physics Analysis of Inductive Power Transfer Systems. In Proceedings of the 2018 IEEE Transportation Electrification Conference and Expo (ITEC), Long Beach, CA, USA, 13–15 June 2018; pp. 130–134. [Google Scholar] [CrossRef]
  4. Zimmer, S.; Helwig, M.; Winkler, A.; Modler, N. One-way vs. two-way coupled simulation: Investigation of thermal management of wireless power transfer modules for electric vehicles. In Proceedings of the 2022 Wireless Power Week (WPW), Bordeaux, France, 5–8 July 2022; pp. 84–89. [Google Scholar] [CrossRef]
  5. Kim, S.; Amirpour, M.; Covic, G.; Bickerton, S. Thermal Characterisation of a Double-D Pad. In Proceedings of the 2019 IEEE PELS Workshop on Emerging Technologies: Wireless Power Transfer (WoW), London, UK, 18–21 June 2019; pp. 1–5. [Google Scholar] [CrossRef]
  6. Zhang, B.; Deng, J.; Wang, W.; Li, L.; Wang, Z.; Wang, S.; Guidi, G. Multi-Objective Thermal Optimization Based on Improved Analytical Thermal Models of a 30 kW IPT System for EVs. IEEE Trans. Transp. Electrif. 2023, 9, 1910–1926. [Google Scholar] [CrossRef]
  7. Helwig, M.; Zimmer, S.; Lucas, P.; Winkler, A.; Modler, N. Multiphysics Investigation of an UltrathinVehicular Wireless Power Transfer Module for Electric Vehicles. Sustainability 2021, 13, 9785. [Google Scholar] [CrossRef]
  8. Amirpour, M.; Kim, S.; Battley, M.P.; Kelly, P.; Bickerton, S.; Covic, G. Coupled electromagnetic–thermal analysis of roadway inductive power transfer pads within a model pavement. Appl. Therm. Eng. 2021, 189, 116710. [Google Scholar] [CrossRef]
  9. Rasekh, N.; Dabiri, S.; Rasekh, N.; Mirsalim, M.; Bahiraei, M. Thermal analysis and electromagnetic characteristics of three single-sided flux pads for wireless power transfer. J. Clean. Prod. 2020, 243, 118561. [Google Scholar] [CrossRef]
  10. Bertoluzzo, M.; Forzan, M.; Di Barba, P.; Mognaschi, M.E.; Sieni, E. Magneto-Thermal Analysis of a Wireless Power Transfer System. In Proceedings of the 2022 IEEE 20th Biennial Conference on Electromagnetic Field Computation (CEFC), Denver, CO, USA, 24–26 October 2022; pp. 1–2. [Google Scholar] [CrossRef]
  11. Evonik Resource Efficiency GmbH. Product Information: ROHACELL® IG-F; Evonik Resource Efficiency GmbH: Hanau-Wolfgang, Germany, 2020. [Google Scholar]
  12. Blinzinger Elektronik GmbH. Data Sheet: Ferrite Core BMF8; Blinzinger Elektronik GmbH: Forchtenberg, Germany, 2020. [Google Scholar]
  13. Digi-Key Electronics Germany GmbH. Data Sheet: PLT32/20/3.2/R-3F36; Digi-Key Electronics Germany GmbH: Munich, Germany, 2020. [Google Scholar]
  14. Huntsman Advanced Materials GmbH. Advanced Materials: Araldite® CW 2243-2 L, Aradur® HY 2966; Huntsman Advanced Materials GmbH: Bergkamen, Germany, 2012. [Google Scholar]
  15. ARCTIC GmbH. Thermal Pad: High Performance Gap Filler; ARCTIC GmbH: Braunschweig, Germany, 2020. [Google Scholar]
  16. DuPont de Nemours, Inc. Betaforce; DuPont de Nemours, Inc.: Wilmington, DE, USA, 2020; Available online: https://www.dupont.materialdatacenter.com (accessed on 30 January 2023).
  17. ACC Silicones Ltd. SILCOTHERM®: Werkstoffe zur Wärmeleitung; ACC Silicones Ltd.: Essen, Germany, 2020. [Google Scholar]
  18. LANXESS Deutschland GmbH. Material Data Sheet: Tepex Dynalite 102-RG600(x)/47%: Roving Glass—PA6 Consolidated Composite Laminate; LANXESS Deutschland GmbH: Leverkusen, Germany, 2020. [Google Scholar]
  19. Ashby, M.F. Composite Materials, Microstructural Design of. In Encyclopedia of Materials: Science and Technology; Elsevier: Amsterdam, The Netherlands, 2001; pp. 1357–1361. [Google Scholar] [CrossRef]
  20. Deutsches Kupferinstitut. Data Sheet: Cu-ETP; Deutsches Kupferinstitut: Düsseldorf, Germany, 2020. [Google Scholar]
  21. Schwarzwälder Metallhandel GmbH. Data Sheet: EN AW 5754 [EN AW-Al Mg3]; Schwarzwälder Metallhandel GmbH: Trossingen, Germany, 2020. [Google Scholar]
  22. Paoli, G.; Biro, O.; Buchgraber, G. Complex representation in nonlinear time harmonic eddy current problems. IEEE Trans. Magn. 1998, 34, 2625–2628. [Google Scholar] [CrossRef]
  23. Steinmetz, C. On the law of hysteresis. Proc. IEEE 1984, 72, 197–221. [Google Scholar] [CrossRef]
  24. Weidmann Technologies Deutschland GmbH. TS2: Fiber Optic Temperature Sensor; Weidmann Technologies Deutschland GmbH: Dresden, Germany, 2020. [Google Scholar]
  25. InfraTec GmbH. VarioCAM® HD Head: Infrarot-Thermografiesystem für den Industrie- und Laboreinsatz; InfraTec GmbH: Dresden, Germany, 2020. [Google Scholar]
  26. Cosinus Messtechnik GmbH. Data Sheet: CWT Mini; Cosinus Messtechnik GmbH: Taufkirchen, Germany, 2020. [Google Scholar]
Figure 1. Exemplary wireless power transfer system.
Figure 1. Exemplary wireless power transfer system.
Eng 04 00049 g001
Figure 2. Comparison of coupling approaches: (a) one-way coupling; (b) two-way coupling.
Figure 2. Comparison of coupling approaches: (a) one-way coupling; (b) two-way coupling.
Eng 04 00049 g002
Figure 3. Proposed CPM design: (a) sandwich design and (b) dimensions (in mm) of the magnetic active components.
Figure 3. Proposed CPM design: (a) sandwich design and (b) dimensions (in mm) of the magnetic active components.
Eng 04 00049 g003
Figure 4. Component testing: experimental setups of (a) coil, (b) ferrites, and (c) shielding with heat sinks (dimensions in mm).
Figure 4. Component testing: experimental setups of (a) coil, (b) ferrites, and (c) shielding with heat sinks (dimensions in mm).
Eng 04 00049 g004
Figure 5. Experimental setup of module testing.
Figure 5. Experimental setup of module testing.
Eng 04 00049 g005
Figure 6. Component-level simulation and experimental validation of the electromagnetic–thermal simulation models’ steady-state behavior. The thermographic images were each taken for comparison after the maximum operation time.
Figure 6. Component-level simulation and experimental validation of the electromagnetic–thermal simulation models’ steady-state behavior. The thermographic images were each taken for comparison after the maximum operation time.
Eng 04 00049 g006
Figure 7. Component-level simulation and experimental validation of the electromagnetic–thermal simulation models’ transient behavior.
Figure 7. Component-level simulation and experimental validation of the electromagnetic–thermal simulation models’ transient behavior.
Eng 04 00049 g007
Figure 8. Module-level simulation and experimental validation of the electromagnetic–thermal simulation models’ transient and steady-state behavior. The thermographic images were taken for comparison after the maximum operation time.
Figure 8. Module-level simulation and experimental validation of the electromagnetic–thermal simulation models’ transient and steady-state behavior. The thermographic images were taken for comparison after the maximum operation time.
Eng 04 00049 g008
Figure 9. Comparison of the numerically calculated power losses up to an operation time of 120 min: (a) coil, (b) BFM8 ferrites, and (c) shielding with heat sinks.
Figure 9. Comparison of the numerically calculated power losses up to an operation time of 120 min: (a) coil, (b) BFM8 ferrites, and (c) shielding with heat sinks.
Eng 04 00049 g009
Table 1. Overview of existing research on coupled simulation approaches for predicting the electromagnetic–thermal behavior of wireless power transfer systems.
Table 1. Overview of existing research on coupled simulation approaches for predicting the electromagnetic–thermal behavior of wireless power transfer systems.
Coupling ApproachLevelComponentAnalyticsExperimentSimulationReference
One-wayModuleGPM xx[5]
One-wayModuleCPMxxx[4]
One-waySystemCPMxxx[6]
One-waySystemGPMxxx[6]
Two-wayModuleCPMxxx[7]
Two-wayModuleGPMxxx[8]
Two-waySystemCPMxxx[7]
Two-waySystemCPMxxx[9]
Two-waySystemCPM x[10]
Two-waySystemGPMx x[3]
x: Information included.
Table 2. Material properties of the used CPM components.
Table 2. Material properties of the used CPM components.
DensityHeat
Capacity
Thermal
Conductivity
Electrical
Conductivity
Relative
Permeability
ρ c p λ σ μ r
( g   c m 3 )( J   k g 1   K 1 )( W   m 1   K 1 )( S   m 1   )(-)
Cu-ETP [20]8.9238538157 × 10 6 1
BFM8 [12]4.87504.20*
3F36 [13]4.87504.202400
IG-F 71 [11]0.0752000.0301
CW 2243-2L [14]1.5818000.761301
Tepex dynalite
102-RG600/47%
1.3512800.2801
Silicone [17]3.2710601
AlMg3 [21]2.6690011723 × 10 6 1
* Calculated from the materials’ B-H curve.
Table 3. Steinmetz coefficients of the used ferrite materials.
Table 3. Steinmetz coefficients of the used ferrite materials.
Ferrite Material C m xy C t 2 C t 1 C t 0
3F36 6.83 × 10 2 1.443.27 8.40 × 10 5 1.10 × 10 2 1.23
BFM8 7.72 × 10 2 1.722.72 1.74 × 10 4 2.93 × 10 2 4.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zimmer, S.; Helwig, M.; Winkler, A.; Modler, N. Comparison of One-Way and Two-Way Coupled Simulation for Thermal Investigation of Vehicular Wireless Power Transfer Modules. Eng 2023, 4, 816-828. https://doi.org/10.3390/eng4010049

AMA Style

Zimmer S, Helwig M, Winkler A, Modler N. Comparison of One-Way and Two-Way Coupled Simulation for Thermal Investigation of Vehicular Wireless Power Transfer Modules. Eng. 2023; 4(1):816-828. https://doi.org/10.3390/eng4010049

Chicago/Turabian Style

Zimmer, Steve, Martin Helwig, Anja Winkler, and Niels Modler. 2023. "Comparison of One-Way and Two-Way Coupled Simulation for Thermal Investigation of Vehicular Wireless Power Transfer Modules" Eng 4, no. 1: 816-828. https://doi.org/10.3390/eng4010049

APA Style

Zimmer, S., Helwig, M., Winkler, A., & Modler, N. (2023). Comparison of One-Way and Two-Way Coupled Simulation for Thermal Investigation of Vehicular Wireless Power Transfer Modules. Eng, 4(1), 816-828. https://doi.org/10.3390/eng4010049

Article Metrics

Back to TopTop