Mixed Finite Element Method for Full-Wave Simulation of Bioelectromagnetism From DC to Microwave Frequencies

Bioelectromagnetism focuses on the study of electromagnetic fields in biological tissues from direct current (DC) to optical frequencies. It is challenging to develop an electromagnetics (EM) simulation method to cover this entire frequency band due to the electrically small/large scattering problem at extremely low/high frequencies. This paper focuses on the band from DC to microwave frequencies in bioelectromagnetism. Its main research objective is to develop a method that can overcome the low frequency breakdown problem at low frequencies (practically DC) and still stay stable at microwave frequencies. Based on the scattered field vector Helmholtz equation, the mixed finite element method (mixed FEM) is developed for the broadband electromagnetic field simulation in biological tissues. By imposing Gauss’ law as the constraint condition, the mixed FEM overcomes the low frequency breakdown problem without resorting to the quasi-static approximation and remains effective and accurate at high frequencies. Extremely low frequency and high frequency numerical results are demonstrated to verify that the mixed FEM is a stable full-wave electromagnetic field simulation method for the full-bandwidth bioelectromagnetism.

[1]  Qing Huo Liu,et al.  3-D MRI-Based Electrical Properties Tomography Using the Volume Integral Equation Method , 2017, IEEE Transactions on Microwave Theory and Techniques.

[2]  P. Kosmas,et al.  Modeling with the FDTD method for microwave breast cancer detection , 2004, IEEE Transactions on Microwave Theory and Techniques.

[3]  H. Gach,et al.  2D & 3D Shepp-Logan Phantom Standards for MRI , 2008, 2008 19th International Conference on Systems Engineering.

[4]  Won Hee Lee,et al.  Electric Field Model of Transcranial Electric Stimulation in Nonhuman Primates: Correspondence to Individual Motor Threshold , 2015, IEEE Transactions on Biomedical Engineering.

[5]  David Isaacson,et al.  Electrical Impedance Tomography , 1999, SIAM Rev..

[6]  Qing Huo Liu,et al.  An Efficient Mixed Finite-Element Time-Domain Method for Complex Electrically Small Problems , 2019, IEEE Transactions on Microwave Theory and Techniques.

[7]  D. Tucker,et al.  EEG source localization: Sensor density and head surface coverage , 2015, Journal of Neuroscience Methods.

[8]  Sara Louie,et al.  Virtual Human Models for Electromagnetic Studies and Their Applications , 2017, IEEE Reviews in Biomedical Engineering.

[9]  Min Han,et al.  Magnetic Induction Tomography , 2015 .

[10]  Feng Liu,et al.  An Efficient Integral-Based Method for Three-Dimensional MR-EPT and the Calculation of the RF-Coil-Induced ${B_z}$ Field , 2018, IEEE Transactions on Biomedical Engineering.

[11]  Xavier Tricoche,et al.  Influence of tissue conductivity anisotropy on EEG/MEG field and return current computation in a realistic head model: A simulation and visualization study using high-resolution finite element modeling , 2006, NeuroImage.

[12]  Qing Huo Liu,et al.  Mixed Spectral-Element Method for Overcoming the Low-Frequency Breakdown Problem in Subsurface EM Exploration , 2017, IEEE Transactions on Geoscience and Remote Sensing.

[13]  M. Hallett Transcranial magnetic stimulation and the human brain , 2000, Nature.

[14]  R. Ilmoniemi,et al.  Magnetoencephalography-theory, instrumentation, and applications to noninvasive studies of the working human brain , 1993 .

[15]  M. Murray,et al.  EEG source imaging , 2004, Clinical Neurophysiology.

[16]  Z. Nie,et al.  VECTOR FINITE ELEMENT ANALYSIS OF MULTICOMPONENT INDUCTION RESPONSE IN ANISOTROPIC FORMATIONS , 2008 .

[17]  Eric Michielssen,et al.  The ICVSIE: A General Purpose Integral Equation Method for Bio-Electromagnetic Analysis , 2018, IEEE Transactions on Biomedical Engineering.

[18]  X. Li,et al.  Confocal microwave imaging for breast cancer detection: localization of tumors in three dimensions , 2002, IEEE Transactions on Biomedical Engineering.

[19]  J. Fermaglich Electric Fields of the Brain: The Neurophysics of EEG , 1982 .

[20]  N. Höfner,et al.  Demonstration of full tensor current density imaging using ultra-low field MRI. , 2019, Magnetic resonance imaging.

[21]  Qingtao Sun,et al.  Mixed Total Field/Scattered Field-Based Discontinuous Galerkin Frequency-Domain Method for Subsurface Sensing , 2019, IEEE Transactions on Geoscience and Remote Sensing.

[22]  C. Gabriel Compilation of the Dielectric Properties of Body Tissues at RF and Microwave Frequencies. , 1996 .