1. Introduction
The analysis of electromagnetic devices remains an important research problem for the electrical engineering community. That class of devices exploits the induction effect for coils located on a magnetic circuit. Nowadays, field-circuit models are very often used to model precisely those devices, which combine electromagnetic field equations and circuit equations of coils. This means that partial differential equations of the magnetic circuit and ordinary differential equations of coils must be solved together. Two classes of magnetic fields are distinguished with respect to the properties of material constituted magnetic circuit: magneto-static, for nonlinear but nonconducting material, and dynamic, for nonlinear and conducting material.
The partial differential equations for magnetic vector potential describing electromagnetic fields are solved numerically on the discretized area of the magnetic circuit. Many approaches are known for that [
1,
2,
3]. Magneto-statics problems lead to a rather high number of algebraic nonlinear equations, magneto-dynamics leads to ordinary nonlinear differential equations. The field-circuit models of electromagnetic devices with magneto-statics fields can be reduced to so-called circuit models of magnetically coupled coils because the unique relations between the linked fluxes and coils currents can be determined. The number of such equations is usually much lower than the number of algebraic equations for magnetic circuit. They form a set of ordinary nonlinear differential equations, which should be solved together with algebraic nonlinear equations of the magneto-statics field. Solving complicates significantly for the field-circuit models accounting for nonlinear and conducted magnetic material. All equations, both magnetic field and coils, constitute one big set of ordinary nonlinear differential equations, which must be solved numerically.
Steady-state solutions are the most interesting ones for engineering applications, and direct steady-state determination for electromagnetic devices remains a research topic. To identify the steady states, the equations are typically solved numerically until the transient component disappears, and only the steady-state components remain [
4,
5,
6,
7]. However, such an approach may be inefficient owing to its high computational time and the creation of unnecessary data. Direct steady-state analysis methods were developed to overcome these challenges. The steady states in electromagnetic devices are usually periodic, which helps to predict the solution. The nonlinearities embedded in their mathematical models make the time-domain approach [
8,
9,
10,
11,
12,
13,
14,
15] preferable over the frequency-domain approach [
16,
17,
18,
19,
20,
21,
22]. Effective algorithms for solving such problems are still being developed.
This article presents an alternative approach for directly determining steady states for field-circuit models of electromagnetic devices, accounting for both nonlinearity and conductivity of magnetic material. The solution combines space and time and is predicted as periodic functions with respect to space variables and time. The finite difference is created using the Discrete Differential Operators (DDOs) with respect to each variable. Those operators relate the value of the first derivative in each discretization point to the values of the function in all discretization points over the period. The DDO with respect to time has been successfully applied to direct steady-state calculations of circuit models of electromagnetic devices in [
14,
15]. The DDO for 1D and 2D problems have been developed and applied in [
23,
24,
25,
26] for analysis of magneto-static field in a transformer window.
These experiences allowed us to extend this methodology to magnetodynamic problems, and field-circuit models were selected as being of great interest to engineering. Such a methodology did not exist (and cannot be found) in the field-circuit model literature. Due to the test nature of the research, a very simple electromagnetic object was selected for analysis. It is an electromagnet coil with a ferromagnetic and conductive cylindrical core when supplied with voltage, which can be described by the field-circuit model using two variables: radius and time. The coil current, magnetic field distribution, flux function and eddy current density in a non-linear and conductive solid magnetic material are calculated with respect to these two variables.
This paper presents in detail the solutions of the steady state equations using the DDOs with respect to space and time. These solutions were found using the minimizing procedure offered by MATLAB software. The resulting coil current and values characterizing the magnetic field in the core are presented in figures, with respect to time over one period. Numerical tests show the sensitivity of these values with respect to the amplitude and pulsation of sinusoidal mono-harmonic voltage supplying the coil.
2. Time-Varying Field-Circuit Analysis by an Explanatory Example
In this article, a simple explanatory example of an axisymmetric system is considered, for which the potential depends only on radius and time. It will be an infinitely long iron rod with a wound winding.
A schematic cross-section of such a magnetic system is shown below,
Figure 1.
In the example shown, potential A depends only on r, so its derivative with respect to z is equal to zero.
The equations for a component of the magnetic potential A in axisymmetric coordinates depend only on radius r, z-coordinate, and time t.
Commonly used in numerical methods, the modification of the equations for
has the form
where:
dla dla
dla , dla
—linear or non-linear magnetization characteristic of isotropic materials
; components of flux density
—average flux of winding’s cross-section area
W—number of conductors in winding’s cross-section area
—number of conductors per winding’s cross-section area
—height of the model in the z-axis.
The area of the magnetic field described by these equations, both space and time, has been discretized. The distribution of discretization points is shown in
Figure 2.
The solution of Equation (1) is predicted in the form of a periodic function of two variables
r and
t approximated by a double Fourier series
This form of the Fourier series ensures the fulfillment of the Neumann condition on the edges of the region, which, given the condition
makes the solution unambiguous. In the discrete version, the solution takes the form
where:
To solve the field-circuit problem in a discrete set of points, discrete differential operators differentiating with respect to r and t for functions and with respect to t for functions and were used. They relate the values of the respective partial derivatives to the function values at the grid points. How to create them is described below.
The first partial derivative with respect to
r has the form:
or in discrete domain
The first partial derivative with respect to
t has the form:
or in discrete domain
where:
,
,
.
The
potential values at the grid points are represented as a vector
The values of other functions at the grid points are arranged similarly.
In the presented approach, derivatives with respect to
r and
t are treated separately and discrete derivative operators are constructed separately for them. If we write the
vector as a matrix
Then the derivative with respect to time will have the operator and will differentiate with respect to the columns i.e., , while the derivative with respect to r will have the operator and will differentiate with respect to the rows i.e., .
2.1. Discrete First-Derivative Operators with Respect to ‘r’
The inverse discrete cosine and sine transforms are used to construct the discrete first-derivative operators with respect to r because, after the first differentiation, the base of the series changes from cosine to sine.
The first-derivative operators with respect to
r are also two, one of the cosine-sine base and second of the sine-cosine base and have the forms, respectively:
where
—frequency matrix.
The
operator can be used to calculate the magnetic flux density given the vector potential
at time
where
,
.
For materials with linear magnetization characteristics, the second derivative of
at time
is calculated as
For non-linear materials, the term
should be replaced with an expression
where
and the vector
has a form analogous to
, and
is a non-linear magnetization curve.
The matrices of cosine and sine transformations are almost unitary except for the element with index (1,1).
However, this does not distort the operation of the second derivative operator with respect to r, which is comprised of two first-derivative operators and , because element (1,1) concerns a constant signal, which is zeroed anyway as a result of multiplication by the frequency matrix . The operators and are, however, singular matrices.
2.2. Discrete Time First-Derivative Operator
A similar discretization as for r was performed for the time coordinate. The first derivative with respect to time can be computed using the discrete derivative operator with respect to time .
The inverse discrete Fourier transform was used to construct the discrete derivative operator
The operator of the first derivative with respect to time
has the form
where
—frequency matrix.
The first derivative of
with respect to time at a distance
is equal to
, while the first derivative of the coil current
and linked flux
are given by
and
, where
,
are vectors:
2.3. Discrete Field-Circuit Equations
2.3.1. Discrete Field-Circuit Equations for a Linear System
In the two-dimensional (time-space) case reduced to a single vector of searched values A, the Kronecker tensor product was used.
If
A1 is an n-by-m matrix and
A2 is a matrix, then the Kronecker tensor product of
A1 and
A2 is a large matrix formed by multiplying by
A2 each element of
A1The vector of discrete time-space values of the vector potential has the form
where the individual parts of A correspond to iron and copper areas and have the form
Vectors of time samples of current and voltage signals of the supply circuit have the following form:
Using the above, we can write the field-circuit equation for linear materials as follows:
where
2.3.2. Discrete Field-Circuit Equation for a Non-Linear System
For nonlinear magnetic materials, the field-circuit equation has the form
In which the vectors and have forms analogous to and . Their element values are computed from the characteristic of the ferromagnetic and copper material for each element of vector and , respectively, and g is the weight for the circuit equation.
The new operator matrices are calculated as follows:
2.4. Methodology of Solving Discrete Equations
The non-linear Equation (14) is solved by minimizing the function
Using MATLAB’s FMINUNC.
In addition, in order to converge the minimization process to the correct minimum value, a gradient and Hessian must also be provided
where:
3. Calculation Results
The aim of the research is to demonstrate the usefulness of the method of discrete differential operators (DDOs) for the qualitative determination of physical phenomena, but also the accuracy of calculating the field values for non-linear and conductive magnetic materials. For testing, the coil is powered by a sinusoidal monoharmonic voltage that changes its amplitude and pulsation. The results obtained with the DDOs method were compared with the results obtained with the FEM method using the Ansys-Maxwell program.
In the discussed example, the system is axisymmetric and infinite in the z-axis. The space-time model applies to a segment with dimensions L = 12 mm L
Cu = 10 mm
and the length of time equal to one period of the supply voltage. The fragment was discretized in the
r axis (
points) and the t axis (
points). The equivalent external circuit (supplying the analyzed part of the coil) consists of a sinusoidal voltage
u with an amplitude of 1.5 V and a small lumped resistance
. The number of turns in this fragment was equal to W = 10. The following material constants for iron were assumed: bulk conductivity
S/m, B-H curve
Figure 3.
Figure 4 shows the coil current obtained by DDO and FEM for a sinusoidal voltage of 50 Hz, the amplitude of which is large enough to saturate the magnetic core, producing a characteristic non-sinusoidal current waveform mainly due to third harmonic.
Figure 5 shows the difference between the currents calculated from the two aforementioned methods resulting mainly from the different grid adopted by Ansys (much dense FEM grid) and the fact that the Ansys-Maxwell program calculated the near-steady state, i.e., the fourth period of the transient state simulation.
The flux function calculated by the operator is shown in
Figure 6. Its changes in time and space could be drawn entirely on the “r-t” plane.
It should be noted that the cross-sections against the radius are nearly sinusoidal in time because the coil is supplied with monoharmonic voltage and the flux function should follow the shape of the voltage.
For small radiuses these functions are damped by eddy currents, which is well known phenomenon.
Also the difference between the calculation result obtained by the DDOs method and FEM is shown in
Figure 7.
The calculation error between the operator and the FEM method is of the order of which means that both methods are comparable.
The magnetic induction waveforms at selected points for several radiuses (given in
Table 1) are presented below.
In
Figure 8 the waveforms are strongly nonlinear due to core saturation and are shifted in time due to eddy currents. Their magnitudes decrease lite bit inside core.
The next figure shows the changing of current density as a function of time and radius and the waveform of the current density at the points given in
Table 1 is shown below (
Figure 9).
In
Figure 10, the current density waveforms were obtained using the DDOs method. It can be observed that the waveforms in the core are in phase opposite to the coil current and are attenuated inside the core, which is a fundamental phenomenon for an alternating magnetic field in a conductive material. The non-linearity of their deformation depends on the radius.
For comparison,
Figure 11 shows the waveform of the current density obtained with the Ansys-Maxwell program in the fourth period of the transient state simulation at the same points
r (
Table 1). The waveforms keep the same character and both are qualitatively and quantitatively identical.
For test purposes, the magnitude of the supply voltage is significantly reduced to keep in linear part of B-H curve.
Figure 12,
Figure 13 and
Figure 14 show the resulting waveforms of coil current, current density and flux density, respectively. All curves are almost sinusoidal but tendencies of current density and flux density for changes with the radius are kept.
The lack of non-linear distortions results from lower magnetic induction values in the range of 0–0.8 T, i.e., in the range of linear changes.
Similar calculations were performed for the voltage excitation with 100 Hz pulsation.
The current and voltage waveforms of the field-excitation coil are shown below (
Figure 15).
The maximum current here is more than 2 times lower than at 50 Hz, therefore the core saturation is lower and the content of higher harmonics is lower.
The flux-function calculated by the DDOs is presented above the ‘r-t’ plane in
Figure 16.
The following figures show the results of the DDOs calculations analogous to those for 50 Hz pulsations.
Figure 17 shows the changing current density as a function of time and radius, it shows a clearly smaller penetration depth compared to the case for 50 Hz (see
Figure 8).
Due to the smaller range of the field impact (stronger skin-effect), other points (given in
Table 2) were selected to present the waveforms
The current density vs. time waveforms for several radius values given in
Table 2 are shown below (
Figure 18).
A smaller field penetration (skin depth) is noticeable here, i.e., a faster decrease in the amplitudes of the induced currents with the distance from the supply winding.
The magnetic flux time waveforms at the same points
r are shown below (
Figure 19).
4. Discussion
Calculations testing the proposed approach showed that the results very well reflect the qualitative and quantitative nature of both the current in the coil and the distribution of magnetic induction as well as eddy currents in the magnetic core. The quantitative assessment was made by comparing the results of the transient simulations with the commercial ANSYS package, confirming the very good quantitative agreement of the results, despite the much sparser rectangular space-time discretization grid for the discrete difference operators used. Calculation time using discrete differential operators was slightly shorter, despite the use of commonly available MATLAB procedures. The obtained calculations of the magnetic field and the supply current using the DDO method show all the characteristic relationships for the time-varying magnetic field in conductive and non-linear materials. The field suppression (skin effect) resulting from the theory and the decrease of the induced currents can be seen with changing the distance from the magnetizing winding. A stronger effect is obtained by increasing the frequency of the supply voltage. The currents induced in the iron are also reversed in phase with respect to the supply current. The field attenuation strongly depends also on the conductivity of the medium. The effect of entering into non-linearity and saturating the core is also visible and consistent with the theory and was quantitatively compared with the results from the almost steady state (obtained after several periods of transient state) in the Ansys program in which the same excitation and circuit were used for power supply.
5. Conclusions
DDO operators have already been used by us to solve a 2D stationary field, e.g., in a transformer window, where the cosine basis worked very well. It ensured the reversibility of the cosine transform operators, allowed the modeling of discontinuous fields at the edges of the region, and easy fulfillment of the Dirichlet and Neumann boundary conditions.
This time we presented an alternative method for the analysis of steady states in electromagnetic systems described by field-circuit models for nonlinear and conductive magnetic media. Discrete space-time differential operators, created in the paper, approximating the partial derivatives appearing in the field equations, and a discrete-time differential operator for magnetizing coils, were used for the analysis. This allowed direct calculations of steady states both in the magnetic core and in the coils, eliminating the calculation of transient states.
The proposed approach was tested on a simple case of a solenoid coil with ferromagnetic and conductive cylindrical core using two variables: the radius and time.