Search

The aim of this paper is to introduce a new mathematical model that simulates myocardial blood perfusion that accounts for multiscale and multiphysics features. Our model incorporates cardiac electrophysiology, active and passive mechanics, hemodynamics, valve modeling, and a multicompartment Darcy model of perfusion. We consider a fully coupled electromechanical model of the left heart that provides input for a fully coupled Navier–Stokes–Darcy model for myocardial perfusion. The fluid dynamics problem is modeled in a left heart geometry that includes large epicardial coronaries, while the multicompartment Darcy model is set in a biventricular myocardium. Using a realistic and detailed cardiac geometry, our simulations demonstrate the biophysical fidelity of our model in describing cardiac perfusion. Specifically, we successfully validate the model reliability by comparing in-silico coronary flow rates and average myocardial blood flow with clinically established values ranges reported in relevant literature. Additionally, we investigate the impact of a regurgitant aortic valve on myocardial perfusion, and our results indicate a reduction in myocardial perfusion due to blood flow taken away by the left ventricle during diastole. To the best of our knowledge, our work represents the first instance where electromechanics, hemodynamics, and perfusion are integrated into a single computational framework. [ABSTRACT FROM AUTHOR]

The aim of this paper is to introduce a new mathematical model that simulates myocardial blood perfusion that accounts for multiscale and multiphysics features. Our model incorporates cardiac electrophysiology, active and passive mechanics, hemodynamics, valve modeling, and a multicompartment Darcy model of perfusion. We consider a fully coupled electromechanical model of the left heart that provides input for a fully coupled Navier–Stokes–Darcy Model for myocardial perfusion. The fluid dynamics problem is modeled in a left heart geometry that includes large epicardial coronaries, while the multicompartment Darcy model is set in a biventricular myocardium. Using a realistic and detailed cardiac geometry, our simulations demonstrate the biophysical fidelity of our model in describing cardiac perfusion. Specifically, we successfully validate the model reliability by comparing in-silico coronary flow rates and average myocardial blood flow with clinically established values ranges reported in relevant literature. Additionally, we investigate the impact of a regurgitant aortic valve on myocardial perfusion, and our results indicate a reduction in myocardial perfusion due to blood flow taken away by the left ventricle during diastole. To the best of our knowledge, our work represents the first instance where electromechanics, hemodynamics, and perfusion are integrated into a single computational framework. [ABSTRACT FROM AUTHOR]

Bicuspid aortic valve (BAV) predisposes to aortic aneurysms with a high prevalence. A first hypothesis for this phenomenon is related to fibrillin deficiency (genetic hypothesis). The present article focused on a complementary, hemodynamic hypothesis stating that it is the peculiar fluid dynamics of blood in the ascending aorta of patients with BAV configurations that leads to aneurysm formation. To corroborate this hypothesis, a parametric study was performed based on numerical simulations of ascending aorta hemodynamics with different configurations of orifice area and valve orientation. The resulting wall shear stress (WSS) distributions and degree of asymmetry of the blood jet were investigated, and surrogate indices introduced. The results showed that WSS was more pronounced in subjects with BAV morphologies, also in the nonstenotic case. In particular, a maximum WSS of 3 Pa was found (vs. 1.5 Pa in subjects with a tricuspid configuration). It is localized at the mid-ascending aorta, the segment more prone to dilate as shown by the index related to maximum WSS (0.869 in BAV vs. 0.322 in tricuspid). Moreover, the asymmetry of the blood flow was found larger for decreasing valve area, the related index at mid-ascending aorta being more than three times higher than that found for tricuspid configuration (0.70 vs. 0.20). Further, WSS and flow asymmetry were higher also at the sinus of Valsalva and sinotubolar junction for a latero-lateral (LL) BAV configuration in keeping with the clinical observation on association between BAV configurations and different aortic aneurysm morphology. These findings may help explain the higher risk of aneurysm formation in BAV patients. The proposed indices will require validation prior to application in clinical settings. [ABSTRACT FROM AUTHOR]

