Conference Agenda

Overview and details of the sessions of this conference. Please select a date or location to show only sessions at that day or location. Please select a single session for detailed view (with abstracts and downloads if available).

 
 
Session Overview
Session
PB-A1: Numerical techniques 2
Time:
Wednesday, 24/May/2023:
1:50pm - 3:20pm

Session Chair: Dr. Zsolt Badics, Tensor Research, LLC, United States of America

Presentations
ID: 314 / PB-A1: 1
Topics: Novel Computational Methods for Machines and Devices
Keywords: Inductors, Equivalent circuit, Nonlinear circuit, Curve fitting

Loss Evaluation of Magnetic Devices using Cauer Equivalent Circuit Identified from Input Impedance

Qiao Liu1, Yuji Uehara2, Naoya Terauchi3, Hajime Igarashi1

1Hokkaido University, Japan; 2Magnetic Device Laboratory Ltd.; 3TAIYO YUDEN CO., LTD.

This article proposed a method for analyzing the inductor based on a Cauer equivalent circuit. The circuit parameters are identified from the input impedance of an inductor. Based on the synthesized Cauer equivalent circuit, it is possible to perform time-domain analysis to evaluate the eddy current losses for arbitrary driving voltage waveforms. Moreover, the magnetic saturation of the core can be considered in the Cauer circuit. It is shown that the copper loss and current response in time domain simulated from synthesized Cauer circuit agree well with the results calculated from time-domain finite element analysis.

PB-A1-1-314.pdf


ID: 191 / PB-A1: 2
Topics: Numerical Techniques, Novel Computational Methods for Machines and Devices
Keywords: Hybrid analytical model (HAM), Magnetic saturation, AC copper loss, Hairpin winding

Hybrid Analytical Model for AC Copper Loss Computation of Hairpin Winding

Yunlu Du1, Yunkai Huang1, Baocheng Guo2, Fei Peng1, Yu Yao1, Jianning Dong3

1Southeast University, China, People's Republic of; 2Nanjing Normal University, China, People's Republic of; 3Delft University of Technology, Netherlands

This paper proposes a hybrid analytical model (HAM) for predicting AC copper losses in hairpin windings. The HAM combines the nonlinear semi-analytical model (SAM) and the one-dimensional (1-D) AC copper loss computation model to solve the inaccurate boundary solutions of the 1-D AC copper loss computation model when magnetic saturation occurs in the stator teeth. This HAM is suitable for a variety of hairpin winding configurations. The finite element (FE) model validates the effectiveness of the HAM. The primary contribution of this research is to present a quick and simple approach for accurately predicting AC copper losses in hairpin windings.

PB-A1-2-191.pdf


ID: 167 / PB-A1: 3
Topics: Numerical Techniques, Novel Computational Methods for Machines and Devices, AI and Machine Learning Technologies
Keywords: Cloud computing, data-driven modeling, finite element analysis, induction motors, machine learning

Robust Development of Induction Motor Active Learning Surrogates

Janne Keränen1, Mikko Tahkola1, Peter Råback2, Álvaro González2, Victor Mukherjee3, Jenni Pippuri-Mäkeläinen1

1VTT Technical Research Centre of Finland, Finland; 2CSC–IT Center for Science, Finland; 3ABB Motors and Generators, Technology Center, Finland

A robust open-source workflow is developed for finite element (FE) data generation for active learning (AL) based surrogate modelling. Machine learning (ML) and FE simulations work iteratively together, ML in a local computer and open-source FE in a containerised cloud. We develop a surrogate model aiming for a lightweight induction machine steady-state torque and total loss prediction as a function of motor frequency, voltage, and slip. Results show how active learning improves the surrogate accuracy, especially in regions where the error was originally the largest. Special attention is given to making the cloud FE solution procedure as robust and fast as possible, e.g., by a special convergence criterion, reliable parallel computation and variable timestep length.

PB-A1-3-167.pdf


