# Mutual Inductance and Magnetic Force Calculations for Bitter Disk Coil (Pancake) with Nonlinear Radial Current and Filamentary Circular Coil with Azimuthal Current

Authors: Slobodan Babic; Cevdet Akyel
Journal: Advances in Electrical Engineering
 (2016)

---

## Abstract

Bitter coils are electromagnets used for the generation of extremely strong magnetic fields superior to 30 T. In this paper we calculate the mutual inductance and the magnetic force between Bitter disk (pancake) coil with the nonlinear radial current and the circular filamentary coil with the azimuthal current. The close form expressed over complete elliptic integrals of the first and second kind as well as Heuman’s Lambda function is obtained for this configuration either for the mutual inductance or for the magnetic force. The results of this method are compared with those obtained by the improved modified filament method for the presented configuration. All results are in an excellent agreement.


---

## Body

## 1. Introduction

In the literature and scientific papers the calculation of the mutual inductance and the magnetic force between ordinary circular coils with the uniform current densities have been given by different methods either analytical or numerical [1–15]. However, there are not a lot references about the mutual inductance and the magnetic force calculation between circular coils in which the current densities are not uniform. The coils with this characteristic are Bitter coils or Bitter solenoids. Bitter coils are used in high magnetic field applications and they differ from the ordinary coils in that they have an inverse radial distribution of current [16, 17]. The interesting point is the calculation of the mutual inductance and the magnetic force between Bitter coils or between one Bitter coil and an ordinary coil. Generally, the calculation of the magnetic force between circular coils is closely related to the calculation of their mutual inductance. Since their mutual energy is equal to the product of their mutual inductance and the currents in the coils, the component of the magnetic force of attraction or repulsion in any direction is equal to the product of the currents multiplied by the differential coefficient of the mutual inductance taken with respect to that coordinate. The mutual inductance and the magnetic force are obtained over multiple integrals with the different kernel functions. Thus, in the calculation of the mutual inductance and the corresponding magnetic force between two coils we need to integrate their kernel functions depending on coil configurations. Hopefully, these kernel functions are Green functions r
-
0.5 and r
-
1.5, where the “r” is the distance between two coils. This integral approach, much easier than the differential approach, leads to relatively simple expressions which are incorporated in these two physical quantities. In this paper we calculate the mutual inductance between the Bitter disk coil with nonlinear radial current and a circular filamentary coil with the azimuthal current (ordinary coil). Coils are coaxial and in air. Either the mutual inductance or the magnetic force is obtained in closed form expressed over complete elliptic integrals of the first and second kind as well as over Heuman’s Lambda function [18, 19]. The results of these calculations will be compared to those obtained by the improved filament method for concerned configuration. The results obtained by these two methods are in an excellent agreement.


## 2. Basic Expressions

The mutual inductance and magnetic force between a disk coil (pancake) with the uniform azimuthal current densityJ
1
=
N
1
I
1
/
(
R
2
-
R
1
) and a filamentary coil with the azimuthal current I
2 (see Figure 1) can be calculated, respectively, by [1, 2](1)
M
=
μ
0
N
1
R
R
2
-
R
1
∫
0
π
∫
R
1
R
2
r
I
cos
⁡
θ
d
r
I
d
θ
r
0
,
(2)
F
=
-
μ
0
N
1
I
1
I
2
R
z
2
-
z
1
R
2
-
R
1
∫
0
π
∫
R
1
R
2
r
I
cos
⁡
θ
d
r
I
d
θ
r
0
3
,where(3)
r
0
=
z
2
-
z
1
2
+
r
I
2
+
R
2
-
2
r
I
R
cos
⁡
θand N
1 is the number of turns of the disk coil (pancake).Figure 1
The thin Bitter disk coil (pancake) and the filamentary circular coil.μ
0
=
4
π
×
10
-
7 H/m is the permeability of free space (vacuum), and r, θ, z are the cylindrical coordinates. Both coils are in air or in a nonmagnetic and nonconducting environment.In [2] the magnetic force is obtained in the close form expressed over the complete elliptic integral of the first kind and Heuman’s Lambda function and the magnetic force in the semianalytical form expressed over complete elliptic integral of the first and second kind, Heuman’s Lambda function plus one simple integral which has to be evaluated by some of the numerical integrations [1].Let us suppose the disk coil be the Bitter coil in which the current density is not uniform [16, 17] and given by(4)
J
1
=
N
1
I
1
r
I
ln
⁡
R
2
/
R
1
.From (1), (2), and (4) we obtain the expressions of the mutual inductance and the magnetic force between the Bitter coil and the filamentary coil as follows:(5)
M
B
=
μ
0
N
1
R
ln
⁡
R
2
/
R
1
∫
0
π
∫
R
1
R
2
cos
⁡
θ
d
r
I
d
θ
r
0
,
(6)
F
B
=
-
μ
0
N
1
I
1
I
2
R
z
2
-
z
1
ln
⁡
R
2
/
R
1
∫
0
π
∫
R
1
R
2
cos
⁡
θ
d
r
I
d
θ
r
0
3
,where(7)
r
0
=
z
2
-
z
1
2
+
r
I
2
+
R
2
-
2
r
I
R
cos
⁡
θ
.


