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
WP1A: Mathematical Modeling #2
Time:
Wednesday, 25/June/2025:
11:05am - 12:45pm

Session Chair: Dr. ZUQI TANG
Session Chair: Prof. Jan Sykulski
Location: Partenope

Poster Room #1

Presentations

A Novel EC Model for CFRP Defect Detection and Classification

mario versaci, giovanni angiulli, francesco carlo morabito

Mediterranea University Reggio Calabria (Italy), Italy

This study adapts a well-known Eddy Current model to detect subsurface defects in CFRP plates, where each defect is modeled as an additive variation of the electrical conductivity tensor. By translating the model into its weak form and applying an energy functional based FEM technique, we obtain 2D maps of $|\mathbf{B}|$ (stable and precise numerical solutions) corresponding to the mechanical stress maps, offering an accurate alternative, economical and effective for non-destructive inspection.

WP1A-1-126.pdf


Consistency and Convergence Analysis of Co-Simulation Techniques for Field-Circuit Coupled Problems

Michael Wiesheu1, Sebastian Schöps1, Idoia Cortes Garcia2

1TU Darmstadt, Germany; 2Eindhoven University of Technology

This work compares different techniques for co-simulation when dealing with field-circuit coupled problems. The convergence behavior of the coupled system is investigated using different time-stepping techniques and coupling methods, such as source and parameter coupling. The speed-up of the methods with regard to accuracy is explored, and challenges which can occur in the context of differential algebraic systems are elaborated.

WP1A-2-210.pdf


Magnetic field calculations with second order accuracy using linear finite elements

Raffaele Albanese1,2, Marco Neri1,2,3, Pasquale Zumbolo1,2

1Università degli studi di Napoli Federico II; 2Consorzio CREATE; 3Università degli studi di Padova

This paper shows an innovative method to evaluate the magnetic field in the vacuum in the presence of ferromagnetic materials with second order accuracy, still using the finite element method (FEM) with linear triangular elements. The procedure is the exten-sion of a technique applied to the solution of magnetohydrodynamic equilibria, based on Helmholtz’s theorem. In this way, post-processing the standard FEM solution, continuity of the magnetic field is ensured by using nodal shape functions for the field compo-nents. The effectiveness of the method is shown in 2D cases for which analytical solutions are available. The possible extension to the calculation of 3D fields is also discussed.

WP1A-3-211.pdf


An Ultra-Fast Numerical Model for AC Losses Evaluation in Superconducting Coils made of HTS tapes

Hocine Menana

University of Lorraine, France, France

In this work, an integrodifferential model in terms of the electric vector potential is combined with the method of images for the electromagnetic field modeling in coils made of high temperature superconducting tapes. The modeling results are compared to measurements achieved on a Rare-earth Barium Copper Oxide (ReBCO) tape-based coil.

WP1A-4-280.pdf


A Framework to Model Skewing of Synchronous Machine for Digital Twin Applications

Victor Mukherjee, Tommi Ryyppo

ABB Oy, Finland

This paper presents a novel framework for modeling rotor skewing in synchronous machines, specifically Permanent Magnet Synchronous Motors (PMSM) and Synchronous Reluctance Motors (SYNRM), suitable for digital twin applications. The method employs a multi-slice approach to approximate the 3D skewing phenomenon, iteratively determining the optimal line current and rotor angle for each slice to achieve target voltage and shaft power. A key feature of this framework is its ability to handle the cross-dependencies between rotor angle, line current, voltage, and power through an adaptive step-size optimization algorithm based on Nelder-Mead. This ensures accurate and stable convergence to the desired operating point, enabling precise simulation of skewed rotor performance, particularly the reduction of torque harmonics. The proposed framework facilitates the development of high-fidelity digital twins for synchronous machines, enabling advanced analysis and optimization of their design and operation.

WP1A-5-365.pdf


AC Losses Modeling in ReBCO Superconducting Coils Using the Volume Integral Method

Sara Fawaz1, Hocine Menana2

1ESME Sudria, France; 2University of Lorraine, France

In this work, different strategies are compared for AC losses modeling in superconducting coils made of Rare-earth Barium Copper Oxide (ReBCO) tapes, using the volume integral method. These strategies explore the possibilities of homogenization of the electromagnetic properties of the tape, in order to reduce computation time. The modeling results are compared to measurements.

WP1A-6-386.pdf


A Novel Energy-Density Based Approach: From Confinement to Radiation in Cylindrical EM Structures

Andrea Pirisi, Riccardo Enrico Zich

Politecnico di Milano, Italy

This work analyzes the transition between quasi-stationary and propagative regimes in electromagnetic wave propagation, considering cylindrical capacitor and inductor configurations as a case study. Starting from the conditions defined by Max Abraham and the subsequent contributions of Schelkunoff—who emphasized energy balance as a crucial variable—the spatio-temporal energy density as a fundamental state variable is here presented. By applying Maxwell's equations through Hertz's method, a rigorous framework to describe energy confinement phenomena in quasi-stationary regimes and their transition to propagative behaviors has been outlined with respect to any cylindrical geometry. The results aim to provide guidelines for designing and optimizing wireless power systems, reconfigurable metamaterials, and advanced electromagnetic devices operating in stratified or anisotropic media.

WP1A-7-390.pdf


