An optimum PML for scattering problems in the time domainModave, Axel ; ; et alin European Physical Journal : Applied physics (in press) In electromagnetic compatibility, scattering problems are defined in an infinite spatial domain, while numerical techniques such as finite element methods require a computational domain that is bounded ... [more ▼] In electromagnetic compatibility, scattering problems are defined in an infinite spatial domain, while numerical techniques such as finite element methods require a computational domain that is bounded. The Perfectly Matched Layer (PML) is widely used to simulate the truncation of the computational domain. However, its performance depends critically on an absorption function. This function is generally tuned by using case-dependent optimization procedures. In this paper, we will present some efficient functions that overcome any tuning. They will be compared using a realistic scattering benchmark solved with the Discontinuous Galerkin method. [less ▲] Detailed reference viewed: 7 (4 ULg) A Frontal Delaunay Quad Mesh Generator Using the L ∞ Norm; ; et al in International Journal for Numerical Methods in Engineering (in press) Detailed reference viewed: 60 (18 ULg) A multirate approach to combine electromagnetic transients and fundamental-frequency simulationsPlumier, Frédéric ; Geuzaine, Christophe ; Van Cutsem, Thierry ![]() Conference (2013, July) Abstract--- A new multi-rate technique is proposed to combine fundamental-frequency simulation (typical of stability studies) with electromagnetic transients simulation. The objective of this hybrid ... [more ▼] Abstract--- A new multi-rate technique is proposed to combine fundamental-frequency simulation (typical of stability studies) with electromagnetic transients simulation. The objective of this hybrid approach is to obtain more accurate simulations than with the fundamental-frequency approximation, while saving computing time by applying the detailed model to a subsystem only, in some neighbourhood of the disturbance. It also allows to remove some limitations of fundamental-frequency simulation, such as the difficulty of simulating unbalanced faults. A relaxation technique is used to iterate between both models with simple interfacing. Preliminary results obtained with a 74-bus test system are presented, together with a comparison with full electromagnetic transients simulation. [less ▲] Detailed reference viewed: 29 (4 ULg) Fifth International Conference on Advanced COmputational Methods in ENgineering (ACOMEN 2011)Béchet, Eric ; ; Geuzaine, Christophe et alin International Journal of Computational & Applied Mathematics (2013), 246 Detailed reference viewed: 23 (9 ULg) An Energy-Based Variational Model of Ferromagnetic Hysteresis for Finite Element Computations; ; Stainier, Laurent et alin Journal of Computational & Applied Mathematics (2013), 246 This paper proposes a macroscopic model for ferromagnetic hysteresis that is well-suited for finite element implementation. The model is readily vectorial and relies on a consistent thermodynamic ... [more ▼] This paper proposes a macroscopic model for ferromagnetic hysteresis that is well-suited for finite element implementation. The model is readily vectorial and relies on a consistent thermodynamic formulation. In particular, the stored magnetic energy and the dissipated energy are known at all times, and not solely after the completion of closed hysteresis loops as is usually the case. The obtained incremental formulation is variationally consistent, i.e., all internal variables follow from the minimization of a thermodynamic potential. [less ▲] Detailed reference viewed: 86 (22 ULg) A PML for convex truncated domains in time-dependent acoustics with a DG-FE discretizationModave, Axel ; ; Delhez, Eric et alin Proceedings of the 11th International Conference on Mathematical and Numerical Aspects of Waves (Waves 2013) (2013, June 03) A new perfectly matched layer (PML) is proposed for convex truncated domains in the context of time-dependent acoustics. With this formulation, the size of the computational domain can be reduced when ... [more ▼] A new perfectly matched layer (PML) is proposed for convex truncated domains in the context of time-dependent acoustics. With this formulation, the size of the computational domain can be reduced when classical shapes of truncated domains are less appropriate. A numerical discretization based on the discontinuous Galerkin method is then described and validated. An example of realistic three-dimensional application is finally proposed. [less ▲] Detailed reference viewed: 14 (4 ULg) Computational Homogenization for Laminated Ferromagnetic Cores in MagnetodynamicsNiyonzima, Innocent ; ; Dular, Patrick et alin IEEE Transactions on Magnetics (2013), 49(5), 2049-2052 In this paper, we investigate the modeling of ferromagnetic multiscale materials. We propose a computational homogenization technique based on the heterogeneous multiscale method (HMM) that includes both ... [more ▼] In this paper, we investigate the modeling of ferromagnetic multiscale materials. We propose a computational homogenization technique based on the heterogeneous multiscale method (HMM) that includes both eddy-current and hysteretic losses at the mesoscale. The HMM comprises: 1) a macroscale problem that captures the slow variations of the overall solution; 2) many mesoscale problems that allow to determine the constitutive law at the macroscale. As application example, a laminated iron core is considered. [less ▲] Detailed reference viewed: 37 (8 ULg)![]() Computing Bounds on the Geometrical Quality of 2D Curvilinear Finite ElementsJohnen, Amaury ; ; et alConference (2013, April 25) The development of high-order computational methods for solving partial differen- tial equations on unstructured grids has been underway for many years. Such meth- ods critically depend on the ... [more ▼] The development of high-order computational methods for solving partial differen- tial equations on unstructured grids has been underway for many years. Such meth- ods critically depend on the availability of high-quality curvilinear meshes, as one badly-shaped element can degrade the solution in the whole domain (J. Shewchuk, “What Is a Good Linear Finite Element? Interpolation, Conditioning, Anisotropy, and Quality Measures”, Preprint, 2002). The usual way of generating curved meshes is to first generate a straight sided mesh and to curve mesh entities that are classified on the boundaries of the domain. The latter operation introduces a “shape-distortion” that should be controlled if we suppose that the straight sided mesh is composed of well-shaped elements. Quality measures allow to quantify to which point an element is well-shaped. They also provide tools to improve the quality of meshes through optimization opera- tions. Many quality measures has been proposed for quadratic triangular finite element. Recently, X. Roca et al. (“Defining Quality Measures for High-Order Planar Triangles and Curved Mesh Generation”, Proceedings of the 20th Interna- tional Meshing Roundtable, 2011) proposed a technique that allows extending any Jacobian based quality measure for linear elements to high-order iso-parametric planar triangles of any interpolation degree. In this work we propose an efficient method to provide accurate bounds on the mag- nitude of the shape distortion of any triangular and quadrangular curved element. The shape distortion is measured with respect to an ideal element, which can e.g. be an equilateral triangle or the element from the original straight-sided mesh. The key feature of the method is that we can adaptively expand functions based on the Jacobian matrix and its determinant in terms of Be ́zier functions. Be ́zier functions have both properties of boundedness and positivity, which allow sharp computation of minimum or maximum of the interpolated functions. [less ▲] Detailed reference viewed: 19 (3 ULg) Evaluation of Shielding Effectiveness in the Time Domain using a DG Method with an Efficient PMLModave, Axel ; Geuzaine, Christophe ; et alPoster (2013, April 24) Development of new electric and electronic devices leads to an increasing importance for solving electromagnetic compatibility problems. In particular, the optimization of radiated emissions is ... [more ▼] Development of new electric and electronic devices leads to an increasing importance for solving electromagnetic compatibility problems. In particular, the optimization of radiated emissions is increasingly done through numerical modeling. In this context, numerical methods allow to evaluate the effectiveness of more and more elaborate shielding enclosures. We propose an application of the nodal Discontinuous Galerkin method to evaluate the effectiveness of composite shieldings in the time domain. This numerical method provides convenient advantages to deal with shielding problems. Firstly, heterogeneous media are naturally included in the discrete formulation. After, the method is famous to allow a massive parallelization for the numerical resolution. Finally, a resolution in the time domain allows to consider transient signals. The numerical scheme is firstly validated using an academic benchmark. In the same way, the PML used to truncate the computational domain is optimized and validated. Applications are then presented for homogeneous and composite shieldings. The preliminary results are in conformity with classical predictive rules. [less ▲] Detailed reference viewed: 25 (4 ULg) Subproblem h-Conform Formulation for Accurate Thin Shell Models Between Conducting and Nonconducting RegionsDang, Quoc Vuong ; Dular, Patrick ; Vazquez Sabariego, Ruth et alin Proceeding of the 9th International Symposium on Electric and Magnetic Fields, EMF 2013 (2013, April 23) A subproblem method (SPM) with h-formulation is developed for correcting the inaccuracies near edges and corners that arise from using thin shell (TS) models to replace thin volume regions by surfaces ... [more ▼] A subproblem method (SPM) with h-formulation is developed for correcting the inaccuracies near edges and corners that arise from using thin shell (TS) models to replace thin volume regions by surfaces. The developed surface-to-volume correction problem is defined as a step of multiple SPs, with inductors and magnetic or conducting regions, some of them being thin. The TS model assumes that the fields in the thin regions are approximated by a priori 1-D analytical distributions along the shell thickness (C. Geuzaine et al., “Dual formulations for the modeling of thin electromagnetic shells using edge elements,” IEEE Trans. Magn., vol. 36, no. 4, pp. 799–802, 2000). Their interior is not meshed and ratherextracted from the studied domain, which is reduced to a zero-thickness double layer with interface conditions (ICs) linked to 1-D analytical distributions that however neglect end and curvature effects. This leads to inaccuracies near edges and corners that increase with the thickness. To cope with these difficulties, the authors have recently proposed a SPM based on the h-formulation for a thin region located between non-conducting regions (Vuong Q. Dang et al., “Subproblem Approach for Thin Shell Dual Finite Element Formulations”, IEEE Trans. Magn., vol. 48, no. 2, pp. 407–410, 2012). The magnetic field h is herein defined in nonconducting regions by means of a magnetic scalar potential , i.e. h = -grad{\phi} , with discontinuities of through the TS. In this paper, the SPM is extended to account for thin regions located between conducting regions or between conducting and nonconducting regions, in the general case of multiply connected regions. In these regions, the potential is not defined anymore on both sides of the TS and the problem has to be expressed in terms of the discontinuities of h, possibly involving on one side only, to be strongly defined via an IC through the TS. In the proposed SP strategy, a reduced problem with only inductors is first solved on a simplified mesh without thin and volume regions. Its solution gives surface sources (SSs) as ICs for added TS regions, and volume sources (VSs) for possible added volume regions. The TS solution is further improved by a volume correction via SSs and VSs that overcome the TS assumptions, respectively suppressing the TS model and adding the volume model. Each SP has its own separate mesh, which increases the computational efficiency. Details on the proposed method will be given in the extended paper, with practical applications. [less ▲] Detailed reference viewed: 26 (3 ULg) A Computational Homogenization Method for the Evaluation of Eddy Current in Nonlinear Soft Magnetic CompositesNiyonzima, Innocent ; Vazquez Sabariego, Ruth ; Dular, Patrick et alin Proceeding of the 9th International Symposium on Electric and Magnetic Fields, EMF 2013 (2013, April 23) The use of the soft magnetic composite (SMC) in electric devices has increased in recent years. These materials made from a metallic powder compacted with a dielectric binder are a good alternative to ... [more ▼] The use of the soft magnetic composite (SMC) in electric devices has increased in recent years. These materials made from a metallic powder compacted with a dielectric binder are a good alternative to laminated ferromagnetic structures as their granular mesoscale structure allows to significantly reduce the eddy current losses. Furthermore unlike the laminated ferromagnetic structures, SMC exhibit isotropic magnetic properties what makes them good candidates for manufacturing machines with 3D flux paths. The isotropy of the thermal conductivity also allows for a more efficient heat dissipation. The use of classical numerical methods such as the finite element method to study the behavior of SMC is computational very expensive. Indeed a very fine mesh would be required in order to capture fine scale variations i.e. variations at level of metallic grains whence the use of multiscale methods for modelling SMC. The application of multiscale method to study the behaviour of SMC is relatively recent. In (A. Bordianu et al “A Multiscale Approach to Predict Classical Losses in Soft Magnetic Composites”, IEEE Trans. Mag., vol. 48, no. 4, 2012.), the authors used a homogenization technique to compute electrical and magnetic constitutive laws on a representative volume element (RVE). These laws were then used in finite element computations. Herein, the RVE has been chosen to account for the grain- grain contact that can occur in a actual SMC structure due to the compaction process and that can lead to the appearance of macroscale eddy currents. In this paper, we will extend the computational homogenization method success- fully used for modelling the behaviour of laminated ferromagnetic cores in mag- netodynamics (I. Niyonzima et al “Computational Homogenization for Laminated Ferromagnetic Cores in Magnetodynamics”, in Proc. of the 15th Biennal Confer- ence on Electromagnetic Field Computation, 2012) to the case of SMC. The method is based on the heterogeneous multiscale method (HMM) and couples two types of problems: a macroscale problem that captures the slow variations of the overall so- lution and many microscale problems that allow to determine the constitutive laws at the macroscale. The choice of RVE will also be discussed. [less ▲] Detailed reference viewed: 50 (10 ULg) Dual Formulations for Accurate Thin Shell Models in a Finite Element Subproblem MethodDang, Quoc Vuong ; Dular, Patrick ; Vazquez Sabariego, Ruth et alin Proceeding of the 19th COMPUMAG Conference on the Computation of Electromagnetic Fields, 2013 (2013, April 01) A subproblem finite with dual finite element magnetostatic and magnetodynamic formulations is developed for correcting the inaccuracies near edges and corners coming from thin shell models, that replace ... [more ▼] A subproblem finite with dual finite element magnetostatic and magnetodynamic formulations is developed for correcting the inaccuracies near edges and corners coming from thin shell models, that replace thin volume regions by surfaces. The surface-to-volume correction problem is defined as one of the multiple subproblems applied to a complete problem, considering successive additions of inductors and magnetic or conducting regions, some of these being thin regions. Each SP requires a proper adapted mesh of its regions, which facilitates meshing and increases computational e ciency. [less ▲] Detailed reference viewed: 21 (3 ULg) Power system dynamic simulation: an iterative multirate approachPlumier, Frédéric ; Geuzaine, Christophe ; Van Cutsem, Thierry ![]() Poster (2013, April) Detailed reference viewed: 22 (2 ULg) Power System Simulation ChallengesAristidou, Petros ; Plumier, Frédéric ; Van Cutsem, Thierry et alScientific conference (2013, March 05) Detailed reference viewed: 21 (2 ULg) Geometrical Validity of Curvilinear Finite ElementsJohnen, Amaury ; ; Geuzaine, Christophe ![]() in Journal of Computational Physics (2013), 233 In this paper, we describe a way to compute accurate bounds on Jacobian de- terminants of curvilinear polynomial finite elements. Our condition enables to guarantee that an element is geometrically valid ... [more ▼] In this paper, we describe a way to compute accurate bounds on Jacobian de- terminants of curvilinear polynomial finite elements. Our condition enables to guarantee that an element is geometrically valid, i.e., that its Jacobian determinant is strictly positive everywhere in its reference domain. It also provides an efficient way to measure the distortion of curvilinear elements. The key feature of the method is to expand the Jacobian determinant using a polynomial basis, built using B ́ezier functions, that has both properties of boundedness and positivity. Numerical results show the sharpness of our estimates. [less ▲] Detailed reference viewed: 69 (20 ULg) Geometrical Validity of High-Order Triangular Finite ElementsJohnen, Amaury ; ; Geuzaine, Christophe ![]() in Engineering with Computers (2012) This paper presents a method to compute accurate bounds on Jacobian determinants of high-order (curvilinear) triangular nite elements. This method can be used to guarantee that a curvilinear triangle is ... [more ▼] This paper presents a method to compute accurate bounds on Jacobian determinants of high-order (curvilinear) triangular nite elements. This method can be used to guarantee that a curvilinear triangle is geometrically valid, i.e., that its Jacobian determinant is strictly positive everywhere in its reference domain. It also provides an e cient way to measure the quality the triangles. The key feature of the method is to expand the Jacobian determinant using a polynomial basis, built using B ezier functions, that has both properties of boundedness and positivity. Numerical results show the sharpness of our estimates. [less ▲] Detailed reference viewed: 68 (14 ULg) Computational Homogenization for Laminated Ferromagnetic Cores in MagnetodynamicsNiyonzima, Innocent ; Vazquez Sabariego, Ruth ; et alin Proceedings of the 15th Biennial IEEE Conference on Electromagnetic Field Computation (CEFC2012) (2012, November) In this paper, we investigate the modeling of fer- romagnetic multiscale materials. We propose a computational homogenization method based on the heterogeneous multiscale method (HMM) that includes eddy ... [more ▼] In this paper, we investigate the modeling of fer- romagnetic multiscale materials. We propose a computational homogenization method based on the heterogeneous multiscale method (HMM) that includes eddy currents and hysteretic losses at the mesoscale. The HMM comprises: 1) a macroscale problem that captures the slow variations of the overall solution; 2) many microscale problems that allow to determine the constitutive law at the macroscale. As application example, a laminated iron core is considered. [less ▲] Detailed reference viewed: 72 (15 ULg) Homogenisation of Windings and Laminations in Time-Domain Finite-Element Modeling of Electrical Machines; Geuzaine, Christophe ; Vazquez Sabariego, Ruth ![]() in Proceedings of the 15th Biennial IEEE Conference on Electromagnetic Field Computation (CEFC2012) (2012, November) This paper deals with time-domain homogenisation of multi-turn windings and laminated cores in 2D and 3D finite- element (FE) modelling of rotating electrical machines. Herein the number of additional ... [more ▼] This paper deals with time-domain homogenisation of multi-turn windings and laminated cores in 2D and 3D finite- element (FE) modelling of rotating electrical machines. Herein the number of additional degrees of freedom (auxiliary field variables) in the homogenised regions can be fixed depending on the extent of the eddy current effects and on the desired accuracy. The homogenisation technique is illustrated and validated by means of a 2D model of a switched reluctance motor. Global quantities, such as the phase impedance, converge very well to those produced by a precise but very expensive 3D model. [less ▲] Detailed reference viewed: 37 (4 ULg) Time-domain finite-element modelling of laminated iron cores – Large skin effect homogenization considering the Jiles-Atherton hysteresis modelVazquez Sabariego, Ruth ; Niyonzima, Innocent ; Geuzaine, Christophe et alin Proceedings of the 15th Biennial IEEE Conference on Electromagnetic Field Computation (CEFC2012) (2012, November) This paper deals with the incorporation of the Jiles- Atherton (J-A) hysteresis model in a time-domain finite-element homogenization technique for laminated iron cores. The separate discretization of each ... [more ▼] This paper deals with the incorporation of the Jiles- Atherton (J-A) hysteresis model in a time-domain finite-element homogenization technique for laminated iron cores. The separate discretization of each lamination is avoided by using dedicated skin-effect basis functions, which also serve to interpolate the J- A hysteretic material law. As validation test case, a stacked ring core surrounded by a toroidal coil is considered. [less ▲] Detailed reference viewed: 44 (6 ULg) An optimum PML for scattering problems in the time domainModave, Axel ; Delhez, Eric ; et alin Proceedings of the 7th European Conference on Numerical Methods in Electromagnetism (NUMELEC2012) (2012, July 05) The Perfectly Matched Layer (PML) is widely used for unbounded problems. However its performances depend critically on an absorption coefficient. This coefficient is generally tuned by using case ... [more ▼] The Perfectly Matched Layer (PML) is widely used for unbounded problems. However its performances depend critically on an absorption coefficient. This coefficient is generally tuned by using case-dependent optimization procedures. In this paper we will present some efficient profiles of the coefficient that overcome every tuning. These profiles and others will be compared by using a scattering benchmark solved with the Discontinuous Galerkin method. [less ▲] Detailed reference viewed: 48 (14 ULg) |
||