ID: 569 / PB-A1: 4
Topics: Numerical Techniques
Keywords: Finite element analysis, Eddy currents, Accelerated computing

(This presentation was moved to PO1)

Conference Organizer

Compumag2023 secretariat

The paper "A Tensor-Analysis-Based Approach to Accelerate 3D Eddy-Current

Finite Element Methods" was moved to PO1.



ID: 136 / PB-A1: 5
Topics: Static and Quasi-Static Fields, Numerical Techniques
Keywords: Domain decomposition, finite-element method, magnetic hysteresis, parallel computing, parallel-in-time method

Parallel-in-Space-and-Time Finite-Element Method for Time-Periodic Magnetic Field Problems with Hysteresis

Yasuhito Takahashi1, Koji Fujiwara1, Takeshi Iwashita2

1Doshisha University, Japan; 2Hokkaido University, Japan

This study investigates the effectiveness of a parallel-in-space-and-time finite-element method (PinSTFEM) in electric machine analyses taking account of magnetic hysteresis. The PinSTFEM is based on a combination of the domain decomposition method and parallel time-periodic explicit-error-correction (PTP-EEC) method. In addition, a magnetization history in each element is appropriately corrected based on time-series data of flux density vectors within the framework of the PTP-EEC method. The parallel performance of the developed method is verified in the 3-dimensional magnetic field analysis of a ring core specimen and a practical interior permanent magnet synchronous motor.

PB-A1-5-136.pdf


ID: 453 / PB-A1: 6
Topics: Numerical Techniques
Keywords: CFIE, regularised CFIE, Calderon preconditioner

The Galerkin method for a regularised combined field integral equation without dual basis functions

Kazuki Niino, Shunpei Yamamoto

Kyoto University, Japan

We propose a formulation of the regularised combined field integral equation (regularised CFIE) and its discretisation only with the Rao-Wilton-Glisson (RWG) method.

The CFIE is a formulation of integral equations, which avoids the so-called ficticious frequencies of integral equations.

The most typical CFIE, which is a linear combination of the electric field integral equation (EFIE) and magnetic field integral equation (MFIE),

is known to be ill-conditioned in terms of computational time.

The regularised CFIE is another formulation of the CFIE to solve the ill-conditioning of the original CFIE

by applying a regularising operator to the part of the EFIE.

In many previous studies the regularising operator is determined based on the Calderon preconditioning.

This regularising operator however takes much more computatonal time than the standard CFIE

since discretising the EFIE with the Calderon preconditioner requires the dual basis function.

In this article we propose a formulation of the regularised CFIE, which can be discretised with the Galerkin method without the dual basis function.

PB-A1-6-453.pdf


ID: 222 / PB-A1: 7
Topics: Numerical Techniques
Keywords: UHVDC transmission lines, UPFEM, ion flow field, rain environment.

Numerical Simulation Method for Ion Flow Field of UHVDC Transmission Lines in Light Rain Environment

Guohua Yue, Zhiye Du, Liancheng Xiu, Hongwei Cai

Wuhan University, China, People's Republic of

The UHVDC transmission line covers a wide area, and the climate along the line is complex and changeable, so rainfall is inevitable. In rainy days, the parameters such as corona onset field, ion mobility and ion recombination coefficient of the conductor will change compared with those in sunny days. In addition, the analysis of the charge capacity and force of suspended droplets, falling raindrops and other particles should be added. Taking the above factors into account and based on the upstream finite element method (UPFEM), this paper proposes a numerical simulation method for the ion flow field of UHVDC transmission lines in light rain environment, and validates it under laboratory conditions. The results show that the total electric field and ion current density have increased to a certain extent compared with sunny weather. The results calculated by this method is close to the measured data in the ±800 kV line in China.

PB-A1-7-222.pdf


ID: 378 / PB-A1: 8
Topics: Numerical Techniques, Multi-Scale Modelling and Homogenization
Keywords: Multi-scale finite element method, numerical simulation, eddy currents

