Custom edge‐element FEM solver and its application to eddy‐current simulation of realistic 2M‐element human brain phantom

Extensive research papers of three-dimensional computational techniques are widely used for the investigation of human brain pathophysiology. Eddy current analyzing could provide an indication of conductivity change within a biological body. A significant obstacle to current trend analyses is the development of a numerically stable and efficiency-finite element scheme that performs well at low frequency and does not require a large number of degrees of freedom. Here, a custom finite element method (FEM) solver based on edge elements is proposed using the weakly coupled theory, which separates the solution into two steps. First, the background field (the magnetic vector potential on each edge) is calculated and stored. Then, the electric scalar potential on each node is obtained by FEM based on Galerkin formulations. Consequently, the electric field and eddy current distribution in the object can be obtained. This solver is more efficient than typical commercial solvers since it reduces the vector eddy current equation to a scalar one, and reduces the meshing domain to just the eddy current region. It can therefore tackle complex eddy current calculations for models with much larger numbers of elements, such as those encountered in eddy current computation in biological tissues. An example is presented with a realistic human brain mesh of 2 million elements. In addition, with this solver, the equivalent magnetic field induced from the excitation coil is applied, and therefore there is no need to mesh the excitation coil. In combination, these significantly increase the efficiency of the solver. Bioelectromagnetics. 39:604-616, 2018. © 2018 Wiley Periodicals, Inc.

[1]  S. Alekseev,et al.  Influence of blood flow and millimeter wave exposure on skin temperature in different thermal models , 2009, Bioelectromagnetics.

[2]  Charbel Farhat,et al.  An Unconventional Domain Decomposition Method for an Efficient Parallel Solution of Large-Scale Finite Element Systems , 1992, SIAM J. Sci. Comput..

[3]  DAWSON TURNER.  The Electrical Resistance of the Blood , 1902, Nature.

[4]  R H Bayford,et al.  Electrical impedance tomography spectroscopy (EITS) for human head imaging. , 2003, Physiological measurement.

[5]  Wuliang Yin,et al.  A Novel Compensation Algorithm for Thickness Measurement Immune to Lift-Off Variations Using Eddy Current Method , 2016, IEEE Transactions on Instrumentation and Measurement.

[6]  A. Kameari,et al.  Calculation of three dimensional eddy current by FEM-BEM coupling method , 1988 .

[7]  John C. Moulder,et al.  Review of Advances in Quantitative Eddy Current Nondestructive Evaluation , 1999 .

[8]  J. S. Cohen,et al.  Differences in metabolite levels upon differentiation of intact neuroblastoma × glioma cells observed by proton NMR spectroscopy , 1983, FEBS letters.

[9]  T. Tenforde,et al.  Theoretical analysis of magnetic field interactions with aortic blood flow. , 1996, Bioelectromagnetics.

[10]  H Buchner,et al.  Current density threshold for the stimulation of neurons in the motor cortex area , 2002, Bioelectromagnetics.

[11]  Joachim Georgii,et al.  A STREAMING APPROACH FOR SPARSE MATRIX PRODUCTS AND ITS APPLICATION IN GALERKIN MULTIGRID METHODS , 2010 .

[12]  B H Brown,et al.  The Sheffield data collection system. , 1987, Clinical physics and physiological measurement : an official journal of the Hospital Physicists' Association, Deutsche Gesellschaft fur Medizinische Physik and the European Federation of Organisations for Medical Physics.

[13]  M. Stuchly,et al.  Influence of human model resolution on computed currents induced in organs by 60-Hz magnetic fields. , 1997, Bioelectromagnetics.

[14]  Jian-Ming Jin,et al.  Finite-element domain decomposition methods for analysis of large-scale electromagnetic problems , 2014 .

[15]  Ulrich Jakobus,et al.  Human exposure assessment in the near field of GSM base-station antennas using a hybrid finite element/method of moments technique , 2003, IEEE Transactions on Biomedical Engineering.

[16]  S. Abdalla,et al.  Low Frequency Dielectric Properties of Human Blood , 2011, IEEE Transactions on NanoBioscience.

