Multi-GPU Accelerated Admittance Method for High-Resolution Human Exposure Evaluation

Objective: A multi-graphics processing unit (GPU) accelerated admittance method solver is presented for solving the induced electric field in high-resolution anatomical models of human body when exposed to external low-frequency magnetic fields. Methods: In the solver, the anatomical model is discretized as a three-dimensional network of admittances. The conjugate orthogonal conjugate gradient (COCG) iterative algorithm is employed to take advantage of the symmetric property of the complex-valued linear system of equations. Compared against the widely used biconjugate gradient stabilized method, the COCG algorithm can reduce the solving time by 3.5 times and reduce the storage requirement by about 40%. The iterative algorithm is then accelerated further by using multiple NVIDIA GPUs. The computations and data transfers between GPUs are overlapped in time by using asynchronous concurrent execution design. The communication overhead is well hidden so that the acceleration is nearly linear with the number of GPU cards. Results: Numerical examples show that our GPU implementation running on four NVIDIA Tesla K20c cards can reach 90 times faster than the CPU implementation running on eight CPU cores (two Intel Xeon E5-2603 processors). Conclusion: The implemented solver is able to solve large dimensional problems efficiently. A whole adult body discretized in 1-mm resolution can be solved in just several minutes. Significance: The high efficiency achieved makes it practical to investigate human exposure involving a large number of cases with a high resolution that meets the requirements of international dosimetry guidelines.

[1]  Akimasa Hirata,et al.  Corrigendum: Computational dosimetry for grounded and ungrounded human models due to contact current (2013 Phys. Med. Biol. 58 5153?72) , 2014 .

[2]  V. Spitzer,et al.  The visible human male: a technical report. , 1996, Journal of the American Medical Informatics Association : JAMIA.

[3]  J G Jefferys,et al.  Experimental neurobiology of epilepsies. , 1994, Current opinion in neurology.

[4]  Roland W. Freund,et al.  An Implementation of the Look-Ahead Lanczos Algorithm for Non-Hermitian Matrices , 1993, SIAM J. Sci. Comput..

[5]  R. W. Lau,et al.  The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues. , 1996, Physics in medicine and biology.

[6]  Sami Ilvonen,et al.  Computational estimation of magnetically induced electric fields in a rotating head , 2009, Physics in medicine and biology.

[7]  Jianping Cheng,et al.  CNMAN: a Chinese adult male voxel phantom constructed from color photographs of a visible anatomical data set. , 2007, Radiation protection dosimetry.

[8]  M. Douglas,et al.  Issues of ICNIRP guidelines when determining compliance with LF exposure limits , 2012, International Symposium on Electromagnetic Compatibility - EMC EUROPE.

[9]  P. Dimbylow,et al.  Development of pregnant female, hybrid voxel-mathematical models and their application to the dosimetry of applied magnetic and electric fields at 50 Hz , 2006, Physics in medicine and biology.

[10]  Hiroshi Kanai,et al.  Impedence Method for Calculation of Power Deposition Patterns in Magnetically Induced Hyperthermia , 1984, IEEE Transactions on Biomedical Engineering.

[11]  Theodoros Samaras,et al.  The effect of head and coil modeling for the calculation of induced electric field during transcranial magnetic stimulation. , 2014, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[12]  Luc Dupré,et al.  An Efficient 3-D Eddy-Current Solver Using an Independent Impedance Method for Transcranial Magnetic Stimulation , 2011, IEEE Transactions on Biomedical Engineering.

[13]  Maria A. Stuchly,et al.  Comparison of magnetically induced elf fields in humans computed by FDTD and scalar potential FD codes , 1996, 1996 Symposium on Antenna Technology and Applied Electromagnetics.

[14]  P. Dimbylow,et al.  Induced current densities from low-frequency magnetic fields in a 2 mm resolution, anatomically realistic model of the body. , 1998, Physics in medicine and biology.

[15]  Teruo Onishi,et al.  Analysis of in situ electric field and specific absorption rate in human models for wireless power transfer system with induction coupling. , 2014, Physics in medicine and biology.

[16]  Akimasa Hirata,et al.  Computational dosimetry for grounded and ungrounded human models due to contact current. , 2013, Physics in medicine and biology.

[17]  A. Hirata,et al.  Improving the computational speed and reducing the staircasing error for simulations of human exposure to low frequency magnetic fields , 2012, International Symposium on Electromagnetic Compatibility - EMC EUROPE.

[18]  A. Hirata,et al.  Reducing the staircasing error in computational dosimetry of low-frequency electromagnetic fields , 2012, Physics in medicine and biology.

