Computational fluid dynamics of blood flow in an idealized left human heart

Fluid-structure interaction analysis of the thromboembolic risk in the left atrial appendage under atrial fibrillation: Effect of hemodynamics and morphological features

BACKGROUND

Complications of atrial fibrillation (AF) include ischemic events originating within the left atrial appendage (LAA), a protrusion of the left atrium with variable morphological characteristics. The role of the patient specific morphology and pathological haemodynamics on the risk of ischemia remains unclear.

METHODS

This work performs a comparative assessment of the hemodynamic parameters among patient-specific LAA morphologies through fluid-structure interaction computational analyses. Three LAA models per each of the four commons patient-specific morphological families (chicken wing, cactus, windsock, and cauliflower) were analysed. Mechanical properties of the tissue were based on experimental uniaxial tests on a young pig's heart. Boundary conditions were imposed based on clinical assessments of filling and emptying volumes. Sinus rhythm and atrial fibrillation operative conditions were simulated and analysed.

RESULTS

For each model, the effect of morphological and functional parameters, such as the number of trabeculae and LAA stroke volume, over the hemodynamics established into the appendage was analysed. Comparison between results obtained in healthy and diseased conditions suggested the introduction of a new parameter to quantify the risk of thrombosis, here called blood stasis factor (BSF). This is defined as the LAA surface area which permanently experiences levels of shear strain rate inferior to a threshold value, set to 5 s-1 (BSF5).

CONCLUSIONS

This work suggests that the current morphological classification is unsuitable to evaluate the probability of thrombus formation. However, hemodynamic parameters easy to determine from clinical examinations, such as normalised stroke volume, LAA orifice flow rate and presence of extensive trabeculations can identify departures from healthy hemodynamics in AF and support a more systematic stratification of the thromboembolic risk.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

The Role of Patient-Specific Morphological Features of the Left Atrial Appendage on the Thromboembolic Risk Under Atrial Fibrillation

Background

A large majority of thrombi causing ischemic complications under atrial fibrillation (AF) originate in the left atrial appendage (LAA), an anatomical structure departing from the left atrium, characterized by a large morphological variability between individuals. This work analyses the hemodynamics simulated for different patient-specific models of LAA by means of computational fluid-structure interaction studies, modeling the effect of the changes in contractility and shape resulting from AF.

Methods

Three operating conditions were analyzed: sinus rhythm, acute atrial fibrillation, and chronic atrial fibrillation. These were simulated on four patient-specific LAA morphologies, each associated with one of the main morphological variants identified from the common classification: chicken wing, cactus, windsock, and cauliflower. Active contractility of the wall muscle was calibrated on the basis of clinical evaluations of the filling and emptying volumes, and boundary conditions were imposed on the fluid to replicate physiological and pathological atrial pressures, typical of the various operating conditions.

Results

The LAA volume and shear strain rates were analyzed over time and space for the different models. Globally, under AF conditions, all models were well aligned in terms of shear strain rate values and predicted levels of risk. Regions of low shear rate, typically associated with a higher risk of a clot, appeared to be promoted by sudden bends and focused at the trabecule and the lobes. These become substantially more pronounced and extended with AF, especially under acute conditions.

Conclusion

This work clarifies the role of active and passive contraction on the healthy hemodynamics in the LAA, analyzing the hemodynamic effect of AF that promotes clot formation. The study indicates that local LAA topological features are more directly associated with a thromboembolic risk than the global shape of the appendage, suggesting that more effective classification criteria should be identified.

Effect of the Alterations in Contractility and Morphology Produced by Atrial Fibrillation on the Thrombosis Potential of the Left Atrial Appendage

Atrial fibrillation (AF) is a common arrhythmia mainly affecting the elderly population, which can lead to serious complications such as stroke, ischaemic attack and vascular dementia. These problems are caused by thrombi which mostly originate in the left atrial appendage (LAA), a small muscular sac protruding from left atrium. The abnormal heart rhythm associated with AF results in alterations in the heart muscle contractions and in some reshaping of the cardiac chambers. This study aims to verify if and how these physiological changes can establish hemodynamic conditions in the LAA promoting thrombus formation, by means of computational fluid dynamic (CFD) analyses. In particular, sinus and fibrillation contractility was replicated by applying wall velocity/motion to models based on healthy and dilated idealized shapes of the left atrium with a common LAA morphology. The models were analyzed and compared in terms of shear strain rate (SSR) and vorticity, which are hemodynamic parameters directly associated with thrombogenicity. The study clearly indicates that the alterations in contractility and morphology associated with AF pathologies play a primary role in establishing hemodynamic conditions which promote higher incidence of ischaemic events, consistently with the clinical evidence. In particular, in the analyzed models, the impairment in contractility determined a decrease in SSR of about 50%, whilst the chamber pathological dilatation contributed to a 30% reduction, indicating increased risk of clot formation. The equivalent rigid wall model was characterized by SSR values about one order of magnitude smaller than in the contractile models, and substantially different vortical behavior, suggesting that analyses based on rigid chambers, although common in the literature, are inadequate to provide realistic results on the LAA hemodynamics.

A parallel non-nested two-level domain decomposition method for simulating blood flows in cerebral artery of stroke patient

Numerical simulation of blood flows in patient-specific arteries can be useful for the understanding of vascular diseases, as well as for surgery planning. In this paper, we simulate blood flows in the full cerebral artery of stroke patients. To accurately resolve the flow in this rather complex geometry with stenosis is challenging and it is also important to obtain the results in a short amount of computing time so that the simulation can be used in pre- and/or post-surgery planning. For this purpose, we introduce a highly scalable, parallel non-nested two-level domain decomposition method for the three-dimensional unsteady incompressible Navier-Stokes equations with an impedance outlet boundary condition. The problem is discretized with a stabilized finite element method on unstructured meshes in space and a fully implicit method in time, and the large nonlinear systems are solved by a preconditioned parallel Newton-Krylov method with a two-level Schwarz method. The key component of the method is a non-nested coarse problem solved using a subset of processor cores and its solution is interpolated to the fine space using radial basis functions. To validate and verify the proposed algorithm and its highly parallel implementation, we consider a case with available clinical data and show that the computed result matches with the measured data. Further numerical experiments indicate that the proposed method works well for realistic geometry and parameters of a full size cerebral artery of an adult stroke patient on a supercomputers with thousands of processor cores.

An image-based computational hemodynamics study of the Systolic Anterior Motion of the mitral valve

Systolic Anterior Motion (SAM) of the mitral valve - often associated with Hypertrophic Obstructive Cardiomyopathy (HOCM) - is a cardiac pathology in which a functional subaortic stenosis is induced during systole by the mitral leaflets partially obstructing the outflow tract of the left ventricle. Its assessment by diagnostic tests is often difficult, possibly underestimating its severity and thus increasing the risk of heart failure. In this paper, we propose a new computational pipeline, based on cardiac cine Magnetic Resonance Imaging (cine-MRI) data, for the assessment of SAM. The pipeline encompasses image processing of the left ventricle and the mitral valve, and numerical investigation of cardiac hemodynamics by means of Computational Fluid Dynamics (CFD) in a moving domain with image-based prescribed displacement. Patient-specific geometry and motion of the left ventricle are considered in view of an Arbitrary Lagrangian-Eulerian approach for CFD, while the reconstructed mitral valve is immersed in the computational domain by means of a resistive method. We assess clinically relevant flow and pressure indicators in a parametric study for different degrees of SAM severity, in order to provide a better quantitative evaluation of the pathological condition. Moreover, we provide specific indications for its possible surgical treatment, i.e. septal myectomy.