[17]  Jian-Ming Jin,et al.  A Nonconformal FEM-DDM With Tree–Cotree Splitting and Improved Transmission Condition for Modeling Subsurface Detection Problems , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[18]  Zhong-Lin Lu,et al.  Magnetic source imaging of the human brain , 1993, Photonics West - Lasers and Applications in Science and Engineering.

[19]  Richard H. Bayford,et al.  Electrical impedance tomography of human brain function using reconstruction algorithms based on the finite element method , 2003, NeuroImage.

[20]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[21]  R H Bayford,et al.  Multi-frequency electrical impedance tomography (EIT) of the adult human head: initial findings in brain tumours, arteriovenous malformations and chronic stroke, development of an analysis method and calibration , 2006, Physiological measurement.

[22]  T. Eibert,et al.  Electromagnetic and thermal analysis for lipid bilayer membranes exposed to RF fields , 1999, IEEE Transactions on Biomedical Engineering.

[23]  Richard H. Bayford,et al.  Three-Dimensional Electrical Impedance Tomography of Human Brain Activity , 2001, NeuroImage.

[24]  Fumio Kojima,et al.  Identification of crack profiles using genetic programming and fuzzy inference , 2001 .

[25]  Nathan Ida,et al.  Introduction to the Finite Element Method , 1997 .

[26]  R H Bayford,et al.  Using the GRID to improve the computation speed of electrical impedance tomography (EIT) reconstruction algorithms. , 2005, Physiological measurement.

[27]  A Tizzard,et al.  Generating accurate finite element meshes for the forward model of the human head in EIT , 2005, Physiological measurement.

[28]  L. Udpa,et al.  Reduced Magnetic Vector Potential Formulation in the Finite Element Analysis of Eddy Current Nondestructive Testing , 2009, IEEE Transactions on Magnetics.

[29]  D. Kaklamani,et al.  Aspects of the Method of Auxiliary Sources (MAS) in computational electromagnetics , 2002 .

[30]  Wuliang Yin,et al.  A Method to Solve the Forward Problem in Magnetic Induction Tomography Based on the Weakly Coupled Field Approximation , 2010, IEEE Transactions on Biomedical Engineering.

[31]  W. Desmet,et al.  A fast frequency sweep approach using Padé approximations for solving Helmholtz finite element models , 2013 .

[32]  Matti Stenroos,et al.  On the accuracy of collocation and Galerkin BEM in the EEG/MEG forward problem , 2012 .

[33]  Wuliang Yin,et al.  Reducing the Lift-Off Effect on Permeability Measurement for Magnetic Plates From Multifrequency Induction Data , 2018, IEEE Transactions on Instrumentation and Measurement.

[34]  O. Bíró.  Edge element formulations of eddy current problems , 1999 .

[35]  Robert Bridson,et al.  A Structural Diagnosis of Some IC Orderings , 2000, SIAM J. Sci. Comput..

[36]  M. A. Roula,et al.  Evaluation of Parallel Accelerators for High Performance Image Reconstruction for Magnetic Induction Tomography , .

[37]  Jian-Ming Jin,et al.  An Accurate and Efficient Finite Element-Boundary Integral Method With GPU Acceleration for 3-D Electromagnetic Analysis , 2014, IEEE Transactions on Antennas and Propagation.

[38]  Jian-Ming Jin,et al.  A New Dual-Primal Domain Decomposition Approach for Finite Element Simulation of 3-D Large-Scale Electromagnetic Problems , 2007, IEEE Transactions on Antennas and Propagation.

[39]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[40]  M A Stuchly,et al.  Factors affecting neural stimulation with magnetic fields. , 1992, Bioelectromagnetics.

[41]  Kan Xu,et al.  Multilevel fast multipole algorithm enhanced by GPU parallel technique for electromagnetic scattering problems , 2010 .

[42]  Riccardo Hertel,et al.  Guided Spin Waves , 2007 .

[43]  John R. Mortarelli.  A Generalization of the Geselowitz Relationship Useful in Impedance Plethysmographic Field Calculations , 1980, IEEE Transactions on Biomedical Engineering.

[44]  N. Hayashi,et al.  ELF electromagnetic environment in power substations. , 1989, Bioelectromagnetics.

[45]  D. Kell,et al.  Dielectric properties of human blood and erythrocytes at radio frequencies (0.2–10 MHz); dependence on cell volume fraction and medium composition , 2005, European Biophysics Journal.