Home Physical Sciences 3-D Electromagnetic field analysis of wireless power transfer system using K computer
Article Open Access

3-D Electromagnetic field analysis of wireless power transfer system using K computer

  • , EMAIL logo , , , and
Published/Copyright: May 30, 2018

Abstract

We analyze the electromagnetic field of a wireless power transfer system using the 3-D parallel finite element method on K computer, which is a super computer in Japan. It is clarified that the electromagnetic field of the wireless power transfer system can be analyzed in a practical time using the parallel computation on K computer, moreover, the accuracy of the loss calculation becomes better as the mesh division of the shield becomes fine.

1 Introduction

The wireless power transfer system is an open circuit and the magnetic leakage flux affects far distance. Thus, it is important to investigate the distant magnetic field of wireless power transfer systems because it is necessary to consider the influence of the magnetic flux on the human body.

In a previous study, we have calculated the distant magnetic field of a wireless power transfer system using the 3-D finite element method [1].

In order to analyze the distant magnetic field and the eddy current loss of the wireless power transfer system accurately, it is necessary to divide the mesh of the wireless power transmitter finely while analyzing the wide area. The 3-D parallel finite element method [2] is effective for such an analysis with long calculation time.

Thus, we analyze a wireless power transfer system in a wide area using the 3-D parallel finite element method on K computer, which is a super computer in Japan.

Consequently, it is clarified the distant magnetic field and the eddy current loss of the shield of the wireless power transfer system at the same time within a practical time.

2 Analysis method

2.1 Fundamental equations of magnetic field

When the magnetic flux varies sinusoidally, the fundamental equations of the electromagnetic field can be written using the magnetic vector potential A and the electric scalar potential φ as follows [3]:

rot(νrotA˙)=J˙0+J˙e(1)
J˙e=σ(jωA˙+gradϕ˙)(2)
divJ˙e=0(3)

where ν is the reluctivity, 0 is the exciting current density, e is the eddy current density, σ is the electrical conductivity, j is the imaginary unit and ω is the angular frequency. The dot over the letter indicates a complex function.

The analysis has been done as a linear problem.

2.2 Electrical loss calculation

The eddy current loss Wed in the aluminum plates of the shield is given as follows:

Wed=12σVc|J˙e|2dv(4)

where Vc is the region of the conductor with the eddy current.

2.3 Parallel computing using DDM and MPI

In this study, the domain decomposition method (DDM) is adopted for the parallel computing. Using the DDM, the analysis domain is divided into multiple subdomains. The parallel computing is performed in distributed memory type parallel computer using message passing interface (MPI) for data communication [4].

2.4 Speed-up and parallel efficiency

The speed-up and parallel efficiency, which represents the performance of parallel computation, are calculated using the parallel number n respectively as follows [5]:

Speedup(ncore)=Elapsedtime(1core)Elapsedtime(ncore)(5)
ParallelEfficiency=Speedup(ncore)n×100(%)(6)

3 Analyzed model and conditions

Figure 1 shows the analyzed model of a wireless power transfer system. The wireless power transmitter shown in Figure 1(a) is consisted of a pair of the ferrite core and the coil. The aluminum plates are located above and below the transmitter as the electromagnetic shield. The iron plate is located in the analyzed area, because it is on the ground in the experiment. The analyzed area (120m × 240m × 240m) shown in Figure 1(d) is a half of the whole region, because of the symmetry in the y-axis direction. The magnetic field at the “measured point” shown in Figures 1(b)and 1(c), is measured in the experiment.

Figure 1 Analyzed model of wireless power transfer system
Figure 1

Analyzed model of wireless power transfer system

Table 1 shows the computing environment. The K computer consists of 8 cores per node, and the main memory has 16 GB per node. In this study, two cores per node are used to calculate the electromagnetic field.

Table 1

Computing environment

ParallelK computer
computer
CPUSPARC64TM VIIIfx 128Gflops (2.0 GHz)
Cores8 cores/node
Memory16 GB/node

Table 2 shows the analysis condition. The analysis has been done using our own program.

Table 2

Analysis condition

Primary coilNumber of turns29
Current (Arms)9.1
Frequency (kHz)92.3
Phase (∘)0.0
Secondary coilNumber of turns30
Current (Arms)9.8
Frequency (kHz)92.3
Phase (∘)−35.0
Ferrite coreRelative permeability3,000
Aluminum plateRelative permeability1.0
Conductivity (S/m)3.8×107

4 Results and discussion

4.1 Elapsed time, speed-up and parallel efficiency

In this section, in order to investigate the characteristics of parallel efficiency, the magnetic field analysis was performed with various number of divisions of the analysis region.

