Magnetic Field of a Hollow Cylinder

© 1998 Waterloo Maple Inc.

NOTE: This worksheet shows how Maple can be used to determine the magnetic screening provided by a hollow cylinder.

Introduction

restart;

A ferromagnetic tube can be used to shield its interior from static magnetic fields. To calculate the screening constant, the tube is presumed to be set up so that its axis is perpendicular to the direction of an outer homogeneous magnetic field H0 [A/m]. The outer radius of the tube is Ro [m], and the inner radius is Ri [m].

[Maple OLE 2.0 Object]

Theory

Under the influence of the outside magnetic field, the ferromagnetic material is magnetized. The resulting magnetic field can be calculated as a superposition of the outside magnetic field and the fields of magnetic line dipoles. The free constants in the equations are calculated with the help of boundary conditions.
The magnetic scalar potential V [A] is defined by H =-grad( V ).
The scalar magnetic potential of a homogeneous magnetic field in polar coordinates (r [m], [Maple Math] [rad]) is given by
[Maple Math] (characteristic constants A [A], B [A/m]).

The magnetic field of a line dipole is given by
[Maple Math] (characteristic constant C [A*m]).

Solution

There are three equations that represent the superposition of a homogeneous magnetic scalar potential and the scalar potential of a line dipole: one for the outside, a second for the wall and a third for the inside of the tube.

The permeabitiy of the material is defined as

[Maple Math]

( [Maple Math] : permeability of vacuum [H/m], mur: relative permeability [-])

Outside:

[Maple Math]

Wall:

[Maple Math]

Inside:

[Maple Math]

On the inside, the potential must be finite so that C3 = 0. On the axis where r = 0 and where [Maple Math] , the potential has the value V3 = V2 = V1 = V0 = A1 = A2 = A3.

Choose V0 = 0. The result on the inside is a homogeneous magnetic field with a field intensity of Hi, so that B3 = Hi.

A3:=0: A2:=0: A1:=0: C3:=0: B3:=Hi:

On the outside, the magnetic field in very big distances is undisturbed by the magnetic field caused by the tube, so that

[Maple Math] ~ [Maple Math] .

This results in B1 = H0.

B1:=H0:

The conditions at the interfaces are that the tangential components of the vector H and the normal components of the vector B (B = [Maple Math] *H) are continuous.

From [Maple Math] in cylinder coordinates, one obtains the following equations:

with(linalg):
Warning, new definition for norm
Warning, new definition for trace

[Maple Math]

h1:=grad(V1, v, coords=cylindrical);

[Maple Math]

h2:=grad(V2, v, coords=cylindrical);

[Maple Math]

h3:=grad(V3,v,coords=cylindrical);

[Maple Math]

At r = Ri, the tangential components of H must be continuous (component of h1, h2 and h3 with 'd/d( [Maple Math] )').

eqHi:=simplify(evalf(subs(r=Ri,h2[2]))=evalf(subs(r=Ri,h3[2])));

[Maple Math]

At r = Ro, the tangential components of H must be continuous.

eqHa:=simplify(evalf(subs(r=Ro,h1[2]))=evalf(subs(r=Ro,h2[2])));

[Maple Math]

At r = Ri, the normal components of B must be continuous (component of h1, h2 and h3 with 'd/d(r)').

eqBi:=simplify(evalf(subs(r=Ri,mu*h2[1]))=evalf(subs(r=Ri,mu0*h3[1])));

[Maple Math]

At r = Ro, the normal components of B must be continuous.

eqBa:=simplify(evalf(subs(r=Ro,mu0*h1[1]))=evalf(subs(r=Ro,mu*h2[1])));

[Maple Math]

Solving these equations for C1, C2, B2 and Hi gives

sol:=simplify(solve({eqHi,eqHa,eqBi,eqBa},{C1,C2,B2,Hi}));

[Maple Math]
[Maple Math]
[Maple Math]
[Maple Math]

assign(sol);

The result for the screening factor 'Screen' (= Hi/H0) is then given by

Screen:=simplify(Hi/H0);

[Maple Math]

Graphical Display of the Solution

For the plots, values are assigned to the parameters. The relative permeability for ferromagnetic materials is a very small value of mur = 6.

Ri:=1: Ro:=2: H0:=10: mur:=6:

For the contour plot, the plotted function must be written in cartesian coordinates so the radius r and the angle phi are substituted by x and y.

[Maple Math]

[Maple Math]

[Maple Math]

The plot a piecewise function is defined.

[Maple Math]

Now the equipotential lines can be plotted. In addition, two circles are included so that the location of the tube is indicated.

with(plottools):
c1 := circle([0,0], Ri, color=red): c2 := circle([0,0], Ro, color=red): p:=plots[contourplot](Pot,x=-3..3,y=-3..3,scaling=CONSTRAINED, color=blue,numpoints=4000):
plots[display](c1,c2,p);

[Maple Plot]

In this case, the screening factor is about 0.56, so that the magnetic field inside the tube is still more than half of the value on the outside.

Screen;

[Maple Math]

Summary

For some applications, magnetic fields need to be shielded. In the preceeding calculation, an illustrative example is given on how this can be accomplished. A ferromagnetic tube (with a very low relative permeability) is placed in a static magnetic field and the screening factor is calculated.