Fiche individuelle
Yvonnick LE MENACH  
Titre  Professeur  
Equipe  Outils et Méthodes Numériques  
Adresse  Université Lille 1 Batiment P2 59655 VILLENEUVED'ASCQ  
Téléphone  +33 (0)320337790  
yvonnick.lemenach@univlille1.fr  
Publications 
ACLI Revue internationale avec comité de lecture 

[1] Study of shaft voltage of a simplified synchronous generator International Journal of Applied Electromagnetics and Mechanics (IJAEM), 08/2018 DARQUES Kévin, TOUNZI Abdelmounaïm, LE MENACH Yvonnick, BEDDEK Karim 
[2] Comparison of Numerical Error Estimators for Eddy Current Problems solved by FEM IEEE Transactions on Magnetics, Vol. 54, N°. 3, pages. 7401204, 03/2018, Abstract TITTARELLI Roberta, LE MENACH Yvonnick, PIRIOU Francis, CREUSE Emmanuel, NICAISE Serge, DUCREUX JeanPierre 
In the domain of field computation with the finite element method, choosing the mesh refinement is an important step to obtain an accurate solution. In order to evaluate the quality of the mesh, a posteriori error estimators are frequently used. In this communication we propose to analyze and to compare residual and equilibrated error estimators for eddy current problems in the case of Aφ and TΩ formulations. The different properties of the estimators will be discussed. 
[3] Numerical Modeling of PCB Planar Inductors: Impact of 3D Modeling on High Frequency Copper Loss Evaluation IET Power Electronics, Vol. 10, N°. 14, pages. 19661974, 11/2017, URL, Abstract TAYLOR Laurent, MARGUERON Xavier, LE MENACH Yvonnick, LE MOIGNE Philippe 
Loss values are key parameters for designing high performance high frequency (HF) magnetic components for power electronics (PE) converters. With the increase of PE switching frequencies, copper losses have to be precisely quantified, ideally until some megahertz. In the literature, many 2D numerical simulations based on finite element analysis (FEA) are performed for such computations. 3D FEA studies of planar components are still limited because of modeling problems, computational resources and computing time. In this paper, quantitative comparisons between 2D and 3D simulation results for planar inductors are achieved focusing on copper loss computation. Results are compared in terms of simulation performances and accuracy. The aim of the paper is to highlight benefits of 2D and 3D FEA simulations in order to
choose the appropriate model according to the studied problem. 
[4] Computation of Magnetic Forces Using Degenerated Airgap Element IEEE Transactions on Magnetics, Vol. 53, N°. 6, pages. 7401304, 06/2017, URL, Abstract SATHYAN Sabin, BELAHCEN Anouar, JUHANI Kataya, HENROTTE François, BENABOU Abdelkader, LE MENACH Yvonnick 
This paper presents an efficient method to calculate the magnetic forces on bodies in contact. The forces are computed through the local application of the virtual work principle on degenerated airgap elements. The results from this method are compared with those from other software and validated with measurements on a permanent magnet setup. Not only is this technique very accurate, but it also reduces the computational burden related to the problematic meshing of thin layers. The implementation of this method in an open source finite element software having facilities for higher order elements and parallel computation unlocks a cost effective and effectual platform for an electromechanical computation of electromechanical devices and magnetic materials. 
[5] Thermal Topology Optimization of a ThreeLayer Laminated Busbar for Power Converters IEEE Transactions on Power Electronics, Vol. 32, N°. 6, pages. 46914699, 06/2017, Abstract PUIGDELLIVOL Oriol, MÉRESSE Damien, LE MENACH Yvonnick, HARMAND Souad, WECXSTEEN JeanFrançois 
This paper focuses on a topology optimization method for laminated busbars in power converters that minimizes the quantity of copper used while keeping the temperature under the allowed limits. Busbars are widely studied for adding their stray inductance to the commutation loop, which causes surge voltage across the power devices. However, the study of heat dissipation is essential to control hotspots in the busbar and preserve the converter components. Current density and temperature are sensitive to shape modifications; hence, topology optimization based on multiphysics simulations is an aspect to be considered when designing prototypes for a good cost performance ratio. The temperature is calculated by an electrothermal twodimensional (2D) finite element method (FEM) superposition approach. Busbar plates are modeled in 2D since the thickness is constant. Furthermore, the different layers are related by the thermal equations reproducing the heat transfers regarding the overlap in the laminated busbar. Simulation results are validated by experimental tests. Comparison with 3D FEM proves the 2D approach to be faster while remaining accurate and a perfect method for topology optimization resolutions, which are very time consuming for threedimensional (3D) geometries. The busbar topology optimization is made by maximizing the energy transfer with the environment and by varying the electric and thermal conductivities of the mesh elements. Optimization leads to more than 50% volume reduction. 
[6] A method coupling modified vector potential A* and homogenization formulations to model short circuits in lamination stacks The European Physical Journal  Applied Physics (EPJ AP), Vol. 75, N°. 3, pages. 11, 08/2016, Abstract ZIANI Smail, HENNERON Thomas, PUIGDELLIVOL Oriol, LE MENACH Yvonnick 
In this paper a method in 2D frequency domain is presented to simulate a laminated iron core with a shortcircuit between several magnetic sheets. The approach consists in coupling homogenization methods and finite element method. The defect is modeled with A* modified vector potential formulation and the rest of the structure with a homogenization method. The coupled method is applied to a lamination stack containing a shortcircuit and compared to the reference, where the A* formulation is applied on the whole domain. Finally, a thermal modeling of lamination stack is presented to study the influence of an insulating defect. 
[7] Estimation of FEM Model Parameters Using Data Assimilation and Its Application to an Electrical Machine IEEE Transactions on Magnetics, Vol. 52, N°. 3, 05/2016, URL BACCHUS Alexandre, TOUNZI Abdelmounaïm, ARGAUD JeanPhilippe, BOURIQUET Bertrand, BIET Mélisande, MACAIRE Ludovic, LE MENACH Yvonnick 
[8] ModelOrder Reduction of Magnetoharmonic Problems Based on POD. Application to Planar Magnetic Components The European Physical Journal  Applied Physics (EPJ AP), Vol. 74, N°. 1, pages. 10903, 04/2016, URL, Abstract TAYLOR Laurent, HENNERON Thomas, MARGUERON Xavier, LE MENACH Yvonnick, LE MOIGNE Philippe 
Predetermination of losses and inductance values in the design phase, is necessary for the development
of highperformance magnetic components for power electronics. Numerical modeling, based on
the Finite Element Method (FEM) can be used to determine the characteristics of a particular component
with a complex geometry in high frequency (HF). These models are very accurate but the computation
time required is high compared to analytical models. The Model Order Reduction (MOR) methods can
be applied to reduce the computation time while maintaining high accuracy. Nowadays, the Proper Orthogonal
Decomposition (POD) is the most popular of MOR approaches. This technique has been applied
to study problems in many elds of engineering. In this paper, the POD method is developed to solve
magnetoharmonic problems in order to study a planar magnetic inductor. 
[9] Finite Element Implementation and Experimental Validation of 2D/3D Magnetic Force Formulas IEEE Transactions on Magnetics, Vol. 52, N°. 2, 03/2016, URL, Abstract BOUGHANMI Walid, HENROTTE François, BENABOU Abdelkader, LE MENACH Yvonnick 
A framework for the finite element implementation of the formulas of local and global electromagnetic forces is presented. The theoretical framework, based on the Lie derivative, is applicable to 2D and 3D problems in the presence of arbitrary materials. The numerical properties of the proposed implementation, in terms of accuracy and convergence, are analyzed. An experimental setup is also presented for which accurate measurements have been carried out, and that can serve as a benchmark for the implementation. 
[10] A posteriori residual error estimators with mixed boundary conditions for quasistatic electromagnetic problem The International Journal for Computation and Mathematics in Electrical and Electronic Engineering (COMPEL), Vol. 34, N°. 3, pages. 724739, 07/2015, URL TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis, NEMITZ Nicolas 
[11] Residual a posteriori estimator for magnetoharmonic potential formulations with global quantities for the source terms IEEE Transactions on Magnetics, Vol. 51, N°. 3, 06/2015, Abstract TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis 
In the modeling of eddy current problems, potential formulations are widely used in recent days. In this paper, the results of residualbased a posteriori error estimators, which evaluate the discretization error in the finiteelement computation, are extended to the case of several kinds of source terms for both A/φ and T/Ω harmonic formulations. The definitions of the estimators are given and some numerical examples are provided to show the behavior of the estimators. 
[12] Finite element mesh adaptation strategies from residual and hierarchical error estimators in eddy current problems IEEE Transactions on Magnetics, Vol. 51, N°. 3, 04/2015, URL, Abstract DULAR Patrick, LE MENACH Yvonnick, TANG Zuqi, CREUSE Emmanuel, PIRIOU Francis 
A strategy of mesh adaptation in eddy current finite element modeling is developed from both residual and hierarchical error estimators. Wished distributions of element sizes of adapted meshes are determined from the elementwise local contributions to the estimators and define constraints for the mesh generator. Uniform distributions of the local error are searched. 
[13] Spacetime residualbased a posteriori estimator for the Αφ formulation in eddy current problems. Magnetics, IEEE Transactions, Vol. 51, N°. 3, 03/2015, Abstract TITTARELLI Roberta, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis, MOREAU Olivier, BOITEAU Olivier 
In this work, an a posteriori residual error estimator is presented for the 3D eddy current problem modeled by the spacetime Aphi potential formulation. It is solved by the finite element method in space and the backward Euler scheme in time. Once the reliability as well as the efficiency of the estimator are established, two numerical tests are proposed: an analytical one in order to validate the theoretical results and a physical one in order to illustrate the performance for a real eddy current problem. 
[14] Comparison of Residual and Hierarchical Finite Element Error Estimators in Eddy Current Problems IEEE Transactions on Magnetics, Vol. 50, N°. 2, pages. 501504, 02/2014, URL, Abstract DULAR Patrick, TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, PIRIOU Francis 
The finite element computation of eddy current problems introduces numerical error. This error can only be estimated. Among all error estimators (EEs) already developed, two estimators, called residual and hierarchical EEs, proven to be reliable and efficient, are theoretically and numerically compared. Both estimators show similar behaviors and locations of the error. 
[15] Residual and equilibrated error estimators for magnetostatic problems solved by finite element method IEEE Transactions on Magnetics, Vol. 49, N°. 12, pages. 5715  5723, 12/2013, URL, Abstract TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis 
In finite element computations, the choice of the mesh is crucial to obtain accurate solutions. In order to evaluate the quality of the mesh, a posteriori error estimators can be used. In this paper, we develop residualbased error estimators for magnetostatic problems with both classical formulations in term of potentials used, as well as the equilibrated error estimator. We compare their behaviors on some numerical applications, to understand the interest of each of them in the remeshing process. 
[16] A posteriori error estimator for harmonic APhi formulation The international journal for computation and mathematics in electrical and electronic engineering, Vol. 32, N°. 4, pages. 1219  1229, 07/2013, URL TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis 
[17] Solution of large stochastic finite element problem – application to ECTNDT IEEE Trans. on. Mag,, Vol. 49, N°. 5, pages. 16051608, 05/2013 BEDDEK Karim, CLENET Stéphane, MOREAU Olivier, LE MENACH Yvonnick 
[18] An Arbitrary Thick Shell Finite Element for EddyCurrent Dual VectorScalar Potential Formulations IEEE Trans. Mag, Vol. 49, N°. 5, pages. 17251728, 05/2013 THOMAS Pierre, LE MENACH Yvonnick 
[19] Test Harness on a Preconditioned Conjugate Gradient Solver on GPUs: An Efficiency Analysis IEEE Trans. Mag, Vol. 49, N°. 5, pages. 17291732, 05/2013 RODRIGUES A. Wendell de O., CHEVALLIER Loïc, LE MENACH Yvonnick, GUYOMARC’H Frédéric 
[20] Stochastic NonDestructive Testing simulation: sensitivity analysis applied to material properties for tube support plate clogging in nuclear power plant steam generator IEEE Transactions on Magnetics, Vol. 49, N°. 9, pages. 18731876, 05/2013 MOREAU Olivier, BEDDEK Karim, CLENET Stéphane, LE MENACH Yvonnick 
[21] Residual based a posteriori error estimators for harmonic A/Phi and T/Omega formulations in eddy current problems IEEE Transactions on Magnetics, Vol. 49, N°. 5, 05/2013, URL TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis 
[22] Residualbased a posteriori estimators for the T/Omega magnetodynamic harmonic formulation of the Maxwell system International Journal of Numerical Analysis and Modeling, Vol. 10, N°. 2, pages. 411429, 02/2013, URL, Abstract CREUSE Emmanuel, NICAISE Serge, TANG Zuqi, LE MENACH Yvonnick, NEMITZ Nicolas, PIRIOU Francis 
In this paper, we focus on an a posteriori residualbased error estimator for the T/Omega
magnetodynamic harmonic formulation of the Maxwell system. Similarly to the A/Phi formulation, the weak continuous and discrete formulations are established, and the wellposedness of
both of them is addressed. Some useful analytical tools are derived. Among them, an adhoc
Helmholtz decomposition for the T/Omega
case is derived, which allows to pertinently split the error.
Consequently, an a posteriori error estimator is obtained, which is proven to be reliable and locally
efficient. Finally, numerical tests confirm the theoretical results. 
[23] Residualbased a posteriori estimators for the A/Phi magnetodynamic harmonic formulation of the Maxwell system Mathematical Models and Methods in Applied Sciences, Vol. 22, N°. 5, pages. 30, 05/2012, URL, Abstract CREUSE Emmanuel, NICAISE Serge, TANG Zuqi, LE MENACH Yvonnick, NEMITZ Nicolas, PIRIOU Francis 
This paper is devoted to the derivation of an a posteriori residualbased error estimator for the APhi magnetodynamic harmonic formulation of the Maxwell system. The weak continuous and discrete formulations are established, and the wellposedness of both of them is addressed. Some useful analytical tools are derived. Among them, an ad hoc Helmholtz decomposition is proven, which allows to pertinently split the error. Consequently, an a posteriori error estimator is obtained, which is proven to be reliable and locally efficient. Finally, numerical tests confirm the theoretical results. 
[24] A threedimensional electromagnetic shell finite element for coupled vectorscalar potential formulations IEEE Trans. mag, Vol. 48, N°. 2, pages. 823826, 02/2012 THOMAS Pierre, LE MENACH Yvonnick 
[25] Automatic MultiGPU Code Generation Applied to Simulation of Electrical Machines IEEE Trans. Mag, Vol. 48, N°. 2, pages. 823826, 02/2012 RODRIGUES A. Wendell de O., GUYOMARC’H Frédéric, DEKEYSER JeanLuc, LE MENACH Yvonnick 
[26] Adaptive method for nonintrusive spectral projection application on a stochastic eddy current NDT problem IEEE Transactions on Magnetics, Vol. 48, N°. 2, pages. 759762, 02/2012, Abstract BEDDEK Karim, CLENET Stéphane, MOREAU Olivier, COSTAN Valentin, LE MENACH Yvonnick, BENABOU Abdelkader 
The NonIntrusive Spectral Projection (NISP) method is widely used for uncertainty quantification in stochastic models. The determination of the expansion of the solution on the polynomial chaos requires the computation of multidimensional integrals. An automatic adaptive algorithm based on nested sparse grids has been developed to evaluate those integrals. The adapted algorithm takes into account the weight of each random variable with respect to the output of the model. To achieve that it constructs anisotropic sparse grid of the mean, leading to a reduction of the number of numerical simulations. Furthermore, the spectral form of the solution is explicitly identified from the constructed quadrature scheme. Numerical results obtained on an industrial application in NDT demonstrate the efficiency of the proposed method. 
[27] Spectral stochastic finite element method for solving 3D stochastic eddy current problems International Journal of Applied Electromagnetics and Mechanics, Vol. 39, N°. 14, pages. 753760, 01/2012, Abstract BEDDEK Karim, CLENET Stéphane, MOREAU Olivier, LE MENACH Yvonnick 
In this paper, a spectral stochastic finite element method (SSFEM) is described and applied to quantify uncertainties in an industrial eddy current non destructive testing (NDT) problem. A new iterative solver based on the fixedpoint method is proposed to solve the large scale linear system. The efficiency of the suggested iterative solver to tackle the SSFEM computational complexity is demonstrated. The SSFEM results and the computation times obtained with 2 and 4 random variables are compared to those given by the anisotropic adaptive sparse grid non intrusive spectral projection (NISP). On this application, the SSFEM dedicated solver is competitive with the adaptive sparse grids strategy. 
[28] High Order Surface Impedance Boundary Conditions with the AΦ Formulation Facta Universitatis: Electronics and Energetics, Vol. 24, N°. 2, pages. 147155, 08/2011 IDA Nathan, LE MENACH Yvonnick, HENNERON Thomas 
[29] 3D stochastic Spectral Finite Element Method in static electromagnetism using vector potential formulation IEEE trans. on Mag., Vol. 47, N°. 5, pages. 12501253, 05/2011, Abstract BEDDEK Karim, LE MENACH Yvonnick, CLENET Stéphane, MOREAU Olivier 
Stochastic spectral finiteelement method can be used to take into account some random aspects in the input data (material characteristic, source terms) involved in static electromagnetism problems. Similarly to the deterministic case, two potential formulations can be used in the stochastic case. The vector potential formulation applied to static problems is developed and compared to the scalar potential one, previously developed. 
[30] The fieldcircuit analysis of the startup operation of the brushless DC motor Computer Applications in electrical Engineering,, pages. 93104, 11/2010 KNYPINSKI L, NOWAK L, RADZIUK K, LE MENACH Yvonnick 
[31] Parallel direct solver for the finite integration technique in electrokinetic problem IEEE Trans mag, Vol. 46, N°. 8, pages. 32693272, 08/2010 TINZEFTE Abdellatif, LE MENACH Yvonnick, KORECKI Julien, GUYOMARC’H Frédéric, PIRIOU Francis 
[32] Hysteresis Phenomenon Implementation in FIT: Validation with Measurements IEEE Transactions on Magnetics, Vol. 46, N°. 8, pages. 32853288, 08/2010, URL, Abstract KORECKI Julien, BENABOU Abdelkader, LE MENACH Yvonnick, DUCREUX JeanPierre, PIRIOU Francis 
The presented work is about the implementation of a vector hysteresis model in a Finite Integration Technique (FIT) calculation code. This is achieved in the magnetostatic case using the vector potential formulation and with the external electric circuit coupling. The used behavior law model is the vectorized JilesAtherton model with the magnetic flux density as input variable. The proposed approach is illustrated by modeling a coil made of Soft Magnetic Composite (SMC) magnetic core. Comparison with measured global quantities is also performed. 
[33] Study of Static and Dynamic Eccentricities of a Synchronous Generator Using 3D FEM IEEE Transactions on Magnetics, Vol. 46, N°. 8, 08/2010, Abstract TANNOIAMAMURA Bruno, LE MENACH Yvonnick, TOUNZI Abdelmounaïm, SADOWSKI Nelson, GUILLOT Eilin 
This paper deals with the study of static and dynamic rotor eccentricities of a synchronous generator using 3D FEM. First, both eccentricity cases are introduced as well as the used approach. Then the models of the structure, based on the magnetic scalar potential formulation and on the magnetic vector potential formulation, are presented. Two eccentricity modeling methods are introduced, the 2D eccentricity and the 3D eccentricity. Results obtained for magnetostatic and magnetodynamic cases related to both defects and their combination are presented and discussed. 
[34] Parallel direct solver for the finite integration technique in electromagnetics COMPEL, Vol. 29, N°. 4, pages. 941949, 07/2010 TINZEFTE Abdellatif, LE MENACH Yvonnick, KORECKI Julien, PIRIOU Francis 
[35] Iterative Solvers For Singular Symmetric Linear Systems in Electromagnetism Low Frequency IEEE Trans Mag, Vol. 45, N°. 3, pages. 1428  1431, 03/2009 TINZEFTE Abdellatif, LE MENACH Yvonnick, PIRIOU Francis 
[36] Study of an electromagnetic gearbox involving two permanent magnet synchronous machines using 3DFEM IEEE Trans Mag, Vol. 44, N°. 11, pages. 4381  4384, 10/2008 AUBERTIN Mathieu, TOUNZI Abdelmounaïm, LE MENACH Yvonnick 
[37] Implementation of a Vector Hysteresis Model in 2D Finite Element Analysis: Study of a RSST with Anisotropic Sample International Journal of Applied Electromagnetics and Mechanics, Vol. 28, N°. 12, pages. 4147, 09/2008, Abstract LEITE Jean Vianei, BENABOU Abdelkader, SADOWSKI Nelson, LE MENACH Yvonnick, CLENET Stéphane 
This work deals with the implementation of a vector hysteresis model in a 2D finite element code with vector potential formulation. The inverse JilesAtherton hysteresis model leads to obtain a differential reluctivity tensor which can be directly applied to the field calculation equations. 
[38] A mixed finite element/meshless natural element method for simulating rotative electromagnetic machines EPJ Appl. Phy, Vol. 43, pages. 197208, 09/2008, Abstract ILLOUL Lounes, LE MENACH Yvonnick, CLENET Stéphane, CHINESTA Francisco 
Efficient numerical models have been already developed in 2D to take into account the movement
of electromagnetic devices with rotating parts in the framework of the Finite Element Method (FEM).
When the movement becomes complex it leads to large mesh distortions.A remeshing step is then required
which increases the computational complexity and can also lead,in some cases,to numerical ripples
on forces and torques due to the ﬁeld projections between old and new meshes.Moreover,remeshing
procedures in 3D remain an open topic.Meshless methods seem an appealing choice for alleviating the
mesh constraints.The Natural Element Method (NEM)which,has known a growing interest in the domain
of mechanics,allows to proceed in the meshless framework,avoiding one of the main drawbacks related
to the vast majority of meshless techniques,as is the imposition of essential boundary conditions.In
this paper,a variant of the NEM,known as constrained natural element method (CNEM)is applied for
simulating electromagnetic machines involving rotating parts.A new mixed strategy combining the ﬁnite
element and the constrained natural element methods is proposed and then tested by using an appropriate
error estimator. 
[39] Comparison between the Mortar Element Method and the polynomial Interpolation Method to model movement in the Finite Element Method IEEE, Vol. 44, N°. 6, pages. 1311317, 06/2008 SHI Xiaodong, LE MENACH Yvonnick, DUCREUX JeanPierre, PIRIOU Francis 
[40] Numerical solutions in primal and dual mesh of magnetostatic problem
solved with the Finite Integration Technique
COMPEL, Vol. 27, N°. 1, pages. 4755, 01/2008, Abstract KORECKI Julien, LE MENACH Yvonnick, DUCREUX JeanPierre, PIRIOU Francis 
A paraître Abstract  To solve the Maxwell’s equations in the magnetosatic case the formulations in terms of vector potential or scalar potential can be used. If to solve the equations we use the Finite Element Method on the primal mesh we obtain two complementary solutions. In this paper, both formulations in terms of potential and using the fields H and B are developed and solved with the help of the Finite Integration Technique. The solutions computed in the primal and dual meshes are analysed and discussed. 
[41] Source Field Computation in NDT Applications IEEE Trans. Mag., Vol. 43, N°. 4, pages. 17851788, 04/2007 HENNERON Thomas, LE MENACH Yvonnick, PIRIOU Francis, CLENET Stéphane, DUCREUX JeanPierre 
[42] COMPARISON OF SLIP SURFACE AND MOVING BAND TECHNIQUES FOR MODELLING MOVEMENT IN 3D WITH FEM COMPEL, Vol. 25, N°. 1, pages. 1730, 01/2006, Abstract SHI Xiaodong, LE MENACH Yvonnick, DUCREUX JeanPierre, PIRIOU Francis 
In this paper the slip surface and the moving band techniques are used to model the rotation of electrical machines in 3D with FEM and the results are compared. For the slip surface technique, several method are explained; such as the lockedstep approach and linear and quadratic interpolation. For this study, the proposed techniques are applied to a permanent magnet synchronous machine. The comparison is carried out at noload for the electromotive force and the cogging torque. Furthermore, the torque is also compared for the shortcircuit case. 
[43] Consideration of the coupling of the magnetic and electric equations with Finite Integration Technique (FIT)
European Physical Journal Applied Physics, Vol. 30, N°. 1, pages. 1721, 04/2005, Abstract ALBERTIER Grégory, LE MENACH Yvonnick, DUCREUX JeanPierre, PIRIOU Francis 
In this paper, the Finite Integration Technique and an approach to solve simultaneously the magnetic and electric circuit equations is presented in magnetostatic case. The developments of the FIT will be compared with those used for the Finite Element Method. As example of application an iron core coil is studied and the results of both approaches are compared 
[44] 3D approaches to determine the end winding inductances of a PMLSM IEEE Trans. Mag, Vol. 38, N°. 2, pages. 989992, 03/2004 TOUNZI Abdelmounaïm, HENNERON Thomas, LE MENACH Yvonnick, ASKOUR Rachid, DUMETZ Eric 
[45] Numerical model to discretize source fields in the 3D finite elements method IEEE Transactions Magnetics, Vol. 36, N°. 4, pages. 676679, 07/2000, Abstract LE MENACH Yvonnick, CLENET Stéphane, PIRIOU Francis 
This communication presents a method to discretize electromagnetic source fields in the Whitney element spaces for static electromagnetic potential formulations. This method requires no finite element solution to calculate this source fields. It is based on the use of a facet or an edge tree. We will test it on electrokinetic, electrostatic and magnetostatic examples 
[46] 3D computation of a claw pole magnet machine using a scalar potential formulation European Physical Journal Applied Physics, Vol. 12, N°. 3, pages. 175182, 01/2000, Abstract LE MENACH Yvonnick, CLENET Stéphane, PIRIOU Francis, LETOURNEAU JeanFrançois, VIAROUGE Philippe 
In this paper, a claw pole permanent magnet machine is studied using a 3D scalar potential formulation. First, the discretization of Maxwell’s equations using Whitney’s elements is presented. Especially, the method to take the winding current distribution into account by means of source fields is detailed. Then, a numerical model of the studied machine is presented. Simulation results (inductance, EMF, torque) are discussed by comparison with experimental ones.

