Search results

Rotational‐translational addition theorems for the vector spheroidal wave functions Ma(i)mn(h; ξ, η, φ) and Na(i)mn(h; ξ, η, φ), i = 1,2,3,4, are derived from those for the corresponding scalar spheroidal wave functions ψ(i)mn(h; ξ, η, φ). A vector spheroidal wave function defined in one spheroidal coordinate system (h; ξ, η, φ) is expressed in terms of a series of vector spheroidal wave functions defined in another spheroidal coordinate system (h′; ξ′, η′, φ′), which is rotated and translated with respect to the first one. These theorems allow a rigorous treatment of boundary value problems relative to time‐harmonic vector field waves in the presence of a system of spheroids with arbitrary orientations. As a special case, general rotational‐translational addition theorems for vector spherical wave functions are also presented.

The finite‐element method can be used for an approximate solution of axisymmetric exterior‐field problems by truncating the unbounded domain, or by applying various techniques of coupling a finite region of interest with the remaining far region, which is properly modelled. In this paper, we propose the solution of axisymmetric exterior‐field problems by using the standard finite‐element method in a bounded, transformed domain obtained by conformal mapping from the original, unbounded one. The transformed functionals have very simple expressions and the exact transforms of the original boundary conditions are used in the transformed domain. Consequently no approximation is introduced in the proposed method and improvements in the accuracy of the solution are obtained as compared with several other methods in common usage, especially with the truncated mesh technique. A few example problems are solved and the presented method is found to be simple and computationally highly efficient. It is particularly recommended for problems with material inhomogeneities and anisotropies within large regions.

A simple and efficient method for the finite‐element solution of three‐dimensional unbounded region field problems is presented in this paper. The proposed technique consists of a global mapping of the original unbounded region onto a bounded domain by applying a standard inversion transformation to the spatial coordinates. Same numerical values of the potential function are assigned to the transformed points. The functional associated to the field problem, which incorporates the boundary conditions, has the same structure in the transformed domain as that in the original one. This allows the implementation of the standard finite‐element method in the bounded transformed domain. The finite‐element solution is obtained on the basis of a complete discretization of the bounded, transformed domain by standard finite elements, with no approximate assumption made for the behaviour of the field at infinity, other than that introduced by the finite‐element idealization. This leads to improved accuracy of the numerical results, compared to those obtained in the original region, for the same number of nodes. Application to three test problems illustrates the high efficiency of the proposed method in terms of both accuracy and computational effort. The technique presented is particularly recommended for exterior‐field problems in the presence of material inhomogeneities and anisotropies.

An exact mathematical solution has been obtained for the quasistationary electromagnetic field of a circular current loop coaxial with a conducting circular cylinder in uniform relative motion with respect to each other. The vector magnetic potential corresponding to a filamentary loop carrying a variable with time current is determined by applying the method of separation of variables. Expressions for the fields outside and inside the cylindrical conductor, the eddy‐current distribution, the power loss and the interaction force are derived directly from the vector potential. The fields due to practical current coils are obtained by integration from the results for the filamentary loop. An approximate simple formula is presented for a loop carrying direct current at high velocities. The analysis performed is relevant to the design and operation of magnetic devices with metallic cores in motion, linear electrical machines and electromagnetic launching systems.

The purpose of this paper is to present three novel techniques aimed at increasing the efficiency of the polarization fixed point method for the solution of nonlinear periodic field problems.

Firstly, the characteristic B‐M resulting from the constitutive relation B‐H is replaced by a relation between the components of the harmonics of the vectors B and M. Secondly, a dynamic overrelaxation method is implemented for the convergence acceleration of the iterative process involved. Thirdly, a modified dynamic overrelaxation method is proposed, where only the relation B‐M between the magnitudes of the field vectors is used.

By approximating the actual characteristic B‐M by the relation between the components of the harmonics of the vectors B and M, the amount of computation required for the field analysis is substantially reduced. The rate of convergence of the iterative process is increased by implementing the proposed dynamic overrelaxation technique, with the convergence being further accelerated by applying the modified dynamic overrelaxation presented. The memory space is also well reduced with respect to existent methods and accurate results for nonlinear fields in a real world structure are obtained utilizing a normal size processor notebook in a time of about one‐half of one minute when no induced currents are considered and of about one minute when eddy currents induced in solid ferromagnetic parts are also fully analyzed.

The originality of the novel techniques presented in the paper consists in the drastic approximations proposed for the material characteristics of the nonlinear ferromagnetic media in the analysis of periodic electromagnetic fields. These techniques are highly efficient and yield accurate numerical results.

The purpose of this paper is to implement the Anderson acceleration for different formulations of eletromagnetic nonlinear problems and analyze the method efficiency and strategies to obtain a fast convergence.

The paper is structured as follows: the general class of fixed point nonlinear problems is shown at first, highlighting the requirements for convergence. The acceleration method is then shown with the associated pseudo-code. Finally, the algorithm is tested on different formulations (finite element, finite element/boundary element) and material properties (nonlinear iron, hysteresis models for laminates). The results in terms of convergence and iterations required are compared to the non-accelerated case.