Table 3 shows the discretization data. In order to investigate the performance of parallel computation, the analysis is performed under the condition that the number of parallel is 600 to 9600. The number of elements of the finite element model is approximately 340-380 million. The number of elements, nodes, and edges of the analysis model increases as the number of parallel increases, because the overlapping elements [4] increase.

Table 3

Discretization data (tetrahedral element)

Number of parallel60012002400
Number of elements351,998,312356,076,245361,138,171
Number of nodes61,893,53663,727,25565,987,137
Number of edges420,245,083428,374,801438,384,168
Number of unknown variables442,377,124451,063,881461,831,437
Number of parallel480072009600
Number of elements368,822,433374,349,636378,617,471
Number of nodes69,451,74071,959,15373,918,884
Number of edges453,641,230464,630,082473,164,131
Number of unknown variables478,262,848490,140,622499,342,885

Figure 2 shows the elapsed time. The calculation time is approximately 400 minutes when number of parallel is 600 and that is approximately 30 minutes when number of parallel is 9600.

Figure 2 Elapsed time
Figure 2

Elapsed time

Figure 3 shows the speed-up. The calculation time of 1 division is calculated from the elapsed times of 600 and 9600 divisions by Amdahl’s law. From Figure 3, we can see that the speed-up close to the ideal can be obtained when number of parallel is 2400 or less. The speed-up falls below the ideal as the number of parallel increases.

Figure 3 Speed-up
Figure 3

Speed-up

Figure 4 shows the parallel efficiency. The parallel efficiency is over 90% when number of parallel is 1200, and that is over 70% even when number of parallel is 9600.

Figure 4 Parallel efficiency
Figure 4

Parallel efficiency

From these results, we confirmed that using parallel computation by K computer, we can analyze the electromagnetic field of the wireless power transfer system in a wide area within a practical time. Moreover, it is confirmed that the electromagnetic field can be calculated with high efficiency even if the number of parallel is increased.

4.2 Mesh division of aluminum plate and analysis accuracy

In this section, we investigated the effects of mesh division of the aluminum plate on the eddy current loss and the distant magnetic field strength. From the results of the previous section, the number of parallel is set to 2400 for the efficient analysis.

Figure 5 shows the aluminum plate and its surrounding meshes. Each figure is an enlarged view of (i) primary side- and (ii) secondary side-aluminum plates shown in Figure 1(a). In order to investigate the effect of the fineness of the aluminum plate mesh on the eddy current loss and the distant magnetic field strength, the analysis is performed under the condition that the number of mesh divisions of the secondary side aluminum plate is 1 to 16. The primary side aluminum plate mesh is divided so as to have the same aspect ratio as the secondary side aluminum plate mesh.

Figure 5 Enlarged view of the aluminum plate and its surrounding meshes
Figure 5

Enlarged view of the aluminum plate and its surrounding meshes

Figure 6 shows the distributions of eddy current loss in the aluminum plate. We can see that the eddy current loss at secondary side is larger than that at primary side. Moreover, we can also see that the eddy current loss is concentrated on the surface of the aluminum plate by dividing the mesh finely.

Figure 6 Distributions of eddy current loss of aluminum plate
Figure 6

Distributions of eddy current loss of aluminum plate

Figure 7 shows the distant magnetic field strength at “measured point” shown in Figure 1 and the eddy current loss of the aluminum plate. The distant magnetic field strength is almost the same regardless of the mesh division of aluminum plate. However, the eddy current loss of the aluminum plate increases as the mesh becomes fine, and finally that converges to a certain value.

Figure 7 Distant magnetic field strength and eddy current loss
Figure 7

Distant magnetic field strength and eddy current loss

Table 4 shows the discretization data and elapsed time. The number of elements and calculation time also becomes to increase as the mesh division of the aluminum plate becomes finer.

Table 4

Discretization data and elapsed time

Number of mesh divisions81216
Number of parallel2,400
Number of elements330,484,691361,138,171361,979,701
Number of nodes60,616,45565,987,13766,458,832
Number of edges401,907,573438,384,168440,406,975
Number of unknown variables417,892,720461,831,437460,144,958
MRTR convergence criterion5.0×10−4
Number of MRTR Iterations18,19764,89877,818
Elapsed time (hours)0.51.82.2

5 Conclusion

In this paper, we analyzed a wireless power transfer system in a wide area using the 3-D parallel finite element method on K computer. Consequently, it is clarified that using parallel computation by K computer, we can analyze the electromagnetic field of the wireless power transfer system in a wide area within a practical time and efficiently calculate even high number of parallel. Moreover, the distant magnetic field strength hardly changes when making the mesh of the aluminum plate finer, however, the eddy current loss of the aluminum plate increases and finally that converges to a certain value.