[47] 3D compatible magnetostatic potential formulations coupled with electrical circuits
COMPEL, Vol. 13, N°. 3, pages. 776786, 01/2000, Abstract LE MENACH Yvonnick, CLENET Stéphane, PIRIOU Francis 
This paper presents a method of coupling magnetostatic potential formulations with electrical circuits for the 3D FEM. Allowing for the current density distribution in stranded conductors, two vectors N and K are introduced. A systematic method of decomposing both vectors in Whitney’s element spaces is suggested. This method can be used for coils with a complex shape. As an example of application a threephase transformer is studied.

[48] Coupling of magnetic and electrical circuit equations using coil shape functions International Journal of Theorical Electrotechnics, pages. 3639, 03/1999, Abstract LE MENACH Yvonnick, CLENET Stéphane, PIRIOU Francis 
In this paper we propose, an approach for coupling magnetostatic and electric circuit equations applicable to the case of stranded conductors. The method based on the calculation of global quantities is presented and applied to the magnetic formulation HsƒÇ. The properties and calculation of source field Hs are also discussed. As an example we have studied a permanent magnet synchronous machine. At no load, the magnetomotive force is determined and compared with experimental results. A calculation when the generator has a resistance load is also presented 
[49] Determination and utilization of the source field in 3D magnetostatic problems IEEE Transactions Magnetics, Vol. 34, N°. 5, pages. 25092512, 09/1998, Abstract LE MENACH Yvonnick, CLENET Stéphane, PIRIOU Francis 
In this paper we study different approaches to introduce the source terms due to the inductors crossed by a uniform current density J 0 in 3D FEM. Two equivalent methods are developed to decompose the current density J in the facet element space with divergence free close to J0. This decomposition is used in the aformulation without gauge condition. Moreover, from the flux facet the source field Hs can be calculated in the edge element space and introduced into the öformulation. As examples of applications we have studied a coil with a complex geometry and a iron core coil 
[50] Two guaranteed equilibrated error estimators for Harmonic formulations in eddy current problems Computers & Mathematics with Applications, A paraître, URL, Abstract CREUSE Emmanuel, LE MENACH Yvonnick, NICAISE Serge, PIRIOU Francis, TITTARELLI Roberta 
In this paper a guaranteed equilibrated error estimator is developed for the 3D harmonic magnetodynamic problem of Maxwell’s system. This system is recasted in the classical Aphi potential formulation and solved by the Finite Element method. The error estimator is built starting from the Aphi numerical solution by a local flux reconstruction technique. Its equivalence with the error in the energy norm is established. A comparison of this estimator with an equilibrated error estimator already developed through a complementary problem points out the advantages and drawbacks of these two estimators. In particular, an analytical benchmark test illustrates the obtained theoretical results and a physical benchmark test shows the efficiency of these two estimators. 
ACLN Revue nationale avec comité de lecture 
[1] Introduction des grandeurs globales dans les problèmes d’électrocinétique résolus avec la F.E.M. et la F.I.T. RIGE, Vol. 11, N°. 23, pages. 185203, 04/2008, Abstract KORECKI Julien, LE MENACH Yvonnick, VERITE JeanClaude, PIRIOU Francis 
Cette communication présente une approche permettant d’imposer le courant ou la différence de potentiel dans le cas d’un problème d’électrocinétique. Pour une formulation en potentiel vecteur, l’approche proposée a été implantée dans un code de calcul basé sur la Technique d’Intégration Finie et dans un code éléments finis. Des simulations sont réalisées pour des cas simples et les résultats sont comparés à ceux obtenus à l’aide de la méthode des éléments finis. 
ACT Conférence internationale avec acte 
[1] An Improved Newton Method Based on choosing Initial Guess Applied to Scalar Potential Formulation in Nonlinear Magnetostatics CEFC 2018, Hangzhou, China, 10/2018 CHERIF Riheb, TANG Zuqi, GUYOMARC’H Frédéric, CHEVALLIER Loïc, LE MENACH Yvonnick 
[2] Application Limits of the Airgap Maxwell Tensor CEFC 2018, 10/2018, Abstract PILE Raphaël, PARENT Guillaume, DEVILLERS Emile, HENNERON Thomas, LE MENACH Yvonnick, LE BESNERAIS Jean, LECOINTE JeanPhilippe 
In an electrical machine, the Maxwell Tensor is widely used to compute global forces or local pressure along a surface in
the air. This communication proposes to highlight the limits of the method with an academic case of slotless stator and rotor.
In particular an analytic demonstration shows the existence of coefficients depending on the geometry and the wavenumber
between the application of the Maxwell Tensor in the airgap and the stator magnetic pressure. 
[3] A Mesh Adaptation Method Based on A Posteriori Error Estimators Applied to Eddy Current Testing CEFC 2018, 10/2018 DEBUSSCHÈRE Lydéric, LE MENACH Yvonnick, CLENET Stéphane, GOURSAUD Benjamin, DUCREUX JeanPierre 
[4] Effect of the end winding on the shaft voltage of a high power alternator EPNC 2018, 06/2018, Abstract DARQUES Kévin, TOUNZI Abdelmounaïm, LE MENACH Yvonnick, BEDDEK Karim 
Shaft voltage is present in all electrical machines and can be useful as a diagnosis variable but first, it is important to quantify the effect of some defects. 2D numerical models have already been used to tackle this point. However, 3D effects due to end windings can have a significant influence on this voltage. This paper deals with the study of a turbogenerator using a 3D finite element model in order to investigate this task. 
[5] An improved starting point for Newton’s method solving 3D nonlinear magnetostatic problems EPNC 2018, Arras, France, 06/2018 CHERIF Riheb, TANG Zuqi, GUYOMARC’H Frédéric, CHEVALLIER Loïc, LE MENACH Yvonnick 
[6] INTEGRATED DRIVES: OPPORTUNITIES, CHALLENGES, CE2I PROJECT PROPOSAL Twentyfifth Symposium on Electromagnetic Phenomena in Nonlinear Circuits (EPNC), 06/2018,, 06/2018, Abstract HARMAND Souad, IDIR Nadir, LECOINTE JeanPhilippe, LEMAIRESEMAIL Betty, LE MENACH Yvonnick, SAUDEMONT Christophe, ROMARY Raphaël, SEMAIL Eric 
Integrated Drives are a target in terms of simplicity in complex systems such as transportation. Anyway the integration of the power component inside the electrical machine induces challenges if the compactness is required. The CE2I project, for “Smart Integrated Drives”, searches to deal simultaneously with the constraints such as better thermal control and modelling, more efficient numerical virtual modelling for cosimulation, exploration of new materials and use of low losses, high frequency widebandgap (WBG) components. 
[7] Harmonicbalance finiteelement homogenisation of laminated iron cores with net circulating currents EMF 2018, Darmstadt, Germany, 04/2018 SABARIEGO Ruth, GYSELINCK Johan, LE MENACH Yvonnick, HENNERON Thomas 
[8] An iterative mesh adaptation method using a size map based on equilibrated error estimator EMF, Darmstadt, 04/2018 DEBUSSCHÈRE Lydéric, LE MENACH Yvonnick, CLENET Stéphane, DUCREUX JeanPierre, GOURSAUD Benjamin 
[9] A method to model the short circuits in the lamination stack Acomen, Gand, 09/2017 ZIANI Smail, HENNERON Thomas, LE MENACH Yvonnick 
[10] Nonlinear Lamination Stacks Studied with Harmonic Balance FEM combined with Homogenization approach Compumag (International Conference on the Computation of Electromagnetic Fields), 03/2017, Abstract ZIANI Smail, HENNERON Thomas, LE MENACH Yvonnick 
Harmonic Balance Finite Element method combined with homogenization method is used to model lamination stacks. The Harmonic Balance gives directly the steadystate solution and the homogenization method reduces the number of unknowns. The numerical model takes into account the nonlinear magnetic behavior and the electric conductivity. The results of the proposed method are compared with those obtained from a classic approach. 
[11] Study of shaft voltage of a simplified synchronous generator ISEM 2017, pages. 2, 03/2017, Abstract DARQUES Kévin, TOUNZI Abdelmounaïm, LE MENACH Yvonnick, BEDDEK Karim 
In large turbogenerators, the analysis of the shaft voltage can be useful to diagnose some defects such as eccentricities or rotor short circuit. To do that, the first step is to quantify the inherent effect of each defect on that voltage. This paper deals with the study of a simplified synchronous generator with a smooth slotless stator using a numerical model in order to well understand the effect of each defect. 
[12] Nonlinear Lamination Stacks Studied with Harmonic Balance FEM Supplied by Magnetic Flux Arising from PWM CEFC 2016, Miami, 11/2016, Abstract ZIANI Smail, HENNERON Thomas, LE MENACH Yvonnick 
The Harmonic Balance combined with the finite element method enables to obtain the steadystate solution of electromagnetic problems. In this communication, we propose to apply this approach to study a laminated iron core submitted to a magnetic flux arising from a Pulse Width Modulation (PWM) voltage. The proposed numerical model takes into account the nonlinear magnetic behaviour, the electric conductivity and a short circuit between steel sheets. 
[13] Computation of Magnetic Forces Using Degenerated Airgap Element CEFC 2016, Miami  EtatsUnis, 11/2016, Abstract SATHYAN Sabin, BELAHCEN Anouar, JUHANI Kataya, HENROTTE François, BENABOU Abdelkader, LE MENACH Yvonnick 
This paper presents an efficient method to calculate the magnetic forces on bodies in contact. The forces are computed through local application of the Virtual Work principle on degenerated airgap elements. The results from this method are compared with those from another software and validated with measurements on a permanent magnet setup. Not only is this technique very accurate, but also reduces the computational burden related to the problematic meshing of thin layers. The implementation of this method in an open source finite element software having facilities for higher order elements and parallel computation unlocks a cost effective and effectual platform for electromechanical computation of electromechanical devices and magnetic materials. 
[14] A coupled method between homogenization and vector potential formulations
to model shortcircuits in lamination stacks MSHOM 2016, Vienna, 09/2016 ZIANI Smail, HENNERON Thomas, LE MENACH Yvonnick 
[15] An approach to model shaft voltage of wound rotor synchronous machines CEFC 2016, pages. 1, 08/2016, Abstract DARQUES Kévin, TOUNZI Abdelmounaïm, LE MENACH Yvonnick, BEDDEK Karim, BIET Mélisande 
In the case of high power synchronous machines, shaft voltage can induce high currents yielding to bearing damages but it can also constitutes a variable to be used to diagnosis some machine defects. Therefore, it is important to well model the phenomenon in order to use it in further diagnosis tool. 
[16] Harmonic Balance Finite Element Method applied for Nonlinear Lamination stack with short circuits EPNC 2016, Helsinki, 06/2016, Abstract ZIANI Smail, MONTIER Laurent, HENNERON Thomas, LE MENACH Yvonnick 
In this paper, the Harmonic Balance Method is used to obtain the steadystate solution of a laminated iron core with an imposed magnetic flux. The numerical model takes into account the nonlinear magnetic behaviour, the electric conductivity and a short circuit between steel sheets. The results of the Harmonic Balance Method are compared with those obtained from a classic approach. 
[17] Geometrical parametric model order reduction of transmission lines EMF 2016, Lyon, France, 04/2016 PIERQUIN Antoine, HENNERON Thomas, LE MENACH Yvonnick, COLAS Frédéric 
[18] A method coupling a modified vector potential A* and a homogenization formulations to
model short circuits in lamination stacks EMF 2016, Lyon, 04/2016 ZIANI Smail, LE MENACH Yvonnick, HENNERON Thomas 
[19] Alternator Rotor Interturn ShortCircuit Identification using FEM Based Learning 9th IFAC SAFEPROCESS 2015 Symposium, 09/2015, Abstract BACCHUS Alexandre, BIET Mélisande, MACAIRE Ludovic, TOUNZI Abdelmounaïm, LE MENACH Yvonnick 
The goal of this paper is to identify rotor interturn shortcircuits of an alternator.
The detection method is based on the analysis of a ux probe signal located in the air gap
of the machine. Previous works have shown that pattern recognition can be applied to detect
such a fault by using the experimental data as prototypes. A new method is developed here by
considering a learning step based on simulation. Therefore the machine is modeled and validated.
A feature selection is made by considering feature correlation and disparity. Finally, k nearest
neighbors is used to classied experimental test data. 
[20] Finite Element Implementation and Experimental Validation of 2D/3D Magnetic Force formulae Compumag 2015, Montreal, Canada, 07/2015 BOUGHANMI Walid, HENROTTE François, BENABOU Abdelkader, LE MENACH Yvonnick 
[21] 3D Modeling of Planar Inductors for Power Electronics Conference on Electromagnetic Field Computation 2014, 05/2014, Abstract TAYLOR Laurent, MARGUERON Xavier, LE MENACH Yvonnick, LE MOIGNE Philippe 
Predetermination of inductances and losses are needed to effectively design magnetic components used in power electronic converters. Their behaviours change with frequency and have to be determined accurately over a wide frequency range, commonly up to 30MHz to account for EMI aspects. In many papers, 2D finite element simulations are widely performed for the calculations of inductances and losses of planar components. These 2D numerical models are generally taken as references to develop analytical models, including the specific case of High Frequency copper losses. 3D studies of such planar components are still limited because they leads to problems related to modeling, computing resources and computing time. In this paper comparisons in terms of performance and precision for 2D and 3D models of planar inductors is presented. The obtained results are compared with impedance measurements on prototypes. 
[22] Error estimation in the computation of induced currents of human body 2013 CIGRE SCC3 & EMFELF Colloquium, Japan, 10/2013 LELONG Thomas, TANG Zuqi, SCORRETTI Riccardo, THOMAS Pierre, LE MENACH Yvonnick, CREUSE Emmanuel 
[23] Comparison of Residual and Hierarchical Finite Element Error
Estimators in Eddy Current Problems COMPUMAG 2013 Budapest, Hungary, 07/2013, URL, Abstract TANG Zuqi, DULAR Patrick, LE MENACH Yvonnick, CREUSE Emmanuel, PIRIOU Francis 
The finite element computation of eddy current
problems gives numerical error. This error cannot be calculated,
but can only be estimated. Among all error estimators already
developed, it is proposed to compare two proven estimators
called residual and hierarchical error estimators. 
[24] Error estimation in the Computation of Induced Current of Human Body in the Case of Low frequency Magnetic Field Excitation COMPUMAG 2013 Budapest, Hungary, 07/2013 LELONG Thomas, TANG Zuqi, SCORRETTI Riccardo, THOMAS Pierre, LE MENACH Yvonnick, CREUSE Emmanuel 
[25] Comparison of supervised classification algorithms combined with feature extraction and selection: Application to a turbogenerator rotor fault detection Diagnostics for Electric Machines, Power electronics and Drives (SDEMPED), 2013 9th IEEE International Symposium on, pages. 558  565, 01/2013, Abstract BACCHUS Alexandre, BIET Mélisande, MACAIRE Ludovic, TOUNZI Abdelmounaïm, LE MENACH Yvonnick 
The goal of this paper consists in applying pattern recognition methods to turbogenerators. Previous works have shown that a monitor, thanks to pattern recognition, is practical on asynchronous machines. This procedure has rarely taken advantage of these methods for turbogenerator. The statistical model has been obtained from harmonics extracted from flux probes and from stator current and voltage. For this purpose, the main way is to build a learning matrix to predict the functional state of a new measurement. Finally, three classifiers have been compared: k Nearest Neighbors, Linear Discriminant Analysis and Support Vector Machines. The best classification result is obtained by Linear Discriminant Analysis combined with Factorial Discriminant Analysis achieving a score of 84.6%. 
[26] Residual based a posteriori error estimator for harmonic APhi and TOmega formulation in eddy current problems CEFC 2012 Oita, Japan, 11/2012, Abstract TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis 
For the eddy current problem, the potential formulations are widely used today. In this communication, residual based a posteriori error estimators are introduced to evaluate the discretization error in the finite element calculation in both case of APhi and TOmega harmonic formulations. An example is carried out to show the behavior of our estimators. 
[27] Test Harness on a Preconditioned Conjugate Gradient Solver on GPUs: An Efficiency Analysis CEFC 2012 International Conference on Electromagnetic Field Computation, Oita, Japan, 11/2012, Abstract RODRIGUES A. Wendell de O., CHEVALLIER Loïc, LE MENACH Yvonnick, GUYOMARC’H Frédéric 
The parallelization of numerical simulation algorithms, i.e., their adaptation to parallel processing architectures, is an aim to reach in order to hinder exorbitant execution times. The parallelism has been imposed at the level of processor architectures and graphics cards are now used for generalpurpose calculation, also known as “GeneralPurpose computation on Graphics Processing Unit (GPGPU)”. The clear benefit is the excellent performance over price ratio. Besides hiding the low level programming, software engineering leads to a faster and more secure application development. This paper presents the real interest of using GPU processors to increase performance of larger problems which concern electrical machines simulation. Indeed, we show that our autogenerated code applied to several models allows achieving speedups of the order of 10x. 
[28] A new preconditioner based on the F.I.T to solve F.E.M problems CEFC 2012 Oita, Japan, 11/2012 TINZEFTE Abdellatif, LE MENACH Yvonnick, GUYOMARC’H Frédéric 
[29] Solution of large stochastic finite element problem – application to ECTNDT CEFC 2012 Oita, Japan, 11/2012 BEDDEK Karim, CLENET Stéphane, MOREAU Olivier, LE MENACH Yvonnick 
[30] An arbitrary thick shell finite element for eddy currents dual vector scalar potential formulations CEFC 2012 Oita, Japan, 11/2012 THOMAS Pierre, LE MENACH Yvonnick 
[31] Impedance Determination in Frequency Domain for Energy Cables by FEM and TLM Proc., INDEL, Banja Luka, 11/2012 NGUYENDUC Q, LE MENACH Yvonnick, CLENET Stéphane, VIZIREANU Darius, COSTAN Valentin 
[32] Damper bars influence on the synchronous machine steady state operating Proc., Momag, Joao PessoaParaiba, 08/2012 TANNOIAMAMURA Bruno, LE MENACH Yvonnick, TOUNZI Abdelmounaïm, SADOWSKI Nelson, GUILLOT Eilin 
[33] A posteriori error estimator for harmonic APhi formulation EPNC 2012 Pula, Croatia, 06/2012, Abstract TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis 
In this paper a residualbased error estimator is proposed to evaluate the numerical error induced by the computation of the electromagnetic systems using a Finite Element Method in the case of the harmonic APhi formulation. This estimator is based on the evaluation of quantities weakly verified by the formulation. Furthermore as this estimator verifies the notions of reliability and efficiency it allows to estimate the qualities of local and global solutions. Two examples of electromagnetic systems with current density induced are used to show the efficiency of the proposed estimator. 
[34] A Nonlinear Model for AC Induced Corrosion AES, Paris, 05/2012 IDA Nathan, LE MENACH Yvonnick, SHAN X, H. PAYER J 
[35] A nonlinear model for AC induced corrosion", Advanced Electromagnetics Symposium, AES, Paris, 04/2012 LE MENACH Yvonnick 
[36] Study of Synchronous Generator Static Eccentricities – FEM Results and Measurements ICEM 2012, Marseille, FRANCE, 01/2012 TANNOIAMAMURA Bruno, LE MENACH Yvonnick, TOUNZI Abdelmounaïm, SADOWSKI Nelson, GUILLOT Eilin 
[37] Eddy current effects on the steady state operating of a nonsalient pole synchronous machine The 15th International Symposium on Applied Electromagnetics and Mechanics, 09/2011, Abstract TANNOIAMAMURA Bruno, LE MENACH Yvonnick, TOUNZI Abdelmounaïm, SADOWSKI Nelson, GUILLOT Eilin 
This paper deals with the study of the effects of the eddy current induced in the
damper bars of a synchronous machine at steady state. A nonsalient pole synchronous
machine, with damper bars, is modelled in magnetostatic and magnetodynamic cases taking
into account the end rings. The time variation of the no load back emf is calculated in both
cases and the results are compared to measurements to highlight the effects of the induced
currents. 
[38] A TIMEDOMAIN IMPLICITSCHEMA DIRECT SOLVER: APPLICATION TO FINITE INTEGRATION SOLUTION
ISEF 2011, 09/2011, Abstract WANG Zifu, LE MENACH Yvonnick, TANG Zuqi, KORECKI Julien, HENNERON Thomas 
In timedomain electromagnetic fields computation, numerical methods (such as Finite Element
Method (FEM), Finite Integration Technique (FIT) [13] and etc.) have been applied. For the time
domain integration solution, explicit and implicit schemas have been widely used.
In comparison with implicit methods, the explicit methods are easier to realize in terms of
computation complexity, however, they are constrained by the stability condition. This condition
may require a small time step and therefore a prohibitive computing time. Another possibility is to
use an implicit schema which ensures the numerical stability. Unfortunately the implicit methods
require equation solved at each time step [4]. As a consequence, despite of a free choice on time
step, the computation time using the full implicit methods increases.
In this paper, fixedpoint explicit calculation is introduced to an implicit schema. This method
combines the two advantages of implicit and explicit methods: no stability condition and no
equation solving. The solver is then applied to a timedomain eddy current problem. Using
orthogonal mesh cells and FIT, the massmatrices in discrete formulations are diagonal. The fixed
point explicit method allows direct calculations without matrix inversion or decomposition.