## 3. Calculation Method

Integrating in (5) and (6) over r
I and θ (substituting θ
=
π
-
β) the mutual inductance and the magnetic force between the Bitter disk and the filamentary circular coil can be expressed, respectively, in an analytical form as follow (see Appendices A and B):(8)
M
B
=
μ
0
N
1
2
ln
⁡
R
2
/
R
1
∑
n
=
1
n
=
2
-
1
n
-
1
T
n
,
(9)
F
B
=
μ
0
N
1
I
1
I
2
ln
⁡
R
2
/
R
1
∑
n
=
1
n
=
2
-
1
n
-
1
S
n
,where(10)
T
n
=
k
n
R
R
n
2
R
R
2
+
z
Q
2
-
R
n
2
-
R
2
-
z
Q
2
K
k
n
+
k
n
R
R
n
R
+
R
n
2
+
z
Q
2
E
k
n
-
π
z
Q
V
n
,
S
n
=
k
n
z
Q
R
R
n
R
R
2
+
z
Q
2
+
R
K
k
n
-
π
2
sgn
⁡
z
Q
V
n
,
V
n
=
1
-
Λ
0
θ
1
n
,
k
n
+
sgn
⁡
R
2
+
z
Q
2
-
R
n
1
-
Λ
0
θ
2
n
,
k
n
,
z
Q
=
z
2
-
z
1
,
k
n
2
=
4
R
R
n
R
+
R
n
2
+
z
Q
2
≤
1
,
m
=
2
R
R
+
z
Q
2
+
R
≤
1
,
θ
1
n
=
arcsin
⁡
z
Q
R
2
+
z
Q
2
+
R
,
θ
2
n
=
arcsin
⁡
1
-
m
1
-
k
n
2
,
k
n
2
≤
m
.Thus, all expressions are obtained in the closed form expressed over complete elliptic integrals of the first and second kind as well as Heuman’s Lambda function [18, 19].


### 3.1. Special Cases and Singularities

Mutual Inductance. In the mutual inductance calculation special cases and singular cases appear for z
Q
=
0 and k
n
2
≠
1 or z
Q = 0 and k
n
2
=
1.


#### 3.1.1.z
Q
=
0 and k
n
2
≠
1

This case is the regular case covered by (8). In this case the term T
n can be given in the following simplified form:(11)
T
n
=
k
n
0
R
+
R
n
R
R
n
R
-
R
n
K
k
n
0
+
R
+
R
n
E
k
n
0
,
(12)
k
n
0
2
=
4
R
R
n
R
+
R
n
2
≤
1
.


#### 3.1.2.z
Q
=
0 and k
n
2
=
1

This case is the singular case (R
=
R
1 or R
=
R
2) for which the term T
n is equal either 4
R
1 or 4
R
2.Magnetic Force. In the magnetic force calculation for z
Q
=
0 its value is zero because the coils are in the same plane.Notation Above. Λ
0 corresponds to Heuman’s Lambda function as defined in [18, 19] and “sign” function returns integer indicating the sign of a number. Functions K and E correspond to complete elliptic integrals of the first and second kind [18, 19].


