A novel 3-D hybrid formulation for time-harmonic eddy current problems in multiply connected domains is presented. The interior problem (in conductive regions) is discretized by the cell method (CM) in terms of magnetic vector potentials a, whereas the exterior problem (in the unbounded air domain) is discretized by the boundary element method (BEM) in terms of reduced scalar potentials {arphi }-{r}. Novel topological constraints are derived from magnetostatic energy conservation by using a decomposition of the magnetic field which minimizes the support and the number of cohomology generators. A fast algebraic procedure to pre-compute source fields to handle efficiently current-driven coils is also proposed. The final matrix system can be solved in a limited number of iterations by transpose-free quasi-minimal residual method with symmetric successive over-relaxation preconditioning. Convergence tests show that numerical results are in a very good agreement with third-order finite element method (FEM) on a 2-D axisymmetric model. The CM-BEM with piecewise uniform approximation shows also to be very effective when analyzing fully 3-D test cases, since second-order FEM accuracy is attained even with coarse mesh refinements, by using, however, a much lower number of degrees of freedom compared to FEM.
A 3-D Hybrid Cell Boundary Element Method for Time-Harmonic Eddy Current Problems on Multiply Connected Domains
Moro, Federico
;
2019
Abstract
A novel 3-D hybrid formulation for time-harmonic eddy current problems in multiply connected domains is presented. The interior problem (in conductive regions) is discretized by the cell method (CM) in terms of magnetic vector potentials a, whereas the exterior problem (in the unbounded air domain) is discretized by the boundary element method (BEM) in terms of reduced scalar potentials {arphi }-{r}. Novel topological constraints are derived from magnetostatic energy conservation by using a decomposition of the magnetic field which minimizes the support and the number of cohomology generators. A fast algebraic procedure to pre-compute source fields to handle efficiently current-driven coils is also proposed. The final matrix system can be solved in a limited number of iterations by transpose-free quasi-minimal residual method with symmetric successive over-relaxation preconditioning. Convergence tests show that numerical results are in a very good agreement with third-order finite element method (FEM) on a 2-D axisymmetric model. The CM-BEM with piecewise uniform approximation shows also to be very effective when analyzing fully 3-D test cases, since second-order FEM accuracy is attained even with coarse mesh refinements, by using, however, a much lower number of degrees of freedom compared to FEM.Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