Multi-Scale Finite Element Method for 3D Eddy Current Problems Based on Distance Regularized Level Set Function

Xinyu Ma1, Nana Duan1, Weijie Xu2, Ruting Tang1, Shuhong Wang1

1Xi'an Jiaotong University, China, People's Republic of; 2State Grid Shaanxi Electric Power Research Institute, China, People's Republic of

3D eddy current calculation plays an important role in the design of electrical equipment. However, due to the special structure of some equipment, numerical calculation often faces great difficulties. For example, in the laminate core for transformers, the overall dimensions of the core can be about 104 times the thickness of the laminate. A large number of elements and nodes will be required using the finite element method in analyzing such equipment with laminates. The computational cost would be extremely high, which cannot be afforded. This paper introduces a multi-scale finite element method (MSFEM) that can effectively reduce computational costs. Using a special interpolation built with an enrichment function, the elements are not restricted by media or geometry. The enrichment function is presented through distance regularization level set evolution, which enables MSFEM to cope with complex 3D model. Taking Team Workshop problem 7 as a numerical example, the accuracy of this method and the saving of computational costs have been proved.

PB-A1-8-378.pdf


ID: 146 / PB-A1: 9
Topics: Numerical Techniques
Keywords: 3D UPFEM, ion flow field, convergence rate, space charge density update

Fast Calculation Method of 3D Ion Flow Field Based on Adaptive Iteration Factor and Half-error Term Correction

Zhiye Du, Guohua Yue, Liancheng Xiu, Gen Li, Jingwen Huang

School of Electrical Engineering and Automation, Wuhan University, China, People's Republic of

The solution speed of 3D ion flow field is generally very slow, which seriously restricts the application of this method in UHVDC transmission line projects. In this paper, based on the 3D upstream finite element method (UPFEM), an adaptive iteration factor for charge density update is introduced, which can automatically control the iteration speed according to the change degree of electric field strength, ensuring the rapid and stable update of charge density. At the same time, to solve the problem that the convergence rate of charge density and electric field strength is inconsistent after the introduction of iteration factor, a half-error term is added to the charge density update formula. This controllable disturbance is only half of the charge density error, so it has little impact on the convergence of the charge density, but can greatly accelerate the convergence of the electric field intensity. It solves the short board effect, balances the convergence of charge density and electric field strength, and improves the overall solution speed of the algorithm. This method is used to calculate the ion flow field between a 3D coaxial cylindrical electrode. The numerical simulation results are compared with the analytical solution. Furthermore, the iterative speed of the traditional method and this method is also compared, which verifies the accuracy and effectiveness of the method proposed in this paper. Finally, the method is applied to a real ± 800kV UHVDC transmission line in China. The results show that the method is in good agreement with the test results and has strong engineering practicability.

PB-A1-9-146.pdf


ID: 383 / PB-A1: 10
Topics: Numerical Techniques
Keywords: Model order reduction (MOR), rotor eccentricity, hybrid POD-interpolation method.

Reduced Order Model of Permanent Magnet Synchronous Motor Considering Rotor Eccentricity

Ze Guo1, Shuai Yan1, Zhuoxiang Ren1,2

1Institute of Electrical Engineering, Chinese Academy of Sciences, China, People's Republic of; 2Group of Electrical and Electronic Engineering of Paris, Sorbonne Universite, CNRS, F-75005 Paris, France

In order to calculate quickly and accurately the magnetic field of the rotating machine when the rotor is eccentric, and study the influence of rotor eccentricity on its operational performance, this paper proposes a hybrid method combing the proper orthogonal decomposition (POD) and the interpolation technique to establish a behavior model of a permanent magnet synchronous motor considering rotor eccentricity. By dividing the parameter space reasonably, the problem of the torque deviation caused by rotor eccentricity submerged in the error caused by POD is avoided. The results show that the proposed model can provide accurate dynamic time-domain response and torque spectrum of the motor, and the simulation time can meet the real-time requirement. It offers an effective method not only for fast design optimization, but also for real-time monitoring and predictive fault diagnosis of the motor.

