scholarly journals A Mathematical Spline-Based Model of Cardiac Left Ventricle Anatomy and Morphology

Author(s):  
Sergei Pravdin

Computer simulation of normal and diseased human heart activity requires a 3D anatomical model of the myocardium, including myofibres. For clinical applications, such a model has to be constructed based on routine methods of cardiac visualisation such as sonography. Symmetrical models are shown to be too rigid, so an analytical non-symmetrical model with enough flexibility is necessary. Based on previously made anatomical models of the left ventricle, we propose a new, much more flexible spline-based analytical model. The model is fully described and verified based on DT-MRI data. We show a way to construct it on the basis of sonography data. To use this model in further physiological simulations, we propose a numerical method to utilise finite differences in solving the reaction-diffusion problem together with an example of scroll wave dynamics simulation.

Author(s):  
Sergei Pravdin

Computer simulation of normal and diseased human heart activity requires a 3D anatomical model of the myocardium, including myofibres. For clinical applications, such a model has to be constructed based on routine methods of cardiac visualisation such as sonography. Symmetrical models are shown to be too rigid, so an analytical non-symmetrical model with enough flexibility is necessary. Based on previously made anatomical models of the left ventricle, we propose a new, much more flexible spline-based analytical model. The model is fully described and verified against DT-MRI data. We show a way to construct it on the basis of sonography data. To use this model in further physiological simulations, we propose a numerical method to utilise finite differences in solving the reaction-diffusion problem together with an example of scroll wave dynamics simulation.


Mathematics ◽  
2020 ◽  
Vol 8 (5) ◽  
pp. 776
Author(s):  
Sergey Pravdin ◽  
Pavel Konovalov ◽  
Hans Dierckx ◽  
Olga Solovyova ◽  
Alexander V. Panfilov

Rotating spiral waves of electrical excitation underlie many dangerous cardiac arrhythmias. The heterogeneity of myocardium is one of the factors that affects the dynamics of such waves. In this paper, we present results of our simulations for scroll wave dynamics in a heterogeneous model of the human left ventricle with analytical anatomically based representation of the geometry and anisotropy. We used a set of 18 coupled differential equations developed by ten Tusscher and Panfilov (TP06 model) which describes human ventricular cells based on their measured biophysical properties. We found that apicobasal heterogeneity dramatically changes the scroll wave dynamics. In the homogeneous model, the scroll wave annihilates at the base, but the moderate heterogeneity causes the wave to move to the apex and then continuously rotates around it. The rotation speed increased with the degree of the heterogeneity. However, for large heterogeneity, we observed formation of additional wavebreaks and the onset of complex spatio-temporal patterns. Transmural heterogeneity did not change the dynamics and decreased the lifetime of the scroll wave with an increase in heterogeneity. Results of our numerical experiments show that the apex may be a preferable location of the scroll wave, which may be important for development of clinical interventions.


Author(s):  
Konstantin S. Ushenin ◽  
Sergei F. Pravdin ◽  
Yuliya S. Alueva ◽  
Tatyana V. Chumarnaya ◽  
Olga E. Solovyova

AbstractThis paper presents first results on the dynamics of filaments of scroll waves of myocardium excitation obtained for personalized models of the left ventricle of the human heart. The paper describes a mathematical model of the left ventricle of the human heart and its electrical activity, numerical methods for the model calculation within the framework of computer implementation, and also the method of personalization of the model according to data of ultrasound examination. We found that regardless of the starting point of wave the filaments of wave drift along a spiral to the attractor located near the apex of the ventricle. The attractor position was essentially different in models constructed from the data of patients without identified pathology and those for patients with an increased left ventricular cavity.


2020 ◽  
Vol 28 (3) ◽  
pp. 147-160
Author(s):  
Andrea Bonito ◽  
Diane Guignard ◽  
Ashley R. Zhang

AbstractWe consider the numerical approximation of the spectral fractional diffusion problem based on the so called Balakrishnan representation. The latter consists of an improper integral approximated via quadratures. At each quadrature point, a reaction–diffusion problem must be approximated and is the method bottle neck. In this work, we propose to reduce the computational cost using a reduced basis strategy allowing for a fast evaluation of the reaction–diffusion problems. The reduced basis does not depend on the fractional power s for 0 < smin ⩽ s ⩽ smax < 1. It is built offline once for all and used online irrespectively of the fractional power. We analyze the reduced basis strategy and show its exponential convergence. The analytical results are illustrated with insightful numerical experiments.


2016 ◽  
Vol 16 (4) ◽  
pp. 609-631 ◽  
Author(s):  
Immanuel Anjam ◽  
Dirk Pauly

AbstractThe results of this contribution are derived in the framework of functional type a posteriori error estimates. The error is measured in a combined norm which takes into account both the primal and dual variables denoted by x and y, respectively. Our first main result is an error equality for all equations of the class ${\mathrm{A}^{*}\mathrm{A}x+x=f}$ or in mixed formulation ${\mathrm{A}^{*}y+x=f}$, ${\mathrm{A}x=y}$, where the exact solution $(x,y)$ is in $D(\mathrm{A})\times D(\mathrm{A}^{*})$. Here ${\mathrm{A}}$ is a linear, densely defined and closed (usually a differential) operator and ${\mathrm{A}^{*}}$ its adjoint. In this paper we deal with very conforming mixed approximations, i.e., we assume that the approximation ${(\tilde{x},\tilde{y})}$ belongs to ${D(\mathrm{A})\times D(\mathrm{A}^{*})}$. In order to obtain the exact global error value of this approximation one only needs the problem data and the mixed approximation itself, i.e., we have the equality$\lvert x-\tilde{x}\rvert^{2}+\lvert\mathrm{A}(x-\tilde{x})\rvert^{2}+\lvert y-% \tilde{y}\rvert^{2}+\lvert\mathrm{A}^{*}(y-\tilde{y})\rvert^{2}=\mathcal{M}(% \tilde{x},\tilde{y}),$where ${\mathcal{M}(\tilde{x},\tilde{y}):=\lvert f-\tilde{x}-\mathrm{A}^{*}\tilde{y}% \rvert^{2}+\lvert\tilde{y}-\mathrm{A}\tilde{x}\rvert^{2}}$ contains only known data. Our second main result is an error estimate for all equations of the class ${\mathrm{A}^{*}\mathrm{A}x+ix=f}$ or in mixed formulation ${\mathrm{A}^{*}y+ix=f}$, ${\mathrm{A}x=y}$, where i is the imaginary unit. For this problem we have the two-sided estimate$\frac{\sqrt{2}}{\sqrt{2}+1}\mathcal{M}_{i}(\tilde{x},\tilde{y})\leq\lvert x-% \tilde{x}\rvert^{2}+\lvert\mathrm{A}(x-\tilde{x})\rvert^{2}+\lvert y-\tilde{y}% \rvert^{2}+\lvert\mathrm{A}^{*}(y-\tilde{y})\rvert^{2}\leq\frac{\sqrt{2}}{% \sqrt{2}-1}\mathcal{M}_{i}(\tilde{x},\tilde{y}),$where ${\mathcal{M}_{i}(\tilde{x},\tilde{y}):=\lvert f-i\tilde{x}-\mathrm{A}^{*}% \tilde{y}\rvert^{2}+\lvert\tilde{y}-\mathrm{A}\tilde{x}\rvert^{2}}$ contains only known data. We will point out a motivation for the study of the latter problems by time discretizations or time-harmonic ansatz of linear partial differential equations and we will present an extensive list of applications including the reaction-diffusion problem and the eddy current problem.


Sign in / Sign up

Export Citation Format

Share Document