[19]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[20]  Jack J. Purdum,et al.  C programming guide , 1983 .

[21]  K. Jokela,et al.  ICNIRP Guidelines GUIDELINES FOR LIMITING EXPOSURE TO TIME-VARYING , 1998 .

[22]  N. Kuster,et al.  Analysis of human brain exposure to low‐frequency magnetic fields: A numerical assessment of spatially averaged electric fields and exposure limits , 2013, Bioelectromagnetics.

[23]  Akimasa Hirata,et al.  Corrigendum: Reducing the staircasing error in computational dosimetry of low-frequency electromagnetic fields , 2012 .

[24]  Michal Okoniewski,et al.  Low Frequency Finite Difference Time Domain (FDTD) for Modeling of Induced Fields in Humans Close to Line Sources , 2000 .

[25]  T. Weiland,et al.  Impact of the displacement current on low-frequency electromagnetic fields computed using high-resolution anatomy models , 2005, Physics in medicine and biology.

[26]  Ling Xia,et al.  An efficient impedance method for induced field evaluation based on a stabilized Bi-conjugate gradient algorithm , 2008, Physics in medicine and biology.

[27]  G. Schmid,et al.  The role of skin conductivity in a low frequency exposure assessment for peripheral nerve tissue according to the ICNIRP 2010 guidelines , 2013, Physics in medicine and biology.

[28]  H. V. D. Vorst,et al.  A Petrov-Galerkin type method for solving Axk=b, where A is symmetric complex , 1990 .

[29]  Jin-Kyu Byun,et al.  Induced Current Calculation in Detailed 3-D Adult and Child Model for the Wireless Power Transfer Frequency Range , 2014, IEEE Transactions on Magnetics.

[30]  H. De Gersem,et al.  Large-Scale Calculations of Low-Frequency Induced Currents in High-Resolution Human Body Models. , 2006, 2006 12th Biennial IEEE Conference on Electromagnetic Field Computation.

[31]  R. W. Lau,et al.  The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz. , 1996, Physics in medicine and biology.

[32]  Gene H. Golub,et al.  Matrix computations , 1983 .

[33]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[34]  Jan De Moerloose,et al.  Application of the finite difference time domain algorithm to quasi‐static field analysis , 1997 .

[35]  Maria A. Stuchly,et al.  Pacemaker interference and low-frequency electric induction in humans by external fields and electrodes , 2000, IEEE Transactions on Biomedical Engineering.

[36]  Niels Kuster,et al.  The Virtual Family—development of surface-based anatomical models of two adults and two children for dosimetric simulations , 2010, Physics in medicine and biology.

[37]  E. Neufeld,et al.  Children and adults exposed to low-frequency magnetic fields at the ICNIRP reference levels: theoretical assessment of the induced electric fields , 2012, Physics in medicine and biology.

[38]  Miriam Leeser,et al.  OpenCL Floating Point Software on Heterogeneous Architectures – Portable or Not? , 2012 .

[39]  C Gabriel,et al.  The dielectric properties of biological tissues: I. Literature survey. , 1996, Physics in medicine and biology.

[40]  Maria A. Stuchly,et al.  Validation of a scalar potential finite-difference code for extremely low frequency magnetic induction modelling using an analytical method , 1996, 1996 Symposium on Antenna Technology and Applied Electromagnetics.

[41]  Akimasa Hirata,et al.  Fast multigrid-based computation of the induced electric field for transcranial magnetic stimulation , 2012, Physics in medicine and biology.

[42]  H. De Gersem,et al.  Large-Scale Calculation of Low-Frequency-Induced Currents in High-Resolution Human Body Models , 2007, IEEE Transactions on Magnetics.

[43]  Joshua R. Smith,et al.  Evaluation of Wireless Resonant Power Transfer Systems With Human Electromagnetic Exposure Limits , 2013 .

[44]  Sang Wook Park,et al.  Calculation Errors of the Electric Field Induced in a Human Body Under Quasi-Static Approximation Conditions , 2013, IEEE Transactions on Microwave Theory and Techniques.

[45]  O. Gandhi,et al.  A 3-D impedance method to calculate power deposition in biological bodies subjected to time varying magnetic fields , 1988, IEEE Transactions on Biomedical Engineering.

[46]  E. Neufeld,et al.  IT’IS Database for Thermal and Electromagnetic Parameters of Biological Tissues , 2012 .

[47]  K G Efthimiadis,et al.  Magnetic stimulation of the spine: the role of tissues and their modelling , 2010, Physics in medicine and biology.