A 3D boundary optimal control for the bidomain-bath system modeling the thoracic shock therapy for cardiac defibrillation

This work is dedicated to study the cardiac defibrillation problem by using an optimal thoracic electroshock treatment. The problem is formulated as an optimal control problem in a 3D domain surrounded by the bath and including the heart. The control corresponds to the thoracic electroshock and the model describing the electrical activity in the heart is the bidomain model. The bidomain model is coupled with the quasi-static Maxwell's equation to consider the effect of an external bathing medium. The existence and uniqueness of a weak solution for the direct problem is assessed as well as the existence of a weak solution for the adjoint problem. The numerical discretization is realized using a finite element method for the spatial discretization and linearly implicit Runge–Kutta methods for the temporal discretization of the partial differential equations. The numerical results are demonstrated for the termination of re-entry waves.