Abstract
This paper presents a fast-multipole surface-charge-simulation method for calculating three-dimensional Laplacian fields in voxel models. This method treats a surface of a voxel that has different inside and outside conductivities as a surface element of the indirect BEM(Boundary Element Method). The main features of the proposed method are as follows. (1) An O(D2) performance in the memory-capacity and operation-cost is provided by applying the diagonal form Fast-Multipole-Method, when the number of voxels is about D3. (2) The boundary matching is imposed by the continuity of the total flux passing through each element, which guarantees the solution globally satisfying the Gauss' law, therefore the solution is globally stabilized. This method is successfully applied to calculate the electric field induced by an applied homogeneous ELF(Extremely Low Frequency) magnetic field in a human head model that has 1mm×1mm×1mm voxel size.
