GB2193322A

GB2193322A - Magnetic field coils

Info

Publication number: GB2193322A
Application number: GB08714434A
Authority: GB
Inventors: Robert Turner; Peter Mansfield; Barry Leonard Walter Chapman
Original assignee: National Research Development Corp UK
Current assignee: National Research Development Corp UK
Priority date: 1986-06-28
Filing date: 1987-06-19
Publication date: 1988-02-03
Anticipated expiration: 2007-06-19
Also published as: EP0252634A2; CA1276979C; GB8615854D0; GB8714434D0; EP0252634B1; US4896129A; JPH0793210B2; GB2193322B; DE3751921T2; EP0252634A3; DE3751921D1; JPS6387712A

Description

GB2193322A 1

SPECIFICATION

Magnetic field coils

The present invention relates to magnetic field coils and more particularly to coils designed to 5 produce a defined field in at least one volume.

In many applications it is important to provide a magnetic field the characteristic of which, for example, in an x-y plane is known. For example, in virtually every application of NIVIR spectros copy and imaging it is essential to use magnetic fields whose spatial variation is controlled. Such fields may be uniform, or they may have a linear variation in either x, y or z (where x, y, z are 10 the co-ordinate axes), or they may be uniform within a certain volume and rapidly changing elsewhere, as in TIVIR. For fine-tuning correction of nearly uniform main fields it is useful to have a set of shim coils, which provide known correction gradients, varying as xy, X2 - y2, yz etc.

It is an object of the present invention to provide coils which are designed in a manner suitable for generating any of the fields mentioned. 15

The coils are formed by conductors, which lie on a surface, and the positions of the conduc tors are derived directly from the specification of the field desired at a further target surface which may be smaller or larger than the former surface. Coils designed according to the principle to be described can provide very large volumes within which the field follows a specified variation, and are thus exceptionally efficient and compact. 20 According to the present invention, there is provided a magnetic coil including a plurality of conductors positioned on a surface or surfaces in which the positions of the conductors are determined by the calculation of the continuous current density required on that surface to generate a specified magnetic field everywhere on a further target surface or at specified points thereon, and calculation of conductor positions on the former surface which provide a good 25 approximation to this current density.

In one class of preferred surfaces, the surface or surfaces are cylinders and the target surface or surfaces are also cylinders.

The cylinder on which the conductors are positioned may, for example, be a cylinder of a non ferromagnetic material such as glass reinforced plastic, or the cylinder may be notional, the coils 30 being merely shaped to conform to a cylindrical or partially cylindrical shape.

In a first preferred embodiment, the cylinder on which the coil conductors are to be mounted is of larger diameter than the diameter of the target cylinder, the coil therefore being designed to produce a known field within the cylinder. The target cylinder may be one fifth to one tenth the diameter of the cylinder or could be zero diameter. 35 In a second preferred embodiment the target cylinder is of larger diameter than the cylinder.

This is particularly useful where a specified field is required outside a coil. A specific case is where a zero field is required outside a coil as is advantageous in NIVIR spectroscopy and imaging.

In a further embodiment the coil design may be obtained from calculations carried out on first 40 and second target cylinders, the first being of smaller diameter and the second of larger diameter. The coil designed by this method can therefore produce defined magnetic fields both within the first cylinder and outside the second cylinder.

In a further embodiment a plurality of coaxial cylindrical or semicylindrical coils can be designed to produce various fields specified on target cylinders in all but one of the intervening 45 spaces, and within the centre cylinder.

In a preferred embodiment the field within the inner cylinder is uniform and outside the outermost cylinder the field is zero. Alternatively, the inner field may be zero with a defined outer field.

The present invention provides a method of producing a coil wound on a surface to produce a 50 desired magnetic field when energised including calculating the continuous current density re quired on the surface to generate a specified magnetic field on the surface of a target surface of different size, calculating conductor positions on the surface to provide a near approximation to this current density and winding the conductors onto the surface in the calculated positions.

Embodiments of the present invention will now be described, by way of example with 55 reference to the accompanying drawings in which:

Figure 1 shows a cylinder or notional cylinder used to explain the principle of the present invention; Figure 2 shows an end view of a squirrel cage r.f. coil with six turns; Figures 3a, b, c and d show the design of a saddle coil for production of an x gradient (one 60 octant only of the coil being shown) and uniformity maps for d = 1.5a and return path at d 2.0a; Figure 4 shows a design of saddle coil (one octant only) for x gradient to be driven in parallel, the return paths not being shown; Figure 5 shows the variation of target magnetic field (Bz) with z for an end corrected saddle 65

2 GB2193322A 2 coil; Figure 6 shows a design of end corrected saddle coil (one octant only); Figure 7 shows a design of longitudinal gradient coil with 12 arcs; and Figure 8 shows the current density required to produce a uniform magnetic field.

The calculation of current density required to generate given field according to the present 5 invention is given as follows:

The crucial observation leading to improved coil designs Is that the z component of the magnetic field produced by the current density F (0, z) flowing on a cylinder of radius a (Figure 1) is given by (see Patent Application 128688).

10 -LO i 00 C0 imo lkz m B (p,,0,,z) - E -CO f dke e k F (k) M(ka)I (kp) (1) z 21r m= -CX) m m v 15 where Fm(k) - 1 f " e -'mo'do 1 f ". e- ikz ' C1 F 0 2w -w -CD 20 Now suppose we are given B, (c, 0, z) for a particular value of c < a that is to say, we are given the field on the further cylinder of radius c. We may invert equation (1) by integrating each side with respect to 0 and z:

B W, k.') - 1- f Ir e-lm'o do 100 e -ikl Z dz B (CIPOPZ) 25 z 21r -ii CO z a k' F M1 (M) K',(kla) Imy(klc) (2) 30 0 m Thus we may write 35 m B z (m,k) Ii a k Kl(ka) I (ke) 0 M m 40 The current density F. (0, z) may then be obtained:

ikz m f dk e F (k) (4) F Tr in 45 1 c -'mo fo dk e -ikz B z (m,k) E e 2wo a m- -(D k Kl(ka) I (ke) 0 m m 50 Hence, given B (o, 0, z) we may obtain the required current F,, (0, z). Note that the other component of the current in the cylinder, F, (0, z) can be obtained using the equation of continuity, which gives 55 aF z 1 W 0 (5a) Wz - - a a 0 or in terms of m and k: 60 kFm(k) - EL Fm(k) - z a 0 (515) Equation (4) has a certain peculiarity. For large k the denominator kk'M(ka)lm(kc) tends to e-k(a-c). 65 3 GB2193322A 3 Hence the integral in equation (4) will fail to converge unless jB, (m, k)l decreases faster than e-l(a-c) as k--oo. In particular this rules out the possibility of creating a field which has the form B,, = Bjzj<cl = 0 elsewhere, since this gives a B. (m, k) which decreases only as 1/k.

A useful class of functions to describe the variation of B,(c, 0, z) is the set 5 g(z) + z (z)n (6a) d 10 where s and n (n even) are integers. For z < d/2 this diverges from zs only by order of (z/d)n, and the Fourier transform of g(z) is well known. For large k this decreases as 17 e7 s'L"x_Y1 p 15 requiring that the characteristic length d of the coil must be such that dsin 7r/n > a. (6b) For many imaging purposes it is sufficient to take n = 6; this gives correction up to sixth order in the z direction. For instance, if a transverse gradient is to be uniform in the z direction, we take s = 0. Then at z = d/2 the magnitude of the gradient has dropped by only 1/64, i.e. 1.5%. 20 Correction to higher order may be achieved by making n larger, but this increases the length of the coil.