In bicuspid aortic valve (BAV) disease, the role of genetic and hemodynamic factors influencing ascending aortic pathology is controversial. To test the effect of BAV geometry on ascending aortic flow, a finite element analysis was undertaken. A surface model of aortic root and ascending aorta was obtained from magnetic resonance images of patients with BAV and tricuspid aortic valve using segmentation facilities of the image processing code Vascular Modeling Toolkit (developed at the Mario Negri Institute). Analytical models of bicuspid (antero-posterior [AP], type 1 and latero-lateral, type 2 commissures) and tricuspid orifices were mathematically defined and turned into a volumetric mesh of linear tetrahedra for computational fluid dynamics simulations. Numerical simulations were performed with the finite element code LifeV. Flow velocity fields were assessed for four levels: aortic annulus, sinus of Valsalva, sinotubular junction, and ascending aorta. Comparison of finite element analysis of bicuspid and tricuspid aortic valve showed different blood flow velocity pattern. Flow in bicuspid configurations showed asymmetrical distribution of velocity field toward the convexity of mid-ascending aorta returning symmetrical in distal ascending aorta. On the contrary, tricuspid flow was symmetrical in each aortic segment. Comparing type 1 BAV with type 2 BAV, more pronounced recirculation zones were noticed in the latter. Finally, we found that in both BAV configurations, maximum wall shear stress is highly localized at the convex portion of the mid-ascending aorta level. Comparison between models showed asymmetrical and higher flow velocity in bicuspid models, in particular in the AP configuration. Asymmetry was more pronounced at the aortic level known to be more exposed to aneurysm formation in bicuspid patients. This supports the hypothesis that hemodynamic factors may contribute to ascending aortic pathology in this subset of patients. [ABSTRACT FROM AUTHOR]

In this work, we performed a computational image-based study of blood dynamics in the whole left heart, both in a healthy subject and in a patient with mitral valve regurgitation. We elaborated multi-series cine-MRI with the aim of reconstructing the geometry and the corresponding motion of left ventricle, left atrium, mitral and aortic valves, and aortic root of the subjects. This allowed us to prescribe such motion to computational blood dynamics simulations where, for the first time, the whole left heart motion of the subject is considered, allowing us to obtain reliable subject-specific information. The final aim is to investigate and compare between the subjects the occurrence of turbulence and the risk of hemolysis and of thrombi formation. In particular, we modeled blood with the Navier–Stokes equations in the arbitrary Lagrangian–Eulerian framework, with a large eddy simulation model to describe the transition to turbulence and a resistive method to manage the valve dynamics, and we used a finite element discretization implemented in an in-house code for the numerical solution. [ABSTRACT FROM AUTHOR]

Bicuspid aortic valve (BAV) is associated with aortic dilatation and aneurysm. Several studies evidenced an eccentric systolic flow in ascending aorta associated with increased wall shear stresses (WSS) and the occurrence of an helical systolic flow. This study seeks to elucidate the connections between jet asymmetry and helical flow in patients with normally functioning BAV and dilated ascending aorta. We performed a computational parametric study by varying, for a patient-specific geometry, the valve area and the flow rate entering the aorta and drawing also a tricuspid valve (TAV). We considered also phase-contrast magnetic resonance imaging of four BAV and TAV patients. Measurement of normalized flow asymmetry index, systolic WSS and of a new index (positive helix fraction, PHF) quantifying the presence of a single a single helical flow were performed. In our computation, BAV cases featured higher values of all indices with respect to TAV in both numerical and imaged-based results. Moreover, all indices increased with decreasing valve area and/or with increasing flow rate. This allowed to separate the BAV and TAV cases with respect to the jet asymmetry, WSS localization and helical flow. Interestingly, these results were obtained without modeling the leaflets. [ABSTRACT FROM AUTHOR]