Acknowledgement

This research used computational resources of the K computer provided by the RIKEN Advanced Institute for Computational Science through the HPCI System Research project (Project ID:hp160188).

References

[1] Kawase Y., Yamaguchi T., Tsukada S., Ota T., Yamamoto T., Distant Magnetic Field Analysis of Wireless Power Transfer System, Proc. of the XXIIth International Conference on Electrical Machines, 2016, 752-756.10.1109/ICELMACH.2016.7732610Search in Google Scholar

[2] Nakano T., Kawase Y., Yamaguchi T., Nakamura M., Parallel Computing of 3-D Eddy-Current Analysis with A-φ Method for Rotating Machines, IEEE Trans. on Magnetics., 2012, 48,2, 975-978.10.1109/TMAG.2011.2172678Search in Google Scholar

[3] Kawase Y., Mori T., Ota T., Magnetic Field Analysis of Coupling Transformers for Electric Vehicle Using 3-D Finite Element Method, IEEE Trans. on Magnetics, 1998, 34, 5, 3186-3189.10.1109/20.717747Search in Google Scholar

[4] Kawase Y., Yamaguchi T., Tsuji T., Tanaka K., Minoshima N., Hattori T., Magnetic Field Analysis of Matrix-Rotor Induction Motor Using Parallel Computing System, IEEE Transaction on Magnetics, 2011, 47, 5, 1062-1065.10.1109/TMAG.2010.2098854Search in Google Scholar

[5] Nakano T., Kawase Y., Yamaguchi T., Nomura S., Parallel Computing of 2-D Finite Element Method for Rotating Machines, Journal of the Japan Society of Applied Electromagnetics and Mechanics, 2011, 19, S157-162.Search in Google Scholar

Received: 2017-11-01
Accepted: 2018-01-16
Published Online: 2018-05-30