The variation B, (c, 0, z) may also be specified point by point, and the Fourier transform B(m,k) computed numerically. It is still necessary to ensure that B(m,k) decrease with k suffici- ently fast for the integral in (4) to converge. 25 Note that most fields of interest have either cylindrical symmetry (m = 0 in equation (1)) or reflection symmetry about the xz or yz plane (m = 1 in equation (1)). Thus we are normally dealing with F a (k) z (0, k) (7) 30 a k K,'(ka) I a (ke) 35 g),r F' (k) = - - 1 B z (1, k) (8) 0 a k K11(ka) Il(ke) Once the continuous current distribution required to generate the specified field has been 40 found, it is necessary to find practical means for approximating such a distribution as closely as possible. In one preferred embodiment these currents are confined to wire arcs. There are two ways of arranging such arcs; (a) equal current in each arc, i.e. series activated, and (b) equal voltage drop down each arc, i.e. parallel activated. It is considerably simpler to calculate the arc positions required in case (a), though case (b) can also be dealt with in an entirely satisfactory 45 way, using an iterative method.

Arc Position Algorithm The basic principle for finding n arc positions for series activation, given the surface current density, is to subdivide the surface of the cylinder into 2n areas, in each of which flows an 50 equal integrated surface current. These areas are of course separated by contours of equal current density. The arc positions are then taken at the odd-numbered contours. In this way each arc carries the current for a surface area on each side of it, closely approximating the effect of a continuous current density. The approximation is valid for magnetic fields calculated at radii less than the cylinder radius minus half the arc separation at that position. 55 In practice this means calculating (o, z) = fz F (0, z) dz - a lo F (0, z) do (9) t 0 c Z- 60 and finding the maximum value of 1,(0j) I, dividing 1. by 2n where n = number of arcs, and setting 4 GB2193322A 4 I 0 (2M - 1).= 1 4, Z) (10) 7n t to determine the arc positions.

In certain cases the are positions may be obtained more easily by integrating the relationship F dz = + a d 10 ( 0, Z) F Z) z which is the differential equation specifying the contours of F directly. This gives, in general F z (, z)F 0 (0,Z) = F 0 (00, z)F z (0, z 0) Q I L) 15 where z. and 0,, are constants.

Specific Cases 20 It is advantageous to consider first the following examples which do not require the full target field approach.

