Magnetic field screens

Abstract: The screen is provided by surrounding the coil producing the magnetic field with a set of electrical conductors the cur­ rents within the conductors being controlled in such a manner that the field is neutralised in a specific region of space, the current distribution within the conductors being determined by calculating the current within a hypothetical superconduct­ ive shield which would have the effect of neutralising the field, the current through the conductors thereby being a substitute for the superconductive shield.

Description: This invention relates to magnetic field screens and has application in NMR imaging apparatus. Current carrying magnet coils are used for a variety of purposes in NMR imaging apparatus. Examples include large electro-magnets designed to provide static magnetic fields to polarise nuclear spins, magnetic field gradient coils which superimpose gradients onto the static polarising field and rf transmitter and receiver coils. In many cases the design of a magnet coil is such as to optimise the magnetic field within a desired volume. However the coil inevitably produces an extraneous magnetic field outside that volume, especially relatively close to the coil. In the case of large bore static electromagnets the high fields they generate will produce undesirably strong extraneous fields at distances outside the magnet corresponding to many magnet diameters. Such magnet systems therefore require much free and unusable space around their installation. Stray magnetic fields may be intolerable in hospitals because of iron structures positioned around the installation site which vitiate the homogeneity of the magnetic field Additionally, electronic equipment may not perform well in an environment which has an extraneous magnetic field. Furthermore, most NMR imaging systems utilise rapidly switched magnetic field gradients in their operation. A major problem epecially where super-conductive magnets are used, is the interaction of the gradient field with the magnet itself. Existing attempts to minimise this interaction include the use of conducting metal screening sleeves. However, induced currents in these sleeves or in the heat shield of the magnet decay with uncontrolled relaxation times which make it difficult or even impossible to implement some of the faster and more efficient NMR imaging methods. This is because the decaying current produces image fields superimposed on the desired gradient field. This uncontrolled time dependence introduces phase artefacts which can completely ruin the image. In order to provide adequate access for patients, and to improve gradient uniformity, it is desirable to maximise the diameter of the magnetic field gradient coils in an NMR imaging machine. However, this often causes the coils to be close to other conductors, either the surfaces of cryogenic vessels (in superconducting magnet systems), electromagnetic coil supports (in resistive magnet systems), or ferromagnetic pole pieces (in ferromagnetic systems). When gradients are switched rapidly, as with many imaging techniques, eddy currents are induced in these conductors which then contribute additional field gradients varying in time and potentially very non-uniform in space. Typical time constants for the decay of the eddy currents vary from a few milliseconds to hundreds of milliseconds, depending on the type of main magnet and the specific coil configuration. The commonest solution to this problem is to tailor the input applied to the amplifiers generating the gradient coil currents in such a way that the gradient fields themselves follow the prescribed time variation. The input voltage and gradient coil currents are characteristically over-driven for the initial part of the on-period of the gradient. But this remedy has a major disadvantage. If the gradient coils are placed close to the coupled surfaces, so that the eddy current field gradients may have the same uniformity as the desired gradient, the gradient coils become very inefficient and a large over capacity in the gradient current amplifiers is required, since the 'reflected' fields will be large and in the opposite sense from the desired fields. If, on the other hand, the gradient coils are reduced in size, to avoid the amplifier capacity problem, then the reflected gradient fields will in general be non-linear over the region of interest. Furthermore, in either case there are likely to be reflected fields from more distant conductors in the main magnet structure, each with its distinct time constant and spatial variation. The only effective solution is in some way to reduce the gradient fields to zero at a finite radius outside the coils, so that no undesired conducting structures may be coupled to them. US Patent No 4,595,899 (Smith) which discloses the features referred to in the preamble of present Claim 14 shows an active screening arrangement in which a first complex inner coil structure comprising a first set of coils is actively screened by a second outer complex coil structure comprising a second set of coils. The arrangement uses geometric similarity in the positioning of the first and second set of coils to achieve cancellation of the field produced by the first inner coil structure. Perfect cancellation of the field produced by the first set of coils is however only achieved at a specific distance external to the second set of coils, the field at other points external to the second set of coils not being reduced to zero. US Patent No 3,671,902 (Westendorp) shows an active screening system comprising a primary coil which is screened by a plurality of outer coils which surround the primary coil. The outer or screening coils serve to cancel the flux created by the inner coil at selected points external to the outer or screening coils. The current carried by the outer or screening coils is calculated to provide the desired active screening effect. The present invention is concerned with cancellation of magnetic fields from a primary coil at all points outside a magnetic screen which encloses the primary coil. The inventive feature is that this is achieved by enclosing the primary coil with a hypothetical superconductive surface in which are induced surface currents. These surface currents and the fields at the surface must obey the boundary condition, namely that the normal component of field B z = 0 and the tangential components of field must satisfy the conditions B x = j y and B y = j x . Knowing B x , B y then defines the surface current distribution necessary to achieve perfect screening. The condition B z = 0 says effectively that the infinitesimal flux d  = B z ds through the area element ds is zero. This alone is a necessary but insufficient condition for perfect screening. Westendorp's method is based on the cancellation of flux through a set of hypothetical circles placed exterior of a shielding coil (col. 3, 1. 34-36). This means that flux is cancelled only at particular regions defined by the hypothetical circles. This also means that Westendorp's equation (col. 3, I. 53) cannot be used to cancel the magnetic field simultaneously at all points beyond the screen. It is clear from the equation on col. 3, I. 53 and the references col. 3, 1. 34, I. 57, I. 72-75; col. 4, I. 6,7, I. 31,32, I. 33-35 and I. 49 that the invention in Westendorp is concerned solely with the cancellation of flux through the hypothetical circles. It is an object of the invention to provide more efficient and effective magnetic field screens for coil geometries useful in NMR imaging. It is a further object of the present invention to provide magnetic field screens for any desired component or components of a magnetic field. The present invention provides a method of designing a resistive screening coil for NMR imaging apparatus as claimed in claim 1. The present invention also provides a screen for a magnetic field produced by an electrical coil implementing the above method, as claimed in claim 3. In one preferred arrangement the conductors of the set are regularly spaced apart from each other. They may be connected electrically in parallel and have different values of resistance in order to produce the desired current distribution. In embodiments of the invention the different values of resistance of the conductors may be produced by different thicknesses of the respective conductors or constructing them with different compositions having appropriate values of electrical resistivity. In alternative preferred arrangements the conductors of the set carry equal currents but are spaced apart from each other by different spacings so as to produce the desired current distribution. It is a further object of the present invention to reduce acoustic vibration in MR gradient coils by using the active screening hereinbefore described. In examples of carrying out the above invention a set of conductors are arranged on a cylindrical former and with appropriate spacing can be fed with equal currents. The spacing of the conductors of the set and the magnitude of the currents are calculated using the currents induced in a flat superconductive screen as the starting point. The calculations described herein represent an analytical formulation of the problem, enabling a fully general calculation of the current density in a cylinder required to cancel outside the cylinder the magnetic fields generated by coils inside. The results obtained are firstly applied to passive shielding, using a thick high-conductivity or superconductive cylindrical tube to solve the reflected fields problem without sacrificing gradient uniformity. The conclusions arrived at can then be applied to active shielding; the calculated current densities within the skin depth of a thick cylinder or in a superconducting cylinder are mimicked using a suitable arrangement of a set of wires supplied with currents of appropriate magnitude. For passive shielding using a normal conducting cylinder to be effective, the skin depth δ in the shield at the lowest frequency represented in the particular gradient switching sequence must be much smaller than the thickness d of the wall of the cylinder. For an echo planar switched gradient, for instance, with a fundamental frequency of 1 kHz, this entails a wall thickness of -10 mm of aluminium. For switched gradients such as that used for slice selection, where there is a non-zero d.c. component of the fields at the cylinder surface, passive shielding is not appropriate. When the criterion of δ/d « 1 is met, the time dependence of reflected fields will be identical to that of the applied field, and it is only the spatial non-uniformity of the net magnetic field which is of concern. Analytic calculation of the ensuing fields enables the necessary corrections to coil spacings to be made. By Lenz's law, a magnetic field screen constructed in accordance with the above design criteria produces a magnetic field which opposes the field generated by the primary magnetic coil that it is designed to screen. For a given current in the primary coil the resultant magnetic field generated within the volume embraced by the coil is reduced and its spatial variation is also changed by the presence of the screen currents thus introducing undesirable variations in the primary field. It is therefore a further object of the invention to provide a screening coil arrangement in which the above disadvantages are overcome. Accordingly, the present invention also provides a screen system as defined in Claim 14. In order that the invention may be more fully understood reference will now be made to the accompanying drawings in which:- Figures 1 to 10 are explanatory of the underlying theory; Figures 11 to 17 illustrate in diagrammatic form various embodiments of the invention; Figure 18 shows a cylindrical co-ordinate system used in calculating the magnetic fields produced by current flow on a cylindrical surface; Figure 19 illustrates saddle-shaped coils used as magnetic gradient field coils; Figure 20 illustrates the configuration of one octant of a set of screening coils in accordance with the calculations described herein; Figure 21 shows curves for optimising the positions of the arcs of the saddle coils of Figure 19; Figure 22 illustrates diagrammatically the magnetic field produced by an unscreened primary coil in the form of a single hoop; Figure 23 illustrates the same hoop with a double screen embodying the invention; Figure 24 is a graph showing the magnetic field at different radial positions that is produced by the Figure 23 arrangement; Figure 25 is a perspective view of a double saddle coil used to produce transverse gradient fields which is screened by a double screen embodying the invention; Figures 26 to 33 illustrate ways of providing varying current distribution over a desired area; Figure 34 shows a parallel band arrangement; Figure 35 shows a series arrangement; Figure 36 shows various arrangements for assisting in reduction of acoustic vibration in MR coils by a gradient wire with (a) a wire pair (b) a wire array (c) a pair of conducting plates used to nullify the main field in the neighbourhood of the gradient coil and (d) a double active screen arrangement; Figure 37 shows winding strategies for screening wires about gradient coils of the (a) standard circular design and (b) saddle design; and Figure 38 illustrates the screening of one coil from the magnetic field of another. The basic theory and ideas are developed from the simple case of an infinite straight wire parallel to an infinite flat conducting sheet. Figure 1 shows a long straight wire carrying a current I. The magnetic field B at a point P which is at a distance r normal to the wire is given by B  = ( µ o 4π ) 2I r If the current is changing with angular frequency ω and the wire is near to an infinite conducting sheet as shown in Figure 2 (also if the current is static and the sheet is superconducting), the magnetic field undergoes distortion at the metal surface. Let us assume no field penetration into the screen, i.e. the angular frequency ω and the metal conductivity δ are sufficiently high, then the boundary conditions for the magnetic field at the surface are B x = 0 B y = B y B z = 0 The details of the field at a point P due to a wire at a distance d from the sheet may be conveniently calculated as shown in Figure 3 using the method of images which assumes a wire at a distance d from the other side of the sheet carrying a current -I. In general, current in the sheet surface is directly related to the tangential field H y . The total field B at P, Figure 4, is given by Equation (1) and may be resolved into components B x and By in the x and y directions respectively B y = ( µ o 4π ) 2l r cos B x = ( µ o 4π ) 2l r sin Taking the image current -I into account, Figure 5, we obtain the total y direction field B TOT at point P for a current I at point 0 and its image current -I at a point 0'. B TOT = B y + B image y We also note from Figure 5 that y = r sin = r' sin and r cos + r cos = 2d If P is on the conducting sheet surface r = r' in which case, see Figure 6. B y = ( µ o 4π ) 4l r cos r 2 = d 2 + y 2 and cos = d/r From Equations (9) and (10) we obtain B y = ( µ o 4π ) 4I d d 2 + y 2 . Consider now the line integral of the magnetic field in and close to the metal surface, Figure 7. By Ampere's theorem we have where j y is the current density per unit area and ds an area element and dl a path element. For a short path I parallel to the surface Hy is constant. The line integral is therefore H y l - H ' y l + (H x + H ' x )dl = J y l where J y is the surface current density per unit length. But in the metal H' y = 0 and H x = H ' x yielding for dl→0 H y = J y The surface current distribution is therefore J y = ( 1 4π ) 4ld d 2 + y 2 = J 1 (d) This function is plotted in Figure 8. To determine the field produced by the surface current density distribution, let us assume we have a surface distribution Jy within a flat metal sheet as in Equation (15). Consider an element of surface dl with current δi, Figure 9. This current is given by δi = ( 1 4π ) 4ld d 2 + y 2 dy The elemental field at point P which is a distance r from the element and distance d from the sheet is δH = δi 2πr , The tangential component of which is:- δH = δi 2πr sin = δi 2π d r 2 = 2ld 2 π ( 1 4π ) dy (d 2 + y 2 ) 2 The total field is which yields B y = I d ( µ o 4π ) Thus the field at point P distance d from the surface current distribution is equivalent to a mirror current of -I at distance 2d from P. The total surface current The results derived above suggest that instead of using a metal plate to screen alternating fields, an active screen comprising a mesh of wires may be used in which a current pattern is externally generating to mimic a desired surface distribution. This situation is shown in Figure 10. Wires at positions y 1 , y 2 ....y n in the y direction all carrying an equal current I are spaced at different discrete intervals corresponding to equal areas A under the J i curve of Figure 8. For equal currents in the wires, the wires must be unequally spaced with positons y n which satisfy the relationship for an even array of 2N wires spread about the y origin. For an odd array of 2N + 1 wires with one wire at the origin we have where n = 1, 2, 3... Alternatively, wires may be equally spaced with different currents chosen such that In either case, since we have arranged to satisfy the original boundary conditions on the surface of a fictitious plate, the magnetic field in the half plane (x,±y) approaches zero. The degree of screening depends ultimately on the number of wires used in the active screening mesh. An example of such an active screen is shown in Figure 11 for a distribution of current carrying conductors corresponding to Figure 10 all carrying an equal current I. The field B on the opposite side of the screen is effectively zero. If two parallel screens are used, each will have primary current distributions of J 1 (d 1 ) and J 1 (d 2 ) given by Equation (15). This is shown in Figure 12. However, each induced current distribution will induce further changes in the distribution in the opposite plate represented by additional terms J n (d n ). This is equivalent to a set of multiple reflections, Figure 13, which correspond to an infinite set of images. When d 1 = d 2 = d, image currents occur at x = ±2nd, n = 1, 2... The total induced surface current in each sheet is the sum of J(d) =∑ n J n (d n ). The boundary conditions at the metal surface ensure that the normal laws of reflection apply. However, when active current screens are used, the reflection laws may be selectively changed to reflect any images using the J n (d)'s corresponding to particular distances d. An example of a gradient active screen is shown in Figure 14 for a single circular hoop of diameter 2a. Let this be screened by an active current mesh in the form of a cylinder of radius a + d. In a metal cylinder, image hoops appear at radii r = b + d, 3b±d, 5b±d, etc. However, since the effect of these distant images diminishes quite rapidly, it is reasonable to approximate a screen with J 1 (d) corresponding to the plane sheet case, Equation (15). Better approximations may be obtained by an iterative numerical approach. Although exact solutions for the surface current in a cylinder exist, when actual wire screens are constructed, the numerical approach is preferred since it automatically takes account of the finite number of wires and their discrete spacing. The screen current is N 2 I 2 = - r 1 r 2 2 · α N 1 I 1 where N 2 is the number of screen wires in the mesh of radius r 2 each carrying a current 12 and N 1 is the number of turns in the primary coil radius r 1 carrying a current I 1 . The factor a is of the order of unity and is chosen to optimise the screening. The whole optimisation procedure is accomplished by a computational process which generalises to applying minimisation of the total field over a limited region of space. Mathematically this is conveniently achieved by a least squares method. For practical ease, it is desirable to have both coil and screen driven from the same current source. Since the total screen current is less than the primary current it will in general be necessary to take parallel combinations of the screen mesh so that the total system may be driven in series. However, a parallel arrangement is also possible in which the screen wires are varied in resistance and/or impedance in such a way that when driven from a voltage source appropriate currents flow. Because of the screening effect the inductance of both versions of the coil system should be low. Sets of screened hoops may be used to construct screened magnets producing uniform magnetic fields. The presence of the screen around one hoop is approximately equivalent to an image hoop with negative current producing an opposing field. For a Helmholtz coil pair the intercoil spacing equals the coil radius a. When screened, however, the spacing must be changed so as to make the second derivative of the field B z with respect to z vanish for the combined coil system. Similarly for a Maxwell pair designed to produce a linear magnetic field gradient, the intercoil spacing is ideally a √3 . This is shown in Figure 15. The two hoops forming a Maxwell pair are screened by a concentric pair of screening meshes offset axially from each other. The combined screen current distribution is also shown. Again, however, when screened coils are used, a new spacing obtains which makes the third derivative of B z with respect to z vanish for the total coil system. The process of optimisation of coil geometry can be done analytically for simple coil structures as discussed above. For more complicated systems such as screened saddle geometry gradient coils, it is preferable simple to find by computational means the position where the chosen derivative vanishes, or is minimised. For NMR imaging sytems using superconducting magnet coils it is convenient to use saddle coils to produce the transverse gradients ∂B z ∂x and ∂B z ∂y . In some imaging techniques, at least one of the gradients can be very large making interaction with the main magnet potentially serious. Figure 16 shows one half of a G x screened gradient coil system. A screened single saddle coil is shown in end view at (a) and plan view at (b). The dotted lines correspond to the screening mesh. To a first approximation, the screen current profile is J 1 (d). Better screening may be obtained by an iterative procedure which minimises the field outside the screen. Figure 17 shows a screened rectangular G x gradient coil set with 1/2 screen. Again, if a d«2a, J 1 (d) may be used as a good approximation for the screen current distribution. For better results other reflections can be included or the iterative procedure used to minimise fields outside the coil. NMR imaging systems require rf coil systems to deliver rf pulses to the specimen and to receive signals induced in the sample. Because of the number of other coil systems required for field gradients, there is always a problem of space. With normal rf coil arrangements the field outside the coil drops off rather slowly. In order to minimise coil interactions which can lower the Q value, rf coil diameters are often chosen to be around 0.5 to 0.7 of the gradient coil diameter. With shielded rf coil designs, it may be possible to utilise more of the available space without loss of performance. A systematic procedure for reducing extraneous magnetic fields outside the active volume of static magnets, field gradient coil systems and rf coils has been described. In NMR imaging, reduction of stray fields in all three types of coil structure is extremely important. The method utilises active magnetic screens and has the advantage that such screens operate independently of frequency down to dc. Some price is paid in terms of reduction of field in the active volume compared with that of the free space value. With time dependent gradients, the price is in general acceptable since for fast NMR imaging schemes employing rapid gradient switching, active coil screening may be the only way in which such imaging schemes may be made to operate in the relatively close confines of an electromagnet. In what has previously been described, iterative and least suqares approximation methods are used to obtain actual screening wire positions. It is possible to obtain these positions directly using analytical methods. If the gradient coils and screens are located on cylindrical formers, it is natural to use cylindrical co-ordinates p, , z in order to retain the symmetry of the system. The z axis is taken to lie along the axis of the cylinder as shown in Figure 18. The vector potential A, is used to describe the magnetic field. This has the components A p , A ø , A z given by A ρ = µ o 4π ∫ J ' ( r ) dv' sin( - ') | r - r '| A  = µ o 4π ∫ J ' ( r ') dv' cos( -') | r - r '| A z = µ o 4π ∫ J z ( r ') dv' | r - r '| where J is the current density and dv' is a volume element corresponding to the position vector r '. There is no current flow in the radial direction in many problems of interest so J has only z and  components. It has been assumed that currents induced in the shield are confined to the surface of a cylinder of radius b. The gradient coils to be shielded are mounted on a cylindrical former or radius a which is coaxial with the shield. The currents can then be written as J = F (z, ) δ (ρ - a) + f (z, ) δ (ρ - b) where F describes the current in the gradient coils, f describes the current induced in the shield and δ is the Dirac delta function. It is possible to derive relationships between F and f from the condition that the radial component of the magnetic field is zero on the surface of the shield. The other constraint which is used is the equation of continuity, which, in the absence of charge accumulation is ▿ . J = 0 The vital stop in the analytic treatment of this system is the use of the Green's function expansion where m is an integer and ρ < (ρ > ) is the lesser (greater) of p and p', and I m (z) and K m (z) are modified Bessel functions. To use equation (30) in equations (25), (26) and (27) it is helpful to define a type of Fourier transform of f and F as follows: where the quantities F m z (k) and F m  (k) are defined in an analogous way. The components of A thus become (e.g. for ρ > b): Similar expressions can be obtained when ρ < a or a < ρ < b. Now let the boundary condition that the radial component of the negative field B ρ = 0 at ρ = b be applied. This is equivalent to 1 b ∂A z ∂ ρ = δ = ∂A  ∂z ρ=b Equations (33) and (34) are now used and terms varying as e im are equated. This gives 2m bk k m (kb) [bI m (kb)f m z (k) + aI m (ka)F m z (k)] = f m  (k)b[I m-1 (kb)K m-1 (kb) + I m+1 (kb)K m+1 (kb)] + F m  (k)a[I m-1 (ka)K m-1 (ka) + I m+1 (ka)K m+1 (ka)] This equation can be simplified using the relations derived from the equations of continuity, 1 b ∂F  ∂ = - ∂F z ∂z for the currents in the shield, and 1 a ∂F  ∂ = - ∂F z ∂z for the currents in the gradient coils. These equations are equivalent to f m  (k) = - kb m f m z (k) and F m  (k) = - ka m F m z (k) The recurrence relations for Bessel functions are also used, from which can be derived the identity I m-1 (Z 1 ) K m-1 (Z 2 ) + I m+1 (z 1 ) K m+1 (z 2 ) + (2m 2 /z 1 z 2 ) I m (z 1 )K m (z 2 ) = - 21 m '(z 1 ) K m '(z 2 ) Here the prime denotes the derivative. Equations (37), (40), (41) and (42) can be combined to obtain the elegant expressions, f m z (k) = - F m z (k) a 2 I m '(ka) b 2 I m '(kb) f m  (k) = - F m  (k) a I m '(ka) b I m '(kb) These identities provide the means of calculating fields due to any combination of currents constrained to flow on the surface of a cylinder inside a conducting shield. This will now be illustrated with an example. For transverse field gradients it is common to use a saddle coil configuration such as shown in Figure 19. Two pairs of saddle coils are shown. One pair extends between -d 1 and -d 2 along the z axis. The other part extends between + d 1 and + d 2 along the z axis. With energisation of the coils as marked a gradient field is produced in the form of a magnetic field along the z axis which has a gradient in the x direction. The field produced by such coils extends widely outside the cylinder on which the coils are wound unless they are shielded by a conductive sleeve. For the standard coil geometry with 120° arcs for the saddle coils. F  (,z) = I{δ(z-d 1 ) + δ(z + d 1 ) - δ(z-d 2 ) - δ (z + d 2 )} {H( + π 3 )(1 - H( - π 3 )) - H( - 2π 3 )(1 - H( + 2π 3 )) - where H(x) is the Heaviside step function. This has the Fourier transform F m  = 2sin(mπ/3) mπ I (cos(kd 1 ) - cos(kd 2 )) x (1 - e imπ ) This is zero for m even or an integer multiple of 3. This leaves non-zero terms for m = 1, 5, 7, 11 etc. The position of the screening arcs can now be determined as follows. First it is necessary to specify the continuous current distribution to, which the arcs approximate as set out in equations (43) and (44). The actual surface current density may then be written and The next step is to determine the stagnation point of this current density, that is, the point S = (O, d) at which f  and f z are both zero, around which the surface current flows. By symmetry this occurs at  = 0, and since f z (0, z) = 0 for all z, it may be found by solving f  (0, z) = 0, by successive approximation. Having found S, the integrated surface current I t crossing the line SP between S and an arbitrary point P on the cylinder (Figure 20) is given by (using the surface version of the divergence theorem). Hence The contours of the induced surface current may now be found by setting I t (, z) = constant. These may be translated into positions of screening arcs in the following way. The total current in the cylinder I t -(0,0) is divided by 2N, where N is the number of screening arcs required, and the M'th arc consists of the contour where I t (,z) = I t (0,0) 2N (2M-1) In practice I t (, z) is calculated over a grid of 50 x 45 points, and the contours are found by linear interpolation between the calculated points. To check that these arcs indeed provide adequate screening of saddle coil fringe fields, the total field can be calculated using the Biot-Savart Law, taking as line elements dl the intervals between successive calculated co-ordinate pairs along each arc. Given a set of saddle coils with radius 0.31 m and arc spacings from the centre d 1 and d 2 of 0.108 m and 0.404 m, the shielding produced at 0.55 m radius, using six screening arcs on a cylinder of radius b = 0.45 m, is as follows: Maximum unscreened field = 0.6 x 10 -6 T/A-turn Maximum screened field (z < 0.5 m) = 0.2 x 10 -7 T/A-turn Maximum screened field (z < 1.0 m) = 0.47 x 10 -7 T/A-turn For comparison, the maximum field at 0.55 m radius produced by an unscreened small scale saddle coil set (a = 0.16 m, d 1 = 0.56 m, d 2 = 0.206 m) is 0.86 x 10 -7 T/A-turn. Figure 20 shows the configuration of one octant of the set of screening coils calculated above. A larger number of screening arcs, or the use of foil rather than wire conductors, will further reduce the fringe fields. Screening a set of saddle coils, the spacing of which has been optimised without the screen present, inevitably reduces the uniformity of the gradients produced. The uniformity may be recovered however by adjusting the arc spacing as follows: The z component of the magnetic field can be derived from equations (34), (35), and (44) to give This becomes, with substition of equations (45) and (46) It is now possible to optimise the gradient linearity by adjusting the arc positions of the saddle coils. The terms for m = 5 are of fifth order in x, y or z, whereas there are terms for m = 1 which are of first order in x and of third order in x, y and z. The optimum choice of the parameters d 1 and d 2 is the one which removes the third order terms. This gives the condition There are now two unknown quantities, D 1 (=d 1 /a) and D 2 (=d 2 /a) and only one constraint, so it is not possible to give a unique choice of parameters. However it can be ensured that each parameter separately satisfies the equation (with α = b/a) for this automatically satisfies equation (54). Values of D 1 and D 2 as a function of a are shown in Figure 21. It is possible to improve on this arrangement by altering D 1 and D 2 slightly (subject to the constraint given by equation (55)) and minimising the fifth order terms. The final choice of values of D 1 and D 2 depend on whether the x or z variation is considered to be more important. The values of D 1 (α) and D 2 (α) shown represent an excellent starting point in the search for the optimum position of the saddle coils. Referring now to Figure 22 there is shown therein a magnetic coil in the form of a single wire hoop 1 of radius a carrying a current + I. In Figure 23 the same wire hoop 1 is surrounded by two active magnetic screens S1 and S2. Each screen comprises a set of electrical current carrying conductors but for simplicity the screens are shown as sections of cylinders. Outer screen S1 is a cylinder of radius b and inner screen S2 is a cylinder of radius c. With appropriate screen current density distributions screens S1 and S2 act together as a flux guide confining the field lines in the manner indicated. The design criteria for the current density distributions are set out below. For a single conducting screen S, the boundary conditions of the magnetic field B(r-r') at the surface of the screen require only that the axial component B z (r-r', z'z') is considered which for a coaxial hoop is independent of azimuthal angle . For perfect screening it is required that I B z (r-a,z) + B s z (r-b,z) = 0 for r > b and for all z, where B z (r-a,z) is the primary hoop field per unit current and B s / z (r-b,z) is the total field generated by the screen. The screen field is the convolution of the surface current density j  (z) with the hoop magnetic field response per unit current. Equation (56) may therefore be written as Note that F m  (q)= J m  (q); F m z (q) = J m z (q) f m  (q) = j m  (q) ; f m z (q) = j m z (q) where q may be taken to be k or z. The spatial Fourier transform of equation (57) gives I B z (r-a,k) = -B z (r-b,k)j  (k) where k is the reciprocal space wave number. Equation (58) may be generalised for two screens with current densities j 1 /  (k) and j 2 /  (k) as set out below. For zero field in the range r ≧ b and the unperturbed hoop field for r ≦ c there is obtained I B z (r-a,k) = -B z (r-b,k)j 1  (k) - B z (r-c,k)j 2  (k) for r ≥ b O = B z (r-b,k)j 1  (k) + B z (r-c,k)j 2  (k) for r ≤ c The current hoop fields are evaluated numerically along the z-axis on the appropriate cylindrical surface using the Biot-Savart Law, then Fourier transformed to k-space. This allows numerical solution of the simultaneous equations (59) and (60) yielding the k-space current densities. These are then inversely transformed to yield the actual screen current density distributions. In this example consider a primary hoop of radius a = 0.5 m shielded by two active screens S1, S2 with radii b = 0.75 m and c = 1.0 m respectively. The hoop current is 1A. Using the computed distributions, the total magnetic field generated by the screen system in the hoop plane (z = 0) is calculated as a function of r. This is shown in Figure 24 which is a graph of a calculated z-component of magnetic field B z (r,0) against radius in the plane of a doubly screened flat current hoop of dimensions as in Figure 23. The unscreened hoop field is shown by the broken curve and is equal to B z for r < 0.75 m. As expected for r < c, the field B z exactly equals that of the unscreened hoop. For r > b, the field B is zero. Between the screens, the field is wholly negative. The total screen currents are given by I 1 = j 1  (0) and I 2 = j 2  (0). In this case I 2 = -I 1 = 0.57 I. The results have been produced by numerical methods using a computer. The inner screen S2 behaves like a complete absorber of the primary field. However, once trapped between the screens, the field is completely internally reflected by S1 and S2. The screen S2 behaves like a perfect one-way mirror. Practical screens having these properties clearly cannot be continuous metal surfaces. Externally driven discrete wire arrays which approximate the calculated continuous surface current densities are used instead. For a common screen current, wires are placed at positions corresponding to equal areas under the current density distribution integrated over discrete intervals, as described in the above patent specifications. Although active screening of a single hoop has been described in the above example, the principles of double screening apply to other geometries, for example flat screens as well as more complex coil geometries, some of which are used to produce linear gradients in NMR imaging. On the basis of the analytical expression for the component of the magnetic field B z which is set out in equation (52), equations (59) and (60) become, for k-space. aII 1 (ka) = -bI 1 (kb)j 1  (k) - cI 1 (kc)j 2  (k) , for r > b O = bK 1 (kb)j 1  (k) + ck 1 (kc)j 2  (k), for r < a which yield the current densities and J 2  (k) = -j 1  (k) b c K 1 (kb) K 1 (kc) For k=0, K 1 (kb)/K 1 (kc) = c/b and I 1 (ka)/I 1 (kb) = a/b from which it is deduced that the total currents I 1 , I 2 flowing in S 1 and S 2 are equal and opposite. From equations (63) and (64) there is obtained = -I 2 If a = c thereby making the primary coil and S 2 coincident, the condition b = a√2 makes I 1 = -I = = I 2 and means that the coil and screens may be placed in series. For discrete wire arrays chosen to approximate the required continuous current distributions, equations (65) and (66) become and N 1 I 1 = -N 2 I 2 where N o I o , N 1 I 1 and N 2 I 2 are the ampere-turns for the primary coil and screens S 1 and S 2 respectively. A more general series coil arrangement is possible by varying both the turns and screen radii in equations (67) and (68). A Maxwell pair of two hoops with opposite currents may be used to generate a z-gradient field. Since double screening produces the free space field for the hoops in the region r < a, the usual coil spacing a 3 obtains for the most linear gradient along the z-axis. Figure 25 shows a simple saddle coil of radius a used to produce a transverse gradient. Let this be shielded by two active cylindrical screens S 1 and S 2 with radii b and c respectively where a ≦ c < b. For a standard saddle geometry with 120° arcs, the primary current is J  (,z) = I{δ(z-d 1 ) + δ(z+d 1 ) - δ(z-d 2 ) - δ(z+d 2 )} × {H(+π/3)[1-H(-π/3)-H(-2π/3)[1-H(+2π/3]} - where H(o) is the Heaviside function. The Fourier transform of equation (69) is J m  (k) = 2sin(mπ/3) mπ I[cos(kd 2 ) - cos(kd 1 )](1-e imπ ) This is zero for m even or an integer of 3. Since there are now many values of m, equations (61) and (62) are generalised to give for two shields J m  (k)aI ' m (ka) =-bI' m (kb)j 1,m  (k) - cI ' m (kc)j 2,m  (k) 0 = bK ' m (kb)j 1,m  (k) + cK ' m (kc)j 2,m  (k) From equations (71) and (72) there is obtained the results that and These results evaluated at k = 0 give for each separate arc pair at ±d 1 and ±d 2 the total azimuthal screen currents I 1 and I 2 . The dominant components of these currents arise from the m = 1 terms and may be simplified by noting that at k = 0, I 1 '(ka)/I 1 '(kb) = 1 and K 1 '(kb)/K 1 '(kb)/K 1 '(kc) = c 2 /b 2 . The z-components of current flowing in the screens may be calculated by noting that div j ] r = R = 0 which expanded gives ∂j z /∂z= - 1 R ∂j  /∂ Fourier transforming equation (76) gives j m z (k) = - m Rk j m  (k) Equation (77) for R = b or c and equations (73) and (74) give on transforming to z-space the actual screen surface densitites. The above results have shown that by introducing a second active screen the spatial response within a primary coil can be made to be independent of the surrounding screens. The inner screen may be positioned to be coincident with the primary coil and still remain the above advantages. While in the above description two active magnetic screens have been employed it is possible to extend the principles of active magnetic screening of coil structures to multiple screens. Use of two or more than two screens has the advantage that the screens can be designed so as not to vitiate or change the character of the magnetic field spatial response from the primary coil structure being screened. This is true even when the inner screen of a two or multiple screen structure is coincident with the primary coil structure. The calculations and analytical expressions presented refer to continuous current density distribitions in the screens. Practical active screens require discrete wire arrays which simulate the continuous current density distributions. Discrete screens also allow exploitation of the selective transmissive and reflective properties of the active magnetic shields. The object of a magnetic screen is to provide a spatial current distribution which mimics that which is theoretically induced in a real and/or fictitious continuous conducting sleeve around a coil structure in which the coil itself is producing a time-dependent magnetic field gradient. Equation 24 describes an arrangement in which wires are equally spaced by with currents chosen again to mimic the induced surface current distribution in a continuous metal screen. Several methods of varying the current in these conductors are:- 1. To include in each conductor a small resistor chosen to give required current. 2. To change the diameter or the shape of the wire so as to affect its resistance in the right manner. 3. Change the composition of the conductor to affect its resistance. These situations are covered in Figures 30 and 31. It will also be clear from the above discussion that when the conductor size is changed as in Figure 31 and in the limit of uninsulated touching wires, we have the situation shown in Figure 27. An alternative arrangement to this is a profiled cross-sectional band of material as shown in Figure 26. Alternative ways of producing this band or its effect are shown in Figures 28 and 29. The thickness t of the band must be chosen such that the electromagnetic penetration depth δ is less than t for the highest frequency present in the current switching waveform. It will be seen that the above conductor arrangements 26 to 31 could be parallel arrangements or bands as in Figure 34 fed with appropriate total current along the edges indicated. However, it is conceivable that one could manufacture multi turn band structures producing a series arrangement as in Figure 35. An alternative approach using a constant standard wire section employs the wire arrangements in Figures 32 and 33. Here the wires are stacked in such a manner as to produce the desired current distribution. If the wires are uninsulated and touching this is an alternative method of producing the equivalent arrangement of Figure 26. However, if the wires are insulated and touching it will readily be seen that all the turns may be in series. A series arrangement will have a much higher inductance but require only a small common current through each turn. The parallel arrangements discussed earlier require the driver circuitry to provide the total field screen current for a one turn arrangement as in Figure 34. In order to obtain images by magnetic resonance it is necessary to switch gradients rapidly. For typical magnetic gradient strengths the current I required is commonly around 150 Amps. These currents, when flowing in the static field B produce a force per unit length of F = I x B on the wires carrying this current. For parts of the gradient coil the field B and the current I may be perpendicular thus maximising this force. The resultant motion of the wires causes acoustic noise which can be very loud if strongly coupled to the coil former. This problem is growing in severity with the use of higher static fields. The situation is further exacerbated by use of larger gradient coils where the wire length is greater over which the force can act. Rapid imaging strategies can also create more noise. Solutions to the noise problem such as embedding the wires in rigid materials like concrete do help by lowering the natural resonances of the coil former and by absorbing some acoustic energy. Surrounding the wires with some soft acoustic absorbing material such as cotton wall can also reduce the noise. But these approaches treat the symptoms rather than the cause. In a further example of the use of the present invention the problem is solved at source by reducing the B field to zero in the vicinity of the wires. This eliminates the force on and hence motion of the wires. As well as solving the noise problem, the lack of motion of the wires removes the possibility of progressive fatiguing and fracture of the conductors. In order to do this the wires on the gradient coil are locally screened magnetically from B. This is achieved by using the principles of active screening. A single screen arrangement is created by arranging the screening wires in series. Different arrangements are shown in Figure 36. With these arrangements the far field, i.e. coordinates x, z for a point P a distance R » a where a is the half separation of the screen pair (Figure 36), is effectively unperturbed. For an infinite straight wire screen the interior and exterior screening fields B c and Bp (R » a) are respectively B c = µ o I πa B p = µ o I 2π a R 4 (z 2 -x 2 ) where R 2 = x 2 + z 2 For small separation B c can be very large and Bp small. This represents a net fall off rate which goes as 1/R 2 . For finite wires the fall off is basically dipolar i.e. goes as 1/R 3 . For parallel infinite sheets the fields are B c = µ o J B p = 0 where J is the current density per unit length. For the arrangements sketched in Figures 36a to 36c respectively the screening efficacy will lie between the cases covered in Equations (78) and (79). Because the screen produces a static field, perturbation of the main field can be eliminated with a shim set. The screen or counter field generator need be active only during the experimental period. To further reduce the extraneous static fields produced by the wire screening arrangement, a second active screen may be used in addition to the first screen. Such a double active screen arrangement is sketched in Figure 36. The wire W is roughly screened by a parallel plate arrangement S 1 . Residual far fields arising from non-cancellation outside the plates is annulled by a second active screen S 2 comprising a set of conductors distributed appropriately. Figure 37 shows two local wire screening strategies for a circular hoop (37a) and a saddle gradient (37b). Using the wire screening arrangement described, the benefits of rapid switching of large gradients may be obtained within a high static magnetic field without acoustic noise. Screening could also be applied to a Hall probe to increase its sensitivity to small field variation by removing the main central field. In addition to the application of active magnetic screening to electromagnet and gradient coil design there are a number of possible applications in rf coil design. For example a fully screened series wound rf coil placed coaxially inside an unscreened coil would have zero mutual inductance. However an NMR sample placed inside the screen coil would sense fields generated by both coils. The coils though coaxial therefore behave electrically as though they were orthogonal. This can have advantages in multi-nuclear irradiation and detection. A schematic example of one such possible configuration is shown in Figure 38 where an inner coil C1 is wound inside an outer coil C2. The inner coil C1 can be screened by providing a suitably wound screening coil CS between the two coils C1 and C2. The screen coil CS is as shown connected in series with coil C1 to pass the same current. It is generally longer than the coils C1 and C2 to provide efficient screening and is wound such that the current in coil CS opposes that in coil C1. The exact positioning of the wires in coil CS is determined by using the above approach.

Claims: A method of designing a resistive screening coil for NMR imaging apparatus for selectively screening the field of a magnetic coil comprising calculating the induced current distribution in a hypothetical continuous superconductive metal surface positioned in the place of wires of the screening coil and then calculating the position and current distribution within the said wires to approximate to said induced current distribution so that the screen acts as a complete reflector of magnetic field. A method as claimed in claim 1 using a further coil in combination with the screening coil, which further coil is designed to selectively reflect desired components of the magnetic field. A screen for a magnetic field produced by an electrical coil implementing the method of claim 1, said screen comprising a set of electrical conductors; said screen further including means for supplying the conductors of the set with electrical currents, wherein the electrical conductors of said screen are so positioned and the electrical currents supplied to the conductors are approximately of such a magnitude to satisfy the conditions such that the magnetic field normal to the screen is zero and the tangential magnetic field components are equal to the respective orthogonal current densities in the screen so that the magnetic field on the side of the screen away from the electrical coil is substantially zero at all points. A screen as claimed in Claim 3 wherein the screen current distribution localised to the surface of the hypothetical superconducting surface is determined by the deconvolution of the magnetic field response function of the unit line elements of that current with the field to be screened; such problems being solved in reciprocal space, which is defined by those co-ordinates conjugate to real space used in an appropriate integral transform. A screen as claimed in claim 3 wherein the conductors of the set are regularly spaced apart from each other. A screen as claimed in claim 3 wherein the conductors of the set are connected electrically in parallel and have different values of resistance in order to produce the desired screen current distribution. A screen as claimed in claim 6 wherein the different values of resistance of the set of conductors are produced by different thicknesses of the respective conductors or by constructing the conductors with different compositions having appropriate values of electrical resistivity. A screen as claimed in claim 3 wherein the conductors of the set carry equal currents but are spaced apart from each other by different spacings in such a manner that the current in one wire is equal to an integrated incremental theoretical superconductive surface current distribution, such incremental values being equal for each wire so as to produce the desired screen current distribution. A screen as claimed in claim 8 wherein the number of screen wires is even and corresponds to the even arc integrated surface current distribution contours. A screen as claimed in claim 8 wherein the number of screen wires is odd and corresponds to the odd arc integrated surface current distribution contours. A screen as claimed in claim 3 wherein the screen surface current distribution Jy for an infinite flat screen, screening an infinite straight line wire, is defined by:- A screen as claimed in claim 3 wherein the primary magnetic field is created by a current F in the electrical coil and in which the current f defines the current in the screen and in which the Fourier transforms of the current components of f are defined in the case of coaxial cylindrical coil and screen geometry as:- and wherein the quantities F m / z (k) and F m /  (k) being the components of the Fourier transform of F are defined in an analogous way. A screen as claimed in claim 12 wherein the relationship for the Fourier components of the currents induced in the screen to those in the gradient coils are defined as f m z (k)= - F m z (k) a 2 I m ' (ka) b 2 I m ' (kb) f m  (k) = - F m  (k) a I m ' (ka) b I m ' (kb) A screen system for a magnetic field created by a coil, the coil being surrounded by two or more active magnetic screening coils through which current is passed, each respective screening coil comprising a set of electrical conductors and means for supplying the conductors of the set with electrical currents characterised in that the system comprises an inner screen and an outer screen, the inner screen lying between the coil and the outer screen, the magnitudes of the electrical currents are such that there is no appreciable magnetic field outside the outer screen and the field within the inner screen substantially corresponds to the field that would be provided by the coil if the screens were not present. A screen as claimed in claim 3 wherein around the screen is a second unscreened coil such that the mutual inductance between the coil and screen and the second coil is zero, which total coil system is supplied with currents at rf frequencies as an NMR transmitter/receiver orthogonal coil system in which the coil and the second coil may be independently tuned to a different frequency.

Inventors: BOWLEY ROGER MALCOLM, CHAPMAN BARRY LEONARD WALTER, MANSFIELD PETER, TURNER ROBERT, BOWLEY, ROGER MALCOLM, CHAPMAN, BARRY LEONARD WALTER, MANSFIELD, PETER, TURNER, ROBERT
Applicants: BTG INT LTD, BTG INTERNATIONAL LIMITED
Application year: 1986
Main classification: H01F27
IPC codes: G01R33, G02B17, A61B5, A61B10, H05K9, H01F7, G12B17, H01F27