Spectral Element-Based Multi-Physical Modeling Framework for Axisymmetric Wireless Power Transfer Systems
Abstract
:1. Introduction
2. Materials and Methods
2.1. Benchmark System
2.2. Global Overview of the SEM
- 1
- An initialization step is performed, in which the nodes where the field variable is approximated are derived. Furthermore, the matrices required to compute integrals and derivatives are derived, which are essential components for the approximation of the governing equations. These aspects are further detailed in Section 2.3.1.
- 2
- The higher-order trial functions employed by the SEM to approximate the solution, achieve the optimal convergence rate in a [−1, 1] square domain. The approximation of the governing equations for such a domain is discussed in Section 2.3.2. In reality, such a domain is scarcely ever observed. Therefore, a mapping is applied, such that the convergence rate is maintained, while the geometrical flexibility of the method is greatly enhanced. In Section 2.3.3, the mapping functions are presented and the weak form of the governing equations under mapping is given for the electromagnetic model. In the case of the thermal model, the latter is given in Section 2.4.1, whereas the mapping functions are equivalent.
- 3
- In order to solve the linear system of equations, the problem is expressed in its equivalent matrix form, which is further detailed in Section 2.3.4 and Section 2.4.1 for the electromagnetic and thermal model, respectively.
- 4
- The boundary conditions are applied, which is the last step before the system of equations is solved. The various boundary conditions occurring in an axisymmetric WPT system are discussed in Section 2.3.5 and Section 2.4.2 for the electromagnetic and thermal model, respectively.
- 5
- The linear system of equations is solved, yielding the solution of the field variable at the nodes in the investigated domain.
- 6
- Finally, various post-processing calculations are performed, e.g., to obtain the various loss components or the average temperature. The post-processing calculations, relevant to axisymmetric WPT systems, are discussed in Section 2.3.6 and Section 2.4.3 for the electromagnetic and thermal model, respectively.
2.3. Electromagnetic Formulation
- 1
- The electrically conductive source regions, i.e., the copper foil windings situated in the primary and secondary core.
- 2
- The electrically conductive non-source regions, i.e., the aluminum housing, which surrounds the design space, and the stainless-steel shaft.
- 3
- The electrically non-conductive non-source regions, i.e., the ferrite cores and air gaps.
2.3.1. Spectral Approximation
Legendre Polynomials
Quadratures
Interpolation
Differentiation
2.3.2. Nodal Galerkin Approximation in the Reference Domain
2.3.3. Mapping Functions for Non-Squared Elements
2.3.4. Nodal Galerkin Approximation in Matrix Form
2.3.5. Boundary Conditions
Continuous Boundary Condition
Dirichlet Boundary Condition
Neumann Boundary Condition
2.3.6. Post-Processing
Magnetic Flux Density
Hysteresis Effect Losses
Joule Effect Losses
Flux-Linkage
2.4. Thermal Formulation
2.4.1. Nodal Galerkin Approximation
2.4.2. Boundary Conditions
Neumann Boundary Condition
Robin Boundary Condition
2.4.3. Post-Processing
2.5. Detailed Overview of the SEM Implementation
- 1
- 2
- 3
- 4
- The sub-matrices and sub-vector are assembled into the global matrix, , and source vector, , of (43). Furthermore, the boundary conditions are imposed according to the approach discussed in Section 2.3.5.
- 5
- The system of equations is solved to obtain the approximate solution at the nodes. In the context of this paper, a direct solver is applied. Alternatively, for problems having a high number of non-zero elements, an iterative solver can be applied. In case the problem is non-linear, a non-linear solver, e.g., the Newton–Raphson or JFNK methods, has to be included to solve the non-linear state.
- 6
- After the solution has been obtained, the post-processing computations of Section 2.3.6 are executed. However, depending on the problem, additional quantities, e.g., the electromagnetic torque, can be added.
2.6. Verification
2.6.1. SEM Model
2.6.2. FEM Model
2.6.3. Method of Comparison
3. Results
4. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
Abbreviations
2D | Two-dimensional |
3D | Three-dimensional |
ac | Alternating current |
DoF | Degree of freedom |
FEM | Finite element method |
JFNK | Jacobian-free Newton–Krylor |
LGL | Legendre–Gauss–Lobatto |
PDEs | Partial differential equations |
SEM | Spectral element method |
SIBCs | Surface impedance boundary conditions |
WPT | Wireless power transfer |
References
- Zhang, C.; Lin, D.; Hui, S.Y.R. Ball-Joint Wireless Power Transfer Systems. IEEE Trans. Power Electron. 2018, 33, 65–72. [Google Scholar] [CrossRef]
- Han, H.; Mao, Z.; Zhu, Q.; Su, M.; Hu, A.P. A 3D Wireless Charging Cylinder With Stable Rotating Magnetic Field for Multi-Load Application. IEEE Access 2019, 7, 35981–35997. [Google Scholar] [CrossRef]
- Maier, M.; Parspour, N. Operation of an Electrical Excited Synchronous Machine by Contactless Energy Transfer to the Rotor. IEEE Trans. Ind. Appl. 2018, 54, 3217–3225. [Google Scholar] [CrossRef]
- Xia, Q.; Yan, L. Application of Wireless Power Transfer Technologies and Intermittent Energy Harvesting for Wireless Sensors in Rotating Machines. Wirel. Power Transf. 2016, 3, 93–104. [Google Scholar] [CrossRef]
- Chen, X.Y.; Jin, J.X.; Zheng, L.H.; Wu, Z.H. A Rotary-Type Contactless Power Transfer System Using HTS Primary. IEEE Trans. Appl. Supercond. 2016, 26, 5501406. [Google Scholar] [CrossRef]
- Qian, W.; Zhang, X.; Fu, Y.; Lu, J.; Bai, H. Applying Normally-Off GaN HEMTs for Coreless High-Frequency Wireless Chargers. CES Trans. Electr. Mach. Syst. 2017, 1, 418–427. [Google Scholar] [CrossRef]
- Covic, G.A.; Boys, J.T. Modern Trends in Inductive Power Transfer for Transportation Applications. IEEE J. Emerg. Sel. Top. Power Electron. 2013, 1, 28–41. [Google Scholar] [CrossRef]
- van den Bossche, A.; Valchev, V.C. Inductors and Transformers for Power Electronics; Taylor & Francis Group: New York, NY, USA, 2005. [Google Scholar]
- Kulkarni, S.V.; Khaparde, S.A. Transformer Engineering; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
- Trevisan, R.; Costanzo, A. A 1-kW Contactless Energy Transfer System Based on a Rotary Transformer for Sealing Rollers. IEEE Trans. Ind. Electron. 2014, 61, 6337–6345. [Google Scholar] [CrossRef]
- Ditze, S.; Endruschat, A.; Schriefer, T.; Rosskopf, A.; Heckel, T. Inductive Power Transfer System with a Rotary Transformer for Contactless Energy Transfer on Rotating Applications. In Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS), Montreal, QC, Canada, 22–25 May 2016; pp. 1622–1625. [Google Scholar] [CrossRef]
- Zhu, X.; Lin, B.; Liu, L.; Luan, Y. Power Transfer Performance and Cutting Force Effects of Contactless Energy Transfer System for Rotary Ultrasonic Grinding. IEEE Trans. Ind. Electron. 2016, 63, 2785–2795. [Google Scholar] [CrossRef]
- Bastiaens, K.; Krop, D.C.J.; Jumayev, S.; Lomonova, E.A. Optimal Design and Comparison of High-Frequency Resonant and Non-Resonant Rotary Transformers. Energies 2020, 13, 929. [Google Scholar] [CrossRef] [Green Version]
- Raminosoa, T.; Wiles, R.H.; Wilkins, J. Novel Rotary Transformer Topology with Improved Power Transfer Capability for High-Speed Applications. IEEE Trans. Ind. Appl. 2020, 56, 277–286. [Google Scholar] [CrossRef]
- Ruviaro, M.; Runcos, F.; Sadowski, N.; Borges, I.M. Analysis and Test Results of a Brushless Doubly Fed Induction Machine with Rotary Transformer. IEEE Trans. Ind. Electron. 2012, 59, 2670–2677. [Google Scholar] [CrossRef]
- Bastiaens, K.; Krop, D.C.J.; Curti, M.; Jumayev, S.; Lomonova, E.A. Coupled Design Approach for an Integrated GaN-Based Wireless Power Transfer Rotating System. Int. J. Appl. Electromagn. Mech. 2021, 64, S51–S64. [Google Scholar] [CrossRef]
- Wang, F.; Feng, T.; Chen, X. Design and Optimisation of a Wireless Power Transfer System for Satellite Application. IET Power Electron. 2019, 12, 2586–2598. [Google Scholar] [CrossRef]
- Godbehere, J.; Hopkins, A.; Yuan, X. Design and Thermal Analysis of a Rotating Transformer. In Proceedings of the IEEE International Electric Machines Drives Conference (IEMDC), San Diego, CA, USA, 12–15 May 2019; pp. 2144–2151. [Google Scholar] [CrossRef]
- Campi, T.; Cruciani, S.; Feliziani, M. Magnetic Shielding of Wireless Power Transfer Systems. In Proceedings of the International Symposium on Electromagnetic Compatibility, Tokyo, Japan, 12–16 May 2014; pp. 422–425. [Google Scholar]
- de Gersem, H. Spectral-Element Method for High-Speed Rotating Cylinders. COMPEL 2009, 28, 730–740. [Google Scholar] [CrossRef]
- De Gersem, H. Combined Spectral-Element, Finite-Element Discretization for Magnetic-Brake Simulation. IEEE Trans. Magn. 2010, 46, 3520–3523. [Google Scholar] [CrossRef]
- Jayasekera, K.A.S.N.; Ciric, I.R. Evaluation of Surface Impedance Models for Axisymmetric Eddy-Current Fields. IEEE Trans. Magn. 2007, 43, 1991–2000. [Google Scholar] [CrossRef]
- Wagner, B.; Renhart, W.; Magele, C. Error Evaluation of Surface Impedance Boundary Conditions with Magnetic Vector Potential Formulation on a Cylindrical Test Problem. IEEE Trans. Magn. 2008, 44, 734–737. [Google Scholar] [CrossRef]
- Stein, A.L.; Kyaw, P.A.; Sullivan, C.R. Wireless Power Transfer Utilizing a High-Q Self-Resonant Structure. IEEE Trans. Power Electron. 2019, 34, 6722–6735. [Google Scholar] [CrossRef]
- Pereira, A.; Sixdenier, F.; Raulet, M.A.; Lefebvre, B.; Burais, N. Comparison Between Numerical and Analytical Methods of AC Resistance Evaluation for Medium-Frequency Transformers: Validation on a Prototype and Thermal Impact Analysis. Can. J. Electr. Comput. Eng. 2017, 40, 101–109. [Google Scholar] [CrossRef]
- Shen, J.; Tang, T.; Want, L.L. Spectral Methods Algorithms, Analysis and Applications; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
- Curti, M.; Van Beek, T.A.; Jansen, J.W.; Gysen, B.L.J.; Lomonova, E.A. General Formulation of the Magnetostatic Field and Temperature Distribution in Electrical Machines Using Spectral Element Analysis. IEEE Trans. Magn. 2018, 54, 1–9. [Google Scholar] [CrossRef]
- Trefethen, L.N. Spectral Methods in Matlab; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2000. [Google Scholar]
- Canuto, C.; Hussaini, M.Y.; Quarteroni, A.; Zang, T.A. Spectral Methods Fundamentals in Single Domains; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
- Curti, M. A Spectral Element Based Framework for Efficient Computation of Electromechanical Problems. Ph.D. Thesis, Eindhoven University of Technology, Eindhoven, The Netherlands, 2019. [Google Scholar]
- Bastiaens, K.; Curti, M.; Krop, D.C.J.; Jumayev, S.; Lomonova, E.A. Spectral Element Method Modeling of Eddy Current Losses in High-Frequency Transformers. Math. Comput. Appl. 2019, 24, 28. [Google Scholar] [CrossRef] [Green Version]
- Curti, M.; Paulides, J.J.H.; Lomonova, E.A. Magnetic Modeling of a Linear Synchronous Machine With the Spectral Element Method. IEEE Trans. Magn. 2017, 53, 8112006. [Google Scholar] [CrossRef] [Green Version]
- Zanella, R.; Nore, C.; Bouillault, F.; Cappanera, L.; Tomas, I.; Mininger, X.; Guermond, J. Study of Magnetoconvection Impact on a Coil Cooling by Ferrofluid with a Spectral/Finite-Element Method. IEEE Trans. Magn. 2018, 54, 4600104. [Google Scholar] [CrossRef]
- Regenauer-Lieb, K.; Veveakis, M.; Poulet, T.; Wellmann, F.; Karrech, A.; Liu, J.; Hauser, J.; Schrank, C.; Gaede, O.; Fusseis, F.; et al. Multiscale Coupling and Multiphysics Approaches in Earth Sciences: Applications. J. Coupled Syst. Multiscale Dyn. 2013, 1, 281–323. [Google Scholar] [CrossRef] [Green Version]
- Zhang, H.; Guo, J.; Lu, J.; Li, F.; Xu, Y.; Downar, T.J. An Assessment of Coupling Algorithms in HTR Simulator TINTE. Nuclear Sci. Eng. 2018, 190, 287–309. [Google Scholar] [CrossRef]
- Price, D.; Radaideh, M.I.; Mui, T.; Katare, M.; Kozlowski, T. Multiphysics Modeling and Validation of Spent Fuel Isotopics Using Coupled Neutronics/Thermal-Hydraulics Simulations. J. Coupled Syst. Multiscale Dyn. 2020, 2020, 2764634. [Google Scholar] [CrossRef]
- Fernández, M.Á.; Moubachir, M. A Newton Method Using Exact Jacobians for Solving Fluid–Structure Coupling. Comput. Struct. 2005, 83, 127–142. [Google Scholar] [CrossRef] [Green Version]
- Ao, X.; Liu, J.; Xia, H.; Yang, Y. A Numerical Study on the Mesoscopic Characteristics of Ti-6Al-4V by Selective Laser Melting. Materials 2022, 15, 2850. [Google Scholar] [CrossRef]
- Yan, D.; Wang, W.; Li, R.; Jiang, S.; Lu, L.; Levtsev, A.; Chen, D. Multi-Physical and Electrochemical Coupling Model for the Protonic Ceramic Fuel Cells with H+/e−/O2− Mixed Conducting Cathodes. Appl. Sci. 2022, 12, 3889. [Google Scholar] [CrossRef]
- Kopriva, D.A. Implementing Spectral Methods for Partial Differential Equations; Springer Science + Business Media B.V.: Dordrecht, The Netherlands, 2009. [Google Scholar]
- Winslow, A. Numerical Calculation of Static Magnetic Fields in an Irregular Triangular Mesh; UCR-7784; University of California: Livermore, CA, USA, 1964. [Google Scholar]
- Touzani, R.; Rappaz, J. Mathematical Models for Eddy Currents and Magnetostatics; Springer Science + Business Media B.V.: Dordrecht, The Netherlands, 2014. [Google Scholar]
- Bermúdez, A.; Gómez, D.; Salgado, P. Mathematical Models and Numerical Simulation in Electromagnetism; Springer International Publishing AG: Cham, Switzerland, 2014. [Google Scholar]
- van de Vosse, F.N.; Minev, P. Spectral Element Methods: Theory and Applications; EUT Report; Department of Mechanical Engineering, Eindhoven University of Technology: Eindhoven, The Netherlands, 1996. [Google Scholar]
- Bertotti, G. Hysteresis in Magnetism; Academic Press, Inc.: San Diego, CA, USA, 1998; ISBN 978-0-12-093270-2. [Google Scholar]
- Altair Engineering, Inc., Meylan, France. User Guide Flux 2019.1. 2018. Available online: https://www.altair.com/flux/ (accessed on 29 March 2022).
- Davis, T.A. UMFPACK User Guide. 2019. Available online: https://fossies.org/linux/SuiteSparse/UMFPACK/Doc/UMFPACK_UserGuide.pdf (accessed on 29 March 2022).
Geometrical Parameters | |||
---|---|---|---|
Parameter | Symbol | Value | Unit |
Offset radius of the geometry | 1.00 | pm | |
Outer radius of the shaft | 1.50 | mm | |
Inner radius of the core | 1.61 | mm | |
Inner radius of the coil area | 2.62 | mm | |
Inner radius of the air gap | 4.13 | mm | |
Outer radius of the coil area | 6.21 | mm | |
Outer radius of the core | 7.30 | mm | |
Inner radius of the housing | 7.50 | mm | |
Radial thickness of the housing | 1.00 | mm | |
Inner radius of the 1st turn | 2.92 | mm | |
Inner radius of the 2nd turn | 3.68 | mm | |
Inner radius of the 3rd turn | 4.95 | mm | |
Inner radius of the 4th turn | 5.74 | mm | |
Radial thickness of the air gap | 5.00 × 10 | mm | |
Base height of the housing | 5.00 × 10 | mm | |
Height of the surrounding air | 5.79 | mm | |
Base height of the core | 1.02 | mm | |
Total height of the coil area | 4.38 | mm | |
Height of the coil | 4.18 | mm | |
Thickness of the foil winding | 1.00 × 10 | mm | |
Physical Quantities | |||
Parameter | Symbol | Value | Unit |
Electrical frequency | f | 1.00 | MHz |
Electrical conductivity of the housing | 3.77 × 10 | Sm | |
Thermal conductivity of the housing | 237 | WmK | |
Electrical conductivity of the shaft | 1.45 × 10 | Sm | |
Thermal conductivity of the shaft | 16.3 | WmK | |
Relative permeability of the core | 2.30 × 10 | - | |
Thermal conductivity of the core | 4.25 | WmK | |
Electrical conductivity of the foil | 5.81 × 10 | Sm | |
Thermal conductivity of the foil | 385 | WmK | |
Thermal conductivity of the insulator | 1.00 | WmK | |
Thermal conductivity of air | 2.57 × 10 | WmK | |
Ambient temperature | 293 | K | |
Convection coefficient | 7.00 | WmK | |
Emissivity coefficient | 4.00 × 10 | - | |
Imposed current (primary side) | 4.95∠0.00 | A | |
Imposed current (secondary side) | 4.95∠90.0 | A | |
Hysteresis-loss coefficient | 10.6 | WsTm | |
Hysteresis-loss coefficient | 1.30 | - | |
Hysteresis-loss coefficient | 2.70 | - |
Quantity | Symbol | FEM | SEM | Unit | Discrepancy [%] |
---|---|---|---|---|---|
Flux-linkage primary coil | 1.271 | 1.270 | Wb-t | 5.540 × 10 | |
Flux-linkage secondary coil | 1.248 | 1.248 | Wb-t | 5.616 × 10 | |
Joule losses primary coil | 189.0 | 189.0 | mW | 2.159 × 10 | |
Joule losses secondary coil | 156.3 | 156.3 | mW | 1.160 × 10 | |
Hysteresis losses primary core | 7.072 | 7.063 | mW | 1.387 × 10 | |
Hysteresis losses secondary core | 37.87 | 37.82 | mW | 1.494 × 10 | |
Eddy current losses shaft | 55.41 | 55.39 | mW | 3.216 × 10 | |
Eddy current losses housing | 14.87 | 14.86 | mW | 7.786 × 10 | |
Temperature primary coil | 61.64 | 61.63 | C | 2.240 × 10 | |
Temperature secondary coil | 68.22 | 68.17 | C | 8.050 × 10 |
Electromagnetic | |||
---|---|---|---|
Quantity | Unit | SEM | FEM |
Number of DoF | - | 1.738 × 10 | 1.488 × 10 |
Number of non-zero elements | - | 3.999 × 10 | 2.347 × 10 |
Sparsity | % | 99.87 | 99.99 |
Thermal | |||
Quantity | Unit | SEM | FEM |
Number of DoF | - | 8.128 × 10 | 1.507 × 10 |
Number of non-zero elements | - | 1.367 × 10 | 2.350 × 10 |
Sparsity | % | 99.79 | 99.99 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Bastiaens, K.; Krop, D.C.J.; Lomonova, E.A. Spectral Element-Based Multi-Physical Modeling Framework for Axisymmetric Wireless Power Transfer Systems. Energies 2022, 15, 3145. https://doi.org/10.3390/en15093145
Bastiaens K, Krop DCJ, Lomonova EA. Spectral Element-Based Multi-Physical Modeling Framework for Axisymmetric Wireless Power Transfer Systems. Energies. 2022; 15(9):3145. https://doi.org/10.3390/en15093145
Chicago/Turabian StyleBastiaens, Koen, Dave C. J. Krop, and Elena A. Lomonova. 2022. "Spectral Element-Based Multi-Physical Modeling Framework for Axisymmetric Wireless Power Transfer Systems" Energies 15, no. 9: 3145. https://doi.org/10.3390/en15093145
APA StyleBastiaens, K., Krop, D. C. J., & Lomonova, E. A. (2022). Spectral Element-Based Multi-Physical Modeling Framework for Axisymmetric Wireless Power Transfer Systems. Energies, 15(9), 3145. https://doi.org/10.3390/en15093145