© 2018 Yoshihiro Kawase et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Articles in the same Issue

  1. Regular Articles
  2. A modified Fermi-Walker derivative for inextensible flows of binormal spherical image
  3. Algebraic aspects of evolution partial differential equation arising in the study of constant elasticity of variance model from financial mathematics
  4. Three-dimensional atom localization via probe absorption in a cascade four-level atomic system
  5. Determination of the energy transitions and half-lives of Rubidium nuclei
  6. Three phase heat and mass transfer model for unsaturated soil freezing process: Part 1 - model development
  7. Three phase heat and mass transfer model for unsaturated soil freezing process: Part 2 - model validation
  8. Mathematical model for thermal and entropy analysis of thermal solar collectors by using Maxwell nanofluids with slip conditions, thermal radiation and variable thermal conductivity
  9. Constructing analytic solutions on the Tricomi equation
  10. Feynman diagrams and rooted maps
  11. New type of chaos synchronization in discrete-time systems: the F-M synchronization
  12. Unsteady flow of fractional Oldroyd-B fluids through rotating annulus
  13. A note on the uniqueness of 2D elastostatic problems formulated by different types of potential functions
  14. On the conservation laws and solutions of a (2+1) dimensional KdV-mKdV equation of mathematical physics
  15. Computational methods and traveling wave solutions for the fourth-order nonlinear Ablowitz-Kaup-Newell-Segur water wave dynamical equation via two methods and its applications
  16. Siewert solutions of transcendental equations, generalized Lambert functions and physical applications
  17. Numerical solution of mixed convection flow of an MHD Jeffery fluid over an exponentially stretching sheet in the presence of thermal radiation and chemical reaction
  18. A new three-dimensional chaotic flow with one stable equilibrium: dynamical properties and complexity analysis
  19. Dynamics of a dry-rebounding drop: observations, simulations, and modeling
  20. Modeling the initial mechanical response and yielding behavior of gelled crude oil
  21. Lie symmetry analysis and conservation laws for the time fractional simplified modified Kawahara equation
  22. Solitary wave solutions of two KdV-type equations
  23. Applying industrial tomography to control and optimization flow systems
  24. Reconstructing time series into a complex network to assess the evolution dynamics of the correlations among energy prices
  25. An optimal solution for software testing case generation based on particle swarm optimization
  26. Optimal system, nonlinear self-adjointness and conservation laws for generalized shallow water wave equation
  27. Alternative methods for solving nonlinear two-point boundary value problems
  28. Global model simulation of OH production in pulsed-DC atmospheric pressure helium-air plasma jets
  29. Experimental investigation on optical vortex tweezers for microbubble trapping
  30. Joint measurements of optical parameters by irradiance scintillation and angle-of-arrival fluctuations
  31. M-polynomials and topological indices of hex-derived networks
  32. Generalized convergence analysis of the fractional order systems
  33. Porous flow characteristics of solution-gas drive in tight oil reservoirs
  34. Complementary wave solutions for the long-short wave resonance model via the extended trial equation method and the generalized Kudryashov method
  35. A Note on Koide’s Doubly Special Parametrization of Quark Masses
  36. On right-angled spherical Artin monoid of type Dn
  37. Gas flow regimes judgement in nanoporous media by digital core analysis
  38. 4 + n-dimensional water and waves on four and eleven-dimensional manifolds
  39. Stabilization and Analytic Approximate Solutions of an Optimal Control Problem
  40. On the equations of electrodynamics in a flat or curved spacetime and a possible interaction energy
  41. New prediction method for transient productivity of fractured five-spot patterns in low permeability reservoirs at high water cut stages
  42. The collinear equilibrium points in the restricted three body problem with triaxial primaries
  43. Detection of the damage threshold of fused silica components and morphologies of repaired damage sites based on the beam deflection method
  44. On the bivariate spectral quasi-linearization method for solving the two-dimensional Bratu problem
  45. Ion acoustic quasi-soliton in an electron-positron-ion plasma with superthermal electrons and positrons
  46. Analysis of projectile motion in view of conformable derivative
  47. Computing multiple ABC index and multiple GA index of some grid graphs
  48. Terahertz pulse imaging: A novel denoising method by combing the ant colony algorithm with the compressive sensing
  49. Characteristics of microscopic pore-throat structure of tight oil reservoirs in Sichuan Basin measured by rate-controlled mercury injection
  50. An activity window model for social interaction structure on Twitter
  51. Transient thermal regime trough the constitutive matrix applied to asynchronous electrical machine using the cell method
  52. On the zagreb polynomials of benzenoid systems
  53. Integrability analysis of the partial differential equation describing the classical bond-pricing model of mathematical finance
  54. The Greek parameters of a continuous arithmetic Asian option pricing model via Laplace Adomian decomposition method
  55. Quantifying the global solar radiation received in Pietermaritzburg, KwaZulu-Natal to motivate the consumption of solar technologies
  56. Sturm-Liouville difference equations having Bessel and hydrogen atom potential type
  57. Study on the response characteristics of oil wells after deep profile control in low permeability fractured reservoirs
  58. Depiction and analysis of a modified theta shaped double negative metamaterial for satellite application
  59. An attempt to geometrize electromagnetism
  60. Structure of traveling wave solutions for some nonlinear models via modified mathematical method
  61. Thermo-convective instability in a rotating ferromagnetic fluid layer with temperature modulation
  62. Construction of new solitary wave solutions of generalized Zakharov-Kuznetsov-Benjamin-Bona-Mahony and simplified modified form of Camassa-Holm equations
  63. Effect of magnetic field and heat source on Upper-convected-maxwell fluid in a porous channel
  64. Physical cues of biomaterials guide stem cell fate of differentiation: The effect of elasticity of cell culture biomaterials
  65. Shooting method analysis in wire coating withdrawing from a bath of Oldroyd 8-constant fluid with temperature dependent viscosity
  66. Rank correlation between centrality metrics in complex networks: an empirical study
  67. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering
  68. Modeling of electric and heat processes in spot resistance welding of cross-wire steel bars
  69. Dynamic characteristics of triaxial active control magnetic bearing with asymmetric structure
  70. Design optimization of an axial-field eddy-current magnetic coupling based on magneto-thermal analytical model
  71. Thermal constitutive matrix applied to asynchronous electrical machine using the cell method
  72. Temperature distribution around thin electroconductive layers created on composite textile substrates
  73. Model of the multipolar engine with decreased cogging torque by asymmetrical distribution of the magnets
  74. Analysis of spatial thermal field in a magnetic bearing
  75. Use of the mathematical model of the ignition system to analyze the spark discharge, including the destruction of spark plug electrodes
  76. Assessment of short/long term electric field strength measurements for a pilot district
  77. Simulation study and experimental results for detection and classification of the transient capacitor inrush current using discrete wavelet transform and artificial intelligence
  78. Magnetic transmission gear finite element simulation with iron pole hysteresis
  79. Pulsed excitation terahertz tomography – multiparametric approach
  80. Low and high frequency model of three phase transformer by frequency response analysis measurement
  81. Multivariable polynomial fitting of controlled single-phase nonlinear load of input current total harmonic distortion
  82. Optimal design of a for middle-low-speed maglev trains
  83. Eddy current modeling in linear and nonlinear multifilamentary composite materials
  84. The visual attention saliency map for movie retrospection
  85. AC/DC current ratio in a current superimposition variable flux reluctance machine
  86. Influence of material uncertainties on the RLC parameters of wound inductors modeled using the finite element method
  87. Cogging force reduction in linear tubular flux switching permanent-magnet machines
  88. Modeling hysteresis curves of La(FeCoSi)13 compound near the transition point with the GRUCAD model
  89. Electro-magneto-hydrodynamic lubrication
  90. 3-D Electromagnetic field analysis of wireless power transfer system using K computer
  91. Simplified simulation technique of rotating, induction heated, calender rolls for study of temperature field control
  92. Design, fabrication and testing of electroadhesive interdigital electrodes
  93. A method to reduce partial discharges in motor windings fed by PWM inverter
  94. Reluctance network lumped mechanical & thermal models for the modeling and predesign of concentrated flux synchronous machine
  95. Special Issue Applications of Nonlinear Dynamics
  96. Study on dynamic characteristics of silo-stock-foundation interaction system under seismic load
  97. Microblog topic evolution computing based on LDA algorithm
  98. Modeling the creep damage effect on the creep crack growth behavior of rotor steel
  99. Neighborhood condition for all fractional (g, f, n′, m)-critical deleted graphs
  100. Chinese open information extraction based on DBMCSS in the field of national information resources
  101. 10.1515/phys-2018-0079
  102. CPW-fed circularly-polarized antenna array with high front-to-back ratio and low-profile
  103. Intelligent Monitoring Network Construction based on the utilization of the Internet of things (IoT) in the Metallurgical Coking Process
  104. Temperature detection technology of power equipment based on Fiber Bragg Grating
  105. Research on a rotational speed control strategy of the mandrel in a rotary steering system
  106. Dynamic load balancing algorithm for large data flow in distributed complex networks
  107. Super-structured photonic crystal fiber Bragg grating biosensor image model based on sparse matrix
  108. Fractal-based techniques for physiological time series: An updated approach
  109. Analysis of the Imaging Characteristics of the KB and KBA X-ray Microscopes at Non-coaxial Grazing Incidence
  110. Application of modified culture Kalman filter in bearing fault diagnosis
  111. Exact solutions and conservation laws for the modified equal width-Burgers equation
  112. On topological properties of block shift and hierarchical hypercube networks
  113. Elastic properties and plane acoustic velocity of cubic Sr2CaMoO6 and Sr2CaWO6 from first-principles calculations
  114. A note on the transmission feasibility problem in networks
  115. Ontology learning algorithm using weak functions
  116. Diagnosis of the power frequency vacuum arc shape based on 2D-PIV
  117. Parametric simulation analysis and reliability of escalator truss
  118. A new algorithm for real economy benefit evaluation based on big data analysis
  119. Synergy analysis of agricultural economic cycle fluctuation based on ant colony algorithm
  120. Multi-level encryption algorithm for user-related information across social networks
  121. Multi-target tracking algorithm in intelligent transportation based on wireless sensor network
  122. Fast recognition method of moving video images based on BP neural networks
  123. Compressed sensing image restoration algorithm based on improved SURF operator
  124. Design of load optimal control algorithm for smart grid based on demand response in different scenarios
  125. Face recognition method based on GA-BP neural network algorithm
  126. Optimal path selection algorithm for mobile beacons in sensor network under non-dense distribution
  127. Localization and recognition algorithm for fuzzy anomaly data in big data networks
  128. Urban road traffic flow control under incidental congestion as a function of accident duration
  129. Optimization design of reconfiguration algorithm for high voltage power distribution network based on ant colony algorithm
  130. Feasibility simulation of aseismic structure design for long-span bridges
  131. Construction of renewable energy supply chain model based on LCA
  132. The tribological properties study of carbon fabric/ epoxy composites reinforced by nano-TiO2 and MWNTs
  133. A text-Image feature mapping algorithm based on transfer learning
  134. Fast recognition algorithm for static traffic sign information
  135. Topical Issue: Clean Energy: Materials, Processes and Energy Generation
  136. An investigation of the melting process of RT-35 filled circular thermal energy storage system
  137. Numerical analysis on the dynamic response of a plate-and-frame membrane humidifier for PEMFC vehicles under various operating conditions
  138. Energy converting layers for thin-film flexible photovoltaic structures
  139. Effect of convection heat transfer on thermal energy storage unit
Downloaded on 1.4.2026 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2018-0044/html
Scroll to top button