The Anderson acceleration provides accelerations up to 75 per cent in the test cases that have been analyzed. For the hysteresis test case, a restart technique is proven to be helpful in analogy to the restarted GMRES technique.

The acceleration that has been suggested in this paper is rarely adopted for the electromagnetic case (it is normally adopted in the electronic simulation case). The procedure is general and works with different magneto-quasi static formulations as shown in the paper. The obtained accelerations allow to reduce the number of iterations required up to 75 per cent in the benchmark cases. The method is also a good candidate in the hysteresis case, where normally the fixed point schemes are preferred to the Newton ones.

A transformer model is used as a benchmark for testing various methods to solve 3D nonlinear periodic eddy current problems. This paper aims to set up a nonlinear magnetic circuit problem to assess the solving procedure of the nonlinear equation system for determining the influence of various special techniques on the convergence of nonlinear iterations and hence the computational time.

Using the T,ϕ-ϕ formulation and the harmonic balance fixed-point approach, two techniques are investigated: the so-called “separate method” and the “combined method” for solving the equation system. When using the finite element method (FEM), the elapsed time for solving a problem is dominated by the conjugate gradient (CG) iteration process. The motivation for treating the equations of the voltage excitations separately from the rest of the equation system is to achieve a better-conditioned matrix system to determine the field quantities and hence a faster convergence of the CG process.

In fact, both methods are suitable for nonlinear computation, and for comparing the final results, the methods are equally good. Applying the combined method, the number of iterations to be executed to achieve a meaningful result is considerably less than using the separated method.

To facilitate a quick analysis, a simplified magnetic circuit model of the 3D problem was generated to assess how the different ways of solutions will affect the full 3D solving process. This investigation of a simple magnetic circuit problem to evaluate the benefits of computational methods provides the basis for considering this formulation in a 3D-FEM code for further investigation.

We present a set of new simple closed‐form analytical formulas for the calculation of the magnetic field produced at any point of space by any solid polyhedral conductor with a uniform current density j. The formulas have been obtained by analytical integration of Ampère's law under the only assumptions that the conductor is bounded by flat surfaces and that j = constant in the conductor. This includes bars, bricks, tetrahedrons, wedges, prisms, trapezoids, pyramids, and polyhedrons in general. The formulas contain no singularities, and can be used for the numerical calculation of the field at any point, including points inside the conductor, or on its surface, edges or corners. The formulas can easily be extended for conductors of infinite length. Extensive numerical tests of the formulas have been performed.

This paper aims to introduce the decomposed harmonic balance finite element method (HBFEM) to decrease the memory requirement in large‐scale computation of the DC‐biasing magnetic field. Harmonic analysis of the flux density and flux distribution was carried out to investigate the DC biased problem in a laminated core model (LCM).

Based on the DC bias test on a LCM, the decomposed HBFEM is applied to accurately calculate the DC‐biasing magnetic field. External electric circuits are coupled with the magnetic field in the harmonic domain. The reluctivity matrix is decomposed and the block Gauss‐Seidel algorithm solves each harmonic solution of magnetic field and exciting current sequentially.

The calculated exciting currents and flux density are compared with that obtained from measurement and time domain finite element analysis, respectively, which demonstrates consistency. The DC bias leads to the significant saturation of the magnetic core and serious distortion of the exciting current. The flux density varies nonlinearly with DC bias excitation.

The harmonic balance method is only applicable in solving the steady state magnetic field. Future improvements in the method are necessary in order to manage the hysteresis effects in magnetic material.

The proposed method to solve the DC biased problem significantly reduces the memory requirement compared to the conventional HBFEM. The decomposed harmonic balance equations are solved efficiently by the block Gauss‐Seidel algorithm combined with the relaxation iterative scheme. An investigation on DC bias phenomena is carried out through the harmonic solution of the magnetic field. The decomposed HBFEM can also be applied to solve 3‐D DC‐biasing magnetic field and eddy current nonlinear problems in a practical power transformer.

In many different engineering fields often there is a need to protect regions from electromagnetic interference. According to static and low-frequency magnetic fields the common strategy bases on using a shield made of conductive or ferromagnetic material. Another screening technique uses solenoids that generate an opposite magnetic field to the external one. The purpose of this paper is to discuss the shielding effect for a magnetic and conducting cylindrical screen rotating in an external static magnetic field.

The magnetic flux density is expressed in terms of the magnetic vector potential. Applying the separation of variables method analytical solutions are obtained for an infinitely long magnetic conducting cylindrical screen rotating in a uniform static transverse magnetic field.

Analytical formulas of the shielding factor for a cylindrical screen of arbitrary conductivity and magnetic permeability are given. A magnetic Reynolds number is found to be an appropriate indication of the change in magnetic field inside the screen. Useful simplified expressions are presented.

This paper treats in a qualitative way the possibility of static magnetic field shielding by using rotating conducting magnetic cylindrical screens. Analytical solutions are given. If the angular velocity is equal to zero or the relative permeability of the shield is equal to one the shielding factor has forms well known from literature.