An oscillation-free fully partitioned scheme for the numerical modeling of cardiac active mechanics

by   Francesco Regazzoni, et al.

In silico models of cardiac electromechanics couple together mathematical models describing different physics. One instance is represented by the model describing the generation of active force, coupled with the one of tissue mechanics. For the numerical solution of the coupled model, partitioned schemes, that foresee the sequential solution of the two subproblems, are often used. However, this approach may be unstable. For this reason, the coupled model is commonly solved as a unique system using Newton type algorithms, at the price, however, of high computational costs. In light of this motivation, in this paper we propose a new numerical scheme, that is numerically stable and accurate, yet within a fully partitioned (i.e. segregated) framework. Specifically, we introduce, with respect to standard segregated scheme, a numerically consistent stabilization term, capable of removing the nonphysical oscillations otherwise present in the numerical solution of the commonly used segregated scheme. Our new method is derived moving from a physics-based analysis on the microscale energetics of the force generation dynamics. By considering a model problem of active mechanics we prove that the proposed scheme is unconditionally absolutely stable (i.e. it is stable for any time step size), unlike the standard segregated scheme, and we also provide an interpretation of the scheme as a fractional step method. We show, by means of several numerical tests, that the proposed stabilization term successfully removes the nonphysical numerical oscillations characterizing the non stabilized segregated scheme solution. Our numerical tests are carried out for several force generation models available in the literature, namely the Niederer-Hunter-Smith model, the model by Land and coworkers, and the mean-field force generation model that we have recently proposed. Finally, we apply the proposed scheme [...]



There are no comments yet.


page 26

page 28


A positivity-preserving and convergent numerical scheme for the binary fluid-surfactant system

In this paper, we develop a first order (in time) numerical scheme for t...

An Energy Stable Linear Diffusive Crank-Nicolson Scheme for the Cahn-Hilliard Gradient Flow

We propose and analyze a linearly stabilized semi-implicit diffusive Cra...

Coupled time discretisation of dynamic damage models at small strains

The dynamic damage model in viscoelastic materials in Kelvin-Voigt rheol...

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...

The Newton Scheme for Deep Learning

We introduce a neural network (NN) strictly governed by Newton's Law, wi...
This week in AI

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