Parallel Newton-Krylov-BDDC and FETI-DP deluxe solvers for implicit time discretizations of the cardiac Bidomain equations

by   Ngoc Mai Monica Huynh, et al.

Two novel parallel Newton-Krylov Balancing Domain Decomposition by Constraints (BDDC) and Dual-Primal Finite Element Tearing and Interconnecting (FETI-DP) solvers are here constructed, analyzed and tested numerically for implicit time discretizations of the three-dimensional Bidomain system of equations. This model represents the most advanced mathematical description of the cardiac bioelectrical activity and it consists of a degenerate system of two non-linear reaction-diffusion partial differential equations (PDEs), coupled with a stiff system of ordinary differential equations (ODEs). A finite element discretization in space and a segregated implicit discretization in time, based on decoupling the PDEs from the ODEs, yields at each time step the solution of a non-linear algebraic system. The Jacobian linear system at each Newton iteration is solved by a Krylov method, accelerated by BDDC or FETI-DP preconditioners, both augmented with the recently introduced deluxe scaling of the dual variables. A polylogarithmic convergence rate bound is proven for the resulting parallel Bidomain solvers. Extensive numerical experiments on linux clusters up to two thousands processors confirm the theoretical estimates, showing that the proposed parallel solvers are scalable and quasi-optimal.



There are no comments yet.


page 5

page 13

page 14


Newton-Krylov-BDDC deluxe solvers for non-symmetric fully implicit time discretizations of the Bidomain model

A novel theoretical convergence rate estimate for a Balancing Domain Dec...

Multirate Partial Differential Equations for the Efficient Simulation of Low-Frequency Problems with Pulsed Excitations

This paper proposes the use of Multirate Partial Differential Equations ...

A BDF2-Semismooth Newton Algorithm for the Numerical Solution of the Bingham Flow with Temperature Dependent Parameters

This paper is devoted to the numerical solution of the non-isothermal in...

Multigrid preconditioners for the mixed finite element dynamical core of the LFRic atmospheric model

Due to the wide separation of time scales in geophysical fluid dynamics,...

A closest point method library for PDEs on surfaces with parallel domain decomposition solvers and preconditioners

The DD-CPM software library provides a set of tools for the discretizati...

Composing Scalable Nonlinear Algebraic Solvers

Most efficient linear solvers use composable algorithmic components, wit...

Parallel matrix-free higher-order finite element solvers for phase-field fracture problems

Phase-field fracture models lead to variational problems that can be wri...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.