Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userElectromechanical modeling of human ventricles with ischemic cardiomyopathy: numerical simulations in sinus rhythm and under arrhythmia
51
0
0.0
(
0
)
Added by
Matteo Salvador
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We developed a novel patient-specific computational model for the numerical simulation of ventricular electromechanics in patients with ischemic cardiomyopathy (ICM). This model reproduces the activity both in sinus rhythm (SR) and in ventricular tachycardia (VT). The presence of scars, grey zones and non-remodeled regions of the myocardium is accounted for by the introduction of a spatially heterogeneous coefficient in the 3D electromechanics model. This 3D electromechanics model is firstly coupled with a 2-element Windkessel afterload model to fit the pressure-volume (PV) loop of a patient-specific left ventricle (LV) with ICM in SR. Then, we employ the coupling with a 0D closed-loop circulation model to analyze a VT circuit over multiple heartbeats on the same LV. We highlight similarities and differences on the solutions obtained by the electrophysiology model and those of the electromechanics model, while considering different scenarios for the circulatory system. We observe that very different parametrizations of the circulation model induce the same hemodynamical considerations for the patient at hand. Specifically, we classify this VT as unstable. We conclude by stressing the importance of combining electrophysiological, mechanical and hemodynamical models to provide relevant clinical indicators in how arrhythmias evolve and can potentially lead to sudden cardiac death.
rate research
Read More
Heart Failure is a major component of healthcare expenditure and a leading cause of mortality worldwide. Despite higher inter-rater variability, endomyocardial biopsy (EMB) is still regarded as the standard technique, used to identify the cause (e.g. ischemic or non-ischemic cardiomyopathy, coronary artery disease, myocardial infarction etc.) of unexplained heart failure. In this paper, we focus on identifying cardiomyopathy as ischemic or non-ischemic. For this, we propose and implement a new unified architecture comprising CNN (inception-V3 model) and bidirectional LSTM (BiLSTM) with self-attention mechanism to predict the ischemic or non-ischemic to classify cardiomyopathy using histopathological images. The proposed model is based on self-attention that implicitly focuses on the information outputted from the hidden layers of BiLSTM. Through our results we demonstrate that this framework carries a high learning capacity and is able to improve the classification performance.
We propose a novel algorithm for the approximation of surface-quasi geostrophic (SQG) flows modeled by a nonlinear partial differential equation coupling transport and fractional diffusion phenomena. The time discretization consists of an explicit strong-stability-preserving three-stage Runge-Kutta method while a flux-corrected-transport (FCT) method coupled with Dunford-Taylor representations of fractional operators is advocated for the space discretization. Standard continuous piecewise linear finite elements are employed and the algorithm does not have restrictions on the mesh structure nor on the computational domain. In the inviscid case, we show that the resulting scheme satisfies a discrete maximum principle property under a standard CFL condition and observe, in practice, its second-order accuracy in space. The algorithm successfully approximates several benchmarks with sharp transitions and fine structures typical of SQG flows. In addition, theoretical Kolmogorov energy decay rates are observed on a freely decaying atmospheric turbulence simulation.
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 [...]
The numerical integration of an analytical function $f(x)$ using a finite set of equidistant points can be performed by quadrature formulas like the Newton-Cotes. Unlike Gaussian quadrature formulas however, higher-order Newton-Cotes formulas are not stable, limiting the usable order of such formulas. Existing work showed that by the use of orthogonal polynomials, stable high-order quadrature formulas with equidistant points can be developed. We improve upon such work by making use of (orthogonal) Gram polynomials and deriving an iterative algorithm, together allowing us to reduce the space-complexity of the original algorithm significantly.
In the paper [Hainaut, D. and Colwell, D.B., A structural model for credit risk with switching processes and synchronous jumps, The European Journal of Finance44(33) (4238):3262-3284], the authors exploit a synchronous-jump regime-switching model to compute the default probability of a publicly-traded company. Here, we first generalize the proposed Levy model to a more general setting of tempered stable processes recently introduced into the finance literature. Based on the singularity of the resulting partial integro-differential operator, we propose a general framework based on strictly positive-definite functions to de-singularize the operator. We then analyze an efficient meshfree collocation method based on radial basis functions to approximate the solution of the corresponding system of partial integro-differential equations arising from the structural credit risk model. We show that under some regularity assumptions, our proposed method naturally de-sinularizes the problem in the tempered stable case. Numerical results of applying the method on some standard examples from the literature confirm the accuracy of our theoretical results and numerical algorithm.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
