Towards Finite Element Modeling of the Accoustics of Human Head

Authors: Maciej Paszynski, Leszek Demkowicz, Jason Kurtz

Abstract:

In this paper, a new formulation for acoustics coupled with linear elasticity is presented. The primary objective of the work is to develop a three dimensional hp adaptive finite element method code destinated for modeling of acoustics of human head. The code will have numerous applications e.g. in designing hearing protection devices for individuals working in high noise environments. The presented work is in the preliminary stage. The variational formulation has been implemented and tested on a sequence of meshes with concentric multi-layer spheres, with material data representing the tissue (the brain), skull and the air. Thus, an efficient solver for coupled elasticity/acoustics problems has been developed, and tested on high contrast material data representing the human head.

Keywords: finite element method, acoustics, coupled problems, biomechanics

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1057433

References:

[1] H. Kopka and P. W. Daly, A Guide to LATEX, 3rd ed. Harlow, England: Addison-Wesley, 1999.
[2] Y.C. Chang and L. Demkowicz, Scattering on a spherical shell problem. Comparison of 3D elasticity and Kirchhoff shell theory results. Computer Assisted Mechanics and Engineering Science, 2:207229, 1995.
[3] L. Demkowicz. Computing with hp Finite Elements. I.One- and Two- Dimensional Elliptic and Maxwell Problems. Chapman & Hall/CRC Press, Taylor and Francis, 2006.
[4] L. Demkowicz, J. Kurtz, D. Pardo, M. Paszy'nski, W. Rachowicz, and A. Zdunek, Computing with hp Finite Elements. II Frontiers: Three- Dimensional Elliptic and Maxwell Problems with Applications. Chapman & Hall/CRC Press, Taylor and Francis, 2007.
[5] P. Geng, J.T. Oden, and R.A. van de Geijn. A parallel multifrontal algorithm and its implementation. Computer Methods in Applied Mechanics and Engineering, 149:289301, 1997.
[6] A. Majda. Compressible Fluid Flow and Systems of Conservation Laws in Several Space Variables, volume 53 of Applied Mathematical Sciences, Springer-Verlag, New York, 1984.
[7] Ch. Michler, L. Demkowicz, J. Kurtz, and D. Pardo. Improving the performance of perfectly matched layers ny means of hp-adaptivity. ICES- Report, 06-17, The University of Texas at Austin, 2006.
[8] MUMPS: a multifrontal massively parallel sparse direct solver, http://www.enseeiht.fr/lima/apo/MUMPS/.
[9] M. Paszy'nski, Performance of Multi Level Parallel Direct Solver for hp Finite Element Method, Lecture Notes in Computer Science, 4967:1303- 1312, 2007.
[10] M. Paszy'nski, D. Pardo, C. Torres-Verdin, L. Demkowicz, and V. Calo, A Multi-Level Direct Substructuring Multi-Frontal Parallel Direct Solver for hp-Finite Element Method, ICES-Report, 07-33, The University of Texas at Austin, 2007
[11] P.-O. Persson and G. Strang. A simple mesh generator in matlab. SIAM Review, 46(2):329 345, 2004.
[12] W. Schroeder, K. Martin and B. Lorensen. The Visualization Toolkit An Object-Oriented Approach To 3D Graphics, 3rd Edition. Kitware, Inc., (http://www.kitware.com).
[13] T. Walsh and L. Demkowicz. A parallel multifrontal solver for hpadaptive finite elements. Technical Report 1, TICAM, The University of Texas at Austin, 1999.
[14] T. Walsh, L. Demkowicz and R. Charles. Boundary element modeling of the external human auditory system. J. Acoust. Soc. Am., 115(3), 2004.
[15] Y. Zhang, C. Bajaj, and B-S Sohn. 3D finite element meshing from imaging data. The special issue of Computer Methods in Applied Me- chanics and Engineering (CMAME) on Unstructured Mesh Generation, 194(48-49):50835106, 2005.