[39] A Three Dimensional Electromagnetic Shell Finite Element for Coupled VectorScalar Potential
Formulations Compumag, Sydney, Australia, 07/2011 THOMAS Pierre, LE MENACH Yvonnick 
[40] Automatic MultiGPU Code Generation Applied to Simulation of Electrical Machines Compumag, Sydney, Australia, 07/2011 RODRIGUES A. Wendell de O., GUYOMARC’H Frédéric, LE MENACH Yvonnick, DEKEYSER JeanLuc 
[41] Proper adaptive method for Non intrusive Spectral Projection – Application on Eddy Current Non Destructive Testing COMPUMAG 2011, 07/2011 BEDDEK Karim, CLENET Stéphane, MOREAU Olivier, COSTAN Valentin, LE MENACH Yvonnick, BENABOU Abdelkader 
[42] Comparison of Residual and Equilibrated Error Estimators for FEM Applied to Magnetostatic Problems COMPUMAG 2011 Sydney, Australie, 07/2011, Abstract TANG Zuqi, LE MENACH Yvonnick, CREUSE Emmanuel, NICAISE Serge, PIRIOU Francis 
In finite element computations, the choice of the mesh is crucial to obtain an accurate solution. In order to evaluate
the quality of the mesh, a posteriori error estimators can be used. In this paper, we analyze and compare the residual and equilibrated error estimators for magnetostatic problems. 
[43] Direct and inverse methods for uncertainty quantification on Non destructive problem ACOMEN’2011, Liège (Belgique), 01/2011, 01/2011 BEDDEK Karim, DE STAELEN Rob, CLENET Stéphane, MOREAU Olivier, LE MENACH Yvonnick 
[44] Proper adaptive method for Non intrusive Spectral Projection – Application on Eddy Current Non Destructive Testing COMPUMAG 2011, 01/2011 BEDDEK Karim, CLENET Stéphane, MOREAU Olivier, COSTAN Valentin, LE MENACH Yvonnick 
[45] Polynomial chaos forward models in Bayesian inference to solve inverse problem ISEM’2011, Naple (Italie), 01/2011 LE MENACH Yvonnick 
[46] Study of Synchronous Generator eccentricities
using analytical approach and FEM ICEM 2010  International Conference on Electrical Machines, Rome, Italy, 09/2010, Abstract TANNOIAMAMURA Bruno, LE MENACH Yvonnick, TOUNZI Abdelmounaïm, SADOWSKI Nelson, GUILLOT Eilin 
This paper deals with the study of rotor turbogenerator
eccentricities using an analytical approach and FEM.
First, we present the cases of eccentricities and the method,
based on the measure of the flux density in the airgap, used to
detect the faults. Both models, the analytical approach and the
numerical one are then introduced. Using these models, static,
dynamic and combined eccentricities are studied at no load. The
results obtained by both approaches are compared and the
limits of the analytical model are underlined relatively to the
effect of the non linearity of the magnetic material
characteristic and the induced eddy currents. 
[47] Study of Interturn Short Circuit in Rotor Windings of a
Synchronous Generator Using FEM CEFC 2010  Conference on Electromagnetic Field Computation, Chicago, USA, 05/2010, Abstract TANNOIAMAMURA Bruno, LE MENACH Yvonnick, TOUNZI Abdelmounaïm, SADOWSKI Nelson, GUILLOT Eilin 
This paper deals with the study of interturn short
circuits in rotor windings of a synchronous generator using FEM.
Results obtained for magnetostatic and magnetodynamic cases
are presented and discussed. 
[48] 3D Stochastic Spectral Finite Element Method in static electromagnetism using vector potential formulation CEFC 2010, Chicago, USA, 05/2010, Abstract BEDDEK Karim, CLENET Stéphane, LE MENACH Yvonnick, MOREAU Olivier 
Stochastic Spectral Finite Element Method (SSFEM) can be used to take into account some random aspects in the input data (material characteristic, source terms) involved in static electromagnetism problems. Similarly to the deterministic case, two potential formulations can be used in the stochastic case. The vector potential formulation applied to static problems is developed and compared to the scalar potential one, previously developed.

