Finite Element Based Solution of Laplace's Equation Applied to Electrical Activity of the Human Body

Authors

  • Zainab T. Baqer Department of Electrical Engineering/ College of Engineering/ University of Baghdad

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

Download data is not yet available.

References

[1] Christopher R. Johnson."Adaptive Finite Element and Local Regularization Methods for the Inverse ECG Problem". Biology Society 17th Annual International Conference, pages 209–210. IEEE Press. Center for Scientific Computing and Imaging, Department of Computer Science, University of Utah, Salt Lake City, Utah 84112 US, 1995.
[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

2010-12-01

How to Cite

Baqer, Z. T. (2010). Finite Element Based Solution of Laplace’s Equation Applied to Electrical Activity of the Human Body. Al-Khwarizmi Engineering Journal, 6(4), 37–51. Retrieved from https://alkej.uobaghdad.edu.iq/index.php/alkej/article/view/511