|
GRADIENTS
WITH INTEGRATED SHIMS - BFG-S SERIES
GRADIENT
DESIGN PRINCIPLES
Eddy
currents are generated by the interaction of the fields external to the gradient
assembly due to the pulsed gradient coils with the external electrically
conductive structures. The
conventional method for dealing with the eddy current problem relies on compensating
for the eddy fields. The
second, preferred, method deals with the problem by eliminating
the eddy currents that generate the eddy fields.
. This technique reduces the
fields outside the pulsed gradient coils substantially to zero. Since the surrounding structure is not exposed to changing
magnetic fields, no eddy currents are induced and no eddy fields are generated.
Pulsed
gradient coils employing this technique are referred to as actively shielded
gradient coils, (ASGCs), because the field that they generate is shielded from
their exterior at all times. This
is achieved by winding an additional set of coils, (the shields), at a radius
larger than the primary windings, (the coils).
The current in the shields is of polarity opposite to that in the coils
and the winding distribution is such that beyond the radius of the shields, at
all points in space, the sum of the field due to the coils and the field due to
the shields is zero. This means
that all of the magnetic flux generated by the coil/shield system that
penetrates the coil returns in the annular space between the coil and the
shield.
It
should be noted that the opposite polarity current in the shield that is
required for cancellation of the coil's external field, creates an internal
field opposite to that of the coil. Thus
the act of shielding the outside makes the ASGC less efficient than an
unshielded pulsed gradient coil, and the closer the shield and coil approach
then the less efficient will be the ASGC. If
they are infinitesimally close to touching then almost equal and opposite
currents in almost identical winding patterns will cancel the outside field but
of course will also almost cancel the internal field.
Then an almost infinite current will be required to produce a finite
field gradient strength. At the
other extreme, as the shield gets farther from the coil, less and less shield
current is required and in the limit it goes to zero leaving the situation as
for an unshielded coil.
COIL
EFFICIENCY
The
most important requirements of an MR experiment that are relevant to the
performance of pulsed gradient coils are: gradient strength, pulse rise time,
gradient linearity and duty cycle. The
specific requirements vary widely from experiment to experiment but there is
always a need to maximize gradient strength, linearity and duty cycle, and
minimize rise time within the constraints of available gradient amplifier power
and cooling capability.
High
gradient strength is required in many applications, including situations where
increased spatial resolution is necessary as in MR microscopy or when very short
imaging times are desired, as is the situation for echo-planar imaging
techniques, or when it is necessary to move through k space in unusual patterns.
High gradients are necessary even for more mundane situations when
minimization of imaging time is desired. For
example, in spin echo imaging the time for a pulse sequence may be shortened by
modifying the phase compensation pulses that occur after the slice select pulse
and before the phase encode pulse. In
such pulses the flat top time may be reduced provided that the gradient strength
is increased to the extent that the integral of gradient strength with respect
to time remains constant.
This
means that whatever the specifications imposed by an MR imaging experiment on
the design of an ASGC, the design must be as efficient as possible.
ASGC efficiency can be expressed in many ways but one of the most
relevant is in terms of the relationship between the energy stored in the
magnetic field, the gradient strength, the diameter of the system and the ratio
of the diameters of the coil and shield. The
reason that the magnetic stored energy is important is that it is related to the
peak power requirement. This occurs
at the end of the ramp up to flat top. For
a linear ramp with respect to time we can write;
Ppeak
= PL
+ PR
(1)
PL
= 2 E / tr
(2)
where:
Ppeak[VA]
= peak power demand
PL[VA]
= maximum reactive power demand
PR[W]
= maximum resistive power demand
E[J]
= magnetic stored energy
tr[s]
= gradient rise time
It can be shown that the magnetic
stored energy is related to the gradient strength and system diameter by the
relation:
(3)
where:
G[T/m]
= gradient strength
s[m]
= radius of shield
0[H/m]
= permeability of free space
= stored energy factor
The stored energy factor
is a dimensionless parameter that
is a function of the ratio of coil radius to shield radius, of the ratio of coil
length to shield radius and the distribution of turns in the coil.
If we write:
0[m]
= radius of coil
ze[m] = half length of coil
= normalized radius of coil
= normalized half length of coil
= normalized turns distribution
then:
It can further be shown that:
(4)
where the approximation improves as
the relative system length
increases.
Equation
3 shows that for geometrically similar systems the magnetic stored energy of the
system is proportional to the square of the gradient strength, e.g. four times
the stored energy for twice the gradient strength and to the fifth power of the
shield radius, e.g. twice the stored energy for a 15 percent increase in the
shield radius. Equation 4 shows
that the stored energy factor
becomes a very sensitive function
of the radius ratio as the normalized coil radius
tends towards unity.
For example, in going from
= 0.8 to
= 0.9 the stored energy factor
increases by a factor of about 2.4.
The
confined space that is usually available for pulsed gradient coil occupation
means that relatively high radius ratios
are the norm.
This, combined with the demand for higher gradient strengths, G, and
shorter pulse rise times, tr, creates a high reactive power demand, PL, as shown
in equation 2. It is usually so high that it dominates the relation between
reactive power and resistive power, PR, in equation 1.
Thus the peak power demand from the gradient amplifiers is essentially
determined by the magnetic stored energy of the system and there is a high
premium on designing an ASGC with minimum stored energy consistent with the
constraints imposed by gradient linearity and duty cycle, the latter being a
function of the resistive power dissipation and cooling efficiency.
COIL
DESIGN OPTIMIZATION
Once
the pulsed gradient coil envelope dimensions have been defined, i.e. the maximum
and minimum OD and ID respectively, three coils and three shields corresponding
to the three pulsed gradient coil axes are then determined according to
tentative conductor thickness dimensions. The
conductor thickness dimensions are supposed to make resistive power dissipation
consistent with duty cycle specifications.
The coils and shields are then maximally separated to minimize the coil
to shield radius ratio
and a tentative coil length
supposed to be consistent with linearity specifications is selected.
Given
a coil length and radius and a coil turns distribution it is possible to obtain
the magnetic scalar potential corresponding to this distribution everywhere
inside an infinitely long cylinder of radius equal to the shield radius subject,
to the boundary condition that the radial field on the cylinder radius is zero,
i.e., there is no penetration of magnetic flux. This involves solving Laplace's equation in cylindrical
coordinates.
The
resulting analytical solution is in terms of Bessel functions.
This solution corresponds to the state of ideal shielding.
Having obtained that magnetic scalar potential solution it is possible to
derive the shield turns distribution from the magnetic scalar potential
distribution at the shield radius. Analytic
expressions involving Bessel functions are also obtained for the source
coefficients corresponding to a spherical harmonic expansion of the magnetic
scalar potential in the central region. This
is used in the process of meeting field gradient linearity specifications.
Although
this distribution is consistent with the magnetic field shielding constraint it
is not in general consistent with the minimum stored energy requirement nor with
the linearity constraints. To
achieve minimum stored energy consistent with constrained gradient linearity
requires the ability to manipulate the distribution of turns within the coil.
To this end the coil turns distribution is modeled as a Fourier series of
appropriate symmetry with unknown amplitudes.
The
solution depends upon finding a set of amplitudes that are consistent with
minimizing stored energy and satisfying the gradient linearity constraint.
The problem is highly non-linear and objective function gradient
evaluations are impossible analytically.
This
precludes the use of many efficient function minimization methods.
The algorithm used to find the minimum energy then is a variant, due to
Brent1,
of the conjugate direction method invented by Powell2.
This method does not require any information concerning the derivatives
with respect to the optimization variables, of the objective function, in this
case the system stored energy; only function evaluations are required.
This alone, only allows a minimum energy solution to be found; to impose
constraints on gradient linearity requires additional treatment. This is achieved by using a modified objective function.
The
objective function i.e., the magnetic stored energy, is modified by a penalty
term that is a function of the extent by which a particular configuration
violates the linearity constraint. The
particular method used is Bertsekis's3
augmented Lagrangian multiplier modification of the exterior penalty method.
The solution then proceeds by a series of unconstrained minimizations of
a modified objective function4.
After
the completion of such a series a solution is obtained that is both of minimum
stored energy and which satisfies the linearity constraints. Constraints other than gradient linearity are also applied
during the solution process. These
include constraints on such parameters as maximum current density and resistive
power dissipation. Having
obtained such a solution the design is further investigated to determine its
thermal and mechanical characteristics. Adjustments
to coil length, conductor thickness, etc., are made and the process is repeated
for as many times as is necessary to obtain an optimal design meeting all design
specifications.
REFERENCES
-
Algorithms for Minimization without
Derivatives, Richard P. Brent, Prentice-Hall Inc., Englewood Cliffs, NJ, 1973.
-
Powell,
M.J.D., An Efficient Method for Finding the Minimum of a Function
of Several Variables without Calculating Derivatives, Comp. J., 7, 155-162,
1965.
-
Constrained Optimization and Lagrange Multiplier Methods, Dimitri P.
Bertsekis, Academic Press, 1982.
-
Nonlinear Programming, Theory, Algorithms, and Applications, Garth P.
McCormick, John Wiley and Sons, New York, 1982.
|