## 3.1. Special Cases and Singularities

Mutual Inductance. In the mutual inductance calculation special cases and singular cases appear for z
Q
=
0 and k
n
2
≠
1 or z
Q = 0 and k
n
2
=
1.


### 3.1.1.z
Q
=
0 and k
n
2
≠
1

This case is the regular case covered by (8). In this case the term T
n can be given in the following simplified form:(11)
T
n
=
k
n
0
R
+
R
n
R
R
n
R
-
R
n
K
k
n
0
+
R
+
R
n
E
k
n
0
,
(12)
k
n
0
2
=
4
R
R
n
R
+
R
n
2
≤
1
.


### 3.1.2.z
Q
=
0 and k
n
2
=
1

This case is the singular case (R
=
R
1 or R
=
R
2) for which the term T
n is equal either 4
R
1 or 4
R
2.Magnetic Force. In the magnetic force calculation for z
Q
=
0 its value is zero because the coils are in the same plane.Notation Above. Λ
0 corresponds to Heuman’s Lambda function as defined in [18, 19] and “sign” function returns integer indicating the sign of a number. Functions K and E correspond to complete elliptic integrals of the first and second kind [18, 19].


## 3.1.1.z
Q
=
0 and k
n
2
≠
1

This case is the regular case covered by (8). In this case the term T
n can be given in the following simplified form:(11)
T
n
=
k
n
0
R
+
R
n
R
R
n
R
-
R
n
K
k
n
0
+
R
+
R
n
E
k
n
0
,
(12)
k
n
0
2
=
4
R
R
n
R
+
R
n
2
≤
1
.


## 3.1.2.z
Q
=
0 and k
n
2
=
1

This case is the singular case (R
=
R
1 or R
=
R
2) for which the term T
n is equal either 4
R
1 or 4
R
2.Magnetic Force. In the magnetic force calculation for z
Q
=
0 its value is zero because the coils are in the same plane.Notation Above. Λ
0 corresponds to Heuman’s Lambda function as defined in [18, 19] and “sign” function returns integer indicating the sign of a number. Functions K and E correspond to complete elliptic integrals of the first and second kind [18, 19].


## 4. Modified Filament Method

In [17] the mutual inductance is calculated between the Bitter coil and the superconducting coil with uniform current density by using the filament method. Here we give the modified formulas for the mutual inductance and the magnetic force for the treated configuration (see Figure 2) using the filament method.Figure 2
Configuration of mesh matrix: Filamentary coil-thin disk coil (pancake).Applying some modification in the mutual inductance calculation given in [17] we deduced the mutual inductance and the magnetic force between the Bitter coil and the filamentary coil as follows: (13)
M
B
F
=
N
1
R
2
-
R
1
2
n
+
1
ln
⁡
R
2
/
R
1
∑
l
=
-
n
l
=
n
M
l
r
I
l
,
F
B
F
=
N
1
I
1
I
2
R
2
-
R
1
2
n
+
1
ln
⁡
R
2
/
R
1
∑
l
=
-
n
l
=
n
F
l
r
I
l
,where(14)
M
l
=
2
μ
0
r
I
l
R
k
l
1
-
k
2
l
2
K
k
l
-
E
k
l
,
F
l
=
-
μ
0
I
1
I
2
c
k
l
4
r
I
l
R
2
-
k
2
l
1
-
k
2
l
E
k
l
-
2
K
k
l
,
r
I
l
=
R
2
+
R
1
2
+
R
2
-
R
1
2
n
+
1
l
l
=
-
n
,
…
,
0
,
…
,
n
,
c
=
z
2
-
z
1
=
z
Q
,
k
2
l
=
4
r
I
l
R
r
I
l
+
R
2
+
c
2
.Expressions (13) will be used to confirm the validity of analytical formulas (8) and (9).


## 5. Examples

