A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs
Centre for Product Design and Manufacturing, Indian Institute of Science, Bangalore, Karnataka 560012, India
arXiv:1311.4533 [cs.CE], (18 Nov 2013)
@article{2013arXiv1311.4533K,
author={Kumara P}, K.},
title={"{A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs}"},
journal={ArXiv e-prints},
archivePrefix={"arXiv"},
eprint={1311.4533},
primaryClass={"cs.CE"},
keywords={Computer Science – Computational Engineering, Finance, and Science, Computer Science – Distributed, Parallel, and Cluster Computing, Computer Science – Mathematical Software, Physics – Computational Physics, Physics – Medical Physics},
year={2013},
month={nov},
adsurl={http://adsabs.harvard.edu/abs/2013arXiv1311.4533K},
adsnote={Provided by the SAO/NASA Astrophysics Data System}
}
In this work, first a Fortran code is developed for three dimensional linear elastostatics using constant boundary elements; the code is based on a MATLAB code developed by the author earlier. Next, the code is parallelized using BLACS, MPI, and ScaLAPACK. Later, the parallelized code is used to demonstrate the usefulness of the Boundary Element Method (BEM) as applied to the realtime computational simulation of biological organs, while focusing on the speed and accuracy offered by BEM. A computer cluster is used in this part of the work. The commercial software package ANSYS is used to obtain the `exact’ solution against which the solution from BEM is compared; analytical solutions, wherever available, are also used to establish the accuracy of BEM. A pig liver is the biological organ considered. Next, instead of the computer cluster, a Graphics Processing Unit (GPU) is used as the parallel hardware. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature. Also, a serial MATLAB code, and both serial and parallel versions of a Fortran code, which can solve three dimensional (3D) linear elastostatic problems using constant boundary elements, are provided as supplementary files that can be freely downloaded.
November 19, 2013 by hgpu
