Jump to content

User:Thypingvinen/sandbox/Halbach

From Wikipedia, the free encyclopedia

The method used to find the field created by the cylinder is mathematically very similar to that used to investigate a uniformly magnetised sphere.[1]

Because of the symmetry of the arrangement along the cylinder's axis, the problem can be treated as two-dimensional. Work in plane-polar coordinates with associated unit vectors and , and let the the cylinder have radial extent . Then the magnetisation in the cylinder walls, which has magnitude , rotates smoothly as

,

while the magnetisation vanishes outside the walls, that is for the bore and surroundings .

By definition, the auxilliary magnetic field strength is related to the magnetisation and magnetic flux density by . Using Gauss' law , this is equivalently

. (1)

Since the problem is static there are no free currents and all time derivatives vanish, so Ampère's law additionally requires , where is the magnetic scalar potential (up to a sign under some definitions). Substituting this back into the previous Equation 1 governing and , we find that we need to solve

, (2)

which has the form of Poisson's equation.

Consider now the boundary conditions at the cylinder-air interfaces and . Integrating over a small loop straddling the boundary and applying Stokes' theorem requires that the parallel component of is continuous. This in turn requires that is continuous across the boundary. (More properly this implies that must differ by a constant across the boundary, but since the physical quantities we are interested in depend on gradients of this potential, we can arbitrarily set the constant to zero for convenience.) To obtain a second set of conditions, integrate Equation 1 across a small volume straddling the boundary and apply the divergence theorem to find

,

where the notation denotes a jump in the quantity across the boundary, and in our case the sign is negative at and positive at . The sign difference is due to the relative orientation of the magnetisation and the surface normal to the part of the integration volume inside the cylinder walls being opposite at the inner and outer boundaries.

In plane-polar coordinates, the divergence of a vector field is given by

. (3)

Similarly, the gradient of a scalar field is given by

. (4)

Combining these two relations, the Laplacian becomes

. (5)

Using Equation 3, the magnetisation divergence in the cylinder walls is

.

Hence Equation 2, which is that we want to solve, becomes by using Equation 5

. (6)

Look for a particular solution of this equation in the cylinder walls. With the benefit of hindsight, consider , because then we have

and also

.

Hence , and comparison with Equation 6 shows that is the appropriate particular solution.

Now consider the homogenous equation for Equation 6, namely . This has the form of Laplace's equation. Through the method of separation of variables, it can be shown that the general homogenous solution whose gradient is periodic in (such that all the physical quantities are single-valued) is given by

,

where the are arbitrary constants. The desired solution will be the sum of the particular and homogenous solutions that satisfies the boundary conditions. Again with the benefit of hindsight, let us set most of the constants to zero immediately and assert that the solution is

,

where now are constants to be determined. If we can choose the constants such that the boundary conditions are satisfied, then by the uniqueness theorem for Poisson's equation, we must have found the solution.

The continuity conditions give

(7)

at the inner boundary and

(8)

at the outer boundary. The potential gradient has non-vanishing radial component in the cylinder walls and in the bore, and so the conditions on the potential derivative become

at the inner boundary and

at the outer boundary. Note that these are identical to Equations 7 and 8, so indeed the guess was consistent. Hence we have and , giving the solution

.

Consequently the magnetic field is given by

,

while the magnetic flux density can then be found everywhere using the previous definition . In the bore, where the magnetisation vanishes, this reduces to . Hence the magnitude of the flux density there is

,

which is independent of position. Similarly, outside the cylinder the magnetisation also vanishes, and since the magnetic field vanishes there, the flux density does too. So indeed the field is uniform inside and zero outside the ideal Halbach cylinder, with a magnitude depending on its physical dimensions.

  1. ^ Fitzpatrick, Richard (2014-06-27). "Uniformly Magnetized Sphere". The University of Texas at Austin. Retrieved 2022-03-28.