Mathematical modelling of active contraction in isolated cardiomyocytes

We investigate the interaction of intracellular calcium spatio-temporal variations with the self-sustained contractions in cardiac myocytes. A consistent mathematical model is presented considering a hyperelastic description of the passive mechanical properties of the cell, combined with an active-strain framework to explain the active shortening of myocytes and its coupling with cytosolic and sarcoplasmic calcium dynamics. A finite element method based on a Taylor-Hood discretization is employed to approximate the nonlinear elasticity equations, whereas the calcium concentration and mechanical activation variables are discretized by piecewise linear finite elements. Several numerical tests illustrate the ability of the model in predicting key experimentally established characteristics including: (i) calcium propagation patterns and contractility, (ii) the influence of boundary conditions and cell shape on the onset of structural and active anisotropy and (iii) the high localized stress distributions at the focal adhesions. Besides, they also highlight the potential of the method in elucidating some important subcellular mechanisms affecting, e.g. cardiac repolarizatio

An adaptive finite element method for the modeling of the equilibrium of red blood cells

International audienceThis contribution is concerned with a the numerical modeling of an isolated red blood cell (RBC), and more generally of phospholipid membranes. We propose an adaptive Eulerian finite element approximation, based on the level set method, of a shape optimization problem arising in the study of RBC's equilibrium. We simulate the equilibrium shapes that minimize the elastic bending energy under prescribed constraints of fixed volume and surface area. An anisotropic mesh adaptation technique is used in the vicinity of the cell's membrane to enhance the robustness of the method. Efficient time and spatial discretizations are considered and implemented. We address in detail the main features of the proposed method and finally we report several numerical experiments in the two-dimensional and the three-dimensional axisymmetric cases. The effectiveness of the numerical method is further demonstrated through numerical comparisons with semi-analytical solutions provided by a reduced order model

A Finite Element Approach For Modeling Biomembranes In Incompressible Power-Law Flow

We present a numerical method to model the dynamics of inextensible biomembranes in a quasi-Newtonian incompressible flow, which better describes hemorheology in the small vasculature. We consider a level set model for the fluid-membrane coupling, while the local inextensibility condition is relaxed by introducing a penalty term. The penalty method is straightforward to implement from any Navier-Stokes/level set solver and allows substantial computational savings over a mixed formulation. A standard Galerkin finite element framework is used with an arbitrarily high order polynomial approximation for better accuracy in computing the bending force. The PDE system is solved using a partitioned strongly coupled scheme based on Crank-Nicolson time integration. Numerical experiments are provided to validate and assess the main features of the method

Numerical Approach Based on the Composition of One-Step Time-Integration Schemes For Highly Deformable Interfaces

In this work, we propose a numerical approach for simulations of large deformations of interfaces in a level set framework. To obtain a fast and viable numerical solution in both time and space, temporal discretization is based on the composition of one-step methods exhibiting higher orders and stability, especially in the case of stiff problems with strongly oscillatory solutions. Numerical results are provided in the case of ordinary and partial differential equations to show the main features and demonstrate the performance of the method. Convergence properties and efficiency in terms of computational cost are also investigated

A Necklace Model for Vesicles Simulations in 2D

International audienceThe aim of this paper is to propose a new numerical model to simulate 2D vesicles interacting with a newtonian fluid. The inextensible membrane is modeled by a chain of circular rigid particles which are maintained in cohesion by using two different type of forces. First, a spring force is imposed between neighboring particles in the chain. Second, in order to model the bending of the membrane, each triplet of successive particles is submitted to an angular force. Numerical simulations of vesicles in shear flow have been run using Finite Element Method and the FreeFem++[1] software. Exploring different ratios of inner and outer viscosities, we recover the well known "Tank-Treading" and "Tumbling" motions predicted by theory and experiments. Moreover, for the first time, 2D simulations of the "Vacillating-Breathing" regime predicted by theory in [2] and observed experimentally in [3] are done without special ingredient like for example thermal fluctuations used in [4]

Mathematical modelling of active contraction in isolated cardiomyocytes

We investigate the interaction of intracellular calcium spatio-temporal variations with the self-sustained contractions in cardiac myocytes. A consistent mathematical model is presented considering a hyperelastic description of the passive mechanical properties of the cell, combined with an active-strain framework to explain the active shortening of myocytes and its coupling with cytosolic and sarcoplasmic calcium dynamics. A finite element method based on a Taylor-Hood discretization is employed to approximate the nonlinear elasticity equations, whereas the calcium concentration and mechanical activation variables are discretized by piecewise linear finite elements. Several numerical tests illustrate the ability of the model in predicting key experimentally established characteristics including: (i) calcium propagation patterns and contractility, (ii) the influence of boundary conditions and cell shape on the onset of structural and active anisotropy and (iii) the high localized stress distributions at the focal adhesions. Besides, they also highlight the potential of the method in elucidating some important subcellular mechanisms affecting, e.g. cardiac repolarization

Reynolds number effect on the dissipation function in wall-bounded flows

International audienceThe evolution with Reynolds number of the dissipation function, normalized by wall variables, is investigated using direct numerical simulation (DNS) databases for incompressible turbulent Poiseuille flow in a plane channel, at friction Reynolds numbers up to Reτ = 2000. DNS results show that the mean part, directly dissipated by the mean flow, reaches a constant value while the turbulent part, converted into turbulent kinetic energy before being dissipated, follows a logarithmic law. This result shows that the logarithmic law of friction can be obtained without any assumption on the mean velocity distribution. The proposed law is in good agreement with experimental results in plane-channel and boundary layer flows