Adaptive BDDC preconditioners for the bidomain model on unstructured ventricular finite element meshes

This study aims to develop and numerically analyze adaptive balancing domain decomposition by constraints (BDDC) preconditioners for unstructured finite element discretizations of the Bidomain model of electrocardiology on patient-specific ventricular geometries. The Bidomain model, one of the most comprehensive mathematical representations of the cardiac bioelectrical activity, consists of a system of an elliptic and a parabolic partial differential equation of reaction-diffusion type. These equations govern the propagation of electrical potentials in the cardiac tissue. They are strongly coupled with a stiff system of ordinary differential equations that describe the evolution of ionic currents across the cardiac cell membrane. Minimizing the computational cost of simulating this bioelectrical activity requires the development of efficient and scalable preconditioners for the linear systems resulting from the model’s discretization. BDDC preconditioners are nonoverlapping domain decomposition algorithms that consist of the solution of concurrent local problems on each subdomain plus a global coarse problem, whose unknowns are vertex and edge/face average values. Adaptive BDDC preconditioners represent an evolution of standard BDDC methods, where the coarse problem is enriched by adding further constraints obtained by solving suitable generalized eigenvalue problems on subdomain edges and faces. The novelty of the present study is to analyze the effectiveness of such adaptive BDDC methods for unstructured finite element discretizations of the Bidomain model on patient-specific left ventricular geometries, using modern high-performance computing parallel architectures. Different refined left ventricular meshes were generated, each incorporating fiber data. The results highlight the efficiency and accuracy of the implemented preconditioners, confirming their optimality and scalability on CPUs architectures.