A Divergence-Free BEM Method to Model Quasi-Static Currents: Application to MRI Coil Design

The modeling of quasi-static optimization problems often involves divergence-free surface current densities. In this paper, a novel technique to implement these currents by using the boundary element method framework is presented. A locally-based characterization of the current density is employed, to render a fully geometry-independent formulation, so that it can be applied to arbitrary shapes. To illustrate the versatility of this approach, we employ it for the design of gradient coils for MRI, providing a solid mathematical framework for this type of problem. 1. INTRODUCTION Many problems in engineering require to determine the spatial distribution of electric currents ∞owing in a conductive surface which satisfles given requirements for the flelds, electromagnetic energy, etc. they produce. The reconstruction of current distribution on the conducting surface subjected to these constraints can be seen as an inverse problem. An appropriate and realistic formulation of this type of problems is presented in this paper, by incorporating a suitable model of the current under search, in terms of the stream function, into the Boundary Element Method (BEM). The use of stream function for the characterization of surface current densities has been widely employed (1). The BEM has been proved to be an excellent tool in the solution of electromagnetic problems (2,3); and the incorporation of the stream function into a numerical computational technique, such us BEM, has also been already considered for the solution of electromagnetic inverse

[1]  D. A. Dunnett Classical Electrodynamics , 2020, Nature.

[2]  Jeng-Tzong Chen,et al.  A Practical Guide to Boundary Element Methods with the Software Library BEMLIB , 2002 .

[3]  R. Bowtell,et al.  Boundary Element Method for an Inverse Problem in Magnetic Resonance Imaging Gradient Coils , 2008 .

[4]  Adib A. Becker,et al.  Numerical solution for an inverse MRI problem using a regularised boundary element method. , 2008 .

[5]  Thomas F. Eibert,et al.  On the calculation of potential integrals for linear source distributions on triangular domains , 1995 .

[6]  R. Holland Finite-Difference Solution of Maxwell's Equations in Generalized Nonorthogonal Coordinates , 1983, IEEE Transactions on Nuclear Science.

[7]  C. Brebbia,et al.  Boundary Element Techniques , 1984 .

[8]  F. París,et al.  Boundary Element Method: Fundamentals and Applications , 1997 .

[9]  S. Pissanetzky Minimum energy MRI gradient coils of general geometry , 1992 .

[10]  Reinhold Ludwig,et al.  A Stream Function Method for Gradient Coil Design , 2005 .

[11]  Use of stream functions for the computation of currents in thin circuits determination of the impedances , 2000 .

[12]  Mattan Kamon,et al.  FASTHENRY: a multipole-accelerated 3-D inductance extraction program , 1994 .

[13]  Dexin Xie,et al.  A Novel Target Field Method for Designing Uniplanar Self-shield Gradient Coils of Fully Open MRI Device , 2007 .

[14]  A New Target Field Method for Optimizing Longitudinal Gradient Coils' Property , 2007 .

[15]  G. Peeren Stream function approach for determining optimal surface currents , 2003 .

[16]  Donglin Zu,et al.  Shimming Permanent Magnet of MRI Scanner , 2007 .

[17]  Simon R. Arridge,et al.  3D Shape Reconstruction in Optical Tomography Using Spherical Harmonics and BEM , 2006 .

[18]  Philip E. Gill,et al.  Practical optimization , 1981 .

[19]  R. Turner,et al.  Gradient coil design: a review of methods. , 1993, Magnetic resonance imaging.

[20]  Manuchehr Soleimani,et al.  A trust region subproblem for 3D electrical impedance tomography inverse problem using experimental data , 2009 .

[21]  Stuart Crozier,et al.  Determining complicated winding patterns for shim coils using stream functions and the target-field method , 2002 .