Finite Element Based Solution of Laplace's Equation Applied to Electrical Activity of the Human Body
Abstract
Computer models are used in the study of electrocardiography to provide insight into physiological phenomena that are difficult to measure in the lab or in a clinical environment.
The electrocardiogram is an important tool for the clinician in that it changes characteristically in a number of pathological conditions. Many illnesses can be detected by this measurement. By simulating the electrical activity of the heart one obtains a quantitative relationship between the electrocardiogram and different anomalies.
Because of the inhomogeneous fibrous structure of the heart and the irregular geometries of the body, finite element method is used for studying the electrical properties of the heart.
This work describes the implementation of the Conjugate Gradient iterative method for the solution of large linear equation systems resulting from the finite element method. A diagonal Jacobi preconditioner is used in order to accelerate the convergence. Gaussian elimination is also implemented and compared with the Precondition Conjugate Gradient (PCG) method and with the iterative method. Different types of matrix storage schemes are implemented such as the Compressed Sparse Row (CSR) to achieve better performance. In order to demonstrate the validity of the finite element analysis, the technique is adopted to solve Laplace's equation that describes the electrical activity of the human body with Dirichlet and Neumann boundary conditions. An automatic mesh generator is built using C++ programming language. Initially a complete finite element program is built to solve Laplace's equation. The same accuracy is obtained using these methods. The results show that the CSR format reduces computation time compared to the order format. The PCG method is better for the solution of large linear system (sparse matrices) than the Gaussian Elimination and back substitution method, while Gaussian elimination is better than iterative method.
Downloads
References
[2] Julian Johnson, "Wireless Electrocardiogram"
[3] Wael Namat Al_Sahar, "The design and implementation of a PC based ECG recording and interpretation system", M.Sc Thesis, University of Technology, Electronic Engineering Dept, 1999.
[4] Pilkington TC, Plonsey R: "Engineering Contributions to Biophysical Electrocardiography", IEEE Press, John Wiley, New York, 1982.
[5] Tirupathi R. Chandrupatla and Ashok D. Belegundu "Introduction to Finite Elements in Engineering", copyright ©1991 by Prentice-Halls Inc.
[6] [Matthew N. O. Sadiku "Numerical Techniques in Electromagnetics"; 2nd edition, © 2001 by CRC Press LLC.
[7] C.O.Moretti, J.B. Cavalcante Neto, T.N. Bittencourt and L.F. Martha,"A Parallel Environment for Three-Dimensional Finite Element Method Analysis", 2000.
[8] O. Skipa, D. Farina, C. Kaltwasser, O.D¨ossel, W. R. Bauer "Fast Interaction ECG Simulation and Optimization-Based Reconstruction of Depolarization" Institute for Biomedical Technique, University of Karlsruhe (TH), Germany Institute for Biophysics, University at Wurzburg, Deutschland [email protected], 2004.
[9] M.Okajima, K.Doniwa, T. Yamana, K. Ohta, N.Suzumura, A. Iwata and K. Kamiya. "Individualized Modeling of the Heart and Torso for Finite Element Method to be used in Forward Calculation of Body Surface Map", IEEE, University, Toyoake, Nagoya Area 470-11, Japan, 1989.
[10] Antonio d' Acierno, Giuseppe De Pietro, Antonio Giordano, "A Parallel Implementation of the Conjugate Gradient Method on the Meiko CS-2", IRSIP _ CNR, Napoli (Italy), 1996.
[11] Ananth Grama, Anshul Gupta, George Karypis and Vipin Kumar "Introduction to Parallel Computing" 2nd edition published 2003.
Downloads
Published
Issue
Section
License
Copyright: Open Access authors retain the copyrights of their papers, and all open access articles are distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution and reproduction in any medium, provided that the original work is properly cited. The use of general descriptive names, trade names, trademarks, and so forth in this publication, even if not specifically identified, does not imply that these names are not protected by the relevant laws and regulations. While the advice and information in this journal are believed to be true and accurate on the date of its going to press, neither the authors, the editors, nor the publisher can accept any legal responsibility for any errors or omissions that may be made. The publisher makes no warranty, express or implied, with respect to the material contained herein.
