A Feasibility Study: Operator Splitting for Solving Anisotropic Problem

Is Bunyamin Suryo, Maureen Clerc


The electroencephalography (EEG) is a non-invasive technique to study electrical brain activity (while the brain is performing a cognitive task). The electrical brain activity is a complex process of electrical propagation because the brain structure is an incredibly complex structure. This complex structure leads to different conductivity properties in terms of its magnitude and orientation, called anisotropic conductivity. Using Maxwell's equations, electrical brain activity has been studied intensively. For simplification, the quasistatic Maxwell’s equations are used to model the electrical brain activity and it leads to deal with a Poisson’s equation. In this research, a feasibility study of using a new method, called Operator Splitting Method (OSM), to solve anisotropic 2-Dimensional (2D) Poisson’s equation is performed. Freeware of the finite element method (FEM), FreeFEM++, is employed to build matrices used in the OSM algorithm. The OSM algorithm which is written in Matlab is then tested to solve anisotropic 2D Laplace’s equation and anisotropic Poisson’s equation with the dipolar source. Afterward, the OSM solutions are validated by using an exact solution and a direct numerical solution. By using L2-Error Norm, the convergence rate of the OSM algorithm is then analyzed. Some numerical experiments have been performed to test the performance of the OSM algorithm. The OSM solution of anisotropic 2D Laplace’s equation coincides with the exact and direct numerical solution of the problem. For anisotropic 2D Poisson’s equation with dipolar source, some similar results have been obtained too. The pattern of the OSM solutions is similar to the pattern of direct numerical solutions of the problem. The results arise a hope to attempt to implement the OSM algorithm for more complex problems such as a realistic human head model.


eeg, numeric; maxwell's equations; anisotropic conductivity

Full Text:



Clerc, M and Papadopoulo, T., Inverse problems in func tional brain imaging, Reference material for the MVA Master 2 class, January 11, 2010.

L. A. Geddes and L. E. Baker, “The specific resistance of biological material—a compendium of data for the biomedical engineer and physiologist,” Medical and bi ological engineering, vol. 5, no. 3, pp. 271–293, 1967.

J. Haueisen, Methods of numerical field calculation for neuromagnetic source localization. Shaker Verlag, 1996.

S. Rush and D. A. Driscoll, “Current distribution in the brain from surface electrodes,” Anesthesia & Analgesia, vol. 47, no. 6, pp. 717–723, 1968.

M. Akhtari, H. Bryant, A. Mamelak, E. Flynn, L. Heller, J. Shih, M. Mandelkem, A. Matlachov, D. Ranken, E. Best, et al., “Conductivities of three-layer live hu man skull,” Brain topography, vol. 14, no. 3, pp. 151– 167, 2002.

P. W. Nicholson, “Specific impedance of cerebral white matter,” Experimental neurology, vol. 13, no. 4, pp. 386–401, 1965.

K. A. Awada, D. R. Jackson, J. T. Williams, D. R. Wilton, S. B. Baumann, and A. C. Papanicolaou, “Com putational aspects of finite element modeling in eeg source localization,” IEEE Transactions on biomedical engineering, vol. 44, no. 8, pp. 736–752, 1997.


DOI: http://dx.doi.org/10.12962/j25807471.v2i2.6397

Creative Commons License
JMES The International Journal of Mechanical Engineering and Sciences by Lembaga Penelitian dan Pengabdian kepada Masyarakat (LPPM) ITS is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.
Based on a work at https://iptek.its.ac.id/index.php/jmes.