Using domain decomposition techniques for the calculation of low-frequency electric current densities in high-resolution 3D human anatomy models

Purpose – Improved numerical calculation techniques for low‐frequency current density distributions within high‐resolution anatomy models caused by ambient electric or magnetic fields or direct contact to potential drops using the finite integration technique (FIT).Design/methodology/approach – The methodology of calculating low‐frequency electromagnetic fields within high‐resolution anatomy models using the FIT is extended by a local grid refinement scheme using a non‐matching‐grid formulation domain. Furthermore, distributed computing techniques are presented. Several numerical examples are analyzed using these techniques.Findings – Numerical simulations of low‐frequency current density distributions may now be performed with a higher accuracy due to an increased local grid resolution in the areas of interest in the human body voxel models when using the presented techniques.Originality/value – The local subgridding approach is introduced to reduce the number of unknowns in the very large‐scale linear a...

[1]  Barbara I. Wohlmuth,et al.  Discretization Methods and Iterative Solvers Based on Domain Decomposition , 2001, Lecture Notes in Computational Science and Engineering.

[2]  C D Werner,et al.  Development of a human body model for numerical calculation of electrical fields. , 2000, Computerized medical imaging and graphics : the official journal of the Computerized Medical Imaging Society.

[3]  Thomas Weiland,et al.  Numerical simulation of low-frequency current density distributions in voxel-based human anatomy models due to ambient electric and magnetic fields , 2004 .

[4]  Riccardo Scorretti,et al.  Electromagnetic fields and human body: a new challenge for the electromagnetic field computation , 2003 .

[5]  V. E. Henson,et al.  BoomerAMG: a parallel algebraic multigrid solver and preconditioner , 2002 .

[6]  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.

[7]  Thomas Weiland,et al.  Simulation of low-frequency fields on high-voltage insulators with light contaminations , 1996 .

[8]  T. Weiland Time Domain Electromagnetic Field Computation with Finite Difference Methods , 1996 .

[9]  Edward K. N. Yung,et al.  TRANSIENT ANALYSIS OFn-TIER GaAs MESFET MATRIX AMPLIFIERS BY THE TLM METHOD , 1996 .

[10]  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.

[11]  M. Stuchly,et al.  Modelling fields induced in humans by 50/60 Hz magnetic fields: reliability of the results and effects of model variations. , 2002, Physics in medicine and biology.

[12]  Thomas Weiland,et al.  3D eddy current computation in the frequency domain regarding the displacement current , 1992 .