A modeling framework for contact, adhesion and mechano-transduction between excitable deformable cells

Cardiac myocytes are the fundamental cells composing the heart muscle. The propagation of electric signals and chemical quantities through them is responsible for their nonlinear contraction and dilatation. In this study, a theoretical model and a finite element formulation are proposed for the simulation of adhesive contact interactions between myocytes across the so-called gap junctions. A multi-field interface constitutive law is proposed for their description, integrating the adhesive and contact mechanical response with their electrophysiological behavior. From the computational point of view, the initial and boundary value problem is formulated as a structure-structure interaction problem, which leads to a straightforward implementation amenable for parallel computations. Numerical tests are conducted on different couples of myocytes, characterized by different shapes related to their stages of growth, capturing the experimental response. The proposed framework is expected to have impact on the understanding how imperfect mechano-transduction could lead to emergent pathological responses.Comment: 31 pages, 17 figure

CFD-based functional imaging for arteries: in vitro validation

De l’imagerie fonctionnelle pour les vaisseaux est développée à partir de données médicales morphologiques (IRM 4D) et hémodynamiques (IRM par contraste de phase dans les plans d’entrée-sortie). Les données fonctionnelles pertinentes (champ de vitesse, frottement pariétal, gradient de pression,…) sont alors calculées en simulant l’écoulement compatible avec les données médicales. On présente les résultats obtenus dans la phase de validation in vitro de cette technique sur un fantôme de crosse aortique

From medical images to individualized cardiac mechanics: A Physiome approach

Cardiac mechanics is a branch of science that deals with forces, kinematics, and material properties of the heart, which is valuable for clinical applications and physiological studies. Although anatomical and biomechanical experiments are necessary to provide the fundamental knowledge of cardiac mechanics, the invasive nature of the procedures limits their further applicability. In consequence, noninvasive alternatives are required, and cardiac images provide an excellent source of subject-specific and in vivo information. Noninvasive and individualized cardiac mechanical studies can be achieved through coupling general physiological models derived from invasive experiments with subject-specific information extracted from medical images. Nevertheless, as data extracted from images are gross, sparse, or noisy, and do not directly provide the information of interest in general, the couplings between models and measurements are complicated inverse problems with numerous issues need to be carefully considered. The goal of this research is to develop a noninvasive framework for studying individualized cardiac mechanics through systematic coupling between cardiac physiological models and medical images according to their respective merits. More specifically, nonlinear state-space filtering frameworks for recovering individualized cardiac deformation and local material parameters of realistic nonlinear constitutive laws have been proposed. To ensure the physiological meaningfulness, clinical relevance, and computational feasibility of the frameworks, five key issues have to be properly addressed, including the cardiac physiological model, the heart representation in the computational environment, the information extraction from cardiac images, the coupling between models and image information, and also the computational complexity. For the cardiac physiological model, a cardiac physiome model tailored for cardiac image analysis has been proposed to provide a macroscopic physiological foundation for the study. For the heart representation, a meshfree method has been adopted to facilitate implementations and spatial accuracy refinements. For the information extraction from cardiac images, a registration method based on free-form deformation has been adopted for robust motion tracking. For the coupling between models and images, state-space filtering has been applied to systematically couple the models with the measurements. For the computational complexity, a mode superposition approach has been adopted to project the system into an equivalent mathematical space with much fewer dimensions for computationally feasible filtering. Experiments were performed on both synthetic and clinical data to verify the proposed frameworks

Integrated Heart - Coupling multiscale and multiphysics models for the simulation of the cardiac function

Mathematical modelling of the human heart and its function can expand our understanding of various cardiac diseases, which remain the most common cause of death in the developed world. Like other physiological systems, the heart can be understood as a complex multiscale system involving interacting phenomena at the molecular, cellular, tissue, and organ levels. This article addresses the numerical modelling of many aspects of heart function, including the interaction of the cardiac electrophysiology system with contractile muscle tissue, the sub-cellular activation-contraction mechanisms, as well as the hemodynamics inside the heart chambers. Resolution of each of these sub-systems requires separate mathematical analysis and specially developed numerical algorithms, which we review in detail. By using specific sub-systems as examples, we also look at systemic stability, and explain for example how physiological concepts such as microscopic force generation in cardiac muscle cells, translate to coupled systems of differential equations, and how their stability properties influence the choice of numerical coupling algorithms. Several numerical examples illustrate three fundamental challenges of developing multiphysics and multiscale numerical models for simulating heart function, namely: (i) the correct upscaling from single-cell models to the entire cardiac muscle, (ii) the proper coupling of electrophysiology and tissue mechanics to simulate electromechanical feedback, and (iii) the stable simulation of ventricular hemodynamics during rapid valve opening and closure

