A macro finite-element formulation for cardiac electrophysiology simulations using hybrid unstructured grids.

Electrical activity in cardiac tissue can be described by the bidomain equations whose solution for large-scale simulations still remains a computational challenge. Therefore, improvements in the discrete formulation of the problem, which decrease computational and/or memory demands are highly desirable. In this study, we propose a novel technique for computing shape functions of finite elements (FEs). The technique generates macro FEs (MFEs) based on the local decomposition of elements into tetrahedral subelements with linear shape functions. Such an approach necessitates the direct use of hybrid meshes (HMs) composed of different types of elements. MFEs are compared to classic standard FEs with respect to accuracy and RAM memory usage under different scenarios of cardiac modeling, including bidomain and monodomain simulations in 2-D and 3-D for simple and complex tissue geometries. In problems with analytical solutions, MFEs displayed the same numerical accuracy of standard linear triangular and tetrahedral elements. In propagation simulations, conduction velocity and activation times agreed very well with those computed with standard FEs. However, MFEs offer a significant decrease in memory requirements. We conclude that HMs composed of MFEs are well suited for solving problems in cardiac computational electrophysiology.

Finite-element simulation of wave propagation in periodic piezoelectric SAW structures.

Many surface acoustic wave (SAW) devices consist of quasiperiodic structures that are designed by successive repetition of a base cell. The precise numerical simulation of such devices, including all physical effects, is currently beyond the capacity of high-end computation. Therefore, we have to restrict the numerical analysis to the periodic substructure. By using the finite-element method (FEM), this can be done by introducing periodic boundary conditions (PBCs) at special artificial boundaries. To be able to describe the complete dispersion behavior of waves, including damping effects, the PBC has to be able to model each mode that can be excited within the periodic structure. Therefore, the condition used for the PBCs must hold for each phase and amplitude difference existing at periodic boundaries. Based on the Floquet theorem, our two newly developed PBC algorithms allow the calculation of both, the phase and the amplitude coefficients of the wave. In the first part of this paper we describe the basic theory of the PBCs. Based on the FEM, we develop two different methods that deliver the same results but have totally different numerical properties and, therefore, allow the use of problem-adapted solvers. Further on, we show how to compute the charge distribution of periodic SAW structures with the aid of the new PBCs. In the second part, we compare the measured and simulated dispersion behavior of waves propagating on periodic SAW structures for two different piezoelectric substrates. Then we compare measured and simulated input admittances of structures similar to SAW resonators.