PB-A1-10-383.pdf


ID: 347 / PB-A1: 11
Topics: Numerical Techniques
Keywords: Computational electromagnetics, Finite element analysis, Forecast uncertainty, Random processes

An Efficient Uncertainty Quantification Method for Large-Scale High-Frequency Electromagnetic Field Problems

Sota Goto1, Amane Takei2, Shigeki Kaneko1, Shinobu Yoshimura1

1The University of Tokyo, Japan; 2University of Miyazaki

This paper describes the development of an integrated method of high-frequency electromagnetic wave analysis using a parallel finite element method and an efficient Uncertainty Quantification (UQ) method, i.e. the Non-Intrusive Polynomial Chaos (NIPC) method, based on the assumption that input parameters are introduced as random fields. First, taking a heat transfer problem with the NIPC method, we investigate the computational complexity of each process appeared in the NIPC method, and clarify that the bottleneck process is the pre-processing eigenvalue analysis. Next, we propose a new method to reduce the computational complexity of the eigenvalue analysis by employing a coarser finite element mesh than an original mesh used in the deterministic analysis. Here we obtain eigenfunctions on the original coarse mesh, and then convert the coarse mesh to the fine mesh. Finally, by integrating the NIPC method with the proposed computational reduction method and the parallel finite element electromagnetic field analysis method, we successfully perform the UQ analysis of large-scale high-frequency electromagnetic wave problem, which has been difficult to achieve so far.

PB-A1-11-347.pdf


ID: 563 / PB-A1: 12
Topics: Numerical Techniques
Keywords: Isogeometric Analysis, Mortar, Refinement

Isogeometric mesh refinement through Non-conforming C0 Isogeometric Mortaring

Maxime Fays, Olivier Chadebec, Brahim Ramdane

Univ. Grenoble Alpes, CNRS, Grenoble INP, France

An original mesh refinement of Isogeometric model through $C^0$ isogeometric mortaring is proposed in this paper to address patch refinement propagation occuring when handling multipatch geometries, a particularly relevent side-effect in electromagnetic problems. The efficency of the proposed method is demonstrated on an academic problem.

PB-A1-12-563.pdf


ID: 254 / PB-A1: 13
Topics: Numerical Techniques
Keywords: Convergence Characteristic, Electric Circuit, Electromagnetic Finite Element Analysis, Solid Conductor

Investigation of Convergence Characteristics of the Iterative Method in Electromagnetic Finite Element Analysis with Electric Circuit

Kazuki SEMBA1, Tomohisa SATO2, Hirokatsu KATAGIRI3, Takashi YAMADA4

1JSOL Corporation, Japan; 2JSOL Corporation, Japan; 3JSOL Corporation, Japan; 4JSOL Corporation, Japan

The effect of AC loss is non-negligible for high-speed rotating machines with their associated PWM control. The same is true for the joule loss in the coil of the induction heating systems. To estimate those losses, solid conductors are used in 3-dimentional electromagnetic finite element analysis (FEA). However, the convergence characteristics of the iterative method deteriorates in FEA with solid conductors. The A-φ method and the application of appropriate initial state for each time step make it possible to overcome this challenge for transient analysis. We present the effectiveness of these two approaches by employing several practical analysis models.

PB-A1-13-254.pdf


ID: 411 / PB-A1: 14
Topics: Numerical Techniques
Keywords: Computational electromagnetics, Newton method

Improving the Newton initial guess for circuit coupled magnetostatic problems

Riheb Cherif1,2, Frédéric Guyomarch2, Zuqi Tang2, Julien Korecki2, Yvonnick Le Menach2

1ESME Research Lab, Lille, France; 2Laboratoire d'Électrotechnique et d'Électronique de Puissance, University of Lille, Arts et Metiers Paris Tech, Centrale Lille, Junia, Lille, France