Volume integral eddy current problem as a Friedrichs’ system

Anthony Torosyan1, Bertrand Bannwarth1, Olivier Chadebec1, Jean-Michel Guichon1, Gérard Meunier1, Jean-René Poirier2

1Univ. Grenoble Alpes, CNRS, Grenoble INP, G2Elab, 38000 Grenoble, France; 2LAPLACE, Université de Toulouse, CNRS, INPT, UPS, Toulouse, France

The paper introduces Discontinuous Galerkin discretization of a volume integral eddy current problem by recasting it into a Friedrichs’ system. A numerical example is given at the end.

WP1A-8-421.pdf


Easy Computation of the Full Helmholtz Decomposition of Tangential Fields on General Multiply Connected Surfaces using Deterministic and Randomized Algorithms

Marc Bakry, Edouard Demaldent

CEA, France

In this work, we construct a full Helmholtz-Hodge decomposition for tangential Hdiv-conforming fields on general multiply connected, eventually open, surfaces using either a singular value decomposition or a randomized algorithm for the computation of the global loop functions. In particular, our method yields explicit coefficients for the global loop functions spanning the harmonic subspace in the RWG basis in the same fashion as for the loop-tree or loop-star decompositions. This is particularly interesting for users who already have a code based on such decompositions as the implementation of our two distinct algorithms can be done in a few lines of code. We give a numerical example in the context of the simulation of the Eddy current-based non-destructive testing of metallic structures using Boundary Integral Equations.

WP1A-9-424.pdf


Exact Modal Decomposition of Axial Currents in a Conducting Cylinder

Nicola Isernia

Università degli Studi di Napoli Federico II, 80125 Napoli, Italy

The modal decomposition of eddy current problems is commonly treated numerically via integral formulations or via the imposition of arbitrary homogeneous boundary conditions at some distance from the sources in case of differential formulations. In this paper the eigenvalues and eigenvectors for axial currents in an infinitely long conducting cylinder surrounded by vacuum are computed analytically, rigorously accounting for the solution of the outer magnetostatic problem. The analytically derived eigenvalues and eigenvectors are validated against an ad hoc numerical model based on an integral formulation, showing that our analytical solution can be used as a reference for other general-purpose numerical models. The procedure highlighted is general and could be extended to different geometries.

WP1A-10-587.pdf


T - A formulation for superconducting HTS tapes: the Good, the Bad and the Ugly

Fabio Freschi1,2, Maurizio Repetto1, Laura Savoldi1, Sofia Viarengo3

1Dipartimento Energia "G. Ferraris", Politecnico di Torino, Italy; 2School of Information Technology and Electrical Engineering, University of Queensland, Brisbane, QLD, Australia; 3Proxima Fusion, Villigen, Switzerland

This paper evaluates the T− A formulation, a hybrid approach for modeling the electromagnetic behavior of high-temperature

superconducting (HTS) tapes. The method leverages a geometric multi-scale approach by combining a 2D discretization for the

superconducting domain, represented by the electric vector potential T , with a 3D discretization for the surrounding domain,

described by the magnetic vector potential A. The T-A formulation efficiently handles the high width-to-thickness ratio of HTS

tapes, reducing the computational complexity of volumetric modeling while avoiding singularities in the governing equations through

resistance-based matrices.

WP1A-11-596.pdf


Virtual Casing Principle in presence of sink/source currents

Nicola Isernia1, Raffaele Sparago2,3, Francisco Javier Artola3, Matthias Hoelzl2, Guglielmo Rubinacci1, Nina Schwarz3, Salvatore Ventre4, Fabio Villone1

1Università degli Studi di Napoli Federico II, 80125 Napoli, Italy; 2Max Planck Institute for Plasma Physics, 85748 Garching bei Munchen, Germany; 3ITER Organization, 13067 St Paul Lez Durance Cedex, France; 4Università degli Studi di Cassino e del Lazio Meridionale, 03043 Cassino, Italy

The possibility of describing the magnetic field due to currents within a finite control volume via a set of equivalent surface current on the boundary of such volume is given by the well-known Virtual Casing Principle. This was illustrated and implemented in the past for the case of solenoidal currents in the volume of interest, i.e. when the currents to represent have a null normal component to the boundary of the control volume. Here we illustrate how to apply and implement the Virtual Casing Principle for the more general case where the currents to represent can cross the boundary of the control volume.

WP1A-12-616.pdf


Modeling and Simulation of Mechanical Energy Harvesting Using the Feedback Preisach Model of Hysteresis

Carmine S. Clemente1, Daniele Davino2, Ciro Visone3, Amr Adly4

1University of Molise, Campobasso, Italy; 2University of Sannio,Benevento, Italy; 3University of Napoli, Napoli, Italy; 4Egypt-Japan University of Science and Technology, Alexandria, Egypt

Over the past few decades, energy harvesting has been a topic of interest. Recently, magnetoelastic materials have been utilized in a broad spectrum of applications including actuation, sensing, vibration damping and mechanical vibration energy harvesting. This paper aims to present an efficient simulation and modeling approach for electrical energy harvesting from mechanical vibrations. The proposed approach utilizes the scalar feedback Preisach-type model of hysteresis cascaded to a scaling polynomial block. Details of the proposed methodology, experimental testing and simulation results are given in the paper.

WP1A-13-389.pdf