Design multiple-layer gradient coils using least-squares finite element method

The design of gradient coils for magnetic resonance imaging is an optimization task in which a specified distribution of the magnetic field inside a region of interest is generated by choosing an optimal distribution of a current density geometrically restricted to specified non-intersecting design surfaces, thereby defining the preferred coil conductor shapes. Instead of boundary integral type methods, which are widely used to design coils, this paper proposes an optimization method for designing multiple layer gradient coils based on a finite element discretization. The topology of the gradient coil is expressed by a scalar stream function. The distribution of the magnetic field inside the computational domain is calculated using the least-squares finite element method. The first-order sensitivity of the objective function is calculated using an adjoint equation method. The numerical operations needed, in order to obtain an effective optimization procedure, are discussed in detail. In order to illustrate the benefit of the proposed optimization method, example gradient coils located on multiple surfaces are computed and characterised.