Incomplete Cholesky Conjugate Gradient Method for the Three- dimensional Forward Problem in Magnetic Induction Tomography Using Finite Element Method
XUAN Yang① WANG Xu①② LIU Cheng’an① YANG Dan② ZHANG Zhimei②
①(Sino-Dutch Biomedical and Information Engineering School, Northeastern University, Shenyang 110004, China) ②(College of Information Science and Engineering, Northeastern University, Shenyang 110004, China)
In 3D forward problem of Magnetic Induction Tomography (MIT), the problems are slow computation speeds and incorrect results due to round-off errors, when calculating the finite element equations with the direct method. Incomplete Cholesky Conjugate Gradient (ICCG) iteration method is used to solve these problems. Round-off errors are compensated by iteration method. An Finite-Element Model (FEM) is built based on the ANSYS software. The FEM equations are solved by the ICCG method. The optimal convergence tolerance value is calculated. Simulation result shows that the ICCG method has advantages in speed and stability compared with direct and Jacobi Conjugate Gradient (JCG) method. The results show that the ICCG method is not affected by meshing perturbation, it can solve the 3D forward problem of MIT correctly.
宣杨,王旭,刘承安,杨丹,张志美. 不完全乔列斯基分解共轭梯度法在磁感应成像三维有限元正问题中的应用[J]. 电子与信息学报, 2016, 38(1): 187-194.
XUAN Yang, WANG Xu, LIU Cheng’an, YANG Dan, ZHANG Zhimei. Incomplete Cholesky Conjugate Gradient Method for the Three- dimensional Forward Problem in Magnetic Induction Tomography Using Finite Element Method. JEIT, 2016, 38(1): 187-194.
WEI H Y and SOLEIMANI M. Electromagnetic tomography for medical and industrial applications: challenges and opportunities[J]. Proceedings of the IEEE, 2013, 101(3): 559-565.
[2]
LU Ma, HUNT Andy, and SOLEIMANI M. Experimental evaluation of conductive flow imaging using magnetic induction tomography[J]. International Journal of Multiphase Flow, 2015, 72: 198-209.
[3]
JIN Gui, SUN Jian, QIN Mingxin, et al. A new method for detecting cerebral hemorrhage in rabbits by magnetic inductive phase shift[J]. Biosensors and Bioelectronics, 2014, 52: 374-378.
WANG L, LIU R G, ZHOU W, et al. Research progress of magnetic induction tomography system and experimental results[J]. Chinese Medical Equipment Journal, 2013, 34(2): 85-88.
[5]
DARRER B J, WATSON J C, BARTLETT P, et al. Toward an automated setup for magnetic induction tomography[J]. IEEE Transactions on Magnetics, 2015, 51(1): 6500104.
[6]
MERWA R, HOLLAUS K, BRANDSTATTER B, et al. Volumetric magnetic induction tomography[J]. Measurement Science and Technology, 2012, 23(5): 055401.
[7]
MERWA R, HOLLAUS K, BRANDSTATTER B, et al. Numerical solution of the general 3D eddy current problem for magnetic induction tomography (spectroscopy)[J]. Physiological Measurement, 2003, 24: 545-554.
[8]
HOLLAUS K, MAC, MERWA R, et al. Numerical simulation of the eddy current problem in magnetic induction tomography for biomedical applications by edge elements[J]. IEEE Transactions on Magnetics, 2004, 40(2): 623-626.
LIU G Q, WANG T, MENG M, et al. Using edge element method to solve the forward problem in magnetic induction tomography[J]. China Medical Engineering, 2006, 25(2): 163-165.
HE W, SONG X D, ZHANG X Y, et al. The edge element method in the forward problem of magnetic induction tomography with homogeneous sinusoidal magnetic excitation[J]. Journal of Chongqing University, 2012, 35(10): 159-164.
KE L, ZHAO L L, DU Q, et al. Simulation and analysis of cerebral hematoma eddy current field in magnetic induction tomography[J]. Journal of System Simulation, 2014, 26(3): 517-522.
[12]
ZHAO Q, CHEN G, HAO J N, et al. Numerical approach for the sensitivity of a high-frequency magnetic induction tomography system based on boundary elements and perturbation method[J]. Measurement Science and Technology, 2013, 4: 074004.
MENG M, JIANG L T, LI S Q, et al. 3-D reconstruction algorithm of magnetic resonance magnetic induction tomography[J]. Chinese Journal of Biomedical Engineering, 2008, 27(5): 650-653.
[14]
SOLEIMANI M, WILLIAM R, and LIONHEART B. Absolute conductivity reconstruction in magnetic induction tomography using a nonlinear method[J]. IEEE Transactions on Medical Imaging, 2006, 25(12): 1521-1530.
[15]
SOLEIMANI1 M, LIONHEART W R B, PEYTON1 A J, et al. A three-dimensional inverse finite-element method applied to experimental eddy-current imaging Data[J]. IEEE Transactions on Magnetics, 2006, 42(5): 1560-1567.
LIN X and LIU Z G. The application of ICCG algorithm in the three dimensional electric field calculation of SF6 tank-type circuit breaker[J]. Proceedings of the CSEE, 2001, 21(2): 21-24.
XIE D X and YANG S Y. Numerical Analysis and Synthesis of Engineering Electromagnetic Field[M]. Beijing: China Machine Press, 2008: 84-85.
[18]
KAMEARI A. Symetric second order edge elements for triangles and terahedra[J]. IEEE Transactions on Magnetics, 1999, 35: 1394-1397.
[19]
GRIFFITHS H, STEWART W R, and GOUGH W. Magnetic induction tomography a measuring system for biological tissues[J]. Annals New York Academy of Sciences, 1999, 20: 335-345.
[20]
HERMANN S, HELMUT K L, and JAVIER R. Magnetic induction tomography: hardware for multi-frequency measurement in biological tissues[J]. Physiological Measurement, 2001, 22: 131-146.