[49] Parallel Sparse Matrix Solver on the GPU Applied to
Simulation of Electrical Machines Compumag 2009 Florianopolis Bresil, 11/2009, Abstract RODRIGUES A. Wendell de O., GUYOMARC’H Frédéric, LE MENACH Yvonnick, DEKEYSER JeanLuc 
The modeling and simulation of electrical systems use a very large
computational algorithms to solve their problems. The Conjugate
Gradient (CG) solver is an effective technique for symmetric positive
definite systems. It is suitable for systems of the form Ax=b, where A
is a known, square, symmetric, positivedefinite sparse matrix, x is a
unknown solution vector and b is a known vector. A and b are
obtained from the discretization of the FiniteElement Method (FEM). 
[50] Study of Static and Dynamic Eccentricities of a Synchronous Generator Using 3D FEM Compumag 2009, Florianopolis, Brésil, 11/2009, Abstract TANNOIAMAMURA Bruno, LE MENACH Yvonnick, TOUNZI Abdelmounaïm, SADOWSKI Nelson, GUILLOT Eilin 
The paper deals with the study of static and dynamic rotor eccentricities of a synchronous generator using 3D FEM. First, both eccentricity cases are introduced as well as the used approach. Then, the model of the structure, based on the vector potential formulation is presented. Results obtained for magnetostatic case related to both defects and their combination are presented and discussed. The saturation effect is highlighted. 
[51] Parallel Direct Solver For The Finite Integration Technique in Electrokinetic Problems Compumag 2009, florianopolis, brésil, 11/2009 TINZEFTE Abdellatif, LE MENACH Yvonnick, KORECKI Julien, PIRIOU Francis, GUYOMARC’H Frédéric 
[52] Parallel direct solver for the finite integration
TECHNIQUE IN ELECTROMAGNETICS ISEF 2009 Arras France, 10/2009, Abstract TINZEFTE Abdellatif, LE MENACH Yvonnick, KORECKI Julien, PIRIOU Francis 
The Finite Integration Technique is used in the different engineering domains as well as to computing the electromagnetic phenomena. This technique is efficient if the mesh is generated by
hexahedron elements. Moreover the matrix system, obtains from a regular mesh, can be exploit to use the parallel direct solver. In fact, in reordering the unknowns by the nested dissection method, it is possible to
construct directly the lower triangular matrix of the Cholesky factorization with many processors without
assembling the matrix system. In this paper, the parallel direct solver is described and the efficiency is shown as function to the number of processors. 
[53] Flux density curves variations using an EF model for turbogenerators diagnostics
ICEM, 09/2008 GUILLOT Eilin, LE MENACH Yvonnick, GUILLOT Yves 
[54] A Tree Construction Technique for the Feecad Electromagnetic 3D System Conf. MOMAG 2008, Florianopolis, Brésil, 09/2008, Abstract BASTOS João Pedro Assumpção, SADOWSKI Nelson, MÜLLER JulianaLuisa, TANNOIAMAMURA Bruno, LE MENACH Yvonnick 
In this paper we present a formulation to solve 3D magnetostatic field called “sourcefield” method. In this method the influence of the imposed coil current is replaced by the magnetic fields created by the coil on the edges of the finite element mesh. This technique is very robust and presents accurate and fast results, but the edge tree construction can be a delicate task. We present, firstly, the sourcefield method and then we propose a way intended to create the edges tree as well as the conductor facets tree, also necessary. 
[55] Analysis of dual formulations for eddy current problems with FIT EPNC 2008, 07/2008, Abstract KORECKI Julien, LE MENACH Yvonnick, PIRIOU Francis, DUCREUX JeanPierre 
This paper presents a method allowing to impose global electrics quantities, in a time step magnetodynamic case, simulated with the Finite Integration Technique. The A and TΩ formulations are developed and two vectors fields are associated in order to impose and compute global quantities, such as currents and voltages. Several simulations, applied to a system with two conductors, were carried out and the results were compared to those obtained with the finite element method. 
[56] Study of an electromagnetic gearbox involving two PM synchronous machines using 3DFEM Intermag’08, Madrid, Espagne, 05/2008, Abstract AUBERTIN Mathieu, TOUNZI Abdelmounaïm, LE MENACH Yvonnick 
Recently, magnetic gearboxes know a real interest and have been the subject of many studies due to their high performances. However, their principle is often issue from the mechanical gearbox one and then they do not offer any degree of freedom. In this paper, we propose a magnetic gearbox constituted of two electromagnetic structures. This introduces two degrees of freedom and then the system can also insure clutch operating. First, the operating of the system is described. Then,
its behavior is studied using a numerical model based on 3DFEM.
Different results given by 3D analysis are presented. 
[57] Iterative Solvers For Singular Symmetric Linear
Systems in Electromagnetism Low Frequency Conference CEFC 08 Athens, 05/2008, Abstract TINZEFTE Abdellatif, LE MENACH Yvonnick, PIRIOU Francis 
Application of different iterative methods (CG,MRTR, SQMR, MINRES) are proposed to solve linear systems from Finite Element Method (FEM) in 3D electromagnetic field analysis, and their convergences are compared. 
[58] Method to connect nonconforming mesh in 3D with the overlapping method CEFC 2008, 05/2008, Abstract KREBS Guillaume, CLENET Stéphane, LE MENACH Yvonnick, LOCMENT Fabrice, PIRIOU Francis 
In this paper, the overlapping method is applied to connect nonconforming tetrahedral meshes. At first the FEM formulation is described. The overlapping method is presented as well as its implementation in a FEM calculation code. At last, applications will be studied and discussed. 
[59] Implementation of a Vector Hysteresis Model in 2D Finite Element Analysis: Study of a RSST with Anisotropic Sample ISEM 2007, Michigan, USA, 09/2007 LEITE Jean Vianei, BENABOU Abdelkader, SADOWSKI Nelson, LE MENACH Yvonnick, CLENET Stéphane, DA SILVA Pedro Armando 
[60] A hybrid movement method to model electrical machines with end winding in 3D Finite Element Method CEFC, 06/2006, Abstract SHI Xiaodong, LE MENACH Yvonnick, ABAKAR Ali, DUCREUX JeanPierre, PIRIOU Francis 
In this paper, a hybrid method is proposed for the simulation of rotating electrical machines taking into account the end winding effect. The method is based on the sliding surface that is decomposed on two parts. In the air gap we use the interpolation method at the second order and in the area of the end winding the Surface Interpolation Method (first order). As application a turbogenerator is studied. 
[61] Study of the coupling between magnetic and electric equations
by using FEM, FIT and Cellmethod EPNC 2006, Maribor (Slovenia), 06/2006, URL, Abstract PARENT Guillaume, LE MENACH Yvonnick, DUCREUX JeanPierre, PIRIOU Francis 
Usually the Finite Element Method (F.E.M) is used to model the electromagnetic devices for low frequency. But other methods can be used as the CellMethod and the Finite Integration Technique (F.I.T.). The presented work is a comparative study about these methods. Their theories are reminded to show that the same source field can be used for all methods. The properties of the mass matrices generated though the three different methods are also discussed. Finally, as an example of application an iron core coil supplied by a voltage source is modelled. 
[62] Sparse Matrix Storage Formats and MatrixVector Products
Compared in Parallel Cluster Computing
EMF2006 ,aussois , france, 06/2006 FISCHBORN Marcos, KUOPENG Patrick, SADOWSKI Nelson, BASTOS João Pedro Assumpção, LE MENACH Yvonnick 
[63] Comparison of Methods to Simulate the Movement of Electrical Machines in 3D Using FEM COMPUMAG, Vol. 3, pages. 80, 06/2005, Abstract SHI Xiaodong, LE MENACH Yvonnick, DUCREUX JeanPierre, PIRIOU Francis 
In this paper three techniques ( the lockedstep approach, the overlapping method and the nodal interpolation method) are compared to model the rotation of electrical machines in 3D using FEM. These methods are associated with a formulation expressed in terms of scalar potential. To extend the possibilities of the interpolation method we propose the Surface Interpolation Method, which allows us to use nonconform meshes in the air gap of the studied machine. These different techniques are used to model a permanent magnet machine and a turbogenerator. The results obtained are compared by calculating the electromotive force and the torque. 
[64] COMPARISON OF DIFFERENT TECHNIQUES TO MODEL THE MOVEMENT IN 3D WITH FEM EPNC, 12/2004, Abstract SHI Xiaodong, LE MENACH Yvonnick, DUCREUX JeanPierre, PIRIOU Francis 
In this paper we compare three techniques to model the movement in 3D FEM using both classical formulations to solve the magnetic equations. To compare the methods a permanent magnet synchronous machine is modelled. The comparison is done at no load on the electromotive force and the cogging torque. 
[65] Error estimation for 3D magnetostatic problems coupled with electrical circuits procceding XVI EPNC, 09/2000 MARQUES Georges, LE MENACH Yvonnick, CLENET Stéphane, PIRIOU Francis 
[66] Stochastic Non Destructive Testing simulation : sensitivity analysis applied to material properties for Tube Support Plate clogging in nuclear power plant steam generator CEFC 2012 Oita, Japan, 11/1990 MOREAU Olivier, BEDDEK Karim, CLENET Stéphane, LE MENACH Yvonnick 
ACN Conférence nationale avec acte 
[1] Technique de génération de maillage par intersection de lignes équipotentielles numelec,Paris, 11/2017 DEBUSSCHÈRE Lydéric, LE MENACH Yvonnick, CLENET Stéphane, CHEVALLIER Loïc, GOURSAUD Benjamin, DUCREUX JeanPierre 
[2] Modélisation d’empilements de tôles par une méthode de couplage entre une formulation A* spectrale et une approche homogénéisée JCGE 2017 (Jeunes Chercheurs en Génie Electrique), 06/2017, Abstract ZIANI Smail, HENNERON Thomas, LE MENACH Yvonnick 
Nous proposons une approche de couplage entre une formulation éléments finis spectrale et une approche d’homogénéisation afin de modéliser un empilement de tôles. L’approche spectrale permet d’obtenir directement le régime permanent, l’homogénéisation permet de réduire le nombre d’inconnues et la formulation classique nous permet de prendre en compte les effets de bord des tôles. Notre approche prend en compte une loi de comportement magnétique non linéaire et une conductivité électrique dans les tôles. Les résultats de notre approche de couplage sont comparés aux résultats d’une approche classique de type éléments finis pas à pas dans le temps. 
[3] Réduction de modèle par POD appliquée à un transformateur planar en régime magnétoharmonique NUMELEC15  St Nazaire, 06/2015, Abstract TAYLOR Laurent, HENNERON Thomas, MARGUERON Xavier, LE MENACH Yvonnick, LE MOIGNE Philippe 
La prédétermination, dans la phase de dimensionnement, des pertes et des valeurs d’inductances est nécessaire à l’élaboration de composants magnétiques performants pour l’électronique de puissance. La modélisation numérique peut être utilisée pour déterminer les caractéristiques propres d’un composant notamment lorsque sa géométrie est "complexe" ou que la fréquence d’utilisation est élevée. Ces modèles sont très précis mais le temps de simulation nécessaire est important par rapport à des modèles analytiques plus simplistes. Les méthodes de réduction de modèles, type POD ou autres, permettent de diminuer le nombre de ces simulations numériques tout en conservant une grande précision. 
[4] Etude expérimentale et modélisation numérique magnétothermique de busbars laminés NUMELEC, 01/2015, Abstract PUIGDELLIVOL Oriol, LE MENACH Yvonnick, HARMAND Souad, MÉRESSE Damien, WECXSTEEN JeanFrançois 
Dans le monde de l’électronique de puissance il est impératif d’avoir une méthode de conception adaptée aux exigences et aux limitations du produit avec une optimisation sur les coûts et les prestations. Dans ce travail, on présente la simulation électromagnétique des busbars dans le but de calculer avec précision l’inductance et la densité de courant et ainsi de prévoir, par un couplage faible, la température de l’ensemble (conducteurs et isolants). 
[5] Modélisation 3D d’Inductances Planar pour l’Électronique de Puissance Symposium de Génie Électrique 2014, 07/2014, URL, Abstract TAYLOR Laurent, MARGUERON Xavier, LE MENACH Yvonnick, LE MOIGNE Philippe 
La prédétermination des valeurs de pertes et d’inductances est nécessaire à l’élaboration de composants magnétiques performants pour l’électronique de puissance. Leurs caractéristiques propres, variant en fonction de la fréquence, doivent être déterminées de manière précise sur une large gamme de fréquences, idéalement jusqu’à 30MHz pour tenir compte des aspects CEM. Dans la littérature, des simulations éléments finis 2D sont largement utilisées pour ce type de calcul. Ces modélisations servent généralement de référence pour le développement de modèles analytiques, notamment pour le cas spécifique des pertes cuivre. L’étude en 3D de ce type de composant est encore très peu développée car elle engendre des problèmes de modélisation, de ressources informatiques et de temps de calcul. Dans cet article la comparaison en termes de performances et de précisions pour des modèles d’inductances planar en 2D et en 3D est présentée. Les résultats obtenus sont comparés à des mesures d’impédances sur des prototypes. 
[6] Imposition d’un courant uniforme dans un conducteur Numelec 2012, Marseille, France, 07/2012, Abstract PIERQUIN Antoine, LE MENACH Yvonnick, ROGER JeanYves, CHEVALLIER Loïc 
Cette communication présente une méthode pour discrétiser
la densité de courant, en imposant une divergence nulle et en minimisant
l’erreur par rapport à la densité de courant exacte et uniforme. 
[7] Comparaison de la SSFEM et de la Projection Spectrale Non Intrusive à l’aide d’un estimateur d’erreur en électromagnétisme 03/2010 BEDDEK Karim, CLENET Stéphane, LE MENACH Yvonnick, MOREAU Olivier 
[8] Modélisation d’une matrice issue de la méthode des éléments finis NUMELEC, 12/2008 TAILLARD Julien, GUYOMARC’H Frédéric, LE MENACH Yvonnick, PIRIOU Francis, DEKEYSER JeanLuc 
AP Autre publication 
[1] La technique d’intégration Finie REE, 12/2006 LE MENACH Yvonnick, DUCREUX JeanPierre 
[2] Les techniques de prise en compte du mouvement des machines électriques tournantes dans la méthode des élément finis REE, 12/2006 LE MENACH Yvonnick, DUCREUX JeanPierre 
HDR Habilitation à diriger des recherches 
[1] Contribution à la modélisation numérique des phénomènes électromagnétiques 3D en basse fréquence Université de Lille, 12/2012, URL LE MENACH Yvonnick 
INV Conférence invité 
[1] Residualbased Error Estimator for Finite Element Method Applied to Quasi Static fields monag, Curitiba, 09/2014 LE MENACH Yvonnick, TANG Zuqi, TITTARELLI Roberta 
[2] Accounting of different movements in the study of electomechanical devices using 3D FEM MOMAG 2008, 09/2008 TOUNZI Abdelmounaïm, LE MENACH Yvonnick, HENNERON Thomas, KREBS Guillaume, CLENET Stéphane 
[3] Numerical solutions of electrokinetic and magnetostatic problems with imposed
global quantities using FEM and FIT MOMAG06 (BeloHorizonte  BRESIL), 07/2006 HENNERON Thomas, KORECKI Julien, LE MENACH Yvonnick, CLENET Stéphane, PIRIOU Francis 
Le L2EP recrute
Dernières actualités
 Thèse, Mario BERMÚDEZ GUEMÁN, 21 Déc. 2018
 Inauguration eCAMPUS, 20 Déc. 2018
 Thèse, Jérôme BUIRE, 14 Déc. 2018
 Thèse, Émile DEVILLERS, 13 Déc. 2018
 Séminiare chaire SBnodesSG, 7 Déc. 2018
 Séminaire, Pr Jiaxin YUAN (Chine), 6 Déc. 2018
 Inauguration projet PANDA, 4 déc. 2018
 L2EP as coordinator of a H2020 project
 Thèse : Bilel ZAIDI, 30 Nov. 2018
 Séminaire, Pr. Alain BOSSAVIT (Supelec), 29 Nov. 2018