This paper aims to propose an Improved Starting Point (ISP-) Newton method applied to vector potential A formulation for circuit coupled magnetostatic problems. These problems usually are known to vary greatly during the time. This impacts the quality of the initial guess of the Newton method and thus increases significantly the computational cost. The Newton method with vector potential formulation A has been analyzed. Numerical examples show the performance of our proposed ISP-Newton method.

PB-A1-14-411.pdf


ID: 454 / PB-A1: 15
Topics: Numerical Techniques
Keywords: Reconfigurable architectures, high-performance computing (HPC), computer architectures, eddy currents

Conceptual Design of FIT Dataflow Machine for 2-D Time Domain Eddy Current Field Simulation

Chenxu Wang, Hideki Kawaguchi

Muroran Institute of Technology, Japan

The authors have been working in the development of dedicated computer based on dataflow architecture to achieve a portable, low power consumption, low cost and green high-performance computing (HPC) technology for electromagnetic fields simulation, which can be effectively used in industry applications. A dataflow architecture machine of 2-D finite integration technique (FIT) based on BiCG-Stab scheme matrix solver for magneto-static fields was proposed in the previous work. In this paper, we discuss a conceptual design of FIT dataflow machine for 2-D time domain eddy current fields simulation.

PB-A1-15-454.pdf


ID: 521 / PB-A1: 16
Topics: Optimization and Design, Numerical Techniques
Keywords: Permanent magnet machines, finite element analysis, finite element methods, magnetic circuit, wind energy

Analyzing the Performance of Radial-Flux PMSGs for Applications in Wind Turbines by using different Types of Windings

EDUARDO ORTIZ1, REYNALDO IRACHETA2

1CINVESTAV, Guadalajara; 2UNISTMO, Mexico

In this paper is presented an improved version of a general methodology for calculating the optimal solution of non-conventional winding distribution for applications in Wind Turbines based on a direct-drive permanent magnet synchronous generator. The proposed methodology is also well suited for different conventional windings in electrical rotating machines. Moreover, the aim of this paper is to compare the performance of wind turbines based on PMSGs with different types of windings such as concentric, integer and distributed. To achieve this, an exhaustive parametric study through 2D and 3D finite element analysis is carried out in order to measure variables such as the copper mass, losses, the power quality and, cogging and ripple torques. Finally, in a comparative table are described the main characteristics of each type of winding and the best performance for the purpose of using in wind turbines.

PB-A1-16-521.pdf


ID: 401 / PB-A1: 17
Topics: Mathematical Modelling and Formulations, Numerical Techniques, Multi-Physics and Coupled Problems
Keywords: Electrical circuit, finite element analysis, mechanical dynamics, motors, three mass system

Finite Element Modelling of Induction Machine Coupled with Three-Mass Mechanical System

Mehrnaz Farzam Far1, Jenni Pippuri-Mäkeläinen1, Janne Keränen1, Anouar Belahcen2, Sampo Laine3, Raine Viitala3

1VTT Technical Research Centre of Finland Ltd, Finland; 2Deptartment of Electrical Engineering and Automation, Aalto University, Finland; 3Department of Mechanical Engineering, Aalto University, Finland

Modern electric drive systems consist of an electric motor drive connected to a mechanical system. Single-mass models are commonly used in electromechanical finite-element analyses to study various aspects of these systems such as the output torque of the electric motor or the transient performance of electric drives. Nevertheless, single-mass models are inherently not sufficient for the analysis of torsional vibrations. To inspect the behavior of the electric drives connected to flexible torsional systems, higher degree-of-freedom mechanical models should be employed. This paper presents a comprehensive finite element model of an induction machine coupled with a three-mass lumped-element model. The basic background for the model and some initial time-stepping results are presented in this paper. The results of the model in the frequency domain will be presented in the full paper.

PB-A1-17-401.pdf