Mohammad Issa, Jean-René Poirier, Ronan Perrussel, Olivier Chadebec and Victor Péron
Thin conducting sheets are used in many electric and electronic devices. Solving numerically the eddy current problems in presence of these thin conductive sheets requires a very…
Abstract
Purpose
Thin conducting sheets are used in many electric and electronic devices. Solving numerically the eddy current problems in presence of these thin conductive sheets requires a very fine mesh which leads to a large system of equations, and it becomes more problematic in case of higher frequencies. The purpose of this paper is to show the numerical pertinence of equivalent models for 3D eddy current problems with a conductive thin layer of small thickness e based on the replacement of the thin layer by its mid-surface with equivalent transmission conditions that satisfy the shielding purpose, and by using an efficient discretization using the boundary element method (BEM) to reduce the computational work.
Design/methodology/approach
These models are solved numerically using the BEM and some numerical experiments are performed to assess the accuracy of the proposed models. The results are validated by comparison with an analytical solution and a numerical solution by the commercial software Comsol.
Findings
The error between the equivalent models and analytical and numerical solutions confirms the theoretical approach. In addition to this accuracy, the computational work is reduced by considering a discretization method that requires only a surface mesh.
Originality/value
Based on a hybrid formulation, the authors present briefly a formal derivation of impedance transmission conditions for 3D thin layers in eddy current problems where non-conductive materials are considered in the interior and the exterior domain of the sheet. BEM is adopted to discretize the problem as there is no need for volume discretization.
Details
Keywords
Christian Vollaire, François Musy and Ronan Perrussel
Propose post processing methods for the edge finite element (FE) method on a tetrahedral mesh. They make it possible to deduce vector values on the vertices from scalar values…
Abstract
Purpose
Propose post processing methods for the edge finite element (FE) method on a tetrahedral mesh. They make it possible to deduce vector values on the vertices from scalar values defined on the edges of the tetrahedra.
Design/methodology/approach
The new proposed techniques are based on a least squares formulation leading to a sparse matrix system to be solved. They are compared in terms of accuracy and CPU time on a FEs formulation for open boundary – frequency domain problems.
Findings
A significant improvement of vector values accuracy on the vertices of the tetrahedra is obtained compared to a classical approach with a very small additional computation time.
Originality/value
This work presents techniques: to obtain the values at the initial nodes of the mesh and not inside the tetrahedra; and to take into account the discontinuity to the interface between two media of different electromagnetic properties.
Details
Keywords
Riccardo Scorretti, Ronan Perrussel, Laurent Morel, Noël Burais and Laurent Nicolas
The classical ϕ‐a formulations for numerical dosimetry of currents induced by extremely low frequency magnetic fields requires that the source field is provided through a vector…
Abstract
Purpose
The classical ϕ‐a formulations for numerical dosimetry of currents induced by extremely low frequency magnetic fields requires that the source field is provided through a vector potential. The purpose of this paper is to present a new formulation t‐b which directly takes the flux density as source term.
Design/methodology/approach
This formulation is implemented through finite element and validated by comparison with analytical solutions. The results obtained by both formulations are compared in the case of an anatomical computational phantom exposed to a vertical uniform field.
Findings
A good agreement between the t‐b formulation and both numerical and analytical computations was found.
Originality/value
This new formulation seems to be more accurate than the ϕ‐a formulation, and is more suited for situations where the magnetic field is known from experimental measurements, as there is no need for a magnetic vector potential.
Details
Keywords
Liang Li, Stéphane Lanteri and Ronan Perrussel
This work is concerned with the development and the numerical investigation of a hybridizable discontinuous Galerkin (HDG) method for the simulation of two‐dimensional…
Abstract
Purpose
This work is concerned with the development and the numerical investigation of a hybridizable discontinuous Galerkin (HDG) method for the simulation of two‐dimensional time‐harmonic electromagnetic wave propagation problems.
Design/methodology/approach
The proposed HDG method for the discretization of the two‐dimensional transverse magnetic Maxwell equations relies on an arbitrary high order nodal interpolation of the electromagnetic field components and is formulated on triangular meshes. In the HDG method, an additional hybrid variable is introduced on the faces of the elements, with which the element‐wise (local) solutions can be defined. A so‐called conservativity condition is imposed on the numerical flux, which can be defined in terms of the hybrid variable, at the interface between neighbouring elements. The linear system of equations for the unknowns associated with the hybrid variable is solved here using a multifrontal sparse LU method. The formulation is given, and the relationship between the considered HDG method and a standard upwind flux‐based DG method is also examined.
Findings
The approximate solutions for both electric and magnetic fields converge with the optimal order of p+1 in L2 norm, when the interpolation order on every element and every interface is p and the sought solution is sufficiently regular. The presented numerical results show the effectiveness of the proposed HDG method, especially when compared with a classical upwind flux‐based DG method.
Originality/value
The work described here is a demonstration of the viability of a HDG formulation for solving the time‐harmonic Maxwell equations through a detailed numerical assessment of accuracy properties and computational performances.