(a) RF coil Here it is well known that a continuous current density 25.25 F z = j2 sin 0 will give a uniform magnetic field in the x direction. An exceptionally good rf coil (Figure

2) can be designed using filamentary currents placed according to the prescription of equation 30 (9). The required current return paths are circular arcs, which can be placed at z = d, sufficiently far from the region of interest for imaging that their effect can be neglected for jzj < d/2. For the straight section (9) gives z) = - a fo I () sino d(p 35 c a = 1 0 COW 40 (10) then gives 1 a 45 -T (2m co,-t m n n 0 12) 50So M = CosM 1/2) 50 Straight line sections placed at these positions, corresponding to equal spacing in the x direction, give a remarkably uniform B,, over a large volume. Generally, if the return path distance d > 2a, the field is uniform to 4% over a cylinder of length d and radius b, where 55 b = a (I Cos-' (1 2n For n = 6, b = 0.8a. 60 By comparison, a typical saddle-type rf coil of similar dimensions, corresponding to n 1 in equation (12), has a volume of similar uniformity of length 0.7d and radius 0.45a. The efficiency of the squirrel cage design, measured as the ratio of B,, to the ampere- turns of wire used,. is smaller by a few per cent than a saddle coil of similar dimensions, though the efficiency of a saddle coil giving an equivalent volume of uniformity is less by a factor of two. Orthogonal 65 GB2193322A 5 arrangement of receiver and transmitter rf coils is achieved easily by rotating the coil design by 7r/2 about the z axis. The squirrel cage design may also be driven in a parallel arrangement, providing a very low effective inductance, and this is useful at high rf frequencies.

(b) Transverse gradient coil: no end-correction An ideal transverse gradient, where 5 BZ x G = p G coso can be achieved using a current density which has m 1 and no variation with z only (cf. 10 equation (8)). Here B ( 1 k 6 (k (14) z 4 15 CL nA F z) a Geoso 2 20 Equation (5a) requires that F IT G cxin(p dz 2p 0 25 -11 F G gino (15) z 2A 0 30 The arc positions may be found using equation (1 l):

Z71 dz add, 35 G --!no - a G coso 2W 2 lic Taking the length of the coil to be d, 40 d 2.2-,A; where oe- = constant, defining CL a contour.

If n arcs are required, the position of the rn'th arc is given by 45 CC50 = d. 2m - 1 16) 2n Normally, to avoid undesirable transverse fields, a pair of such coils would be combined to provide symmetry about the z = 0 plane, as shown in Figure 3. Here the two coils are 50 overlapped to minimise the separation between adjacent wire segments, and hence to maximise the radius of gradient uniformity. The current return paths are taken simply as circular arcs at z = 1.25d. Typically, for six-arc coils, with d/a = 2, the volume over which the gradient is uniform to within 5% is a cylinder of radius 0.8a and length 2a. Increasing the number of arcs will increase the total inductance, if they are driven in series, but will also increase the radius of 55 the uniform gradient region, The wires may alternatively be moved into positions where they can be driven in parallel, the current being then determined by the relative length of the are. This is an efficient arrangement for maximising the volume of uniform gradient. An iterative program has been written to achieve this, a typical coil design being shown in Figure 4. 60 (c) Transverse gradient coil: with end corrections For this coil design it is necessary to specify the desired field B, more carefully to. include explicitly the effect of the current return paths. A combination of the functions 9(z) described in equation (6a) is useful:

6 GB2193322A 6 d d, B z (C'O'Z) z c CORb 57 + + (z 12 - dl d 2 5 This gives a variation with z shown in Figure 5. For Izi < dj/2 the uniformity is excellent along the z axis.

Use of equation (2) then gives, for m 1 d, d, -d, k/2 10 B_(rr,k) = G c [2e sin(n/6 + d k! '13 j 1, 2) 19 W 15 2e- d 2 k12 sin (Tr/6 +dk31+ -d A k_.-d 2 k 2 2d (18) and hence 20 -G d, d_e F kk) 0 61ja(d 2 dj) k K i I(Ka)i (ke) 09) 25 -d,k/2 3 -d 2 k/2 -d,k -d,k sin(-!+ d ky) -2e s i n x [2e d k - 6 1 2 2 2 30 The integral in equation (4) converges for d i gin 1r/6 35 i.e.

d, 2a 40 Note that F-(k -(k) K 1 0 (20) 45 The arcs may be found as before, by setting 2 rn. - 1 Z) (211) 2n 50 with the results shown, for n = 6, in Figure 6. Since for z 0 and 7r/2 the wires are relatively widely spaced, it is preferable to have a larger number of turns n with this design. Given d, = 2a, d. = 3a and n t: 10, the coil gives a gradient which is uniform to 2% within a cylinder of radius 0.8a and length 2a.

(d) Longitudinal gradient coil 55 This coil is a distributed-arc version of the well known Maxwell pair arrangement. The target field Q, is taken as z) = B z (C 09 Z) G (22) 60 + (Z)C, z Then 7 GB2193322A 7 B (0, k) (2edk/2 cOs (,! + dk v3 e-dk) Z 3 3 2 (23) 5 and 1 -ikz(2,-dK/2 costr/3 + dk e-dk) 2 0, 68 a -- k KI(ka) 1 (kc) 10 0 C (24) The shortest coil achievable, using the given functional form of target field, has a length

2a'%/3, which is double that of a conventional Maxwell pair, but a twelve turn version of this coil (Figure 7) gives a linear z gradient which is uniform to 5% within a cylinder of radius 0.7a and 15 length 2.0a, a volume six times larger than a comparable Maxwell pair.

(e) Uniform field coil

The algorithm described lends itself simply to the generation of highly uniform main magnetic fields. In order to provide correction to order (z/d)12 it is necessary to have a coil length d = a/sin 7r/12 = 3.86a. In return for this somewhat long configuration it is possible to achieve a 20 magnetic field which is uniform to 0.2 ppm over 100% of the coil diameter and out to z = a.

The current density required (Figure 8) can be provided using a large number of turns of wire, since for this application a high inductance causes no problems. Here B 25 B_(c, (p, z) = 0 1 - 1 (25) + (zff d and thus 30 B 1Td dK i n (2F-1) iT 12 (2s-, iT] B (O,k) gin[(2!',1)7 + dk coZ 6 -,=I - i r_ 12 (26) 35 F(O,z) =' F(o,z) is found using equations (3) and (4).

(f) Shim coils The normal complement of shim coils provided for modern high quality magnets for spectros- copy or imaging consists of the following: 40 (i) X gradient (ii) Y gradient (iii) Z gradient 45 (iv) XY gradient (V) XZ gradient 50 (Vi) YZ gradient 2 (Vii) z gradient (Viii) Z3 gradient 55 4 (ix) Z gradient W X 2_Y2 gradient 60 Coils specified as in sections (b), (c) and (d) can supply the first order gradients over a large usable volume. In order to provide the remaining gradients using distributed windings, the Aesired magnetic fields may be specified as follows, with the appropriate azimuthal coefficent numbers m listed: 65 8 GB2193322A 8 2SY- B (c'O'Z) c 2 sin20 0 (z) z m 2 cj XZ B Z(C#ORZ) cc c cosq) h(z) 5 M YZ B z (CPO$Z) a c -,in(p h(z) m = 1 10 z 2 B Z(COM) C, Z 2 g(z) in = 0 Q1 z 3 Bz(CROOZ) cc Z3 g(z) 15 m 0 Z4 B z(c,O.Z) z 4 g(z) M 0 20 X 2 -Y 2 B Z(CROOZ) C 2' cos20 g(z) m +2 where - g(z) 1 (z)6. and h(z) 25 1 - + d + The current densities required to generate these fields, and thus the arc positions approximat ing them, can be derived as shown previously. In general these designs can be realised using - 30 many turns of wire, since shim coils are normally not switched, and therefore may be of high inductance. The volume over which these coils give the specified field to within 1.5% is then normally a cylinder of length 2a and radius a.

When specific fields are required in more than one region, the above principles can be applied to produce integrated coil design. 35 For example, consider the case of a screened coil system, where a target field B(r,z) is required in some interior region (r < a) and a second target field of B=0 is required in the exterior region r > b. Two cylindrical surface currents are then necessary, at radii a and b, where b > a.

In this approach we generalise Equation 1 to 40 0 = [bI'(kb)j"ln(k) + aI'(ka)j2,m (k)]K (kr) (27a) m m m 45 kB m (k,r) [bK(kb)j"m(k) + aK(ka)j2 'm(k)JI (kr) (271b) z (t M 50 where B,(k.r) is the Fourier transform of the desired target field jJ, is the surface current density on a cylinder radius b, and j2 is the surface current density on a cylinder of radius a. The Fourier

0 transforms of the resulting current densities can be written explicitly as:

- -- 55 B m (k- j 1 m (k) 1 (28) 0 I(kD) bk I (k.-)[ K' (kb) m K(ka) 0 m m i 1 (ka) m 60 2 11 (kb) 9 (29) j 'mk'k). -b- m 1 m (K) 0 a 11 (ka). j 0 m 9 GB2193322A 9 For a uniform field B along the z-axis, a possible target function is

B (z,r,O) - f(z)g(r,o) - B11 (.A)' [1 - exp -( Z),011-1 (30) Z p q 5 where p and q are assignable constants, 1 m assignable integers and g(r, O) = B is in this case independent of r and 0. By choosing 1 = 3 we ensure that the far field varies as 1 /z. Other types of coil may require specific r and 0 dependence as reflected in g(r,o) together with a different value for 1. For a short solenoid, the axial variation of field would go as z '+m. The constant p is approximately equal to coil radius and q the coil half- length. The procedure is 10 therefore to evaluate f(k) from Equation (30). Equations (27a and b) are then solved for jl(k) and 0 j2(k). These current densities are then inversely transformed to yield the actual current density 0 distributions. For a uniform field gradient along the z-axis, the target function is simply zf(z) with appropriate adjustments in 1 and m.

Although the Bessel functions are used as single hoop response functions, Equation (27) may 15 be generalised to any coordinate system and arbitrary coil and screen geometry. In this case we may write Equation (30) generally as Z 1 0 = B 0 l (k) (k) + B 0 2 (k) j2 (k) (31a) 20 B i ( (k) (31b) )=-Bl Wjl(k), - B' (k) 0 '1 0 2 j2- where Ba,,,(k are the k-space unit current response functions for a line or loop element 25 0, calculated from the Biot-Savart Law and evaluated along a path in real space externally (e) or internally (i) with respect to the loop element.

Further screens may be added. For n active screens a generalised n x n matrix equation is obtained:

30 BM = -3 0 (k)j(k) (32) where 8,(k) is the generalised transfer matrix for the set of screens. B(k) is the specified field or target response vector and j(k) the current density vector. For finite screen and/or coil 35 currents, we also require detB.(k)91-0. So far we have considered specification of the target response vector along one coordinate axis. For specification along two axes B(k) or j(k) and B.(k) become 2nd and 3rd rank tensors respectively. It is therefore possible to completely specify the magnetic fields in n simply connected regions with n closed active screens. For open screens, the target fields can be approximated only over confined regions dependent on the required 40 degree of approximation and the geometry of the screens.

The above description is mainly concerned with screening fields, so that beyond a given radius, the field vanishes. However from the analysis presented the reverse case may also have application. For example, generation of magnetic fields of known spatial dependence external to the coil system structure could be useful in surface coil design and also in the design of 45 magnets, gradient coils and RF coils in oil borehole NIVIR spectrometers.

Claims

1. A magnetic coil including a plurality of conductors positioned on a surface or surfaces in which the positions of the conductors are determined by the calculation of the continuous 50 current density required on that surface to generate a specified magnetic field everywhere on a further target surface or at specified points thereon, and calculation of conductor positions on the former surface which provide a good approximation to this current density.

2. A magnetic coil as claimed in Claim 1 in which the surface or surfaces are cylinders and the target surface or surfaces are also cylinders. 55

3. A magnetic coil as claimed in Claim 2 in which the cylinder on which the conductors are positioned is a cylinder of a non-ferromagnetic material, or the cylinder may be notional, the coils being merely shaped to conform to a cylindrical or partially cylindrical shape.

4. A magnetic coil as claimed in Claim 2 in which the cylinder on which the coil wires are to be mounted is of larger diameter than the diameter of the target cylinder, the coil therefore 60 being designed to produce a known field within the cylinder.

5. A magnetic coil as claimed in Claim 4 in which the target cylinder is one fifth to one tenth the diameter of the cylinder.

6. A magnetic coil as claimed in Claim 4 in which the target cylinder is of zero diameter.

7. A magnetic coil as claimed in Claim 4 in which the target cylinder is of larger diameter than 65 GB2193322A 10 the cylinder.

8. A magnetic coil as claimed in Claim 2 in which the coil design is obtained from calculations carried out on first and second target cylinders, the first being of smaller diameter and the second of larger diameter.

9. A magnetic coil as claimed in Claim 8 in which a plurality of coaxial cylindrical or semi- 5 cylindrical coils are designed to produce various fields specified on target cylinders in the intervening spaces, and within the central cylinder.

10. A magnetic coil as claimed in Claim 9 in which the field within the inner cylinder is uniform and outside the outermost cylinder the field is zero.

11. A magnetic coil as claimed in Claim 9 in which the inner field is zero with a defined outer 10 field. 1V

12. A method of producing a coil wound on a surface to produce a desired magnetic field I when energised including calculating the continuous current density required on the surface to;, generate a specified magnetic field on the surface of a target surface of different size, calculating conductor positions on the former surface to provide a near approximation to this current 15 density and winding the conductors onto the former surface in the calculated positions.

Published 1988 at The Patent Office, State House, 66/71 High Holborn, London WC1R 4TP. Further copies may be obtained from The Patent Office, Sales Branch, St Mary Cray, Orpington, Kent BR5 3RD. Printed by Burgess & Son (Abingdon) Ltd. Con. 1/87.

V1 i