To verify the validity of the new formulas, we apply them to the following set of examples.Example 1.
Calculate the mutual inductance and the magnetic force between the Bitter disk coil and a filamentary circular coil with the following dimensions and the number of turns. Currents in coils are unit.
Filamentary Circular Coil. R
=
20 mm, z
1
=
0 mm.
Bitter Disk Coil. R
1
=
40 mm, R
2
=
60 mm, z
2
=
50 mm, N
1
=
100.
Applying (8) and (9) the mutual inductance and the magnetic force are, respectively,(15)
M
B
=
521.472202513
09271
n
H
,
F
B
=
-
16.016662833788
816
μ
N
.The computational time is about 0.037677 seconds.
By using the modified filament, (13) and (14), we obtain(16)
M
B
F
=
521.472202513
1178
n
H
,
F
B
F
=
-
16.016662833788
95
μ
N
.The number of subdivision was n
=
10000
00 and the computational time about 0.362069 seconds.Example 2.
Calculate the mutual inductance between a disk and a filamentary circular coil with the following dimensions and the number of turns.
Filamentary Circular Coil. R
=
20 mm, z
1
=
0 mm.
Disk Coil. R
1
=
40 mm, R
2
=
60 mm, z
2
=
0 mm, N
1
=
100.
Applying (8) and (9) the mutual inductance and the magnetic force are, respectively,(17)
M
B
=
1.7442876932070
914
μ
H
,
F
B
=
0
N
.The computational time is about 0.020657 seconds.
By using the modified filament, (13) and (14), we obtain(18)
M
B
F
=
1.7442876932070
58
μ
H
,
F
B
F
=
0
N
.The number of subdivision was n
=
20000 and the computational time about 0.030564 seconds.Example 3.
Calculate the mutual inductance between a disk and a filamentary circular coil with the following dimensions and the number of turns.
Filamentary Circular Coil. R
=
40 mm, z
1
=
0 mm.
Disk Coil. R
1
=
40 mm, R
2
=
60 mm, z
2
=
0 mm, N
1
=
100.
This case is the singular case becauseR
=
R
1.
Applying (8), (9), and (11) the mutual inductance and the magnetic force are, respectively,(19)
M
B
=
10.93167
7144024915
μ
H
,
F
B
=
0
N
.The computational time is about 0.020223 seconds.
By using the modified filament, (13) and (14), we obtain(20)
M
B
F
=
10.93167
597572536
μ
H
,
F
B
F
=
0
N
.The number of subdivision was n
=
10000
00 and the computational time about 0.513017 seconds.Example 4.
Calculate the mutual inductance between a disk and a filamentary circular coil with the following dimensions and the number of turns.
Filamentary Circular Coil. R
=
60 mm, z
1
=
0 mm.
Disk Coil. R
1
=
40 mm, R
2
=
60 mm, z
2
=
0 mm, N
1
=
100.
This case is also the singular case for whichR
=
R
2.
Applying (8), (9), and (11) the mutual inductance and the magnetic force are, respectively,(21)
M
B
=
14.06205
6002507577
μ
H
,
F
B
=
0
N
.The computational time is about 0.03929 seconds.
By using the modified filament, (13) and (14), we obtain(22)
M
B
F
=
14.06205
49882661
μ
H
,
F
B
F
=
0
N
.The number of subdivision was n
=
10000
00 and the computational time about 0.487595 seconds.In comparative calculations the significant figures which agree are bolded.Thus, with these examples we confirmed the validity of the present analytical method for calculating the mutual inductance and the magnetic force between the Bitter disk with the nonlinear radial current distribution and an ordinary filamentary circular coil. All calculations were executed in Matlab programming.


## 6. Conclusion

The new accurate mutual inductance and magnetic force formulas for the system of the Bitter disk with the nonlinear radial current and the filamentary circular with the azimuthal in air are derived and presented in this paper. All expressions for either the mutual inductance or the magnetic force are obtained in the close form expressed over complete elliptic integrals of the first and second kind and Heuman’s Lambda function. Also we gave in this paper the improved formulas for the mutual inductance and the magnetic force between treated coils by using the filament method. All singular cases are obtained in a closed form. The presented method can be used to calculate these important electrical quantities such as the mutual inductance and the magnetic force for coils combinations comprising Bitter coils of rectangular cross section and the circular coils with uniform current densities.



---

*Source: 3654021-2016-09-18.xml*