Numerical solution for an inverse MRI problem using a regularised boundary element method

We investigate the reconstruction of a divergence-free surface current distribution from knowledge of the magnetic flux density in a prescribed region of interest in the framework of static electromagnetism. This inverse problem is motivated by the design of gradient coils used in magnetic resonance...

Full description

Saved in:
Bibliographic Details
Published in:Engineering analysis with boundary elements Vol. 32; no. 8; pp. 658 - 675
Main Authors: Marin, Liviu, Power, Henry, Bowtell, Richard W., Cobos Sanchez, Clemente, Becker, Adib A., Glover, Paul, Jones, Arthur
Format: Journal Article
Language:English
Published: Elsevier Ltd 01.08.2008
Subjects:
Boundary element method
Coil design
Coils
Current density
Density
Design engineering
Divergence-free boundary element method (BEM)
Electromagnetism
Inverse problem
Magnetic resonance imaging (MRI)
Mathematical analysis
Mathematical models
Regularisation
Streams
Magnetic resonance imaging (MRI)
Electromagnetism
Coil design
Regularisation
Divergence-free boundary element method (BEM)
Inverse problem
ISSN:0955-7997, 1873-197X
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We investigate the reconstruction of a divergence-free surface current distribution from knowledge of the magnetic flux density in a prescribed region of interest in the framework of static electromagnetism. This inverse problem is motivated by the design of gradient coils used in magnetic resonance imaging (MRI) and is formulated using its corresponding integral representation according to potential theory. A constant boundary element method (BEM) which satisfies the continuity equation for the current density, i.e. divergence-free BEM, and was originally proposed by Lemdiasov and Ludwig [A stream function method for gradient coil design. Concepts Magn Reson B Magn Reson Eng 2005;26B:67–80], is presented based on geometrical arguments with respect to the linear (flat) triangular boundary elements employed in order to emphasise its possible extension to further higher-order divergence-free interpolations. Since the discretised BEM system is ill-posed and hence the associated least-squares solution may be inaccurate and/or physically meaningless, the Tikhonov regularisation method is employed in order to retrieve accurate and physically correct solutions. A rigorous numerical approximation for the calculation the magnetic energy, which reduces the errors induced by employing the approach of Lemdiasov and Ludwig [A stream function method for gradient coil design. Concepts Magn Reson B Magn Reson Eng 2005;26B:67–80], is also proposed.
Bibliography:ObjectType-Article-2
SourceType-Scholarly Journals-1
ObjectType-Feature-1
content type line 23
ISSN:0955-7997
1873-197X
DOI:10.1016/j.enganabound.2007.11.015