Real-time 3D reconstruction of non-rigid shapes with a single moving camera

© . This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/This paper describes a real-time sequential method to simultaneously recover the camera motion and the 3D shape of deformable objects from a calibrated monocular video. For this purpose, we consider the Navier-Cauchy equations used in 3D linear elasticity and solved by finite elements, to model the time-varying shape per frame. These equations are embedded in an extended Kalman filter, resulting in sequential Bayesian estimation approach. We represent the shape, with unknown material properties, as a combination of elastic elements whose nodal points correspond to salient points in the image. The global rigidity of the shape is encoded by a stiffness matrix, computed after assembling each of these elements. With this piecewise model, we can linearly relate the 3D displacements with the 3D acting forces that cause the object deformation, assumed to be normally distributed. While standard finite-element-method techniques require imposing boundary conditions to solve the resulting linear system, in this work we eliminate this requirement by modeling the compliance matrix with a generalized pseudoinverse that enforces a pre-fixed rank. Our framework also ensures surface continuity without the need for a post-processing step to stitch all the piecewise reconstructions into a global smooth shape. We present experimental results using both synthetic and real videos for different scenarios ranging from isometric to elastic deformations. We also show the consistency of the estimation with respect to 3D ground truth data, include several experiments assessing robustness against artifacts and finally, provide an experimental validation of our performance in real time at frame rate for small mapsPeer ReviewedPostprint (author's final draft

Isogeometric Kirchhoff–Love shell formulations for general hyperelastic materials

We present formulations for compressible and incompressible hyperelastic thin shells which can use general 3D constitutive models. The necessary plane stress condition is enforced analytically for incompressible materials and iteratively for compressible materials. The thickness stretch is statically condensed and the shell kinematics are completely described by the first and second fundamental forms of the midsurface. We use C1-continuous isogeometric discretizations to build the numerical models. Numerical tests, including structural dynamics simulations of a bioprosthetic heart valve, show the good performance and applicability of the presented methods

Spatio-Temporal Tensor Decomposition of a Polyaffine Motion Model for a Better Analysis of Pathological Left Ventricular Dynamics

International audienceGiven that heart disease can cause abnormal motion dynamics over the cardiac cycle, which can then affect cardiac function, understanding and quantifying cardiac motion can provide insight for clinicians to aid in diagnosis, therapy planning, as well as to determine the prognosis for a given patient. The goal of this paper is to extract population-specific cardiac motion patterns from 3D displacements in order to firstly identify the mean motion behaviour in a population and secondly to describe pathology-specific motion patterns in terms of the spatial and temporal aspects of the motion. Since there are common motion patterns observed in patients suffering from the same condition, extracting these patterns can lead towards a better understanding of a disease. Quantifying cardiac motion at a population level is not a simple task since images can vary widely in terms of image quality, size, resolution and pose. To overcome this, we analyse the parameters obtained from a cardiac-specific Polyaffine motion tracking algorithm, which are aligned both spatially and temporally to a common reference space. Once all parameters are aligned, different subjects and different populations can be compared and analysed in the space of Polyaffine transformations by projecting the transformations to a reduced-order subspace in which dominant motion patterns in each population can be extracted and analysed. Using tensor decomposition allows the spatial and temporal aspects to be decoupled in order to study the different components individually. The proposed method was validated on healthy volunteers and Tetralogy of Fallot patients according to known spatial andtemporal behaviour for each population. A key advantage of the proposed method is the ability to regenerate motion sequences from the respective models, thus the models can be visualised in terms of the full motion, which allows for better understanding of the motion dynamics of different populations